Producing DFT/MM enzyme reaction trajectories from SCC-DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction

Multidimensional Free Energy and Accelerated Quantum Library Methods Provide a Gateway to Glycoenzyme Conformational, Electronic, and Reaction Mechanisms

Enzyme reactions are complex to simulate accurately, and none more so than glycoenzymes (glycosyltransferase and glycosidases). A rigorous sampling of the protein frame and the conformationally plural carbohydrate substrate coupled with an unbiased treatment of the electron dynamics is needed to discover the true reaction landscapes. Here, we demonstrate the effectiveness of two computational methods ported in libraries that we have developed. The first is a flat histogram free energy method called FEARCF capable of multidimensional sampling and rapidly converging to a complete coverage of phase space. The second, the Quantum Supercharger Library (QSL), is a method that accelerates the computation of the ab initio electronic wave function as well as the integral derivatives on graphical processing units (GPUs). These QSL accelerated computations form the core components needed for direct quantum dynamics and QM/MM dynamics when coupled with legacy codes such as GAMESS and NWCHEM, making state of the art hyper-parallel electronic computations in chemistry and chemical biology possible. The combination of QSL (acceleration of ab initio QM computation) and FEARCF (multidimensional hyper-parallel reaction dynamics) makes the simulation of ab initio QM/MM reaction dynamics of enzyme catalysis feasible. Enzymes that process carbohydrates pose an added challenge as their pyranose ring substrates span multidimensional conformational space whose sampling is an intimate function of the catalytic mechanism. Here, we use the pairing of FEARCF and QSL to simulate the catalytic effect of the glycoenzyme β-N-acetylglucosamine transferase (OGT). The reaction mechanism is discovered from a variable three bond reaction surface using SCCDFTB. The role of the OGT in distorting the pyranose ring of β-N-acetylglucosamine (GlcNAc) away from the equilibrium 4C1 chair conformation toward the E3 envelope needed for the transition state is discovered from its pucker free energy hypersurfaces (or free energy volume, FEV). A complete GlcNAc ring pucker HF 6-31g FEV is constructed from ab initio QM dynamics in vacuum and ab initio QM/MM dynamics in the OGT catalytic domain. The OGT is shown to clearly lower the pathway toward the transition state E3 ring conformer as well as stabilize it by 1.63 kcal/mol. Illustrated here is the use of QSL accelerated ab initio QM/MM dynamics that thoroughly explores carbohydrate catalyzed reactions through a FEARCF multidimensional sampling of the interdependence between reaction and conformational space. This demonstrates how experimentally inaccessible molecular and electronic mechanisms that underpin enzyme catalysis can be discovered by directly modeling the dynamics of these complex reactions.

Exploration of the cofactor specificity of wild-type phosphite dehydrogenase and its mutant using molecular dynamics simulations

Phosphite dehydrogenase (Pdh) catalyzes the NAD-dependent oxidation of phosphite to phosphate with the formation of NADH. It can be used in several bioorthogonal systems for metabolic control and related applications, for example, bioelectricity. At present, NAD has poor stability at high concentrations and costs are expensive. Implementation of a non-natural cofactor alternative to the ubiquitous redox cofactor nicotinamide adenosine dinucleotide (NAD) is of great scientific and biotechnological interest. Several Pdhs have been engineered to favor a smaller-sized NAD analogue with a cheaper price and better thermal stability, namely, nicotinamide cytosine dinucleotide (NCD). However, the conformational changes of two cofactors binding to Pdh remain unknown. In this study, five molecular dynamics (MD) simulations were performed to exploit the different cofactors binding to wild-type (WT) Pdh and mutant-type (MT) Pdh (I151R/P176E/M207A). The results were as follows: First, compared with WT Pdh, the cofactor-binding pocket of mutant Pdh became smaller, which may favor a smaller-sized NCD. Second, secondary structure analysis showed that the alpha helices in residues 151–207 partly disappeared in mutant Pdh binding to NAD or NCD. Our theoretical results may provide a basis for further studies on the Pdh family.

Phosphite dehydrogenase (Pdh) catalyzes the NAD-dependent oxidation of phosphite to phosphate with the formation of NADH.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.