Free energy of burying hydrophobic residues in the interface between protein subunits

We have obtained an experimental estimate of the free energy change associated with variations at the interface between protein subunits, a subject that has raised considerable interest since the concept of accessible surface area was introduced by Lee and Richards [Lee, B. & Richards, F. M. (1971) J. Mol. Biol. 55, 379–400]. We determined by analytical ultracentrifugation the dimer–tetramer equilibrium constant of five single and three double mutants of human Hb. One mutation is at the stationary α1β1 interface, and all of the others are at the sliding α1β2 interface where cleavage of the tetramer into dimers and ligand-linked allosteric changes are known to occur. A surprisingly good linear correlation between the change in the free energy of association of the mutants and the change in buried hydrophobic surface area was obtained, after corrections for the energetic cost of losing steric complementarity at the αβ dimer interface. The slope yields an interface stabilization free energy of −15 ± 1.2 cal/mol upon burial of 1 Å2 of hydrophobic surface, in very good agreement with the theoretical estimate given by Eisenberg and McLachlan [Eisenberg, D. & McLachlan, A. D. (1986) Nature (London) 319, 199–203].

Bond orientational order, molecular motion, and free energy of high-density DNA mesophases.

By equilibrating condensed DNA arrays against reservoirs of known osmotic stress and examining them with several structural probes, it has been possible to achieve a detailed thermodynamic and structural characterization of the change between two distinct regions on the liquid-crystalline phase diagram: (i) a higher density hexagonally packed region with long-range bond orientational order in the plane perpendicular to the average molecular direction and (ii) a lower density cholesteric region with fluid-like positional order. X-ray scattering on highly ordered DNA arrays at high density and with the helical axis oriented parallel to the incoming beam showed a sixfold azimuthal modulation of the first-order diffraction peak that reflects the macroscopic bond-orientational order. Transition to the less-dense cholesteric phase through osmotically controlled swelling shows the loss of this bond orientational order, which had been expected from the change in optical birefringence patterns and which is consistent with a rapid onset of molecular positional disorder. This change in order was previously inferred from intermolecular force measurements and is now confirmed by 31P NMR. Controlled reversible swelling and compaction under osmotic stress, spanning a range of densities between approximately 120 mg/ml to approximately 600 mg/ml, allow measurement of the free-energy changes throughout each phase and at the phase transition, essential information for theories of liquid-crystalline states.

Bioenergetic scaling: metabolic design and body-size constraints in mammals.

The cytosolic phosphorylation ratio ([ATP]/[ADP][P(i)]) in the mammalian heart was found to be inversely related to body mass with an exponent of -0.30 (r = 0.999). This exponent is similar to -0.25 calculated for the mass-specific O2 consumption. The inverse of cytosolic free [ADP], the Gibbs energy of ATP hydrolysis (delta G'ATP), and the efficiency of ATP production (energy captured in forming 3 mol of ATP per cycle along the mitochondrial respiratory chain from NADH to 1/2 O2) were all found to scale with body mass with a negative exponent. On the basis of scaling of the phosphorylation ratio and free cytosolic [ADP], we propose that the myocardium and other tissues of small mammals represent a metabolic system with a higher driving potential (a higher delta G'ATP from the higher [ATP]/[ADP][P(i)]) and a higher kinetic gain [(delta V/Vmax)/delta [ADP]] where small changes in free [ADP] produce large changes in steady-state rates of O2 consumption. From the inverse relationship between mitochondrial efficiency and body size we calculate that tissues of small mammals are more efficient than those of large mammals in converting energy from the oxidation of foodstuffs to the bond energy of ATP. A higher efficiency also indicates that mitochondrial electron transport is not the major site for higher heat production in small mammals. We further propose that the lower limit of about 2 g for adult endotherm body size (bumblebee-bat, Estrucan shrew, and hummingbird) may be set by the thermodynamics of the electron transport chain. The upper limit for body size (100,000-kg adult blue whale) may relate to a minimum delta G'ATP of approximately 55 kJ/mol for a cytoplasmic phosphorylation ratio of 12,000 M-1.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Interactions between tRNA identity nucleotides and their recognition sites in glutaminyl-tRNA synthetase determine the cognate amino acid affinity of the enzyme.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs both ensure accurate RNA recognition and prevent the binding of noncognate substrates. Here we show for Escherichia coli glutaminyl-tRNA synthetase (GlnRS; EC 6.1.1.18) that the accuracy of tRNA recognition also determines the efficiency of cognate amino acid recognition. Steady-state kinetics revealed that interactions between tRNA identity nucleotides and their recognition sites in the enzyme modulate the amino acid affinity of GlnRS. Perturbation of any of the protein-RNA interactions through mutation of either component led to considerable changes in glutamine affinity with the most marked effects seen at the discriminator base, the 10:25 base pair, and the anticodon. Reexamination of the identity set of tRNA(Gln) in the light of these results indicates that its constituents can be differentiated based upon biochemical function and their contribution to the apparent Gibbs' free energy of tRNA binding. Interactions with the acceptor stem act as strong determinants of tRNA specificity, with the discriminator base positioning the 3' end. The 10:25 base pair and U35 are apparently the major binding sites to GlnRS, with G36 contributing both to binding and recognition. Furthermore, we show that E. coli tryptophanyl-tRNA synthetase also displays tRNA-dependent changes in tryptophan affinity when charging a noncognate tRNA. The ability of tRNA to optimize amino acid recognition reveals a novel mechanism for maintaining translational fidelity and also provides a strong basis for the coevolution of tRNAs and their cognate synthetases.