Showing papers by "Wilfred F. van Gunsteren published in 2012"

On developing coarse-grained models for biomolecular simulation: a review.

Abstract: So-called coarse-grained models are a popular type of model for accessing long time scales in simulations of biomolecular processes. Such models are coarse-grained with respect to atomic models. But any modelling of processes or substances involves coarse-graining, i.e. the elimination of non-essential degrees of freedom and interactions from a more fine-grained level of modelling. The basic ingredients of developing coarse-grained models based on the properties of fine-grained models are reviewed, together with the conditions that must be satisfied in order to preserve the correct physical mechanisms in the coarse-graining process. This overview should help the reader to determine how realistic a coarse-grained model of a biomolecular system is, i.e. whether it reflects the underlying physical mechanisms or merely provides a set of pretty pictures of the process or substances of interest.

Recent advances in computational actinoid chemistry.

Abstract: We briefly review advances in computational actinoid (An) chemistry during the past ten years in regard to two issues: the geometrical and electronic structures, and reactions. The former addresses the An–O, An–C, and M–An (M is a metal atom including An) bonds in the actinoid molecular systems, including actinoid oxo and oxide species, actinoid–carbenoid, dinuclear and diatomic systems, and the latter the hydration and ligand exchange, the disproportionation, the oxidation, the reduction of uranyl, hydroamination, and the photolysis of uranium azide. Concerning their relevance to the electronic structures and reactions of actinoids and their importance in the development of an advanced nuclear fuel cycle, we also mentioned the work on actinoid carbides and nitrides, which have been proposed to be candidates of the next generation of nuclear fuel, and the oxidation of PuOx, which is important to understand the speciation of actinoids in the environment, followed by a brief discussion on the urgent need for a heavier involvement of computational actinoid chemistry in developing advanced reprocessing protocols of spent nuclear fuel. The paper is concluded with an outlook.

New functionalities in the GROMOS biomolecular simulation software

Abstract: Since the most recent description of the functionalities of the GROMOS software for biomolecular simulation in 2005 many new functions have been implemented. In this article, the new functionalities that involve modified forces in a molecular dynamics (MD) simulation are described: the treatment of electronic polarizability, an implicit surface area and internal volume solvation term to calculate interatomic forces, functions for the GROMOS coarse-grained supramolecular force field, a multiplicative switching function for nonbonded interactions, adiabatic decoupling of a number of degrees of freedom with temperature or force scaling to enhance sampling, and nonequilibrium MD to calculate the dielectric permittivity or viscosity. Examples that illustrate the use of these functionalities are given. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2012

Solvating atomic level fine-grained proteins in supra-molecular level coarse-grained water for molecular dynamics simulations

Abstract: Simulation of the dynamics of a protein in aqueous solution using an atomic model for both the protein and the many water molecules is still computationally extremely demanding considering the time scale of protein motions. The use of supra-atomic or supra-molecular coarse-grained (CG) models may enhance the computational efficiency, but inevitably at the cost of reduced accuracy. Coarse-graining solvent degrees of freedom is likely to yield a favourable balance between reduced accuracy and enhanced computational speed. Here, the use of a supra-molecular coarse-grained water model that largely preserves the thermodynamic and dielectric properties of atomic level fine-grained (FG) water in molecular dynamics simulations of an atomic model for four proteins is investigated. The results of using an FG, a CG, an implicit, or a vacuum solvent environment of the four proteins are compared, and for hen egg-white lysozyme a comparison to NMR data is made. The mixed-grained simulations do not show large differences compared to the FG atomic level simulations, apart from an increased tendency to form hydrogen bonds between long side chains, which is due to the reduced ability of the supra-molecular CG beads that represent five FG water molecules to make solvent-protein hydrogen bonds. But, the mixed-grained simulations are at least an order of magnitude faster than the atomic level ones.

On the effect of a variation of the force field, spatial boundary condition and size of the QM region in QM/MM MD simulations

Abstract: During the past years, the use of combined quantum–classical, QM/MM, methods for the study of complex biomolecular processes, such as enzymatic reactions and photocycles, has increased considerably. The quality of the results obtained from QM/MM calculations is largely dependent on five aspects to be considered when setting up a molecular model: the QM Hamiltonian, the MM Hamiltonian or force field, the boundary and coupling between the QM and MM regions, the size of the QM region and the boundary condition for the MM region. In this study, we systematically investigate the influence of a variation of the molecular mechanics force field and the size of the QM region in QM/MM MD simulations on properties of the photoactive part of the blue light photoreceptor protein AppA. For comparison, we additionally performed classical MD simulations and studied the effect of a variation of the type of spatial boundary condition. The classical boundary conditions and the force field used in a QM/MM MD simulation are shown to have non-neglegible effects upon the structural and energetic properties of the protein which makes it advisable to minimize computational artifacts in QM/MM MD simulations by application of periodic boundary conditions and a thermodynamically calibrated force field. A comparison of the structural and energetic properties of MD simulations starting from two alternative, different X-ray structures for the blue light utilizing flavin protein in its dark state indicates a slight preference of the two force fields used for the so-called Anderson structure over the Jung structure. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2012

Structural Effects of an Atomic-Level Layer of Water Molecules around Proteins Solvated in Supra-Molecular Coarse-Grained Water

Abstract: Atomistic molecular dynamics simulations of proteins in aqueous solution are still limited to the multinanosecond time scale and multinanometer range by computational cost. Combining atomic solutes with a supra-molecular solvent model in hybrid fine-grained/coarse-grained (FG/CG) simulations allows atomic detail in the region of interest while being computationally more efficient. A recent comparison of the properties of four proteins in CG water versus FG water showed the preservation of the secondary and tertiary structure with a computational speed-up of at least an order of magnitude. However, an increased occurrence of hydrogen bonds between side chains was observed due to a lack of hydrogen-bonding partners in the supra-molecular solvent. Here, the introduction of a FG water layer around the protein to recover the hydrogen-bonding pattern of the atomistic simulations is studied. Three layer thicknesses of 0.2, 0.4, and 0.8 nm are considered. A layer thickness of 0.8 nm is found sufficient to recover...

Mixing coarse-grained and fine-grained water in molecular dynamics simulations of a single system

Abstract: The use of a supra-molecular coarse-grained (CG) model for liquid water as solvent in molecular dynamics simulations of biomolecules represented at the fine-grained (FG) atomic level of modelling may reduce the computational effort by one or two orders of magnitude. However, even if the pure FG model and the pure CG model represent the properties of the particular substance of interest rather well, their application in a hybrid FG/CG system containing varying ratios of FG versus CG particles is highly non-trivial, because it requires an appropriate balance between FG-FG, FG-CG, and CG-CG energies, and FG and CG entropies. Here, the properties of liquid water are used to calibrate the FG-CG interactions for the simple-point-charge water model at the FG level and a recently proposed supra-molecular water model at the CG level that represents five water molecules by one CG bead containing two interaction sites. Only two parameters are needed to reproduce different thermodynamic and dielectric properties of liquid water at physiological temperature and pressure for various mole fractions of CG water in FG water. The parametrisation strategy for the FG-CG interactions is simple and can be easily transferred to interactions between atomistic biomolecules and CG water.

Thirty-five years of biomolecular simulation: development of methodology, force fields and software

Abstract: Computer simulation of biomolecular systems has become a standard research instrument for the investigation of biomolecular processes at the atomic level of modelling and interpretation. A bird's e...

Assessment of enveloping distribution sampling to calculate relative free enthalpies of binding for eight netropsin-DNA duplex complexes in aqueous solution.

Abstract: The performance of enveloping distribution sampling (EDS) simulations to estimate free enthalpy differences associated with seven alchemical transformations of A-T into G-C base pairs at the netropsin binding site in the minor groove of a 13-base pair DNA duplex in aqueous solution is evaluated. It is demonstrated that sufficient sampling can be achieved with a two-state EDS Hamiltonian even for large perturbations such as the simultaneous transformation of up to three A-T into three G-C base pairs. The two parameters required to define the EDS reference state Hamiltonian are obtained automatically using a modified version of a scheme presented in earlier work. The sensitivity of the configurational sampling to a variation of these parameters is investigated in detail. Although for relatively small perturbations, that is, one base pair, the free enthalpy estimate depends only weakly on the EDS parameters, the sensitivity is stronger for the largest perturbation. Yet, EDS offers various convenient measures to evaluate the degree of sampling and thus the reliability of the free enthalpy estimate and appears to be an efficient alternative to the conventional thermodynamic integration methodology to obtain free energy differences for molecular systems.

Free enthalpies of replacing water molecules in protein binding pockets

Abstract: Water molecules in the binding pocket of a protein and their role in ligand binding have increasingly raised interest in recent years. Displacement of such water molecules by ligand atoms can be either favourable or unfavourable for ligand binding depending on the change in free enthalpy. In this study, we investigate the displacement of water molecules by an apolar probe in the binding pocket of two proteins, cyclin-dependent kinase 2 and tRNA-guanine transglycosylase, using the method of enveloping distribution sampling (EDS) to obtain free enthalpy differences. In both cases, a ligand core is placed inside the respective pocket and the remaining water molecules are converted to apolar probes, both individually and in pairs. The free enthalpy difference between a water molecule and a CH3 group at the same location in the pocket in comparison to their presence in bulk solution calculated from EDS molecular dynamics simulations corresponds to the binding free enthalpy of CH3 at this location. From the free enthalpy difference and the enthalpy difference, the entropic contribution of the displacement can be obtained too. The overlay of the resulting occupancy volumes of the water molecules with crystal structures of analogous ligands shows qualitative correlation between experimentally measured inhibition constants and the calculated free enthalpy differences. Thus, such an EDS analysis of the water molecules in the binding pocket may give valuable insight for potency optimization in drug design.

On the calculation of 3 J αβ-coupling constants for side chains in proteins

Abstract: Structural knowledge about proteins is mainly derived from values of observables, measurable in NMR spectroscopic or X-ray diffraction experiments, i.e. absorbed or scattered intensities, through theoretically derived relationships between structural quantities such as atom positions or torsional angles on the one hand and observable quantities such as squared structure factor amplitudes, NOE intensities or 3 J-coupling constants on the other. The standardly used relation connecting 3 J-couplings to torsional angles is the Karplus relation, which is used in protein structure refinement as well as in the evaluation of simulated properties of proteins. The accuracy of the simple and generalised Karplus relations is investigated using side-chain structural and 3 J αβ-coupling data for three different proteins, Plastocyanin, Lysozyme, and FKBP, for which such data are available. The results show that the widely used Karplus relations are only a rough estimate for the relation between 3 J αβ-couplings and the corresponding χ1-angle in proteins.

Temperature dependence of the dielectric permittivity of acetic acid, propionic acid and their methyl esters: a molecular dynamics simulation study.

Abstract: For most liquids, the static relative dielectric permittivity is a decreasing function of temperature, because enhanced thermal motion reduces the ability of the molecular dipoles to orient under the effect of an external electric field. Monocarboxylic fatty acids ranging from acetic to octanoic acid represent an exception to this general rule. Close to room temperature, their dielectric permittivity increases slightly with increasing temperature. Herein, the causes for this anomaly are investigated based on molecular dynamics simulations of acetic and propionic acids at different temperatures in the interval 283-363 K, using the GROMOS 53A6(OXY) force field. The corresponding methyl esters are also considered for comparison. The dielectric permittivity is calculated using either the box-dipole fluctuation (BDF) or the external electric field (EEF) methods. The normal and anomalous temperature dependences of the permittivity for the esters and acids, respectively, are reproduced. Furthermore, in the EEF approach, the response of the acids to an applied field of increasing strength is found to present two successive linear regimes before reaching saturation. The low-field permittivity e, comparable to that obtained using the BDF approach, increases with increasing temperature. The higher-field permittivity e' is slightly larger, and decreases with increasing temperature. Further analyses of the simulations in terms of radial distribution functions, hydrogen-bonded structures, and diffusion properties suggest that increasing the temperature or the applied field strength both promote a relative population shift from cyclic (mainly dimeric) to extended (chain-like) hydrogen-bonded structures. The lower effective dipole moment associated with the former structures compared to the latter ones provides an explanation for the peculiar dielectric properties of the two acids compared to their methyl esters.

Implicit Solvation Parameters Derived from Explicit Water Forces in Large-Scale Molecular Dynamics Simulations.

Abstract: Implicit solvation is a mean force approach to model solvent forces acting on a solute molecule. It is frequently used in molecular simulations to reduce the computational cost of solvent treatment. In the first instance, the free energy of solvation and the associated solvent–solute forces can be approximated by a function of the solvent-accessible surface area (SASA) of the solute and differentiated by an atom–specific solvation parameter σiSASA. A procedure for the determination of values for the σiSASA parameters through matching of explicit and implicit solvation forces is proposed. Using the results of Molecular Dynamics simulations of 188 topologically diverse protein structures in water and in implicit solvent, values for the σiSASA parameters for atom types i of the standard amino acids in the GROMOS force field have been determined. A simplified representation based on groups of atom types σgSASA was obtained via partitioning of the atom–type σiSASA distributions by dynamic programming. Three groups of atom types with well separated parameter ranges were obtained, and their performance in implicit versus explicit simulations was assessed. The solvent forces are available at http://mathbio.nimr.mrc.ac.uk/wiki/Solvent_Forces.

Coarse-grained models for the solvents dimethyl sulfoxide, chloroform, and methanol.

Abstract: The time- and length-scale accessible to molecular dynamics simulations of biomolecular systems using atomic-level (AL) models is most limited by the calculation of the solvent-solvent interactions, which comprise the majority of the interactions and yet are seldom of specific interest. Coarse-graining (CG), in which multiple solvent molecules are subsumed into a single bead, provides a means of overcoming this limitation without resorting to implicit solvation models, which basically misrepresent the hydrophobic effect. Most existing CG models, however, do not explicitly include electrostatic interactions, and thus fail to reproduce important properties of the solvent such as dielectric screening. Moreover, CG models for one type of solvent molecule are seldom compatible with those for other solvents. Here, we develop polarizable CG models for the solvents dimethyl sulfoxide, chloroform, and methanol that are compatible with an existing CG model for water. The inclusion of polarizability greatly improves the reproduction of thermodynamic data measured experimentally and calculated from AL simulations for both the pure liquids and binary mixtures.

Interfacing the GROMOS (bio)molecular simulation software to quantum‐chemical program packages

Abstract: The newly implemented quantum-chemical/molecular-mechanical (QM/MM) functionality of the Groningen molecular simulation (GROMOS) software for (bio)molecular simulation is described. The implementation scheme is based on direct coupling of the GROMOS C++ software to executables of the quantum-chemical program packages MNDO and TURBOMOLE, allowing for an independent further development of these packages. The new functions are validated for different test systems using program and model testing techniques. The effect of truncating the QM/MM electrostatic interactions at various QM/MM cutoff radii is discussed and the application of semiempirical versus density-functional Hamiltonians for a solute molecule in aqueous solution is compared. © 2012 Wiley Periodicals, Inc.

Helical Content of a β3‐Octapeptide in Methanol: Molecular Dynamics Simulations Explain a Seeming Discrepancy between Conclusions Derived from CD and NMR Data

Abstract: Connecting experimental observables with the underlying conformational ensemble is a long-standing problem in the structure determination of biomolecules. The simulations described in this article attempt to resolve a seeming discrepancy between the conformational features derived from measured NOE intensities, (3)J-coupling constants, and circular dichroism (CD) spectra for two β-peptides differing in a linker between two side-chains. Although both peptides are very similar in terms of the r(-6) averaged distances between atom pairs involved in the observed NOEs, the molecular dynamics trajectories suggest why the CD spectra show a greater 3(14)-helical propensity for the linked, cyclic peptide than for the linear one, whereas slightly more NMR NOE peaks are observed and assigned for the latter. The nine 100 ns unrestrained simulations show better agreement with the observed experimental data than the single conformations derived from the published NMR structures by additional energy minimization with the GROMOS force field. They show why the seemingly contradictory quantities obtained by NMR and CD spectroscopy can arise from a single conformational ensemble.

Reoptimized interaction parameters for the peptide-backbone model compound N-methylacetamide in the GROMOS force field: influence on the folding properties of two beta-peptides in methanol.

Abstract: Considering N-methylacetamide (NMA) as a model compound, new interaction parameters are developed for the amide function in the GROMOS force field that are compatible with the recently derived 53A6OXY parameter set for oxygen-containing chemical functions. The resulting set, referred to as 53A6OXY+A, represents an improvement over earlier GROMOS force-field versions in the context of the pure-liquid properties of NMA, including the density, heat of vaporization, dielectric permittivity, self-diffusion constant and viscosity, as well as in terms of the Gibbs hydration free energy of this molecule. Assuming that NMA represents an adequate model compound for the backbone of peptides, 53A6OXY+A may be expected to also provide an improved description of polypeptide chains. As an initial test, simulations are reported for two beta-peptides characterized by very different folding properties in methanol. For these systems, earlier force-field versions provided good agreement with experimental NMR data, and the test shows that the improved description achieved in the context of NMA is not accompanied by any deterioration in the representation of the conformational properties of these peptides. (c) 2012 Wiley Periodicals, Inc.

Using enveloping distribution sampling to compute the free enthalpy difference between right- and left-handed helices of a β-peptide in solution

Abstract: Recently, the method of enveloping distribution sampling (EDS) to efficiently obtain free enthalpy differences between different molecular systems from a single simulation has been generalized to compute free enthalpy differences between different conformations of a system [Z. X. Lin, H. Y. Liu, S. Riniker, and W. F. van Gunsteren, J. Chem. Theory Comput. 7, 3884 (2011)]10.1021/ct200623b. However, the efficiency of EDS in this case is hampered if the parts of the conformational space relevant to the two end states or conformations are far apart and the conformational diffusion from one state to the other is slow. This leads to slow convergence of the EDS parameter values and free enthalpy differences. In the present work, we apply the EDS methodology to a challenging case, i.e., to calculate the free enthalpy difference between a right-handed 2.710/12-helix and a left-handed 314-helix of a hexa-β-peptide in solution from a single simulation. No transition between the two helices was detected in a standard...

Validation of the GROMOS 54A7 Force Field Regarding Mixed α/β‐Peptide Molecules

Abstract: A molecular-dynamics (MD) simulation study of two heptapeptides containing α- and β-amino acid residues is presented. According to NMR experiments, the two peptides differ in dominant fold when solvated in MeOH: peptide 3 adopts predominantly β-hairpin-like conformations, while peptide 8 adopts a 14/15-helical fold. The MD simulations largely reproduce the experimental data. Application of NOE atomatom distance restraining improves the agreement with experimental data, but reduces the conformational sampling. Peptide 3 shows a variety of conformations, while still agreeing with the NOE and 3J-coupling data, whereas the conformational ensemble of peptide 8 is dominated by one helical conformation. The results confirm the suitability of the GROMOS 54A7 force field for simulation or structure refinement of mixed α/β-peptides in MeOH.

Influence of 63Ser phosphorylation and dephosphorylation on the structure of the stathmin helical nucleation sequence: a molecular dynamics study.

Abstract: Phosphorylation is an important mechanism regulating protein-protein interactions involving intrinsically disordered protein regions. Stathmin, an archetypical example of an intrinsically disordered protein, is a key regulator of microtubule dynamics in which phosphorylation of 63Ser within the helical nucleation sequence strongly down-regulates the tubulin binding and microtubule destabilizing activities of the protein. Experimental studies on a peptide encompassing the 19-residue helical nucleation sequence of stathmin (residues 55-73) indicate that phosphorylation of 63Ser destabilizes the peptide's secondary structure by disrupting the salt bridges supporting its helical conformation. In order to investigate this hypothesis at atomic resolution, we performed molecular dynamics simulations of nonphosphorylated and phosphorylated stathmin-[55-73] at room temperature and pressure, neutral pH, and explicit solvation using the recently released GROMOS force field 54A7. In the simulations of nonphosphorylated stathmin-[55-73] emerged salt bridges associated with helical configurations. In the simulations of 63Ser phosphorylated stathmin-[55-73] these configurations dispersed and were replaced by a proliferation of salt bridges yielding disordered configurations. The transformation of the salt bridges was accompanied by emergence of numerous interactions between main and side chains, involving notably the oxygen atoms of the phosphorylated 63Ser. The loss of helical structure induced by phosphorylation is reversible, however, as a final simulation showed. The results extend the hypothesis of salt bridge derangement suggested by experimental observations of the stathmin nucleation sequence, providing new insights into regulation of intrinsically disordered protein systems mediated by phosphorylation.

Molecular dynamics simulation of thionated hen egg white lysozyme

Abstract: Understanding of the driving forces of protein folding is a complex challenge because different types of interactions play a varying role. To investigate the role of hydrogen bonding involving the backbone, the effect of thio substitutions in a protein, hen egg white lysozyme (HEWL), was investigated through molecular dynamics simulations of native as well as partly (only residues in loops) and fully thionated HEWL using the GROMOS 54A7 force field. The results of the three simulations show that the structural properties of fully thionated HEWL clearly differ from those of the native protein, while for partly thionated HEWL they only changed slightly compared with native HEWL. The analysis of the torsional-angle distributions and hydrogen bonds in the backbone suggests that the α-helical segments of native HEWL tend to show a propensity to convert to 310-helical geometry in fully thionated HEWL. A comparison of the simulated quantities with experimental NMR data such as nuclear overhauser effect (NOE) atom–atom distance bounds and 3JHNHα-couplings measured for native HEWL illustrates that the information content of these quantities with respect to the structural changes induced by thionation of the protein backbone is rather limited.

Recent Advances in Computational Actinoid Chemistry

Abstract: We briefly review advances in computational actinoid (An) chemistry during the past ten years in regard to two issues: the geometrical and electronic structures, and reactions. The former addresses the An–O, An–C, and M–An (M is a metal atom including An) bonds in the actinoid molecular systems, including actinoid oxo and oxide species, actinoid–carbenoid, dinuclear and diatomic systems, and the latter the hydration and ligand exchange, the disproportionation, the oxidation, the reduction of uranyl, hydroamination, and the photolysis of uranium azide. Concerning their relevance to the electronic structures and reactions of actinoids and their importance in the development of an advanced nuclear fuel cycle, we also mentioned the work on actinoid carbides and nitrides, which have been proposed to be candidates of the next generation of nuclear fuel, and the oxidation of PuOx, which is important to understand the speciation of actinoids in the environment, followed by a brief discussion on the urgent need for a heavier involvement of computational actinoid chemistry in developing advanced reprocessing protocols of spent nuclear fuel. The paper is concluded with an outlook.

An improved structural characterisation of reduced French bean plastocyanin based on NMR data and local-elevation molecular dynamics simulation

Abstract: Deriving structural information about a protein from NMR experimental data is still a non-trivial challenge to computational biochemistry. This is because of the low ratio of the number of independent observables to the number of molecular degrees of freedom, the approximations involved in the different relationships between particular observable quantities and molecular conformation, and the averaged character of the experimental data. For example, protein 3J-coupling data are seldom used for structure refinement because of the multiple-valuedness and limited accuracy of the Karplus relationship linking a 3J-coupling to a torsional angle. Moreover, sampling of the large conformational space is still problematic. Using the 99-residue protein plastocyanin as an example we investigated whether use of a thermodynamically calibrated force field, inclusion of solvent degrees of freedom, and application of adaptive local-elevation sampling that accounts for conformational averaging produces a more realistic representation of the ensemble of protein conformations than standard single-structure refinement in a non-explicit solvent using restraints that do not account for averaging and are partly based on non-observed data. Yielding better agreement with observed experimental data, the protein conformational ensemble is less restricted than when using standard single-structure refinement techniques, which are likely to yield a picture of the protein which is too rigid.

Ester-linked hen egg white lysozyme shows a compact fold in a molecular dynamics simulation - possible causes and sensitivity of experimentally observable quantities to structural changes maintaining this compact fold

Abstract: Prediction and understanding of the folding and stability of the 3D structure of proteins is still a challenge The different atomic interactions, such as non polar contacts and hydrogen bonding, are known but their exact relative weights and roles when contributing to protein folding and stability are not identified Initiated by a previous molecular dynamics simulation of fully ester-linked hen egg white lysozyme (HEWL), which showed a more compact fold of the ester-linked molecule compared to the native one, three variants of this protein are analyzed in the present study These are 129-residue native HEWL, partly ester-linked HEWL, in which only 34 peptide linkages that are not involved in the helical or β-strand parts of native HEWL were replaced by ester linkages, and fully (126 residues) ester-linked HEWL Native and partly ester-linked HEWL showed comparable behaviour, whereas fully ester-linked HEWL could not maintain the native secondary structure of HEWL in the simulation and adopted a more compact fold The conformational changes were analyzed by comparing simulation averaged values of quantities that can be measured by NMR, such as (1)H-(15)N backbone order parameters, residual dipolar couplings, proton-proton NOE distances and (3)J-couplings with the corresponding values derived from experimental NMR data for native HEWL The information content of the latter appeared to be insufficient to detect the local conformational rearrangements upon esterification of the loop regions of the protein For fully ester-linked HEWL, a significantly reduced agreement was observed Upon esterification, the backbone-side chain and side chain-side chain hydrogen-bonding pattern of HEWL changes to maintain its compactness and thus the structural stability of the ester-linked lysozymes

Molecular dynamics simulation of the last step of a catalytic cycle: Product release from the active site of the enzyme chorismate mutase from Mycobacterium tuberculosis

Abstract: The protein chorismate mutase MtCM from Mycobacterium tuberculosis catalyzes one of the few pericyclic reactions known in biology: the transformation of chorismate to prephenate Chorismate mutases have been widely studied experimentally and computationally to elucidate the transition state of the enzyme catalyzed reaction and the origin of the high catalytic rate However, studies about substrate entry and product exit to and from the highly occluded active site of the enzyme have to our knowledge not been performed on this enzyme Crystallographic data suggest a possible substrate entry gate, that involves a slight opening of the enzyme for the substrate to access the active site Using multiple molecular dynamics simulations, we investigate the natural dynamic process of the product exiting from the binding pocket of MtCM We identify a dominant exit pathway, which is in agreement with the gate proposed from the available crystallographic data Helices H2 and H4 move apart from each other which enables the product to exit from the active site Interestingly, in almost all exit trajectories, two residues arginine 72 and arginine 134, which participate in the burying of the active site, are accompanying the product on its exit journey from the catalytic site

Exploring the properties of small molecule protein binding via molecular simulations: the TRSH–p53 core domain complex

Abstract: Molecular dynamics simulations have been performed to investigate the binding of tris(hydroxymethyl)-aminomethane to the surface of the core domain of the mouse cellular tumor antigen p53 employing the GROMOS 45A4 and 53A6 force field parameter sets. A close investigation of the crystal structure reported by Ho et al. revealed that the protonated form is bound to the protein, i.e. a tris(hydroxymethyl)-methylammonium ion (TRSH). Molecular Dynamics (MD) simulations indicate that the p53 protein gains stability upon binding the ligand. In addition to MD simulations of the p53 protein with and without the TRSH compound, thermodynamic integration was utilised to estimate the free enthalpy of binding of the TRSH–p53 complex, which was estimated to be −49 and −54 kJ mol−1 utilising the 45A4 and 53A6 force fields, respectively.

Test of a method for sampling the internal degrees of freedom of a flexible solute molecule based on adiabatic decoupling and temperature or force scaling

Abstract: Simulation of the folding equilibrium of a polypeptide in solution is a computational challenge. Standard molecular dynamics (MD) simulations of such systems cover hundreds of nanoseconds, which is barely sufficient to obtain converged ensemble averages for properties that depend both on folded and unfolded peptide conformations. If one is not interested in dynamical properties of the solute, techniques to enhance the conformational sampling can be used to obtain the equilibrium properties more efficiently. Here the effect on particular equilibrium properties at 298 K of adiabatically decoupling the motion a β-hepta-peptide from the motion of the solvent and subsequently up-scaling its temperature or down-scaling the forces acting on it is investigated. The ensemble averages and rate of convergence are compared to those for standard MD simulations at two different temperatures and a simulation in which the temperature of the solute is increased to 340 K while keeping the solvent at 298 K. Adiabatic decoup...

Calculation of the relative free energy of oxidation of Azurin at pH 5 and pH 9.

Abstract: Free energy calculations are described for the small copper-containing redox protein Azurin from Pseudomonas aeruginosa. A thermodynamic cycle connecting the reduced and oxidized states at pH 5 and pH 9 is considered, allowing for an assessment of convergence in terms of hysteresis and cycle closure. Previously published thermodynamic integration (TI) data is compared to Hamiltonian replica exchange TI (RE-TI) simulations using different simulation setups. The effects of varying simulation length, initial structure, position restraints on particular atoms, and the strength of temperature coupling are studied. Although the overall simulation times are too short to observe an experimentally described peptide plane rotation, it is found that RE-TI simulations do stimulate the distribution of conformational changes over the relevant values of the TI coupling parameter λ. This results in significantly improved values for hysteresis and cycle closure when compared to regular TI. © 2012 Wiley Periodicals, Inc.