|Home | About | Journals | Submit | Contact Us | Français|
Single molecules in a tunnel junction can now be interrogated reliably using chemically-functionalized electrodes. Monitoring stochastic bonding fluctuations between a ligand bound to one electrode and its target bound to a second electrode (“tethered molecule-pair” configuration) gives insight into the nature of the intermolecular bonding at a single molecule-pair level, and defines the requirements for reproducible tunneling data. Simulations show that there is an instability in the tunnel gap at large currents, and this results in a multiplicity of contacts with a corresponding spread in the measured currents. At small currents (i.e. large gaps) the gap is stable, and functionalizing a pair of electrodes with recognition reagents (the “free analyte” configuration) can generate a distinct tunneling signal when an analyte molecule is trapped in the gap. This opens up a new interface between chemistry and electronics with immediate implications for rapid sequencing of single DNA molecules.
Reading chemical identity at the single atom or single molecule level has long been a goal of the electron tunneling community. Recognition tunneling achieves that goal, serves as a remarkable probe of chemical bonds at the single molecule level and enables a new type of single molecule electronic sensor.
Atomic identification of sites on clean semiconductor surfaces in ultrahigh vacuum (UHV) was achieved by atomic resolution tunneling spectroscopy in the mid 1980's1. This approach has been extended to a number of molecular systems,2 including DNA.3, 4 Potential distributions in UHV are reproducible enough for precise molecular identification in many cases. In particular, the Hipps-Mazur lab has shown that there is a strong correlation between the bias at which features are observed in the tunneling spectroscopy and the reduction potential of the molecules.2
The situation is much more difficult in solution where solvent reorganization and an unknown potential distribution (owing to poorly characterized screening effects) complicate the analysis,5 though the problem can be simplified when chemically well-defined contacts are used between the target molecules and the metal electrodes.6, 7
In an ingenious experiment, Tao showed that tunneling transmission could be gated by controlling the potential of a redox active molecule relative to a standard reference electrode.8 This solves the problem of an unknown potential distribution, and single molecules can be identified, but only if they contain a redox-active species with a distinctive reduction or oxidation feature at an accessible potential. A completely new approach was introduced by Ohshiro and Umezawa in 2006.9 They demonstrated that STM images of molecular adsorbates showed chemical contrast if the scanning probe was functionalized with a molecule that bound to specific targets on the surface. In their work, a gold STM probe was functionalized with a thiolated DNA base, and images were made of monolayers of DNA bases, or peptide nucleic acid (PNA, a DNA analog with a peptide backbone) bonded to a gold substrate. Enhanced contrast was found when the base on the probe was the Watson-Crick hydrogen bonding complement of the bases on the substrate. Analysis of the interaction between the functionalized STM probe and the target molecules shows that it is dominated by mechanical adhesion between the probe and target molecules, with essentially no contribution by a “chemically-sensitive” tunnel current.10 It appears that the STM probe sticks to its complementary target for a longer time as the tip is scanned, resulting in more charge transfer between probe and substrate. Thus the tunnel current serves as a method for monitoring mechanical adhesion, and the method more closely resembles friction-force chemical sensing11 or single-molecule force spectroscopy12 as usually implemented with an atomic force microscope. Nonetheless, a variant of the Ohshiro-Umezawa method does permit identification of molecular species on a surface13 and has been used to “read” the composition of DNA oligomers, albeit with a resolution limited to blocks of about ten bases, owing to the strong mechanical interactions between the probe and surface.14
We have recently refined the method further, and it is now at the point where true single molecule reads are possible with the exquisite degree of recognition that comes from matching sets of hydrogen-bond donors and acceptors on probe and target.15 Rather than move the probe relative to the substrate, we fix a tunnel gap and monitor the stochastic fluctuations in bonding as the target molecule binds, unbinds and rebinds the probe molecule (Figure 1A). Such stochastic switching has been observed before as a consequence of the lability of thiol contacts between a molecule and a gold surface,16 and was used as the basis of an ingenious method for measuring single-molecule conductivity introduced by Nichols and Haiss.17, 18 This gives rise to a characteristic “telegraph” noise in the tunnel current as illustrated in Figure 1B. The gap is set to an initial “baseline” conductance, GBL, the current increasing by an amount corresponding to the “on conductance”, ΔGON when a molecular bridge spans the gap. In the first implementation of “Recognition Tunneling” (which we call the “tethered molecule-pair” configuration, Figure 1A, B) we have measured the conductance of Watson-Crick base pairs using a base bonded to a probe and a nucleoside bonded to a surface, and monitoring the telegraph noise generated both by breaking of the hydrogen bonds and by breaking of the gold-thiol contacts used to hold the molecules to electrodes.15 These measurements permit direct probing of bonds between single molecule pairs, measuring both their lifetimes and their electronic transport properties.
In a further refinement of Recognition Tunneling (which we call the “free analyte” configuration, Figure 1C), we have functionalized a pair of tunneling electrodes with molecules that present both hydrogen bond donors and acceptors, setting the gap to be large enough for target molecules to diffuse into it and bridge the sensing molecules on the surface to the sensing molecules on the probe via a network of hydrogen bonds. Differences in the pattern of hydrogen bonding lead to different levels of tunnel current as targets bind transiently in the gap (Figure 1D). These signals can be distinctive enough to allow the identity of the four DNA bases to be read with a confidence level that is greater than 0.6 on a single molecule read.19 This technique could be adapted to read small drug molecules, metabolites, amino acids, small peptides and so on.
It is the purpose of this review to provide a simple background to Recognition Tunneling, an outline of hydrogen-bond mediated tunneling and a description of experimental methods, as well as a summary of the results obtained to date. We provide, in addition, new modeling that we have carried out to understand the way in which the telegraph noise distribution changes with the size of the tunnel gap.
We begin with a brief overview of the phenomenon of electron tunneling. Technical details are available in a number of theoretical reviews.7, 20-24. Figure 2A shows a pair of metal electrodes, separated by a gap, L. The electrons that transport current lie near the highest occupied state at the Fermi energy, EF. The potential barrier that retains electrons within the metal is V volts so the work function, ϕ, is given by ϕ = V – EF. In the absence of a state in the gap (i.e., without the state represented by the thick line ΔE above the Fermi energy) the tunnel conductance is given approximately by
where G0is the quantum of conductance, μS, ϕ is the work function in units of eV and L is the tunnel gap in Å. Thus, when an automatically-sharp point contact is made, the conductance is approximately 77.5 μS, falling by a factor of 10 for each Ångstrom that the gap is opened for a metal like gold (ϕ ≈ 5eV). This very rapid decay with distance is what makes tunneling such a localized probe and therefore capable of reading single bases in a DNA polymer.
If an electronic eigenstate exists in the gap, then the effective barrier is reduced to approximately ΔE (thick solid line in Figure 2A). β, the inverse electronic decay length is then given by if this state connects the left and right electrodes. A real molecule spanning the gap will consist of molecular orbitals often well described in terms of linear combinations of the atomic orbitals of the constituent atoms and the hopping matrix elements between adjacent orbitals. Figure 2B is a schematic representation of this situation showing atomic states closest to the Fermi energy (E1 - E4) with hopping matrix elements HL1, H12....H34, H4R. Thus, the states of the system, ψn, are molecular orbitals, as modified by interactions with the left (HL1) and right (H4R) electrodes.
Quantum mechanics requires that the probability amplitude for every path that connects the left electrode to the right electrode be added. The paths are all those for which molecular orbitals interact with both electrodes. The probability that an electron is propagated to the opposite electrode is then calculated from the square of this sum of probability amplitudes. Mathematically, the transmission is calculated with a Green's function:
Here, L and R represent states at the energy E on the left and right electrodes and E0n is the eigenenergy of the nth eigenstate of the system and ε is an infinitesimal. The overall conductance is calculated from , where T is the transmission function and the integral extends over energies above and below the metal Fermi-levels by an amount equal to the applied bias “V.” The transmission function involves an (L,R) matrix element of |U G U|2, where U is the “perturbation” connecting the L and R electrodes to the molecules wired between them.
Equation 2 casts the result as a two terminal Büttiker expression25 for the current in terms of a transmission function. The means to calculating the transmission function is through the quantum mechanical Green's function. The earliest development the quantum mechanical formulations in terms of local (atomic) orbitals was constructed by Caroli et al.26 There have been many applications of this formulation and it has been developed using a variety of local orbital techniques.
The key ingredient in all formulations is the Green's function operator, (E - H)-1, where H is the electronic Hamiltonian. In the time domain, this operator is the time evolution operator e–iHt/ħ. It contains all the information concerning the coherent propagation of the electrons through the molecules. The Hamiltonian H describes the left-metal, molecule(s), right-metal and the couplings, U, between the left-metal and molecule(s) and between the right-metal and molecule(s). In its simplest form the propagation through the molecule contains elements of the Green's function of the form shown in Equation 2.
The full theory can be cast as a scattering problem – an electron comes in (say) from the left-metal, is scattered by the molecule(s) manifest by tunneling, and finally enters the right-metal. This scattering formulation treats the problem as three non-interacting systems (with net Hamiltonian H0) and with an interaction/coupling perturbation U that connects the three systems and allows transmission. The complete system Hamiltonian H is H0 + U. The coupling U is included to all orders where the interacting Greens function G(E)=(E-H)-1 is written in terms of the non-interacting Green's functions, G0(E) = (E-H0)-1. Dyson's equation relates the two, G=G0+G0UG. The scattering results in transitions from one metal to the next by a Fermi's golden rule formula where the effective perturbation is UG(E)U, rather than just U as in a textbook formula for Fermi's golden rule. This is the approach taken by Tomfohr et al.27 where a full discussion of the theory and how it is implemented in local orbital density functional theory can be found.
The presence of G(E) in the transmission includes the energy levels of the molecule and naturally introduces the possibility of resonances and multiple reflections. Keldysh-like28 formulations of non-equilibrium Green's functions go further in that they incorporate effects of the self-consistent shifts of the atomic orbitals that are produced by the electric field introduced by the bias.29, 30
The form of the Green's function (Equation 2) has some important experimental implications. If the connection is severed (U=0), the current drops to zero. The Green's function of the molecule is dressed by the contacts broadening out the electronic levels. Under most conditions, the Green's function propagator between the two contacts decays exponentially, which is a signature of electron tunneling. The key point here is that the states that mediate tunneling are those that overlap both the left and right electrodes. Any “break” in the quantum coupling between the electrodes will not allow the two terms in the numerator of equation 2 to be non-zero simultaneously. This has two consequences. The first is that molecules must be coupled to the electrodes, and the most robust way to do this direct chemical bonding between the ends of the molecule and the electrodes.6 The second consequence is that the coupling strengths throughout the system (the Hmn) must be large enough not to be disrupted by thermal fluctuations.
The discussion that follows equation 1 appears to imply that β → 0 as ΔE → 0 (i.e., ϕ → 0 in equation 1). In solid-state systems, this is one condition for resonant tunneling (the other being tight coupling to the electrodes 31). In polarizable molecular systems the corresponding process is electrochemical reduction, the charge being trapped on an acceptor level by molecular and solvent reorganization. This is generally an undesirable outcome: Many redox processes are irreversible as a consequence of the generation of reactive species that, in turn, can alter the readout system itself by reacting with recognition molecules. Secondly, redox processes on single molecules can be quite heterogeneous if the environment is heterogeneous.32
Given that redox processes are to be avoided, how close to the Fermi level can a state of the combined electrode-molecule-electrode system be? There are two ways to address this, both leading to the same answer. Acceptor states on the molecule will be distributed over a range of energies comparable to thermal energy, kBT. Thus we require |E — E0n| > kBT to avoid redox processes. However, we also require that the ‘weakest link’ in the system, for quantum coherence to be maintained across the gap. Since the smallest energy splitting (2ΔE if we assume the Fermi level is in the middle of the molecular gap so that the closest orbital is shifted by ΔE from the Fermi level) is approximately equal to , we conclude that we require ΔE >> kBT. Thus β cannot be zero.
The foregoing discussion allows us to estimate the largest distance over which single molecule recognition tunneling reads may be made as follows: Let us require . Let us assume a minimum read time of 1ms and that our electronics has an input capacitance of 1 pF (it is difficult to do better than this). Then the current-to-voltage conversion resistor cannot be greater than 1 GΩ (so that RC = 1ms) so the Johnson noise (given by ) in a bandwidth, Δf = 1 kHz is about 0.1 pA. Assuming a maximum applied bias of 0.5 V (less than this will often be required to avoid redox processes) the smallest useful gap conductance is then 0.2 pS. According to equation 1, this corresponds to a value of βL=19.7 which, β ≥ 0.5 Å-1 yields a maximum value for L of about 4 nm. This is the maximum possible value with energy levels in just the right place with respect to the Fermi level. In less favorable situatios, this maximum distace will be smaller. Thus Recognition Tunneling is restricted to combinations of analytes and reading reagents that are smaller than 4 nm. Taking the target to be no more than 3 nm diameter would limit the upper read size to proteins of no more than 104 Daltons (assuming that the electronic structure was just right).
Hydrogen bond strengths lie between those of covalent and Van der Waal's bonds.33 With bond strengths in solution of just a few times kBT, hydrogen bonds make an ideal “molecular Velcro” for making transient contacts to analytes for applications like DNA sequencing. Can reasonable tunnel currents be sustained through such “weak” bonds? In fact, this is not the right way to ask the question, for the bond strength is related to the overall lowering of electronic energy in occupied states, whereas the tunneling rate will be determined by the energy difference between the Fermi energy of the metal and the state closest to it that spans the tunnel gap (c.f. equation 2). We need to determine, for a hydrogen bonded system, the β to use in equation 1.
One simple way to determine an upper limit on β is with complex band structure.34 This technique uses standard (fast) codes for determining the band structure of a periodic system, but uses a complex wavevector, keeping track of the results for both the real and imaginary components. States with an imaginary wavevector are forbidden in truly infinite systems but are allowed for finite (e.g. molecular) systems. The magnitude of the imaginary wavevector is the inverse decay constant over 2, β/2, for that state (β/2 being a real variable in a decaying wave-function and β the decay constant for electron density). Assuming that the state that mediates tunneling at the Fermi energy lies halfway between the HOMO and the LUMO, the value of β can be read directly from the imaginary part of the complex band structure. A very simple (and chemically unrealistic) model of an H-bonded system is shown as the line of planar water molecules on the left of Figure 3A. They are arranged in an infinite lattice of lattice constant a, and the band structure calculated with local orbitals is shown to the right.35 The states with imaginary wavevector are plotted to the right, and the value of β midway between the HOMO and LUMO is 0.99 Å-1 (9.9 nm-1). This is very close to the value found in sigma-bonded systems like alkane chains. Thus, this hydrogen bonded system, despite being held together by “weak” hydrogen bonds, propagates electrons as well as an alkane chain. This is not a new result. Studies of optically-excited electron transfer mediated by π-bonds, hydrogen bonds and σ-bonds show that the transfer rate via hydrogen bonds is greater than that via σ-bonds.36
The value of 0.99 Å-1 for β for the water chain is calculated for the whole system comprising a unit cell of the lattice in an infinite chain (Figure 3A). When the β strictly associated with the H-bond itself is calculated from the conductance decay as just the hydrogen bonds are stretched, a much larger value (of about 3 Å-1) is calculated.35 Thus the current should be very sensitive to straining of these bonds.
In reality, the Fermi level alignment is rarely exactly halfway between the HOMO and LUMO. Exact Green's function calculations are required for the entire molecular system as connected to semi-infinite slab electrodes. Structures for the two Watson-Crick base pairs, a G:T “wobble” base pair and a triply-bonded 2-aminoadenine:thymine base pair are shown in Figure 3B. The yellow atoms are the sulfur atoms that were place in a three-fold hollow on Au(111) surfaces (all structures were relaxed prior to calculation of the Green's function).37 The calculated density of states for each type of atom (color coded as explained in the Figure caption) is shown in Figure 3C. The solid lines are projected for atoms in the purines and the dashed lines are for atoms in the pyrimidines. The Fermi level lies close to the LUMO where states are dominated by the pyrimidines. In consequence, the effective β for these molecular systems is quite small and the calculated tunnel conductances quite high.37
Haiss and Nichols pioneered the monitoring of stochastic fluctuations in molecular junctions as a method for determining the single molecule conductance.17,18 In this technique, a molecule with two reactive end groups (thiols) is trapped between two (gold) electrodes. Stochastic fluctuations of the metal-molecule interface16 result in a train of “on”-“off” switching events in the tunnel current (Figure 1B). If a single molecule is trapped, then this current has just two levels (“on” and “off”) and so resembles the telegrapher's Morse code, hence “telegraph noise”. Surprisingly, this method has not been widely adopted yet. It has some advantages over the more popular “break junction” methods.38, 39 Firstly, the signal identifies junctions containing single molecules unambiguously because switching of more than one molecule gives rise to more than two levels in the current. In the case of break junction data, the number of molecules in the tunnel gap has to be inferred statistically. Secondly, the junction remains static, avoiding the constant reforming that occurs with break junctions. And finally, because of the separation of the two electrodes, it is possible to study different molecules on the probe and substrate, desirable for the “tethered molecule-pair” version of Recognition Tunneling. These would be rapidly scrambled in a break junction measurement.
Haiss and Nichols reported that the molecular conductances measured with telegraph noise varies with the gap, and can become multivalued at small gaps. 17, 18 For these reasons, we have found it convenient to make a scatter plot of ΔGON vs. GBL and data for three combinations of basepairing have been presented in this form elsewhere. 15 We have now measured several other types of molecular junction (Chang et al., unpublished) and use some of this data here to show what appear to be rather universal characteristics. Figure 4 shows data for a guanine-cytidine junction (A) and a junction spanned by just an octanedithiol (B). The guanine-cytidine junction data are representative of the two other base-paired systems studied to date. 15 In all these systems, there is a region where ΔGON increases approximately linearly with GBL (i.e., exponentially with decreasing gap) at very small conductances (i.e., large gaps) as shown by the region enclosed by the green box on the lower left corners of Figs. 3A and 3B. As GBL is made larger, measured values of ΔGON spread out over an increasing range of values between a lower “constant” plateau (red dashed boxes) and an upper bound that increase with GBL. Understanding this behavior is critical for the design of Recognition Tunneling readers.
We have simulated a junction containing octanedithiol, using the calculated interaction forces in the junctions to predict the types of structure that are likely to arise as the molecule and electrodes influence each other. As we show below, gaps of different sizes are spanned by configurations that involve distortions of both the molecular bonds and the bonds in the metal contact, as well as the molecule-metal bond. We then calculated the electronic conductance for a number of these optimized structures.
Total energy DFT calculations were carried out using the VASP code,40 within generalized gradient approximation PW9141. We have studied the system consisting of a single octanedithiol molecule sandwiched between two gold electrodes represented by two asymmetric surfaces having a~(3×3) periodicity in the xy-plane. We also imposed periodic boundary conditions in the direction perpendicular to the surface, artificially joining the last two layers of the system. The stretching of the system is simulated by increasing the distance between the two limiting layers by steps of 0.25Å. After each step, the system was allowed to relax toward its configuration of minimum energy; in this relaxation, only the atoms located in the last two layers remain fixed. The convergence criteria for total energy and force were 10-6 eV and 10-2 eV/Å, respectively. Plane-wave functions were expanded on a mesh of 12 Monkhorst-Pack (MP) k-points on the two-dimensional Brillouin zone with a plane-wave cutoff of 290 eV.
Subsequently, we used the optimized atomic structures to calculate the current flowing through molecular junctions. The conductance calculations were carried out using the Green's function formalism42 implemented in a local basis set DFT code (Fireball43, 44) using 64 MP k-points. The results of these simulations are shown in Figure 4C.
The calculated conductances are about an order of magnitude higher than the measured conductances but the general shape of the ΔGON vs. GBL simulated plot (Figure 4C) is strikingly similar to the measured data (Figure 4B): (1) At low values of GBL the simulated data shows a region where ΔGON increases linearly with GBL. (2) Above a threshold value of GBL, data for ΔGON is multivalued, with the maximum value rising with increasing GBL at about he rate observed in the experimental data.
The simulations were carried out for number of contact geometries, only a representative selection of which (A-D) are presented here with their structures illustrated in Figure 5. An important result is that the tunneling current is insensitive to contact geometries at the largest gaps. Figure 4C also plots the “vacuum gap”, ZVAC between the gold atom attached to the molecule and the remainder of the gold cluster (Figure 5E). The region of monotonic dependence of ΔGON vs. GBL corresponds to the “non-contact” region where ZVACvaries rapidly with GBL Figure 5E shows an example of how contact geometry can vary spontaneously over the course of repeated making the breaking of contacts: the lower contact reorganizes on breaking of the top contact. Note how this simulation (and others like it) show clearly that it is an Au-Au bond that breaks and not an Au-sulfur bond. This result is predicted on the (usually-assumed) dissociative adsorption of SH on Au, but there is some evidence the hydrogen remains on the surface.45 The same picture is recapitulated in our simulations of non-dissociative adsorption, save that it is the SH-Au bond that is the weak link that stretches (as opposed to the (Au-S)-S bond), however the (Au-S)-S bond is found to be more stable than the (Au-H)-S bond by about 1 eV in our calculations.
This region of monotonic ΔGON vs. GBL (“non-contact”) sensing may require a flexible molecular junction, for such a region is absent in data obtained with a stiff molecule (dibenzenedithiol -Figure 4D). The DFT calculations capture this behavior for this molecule too. The molecule jumps directly from a non-bonded to a bonded configuration with no intermediate stretching (Hapala and Jelinik, in preparation).
Finally, this analysis suggests that an “equilibrium” configuration occurs just when a molecule first jumps into contact. In a probe of finite radius, we expect that many such configurations arise, and that these correspond to the “plateaus” observed in the plots of ΔGON vs. GBL (red dashed boxes in Figures 4A and 4B). For DNA base-pairs, the average values of conductance obtained from these plateaus lie within a factor two of the predicted values, and follow the theoretically-predicted order of conductances.15
Here, we want to draw a distinction between “electronic bond breaking” (i.e., shutting off tunnel transport) and “mechanical bond breaking (i.e., completely separating the bonding partners). The first observation of these bonding fluctuations by scanning tunneling microscopy 16 showed that the conductance fluctuated while the molecule remained in place. This implies that small fluctuations in bonding can disrupt the electronic coupling essential for tunneling while leaving the molecule still “attached” to the surface. Despite this caveat, the distributions of “on” times clearly contain chemical information, as shown in Figure 6. Figure 6A shows distributions (left sides of the panels) for guanine interacting with cytosine (G:C), 2-aminoadenine interacting with thymine (2AA:T) and adenine interacting with thymine (A:T – structures are shown in Figure 6B). The lifetime distributions are clearly bimodal, with one peak near 45 ms (S-Au) and another near 8 ms (HB). Contrasting these distributions to that measured for an octanedithiol (Figure 6C) shows that the 45 ms features are a consequence of fluctuations of the (S-Au)-Au bonds, so that the 8 ms features must be associated with the hydrogen bonded complexes. The differences in lifetime between the two- and three-hydrogen bonded complexes is small (and the distributions are very wide) but the relative frequency of hydrogen bond breaking is clearly larger relative to (S-Au)-S breaking in the A:T complex than in the G:C and 2AA:T complexes.
The lifetime histograms are resolved as a function of GBL on the right-hand side of Figure 6A. Interestingly, there is little evidence of (S-Au)-Au bond breaking at small GBL in the “non-contact” region. Hydrogen bond breaking dominates initially, but then falls off rapidly with increasing GBL . At small gaps (S-Au)-Au bond breaking appears to be the dominant mode, possibly because of the fact that (S-Au)-Au bond is an intrinsic property of the Au surface. In contrast, the hydrogen bonded components appear to require space to accommodate fluctuations. The intrinsic lifetime of these bonds in a confined space may be very long.
These data clearly demonstrate the ability of this technique to extract bonding information on a single molecule basis. Recent (unpublished) studies of the bonding lifetimes that use force spectroscopy to estimate the off-rate in nanojunctions suggest that the true off rate is very slow. If this is indeed that case, then the system remains “bonded” for long times, with the telegraph noise reporting local fluctuations that interrupt tunnel current, consistent with what was observed by scanning tunneling microscopy.16 The exact nature of these fluctuations is unknown, but it is clear that tunneling will report on even small (atomic-scale) disruptions of the bonding. For example, a motion of just 1 Å will reduce current by a factor 20x.35 On the other hand, the long range electrostatic interactions between charged sites will still contribute significantly to bond strength at much larger separations. Therefore, it is possible for the tunnel current to drop below the background noise while a significant bonding energy remains. This could explain one of the experimental mysteries of Recognition Tunneling. The stability of the tunnel junctions is quite remarkable even in the absence of servo control, currents remaining stable for periods of up to ten seconds. Perhaps the junctions are held together by bonds that are stable on time scales that are much longer than the timescale of the telegraph noise fluctuations.
The foregoing discussion of the “tethered molecule-pair” configuration (Figure 1A) sets the stage for design of a new type of analytical system that traps analytes, the “free analyte” configuration, illustrated in Figure 1B. We recently implemented such a scheme for recognition of DNA nucleosides19 using 4-mercapto benzoic acid as the recognition molecule (here R1=R2 in Figure 1C). In organic solvent (trichlorobenzene) the benzoic acid is neutral, containing a hydrogen bond donor and a hydrogen bond acceptor. Pairs of this reagent can form four distinctive hydrogen bonded complexes with all four DNA bases as shown in Figure 7 A-D.19 The choice of gap size is critical. Any larger, and the smaller bases (like T) will not be trapped on both sides. On the other hand, smaller gaps offer the possibility of multiple binding arrangements for larger bases (like A). At this large gap value (2.5 nm, calculated from measured tunnel decay rates13) the background tunnel current is just 6 pA (at a bias of 0.5V) so that the small gap conductance (12 pS) ensures that we are in the “non-contact” regime where contact geometry variations in conductance appear to be insignificant. In the absence of analyte, the signal (Figure 1E) is also free of the noise spikes that occur in smaller junctions where the benzoic acids molecules can bond with each other across the junction. When nucleosides are injected into the solvent that surrounds the tunnel junction, current spikes are seen almost immediately (Figure 1F). These events often show the telegraph noise (inset) characteristic of binding and unbinding events. This particular reading reagent gives distinctive current reads for each of the four bases, and experimentally measured current distributions are shown in Figure 7G. Reads are not completely separated, but the overlap is such that a correct assignment can be made at the 0.6 confidence level (or better) on the first read. We have also shown that mixtures of bases give clearly resolved signals.19
It is instructive to probe the role of the recognition reagents by comparing reads with, and without them. Figure 8A shows a distribution of currents measured in a somewhat smaller gap (baseline conductance = 20 pS) with a bare gold substrate and probe. The target is deoxyguanosine (the gap had to be decreased somewhat in order to observe current spikes with no functional groups on the electrodes). When one electrode in the gap is functionalized (and set to a baseline conductance of 12 pS) the resultant distribution is dramatically narrowed (Figure 8B). A second small peak at twice the current of the first peak appears, a characteristic feature of junctions like this when two molecule span the gap.7 Functionalizing both electrodes does not decrease the width of the distribution any further. It appears that the target is trapped by the hydrogen bonds, eliminating motional broadening of the tunnel current distribution, but this is not the correct interpretation of the narrowing. This can be seen from the measured distributions of on-state lifetimes shown in Figure 8C. The distributions are very similar for bare electrodes (blue line) and one functionalized electrode (red). Thus the trapping time of molecules in the gap does not depend strongly on the chemical functionalization of the gap. Similar results were obtained for the other three bases.19 It is known that the bases interact quite strongly with gold via the lone-pair electrons on the imines and amines, so the reads taken with bare electrodes almost certainly reflect transient binding of the targets to the gold electrodes. Thus the observed narrowing of the current distribution on functionalizing one of the electrodes reflects a reduction in the range of static bound conformations. Tunneling reads require an interaction between the target molecules and the electrodes. Recognition Tunneling makes this a controlled interaction.
Clearly, Recognition Tunneling is a powerful new tool for investigating bonding in a small junction. It is clearly capable of single molecule detection and might enable a new class of single molecule sensors. As far as DNA sequencing goes, three questions remain: (1) Can the bases be distinguished with reads made in aqueous buffered electrolyte? (2) Can single base resolution be obtained in the context of a DNA molecule where each base is separated by about 0.3 nm? (3) Can chemical modifications of the bases (important as epigenetic coding) be read? Preliminary work in our laboratory suggests that the answer to all three questions is positive. Electrochemical leakage in aqueous buffer is not a problem with good insulation46 of tunneling probes. Signal levels appear to be even larger that those found in organic solvent, despite the expectation that hydrogen bonding between the recognition elements will be weaker in water because of competition from water molecules.
Evidence for single base reads is seen in long runs of telegraph noise as small DNA oligomers with an alternating base sequence are trapped in the gap (Huang, manuscript in preparation). If more than one type of base contributed to the signal, this two level switching would be hard to account for. Finally, experiments with deoxycytidinemonophosphate and 5-methylcytidine monophosphate show that there are clear differences in the size of the tunneling signals. It appears that hydration is playing a role here, because no difference was observed between deoxycytidine and 5-methyl-deoxycytidine measured in organic solvent.19 Methylation has long been known to affect the melting temperature of double helical DNA, reflecting its effects on the bonding between helices in water.47
Recognition tunneling is not limited to nucleic acid targets (DNA and RNA). Half of the amino acid residues have hydrogen-bonding sites in addition to those on the peptide backbone itself. Five (asn, glu, gln, his, arg) have at least two or more sites while another five (lys, ser, thr, try, trp) have one other site available. This raises the possibility of identifying specific peptide sequences in protein biomarkers using Recognition Tunneling, thus avoiding immunoassays that require antibodies for recognition and detection.
We acknowledge financial support from the DNA sequencing technology program of the National Institute of Human Genome Research (HG004378).
Molecule monolayers were formed by adsorption of thiolated molecules onto flame-annealed Au(111) substrates.48 STM imaging showed that 5’-thiolated-nucleosides normally form disordered monolayers. The thickness of the adlayers was measured by ellipsometry and found to be consistent with a monolayer. Fourier transform infrared spectroscopy (FTIR) was used to demonstrate the presence of the desired compounds and confirm an upright orientation on the gold substrates. The formation of appropriately hydrogen-bonded complexes on these surfaces was confirmed with surface plasmon resonance studies of the binding of hydrogen-bonded targets, and also with studies of the adhesion of functionalized atomic force microscope (AFM) probes.
STM tips were made from 0.25 mm diameter gold wire (purity 99.999% from Alfa Aesar) by either AC (30V, 4.2 KHz) or DC (2V) electrochemical etching in mixed solutions of concentrated HCl and ethanol (50:50, v/v). The gold tips were immersed in Piranha Solution (3:1 H2SO4:H2O2 - 30% by volume – use caution as this solution is extremely reactive and will explode on contact with organics for 30 s, rinsed with DI water and dried in an N2 stream prior to use or modification. A typical procedure for functionalization is as follows: tips were immersed in a solution of the probe molecules for periods ranging from 2 hours to overnight, then rinsed with clean solvent repeatedly and blown dry in an N2 stream. Tips for measurements in aqueous solution were coated with wax46 or high density polyethylene.49 High density polyethylene coatings are required when the tip functionalization procedure uses a solvent that dissolves wax insulation. We know of no analytical method sensitive enough to test for successful functionalization of the probes, so tunneling measurements are first carried out with unfunctionalized probes to provide a bench mark for recognizing unsuccessful funtionalization. We currently succeed in functionalizing the majority of the probes we prepare.
Measurements were carried out on a PicoSPM (Agilent, Chandler). The STM tip movement was controlled by custom labview programs and the data were recorded by a Yokogawa digital oscilloscope. 1,2,4-trichlorobenzene (TCB) and aqueous buffered electrolytes solution (0.1x PBS) were both used as working solutions. We performed i-z (STM break-junction) and i-t experiments (STM fixed junction) measurements. For i-z measurements, tunneling current (i) at fixed bias was recorded versus distance (z) while withdrawing tip from the surface in z direction.10
For i-t measurement, tunneling current (i) at fixed bias was recorded versus time (t) while holding tip at a fixed distance above the surface. We developed two methods for these “fixed gap” measurements: 1. Open loop without servo control.15 2. Closed loop with weak servo control.19 For the open loop method, the junction is repeatedly reset by making ever smaller steps in the PZT voltage to avoid PZT creep. The junctions are found to be stable (no detectable shift in the background current) over periods of several seconds, perhaps as a result of the bonding process discussed above.
For the “closed loop” method, the STM servo gains were set low so that only spikes of long duration were affected by the action of the current-control servo. We characterized the response time of the servo by recording thermal noise spectra with, and without servo control.19 The closed-loop method does cause some distortion in signals from bound states of long duration, but it is much easier to implement.