PubMed | Spain; Biomol Informatics SL, Autonomous University of Madrid and Molecular Modelling Group
Type: | Journal: Advances in protein chemistry and structural biology | Year: 2015
Quantum mechanics/molecular mechanics (QM/MM) methods are excellent tools for the modeling of biomolecular reactions. Recently, we have implemented a new QM/MM method (Fireball/Amber), which combines an efficient density functional theory method (Fireball) and a well-recognized molecular dynamics package (Amber), offering an excellent balance between accuracy and sampling capabilities. Here, we present a detailed explanation of the Fireball method and Fireball/Amber implementation. We also discuss how this tool can be used to analyze reactions in biomolecules using steered molecular dynamics simulations. The potential of this approach is shown by the analysis of a reaction catalyzed by the enzyme triose-phosphate isomerase (TIM). The conformational space and energetic landscape for this reaction are analyzed without a priori assumptions about the protonation states of the different residues during the reaction. The results offer a detailed description of the reaction and reveal some new features of the catalytic mechanism. In particular, we find a new reaction mechanism that is characterized by the intramolecular proton transfer from O1 to O2 and the simultaneous proton transfer from Glu 165 to C2.
PubMed | Spain; Biomol Informatics SL, Copenhagen University and Molecular Modelling Group
Type: Journal Article | Journal: PloS one | Year: 2015
The continued development and utility of molecular dynamics simulations requires improvements in both the physical models used (force fields) and in our ability to sample the Boltzmann distribution of these models. Recent developments in both areas have made available multi-microsecond simulations of two proteins, ubiquitin and Protein G, using a number of different force fields. Although these force fields mostly share a common mathematical form, they differ in their parameters and in the philosophy by which these were derived, and previous analyses showed varying levels of agreement with experimental NMR data. To complement the comparison to experiments, we have performed a structural analysis of and comparison between these simulations, thereby providing insight into the relationship between force-field parameterization, the resulting ensemble of conformations and the agreement with experiments. In particular, our results show that, at a coarse level, many of the motional properties are preserved across several, though not all, force fields. At a finer level of detail, however, there are distinct differences in both the structure and dynamics of the two proteins, which can, together with comparison with experimental data, help to select force fields for simulations of proteins. A noteworthy observation is that force fields that have been reparameterized and improved to provide a more accurate energetic description of the balance between helical and coil structures are difficult to distinguish from their unbalanced counterparts in these simulations. This observation implies that simulations of stable, folded proteins, even those reaching 10 microseconds in length, may provide relatively little information that can be used to modify torsion parameters to achieve an accurate balance between different secondary structural elements.