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

Abstract

The SCC-DFTB/MIO/CHARMM free energy surface for a glycosyltransferase, TcTS, is benchmarked against a DFT/MM reaction trajectory using the same CHARMM MM force field ported to the NWChem package. The popular B3LYP functional, against which the MIO parameter set was parameterized is used to optimize TS structures and run DFT reaction dynamics. A novel approach was used to generate reaction forces from a SCC-DFTB/MIO/CHARMM reaction surface to drive B3LYP/6-31G/MM and B3LYP/6-31G(d)/MM reaction trajectories. Although TS structures compare favorably, differences stemming primarily from a minimal basis set approximation prevented a successful 6-31G(d) FEARCF reaction dynamics trajectory. None the less, the dynamic evolution of the B3LYP/6-31G/MM-computed electron density provided an opportunity to perform NBO analysis along the reaction trajectory. Here, we illustrate that a successful ab initio reaction trajectory is computationally accessible when the underlying potential energy function of the semi-empirical method used to produce driving forces is sufficiently close to the ab initio potential. © 2017 Wiley Periodicals, Inc.

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

“…All residues were assigned in their standard protonation (pH = 7), and all missing hydrogen atoms were generated using Discovery Studio 4.0 client software. NCD and NAD were then optimized according to the density functional theory at the B3LYP/6-31G* level 23 by using Gaussian 09 software. 24 …”

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

“…The implementation of the DFTB method in codes widely used in hybrid DFT/MM calculations has considerably facilitated access to this method for such hybrid studies. Very different enzymatic mechanisms have been explored, such as proton transfer reactions or proton storage [241,242], histone methylation [243], C-terminal residue cleavage [244], amide hydrolysis [245], glycosylation/deglycosylation [246,247], inactivation of a new tuberculosis target [248], hydrolysis of organophosphorus [51] or proton-coupled electron transfer reactions [249]. One can also cite DFTB studies aimed at investigating substrate promiscuity [250], ion binding and transport by membrane proteins [251], proton distribution over multiple binding sites of a membrane protein [252] or evaluating the pKa of protein residues [253].…”

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

“…Producing DFT/MM enzyme reaction trajectories from SCC-DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction. J Comput Chem 2017 , 38, 1789–1798 …”

“…The vast and complex conformational space that make glycans important functional units in cellular biology presents a computational modeling challenge to sufficiently sample hemiacetal phase space. ,, While the cellular processing of glycans via glycoenzymes (glycosyltransferases and glycosidases) are known drug targets of a wide range of communicable and noncommunicable diseases, unpacking the mechanistic details and designing new drugs requires commensurate ab initio level QM/MM tools to accurately simulate glycoenzyme catalytic action . The epimeric choices leading to variation in cyclic pyranose monomers, the ability of the pyranose rings to pucker into 38 distinct conformations, and the 5 or more functional groups bonded to the ring carbons that can rotate into distinct orientations (cis, trans, gauche) makes carbohydrates ideal for data storage and information transfer.…”

“…Using the comprehensive sampling of reaction space, the competition between substitution and elimination pathways in the Trypanosoma cruzi trans-sialidase (TcTS) catalyzed reaction and TcTS’s suppression of its side reactions to yield a single product were discovered through a CHARMM implementation of FEARCF to compute a SCCDFTB/MM free energy surface W (ξ) or more specifically here, a reaction surface . The underlying electronic dynamics were discovered using a NWChem implementation of FEARCF to compute the B3LYP/6-31G/MM dynamics employing driving forces from the SCC-DFTB/MM W (ξ) for the TcTS glycosylation reaction …”