Multiobjective Evolutionary Strategy for Improving Semiempirical Hamiltonians in the Study of Enzymatic Reactions at the QM/MM Level of Theory

Quantum mechanics/molecular mechanics (QM/MM) simulations are crucial for understanding enzymatic reactions, but their accuracy depends heavily on the quantum-mechanical method used. Semiempirical methods offer computational efficiency but often struggle with accuracy in complex systems. This work presents a novel multiobjective evolutionary strategy for optimizing semiempirical Hamiltonians, specifically designed to enhance their performance in enzymatic QM/MM simulations while remaining broadly applicable to condensed-phase systems. Our methodology combines automated parameter optimization, targeting ab initio or density functional theory (DFT)-reference potential energy surfaces, atomic charges, and gradients, with comprehensive validation through minimum free energy path (MFEP) calculations. To demonstrate its effectiveness, we applied our approach to improve the GFN2-xTB Hamiltonian using two enzymatic systems that involve hydride transfer reactions where the activation energy barrier is severely underestimated: Crotonyl-CoA carboxylase/reductase (CCR) and dihydrofolate reductase (DHFR). The optimized parameters showed significant improvements in reproducing potential and free energy surfaces, closely matching higher-level DFT calculations. Through an efficient two-stage optimization process, we first developed parameters for CCR using reaction path data, then refined these parameters for DHFR by incorporating a targeted set of additional training geometries. This strategic approach minimized the computational cost while achieving accurate descriptions of both systems, as validated through QM/MM simulations using the Adaptive String Method (ASM). Our method represents an efficient approach for optimizing semiempirical methods to study larger systems and longer time scales, with potential applications in enzymatic reaction mechanism studies, drug design, and enzyme engineering.

Authors

Institutions

Publication Details

Journal
Journal of Chemical Theory and Computation
Published
2025-05-07
DOI
https://doi.org/10.1021/acs.jctc.5c00247
Citations
1
Primary Topic
Protein Structure and Dynamics
Type
article
Field-Weighted Citation Impact
0.81

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Multiobjective Evolutionary Strategy for Improving Semiempirical Hamiltonians in the Study of Enzymatic Reactions at the QM/MM Level of Theory

José Luís Velázquez‐Libera, Rodrigo Recabarren, Yamisleydi Salgueiro, J. Javier Ruiz‐Pernía et al.
1 citations
Journal of Chemical Theory and Computation
Protein Structure and Dynamics
0.81
article

Multiobjective Evolutionary Strategy for Improving Semiempirical Hamiltonians in the Study of Enzymatic Reactions at the QM/MM Level of Theory

José Luís Velázquez‐Libera, Rodrigo Recabarren, Yamisleydi Salgueiro, J. Javier Ruiz‐Pernía, Esteban Vöhringer‐Martinez, Iñaki Tuñón, Julio Caballero
article en
1 citations

Abstract

Quantum mechanics/molecular mechanics (QM/MM) simulations are crucial for understanding enzymatic reactions, but their accuracy depends heavily on the quantum-mechanical method used. Semiempirical methods offer computational efficiency but often struggle with accuracy in complex systems. This work presents a novel multiobjective evolutionary strategy for optimizing semiempirical Hamiltonians, specifically designed to enhance their performance in enzymatic QM/MM simulations while remaining broadly applicable to condensed-phase systems. Our methodology combines automated parameter optimization, targeting ab initio or density functional theory (DFT)-reference potential energy surfaces, atomic charges, and gradients, with comprehensive validation through minimum free energy path (MFEP) calculations. To demonstrate its effectiveness, we applied our approach to improve the GFN2-xTB Hamiltonian using two enzymatic systems that involve hydride transfer reactions where the activation energy barrier is severely underestimated: Crotonyl-CoA carboxylase/reductase (CCR) and dihydrofolate reductase (DHFR). The optimized parameters showed significant improvements in reproducing potential and free energy surfaces, closely matching higher-level DFT calculations. Through an efficient two-stage optimization process, we first developed parameters for CCR using reaction path data, then refined these parameters for DHFR by incorporating a targeted set of additional training geometries. This strategic approach minimized the computational cost while achieving accurate descriptions of both systems, as validated through QM/MM simulations using the Adaptive String Method (ASM). Our method represents an efficient approach for optimizing semiempirical methods to study larger systems and longer time scales, with potential applications in enzymatic reaction mechanism studies, drug design, and enzyme engineering.

Journal of Chemical Theory and ComputationVol. 21(10)
University of Talca (CL), Universitat de València (ES), University of Concepción (CL), Institute of Industrial Engineering (AU)
Ministerio de Ciencia, Innovación y Universidades, Max-Planck-Gesellschaft, European Regional Development Fund, Fondo Nacional de Desarrollo Científico y Tecnológico
Zero hunger
Openalex Percentile: Top 28%
Protein Structure and Dynamics
0.81
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.