Diphtheria toxin catalyzed hydrolysis of NAD(+): molecular dynamics study of enzyme-bound substrate, transition state, and inhibitor.

The mechanism of the diphtheria toxin-catalyzed hydrolysis of NAD(+) was investigated by quantum chemical calculations and molecular dynamics simulations. Several effects that could explain the 6000-fold rate acceleration (Delta Delta G(++) approximately 5 kcal/mol) by the enzyme were considered. First, the carboxamide arm of the enzyme-bound NAD(+) adopts a trans conformation while the most stable conformation is cis. The most stable conformation for the nicotinamide product has the amide carbonyl trans. The activation energy for the cleavage of the ribosidic bond is reduced by 2 kcal/mol due to the relaxation of this ground state conformational stress in the transition state. Second, molecular dynamics simulations to the nanosecond time range revealed that the carboxylate of Glu148 forms a hydrogen bond to the substrate's 2' hydroxyl group in E.S (approximately 17% of the time) and E.TS (approximately 57% of the time) complexes. This interaction is not seen in crystal structures. The ApUp inhibitor is held more tightly by the enzyme than the transition state and the substrate. Analysis of correlated motions reveals differences in the pattern of anticorrelated motions for protein backbone atoms when the transition state occupies the active site as compared to the E.NAD(+) complex.     

In silico model for P-glycoprotein substrate prediction: insights from molecular dynamics and in vitro studies

P-glycoprotein (P-gp) is a plasma membrane efflux transporter belonging to ATP-binding cassette superfamily, responsible for multidrug resistance in tumor cells. Over-expression of P-gp in cancer cells limits the efficacy of many anticancer drugs. A clear understanding of P-gp substrate binding will be advantageous in early drug discovery process. However, substrate poly-specificity of P-gp is a limiting factor in rational drug design. In this investigation, we report a dynamic trans-membrane model of P-gp that accurately identified the substrate binding residues of known anticancer agents. The study included homology modeling of human P-gp based on the crystal structure of C. elegans P-gp, molecular docking, molecular dynamics analyses and binding free energy calculations. The model was further utilized to speculate substrate propensity of in-house anticancer compounds. The model demonstrated promising results with one anticancer compound (NSC745689). As per our observations, the molecule could be a potential lead for anticancer agents devoid of P-gp mediated multiple drug resistance. The in silico results were further validated experimentally using Caco-2 cell lines studies, where NSC745689 exhibited poor permeability (P _app 1.03 ± 0.16 × 10^?6 cm/s) and low efflux ratio of 0.26.     

Ab initio molecular dynamics calculations to study catalysis

The modern versions of the density functional theory (DFT), especially those using the generalized gradient approximation (GGA), have reached (almost) chemical accuracy and thus can be applied to study problems of real chemical interest such as catalysis. The important equations for the DFT, the local density approximation (LDA), and GGA are given. The full-potential linearized augmented plane wave method (LAPW) is used to check the accuracy of GGA in solids. The basic concepts of the ab initio molecular dynamics (MD) method by Car and Parrinello and its recent implementation using the projector augmented Wave (PAW) method which use a similar augmentation as LAPW are described. PAW applications to ferrocene and beryllocene are summarized, which illustrate that vibrational or fluxional behavior are well described. Sodalite, an aluminosilicate, is discussed as a generic zeolite in comparison with gmelinite. A study of the dynamics of such a system allows the determination of, e.g., the proton stretch vibrations which can be related to infrared spectra. This is illustrated for the OH stretch vibration of the acid site in silicon-rich sodalite. With this methodology, we are able to study the interaction of methanol trapped inside the cage structure of silicon-rich sodalite and to gain new insight into crucial steps of catalytic reactions, namely, the hydrogen-bonding and the possible protonation in this system, or a proton-exchange reaction. The strategies for parallelizing the PAW code are outlined. © 1997 John Wiley & Sons, Inc.     

Exploring phase-transfer catalysis with molecular dynamics and 3D/4D quantitative structure-selectivity relationships

Quantitative Structure-Selectivity Relationships (QSSR) are developed for a library of 40 phase-transfer asymmetric catalysts, based around quaternary ammonium salts, using Comparative Molecular Field Analysis (CoMFA) and closely related variants. Due to the flexibility of these catalysts, we use molecular dynamics (MD) with an implicit Generalized Born solvent model to explore their conformational space. Comparison with crystal data indicates that relevant conformations are obtained and that, furthermore, the correct biphenyl twist conformation is predicted, as illustrated by the superiority of the resulting model (leave-one-out q(2) = 0.78) compared to a random choice of low-energy conformations for each catalyst (average q(2) = 0.22). We extend this model by incorporating the MD trajectory directly into a 4D QSSR and by Boltzmann-weighting the contribution of selected minimized conformations, which we refer to as '3.5D' QSSR. The latter method improves on the predictive ability of the 3D QSSR (leave-one-out q(2) = 0.83), as confirmed by repeated training/test splits.     

A molecular dynamics simulation study of peptide deformylase from Leptospira interrogans complex：Exploring the closing mechanism of the substrate pocket

To explore the closing mechanism of the substrate pocket,we perform a 16,000 ps molecular dynamics simulation separately on the ligand-free and actinonin-bound peptide deformylase from Leptospira interrogans.Our results show that the CD-loop, hydrophilic inhibitor and hydrophobic cluster are necessary for the formation of semi-open conformation,and Tyr71 plays an important role in mediating the movements of CD-loop.The average MD structure of the actinonin-bound LiPDF complex approaches to the crystal structure.These are consistent with experiment very well.     

Structural explanation for the tunable substrate specificity of an <Emphasis Type="Italic">E. coli</Emphasis> nucleoside hydrolase: insights from molecular dynamics simulations

Parasitic protozoa rely on nucleoside hydrolases that play key roles in the purine salvage pathway by catalyzing the hydrolytic cleavage of the N-glycosidic bond that connects nucleobases to ribose sugars. Cytidine–uridine nucleoside hydrolase (CU–NH) is generally specific toward pyrimidine nucleosides; however, previous work has shown that replacing two active site residues with Tyr, specifically the Thr223Tyr and Gln227Tyr mutations, allows CU–NH to process inosine. The current study uses molecular dynamics (MD) simulations to gain atomic-level insight into the activity of wild-type and mutant E. coli CU–NH toward inosine. By examining systems that differ in the identity and protonation states of active site catalytic residues, key enzyme-substrate interactions that dictate the substrate specificity of CU–NH are identified. Regardless of the wild-type or mutant CU–NH considered, our calculations suggest that inosine binding is facilitated by interactions of the ribose moiety with active site residues and Ca²⁺, and π-interactions between two His residues (His82 and His239) and the nucleobase. However, the lack of observed activity toward inosine for wild-type CU–NH is explained by no residue being correctly aligned to stabilize the departing nucleobase. In contrast, a hydrogen-bonding network between hypoxanthine and a newly identified general acid (Asp15) is present when the two Tyr mutations are engineered into the active site. Investigation of the single CU–NH mutants reveals that this hydrogen-bonding network is only maintained when both Tyr mutations are present due to a π-interaction between the residues. These results rationalize previous experiments that show the single Tyr mutants are unable to efficiently hydrolyze inosine and explain how the Tyr residues work synergistically in the double mutant to stabilize the nucleobase leaving group during hydrolysis. Overall, our simulations provide a structural explanation for the substrate specificity of nucleoside hydrolases, which may be used to rationally develop new treatments for kinetoplastid diseases.

Graphical Abstract

Open image in new window

     

Ground state structures and excited state dynamics of pyrrole-water complexes: ab initio excited state molecular dynamics simulations

Structures of the ground state pyrrole-(H2O)n clusters are investigated using ab initio calculations. The charge-transfer driven femtosecond scale dynamics are studied with excited state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field method for pyrrole-(H2O)n clusters. Upon the excitation of these clusters, the charge density is located over the farthest water molecule which is repelled by the depleted pi-electron cloud of pyrrole ring, resulting in a highly polarized complex. For pyrrole-(H2O), the charge transfer is maximized (up to 0.34 a.u.) around approximately 100 fs and then oscillates. For pyrrole-(H2O)2, the initial charge transfer occurs through the space between the pyrrole and the pi H-bonded water molecule and then the charge transfer takes place from this water molecule to the sigma H-bonded water molecule. The total charge transfer from the pyrrole to the water molecules is maximized (up to 0.53 a.u.) around approximately 100 fs.     

Promoting vibrations in human purine nucleoside phosphorylase. A molecular dynamics and hybrid quantum mechanical/molecular mechanical study

Crystallographic studies of human purine nucleoside phosphorylase (hPNP) with several transition-state (TS) analogues in the immucillin family showed an unusual geometric arrangement of the atoms O-5', O-4', and O(P), the nucleophilic phosphate oxygen, lying in a close three-oxygen stack. These observations were corroborated by extensive experimental kinetic isotope effect analysis. We propose that protein-facilitated dynamic modes in hPNP cause this stack, centered on the ribosyl O-4' oxygen, to squeeze together and push electrons toward the purine ring, stabilizing the oxacarbenium character of the TS. As the N-ribosidic bond is cleaved during the reaction, the pK(a) values of N-7 and O-6 increase by the electron density expelled by the oxygen-stack compression toward the purine ring. Increased electron density in the purine ring improves electrostatic interactions with nearby residues and facilitates the abstraction of a proton from a solvent proton or an unidentified general acid, making the purine a better leaving group, and accelerating catalysis. Classical and mixed quantum/classical molecular dynamics (MD) simulations of the Michaelis complex of hPNP with the substrates guanosine and phosphate were performed to assess the existence of protein-promoting vibrations (PPVs). Analogous simulations were performed for the substrates in aqueous solution. In the catalytic site, the O-5', O-4', and O(P) oxygens vibrate at frequencies of ca. 125 and 465 cm(-1), as opposed to 285 cm(-1) in the absence of hPNP. The hybrid quantum mechanical/molecular mechanical method was used to assess whether this enzymatic vibration pushing the oxygens together is coupled to the reaction coordinate, and thus has a direct positive impact on catalysis. The potential energy surface for the phosphorolysis reaction for several snapshots taken from the classical MD simulation showed substantial differences in oxygen compression. Our calculations showed the existence of PPVs coupled to the reaction coordinate, which effect electronic alterations in the active site by pushing the three oxygen centers together in proximity, and accelerate substrate turnover in the phosphorolysis reaction catalyzed by hPNP.     

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.     

Glycoside hydrolase stabilization of transition state charge: new directions for inhibitor design

Carbasugars are structural mimics of naturally occurring carbohydrates that can interact with and inhibit enzymes involved in carbohydrate processing. In particular, carbasugars have attracted attention as inhibitors of glycoside hydrolases (GHs) and as therapeutic leads in several disease areas. However, it is unclear how the carbasugars are recognized and processed by GHs. Here, we report the synthesis of three carbasugar isotopologues and provide a detailed transition state (TS) analysis for the formation of the initial GH-carbasugar covalent intermediate, as well as for hydrolysis of this intermediate, using a combination of experimentally measured kinetic isotope effects and hybrid QM/MM calculations. We find that the α-galactosidase from Thermotoga maritima effectively stabilizes TS charge development on a remote C5-allylic center acting in concert with the reacting carbasugar, and catalysis proceeds via an exploded, or loose, S_N2 transition state with no discrete enzyme-bound cationic intermediate. We conclude that, in complement to what we know about the TS structures of enzyme-natural substrate complexes, knowledge of the TS structures of enzymes reacting with non-natural carbasugar substrates shows that GHs can stabilize a wider range of positively charged TS structures than previously thought. Furthermore, this enhanced understanding will enable the design of new carbasugar GH transition state analogues to be used as, for example, chemical biology tools and pharmaceutical lead compounds.

Positive charge stabilized on remote C5-allylic center with catalysis occurring via a loose S_N2 transition state.     

Computer simulations of trypanosomal nucleoside hydrolase: determination of the protonation state of the bound transition-state analogue

Inosine-uridine nucleoside hydrolase (IU-NH) catalyzes the hydrolysis of nucleosides into base and ribose moieties via a ribooxocarbenium ion transition state, which has been characterized using kinetic isotope effects. Protozoan parasites lack de novo purine and pyrimidine biosynthesis and depend on the purine salvage from the host. Vern Schramm and co-workers characterized p-aminophenyliminoribitol (pAPIR) to be a potent inhibitor of IU-NH from Crithidia fasciculata with K(d) of 30 nM. The cyclic amine function of the iminoribitol ring can be either protonated (pAPIRH(+)) or unprotonated (pAPIR). pAPIRH(+) resembles the charge and geometry of the ribooxocarbenium ion transition state and can be looked upon as a transition-state analogue inhibitor; however, it is known that the pAPIR species is initially bound to the enzyme. We have characterized the pAPIRH(+) species as resident of the active site using ab initio calculations and molecular dynamics simulations. This is a novel use of molecular dynamics to investigate the protonation state of the bound ligand to the active site. Nanosecond molecular dynamics simulations reveal a short hydrogen-bonding network between pAPIRH(+)-O2'-Asp14-His241 triad, which is not seen in the crystal structure. Other features discussed are: hydrogen bonding between pAPIRH(+) and Asn168, unusual geometry of the iminoribitol ring, and hydrophobic interactions.     

The glass transition and the distribution of voids in room-temperature ionic liquids: a molecular dynamics study

The glass transition in prototypical room temperature ionic liquids has been investigated by molecular dynamics simulations based on an Amber-like empirical force field. Samples of [C(4)mim][PF(6)], [C(4)mim][Tf(2)N], and [C(3)mim][Tf(2)N] have been quenched from the liquid phase at T = 500 to a glassy state at T ～ 0 K in discontinuous steps of 20 K every 1.2 ns. The glass temperature estimated by simulation (T(g) = 209 K for [C(4)mim][PF(6)], T(g) = 204 K for [C(4)mim][Tf(2)N], and T(g) = 196 K for [C(3)mim][Tf(2)N]) agrees semi-quantitatively with the experimental values (T(g) = 193÷196 K for [C(4)mim][PF(6)], T(g) = 186÷189 K for [C(4)mim][Tf(2)N], and T(g) = 183 K for [C(3)mim][Tf(2)N]). A model electron density is introduced to identify voids in the system. The temperature dependence of the size distribution of voids provided by simulation reproduce well the experimental results of positron annihilation lifetime spectroscopy reported in G. Dlubek, Y. Yu, R. Krause-Rehberg, W. Beichel, S. Bulut, N. Pogodina, I. Krossing, and Ch. Friedrich, J. Chem. Phys. 133, 124502 (2010), with only one free parameter needed to fit the experimental data.     

Structure conversions of cellulose IIII crystal models in solution state: a molecular dynamics study

This paper re-examines our previous molecular dynamics (MD) study on cellulose III_I crystal models with finite dimensions solvated in explicit water molecules. Eight crystal models, differing in a constituent lattice plane and dimensions, were studied. One calculation allowed for O–H and C–H bond stretching, and had a small time step of 0.5 fs. The other calculation adopted non-scaling factors of the 1–4 non-bonded interactions. As in our previous study, in the former MD calculations, six of the eight crystal models exhibited structure conversion with cooperative chain slippages generated by a progressive fiber bend. This converted the initial non-staggered chain packing of cellulose III_I into a near one-quarter staggering and gave the crystal model a triclinic-like configuration. In contrast, in the non-1–4 scaling MD calculations, all of the eight crystal models retained the initial cellulose III_I crystal structure. Another series of non-1–4 scaling MD calculations were performed for the four crystal models containing chains with a degree of polymerization (DP) of 40 at 370 K, which simulated hot water treatment to convert cellulose III_I to Iβ. Some of the hydroxymethyl groups irreversibly rotated from gt into tg conformation. This accompanied exchange of the intrasheet hydrogen bonding scheme along the (1 ?1 0) lattice plane from O2–O6 to O3–O6. The original corrugated (1 ?1 0) chain sheet was partly converted into a cellulose I-like flat chain sheet.     

Exploring and characterizing the folding processes of Lys- and Arg-containing Ala-based peptides: a molecular dynamics study

In this study, molecular dynamics simulations were carried out on Lys- and Arg-containing Ala-based peptides (i.e. Ace-(AAAAK)(n)A-NH(2) and Ace-(AAAAR)(n)A-NH(2), where n=1-4), in order to explore and characterize their folding processes. For the oligopeptides, the evolution of α-helical structure with regard to the whole conformation, as well as to each residue was investigated, and the helix-forming propensities were characterized. On the basis of the helicity curves, representing the alteration of average helicity as a function of time, the typical time values describing the folding processes and subprocesses were identified. In the case of each peptide, the evolution and role of helix-stabilizing, non-local and side-chain-to-backbone H-bonds were examined. The appearing i←i+4 H-bonds pointed out the role of these interactions in the stabilization of α-helical conformations, while the occurring i←i+3 H-bonds indicated the presence of β-turn or 3(10)-helical structures. Studying the formation and role of non-local and side-chain-to-backbone H-bonds led to the observation that these types of interactions produced an effect on the evolution of helical conformations, as well as on the folding processes.     

Structure of transition state in catalysis of chain transfer to a monomer

Translated from Teoreticheskaya Éksperimental'naya Khimiya, Vol. 26, No. 1, p. 128, January–February, 1990.     

Computational study of ketosteroid isomerase: insights from molecular dynamics simulation of enzyme bound substrate and intermediate

Delta(5)-3-Ketosteroid Isomerase (KSI) catalyzes the isomerization of 5,6-unsaturated ketosteroids to their 4,5-unsaturated isomers at a rate approaching the diffusion limit. The isomerization reaction follows a two-step general acid-base mechanism starting with Asp38-CO(2)(-) mediated proton abstraction from a sp(3)-hybridized carbon atom, alpha to carbonyl group, providing a dienolate intermediate. In the second step, Asp38-CO(2)H protonates the C6 of the intermediate providing a 4,5-unsaturated ketosteroid. The details of the mechanism have been highly controversial despite several experimental and computational studies of this enzyme. The general acid-base catalysis has been proposed to involve either a catalytic diad or a cooperative hydrogen bond mechanism. In this paper, we report our results from the 1.5 nanosecond molecular dynamics (MD) simulation of enzyme bound natural substrate (E.S) and enzyme bound intermediate (E.In) solvated in a TIP3P water box. The final coordinates from our MD simulation strongly support the cooperative hydrogen bond mechanism. The MD simulation of E.S and E.In shows that both Tyr14 and Asp99 are hydrogen bonded to the O3 of the substrate or intermediate. The average hydrogen bonding distance between Tyr14-OH and O3 becomes shorter and exhibits less fluctuation on E.S --> E.In. We also observe dynamic motions of water moving in and out of the active site in the E.S structures. This free movement of water disappears in the E.In structures. The active site is shielded by hydrophobic residues, which come together and squeeze out the waters from the active site in the E.In complex.     

Molecular dynamics study of taxadiene synthase catalysis

Molecular dynamics (MD) simulations have been performed to study the dynamic behavior of noncovalent enzyme carbocation complexes involved in the cyclization of geranylgeranyl diphosphate to taxadiene catalyzed by taxadiene synthase (TXS). Taxadiene and the observed four side products originate from the deprotonation of carbocation intermediates. The MD simulations of the TXS carbocation complexes provide insights into potential deprotonation mechanisms of such carbocations. The MD results do not support a previous hypothesis that carbocation tumbling is a key factor in the deprotonation of the carbocations by pyrophosphate. Instead water bridges are identified which may allow the formation of side products via multiple proton transfer reactions. A novel reaction path for taxadiene formation is proposed on the basis of the simulations. © 2018 Wiley Periodicals, Inc.     

On the nature of the transition state in catechol O-methyltransferase. A complementary study based on molecular dynamics and potential energy surface explorations

The way in which enzymes influence the rate of chemical processes is still a question of debate. The protein promotes the catalysis of biochemical processes by lowering the free energy barrier in comparison with the reference uncatalyzed reaction in solution. In this article we are reporting static and dynamic aspects of the enzyme catalysis in a bimolecular reaction, namely a methyl transfer from S-adenosylmethionine to the hydroxylate oxygen of a substituted catechol catalyzed by catechol O-methyltransferase. From QM/MM optimizations, we will first analyze the participation of the environment on the transition vector. The study of molecular dynamics trajectories will allow us to estimate the transmission coefficient from a previously localized transition state as the maximum in the potential of mean force profile. The analysis of the reactive and nonreactive trajectories in the enzyme environment and in solution will also allow studying the geometrical and electronic changes, with special attention to the chemical system movements and the coupling with the environment. The main result, coming from both analyses, is the approximation of the magnesium cation to the nucleophilic and the hydroxyl group of the catecholate as a result of a general movement of the protein, stabilizing in this way the transition state. Consequently, the free energy barrier of the enzyme reaction is dramatically decreased with respect to the reaction in solution.     

Glass transition in single poly(ethylene oxide) chain: A molecular dynamics simulation study

Local dynamics of single poly(ethylene oxide) chain in various environments (bulk, film, and isolated systems) has been characterized by the reorientation functions of various backbone bond vectors. Within any observation time, the variations of these reorientation functions with the temperature can be well described by the Kohlrausch?Williams?Watts (KWW) like equation, in which the fitted temperature parameter is identified as the glass transition temperature (T _g). The so‐obtained T _g for that polymer faithfully reveals the effects of the observation time, chain flexibility and vector range on the local dynamics. Furthermore, it is found that the KWW like relation is also applicable to the temperature‐dependence of the fraction of frozen atoms or torsions defined by the trajectory radii of gyration or the conformational transitions. Consequently, different motions lead to different values of T _g for the same system. Despite all, the consistent trend can be yielded, namely, T _g (bulk) > T _g (film) > T _g (isolated), which captures the effects of free surfaces on enhanced dynamics. In addition, dynamics heterogeneity in the systems can be quantitatively revealed. The newly proposed method holds a bright promise to predict the T _g values of complex polymers especially for comparisons. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55 , 178–188     

Global dynamics and transition state theories: Comparative study of reaction rate constants for gas-phase chemical reactions

In this review article, we present a systematic comparison of the theoretical rate constants for a range of bimolecular reactions that are calculated by using three different classes of theoretical methods: quantum dynamics (QD), quasi-classical trajectory (QCT), and transition state theory (TST) approaches. The study shows that the difference of rate constants between TST results and those of the global dynamics methods (QD and QCT) are seen to be related to a number of factors including the number of degrees-of-freedom (DOF), the density of states at transition state (TS), etc. For reactions with more DOF and higher density of states at the TS, it is found that the rate constants from TST calculations are systematically higher than those obtained from global dynamics calculations, indicating large recrossing effect for these systems. The physical insight of this phenomenon is elucidated in the present review.