Antigen binding forces of individually addressed single-chain Fv antibody molecules

Antibody single-chain Fv fragment (scFv) molecules that are specific for fluorescein have been engineered with a C-terminal cysteine for a directed immobilization on a flat gold surface. Individual scFv molecules can be identified by atomic force microscopy. For selected molecules the antigen binding forces are then determined by using a tip modified with covalently immobilized antigen. An scFv mutant of 12% lower free energy for ligand binding exhibits a statistically significant 20% lower binding force. This strategy of covalent immobilization and measuring well separated single molecules allows the characterization of ligand binding forces in molecular repertoires at the single molecule level and will provide a deeper insight into biorecognition processes.

Structure-based conformational preferences of amino acids

Proteins can be very tolerant to amino acid substitution, even within their core. Understanding the factors responsible for this behavior is of critical importance for protein engineering and design. Mutations in proteins have been quantified in terms of the changes in stability they induce. For example, guest residues in specific secondary structures have been used as probes of conformational preferences of amino acids, yielding propensity scales. Predicting these amino acid propensities would be a good test of any new potential energy functions used to mimic protein stability. We have recently developed a protein design procedure that optimizes whole sequences for a given target conformation based on the knowledge of the template backbone and on a semiempirical potential energy function. This energy function is purely physical, including steric interactions based on a Lennard-Jones potential, electrostatics based on a Coulomb potential, and hydrophobicity in the form of an environment free energy based on accessible surface area and interatomic contact areas. Sequences designed by this procedure for 10 different proteins were analyzed to extract conformational preferences for amino acids. The resulting structure-based propensity scales show significant agreements with experimental propensity scale values, both for α-helices and β-sheets. These results indicate that amino acid conformational preferences are a natural consequence of the potential energy we use. This confirms the accuracy of our potential and indicates that such preferences should not be added as a design criterion.

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78°) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 μM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

Semirational design of active tumor suppressor p53 DNA binding domain with enhanced stability

We have designed a p53 DNA binding domain that has virtually the same binding affinity for the gadd45 promoter as does wild-type protein but is considerably more stable. The design strategy was based on molecular evolution of the protein domain. Naturally occurring amino acid substitutions were identified by comparing the sequences of p53 homologues from 23 species, introducing them into wild-type human p53, and measuring the changes in stability. The most stable substitutions were combined in a multiple mutant. The advantage of this strategy is that, by substituting with naturally occurring residues, the function is likely to be unimpaired. All point mutants bind the consensus DNA sequence. The changes in stability ranged from +1.27 (less stable Q165K) to −1.49 (more stable N239Y) kcal mol−1, respectively. The changes in free energy of unfolding on mutation are additive. Of interest, the two most stable mutants (N239Y and N268D) have been known to act as suppressors and restored the activity of two of the most common tumorigenic mutants. Of the 20 single mutants, 10 are cancer-associated, though their frequency of occurrence is extremely low: A129D, Q165K, Q167E, and D148E are less stable and M133L, V203A and N239Y are more stable whereas the rest are neutral. The quadruple mutant (M133LV203AN239YN268D), which is stabilized by 2.65 kcal mol−1 and Tm raised by 5.6°C is of potential interest for trials in vivo.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 × 106 M−1, two orders of magnitude higher than that for biotin (1 × 104 M−1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 × 104 M−1) was lower than predicted (2.9 × 105 M−1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Peptide interfacial adsorption is kinetically limited by the thermodynamic stability of self association

We present a study of the adsorption of two peptides at the octane–water interface. The first peptide, Lac21, exists in mixed monomer–tetramer equilibrium in bulk solution with an appreciable monomer concentration. The second peptide, Lac28, exists as a tetramer in solution, with minimal exposed hydrophobic surface. A kinetic limitation to interfacial adsorption exists for Lac28 at moderate to high surface coverage that is not observed for Lac21. We estimate the potential energy barrier for Lac28 adsorption to be 42 kJ/mol and show that this is comparable to the expected free energy barrier for tetramer dissociation. This finding suggests that, at moderate to high surface coverage, adsorption is kinetically limited by the availability of interfacially active monomeric “domains” in the subinterfacial region. We also show how the commonly used empirical equation for protein adsorption dynamics can be used to estimate the potential energy barrier for adsorption. Such an approach is shown to be consistent with a formal description of diffusion–adsorption, provided a large potential energy barrier exists. This work demonstrates that the dynamics of interfacial adsorption depend on protein thermodynamic stability, and hence structure, in a quantifiable way.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

Adenosine diphosphate glucose pyrophosphatase: A plastidial phosphodiesterase that prevents starch biosynthesis

A distinct phosphodiesterasic activity (EC 3.1.4) was found in both mono- and dicotyledonous plants that catalyzes the hydrolytic breakdown of ADPglucose (ADPG) to produce equimolar amounts of glucose-1-phosphate and AMP. The enzyme responsible for this activity, referred to as ADPG pyrophosphatase (AGPPase), was purified over 1,100-fold from barley leaves and subjected to biochemical characterization. The calculated Keq′ (modified equilibrium constant) value for the ADPG hydrolytic reaction at pH 7.0 and 25°C is 110, and its standard-state free-energy change value (ΔG′) is −2.9 kcal/mol (1 kcal = 4.18 kJ). Kinetic analyses showed that, although AGPPase can hydrolyze several low-molecular weight phosphodiester bond-containing compounds, ADPG proved to be the best substrate (Km = 0.5 mM). Pi and phosphorylated compounds such as 3-phosphoglycerate, PPi, ATP, ADP, NADP+, and AMP are inhibitors of AGPPase. Subcellular localization studies revealed that AGPPase is localized exclusively in the plastidial compartment of cultured cells of sycamore (Acer pseudoplatanus L.), whereas it occurs both inside and outside the plastid in barley endosperm. In this paper, evidence is presented that shows that AGPPase, whose activity declines concomitantly with the accumulation of starch during development of sink organs, competes with starch synthase (ADPG:1,4-α-d-glucan 4-α-d-glucosyltransferase; EC 2.4.1.21) for ADPG, thus markedly blocking the starch biosynthesis.

Native protein sequences are close to optimal for their structures

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.