A microiterative intrinsic reaction coordinate method for large QM/MM systems.

Area of Science:

• Computational Chemistry
• Theoretical Chemistry
• Chemical Dynamics

Background:

• Intrinsic Reaction Coordinate (IRC) computations are crucial for understanding chemical reaction pathways.
• Standard IRC methods are often computationally prohibitive for large molecular systems, especially those studied with QM/MM.
• A need exists for efficient IRC methods applicable to complex, large-scale chemical systems.

Purpose of the Study:

• To develop and implement an efficient intrinsic reaction coordinate (IRC) computational strategy for large systems.
• To adapt existing IRC algorithms for use with Quantum Mechanics/Molecular Mechanics (QM/MM) methods.
• To demonstrate the broad applicability of the developed method for studying enzymatic reactions.

Main Methods:

• A novel strategy analogous to microiterative transition state optimization is employed.
• IRC equations are applied to a core reaction region, with remaining degrees of freedom relaxed iteratively.
• The implementation integrates with stabilized Euler, local quadratic approximation, and Hessian predictor-corrector IRC algorithms.
• The method was tested using QM/MM setups for enzymatic reactions.

Main Results:

• The developed method successfully performs IRC computations on large systems using QM/MM.
• Validation on small gas-phase systems shows agreement with standard IRC procedures.
• The approach demonstrates broad applicability, successfully applied to two enzymatic reactions.

Conclusions:

• The new IRC strategy effectively extends the application of IRC computations to large molecular systems.
• This development facilitates theoretical studies of complex chemical reactions, including those in biological systems.
• The method offers a valuable tool for computational chemists and biochemists studying reaction mechanisms.