Dataset Information

Enhancing Paradynamics for QM/MM Sampling of Enzymatic Reactions.

ABSTRACT: Despite the enormous increase in computer power, it is still extremely challenging to obtain computationally converging sampling of ab initio QM/MM (QM(ai)/MM) free energy surfaces in condensed phases. The sampling problem can be significantly reduced by the use of the reference potential paradynamics (PD) approach, but even this approach still requires major computer time in studies of enzymatic reactions. To further reduce the sampling problem we developed here a new PD version where we use an empirical valence bond reference potential that has a minimum rather than a maximum at the transition state region of the target potential (this is accomplished conveniently by shifting the EVB of the product state). Hence, we can map the TS region in a more efficient way. Here, we introduce and validate the inverted EVB PD approach. The validation involves the study of the S(N)2 step of the reaction catalyzed by haloakene dehalogenase (DhlA) and the GTP hydrolysis in the RasGAP system. In addition, we have also studied the corresponding reaction in water for each of the systems described here and the reaction involving trimethylsulfonium and dimethylamine in solution. The results are encouraging and the new strategy appears to provide a powerful way of evaluating QM(ai)/MM activation free energies.

SUBMITTER: Lameira J

PROVIDER: S-EPMC5380004 | biostudies-literature | 2016 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Enhancing Paradynamics for QM/MM Sampling of Enzymatic Reactions.

Lameira Jerônimo J Kupchencko Ilya I Warshel Arieh A

The journal of physical chemistry. B 20160229 9

Despite the enormous increase in computer power, it is still extremely challenging to obtain computationally converging sampling of ab initio QM/MM (QM(ai)/MM) free energy surfaces in condensed phases. The sampling problem can be significantly reduced by the use of the reference potential paradynamics (PD) approach, but even this approach still requires major computer time in studies of enzymatic reactions. To further reduce the sampling problem we developed here a new PD version where we use an ...[more]

PMID: 26866994

Similar Datasets

Project description:A combined quantum mechanical/molecular mechanical (QM/MM) potential function is used in a thermodynamic integration approach to calculate the pK(a) of residue 66 in two mutants (V66E, V66D) of Staphylococal nuclease relative to solution. Despite the similarity in chemical nature and experimentally measured pK(a) of the two buried titritable residues, the behaviors of the two mutants and the computed pK(a) values vary greatly in the simulations. For Glu66, the side chain is consistently observed to spontaneously flip out from the protein interior during titration, and the overall protein structure remains stable throughout the simulations. The computed pK(a) shifts using conventional sampling techniques with multiple nanoseconds per lambda window (Set A and B) are generally close to the experimental value, therefore indicating that large-scale conformational rearrangements are not as important for V66E as suggested by the recent study of Warshel and co-worker. For Asp66, by contrast, flipping of the shorter side chain is not sufficient for getting adequate solvent stabilization of the ionized state. As a result, more complex behaviors such as partial unfolding of a nearby beta-sheet region is observed, and the computed pK(a) shift is substantially higher than the experimental value unless Asp66 is biased to adopt the similar configurations as Glu66 in the V66E simulations. Collectively, these studies suggest that the lack of electronic polarization is not expected to be the dominant source of error in microscopic pK(a) shift calculations, while the need of enhanced sampling is more compelling for predicting the pK(a) of buried residues. Furthermore, the comparison between V66E and V66D also highlights that the microscopic interpretation of similar apparent pK(a) values and effective "dielectric constants" of proteins can vary greatly in terms of the residues that make key contributions and the scale of structural/hydration response to titration, the latter of which is difficult to predict a priori. Perturbative analyses of interactions that contribute to the titration free energy point to mutants that can be used to verify the microscopic mechanisms of titration in V66E/D SNase proteins.

Dataset Information

Enhancing Paradynamics for QM/MM Sampling of Enzymatic Reactions.

Publications

Enhancing Paradynamics for QM/MM Sampling of Enzymatic Reactions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets