Binding affinity is measured by dissociation constant, Kd, which uses concentration as units. The universal concentration units facilitate direct comparison of affinities for different binding events. However, Kd is a thermodynamic parameter, which lacks kinetic information of a binding event. In addition, Kd does not reveal the mechanical property of the binding, which emerges as a critical element for many physiologically significant processes such as DNA replication, RNA transcription, and protein translation. Here we propose a new parameter, mechanical affinity, to delineate kinetic and mechanical features of a binding event. The mechanical affinity is equivalent to the work required to dissemble the chemical binding between a ligand and a receptor. During this process, it must cover dissipated heat that originates from the relative movement between a ligand and a receptor. Because dissipated heat varies with unfolding direction or rate of mechanical perturbation, the mechanical affinity is a function of these two variables. Screening of chemicals using rupture force of a ligand-receptor complex or mechanical affinity is discussed at the end of this review. The interrogation on the mechanical interaction between a ligand and a receptor provides a new perspective not available in conventional thermodynamic evaluation of binding processes.
Binding between two components represents a universal yet pivotal step for numerous processes. In the human body, binding between antibody and antigen initiates immune responses. For cells, binding between a receptor and a ligand is the first step to trigger signal transduction that involves a multitude of biochemical reactions, most of which occur after an enzyme is bound with a substrate. In materials synthesis, binding of reactants to catalysts produces organic and inorganic compounds. Self-assembly is a general strategy by which many materials are prepared. Association or binding between materials with similar properties is a driving force in the self-assembly. In biosensing, binding between an analyte and a receptor constitutes the recognition stage that generates signals to be amplified in subsequent steps.
Dissociation constant, Kd, measures chemical binding affinity. It depicts the concentration at which half of the binding sites are occupied by ligands. Binding with higher affinities suggests that fewer ligands (lower concentration) are required to occupy the binding sites, which yields a smaller Kd value. The clear-cut physical significance with universal concentration units in Kd renders the comparison between different binding events obvious and intuitive. In addition, measurement of Kd is straightforward. A series of ligand concentrations is all that is required to obtain a Kd value.
Binding is a dynamic process in which association and dissociation proceed simultaneously. Although Kd provides a measurement of binding strength, it does not reveal dynamic nature in a binding, which is better described by the association (kon) and the dissociation (koff) rate constants. Interpretation of these rate constants is not straightforward as different units may exist. It is even more problematic to compare kon or koff between the binding events that have different reaction orders. To collect dynamic information of binding events, it is desirable to adopt a universal variable, similar to the Kd, that can compare different binding processes.
In this review, we propose such a variable, mechanical affinity, to represent both energetic and dynamic aspects of binding events that are essential for many vital biological processes such as enzymatic actions of motor proteins and rolling and adhesion of leukocytes during inflammation. First, we will discuss dissociation constant, Kd, that reflects the thermodynamic nature of a binding process. Conventional methods to measure this variable will be introduced next. In a subsequent section, we will elaborate the concept of mechanical stability and mechanical affinity. These two variables will be compared to the chemical affinity that has been commonly used to define the strength of binding. In the final section of this review, we propose screening methods that exploit mechanical property as a new metrics to evaluate binding processes.
Kd as a common variable to measure chemical affinity
For a simple two-component reaction, A+B→C, the reaction equilibrium constant for association (Ka) is expressed as Ka=[C]eq/([A]eq[B]eq), where [A]eq, [B]eq, and [C]eq represent the concentrations of specific species at equilibrium. The dissociation constant, Kd, is a reciprocal of Ka. As a thermodynamic variable, Kd is related to free energy change of binding, ΔGbinding (or ΔGaffinity), by the expression lnKd=-ΔGbinding/RT, where T is absolute temperature and R is the molar gas constant.
Several methods exist for Kd determination. The most widely used method is based on the Langmuir isotherm approach (Langmuir 1918). In this method, as shown in Figure 1, the fraction of ligand-bound receptors increases with increasing ligand concentration until a plateau is reached, indicating that binding sites are saturated. The Langmuir model [Eq. (1)] (Langmuir 1918) is based on the assumption that all binding sites are equivalent and that the binding ability of molecules is not affected by the occupancy of nearby sites.
where C is the concentration of adsorbate (or ligand), Γ is the amount of adsorbed ligand, and Γmax is the maximum amount of adsorbate on a surface.
Different techniques have been employed to measure Kd of the interaction between a ligand and a receptor (Connors 1987). Among label-free approaches, UV-visible absorption spectroscopy is the simplest and commonly used technique (Xia et al. 2010). Isothermal titration calorimetry not only can measure Kd of the binding but also reveal thermodynamic parameters associated with the binding (Jelesarov and Bosshard 1999, Belozerova and Levicky 2012). Recently, surface plasmon resonance has become a popular optical technique to probe binding on a metal surface (Homola 2003, Karlsson 2004, Mayer and Hafner 2011). It provides quantitative information on the binding affinity as well as the binding kinetics. All these experiments require increasing concentrations of either a ligand or a receptor (Figure 1). The other challenge in these methods is their low sensitivity. Owing to the high sensitivity and broad dynamic range, fluorescence-based approaches have become useful techniques to monitor binding interactions. Among different fluorescence methods, fluorescence polarization anisotropy (Laurence 1952, Jameson and Sawyer 1995), fluorescence resonance energy transfer (FRET) (De Cian et al. 2007), and laser-induced fluorescence capillary electrophoresis (German et al. 1998, Yangyuoru et al. 2012) have been extensively used. However, labeling with a fluorophore may interfere with the binding. In addition, fluorescence signal from environment presents noise to the measurements. To address this issue, recently Kd has been determined by mechanical force that is required to dissemble a ligand-receptor complex (Koirala et al. 2011a, Nguyen et al. 2011, Paramanathan et al. 2012, Yangyuoru et al. 2012).
Mechanical properties of the binding
Mechanical force and its measurement
Mechanical force is a universal parameter that can be used to characterize a wide range of systems ranging from chemical bonds, intermolecular interactions, macromolecular structures, and many biological and biochemical processes. The basic interaction between individual atoms is described in covalent or noncovalent bonds, whose strength can be characterized with force. Change in the conformation of a macromolecule is often accompanied by a variation of the tension sustained by that molecule. Many vital biological processes can be monitored and characterized in terms of mechanical forces. For example, the biophysical aspects of the attachment and the rolling of leukocytes along the wall of blood capillaries during inflammation involve a complex balance of forces arising from blood-flow-induced hydrodynamic shearing effects and the adhesive forces between the leukocytes and the capillary wall (Konstantopoulos et al. 1998, Simon and Goldsmith 2002, Simon and Green 2005). Biochemical actions of many proteins, such as DNA (Wuite et al. 2000) or RNA polymerase (Yin et al. 1995, Davenport et al. 2000), exonuclease (Perkins et al. 2003), myosin (Finer et al. 1994), and kinesin (Visscher et al. 1999), generate forces up to tens of piconewtons. Moreover, activities of those proteins are often influenced by an externally applied force (Davenport et al. 2000). In ligand-receptor binding, the strength of the interaction between two binding partners can be measured in terms of the mechanical force that is required to destroy the interaction from a specific direction (Ainavarapu et al. 2005, Koirala et al. 2011a, Nguyen et al. 2011, Paramanathan et al. 2012, Yangyuoru et al. 2012). In such a case, force is applied at a constant rate to overcome the energy barrier. The dynamics of the binding process can be retrieved by measuring the unbinding rate at a constant force or by measuring rupture force as a function of the force loading rate.
Despite its importance, there exists a lack of ensemble-average methods to interrogate force in solution. Hydrodynamic forces from a flowing fluid have been used to monitor the mechanical aspect of biologically relevant processes such as adhesion and rolling of white blood cells inside blood capillaries during inflammation (Dong and Lei 2000, Bianchi et al. 2013). However, such methods have difficulties to explore the roles of individual binding units at the molecular level. As a result, only little information has been collected to understand the mechanochemistry (Keller and Bustamante 2000), an emerging field to study the coupling between chemical energy and mechanical energy. Single-molecule techniques, however, start to address this problem. Compared to bulk assays where ensemble information is obtained, single-molecule experiments produce stochastic and discrete signals that can be statistically analyzed to reveal the property of subgroups in a population and the energetics of a particular reaction trajectory in a transition process (Koirala et al. 2013, Yu and Mao 2013). Among many tools to manipulate single molecules, force-based techniques, such as atomic force microscopy (AFM) (Binnig et al. 1986, Rief et al. 1997), optical tweezers (Ashkin et al. 1986, Smith et al. 1996, Visscher et al. 1999), microneedles (Kishino and Yanagida 1988), and magnetic tweezers (Smith et al. 1992, Strick et al. 1996), provide a unique capability to investigate the force generated by or applied to individual molecules. Details of each method are reviewed in recent publications (Bustamante et al. 2000, Neuman and Nagy 2008, De Vlaminck and Dekker 2012). Here, we present a brief introduction of these methods.
AFM is commercially available. It uses a microfabricated cantilever to attach a biomolecule through one of its ends (Figure 2A). The other end of the biomolecule is immobilized onto a substrate surface. The displacement between the cantilever and the surface is often controlled by a piezoelectric device. Force signal, which ranges from 10-11 to 10-7 N, in the biomolecule is measured by the deflection of the cantilever. This method has nanometer spatial resolution and piconewton force resolution. The stiffness of a cantilever ranges from 0.001 to 100 Nm-1. It is noteworthy that the stiffer the cantilever, the lower the force sensitivity. As the dynamic range of the force measurement is broad, this method has been widely used to investigate the strength of intramolecular and intermolecular interactions in biomolecules.
Microneedles are usually 50–500 μm long and 0.1–1 μm in diameter. Because of the lower stiffness (10-6 to 1 Nm-1), this method has an advantage over AFM to study delicate systems. Setup is similar to AFM in which a biomolecule is sandwiched between the tip of a microneedle and a substrate surface (Figure 2B). Force signal can be measured by observing the displacement of a bendable microneedle via imaging the microneedle itself or by using a chemically etched optical fiber that projects a reflected light from its tip onto a photodiode. This method can measure force from 10-12 to 10-10 N with a minimum distance of 10-9 m. The instrument is not commercially available.
Magnetic tweezers use permanent or electromagnets to trap magnetic microparticles (Figure 2C). Biomolecules can be tethered between magnetic particles and a nonmagnetic surface. This method can measure force from 10-14 to 10-11 N with a minimum distance of few nanometers. Because it can trap multiple particles, throughput can be higher compared to other force-based single-molecule approaches. It has a unique capability to rotate magnetic particles, thereby allowing torque measurements of torsionally constrained biomolecules. Recently, various types of magnetic tweezers have been developed (Yan et al. 2004, Lipfert et al. 2011), however, the major disadvantage associated with these methods lies in the fact that force is not measured directly (Bustamante et al. 2000, Neuman and Nagy 2008, De Vlaminck and Dekker 2012).
Optical tweezers can trap dielectric microparticles at a laser focus (Figure 2D). Trap stiffness of a laser trap (10-10 to 10-3 Nm-1) is much smaller than that of AFM cantilevers, which allows a better force resolution in the range of 10-13 to 10-10 N. It has a spatial resolution of 10-9–10-10 m, which allows resolving single base pairs in DNA. The force exerted on a microparticle can be flexibly controlled by parameters such as laser power, particle size, and difference in the refractive index between the particle and the trapping medium. In single-beam optical tweezers, a biomolecule is usually tethered between a trapped particle at the laser focus and a surface of a substrate or a particle held at the tip of a micropipette by suction (Kellermayer et al. 1997, Wang et al. 1997, Cecconi et al. 2005, Moffitt et al. 2008). In dual-beam optical tweezers, the molecule is often tethered between two trapped particles at two separate laser foci (Moffitt et al. 2008, Woodside et al. 2008, Yu et al. 2009). The major disadvantage of this method is laser-induced damage of biomolecules.
Using these single-molecule tools, different approaches can be adopted to investigate the mechanical property of a biomolecule and its interaction with a ligand. In one setup, two binding partners can be immobilized to two surfaces separately. After binding between these two partners is accomplished by bringing two surfaces together, moving one surface away from the other increases the mechanical tension. The increased tension eventually destroys the binding interaction. However, the major challenges in this approach are the precise measurement of the distance change during the unfolding and the reversibility of the process. In another approach, an intramolecularly folded biomolecule is first tethered between two surfaces using affinity or covalent linkages (Figure 2). As one surface is moved away from another by a computer-controlled motor, the tension in the biomolecule increases until unfolding occurs, which is manifested by a sudden change in the tension or the end-to-end distance in a force-extension (F-X) curve recorded for this process. During the opposite process in which the end-to-end distance becomes shorter, the biomolecule refolds. The same processes can be performed in the presence of the ligand to evaluate ligand-receptor interactions.
Mechanical unfolding or refolding processes present a completely different mechanism compared to the temperature- or chemical denaturant-mediated unfolding or refolding. Whereas the force-induced processes have localized effects on a molecule, the temperature- or chemical-assisted events are global in nature. The localized force effect suggests that the direction of unfolding is rather critical for the process. However, it is expected that unfolding/refolding should be an ensemble average of all unfolding directions that are presented in the temperature- or the chemical-assisted transition processes owing to their global nature.
From energy perspective, force (F) reduces the free energy of a system by a factor of F×x, where x is the unfolding distance of a biomolecule. As shown in Figure 3, such a reduction brings down the energy barrier of a two-state folding/unfolding system by F×x†, where x† is the transition distance, or the distance between the folded structure and the transition state. With application of a specific force (Frupture), the energy barrier is reduced to a level within reach of the thermal energy of the environment (kBT, where kB is the Boltzmann constant and T is absolute temperature). This leads to unfolding of the structure. Assuming x† does not change during the application of force, the magnitude of Frupture is inversely dependent on the transition distance with a larger distance corresponding to a smaller Frupture value (Evans 2001). The relationship between Frupture and x† can be understood by the elasticity of a molecule. A soft and flexible molecule requires a long distance to be unfolded (larger x†). The force (Frupture) along this long distance is expected to be small. Moreover, a rigid but fragile molecule only needs a little stretch in its structure (smaller x†) to induce unfolding. However, the force leading to this little perturbation (Frupture) is expected to be significantly higher than the previous case.
Because Frupture alone does not fully account for the mechanical perturbation defined by F×x†, there is a caveat to use this parameter to compare the mechanical stability of folded molecules. Comparison of x† among DNA secondary species, such as hairpins (Woodside et al. 2006), DNA G-quadruplexes (Yu et al. 2009, Yu et al. 2012a,b, Koirala et al. 2011a), DNA i-motifs (Dhakal et al. 2010), and their intermediates (Dhakal et al. 2012, Koirala et al. 2012), has suggested that x† varies with different molecules. Within the same molecule, x† is anisotropic among different unfolding trajectories (Yu et al. 2012a). Therefore, mechanical aspects of unfolding can be described accurately only when both Frupture and x† are taken into account. The other caution to use Frupture for the comparison is that magnitude of Frupture is dependent on the force loading rate. The dependence is rather pronounced for processes with slow unfolding rates (Evans and Ritchie 1997, Strunz et al. 1999). In such a case, Frupture for reversible folding and unfolding processes occurs in a timescale often not reachable in an experiment. Although Frupture is straightforward to understand and has universal units in newtons, the above arguments have made Frupture a not-so-ideal parameter to compare different binding processes.
We propose that unfolding work, Wreceptor-unfold= Frupture×Δx, where Δx is the difference in the end-to-end distance between a folded and an unfolded structure along a specific unfolding trajectory (see shaded region in the F-X curve in Figure 4, left panel), is a suitable parameter to account for mechanical properties of an unfolding event. With the single-molecule techniques described above (see Figure 2), Frupture and Δx can be directly measured from the F-X curves (see Figure 4, left panel, for example). Frupture and Δx provide the information on the mechanical stability and the size of the folded biomolecule, respectively.
Wreceptor-unfold has universal energy units, J/mol. Therefore, it can be used among different systems to facilitate direct comparison. This parameter not only considers the work to disrupt the folded structure per se, it also accounts for the work to extend the folded structure from 0 pN to Frupture and the work to stretch the unfolded structure from 0 pN to Frupture (Liphardt et al. 2001). These works represent an actual mechanical stability of a folded structure, which must resist all external perturbations. Because Frupture is a kinetic parameter that is dependent on the force loading rate (see above), Wreceptor-unfold is expected to have a similar feature. In fact, it has been shown that Frupture and Wreceptor-unfold are dependent on the force loading rate in single molecular force spectroscopy experiments (Evans and Ritchie 1997, Rief et al. 1997). Wreceptor-unfold is anisotropic because different unfolding trajectories often suggest different interactions between a receptor and surrounding solvent molecules. This property has been confirmed in the mechanical unfolding of a human telomeric G-quadruplex from different trajectories, which is made available by a click-chemistry modification of specific nucleotides in human telomeric DNA sequences (Yu et al. 2012b, Dhakal et al. 2013).
With a ligand bound to a molecule, the folding and unfolding energetic landscape of the molecule changes. Due to the stabilization effect of the ligand, the free energy of the molecule bound with ligand decreases, whereas that of the unfolded state with a dissociated ligand does not vary significantly (Figure 3). Assuming the energy for a transition barrier does not vary significantly in the presence of a ligand, the transition barrier to unfold the ligand-bound receptor increases.
On the other hand, the transition distance (x†) of a ligand-bound molecule is likely different from that of a free, folded molecule. Retreatment of our recent data on an ATP-bound DNA aptamer (Yangyuoru et al. 2012) shows that the transition distance becomes shorter with respect to that of a free DNA aptamer (Figure 5). Because rupture force (Frupture) is related to the magnitude of x† by Frupture×x† (see Figure 3), it is possible that Frupture shows a marginal increase or even a decrease when x† increases significantly. Such a variation in Frupture may hold even though energy barrier increases as a result of ligand binding. Therefore, Frupture may not accurately represent the mechanical property of a ligand-bound receptor, which is similar to the case for a folded molecule without a bound ligand as discussed above.
Instead, we propose that the work to unfold a ligand-receptor complex, Wcomplex-unfold, can better reflect the mechanical property of the complex. This work can be dissected into two components,
where Wreceptor-unfold is the work to unfold a free, folded receptor as defined above, and Waffinity corresponds to the work required to dissemble the ligand-receptor complex, which is equal to the mechanical affinity. Measurement of mechanical affinity is straightforward from Eq. (2): It is equivalent to the difference between the mechanical work to unfold a free receptor and that to dissemble a receptor-ligand complex, Waffinity=Wcomplex-unfold-Wreceptor-unfold.
Similar to Wreceptor-unfold, Wcomplex-unfold is anisotropic and dependent on the force loading rate. This is because energy required to strip off a ligand from a receptor is likely different among mechanical perturbations applied from different directions. Figure 6 shows the disassembly of a complex in which a ligand, such as histone, is wrapped by a DNA (a receptor). In the first case (Figure 6, top panel), the disruption of the ligand-receptor complex is achieved by grabbing a DNA strand at two positions in the middle of the DNA strand, thereby unwinding the rest of the strand through a ‘sliding’ action (the ‘slide open’ or ‘unwind open’ event). In the second scenario (Figure 6, bottom panel), the ligand-bound DNA is pulled by grabbing the two termini of the DNA strand. In this geometry, the physical contact between a DNA and a ligand becomes tighter until the ligand is popped out from the top or bottom of the tightened DNA strand (the ‘pop open’ process) (Brower-Toland et al. 2002). From a perspective of mechanics, the former process generates less friction between the DNA and the ligand compared to the latter. Such a difference will be reflected in the unfolding work, which must overcome all friction before a ligand is released from a receptor. Because friction originates from a relative movement between a ligand and a receptor during unfolding, the dissipated heat accompanied by frictional is a function of force loading rate: The faster the loading rate, the more the frictional heat. As dissipated heat is an inherent component in Waffinity, this example clearly shows that Waffinity is dynamic and has a localized feature.
Relationship between mechanical affinity and chemical affinity
Chemical affinity between a ligand and a receptor is attributed to intermolecular forces. From free-energy perspective, it is equivalent to the change in free energy of the binding between a ligand and a receptor (ΔGaffinity=ΔGbinding; it is a negative value). Using a Hess-like cycle (Koirala et al. 2011a), ΔGaffinity can be calculated as the difference between the change in free energy of unfolding a free receptor (ΔGreceptor-unfold, a smaller positive value) and that of unfolding a ligand-receptor complex (ΔGcomplex-unfold, a bigger positive value),
From mechanical energy perspective, the work to dissemble a ligand-receptor complex [Waffinity; see Eq. (2)] must overcome the intermolecular forces that constitute the chemical affinity (-ΔGaffinity). In addition, Waffinity should also cover the dissipated heat (Wfriction) originated from the relative motion between a ligand and a receptor during the disassembly process. Thus, the mechanical affinity can be expressed by the following equation,
which represents the first law of thermodynamics for ligand-receptor interactions from a mechanical perspective. For a reversible process in which the force loading rate is infinitely slow, the chemical affinity (ΔGaffinity) and mechanical affinity (Waffinity) are equal (ΔGaffinity=-Waffinity) because dissipated heat (Wfriction) is negligible under such a situation.
It is noteworthy that ΔGaffinity, as a state function like ΔGunfold, should not vary with unfolding trajectories or loading rates. The isotropic nature of ΔGunfold has been confirmed in the unfolding of a DNA G-quadruplex from different directions (Yu et al. 2012b, Dhakal et al. 2013). In another set of experiments, by varying the force loading rate during unfolding or refolding processes, ΔGunfold can be retrieved from the work at which the unfolding and the refolding share the same probability (Crooks fluctuation theorem; Crooks 1999, Collin et al. 2005). Extending this trend, ΔGunfold is equivalent to the work when work histograms of the unfolding and the refolding superimpose, which occurs when folding and unfolding become reversible under infinitely slow force loading rate.
However, such a slow rate cannot be achieved in actual experiments within a reasonable time frame. In addition, mechanical unfolding of a bound complex is often not reversible as refolding of a receptor must occur prior to the binding of the ligand, whereas unfolding of a ligand-bound complex can happen cooperatively (simultaneously). Furthermore, refolding may not always be measured. For the slow refolding process, it occurs at a low-force region. The compromised signal-to-noise ratio in this region prevents a reliable measurement of the refolding force (Frefold) and refolding distance (Δx).
In practice, ΔGunfold can be retrieved from the mechanical unfolding work by using Jarzynski equality theorem (Jarzynski 1997, Liphardt et al. 2002) for nonequilibrium systems. To do this, the hysteresis area between the stretching and relaxing F-X curves (see Figure 4), which is equivalent to the work done during the unfolding process, is calculated from hundreds of F-X curves. At the single-molecule level at which unfolding work is comparable to thermal energy (kBT), there exists a small probability that contribution from thermal bath leads to an unfolding work smaller than ΔGunfold. By increasing weighing factors of these smaller works with respect to other works, ΔGunfold can be retrieved by the following equation (Jarzynski 1997, Liphardt et al. 2002),
where N is the total number of F-X curves and W is the nonequilibrium work done to unfold a structure. In principle, chemical affinity, ΔGaffinity, can be retrieved from mechanical affinity, Waffinity, using this Jarzynski approach.
Ligand screening based on mechanical properties
Currently, almost all drug screening is based on thermodynamic binding (or chemical affinity) between target receptors and small chemicals (ligands). Although this approach is highly effective to identify drug leads, it does not have the capability to reveal binding kinetics and completely ignores mechanical aspects of a binding.
Kinetics is a rather important aspect in receptor-ligand interactions. Many times, a transient binding between a ligand and a receptor is sufficient to lead to subsequent effects in cells. However, tight binding is not always desirable as it may cause unwanted long-term side effects due to the slow dissociation of a ligand-receptor complex. The mechanical property of a binding between a receptor and a ligand is a completely different concept in that it can profile kinetics as well as mechanics of the binding. Such a kinetic property is inherent for many critical processes orchestrated by motor enzymes such as RNA and DNA polymerases, ribosomes, and DNA/RNA helicases. It also plays a vital role in many biological processes such as attachment and rolling of the leukocytes in blood tubes during inflammation. For all these processes, mechanical properties of a receptor-ligand interaction are physiologically as important as thermodynamic features. In protein translation for example, there is a tug-of-war between the mechanical stability of RNA secondary structures and the stall force of ribosome, which represents a maximum load force as ribosome moves along the RNA template. Any potential secondary structures in mRNA must be resolved before a ribosome can continuously synthesize proteins. When a ligand binds to an RNA secondary structure, the mechanical affinity adds to the mechanical stability of a free RNA structure [Eq. (2)], increasing the mechanical work a ribosome must spend to clear the roadblocks.
However, depending on the change in the transition distance (x†) upon ligand binding, the Frupture of the RNA-ligand complex may not exhibit a significant difference compared to those of the RNA structure without ligand. If this Frupture is within reach of a ribosome, i.e., Frupture<stall force of the ribosome, the RNA structure bound with ligand may not be effective to block the translation. Therefore, for the purpose of blocking protein translation, screening of ligands based on chemical or mechanical affinity may not always produce desirable effects. Instead, the screening based on the Frupture of ligand-bound RNA secondary structures can be more effective.
One major advantage of using Frupture, instead of chemical affinity, in the drug screening lies in the fact that the method can identify ligands with weak chemical affinities or fast koff rates. Current chemical-affinity-based screening prefers ligands with strong thermodynamic binding capability. However, experiments have shown that ligands with weak binding can have important physiological effects (Maerkl and Quake 2007). Moreover, in many surface-based heterogeneous screening, bindings with fast koff rate are often overlooked (Maerkl and Quake 2007), even though they have high chemical affinities. The fast koff issue can be addressed as mechanical screening has submicrosecond temporal resolution, and is homogeneous due to its single molecular nature. The weak binding can be evaluated by measuring the mechanical affinity, which can be readily obtained from the area underneath the rupture events as shown in the F-X curves (Figure 4). In fact, binding as weak as Kd=42 μm has been successfully measured recently (Koirala et al. 2011a). In the case where binding affinity is within a few kBT, spontaneous rupture of an interaction becomes significant. To stabilize the binding, often times many such weak interactions are present. Whereas it is difficult for ensemble-based approaches to characterize these multiple events, single-molecule methods described here offer a unique opportunity to probe these interactions.
However, caution must be given to consider the change in x† in the mechanical-affinity-based evaluation. As discussed above, if a ligand increases the transition distance upon binding, Frupture may not show an obvious increase. In such a case, force loading rate can be varied to search for a particular x† that can still yield a detectable change in Frupture (Evans 2001). For binding processes that lead to the evolution of different conformations, comparison of both Fruptrue and Δx between different populations could be useful to probe these successive events.
Another advantage of mechanical-property-based evaluation is that it can directly assess binding from mechanical perturbations applied from different directions. For physiologically relevant nucleic acid structures, mechanical perturbations from motor proteins mostly come from the 5′ to 3′ direction. However, when tertiary interaction exists between a protein and a nucleic acid structure, the perturbation is no longer necessarily from this direction. Using a click chemistry modification approach, the desired mechanical unfolding direction can be generated for biomacromolecules (Yu et al. 2012b, Dhakal et al. 2013). Such a direction-specific evaluation helps to evaluate the effect of a ligand on a receptor from physiologically important directions that are often blurred in chemical-affinity-based approaches due to their ensemble average nature.
A drawback for ligand screening using mechanical properties is the throughput of an assay. Optical tweezers and AFM can only evaluate one molecule a time. Although throughput can be increased by using microfluidic channels in which each channel can contain a different ligand (Koirala et al. 2011b), the point-by-point data collection still limits the throughput of the method. In AFM, it is possible to perform multiplexing screening by incorporation of multiple cantilevers (Fritz et al. 2000). Magnetic tweezes may provide an alternative approach with increased throughput (Danilowicz et al. 2009). However, the highest force it can apply may not be enough to dissemble a receptor-ligand complex.
Recently, Severin et al. (2011) employed an elegant strategy to evaluate DNA-binding ligands by placing DNA molecules between two planes followed by unfolding of the DNA-ligand complex after moving apart these two surfaces (Figure 7). By observing the fluorescence signal left at one of the surfaces, the effect of ligand binding can be quantified. This method does not require significant investment in instrumentation while providing a high-throughput capacity as multiple ligands can be evaluated in addressable wells on a chip.
Conclusions and outlook
With the invention of force-based single-molecule techniques such as optical tweezers and AFM, now it is possible to evaluate mechanical properties of the binding between a ligand and a receptor, which are largely overlooked in thermodynamic-based methods. We propose here that by using a universal parameter, mechanical affinity, binding can be conveniently evaluated from the unfolding work of a bound complex. The mechanical affinity is anisotropic and reflects the dynamic nature of the interaction between a ligand and a receptor. By exploiting mechanical properties such as rupture force, transition distance, and mechanical affinity, new ligand-screening methods can be designed to evaluate ligand-receptor complex from a specific unfolding direction. Compared to current screening approaches that are based on isotropic thermodynamic properties, these newly proposed methods are physiologically more relevant, especially in biological processes that involve motor proteins.
About the authors
Deepak Koirala was born in 1982 and is a citizen of Nepal. From 1999 to 2005, he studied Chemistry and received his BS and MS degrees from Tribhuvan University, Kathmandu, Nepal. After 2 years spent in teaching, he joined the Department of Chemistry and Biochemistry at Kent State University in August 2008 as a graduate student. Currently, he is working in Dr. Hanbin Mao’s research laboratory pursuing a doctorate degree. His research involves studies on mechanochemistry and transition dynamics of DNA secondary structures such as G-quadruplexes and their binding interactions with small-molecule ligands at single-molecule level. He is also working on the development of highly sensitive and high-throughput nanomechanical biosensing platforms.
Philip Yangyuoru is a native of Ghana. He received his Bachelor’s degree in Chemistry from the University of Ghana in 2002 and MS in Chemistry from Illinois State University under Dr. C.F. Shaw III in 2007, studying the interaction of gold drugs and glutathione by electrospray ionization mass spectrometry (ESI-MS). Currently, he is a doctoral student at the Department of Chemistry and Biochemistry, Kent State University, and working with Dr. Mao’s research group. His research focuses on investigating the mechanical stability, kinetics, and thermodynamics of nucleic acid structures and their interactions with small-molecule ligands in the bulk and at the single-molecule levels using capillary electrophoresis and laser tweezers, respectively. He is also interested in the use of DNA aptamers to develop force-based biosensors.
Hanbin Mao completed his BS degree at the West China University of Medical Sciences (now Sichuan University, China) in 1995 and received his PhD degree in Analytical Chemistry from Texas A&M University in 2003 focusing on lab-on-a-chip systems. He spent the following 2 years as a postdoctoral associate at the University of California at Berkeley and the Laurence Berkeley National Lab, where he studied biophysical aspects of biomolecular structures at the single molecular level. In 2005, he started his professional academic career at the Department of Chemistry and Biochemistry, Kent State University. His main research interests are the structure and transition dynamics of non-B DNA species, as well as the development of highly sensitive biosensors.
We thank the National Science Foundation (NSF) CHE-1026532 for financial support.
Ainavarapu, S. R. K.; Li, L.; Badilla, C. L.; Fernandez, J. M. Ligand binding modulates the mechanical stability of dihydrofolate reductase. Biophys. J. 2005, 89, 3337–3344.Search in Google Scholar
Ashkin, A.; Dziedzic, J. M.; Bjorkholm, J. E.; Chu, S. Observation of a single-beam gradient force optical trap for dielectric particles. Optics Lett. 1986, 11, 288.Search in Google Scholar
Belozerova, I.; Levicky, R. Melting thermodynamics of reversible DNA/ligand complexes at interfaces. J. Am. Chem. Soc. 2012, 134, 18667–18676.Search in Google Scholar
Bianchi, E.; Molteni, R.; Pardi, R.; Dubini, G. Microfluidics for in vitro biomimetic shear stress-dependent leukocyte adhesion assays. J. Biomech. 2013, 46, 276–283.Search in Google Scholar
Brower-Toland, B. D.; Smith, C. L.; Yeh, R. C.; Lis, J. T.; Peterson, C. L.; Wang, M. D. Mechanical disruption of individual nucleosomes reveals a reversible multistage release of DNA. Proc. Natl. Acad. Sci. USA 2002, 99, 1960–1965.Search in Google Scholar
Bustamante, C.; Macosko, J. C.; Wuite, G. J. L. Grabbing the cat by the tail: manipulating molecules one by one. Nat. Rev. Mol. Cell Biol. 2000, 1, 130–136.Search in Google Scholar
Cecconi, C.; Shank, E. A.; Bustamante, C.; Marqusee, S. Direct observation of the three-state folding of a single protein molecule. Science 2005, 309, 2057–2060.Search in Google Scholar
Collin, D.; Ritort, F.; Jarzynski, C.; Smith, S. B.; Tinoco, I. J.; Bustamante, C. Verification of the Crooks fluctuation theorem and recovery of RNA folding free energies. Nature 2005, 437, 231–234.Search in Google Scholar
Connors, K. A. Binding Constants. John Wiley & Sons: New York, 1987.Search in Google Scholar
Crooks, G. E. Entropy production fluctuation theorem and the nonequilibrium work relation for free-energy differences. Phys. Rev. E 1999, 60, 2721–2726.10.1103/PhysRevE.60.2721Search in Google Scholar PubMed
Danilowicz, C.; Lee, C. H.; Kim, K.; Hatch, K.; Coljee, V. W.; Kleckner, N.; Prentiss, M. Single molecule detection of direct, homologous, DNA/DNA pairing. Proc. Natl. Acad. Sci. USA 2009, 106, 19824–19829.Search in Google Scholar
Davenport, R. J.; Wuite, G. J. L.; Landick, R.; Bustamante, C. Single-molecule study of transcriptional pausing and arrest by E. coli RNA polymerase. Science 2000, 287, 2497–2500.Search in Google Scholar
De Cian, A.; Guittat, L.; Kaiser, M.; Sacca, B.; Amrane, S.; Bourdoncle, A.; Alberti, P.; Teulade-Fichou, M. P.; Lacroix, L.; Mergny, J. L. Fluorescence-based melting assays for studying quadruplex ligands. Methods 2007, 42, 183–195.Search in Google Scholar
De Vlaminck, I.; Dekker, C. Recent advances in magnetic tweezers. Annu. Rev. Biophys. 2012, 41, 453–72.Search in Google Scholar
Dhakal, S.; Schonhoft, J. D.; Koirala, D.; Yu, Z.; Basu, S.; Mao, H. Coexistence of an ILPR i-motif and a partially folded structure with comparable mechanical stability revealed at the single-molecule level. J. Am. Chem. Soc. 2010, 132, 8991–8997.Search in Google Scholar
Dhakal, S.; Lafontaine, J. L.; Yu, Z.; Koirala, D.; Mao, H. Intramolecular folding in human ILPR fragment with three C-rich repeats. PLoS One 2012, 7, e39271.10.1371/journal.pone.0039271Search in Google Scholar PubMed PubMed Central
Dhakal, S.; Cui, Y.; Koirala, D.; Ghimire, C.; Kushwaha, S.; Yu, Z.; Yangyuoru, P. M.; Mao, H. Structural and mechanical properties of individual human telomeric G-quadruplexes in molecularly crowded solutions. Nucleic Acids Res. 2013, 41, 3915–3923.Search in Google Scholar
Dong, C.; Lei, X. X. Biomechanics of cell rolling: shear flow, cell-surface adhesion, and cell deformability. J. Biomech. 2000, 33, 35–43.Search in Google Scholar
Evans, E. Probing the relation between force-lifetime-and chemistry in single molecular bonds. Annu. Rev. Biophys. Biomol. Struct. 2001, 30, 105–28.Search in Google Scholar
Evans, E.; Ritchie, K. Dynamic strength of molecular adhesion bonds. Biophys. J. 1997, 72, 1541–1555.Search in Google Scholar
Finer, J. T.; Simmons, R. M.; Spudich, J. A. Single myosin molecule mechanics: piconewton forces and nanometer steps. Nature 1994, 368, 113–119.Search in Google Scholar
Fritz, J.; Baller, M. K.; Lang, H. P.; Rothuizen, H.; Vettiger, P.; Meyer, E.; Güntherodt, H.-J.; Gerber, C.; Gimzewski, J. K. Translating biomolecular recognition into nanomechanics. Science 2000, 288, 316–318.Search in Google Scholar
German, I.; Buchanan, D. D.; Kennedy, R. T. Aptamers as ligands in affinity probe capillary electrophoresis. Anal. Chem. 1998, 70, 4540–4545.Search in Google Scholar
Homola, J. Present and future of surface plasmon resonance biosensors. Anal. Bioanal. Chem. 2003, 377, 528–539.Search in Google Scholar
Jameson, D. M.; Sawyer, W. H. Fluorescence anisotropy applied to biomolecular interactions. Methods Enzymol. 1995, 246, 283–300.Search in Google Scholar
Jarzynski, C. Nonequilibrium equality for free energy differences. Phys. Rev. Lett. 1997, 78, 2690–2693.Search in Google Scholar
Jelesarov, I.; Bosshard, H. R. Isothermal titration calorimetry and differential scanning calorimetry as complementary tools to investigate the energetics of biomolecular recognition. J. Mol. Recognit. 1999, 12, 3–18.Search in Google Scholar
Karlsson, R. SPR for molecular interaction analysis: a review of emerging application areas. J. Mol. Recognit. 2004, 17, 151–61.Search in Google Scholar
Keller, D.; Bustamante, C. The mechanochemistry of molecular motors. Biophys. J. 2000, 78, 541–556.Search in Google Scholar
Kellermayer, M. S.; Smith, S. B.; Granzier, H. L.; Bustamante, C. Folding-unfolding transitions in single titin molecules characterized with laser tweezers. Science 1997, 276, 1112–1116.Search in Google Scholar
Koirala, D.; Dhakal, S.; Ashbridge, B.; Sannohe, Y.; Rodriguez, R.; Sugiyama, H.; Balasubramanian, S.; Mao, H. A single-molecule platform for investigation of interactions between G-quadruplexes and small-molecule ligands. Nat. Chem. 2011a, 3, 782–787.Search in Google Scholar
Koirala, D.; Yu, Z.; Dhakal, S.; Mao, H. Detection of single nucleotide polymorphism using tension-dependent stochastic behavior of a single-molecule template. J. Am. Chem. Soc. 2011b, 133, 9988–9991.Search in Google Scholar
Koirala, D.; Mashimo, T.; Sannohe, Y.; Yu, Z.; Mao, H.; Sugiyama, H. Intramolecular folding in three tandem guanine repeats of human telomeric DNA. Chem. Commun. 2012, 48, 2006–2008.Search in Google Scholar
Koirala, D.; Ghimire, C.; Bohrer, C.; Sannohe, Y.; Sugiyama, H.; Mao, H. Long-loop G-quadruplexes are misfolded population minorities with fast transition kinetics in human telomeric sequences. J. Am. Chem. Soc. 2013, 135, 2235–2241.Search in Google Scholar
Konstantopoulos, K.; Kukreti, S.; McIntire, L. V. Biomechanics of cell interactions in shear fields. Adv. Drug Deliv. Rev. 1998, 33, 141–164.Search in Google Scholar
Langmuir, I. The adsorption of gases on plane surfaces of glass, mica and platinum. J. Am. Chem. Soc. 1918, 40, 1361–1403.Search in Google Scholar
Laurence, D. J. R. A study of the adsorption of dyes on bovine serum albumin by the method of polarization of fluorescence. Biochem. J. 1952, 51, 168–180.Search in Google Scholar
Li, P. T. X.; Collin, D.; Smith, S. B.; Bustamante, C.; Tinoco, I. Jr. Probing the mechanical folding kinetics of TAR RNA by hopping, force-jump, and force-ramp methods. Biophys. J. 2006, 90, 250–260.Search in Google Scholar
Lipfert, J.; Wiggin, M.; Kerssemakers, J. W.; Pedaci, F.; Dekker, N. H. Freely orbiting magnetic tweezers to directly monitor changes in the twist of nucleic acids. Nat. Commun. 2011, 2, 439–448.Search in Google Scholar
Liphardt, J.; Onoa, B.; Smith, S. B.; Tinoco, I. Jr.; Bustamante, C. Reversible unfolding of single RNA molecules by mechanical force. Science 2001, 292, 733–737.Search in Google Scholar
Liphardt, J.; Dumont, S.; Smith, S. B.; Tinoco, I. Jr.; Bustamante, C. Equilibrium information from nonequilibrium measurements in an experimental test of Jarzynski’s equality. Science 2002, 296, 1832–1835.Search in Google Scholar
Maerkl, S. J.; Quake, S. R. A systems approach to measuring the binding energy landscapes of transcription factors. Science 2007, 315, 233–237.Search in Google Scholar
Mayer, K. M.; Hafner, J. H. Localized surface plasmon resonance sensors. Chem. Rev. 2011, 111, 3828–3857.Search in Google Scholar
Moffitt, J. R.; Chemla, Y. R.; Smith, S. B.; Bustamante, C. Recent advances in optical tweezers. Annu. Rev. Biochem. 2008, 77, 205–228.10.1146/annurev.biochem.77.043007.090225Search in Google Scholar PubMed
Neuman, K. C.; Nagy, A. Single-molecule force spectroscopy: optical tweezers, magnetic tweezers and atomic force microscopy. Nat. Methods 2008, 5, 491–505.Search in Google Scholar
Nguyen, T.-H.; Steinbock, L. J.; Butt, H.-J.; Helm, M.; Berger, R. D. Measuring single small molecule binding via rupture forces of a split aptamer. J. Am. Chem. Soc. 2011, 133, 2025–2027.Search in Google Scholar
Paramanathan, T.; Vladescu, I.; McCauley, M. J.; Rouzina, I.; Williams, M. C. Force spectroscopy reveals the DNA structural dynamics that govern the slow binding of Actinomycin D. Nucleic Acids Res. 2012, 40, 4925–4932.Search in Google Scholar
Perkins, T. T.; Dalal, R. V.; Mitsis, P. G.; Block, S. M. Sequence-dependent pausing of single lambda exonuclease molecules. Science 2003, 301, 1914–1918.Search in Google Scholar
Rief, M.; Gautel, M.; Oesterhelt, F.; Fernandez, J. M.; Gaub, H. E. Reversible unfolding of individual titin immunoglobulin domains by AFM. Science 1997, 276, 1109–1112.Search in Google Scholar
Severin, P. M. D.; Ho, D.; Gaub, H. E. A high throughput molecular force assay for protein DNA interactions.Lab Chip 2011, 11, 856–862.Search in Google Scholar
Simon, S. I.; Goldsmith, H. L. Leukocyte adhesion dynamics in shear flow. Ann. Biomed. Eng. 2002, 30, 315–332.Search in Google Scholar
Simon, S. I.; Green, C. E. Molecular mechanics and dynamics of leukocyte recruitment during inflammation. Annu. Rev. Biomed. Eng. 2005, 7, 151–185.Search in Google Scholar
Smith, S. B.; Finzi, L.; Bustamante, C. Direct mechanical measurements of the elasticity of single DNA molecules by using magnetic beads. Science 1992, 258, 1122–1126.Search in Google Scholar
Smith, S. B.; Cui, Y. J.; Bustamante, C. Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules. Science 1996, 271, 795–799.Search in Google Scholar
Strick, T. R.; Allemand, J.-F.; Bensimon, D.; Bensimon, A.; Croquette, V. The elasticity of a single supercoiled DNA molecule. Science 1996, 271, 1835–1837.Search in Google Scholar
Strunz, T.; Oroszlan, K.; Schäfer, R.; Güntherodt, H. J. Dynamic force spectroscopy of single DNA molecules. Proc. Natl. Acad. Sci. USA 1999, 96, 11277–11282.Search in Google Scholar
Visscher, K.; Schnitzer, M. J.; Block, S. M. Single kinesin molecules studied with a molecular force clamp. Nature 1999, 400, 184–189.Search in Google Scholar
Wang, M. D.; Yin, H.; Landick, R.; Gelles, J.; Block, S. M. Stretching DNA with optical tweezers. Biophys. J. 1997, 72, 1335–1346.Search in Google Scholar
Woodside, M. T.; Behnke-Parks, W. M.; Larizadeh, K.; Travers, K.; Herschlag, D.; Block, S. M. Nanomechanical measurements of the sequence-dependent folding landscapes of single nucleic acid hairpins. Proc. Natl. Acad. Sci. USA 2006, 103, 6190–6195.Search in Google Scholar
Woodside, M. T.; García-García, C.; Block, S. M. Folding and unfolding single RNA molecules under tension. Curr. Opin. Chem. Biol. 2008, 12, 640–646.Search in Google Scholar
Wuite, G. J. L.; Smith, S. B.; Young, M.; Keller, D.; Bustamante, C. Single-molecule studies of the effect of template tension on T7 DNA polymerase activity. Nature 2000, 404, 103–106.Search in Google Scholar
Xia, F.; Zuo, X.; Yang, R.; Xiao, Y.; Kang, D.; Vallee-Belisle, A.; Gong, X.; Yuen, J. D.; Hsu, B. B.; Heeger, A. J.; Plaxco, K. W. Colorimetric detection of DNA, small molecules, proteins, and ions using unmodified gold nanoparticles and conjugated polyelectrolytes. Proc. Natl. Acad. Sci. USA 2010, 107, 10837–10841.Search in Google Scholar
Yangyuoru, P. M.; Dhakal, S.; Yu, Z.; Koirala, D.; Mwongela, S. M.; Mao, H. Single-molecule measurements of the binding between small molecules and DNA aptamers. Anal. Chem. 2012, 84, 5298–5303.Search in Google Scholar
Yin, H.; Wang, M. D.; Svoboda, K.; Landick, R.; Block, S. M.; Gelles, J. Transcription against an applied force. Science 1995, 270, 1653–1657.Search in Google Scholar
Yu, Z.; Mao, H. Non-B DNA Structures show diverse conformations and complex transition kinetics comparable to RNA or proteins – a perspective from mechanical unfolding and refolding experiments. Chem. Rec. 2013, 13, 102–116.Search in Google Scholar
Yu, Z.; Schonhoft, J. D.; Dhakal, S.; Bajracharya, R.; Hegde, R.; Basu, S.; Mao, H. ILPR G-Quadruplexes formed in seconds demonstrate high mechanical stabilities. J. Am. Chem. Soc. 2009, 131, 1876–1882.Search in Google Scholar
Yu, Z.; Gaerig, V.; Cui, Y.; Kang, H.; Gokhale, V.; Zhao, Y.; Hurley, L. H.; Mao, H. The tertiary DNA structure in the single-stranded hTERT promoter fragment unfolds and refolds by parallel pathways via cooperative or sequential events. J. Am. Chem. Soc. 2012a, 134, 5157–5164.Search in Google Scholar
Yu, Z.; Koirala, D.; Cui, Y.; Easterling, L. F.; Zhao, Y.; Mao, H. Click chemistry assisted single-molecule fingerprinting reveals a 3D biomolecular folding funnel. J. Am. Chem. Soc. 2012b, 134, 12338–12341.Search in Google Scholar
©2013 by Walter de Gruyter Berlin Boston
This article is distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.