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.

Low dimensionality of supraspinally induced force fields

Recent experiments using electrical and N-methyl-d-aspartate microstimulation of the spinal cord gray matter and cutaneous stimulation of the hindlimb of spinalized frogs have provided evidence for a modular organization of the frog’s spinal cord circuitry. A “module” is a functional unit in the spinal cord circuitry that generates a specific motor output by imposing a specific pattern of muscle activation. The output of a module can be characterized as a force field: the collection of the isometric forces generated at the ankle over different locations in the leg’s workspace. Different modules can be combined independently so that their force fields linearly sum. The goal of this study was to ascertain whether the force fields generated by the activation of supraspinal structures could result from combinations of a small number of modules. We recorded a set of force fields generated by the electrical stimulation of the vestibular nerve in seven frogs, and we performed a principal component analysis to study the dimensionality of this set. We found that 94% of the total variation of the data is explained by the first five principal components, a result that indicates that the dimensionality of the set of fields evoked by vestibular stimulation is low. This result is compatible with the hypothesis that vestibular fields are generated by combinations of a small number of spinal modules.

Stability of nonrelativistic quantum mechanical matter coupled to the (ultraviolet cutoff) radiation field

We announce a proof of H-stability for the quantized radiation field, with ultraviolet cutoff, coupled to arbitrarily many non-relativistic quantized electrons and static nuclei. Our result holds for arbitrary atomic numbers and fine structure constant. We also announce bounds for the energy of many electrons and nuclei in a classical vector potential and for the eigenvalue sum of a one-electron Pauli Hamiltonian with magnetic field.

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.

Cosmic velocity fields and their interpretation

We review the current status of our knowledge of cosmic velocity fields, on both small and large scales. A new statistic is described that characterizes the incoherent, thermal component of the velocity field on scales less than 2h−1 Mpc (h is H0/100 km·s−1·Mpc−1, where H0 is the Hubble constant and 1 Mpc = 3.09 × 1022 m) and smaller. The derived velocity is found to be quite stable across different catalogs and is of remarkably low amplitude, consistent with an effective Ω ∼ 0.15 on this scale. We advocate the use of this statistic as a standard diagnostic of the small-scale kinetic energy of the galaxy distribution. The analysis of large-scale flows probes the velocity field on scales of 10–60 h−1 Mpc and should be adequately described by linear perturbation theory. Recent work has focused on the comparison of gravity or density fields derived from whole-sky redshift surveys of galaxies [e.g., the Infrared Astronomical Satellite (IRAS)] with velocity fields derived from a variety of sources. All the algorithms that directly compare the gravity and velocity fields suggest low values of the density parameter, while the POTENT analysis, using the same data but comparing the derived IRAS galaxy density field with the Mark-III derived matter density field, leads to much higher estimates of the inferred density. Since the IRAS and Mark-III fields are not fully consistent with each other, the present discrepancies might result from the very different weighting applied to the data in the competing methods.

Generation of a high-titer retroviral vector capable of expressing high levels of the human beta-globin gene.

Retrovirus-mediated gene transfer into hematopoietic cells may provide a means of treating both inherited and acquired diseases involving hematopoietic cells. Implementation of this approach for disorders resulting from mutations affecting the beta-globin gene (e.g., beta-thalassemia and sickle cell anemia), however, has been hampered by the inability to generate recombinant viruses able to efficiently and faithfully transmit the necessary sequences for appropriate gene expression. We have addressed this problem by carefully examining the interactions between retroviral and beta-globin gene sequences which affect vector transmission, stability, and expression. First, we examined the transmission properties of a large number of different recombinant proviral genomes which vary both in the precise nature of vector, beta-globin structural gene, and locus control region (LCR) core sequences incorporated and in the placement and orientation of those sequences. Through this analysis, we identified one specific vector, termed M beta 6L, which carries both the human beta-globin gene and core elements HS2, HS3, and HS4 from the LCR and faithfully transmits recombinant proviral sequences to cells with titers greater than 10(6) per ml. Populations of murine erythroleukemia (MEL) cells transduced by this virus expressed levels of human beta-globin transcript which, on a per gene copy basis, were 78% of the levels detected in an MEL-derived cell line, Hu11, which carries human chromosome 11, the site of the beta-globin locus. Analysis of individual transduced MEL cell clones, however, indicated that, while expression was detected in every clone tested (n = 17), the levels of human beta-globin treatment varied between 4% and 146% of the levels in Hu11. This clonal variation in expression levels suggests that small beta-globin LCR sequences may not provide for as strict chromosomal position-independent expression of beta-globin as previously suspected, at least in the context of retrovirus-mediated gene transfer.