Biomedical subjects

T Simonson

Publications and source records attributed to T Simonson.

At least 19 recordsLinked to original sources

PubMed2001

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.

Aspartate-tRNA Ligase↗

PubMed2001

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.

Computer Simulation↗

PubMed2001

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.

Amino Acid Substitution↗

PubMed2001

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.

Acid-Base Equilibrium↗

PubMed1999

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.

Journal Article↗

PubMed1999

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.

Amino Acid Sequence↗

PubMed1999

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.

Amino Acid Sequence↗

PubMed1998

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.

Computer Simulation↗

PubMed1998

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.

Adenosine Triphosphate↗

PubMed1998

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.

Algorithms↗

PubMed1998

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.

Amino Acid Sequence↗

PubMed1997

MRI in neuro-Behçet's disease.

Our purpose was to characterise specific MRI findings and to determine their value in neuro-Behçet's disease. We examined 17 patients (14 men, 3 women) with neuro-Behçet's disease using T1- and T2-weighted spin-echo images and contrast-enhanced images at 0.5 T. There were 13 patients (76.5%) who had single or multiple lesions. Most of these were in the basal ganglia, brain stem or deep white matter region, giving high signal on T2-weighted images and isointense or low signal on T1-weighted images. In 3 cases (17.6%) there was linear high signal along the posterior limb of the internal capsule on T2-weighted images. This was considered as a potential differentiating feature of neuro-Behçet's disease. Contrast-enhancement was seen in 17 lesions in 7 patients.

Adult↗

PubMed1997

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.

Computer Simulation↗

PubMed1996

MRI of brucella polyneuritis in a child.

MRI of a rare case of polyneuritis following brucellosis is presented. An 8-year-old boy with a history of brucellosis showed marked enlargement of the origin of the left fifth cranial nerve and contrast enhancement of the fifth cranial nerves and the Gasserian ganglion in Meckel's caves. A well-defined low-signal cystic mass showing peripheral contrast enhancement lay anteromedial to the left seventh and eighth cranial nerves. The seventh and eighth nerves showed marked enhancement. The prechiasmatic portion of the right optic nerve showed fusiform enlargement with marked contrast enhancement of the nerve and right optic tract. MRI immediately after treatment showed no enhancement of the cranial nerves. Enlargement of the origin of the left fifth cranial nerve persisted, but there was no enlargement of the optic nerve. The mass near the left seventh and eighth cranial nerves was still present with no obvious morphological change, but no longer showed peripheral enhancement.

Brucellosis↗

PubMed1996

Polar fluctuations in proteins: molecular-dynamic studies of cytochrome c in aqueous solution.

The equilibrium fluctuations (the polar fluctuations), of yeast cytochrome c are studied using nanosecond molecular-dynamic simulations in a spherical droplet of water, with a particular emphasis on the fluctuations of the total dipole moment, which determine the average relative permittivity. These fluctuations follow a simple probability distribution, predicted by continuum electrostatics, and already observed in simulations of several polar liquids. An important component consists of diffusive, mutually independent, motions of the charged side chains at the protein surface. A quasiharmonic normal mode analysis of the trajectory shows that while motions covering a large range of frequencies contribute to the polar fluctuations, the four lowest frequency modes account for 50% of the overall static relative permittivity of ca. 25. The fluctuations of the protein bulk, i.e. parts other than the charged side chains, are distributed over a larger number of modes. Modes up to at least 60 cm-1 contribute to the average relative permittivity of the protein interior of ca. 4. The water surrounding the protein, despite the structural perturbation represented by the protein, has fluctuations similar to pure water, consistent with the idea of a linear solvent response to the protein charges. The relationship between the microscopic fluctuations seen in the simulations and simple continuum models is discussed.

Cytochrome c Group↗

PubMed1995

Internal and interfacial dielectric properties of cytochrome c from molecular dynamics in aqueous solution.

The dielectric properties of proteins are central to their stability and activity. We use the Fröhlich-Kirkwood theory of dielectrics to analyze two 1-ns molecular dynamics simulations of ferro- and ferricytochrome c in spherical droplets of 1400 water molecules. Protein and solvent are idealized as a series of concentric, spherical, dielectric media. Analysis results depend strongly on the treatment of the charged protein side chains at the protein/solvent interface. If charged side chains are viewed as part of the protein medium, then the protein dipole fluctuations are dominated by large, mutually uncorrelated, anisotropic, motions of the charged side chains. It is then incorrect to view the protein region as a single, homogeneous dielectric material. If one does take this view, estimates of the protein "dielectric constant" vary from 16 to 37, depending on the exact choice of model parameters. In contrast, if the charged portions of the charged side chains are viewed as part of the solvent medium, then theory and simulation are consistent: the protein dipole fluctuations excluding charged side chains are roughly those of a homogeneous, isotropic dielectric medium, with a dielectric constant of 4.7 +/- 1.0 (ferro) or 3.4 +/- 1.0 (ferri), in agreement with powder experiments. Statistical uncertainty and sensitivity to model parameters are small. Analysis of the radial dependence of the dipole fluctuations suggests that the inner half of the protein has a somewhat lower dielectric constant of 1.5-2, consistent with its biological function in electron transfer. These results suggest that Poisson-Boltzmann models could treat the protein bulk as a low-dielectric medium and the charged surface groups as part of the solvent region.

Cytochrome c Group↗

PubMed1995

Proline cis-trans isomerization in staphylococcal nuclease: multi-substrate free energy perturbation calculations.

Staphylococcal nuclease A exists in two folded forms that differ in the isomerization state of the Lys 116-Pro 117 peptide bond. The dominant form (90% occupancy) adopts a cis peptide bond, which is observed in the crystal structure. NMR studies show that the relatively small difference in free energy between the cis and trans forms (delta Gcis-->trans approximately 1.2 kcal/mol) results from large and nearly compensating differences in enthalpy and entropy (delta Hcis-->trans approximately delta TScis-->trans approximately 10 kcal/mol). There is evidence from X-ray crystal structures that the structural differences between the cis and the trans forms of nuclease are confined to the conformation of residues 112-117, a solvated protein loop. Here, we obtain a thermodynamic and structural description of the conformational equilibrium of this protein loop through an exhaustive conformational search that identified several substates followed by free energy simulations between the substrates. By partitioning the search into conformational substates, we overcame the multiple minima problem in this particular case and obtained precise and reproducible free energy values. The protein and water environment was implicitly modeled by appropriately chosen nonbonded terms between the explicitly treated loop and the rest of the protein. These simulations correctly predicted a small free energy difference between the cis and trans forms composed of larger, compensating differences in enthalpy and entropy. The structural predictions of these simulations were qualitatively consistent with known X-ray structures of nuclease variants and yield a model of the unknown minor trans conformation.

Amino Acid Sequence↗

PubMed1995

The role of magnetic resonance angiography in head and neck surgery.

Magnetic resonance angiography (MRA) is a recently developed, noninvasive vascular imaging technique. The authors of this investigation assessed the diagnostic value of MRA, along with its influence on therapeutic decisions, in 11 patients with a variety of head and neck disorders. In 5 patients, MRA diagnosed or ruled out an intrinsic vascular lesion. MRA was used to evaluate 5 of 8 patients with cancer for evidence of direct tumor involvement of vascular structures. Other uses of MRA included preoperative determination of tumor vascularity and delineation of anatomic relationships between normal vessels and head and neck pathology. Overall, MRA results guided management in 10 patients, and in some cases it determined the extent of surgical intervention. Because MRA is safer and more practical than traditional angiography, the authors recommend more frequent use of this imaging technique in the practice of head and neck surgery.

Adult↗