MoD-QM/MM Structural Refinement Method: Characterization of Hydrogen Bonding in the Oxytricha nova G-Quadruplex

Area of Science:

• Computational Chemistry
• Biophysics
• Molecular Modeling

Background:

• Accurate structural and electronic descriptions of extended systems are crucial in computational chemistry.
• Existing hybrid QM/MM methods often face limitations in self-consistency and accuracy for complex systems.
• Refining geometries and electrostatic potentials is essential for reliable molecular simulations.

Purpose of the Study:

• Introduce a generalized Moving-Domain Quantum Mechanics/Molecular Mechanics (MoD-QM/MM) method for self-consistent structural refinement.
• Enhance the ab initio quality description of extended systems' geometry and molecular electrostatic potential.
• Improve the accuracy of NMR chemical shift predictions for complex biological molecules.

Main Methods:

• Developed a generalized MoD-QM/MM protocol for iterative optimization of QM layers within a molecular environment.
• Applied the method to model systems for full QM optimization benchmarking.
• Refined the hydrogen bonding geometry of Oxytricha nova guanine quadruplex using the MoD-QM/MM approach.
• Performed (1)H NMR chemical shift calculations using the gauge-independent atomic orbital (GIAO) method.

Main Results:

• The generalized MoD-QM/MM method provides a self-consistent protocol for structural refinement of extended systems.
• Benchmarking on model systems confirmed the method's ability to achieve ab initio quality descriptions.
• Refinement of the guanine quadruplex hydrogen bonding geometry showed significant improvements.
• Solvated MoD-QM/MM structures yielded NMR simulations in much better agreement with experimental data compared to X-ray or classical MM models.

Conclusions:

• The generalized MoD-QM/MM method offers a robust computational protocol for accurate structural and electronic characterization of extended systems.
• Explicit solvent molecular dynamics simulations combined with MoD-QM/MM significantly enhance NMR prediction accuracy.
• This approach represents a substantial advancement over traditional methods relying solely on X-ray structures or classical force fields for NMR analysis.