Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.

A Specific Point Mutant at Position 1 of the Influenza Hemagglutinin Fusion Peptide Displays a Hemifusion Phenotype

We showed previously that substitution of the first residue of the influenza hemagglutinin (HA) fusion peptide Gly1 with Glu abolishes fusion activity. In the present study we asked whether this striking phenotype was due to the charge or side-chain volume of the substituted Glu. To do this we generated and characterized six mutants with substitutions at position 1: Gly1 to Ala, Ser, Val, Glu, Gln, or Lys. We found the following. All mutants were expressed at the cell surface, could be cleaved from the precursor (HA0) to the fusion permissive form (HA1-S-S-HA2), bound antibodies against the major antigenic site, bound red blood cells, and changed conformation at low pH. Only Gly, Ala, and Ser supported lipid mixing during fusion with red blood cells. Only Gly and Ala supported content mixing. Ser HA, therefore, displayed a hemifusion phenotype. The hemifusion phenotype of Ser HA was confirmed by electrophysiological studies. Our findings indicate that the first residue of the HA fusion peptide must be small (e.g., Gly, Ala, or Ser) to promote lipid mixing and must be small and apolar (e.g., Gly or Ala) to support both lipid and content mixing. The finding that Val HA displays no fusion activity underscores the idea that hydrophobicity is not the sole factor dictating fusion peptide function. The surprising finding that Ser HA displays hemifusion suggests that the HA ectodomain functions not only in the first stage of fusion, lipid mixing, but also, either directly or indirectly, in the second stage of fusion, content mixing.

Crystal structure of the HIV-1 integrase catalytic core and C-terminal domains: A model for viral DNA binding

Insolubility of full-length HIV-1 integrase (IN) limited previous structure analyses to individual domains. By introducing five point mutations, we engineered a more soluble IN that allowed us to generate multidomain HIV-1 IN crystals. The first multidomain HIV-1 IN structure is reported. It incorporates the catalytic core and C-terminal domains (residues 52–288). The structure resolved to 2.8 Å is a Y-shaped dimer. Within the dimer, the catalytic core domains form the only dimer interface, and the C-terminal domains are located 55 Å apart. A 26-aa α-helix, α6, links the C-terminal domain to the catalytic core. A kink in one of the two α6 helices occurs near a known proteolytic site, suggesting that it may act as a flexible elbow to reorient the domains during the integration process. Two proteins that bind DNA in a sequence-independent manner are structurally homologous to the HIV-1 IN C-terminal domain, suggesting a similar protein–DNA interaction in which the IN C-terminal domain may serve to bind, bend, and orient viral DNA during integration. A strip of positively charged amino acids contributed by both monomers emerges from each active site of the dimer, suggesting a minimally dimeric platform for binding each viral DNA end. The crystal structure of the isolated catalytic core domain (residues 52–210), independently determined at 1.6-Å resolution, is identical to the core domain within the two-domain 52–288 structure.

Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations

We describe and test a Markov chain model of microsatellite evolution that can explain the different distributions of microsatellite lengths across different organisms and repeat motifs. Two key features of this model are the dependence of mutation rates on microsatellite length and a mutation process that includes both strand slippage and point mutation events. We compute the stationary distribution of allele lengths under this model and use it to fit DNA data for di-, tri-, and tetranucleotide repeats in humans, mice, fruit flies, and yeast. The best fit results lead to slippage rate estimates that are highest in mice, followed by humans, then yeast, and then fruit flies. Within each organism, the estimates are highest in di-, then tri-, and then tetranucleotide repeats. Our estimates are consistent with experimentally determined mutation rates from other studies. The results suggest that the different length distributions among organisms and repeat motifs can be explained by a simple difference in slippage rates and that selective constraints on length need not be imposed.

Point mutant mice with hypersensitive α4 nicotinic receptors show dopaminergic deficits and increased anxiety

Knock-in mice were generated that harbored a leucine-to-serine mutation in the α4 nicotinic receptor near the gate in the channel pore. Mice with intact expression of this hypersensitive receptor display dominant neonatal lethality. These mice have a severe deficit of dopaminergic neurons in the substantia nigra, possibly because the hypersensitive receptors are continuously activated by normal extracellular choline concentrations. A strain that retains the neo selection cassette in an intron has reduced expression of the hypersensitive receptor and is viable and fertile. The viable mice display increased anxiety, poor motor learning, excessive ambulation that is eliminated by very low levels of nicotine, and a reduction of nigrostriatal dopaminergic function upon aging. These knock-in mice provide useful insights into the pathophysiology of sustained nicotinic receptor activation and may provide a model for Parkinson's disease.

Creating a Web-accessible, point-of-care, team-based information system (PoinTIS): the librarian as publisher*

The Internet has created new opportunities for librarians to develop information systems that are readily accessible at the point of care. This paper describes the multiyear process used to justify, fund, design, develop, promote, and evaluate a rehabilitation prototype of a point-of-care, team-based information system (PoinTIS) and train health care providers to use this prototype for their spinal cord injury and traumatic brain injury patient care and education activities. PoinTIS is a successful model for librarians in the twenty-first century to serve as publishers of information created or used by their parent organizations and to respond to the opportunities for information dissemination provided by recent technological advances.

Resolution of Holliday junctions in genetic recombination: RuvC protein nicks DNA at the point of strand exchange.

The RuvC protein of Escherichia coli catalyzes the resolution of recombination intermediates during genetic recombination and the recombinational repair of damaged DNA. Resolution involves specific recognition of the Holliday structure to form a complex that exhibits twofold symmetry with the DNA in an open configuration. Cleavage occurs when strands of like polarity are nicked at the sequence 5'-WTT decreases S-3' (where W is A or T and S is G or C). To determine whether the cleavage site needs to be located at, or close to, the point at which DNA strands exchange partners, Holliday structures were constructed with the junction points at defined sites within this sequence. We found that the efficiency of resolution was optimal when the cleavage site was coincident with the position of DNA strand exchange. In these studies, junction targeting was achieved by incorporating uncharged methyl phosphonates into the DNA backbone, providing further evidence for the importance of charge-charge repulsions in determining DNA structure.

Ras-catalyzed hydrolysis of GTP: a new perspective from model studies.

Despite the biological and medical importance of signal transduction via Ras proteins and despite considerable kinetic and structural studies of wild-type and mutant Ras proteins, the mechanism of Ras-catalyzed GTP hydrolysis remains controversial. We take a different approach to this problem: the uncatalyzed hydrolysis of GTP is analyzed, and the understanding derived is applied to the Ras-catalyzed reaction. Evaluation of previous mechanistic proposals from this chemical perspective suggests that proton abstraction from the attacking water by a general base and stabilization of charge development on the gamma-phosphoryl oxygen atoms would not be catalytic. Rather, this analysis focuses attention on the GDP leaving group, including the beta-gamma bridge oxygen of GTP, the atom that undergoes the largest change in charge in going from the ground state to the transition state. This leads to a new catalytic proposal in which a hydrogen bond from the backbone amide of Gly-13 to this bridge oxygen is strengthened in the transition state relative to the ground state, within an active site that provides a template complementary to the transition state. Strengthened transition state interactions of the active site lysine, Lys-16, with the beta-nonbridging phosphoryl oxygens and a network of interactions that positions the nucleophilic water molecule and gamma-phosphoryl group with respect to one another may also contribute to catalysis. It is speculated that a significant fraction of the GAP-activated GTPase activity of Ras arises from an additional interaction of the beta-gamma bridge oxygen with an Arg side chain that is provided in trans by GAP. The conclusions for Ras and related G proteins are expected to apply more widely to other enzymes that catalyze phosphoryl (-PO(3)2-) transfer, including kinases and phosphatases.

Evolutionarily conserved Galphabetagamma binding surfaces support a model of the G protein-receptor complex.

The pivotal role of G proteins in sensory, hormonal, inflammatory, and proliferative responses has provoked intense interest in understanding how they interact with their receptors and effectors. Nonetheless, the locations of the receptors and effector binding sites remain poorly characterized, although nearly complete structures of the alphabetagamma heterotrimeric complex are available. Here we apply evolutionary trace (ET) analysis [Lichtarge, O., Bourne, H. R. & Cohen, F. E. (1996) J. Mol. Biol. 257, 342-358] to propose plausible locations for these sites. On each subunit, ET identifies evolutionarily selected surfaces composed of residues that do not vary within functional subgroups and that form spatial clusters. Four clusters correctly identify subunit interfaces, and additional clusters on Galpha point to likely receptor or effector binding sites. Our results implicate the conformationally variable region of Galpha in an effector binding role. Furthermore the range of predicted interactions between the receptor and Galphabetagamma, is sufficiently limited that we can build a low resolution and testable model of the receptor-G protein complex.

Direct measurement of salt-bridge solvation energies using a peptide model system: implications for protein stability.

The solvation energies of salt bridges formed between the terminal carboxyl of the host pentapeptide AcWL- X-LL and the side chains of Arg or Lys in the guest (X) position have been measured. The energies were derived from octanol-to-buffer transfer free energies determined between pH 1 and pH 9. 13C NMR measurements show that the salt bridges form in the octanol phase, but not in the buffer phase, when the side chains and the terminal carboxyl group are charged. The free energy of salt-bridge formation in octanol is approximately -4 kcal/mol (1 cal = 4.184 J), which is equal to or slightly larger than the sum of the solvation energies of noninteracting pairs of charged side chains. This is about one-half the free energy that would result from replacing a charge pair in octanol with a pair of hydrophobic residues of moderate size. Therefore, salt bridging in octanol can change the favorable aqueous solvation energy of a pair of oppositely charged residues to neutral or slightly unfavorable but cannot provide the same free energy decrease as hydrophobic residues. This is consistent with recent computational and experimental studies of protein stability.