Tibetans exhibit lower hemoglobin concentration and decreased heart response to hypoxia during poikilocapnia at intermediate altitude relative to Han Chinese.

Background: High-altitude populations exhibit distinct cellular, respiratory, and cardiovascular phenotypes, some of which provide adaptive advantages to hypoxic conditions compared to populations with sea-level ancestry. Studies performed in populations with a history of high-altitude residence, such as Tibetans, support the idea that many of these phenotypes may be shaped by genomic features that have been positively selected for throughout generations. We hypothesize that such traits observed in Tibetans at high altitude also occur in Tibetans living at intermediate altitude, even in the absence of severe sustained hypoxia. Methodology: We studied individuals of high-altitude ancestry (Tibetans, n = 17 females; n = 12 males) and sea-level ancestry (Han Chinese, n = 6 females; n = 10 males), both who had been living at ∼1300 m (∼4327 ft) for at least 18 months. We measured hemoglobin concentration ([Hb]), hypoxic ventilatory response (HVR), and hypoxic heart rate response (HHRR) with end-tidal CO2 (PetCO2) held constant (isocapnia) or allowed to decrease with hypoxic hyperventilation (poikilocapnia). We also quantified the contribution of CO2 on ventilation and heart rate by calculating the differences of isocapnic versus poikilocapnic hypoxic conditions (Δ V˙I/ΔPetCO2 and ΔHR/ΔPetCO2, respectively). Results: Male Tibetans had lower [Hb] compared to Han Chinese males (p < 0.05), consistent with reports for individuals from these populations living at high altitude and sea level. Measurements of ventilation (resting ventilation, HVR, and PetCO2) were similar for both groups. Heart rate responses to hypoxia were similar in both groups during isocapnia; however, HHRR in poikilocapnia was reduced in the Tibetan group (p < 0.03), and the heart rate response to CO2 in hypoxia was lower in Tibetans relative to Han Chinese (p < 0.01). Conclusion: These results suggest that Tibetans living at intermediate altitude have blunted cardiac responses in the context of hypoxia. Hence, only some of the phenotypes observed in Tibetans living at high altitude are observed in Tibetans living at intermediate altitude. Whereas blunted cardiac responses to hypoxia is revealed at intermediate altitudes, manifestation of other physiological adaptations to high altitude may require exposure to more severe levels of hypoxia.

Sea-level haemoglobin concentration is associated with greater exercise capacity in Tibetan males at 4200 m.

NEW FINDINGS: What is the topic of this review? Recent developments link relatively lower hemoglobin concentration in Tibetans at high altitude to exercise capacity and components of oxygen transport. What advances does it highlight? Haemoglobin concentration (ranging from 15.2 to 22.9 g dl(-1) ) in Tibetan males was negatively associated with peak oxygen (O2 ) uptake per kilogram, cardiac output and muscle O2 diffusion conductance. Most variance in the peak O2 uptake per kilogram of Tibetan males was attributed to cardiac output, muscle diffusional conductance and arterial partial pressure of CO2 . The mechanisms underlying these differences in oxygen transport in Tibetans require additional analyses. Despite residence at >4000 m above sea level, many Tibetan highlanders, unlike Andean counterparts and lowlanders at altitude, exhibit haemoglobin concentration ([Hb]) within the typical sea-level range. Genetic adaptations in Tibetans are associated with this relatively low [Hb], yet the functional relevance of the lower [Hb] remains unknown. To address this, we examined each major step of the oxygen transport cascade [ventilation (VE), cardiac output (QT) and diffusional conductance in lung (DL) and muscle (DM)] in Tibetan males at maximal exercise on a cycle ergometer. Ranging from 15.2 to 22.9 g dl(-1) , [Hb] was negatively associated with peak O2 uptake per kilogram (r = -0.45, P < 0.05) and both cardiac output (QT/kg: r = -0.54, P < 0.02) and muscle O2 diffusion conductance (DM/kg: r = -0.44, P < 0.05) but not ventilation, arterial partial pressure of O2 or pulmonary diffusing capacity. Most variance in peak O2 uptake per kilogram was attributed to QT, DM and arterial partial pressure of CO2 (r(2)  = 0.90). In summary, lack of polycythaemia in Tibetans is associated with increased exercise capacity, which is explained by elevated cardiac, muscle and, to a small extent, ventilatory responses rather than pulmonary gas exchange. Whether lower [Hb] is the cause or result of these changes in O2 transport or is causally unrelated will require additional study.

Adaptive genetic changes related to haemoglobin concentration in native high-altitude Tibetans.

NEW FINDINGS: What is the topic of this review? Tibetans have genetic adaptations that are hypothesized to underlie the distinct set of traits they exhibit at altitude. What advances does it highlight? Several adaptive signatures in the same genomic regions have been identified among Tibetan populations resident throughout the Qinghai-Tibetan Plateau. Many highland Tibetans exhibit a haemoglobin concentration within the range expected at sea level, and this trait is associated with putatively adaptive regions harbouring the hypoxia-inducible factor pathway genes EGLN1, EPAS1 and PPARA. Precise functional variants at adaptive loci and relationships to physiological traits, beyond haemoglobin concentration, are currently being examined in this population. Some native Tibetan, Andean and Ethiopian populations have lived at altitudes ranging from 3000 to >4000 m above sea level for hundreds of generations and exhibit distinct combinations of traits at altitude. It was long hypothesized that genetic factors contribute to adaptive differences in these populations, and recent advances in genomics provide evidence that some of the strongest signatures of positive selection in humans are those identified in Tibetans. Many of the top adaptive genomic regions highlighted thus far harbour genes related to hypoxia sensing and response. Putatively adaptive copies of three hypoxia-inducible factor pathway genes, EPAS1, EGLN1 and PPARA, are associated with sea-level range, rather than elevated, haemoglobin concentration observed in many Tibetans at high altitude, and recent studies provide insight into some of the precise adaptive variants, timing of adaptive events and functional roles. While several studies in highland Tibetans have converged on a few hypoxia-inducible factor pathway genes, additional candidates have been reported in independent studies of Tibetans located throughout the Qinghai-Tibetan Plateau. Various aspects of adaptive significance have yet to be identified, integrated, and fully explored. Given the rapid technological advances and interdisciplinary efforts in genomics, physiology and molecular biology, careful examination of Tibetans and comparisons with other distinctively adapted highland populations will provide valuable insight into evolutionary processes and models for both basic and clinical research.

Low haemoglobin concentration in Tibetan males is associated with greater high-altitude exercise capacity.

Tibetans living at high altitude have adapted genetically such that many display a low erythropoietic response, resulting in near sea-level haemoglobin (Hb) concentration. We hypothesized that absence of the erythropoietic response would be associated with greater exercise capacity compared to those with high [Hb] as a result of beneficial changes in oxygen transport. We measured, in 21 Tibetan males with [Hb] ranging from 15.2 g dl(-1) to 22.9 g dl(-1) (9.4 mmol l(-1) to 14.2 mmol l(-1) ), [Hb], ventilation, volumes of O2 and CO2 utilized at peak exercise (V̇O2 and V̇CO2), heart rate, cardiac output and arterial blood gas variables at peak exercise on a cycle ergometer at ∼4200 m. Lung and muscle O2 diffusional conductances were computed from these measurements. [Hb] was related (negatively) to V̇O2 kg(-1) (r = -0.45, P< 0.05), cardiac output kg(-1) (QT kg(-1) , r = -0.54, P < 0.02), and O2 diffusion capacity in muscle (DM kg(-1) , r = -0.44, P<0.05), but was unrelated to ventilation, arterial partial pressure of O2 (PaO2) or pulmonary diffusing capacity. Using multiple linear regression, variance in peak V̇O2 kg(-1) was primarily attributed to QT, DM, and PCO2 (R(2) = 0.88). However, variance in pulmonary gas exchange played essentially no role in determining peak V̇O2. These results (1) show higher exercise capacity in Tibetans without the erythropoietic response, supported mostly by cardiac and muscle O2 transport capacity and ventilation rather than pulmonary adaptations, and (2) support the emerging hypothesis that the polycythaemia of altitude, normally a beneficial response to low cellular PO2, may become maladaptive if excessively elevated under chronic hypoxia. The cause and effect relationships among [Hb], QT, DM, and PCO2 remain to be elucidated.

Increased blood-oxygen binding affinity in Tibetan and Han Chinese residents at 4200 m.

High-altitude natives are challenged by hypoxia, and a potential compensatory mechanism could be reduced blood oxygen-binding affinity (P50), as seen in several high-altitude mammalian species. In 21 Qinghai Tibetan and nine Han Chinese men, all resident at 4200 m, standard P50 was calculated from measurements of arterial PO2 and forehead oximeter oxygen saturation, which was validated in a separate examination of 13 healthy subjects residing at sea level. In both Tibetans and Han Chinese, standard P50 was 24.5 ± 1.4 and 24.5 ± 2.0 mmHg, respectively, and was lower than in the sea-level subjects (26.2 ± 0.6 mmHg, P < 0.01). There was no relationship between P50 and haemoglobin concentration (the latter ranging from 15.2 to 22.9 g dl(-1) in Tibetans). During peak exercise, P50 was not associated with alveolar-arterial PO2 difference or peak O2 uptake per kilogram. There appears to be no apparent benefit of a lower P50 in this adult high-altitude Tibetan population.

Comparison of motor unit action potential characteristics and hand dominance using monopolar needle electrodes in the abductor pollicis brevis and abductor digiti minimi muscles.

This study examined specific electrical characteristics of voluntary single motor unit action potentials (SMUAPs): amplitude, duration, phase change, and rate of rise. These characteristics, which were detected from two intrinsic muscles of the hand--the abductor pollicis brevis and the abductor digiti minimi--were compared to hand dominance. Forty subjects participated in the study. Five characteristics were detected from each muscle using the quadrant technique while the subject produced a minimal isometric contraction. Based on results of an ANOVA [two-factor with replication] test, our study revealed no significant difference between muscles in the dominant and non-dominant hands. Descriptive statistics for each muscle characteristic are presented. This study has identified parameters for SMUAP characteristics detected in non-impaired individuals ranging in age from 20 to 43 years. The normative parameters serve as a valuable base from which one may examine potential neuronal damage from cumulative trauma disorders.

Dielectric relaxation in an enzyme active site: molecular dynamics simulations interpreted with a macroscopic continuum model.

Dielectric relaxation plays an important role in many chemical processes in proteins, including acid-base titration, ligand binding, and charge transfer reactions. Its complexity makes experimental characterization difficult, and so, theoretical approaches are valuable. The comparison of molecular dynamics free energy simulations with simpler models such as a dielectric continuum model is especially useful for obtaining qualitative insights. We have analyzed a charge insertion process that models deprotonation or mutation of an important side chain in the active site of the enzyme aspartyl-tRNA synthetase. Complexes with the substrate aspartate and the analogue asparagine were studied. The resulting dielectric relaxation was found to involve both ligand and side chain rearrangements in the active site and to account for a large part of the overall charging free energy. With the continuum model, charge insertion is performed along a two-step pathway: insertion into a static environment, followed by relaxation of the environment. These correspond to different physical processes and require different protein dielectric constants. A low value of approximately 1 is needed for the static step, consistent with the parametrization of the molecular mechanics charge set used. A value of 3-6 (depending on the exact insertion site and the nature of the ligand) is needed to describe the dielectric relaxation step. This moderate value indicates that, for this system, the local protein polarizability in the active site is within at most a factor of 2 of that expected at nonspecific positions in a protein interior.

Protein molecular dynamics with the generalized Born/ACE solvent model.

Implicit solvent models are increasingly important for the study of proteins in aqueous solution. Here, the generalized Born (GB) solvent polarization model as implemented in the analytical ACE potential [Schaefer and Karplus (1996) J Phys Chem 100:1578] is used to perform molecular dynamics simulations of two small, homologous proteins: the immunoglobulin-binding domain of streptococcal protein G and the Ras binding domain of Raf. Several model parameterizations are compared through more than 60 ns of simulation. Results are compared with two simpler solvent models-an accessible surface area model and a distant-dependent dielectric model, with finite-difference Poisson calculations, with existing explicit solvent simulations, and with experimental data. The simpler models yield stable but distorted structures. The best GB/ACE implementation uses a set of atomic Voronoi volumes reported recently, obtained by averaging over a large database of crystallographic protein structures. A 20% reduction is applied to the volumes, compensating in an average sense for an excessive de-screening of individual charges inherent in the ACE self-energy and for an undersolvation of dipolar groups inherent in the GB screening function. This GB/ACE parameterization yields stable trajectories on the 0.5-1-ns time scale that deviate moderately (approximately 1.5-2.5 A) from the X-ray structure, reproduce approximately the surface distribution of charged, polar, and hydrophobic groups, and reproduce accurately backbone flexibility as measured by amide NMR-order parameters. Over longer time scales (1.5-3 ns), some of the protein G runs escape from the native energy basin and deviate strongly (3 A) from the native structure. The conformations sampled during the transition out of the native energy basin are overstabilized by the GB/ACE solvation model, as compared with a numerical treatment of the full dielectric continuum model.

Brain MR imaging in the evaluation of chronic headache in patients without other neurologic symptoms.

RATIONALE AND OBJECTIVES: The authors investigated the use of magnetic resonance (MR) imaging of the brain in adult patients with a primary complaint of chronic headache and no other neurologic symptoms or findings and determined the yield and MR predictors of major abnormalities in these patients. MATERIALS AND METHODS: The medical records and MR images of 402 adult patients with chronic headache were retrospectively reviewed. All patients had been evaluated and referred by the neurology service. The findings were categorized as either negative or positive for major abnormality. Multivariate analysis with a linear logistic regression technique was performed on the clinical data, which included patient age, patient sex, and headache type. RESULTS: Major abnormalities were found in 15 patients (3.7%), consisting of seven women (2.4%) and eight men (6.9%). Major abnormalities were found in 0.6% of those with migraine headaches, 1.4% with tension headaches, none with mixed migraine and tension headaches, 14.1% with atypical headaches, and 3.8% with other types of headaches. Multivariate analysis showed that the atypical headache type was the most significant predictor of major abnormality. CONCLUSION: The yield of major abnormalities found with brain MR imaging in patients with isolated chronic headache is low. However, those patients with atypical headaches have a higher yield of major abnormalities and may benefit from imaging.

Factors influencing visualization of vertebral metastases on MR imaging versus bone scintigraphy.

OBJECTIVE: The purpose of this study was to investigate whether the location and size of vertebral body metastases influence the difference in detection rates between MR imaging and bone scintigraphy. MATERIALS AND METHODS: We retrospectively evaluated the vertebral body lesions detected on MR imaging in 74 patients with known widely disseminated metastatic disease. Three radiologists independently reviewed the MR images and bone scintigraphs. MR imaging findings included lesion size and its spatial relationship to the bony cortex (intramedullary, subcortical, and transcortical) and results were correlated with those of planar technetium 99m bone scintigraphy. RESULTS: Findings on bone scans were negative for all intramedullary lesions without cortical involvement shown on MR imaging, regardless of their size. Findings on bone scans (71.3% for transcortical and 33.8% for subcortical) were frequently positive for lesions with cortical involvement (trans- or subcortical), and the probability of positive findings on bone scans was also influenced by the lesion size. Statistical analysis showed a positive correlation among cortical involvement, lesion size, and positive findings on bone scintigraphy (p < 0.0001). CONCLUSION: Location (the presence of cortical bone involvement on MR imaging) and size of the vertebral body metastases appear to be important contributing factors to the difference in detection rates between MR imaging and bone scintigraphy. Cortical involvement is likely the cause of positive findings on bone scans. Early vertebral metastases tend to be small and located in the medullary cavity without cortical involvement, and therefore, findings may be positive on MR images but negative on bone scans.

Macromolecular electrostatics: continuum models and their growing pains.

Theoretical understanding of macromolecular electrostatics has advanced substantially over the past year. Continuum models have given promising results for calculating protein-ligand binding free energy differences, as well as pK(a)s and redox properties, particularly with explicit treatment of multiple conformers. Generalized Born and other techniques have led to the first molecular dynamics simulations of proteins and RNA with continuum solvent. Continuum and microscopic descriptions of dielectric relaxation have been critically compared.

Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase.

Specific amino acid binding by aminoacyl-tRNA synthetases (aaRS) is necessary for correct translation of the genetic code. Engineering a modified specificity into aminoacyl-tRNA synthetases has been proposed as a means to incorporate artificial amino acid residues into proteins in vivo. In a previous paper, the binding to aspartyl-tRNA synthetase of the substrate Asp and the analogue Asn were compared by molecular dynamics free energy simulations. Molecular dynamics combined with Poisson-Boltzmann free energy calculations represent a less expensive approach, suitable for examining multiple active site mutations in an engineering effort. Here, Poisson-Boltzmann free energy calculations for aspartyl-tRNA synthetase are first validated by their ability to reproduce selected molecular dynamics binding free energy differences, then used to examine the possibility of Asn binding to native and mutant aspartyl-tRNA synthetase. A component analysis of the Poisson-Boltzmann free energies is employed to identify specific interactions that determine the binding affinities. The combined use of molecular dynamics free energy simulations to study one binding process thoroughly, followed by molecular dynamics and Poisson-Boltzmann free energy calculations to study a series of related ligands or mutations is proposed as a paradigm for protein or ligand design. The binding of Asn in an alternate, "head-to-tail" orientation observed in the homologous asparagine synthetase is analyzed, and found to be more stable than the "Asp-like" orientation studied earlier. The new orientation is probably unsuitable for catalysis. A conserved active site lysine (Lys198 in Escherichia coli) that recognizes the Asp side-chain is changed to a leucine residue, found at the corresponding position in asparaginyl-tRNA synthetase. It is interesting that the binding of Asp is calculated to increase slightly (rather than to decrease), while that of Asn is calculated, as expected, to increase strongly, to the same level as Asp binding. Insight into the origin of these changes is provided by the component analyses. The double mutation (K198L,D233E) has a similar effect, while the triple mutation (K198L,Q199E,D233E) reduces Asp binding strongly. No binding measurements are available, but the three mutants are known to have no ability to adenylate Asn, despite the "Asp-like" binding affinities calculated here. In molecular dynamics simulations of all three mutants, the Asn ligand backbone shifts by 1-2 A compared to the experimental Asp:AspRS complex, and significant side-chain rearrangements occur around the pocket. These could reduce the ATP binding constant and/or the adenylation reaction rate, explaining the lack of catalytic activity in these complexes. Finally, Asn binding to AspRS with neutral K198 or charged H449 is considered, and shown to be less favorable than with the charged K198 and neutral H449 used in the analysis.

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding.

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.

Molecular dynamics simulations of the Ras:Raf and Rap:Raf complexes.

The protein Raf is an immediate downstream target of Ras in the MAP kinase signalling pathway. The complex of Ras with the Ras-binding domain (RBD) of Raf has been modelled by homology to the (E30D,K31E)-Rap1A:RBD complex, and both have been subjected to multiple molecular dynamics simulations in solution. While both complexes are stable, several rearrangements occur in the Ras:RBD simulations: the RBD loop 100-109 moves closer to Ras, Arg73 in the RBD moves towards Ras to form a salt bridge with Ras-Asp33, and Loop 4 of the Ras switch II region shifts upwards toward the RBD. The Ras:RBD interactions (including the RBD-Arg73 interaction) are consistent with available NMR and mutagenesis data on the Ras: RBD complex in solution. The Ras switch II region does not interact directly with the RBD, although indirect interactions exist through the effector domain and bridging water molecules. No large-scale RBD motion is seen in the Ras:RBD complex, compared to the Rap:RBD complex, to suggest an allosteric activation of Raf by Ras. This may be because the Raf kinase domain (whose structure is unknown) is not included in the model.

Implicit solvent models.

Implicit solvent models for biomolecular simulations are reviewed and their underlying statistical mechanical basis is discussed. The fundamental quantity that implicit models seek to approximate is the solute potential of mean force, which determines the statistical weight of solute conformations, and which is obtained by averaging over the solvent degrees of freedom. It is possible to express the total free energy as the reversible work performed in two successive steps. First, the solute is inserted in the solvent with zero atomic partial charges; second, the atomic partial charges of the solute are switched from zero to their full values. Consequently, the total solvation free energy corresponds to a sum of non-polar and electrostatic contributions. These two contributions are often approximated by simple geometrical models (such as solvent exposed area models) and by macroscopic continuum electrostatics, respectively. One powerful route is to approximate the average solvent density distribution around the solute, i.e. the solute-solvent density correlation functions, as in statistical mechanical integral equations. Recent progress with semi-analytical approximations makes continuum electrostatics treatments very efficient. Still more efficient are fully empirical, knowledge-based models, whose relation to explicit solvent treatments is not fully resolved, however. Continuum models that treat both solute and solvent as dielectric continua are also discussed, and the relation between the solute fluctuations and its macroscopic dielectric constant(s) clarified.

Crystallography & NMR system: A new software suite for macromolecular structure determination.

A new software suite, called Crystallography & NMR System (CNS), has been developed for macromolecular structure determination by X-ray crystallography or solution nuclear magnetic resonance (NMR) spectroscopy. In contrast to existing structure-determination programs, the architecture of CNS is highly flexible, allowing for extension to other structure-determination methods, such as electron microscopy and solid-state NMR spectroscopy. CNS has a hierarchical structure: a high-level hypertext markup language (HTML) user interface, task-oriented user input files, module files, a symbolic structure-determination language (CNS language), and low-level source code. Each layer is accessible to the user. The novice user may just use the HTML interface, while the more advanced user may use any of the other layers. The source code will be distributed, thus source-code modification is possible. The CNS language is sufficiently powerful and flexible that many new algorithms can be easily implemented in the CNS language without changes to the source code. The CNS language allows the user to perform operations on data structures, such as structure factors, electron-density maps, and atomic properties. The power of the CNS language has been demonstrated by the implementation of a comprehensive set of crystallographic procedures for phasing, density modification and refinement. User-friendly task-oriented input files are available for nearly all aspects of macromolecular structure determination by X-ray crystallography and solution NMR.

Engineering an Mg2+ site to replace a structurally conserved arginine in the catalytic center of histidyl-tRNA synthetase by computer experiments.

Histidyl-tRNA synthetase (HisRS) differs from other class II aminoacyl-tRNA synthetases (aaRS) in that it harbors an arginine at a position where the others bind a catalytic Mg2+ ion. In computer experiments, four mutants of HisRS from Escherichia coli were engineered by removing the arginine and introducing a Mg2+ ion and residues from seryl-tRNA synthetase (SerRS) that are involved in Mg2+ binding. The mutants recreate an active site carboxylate pair conserved in other class II aaRSs, in two possible orders: Glu-Asp or Asp-Glu, replacing Glu-Thr in native HisRS. The mutants were simulated by molecular dynamics in complex with histidyl-adenylate. As controls, the native HisRS was simulated in complexes with histidine, histidyl-adenylate, and histidinol. The native structures sampled were in good agreement with experimental structures and biochemical data. The two mutants with the Glu-Asp sequence showed significant differences in active site structure and Mg2+ coordination from SerRS. The others were more similar to SerRS, and one of them was analyzed further through simulations in complex with histidine, and His+ATP. The latter complex sampled two Mg2+ positions, depending on the conformation of a loop anchoring the second carboxylate. The lowest energy conformation led to an active site geometry very similar to SerRS, with the principal Mg2+ bridging the alpha- and beta-phosphates, the first carboxylate (Asp) coordinating the ion through a water molecule, and the second (Glu) coordinating it directly. This mutant is expected to be catalytically active and suggests a basis for the previously unexplained conservation of the active site Asp-Glu pair in class II aaRSs other than HisRS.

Conformation of the Ras-binding domain of Raf studied by molecular dynamics and free energy simulations.

Recognition of Ras by its downstream target Raf is mediated by a Ras-recognition region in the Ras-binding domain (RBD) of Raf. Residues 78-89 in this region occupy two different conformations in the ensemble of NMR solution structures of the RBD: a fully alpha-helical one, and one where 87-90 form a type IV beta-turn. Molecular dynamics simulations of the RBD in solution were performed to explore the stability of these and other possible conformations of both the wild-type RBD and the R89K mutant, which does not bind Ras. The simulations sample a fully helical conformation for residues 78-89 similar to the NMR helical structures, a conformation where 85-89 form a 3(10)-helical turn, and a conformation where 87-90 form a type I beta-turn, whose free energies are all within 0.3 kcal/mol of each other. NOE patterns and H(alpha) chemical shifts from the simulations are in reasonable agreement with experiment. The NMR turn structure is calculated to be 3 kcal/mol higher than the three above conformations. In a simulation with the same implicit solvent model used in the NMR structure generation, the turn conformation relaxes into the fully helical conformation, illustrating possible structural artifacts introduced by the implicit solvent model. With the Raf R89K mutant, simulations sample a fully helical and a turn conformation, the turn being 0.9 kcal/mol more stable. Thus, the mutation affects the population of RBD conformations, and this is expected to affect Ras binding. For example, if the fully helical conformation of residues 78-89 is required for binding, its free energy increase in R89K will increase the binding free energy by about 0.6 kcal/mol.

Specific amino acid recognition by aspartyl-tRNA synthetase studied by free energy simulations.

Specific amino acid binding by aminoacyl-tRNA synthetases is necessary for correct translation of the genetic code. To obtain insight into the origin of the specificity, the binding to aspartyl-tRNA synthetase (AspRS) of the negatively charged substrate aspartic acid and the neutral analogue asparagine are compared by use of molecular dynamics and free energy simulations. Simulations of the Asn-AspRS complex show that although Asn cannot bind in the same position as Asp, several possible positions exist 1.5 to 2 A away from the Asp site. The binding free energy of Asn in three of these positions was compared to that of Asp through alchemical free energy simulations, in which Asp is gradually mutated ito Asn in the complex with the enzyme. To correctly account for the electrostatic interactions in the system (including bulk solvent), a recently developed hybrid approach was used, in which the region of the mutation site is treated microscopically, whereas distant protein and solvent are treated by continuum electrostatics. Seven free energy simulations were performed in the protein and two in solution. The various Asn positions and orientations sampled at the Asn endpoints of the protein simulations yielded very similar free energy differences. The calculated Asp-->Asn free energy change is 79.8(+/-1.5) kcal/mol in solution and 95.1(+/-2.8) kcal/mol in the complex with the protein. Thus, the substrate Asp is predicted to bind much more strongly than Asn, with a binding free energy difference of 15.3 kcal/mol. This implies that erroneous binding of Asn by AspRS is highly improbable, and cannot account for any errors in the translation of the genetic code. Almost all of the protein contributions to the Asp versus Asn binding free energy difference arise from an arginine and a lysine residue that hydrogen bond to the substrate carboxylate group and an Asp and a Glu that hydrogen bond to these; all four amino acid residues are completely conserved in AspRSs. The protein effectively "solvates" the Asp side-chain more strongly than water does. The simulations are analyzed to determine the interactions that Asn is able to make in the binding pocket, and which sequence differences between AspRS and the highly homologous AsnRS are important for modifying the amino acid specificity. A double or triple mutation of AspRS that could make it specific for Asn is proposed, and supported by preliminary simulations of a mutant complex.

Redox properties of cytochrome c: novel linear response and hybrid continuum-microscopic methodologies.

Redox properties of yeast cytochrome c are estimated using molecular dynamics combined with a simple linear response approximation, as well as a hybrid continuum-molecular dynamics (COMD) approach. In both approaches, the free energy associated with an electrostatic perturbation (a redox electron) is separated into its relaxation and static (non-relaxation) components. The static component is calculated from the molecular dynamics simulation. The relaxation component is then calculated with a linear response approximation, either from the molecular dynamics, or from a separate continuum calculation. This latter hybrid approach exploits the relative robustness of continuum models for dealing with large perturbations, while avoiding some of their limitations. It is quite general, and could be applied for example to pKa calculations.