Reduction pathway of glutaredoxin 1 investigated with QM/MM molecular dynamics using a neural network correction

Abstract

Glutaredoxins are small enzymes that catalyze the oxidation and reduction of protein disulfide bonds by the thiol-disulfide exchange mechanism. They have either one or two cysteines in their active site, resulting in different catalytic reaction cycles that have been investigated in many experimental studies. However, the exact mechanisms are not yet fully known, and to our knowledge, no theoretical studies have been performed to elucidate the underlying mechanism. In this study, we investigated a proposed mechanism for the reduction of the disulfide bond in the protein HMA4n by a mutated monothiol Homo sapiens glutaredoxin (HsGrx1) and the co-substrate glutathione (GSH). The catalytic cycle involves three successive thiol-disulfide exchanges that occur between the molecules. To estimate the regioselectivity of the different attacks,classical molecular dynamics simulations were performed and the trajectories analyzed regarding the sulfur--sulfur distances and the attack angles between the sulfurs. The free energy profile of each reaction was obtained with hybrid quantum mechanical/molecular mechanical metadynamics simulations. Since this required extensive phase space sampling, the semi-empirical density functional tight-binding (DFTB) method was used to describe the reactive cysteines. For an accurate description, we used specific reaction parameters fitted to B3LYP energies of the thiol-disulfide exchange and a machine learned energy correction that was trained on CCSD(T) energies of thiol-disulfide exchanges. Our calculations show the same regiospecifity as observed in the experiment and the obtained barrier heights are about 12 and 20~kcal/mol for the different reaction steps, which confirms the proposed pathway.

Modern Alchemical Free Energy Methods for Drug Discovery Explained

“…Emerging technologies that promise to improve the accuracy and precision of AFE calculations in drug discovery include the development of end-state ensemble reservoirs that can be used within networks to ensure consistent end states as well as continued evolution of sampling methods. Finally, there are a number of exciting advances in the development of new force fields that promise higher levels of accuracy, most notably classical polarizable force fields − and machine learning potentials (MLPs). − Of particular promise are the new methods that combine fast, approximate quantum-mechanical (QM) models with MLP corrections to achieve high accuracy (QM/Δ-MLPs). ,− These methods are “universal” in the sense that unlike molecular mechanical force fields (including polarizable force fields), they do not assume a predetermined bonding topology and are able to accurately model different tautomers and protonation states. This is highly significant, as 30% of the compounds in vendor databases and 21% drug databases have potential tautomers , and it has been estimated that up to 95% of drug molecules contain ionizable groups …”

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

“…In the present work, we develop a Quantum Deep-learning Potential Interaction (QDπ) model that uses a fast third-order self-consistent density-functional tight-binding (DFTB3/3OB) model , that is corrected to a quantitatively high-level of accuracy through a range-corrected deep-learning potential (DPRc). , In this way, the QDπ model developed here is the form of a quantum mechanical/machine learning potential correction (QM/Δ-MLP). ,− The use of DFTB3 as a robust QM base model has several important advantages. First, it provides a reasonable description of the conformational potential energy landscape, greatly reducing the requirement to explicitly train the MLP to avoid inaccessible high-energy regions.…”