| 1. |
- Oanca, Gabriel, et al.
(författare)
-
Efficient Empirical Valence Bond Simulations with GROMACS
- 2023
-
Ingår i: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 19:17, s. 6037-6045
-
Tidskriftsartikel (refereegranskat)abstract
- We describe a protocol to perform empirical valence bond (EVB) simulations using GROMACS software. EVB is a fast and reliable method that allows one to determine the reaction free-energy profiles in complex systems, such as enzymes, by employing classical force fields to represent a chemical reaction. Therefore, running EVB simulations is basically as fast as any classical molecular dynamics simulation, and the method uses standard free-energy calculations to map the free-energy change along a given reaction path. To exemplify and validate our EVB implementation, we replicated two cases of our earlier enzyme simulations. One of these addresses the decomposition of the activation free energy into its enthalpic and entropic components, and the other is focused on calculating the overall catalytic effect of the enzyme compared to the same reaction in water. These two examples give virtually identical results to those obtained with programs that were specifically designed for EVB simulations and show that the GROMACS implementation is robust and can be used for very large systems.
|
|
| 2. |
- Oanca, Gabriel, et al.
(författare)
-
Insights into enzyme point mutation effect by molecular simulation : phenylethylamine oxidation catalyzed by monoamine oxidase A
- 2016
-
Ingår i: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 18:19, s. 13346-13356
-
Tidskriftsartikel (refereegranskat)abstract
- The I335Y point mutation effect on the kinetics of phenylethylamine decomposition catalyzed by monoamine oxidase A was elucidated by means of molecular simulation. The established empirical valence bond methodology was used in conjunction with the free energy perturbation sampling technique and a classical force field representing the state of reactants and products. The methodology allows for the simulation of chemical reactions, in the present case the breaking of the alpha-C-H bond in a phenylethylamine substrate and the subsequent hydrogen transfer to the flavin cofactor, resulting in the formation of the N-H bond on flavin. The empirical parameters were calibrated against the experimental data for the simulated reaction in a wild type protein and then used for the calculation of the reaction free energy profile in the I335Y mutant. In very good agreement with the measured kinetic data, mutation increases the free energy barrier for the rate limiting step by slightly more than 1 kcal mol(-1) and consequently decreases the rate constant by about an order of magnitude. The magnitude of the computed effect slightly varies with simulation settings, but always remains in reasonable agreement with the experiment. Analysis of trajectories reveals a major change in the interaction between phenyl rings of the substrate and the neighboring Phe352 residue upon the I335Y mutation due to the increased local polarity, leading to an attenuated quadrupole interaction between the rings and destabilization of the transition state. Additionally, the increased local polarity in the mutant allows for a larger number of water molecules to be present near the active site, effectively shielding the catalytic effect of the enzyme and contributing to the increased barrier.
|
|
| 3. |
- Oanca, Gabriel, et al.
(författare)
-
Why Do Empirical Valence Bond Simulations Yield Accurate Arrhenius Plots?
- 2024
-
Ingår i: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 20:6, s. 2582-2591
-
Tidskriftsartikel (refereegranskat)abstract
- Computer simulations of the temperature dependence of enzyme reactions using the empirical valence bond (EVB) method have proven to give very accurate results in terms of the thermodynamic activation parameters. Here, we analyze the reasons for why such simulations are able to correctly capture activation enthalpies and entropies and how sensitive these quantities are to parametrization of the reactive potential energy function. We examine first the solution reference reaction for the enzyme ketosteroid isomerase, which corresponds to the acetate catalyzed deprotonation of the steroid in water. The experimentally determined activation parameters for this reaction turn out to be remarkably well reproduced by the calculations. By modifying the EVB potential so that the activation and reaction free energies become significantly shifted, we show that the activation entropy is basically invariant to such changes and that ΔS⧧ is instead determined by the specific mixture of the underlying force fields in the transition state region. The coefficients of this mixture do not change appreciably when the EVB potential is modified within reasonable limits, and hence, the estimate of ΔS⧧ becomes very robust. This is further verified by examining a more complex concerted hydride and proton transfer reaction in the enzyme hydroxybutyrate dehydrogenase.
|
|
