Atomistic Simulation of Adiabatic Reactive Processes Based on Multi-State Potential Energy Surfaces

Area of Science:

• Computational chemistry
• Biophysics
• Chemical kinetics

Background:

• Molecular dynamics (MD) simulations are crucial for studying chemical reactions at the atomic level.
• Existing MD methods can be computationally intensive for reactive processes.
• The adiabatic reactive molecular dynamics (ARMD) method offers a computationally efficient framework.

Purpose of the Study:

• To generalize the ARMD method for arbitrary reactive processes.
• To apply the generalized ARMD to study NO rebinding in myoglobin and conformational transitions in neuroglobin.
• To validate the ARMD method's accuracy and explore its capabilities.

Main Methods:

• Generalization of the ARMD method to handle arbitrary two-state reactive processes.
• Treatment of reactants and products using an atomistic force field.
• Application to myoglobin-NO (MbNO) kinetics and neuroglobin (Ngb) conformational changes.

Main Results:

• The generalized ARMD method successfully reproduced the nonexponential kinetics of NO rebinding to myoglobin (MbNO).
• The study investigated the influence of potential energy surface separation on MbNO kinetics.
• The method's functionality was explored through application to neuroglobin conformational transitions.

Conclusions:

• The generalized ARMD method is a versatile and accurate tool for studying complex chemical reactions and conformational dynamics.
• ARMD provides a computationally efficient approach for simulating biomolecular systems.
• The findings validate ARMD for kinetic studies and conformational analysis in proteins.