Supra-Atomic Coarse-Grained GROMOS Force Field for Aliphatic Hydrocarbons in the Liquid Phase.

A supra-atomic coarse-grained (CG) force field for liquid n-alkanes is presented. The model was calibrated using experimental thermodynamic data and structural as well as energetic properties for 14 n-alkanes as obtained from atomistic fine-grained (FG) simulations of the corresponding hydrocarbons using the GROMOS 45A3 biomolecular force field. A variation of the nonbonded force-field parameters obtained from mapping the FG interactions onto the CG degrees of freedom to fit the density and heat of vaporization to experimental values turned out to be mandatory for a correct reproduction of these data by the CG model, while the bonded force-field parameters for the CG model could be obtained from a Boltzmann-weighted fit with some variations with respect to the corresponding properties from the FG simulations mapped onto the CG degrees of freedom. The model presents 6 different CG bead types, for bead sizes from 2 to 4 distinguishing between terminal and nonterminal beads within an alkane chain (end or middle). It contains different nonbonded Lennard-Jones parameters for the interaction of CG alkanes with CG water. The CG alkane model was further tested by comparing predictions of the excess free energy, the self-diffusion constant, surface tension, isothermal compressibility, heat capacity, thermal expansion coefficient, and shear viscosity for n-alkanes to experimental values. The CG model offers a thermodynamically calibrated basis for the development of CG models of lipids.

The key to predicting the stability of protein mutants lies in an accurate description and proper configurational sampling of the folded and denatured states.

BACKGROUND: The contribution of particular hydrogen bonds to the stability of a protein fold can be investigated experimentally as well as computationally by the construction of protein mutants which lack particular hydrogen-bond donors or acceptors with a subsequent determination of their structural stability. However, the comparison of experimental data with computational results is not straightforward. One of the difficulties is related to the representation of the unfolded state conformation. METHODS: A series of molecular dynamics simulations of the 34-residue WW domain of protein Pin1 and 20 amide-to-ester mutants started from the X-ray crystal structure and the NMR solution structure are analysed in terms of backbone-backbone hydrogen bonding and differences in free enthalpies of folding in order to provide a structural interpretation of the experimental data available. RESULTS: The contribution of the different ß-sheet hydrogen bonds to the relative stability of the mutants with respect to wild type cannot be directly inferred from experimental thermal denaturation temperatures or free enthalpies of chaotrope denaturation for the different mutants, because some ß-sheet hydrogen bonds show sizeable variation in occurrence between the different mutants. CONCLUSIONS: A proper representation of unfolded state conformations appears to be essential for an adequate description of relative stabilities of protein mutants. GENERAL SIGNIFICANCE: The simulations may be used to link the structural Boltzmann ensembles to relative free enthalpies of folding between mutants and wild-type protein and show that unfolded conformations have to be treated with a sufficient level of detail in free energy calculations of protein stability. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Exploration of swapping enzymatic function between two proteins: a simulation study of chorismate mutase and isochorismate pyruvate lyase.

The enzyme chorismate mutase EcCM from Escherichia coli catalyzes one of the few pericyclic reactions in biology, the transformation of chorismate to prephenate. The isochorismate pyruvate lyase PchB from Pseudomonas aeroginosa catalyzes another pericyclic reaction, the isochorismate to salicylate transformation. Interestingly, PchB possesses weak chorismate mutase activity as well thus being able to catalyze two distinct pericyclic reactions in a single active site. EcCM and PchB possess very similar folds, despite their low sequence identity. Using molecular dynamics simulations of four combinations of the two enzymes (EcCM and PchB) with the two substrates (chorismate and isochorismate) we show that the electrostatic field due to EcCM at atoms of chorismate favors the chorismate to prephenate transition and that, analogously, the electrostatic field due to PchB at atoms of isochorismate favors the isochorismate to salicylate transition. The largest differences between EcCM and PchB in electrostatic field strengths at atoms of the substrates are found to be due to residue side chains at distances between 0.6 and 0.8 nm from particular substrate atoms. Both enzymes tend to bring their non-native substrate in the same conformation as their native substrate. EcCM and to a lower extent PchB fail in influencing the forces on and conformations of the substrate such as to favor the other chemical reaction (isochorismate pyruvate lyase activity for EcCM and chorismate mutase activity for PchB). These observations might explain the difficulty of engineering isochorismate pyruvate lyase activity in EcCM by solely mutating active site residues.

Structure of hen egg-white lysozyme solvated in TFE/water: a molecular dynamics simulation study based on NMR data.

Various experimental studies of hen egg white lysozyme (HEWL) in water and TFE/water clearly indicate structural differences between the native state and TFE state of HEWL, e.g. the helical content of the protein in the TFE state is much higher than in the native state. However, the available detailed NMR studies were not sufficient to determine fully a structure of HEWL in the TFE state. Different molecular dynamics (MD) simulations, i.e. at room temperature, at increased temperature and using proton-proton distance restraints derived from NMR NOE data, have been used to generate configurational ensembles corresponding to the TFE state of HEWL. The configurational ensemble obtained at room temperature using atom-atom distance restraints measured for HEWL in TFE/water solution satisfies the experimental data and has the lowest protein energy. In this ensemble residues 50-58, which are part of the ß-sheet in native HEWL, adopt fluctuating α-helical secondary structure.

Multi-resolution simulation of biomolecular systems: a review of methodological issues.

Theoretical-computational modeling with an eye to explaining experimental observations in regard to a particular chemical phenomenon or process requires choices concerning essential degrees of freedom and types of interactions and the generation of a Boltzmann ensemble or trajectories of configurations. Depending on the degrees of freedom that are essential to the process of interest, for example, electronic or nuclear versus atomic, molecular or supra-molecular, quantum- or classical-mechanical equations of motion are to be used. In multi-resolution simulation, various levels of resolution, for example, electronic, atomic, supra-atomic or supra-molecular, are combined in one model. This allows an enhancement of the computational efficiency, while maintaining sufficient detail with respect to particular degrees of freedom. The basic challenges and choices with respect to multi-resolution modeling are reviewed and as an illustration the differential catalytic properties of two enzymes with similar folds but different substrates with respect to these substrates are explored using multi-resolution simulation at the electronic, atomic and supra-molecular levels of resolution.

Structural effects of an atomic-level layer of water molecules around proteins solvated in supra-molecular coarse-grained water.

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 the behavior of the proteins in the atomistic simulations, whereas the hybrid simulation is still three times more efficient than the atomistic one and the cutoff radius for nonbonded interactions could be increased from 1.4 to 2.0 nm.

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

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 calculation of ³Jαß-coupling constants for side chains in proteins.

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.

Molecular dynamics simulation of thionated hen egg white lysozyme.

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 3(10)-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 (3)J((H)(N)(H)(α))-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.

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.

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.

Biomolecular structure refinement using the GROMOS simulation software.

For the understanding of cellular processes the molecular structure of biomolecules has to be accurately determined. Initial models can be significantly improved by structure refinement techniques. Here, we present the refinement methods and analysis techniques implemented in the GROMOS software for biomolecular simulation. The methodology and some implementation details of the computation of NMR NOE data, (3)J-couplings and residual dipolar couplings, X-ray scattering intensities from crystals and solutions and neutron scattering intensities used in GROMOS is described and refinement strategies and concepts are discussed using example applications. The GROMOS software allows structure refinement combining different types of experimental data with different types of restraining functions, while using a variety of methods to enhance conformational searching and sampling and the thermodynamically calibrated GROMOS force field for biomolecular simulation.

Definition and testing of the GROMOS force-field versions 54A7 and 54B7.

New parameter sets of the GROMOS biomolecular force field, 54A7 and 54B7, are introduced. These parameter sets summarise some previously published force field modifications: The 53A6 helical propensities are corrected through new φ/ψ torsional angle terms and a modification of the N-H, C=O repulsion, a new atom type for a charged -CH(3) in the choline moiety is added, the Na(+) and Cl(-) ions are modified to reproduce the free energy of hydration, and additional improper torsional angle types for free energy calculations involving a chirality change are introduced. The new helical propensity modification is tested using the benchmark proteins hen egg-white lysozyme, fox1 RNA binding domain, chorismate mutase and the GCN4-p1 peptide. The stability of the proteins is improved in comparison with the 53A6 force field, and good agreement with a range of primary experimental data is obtained.

GROMOS++ Software for the Analysis of Biomolecular Simulation Trajectories.

GROMOS++ is a set of C++ programs for pre- and postprocessing of molecular dynamics simulation trajectories and as such is part of the GROningen MOlecular Simulation software for (bio)molecular simulation. It contains more than 70 programs that can be used to prepare data for the production of molecular simulation trajectories and to analyze these. These programs are reviewed and the various structural, dynamic, and thermodynamic quantities that can be analyzed using time series, correlation functions, and distributions are described together with technical aspects of their implementation in GROMOS. A few examples of the use of GROMOS++ for the analysis of MD trajectories are given. A full list of all GROMOS++ programs, together with an indication of their capabilities, is given in the Appendix .