Molecular dynamics simulations of the first steps of the reaction catalyzed by HIV-1 protease

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.     

Intramembrane proton binding site linked to activation of bacterial pentameric ion channel

Prokaryotic orthologs of eukaryotic Cys-loop receptor channels recently emerged as structural and mechanistic surrogates to investigate this superfamily of intercellular signaling proteins. Here, we examine proton activation of the prokaryotic ortholog GLIC using patch clamp electrophysiology, mutagenesis, and molecular dynamics (MD) simulations. Whole-cell current recordings from human embryonic kidney (HEK) 293 cells expressing GLIC show half-maximal activation at pH 6, close to the pK(a) of histidine, implicating the three native His residues in proton sensing linked to activation. The mutation H235F abolishes proton activation, H277Y is without effect, and all nine mutations of His-127 prevent expression on the cell surface. In the GLIC crystal structure, His-235 on transmembrane (TM) α-helix 2, hydrogen bonds to the main chain carbonyl oxygen of Ile-259 on TM α-helix 3. MD simulations show that when His-235 is protonated, the hydrogen bond persists, and the channel remains in the open conformation, whereas when His-235 is deprotonated, the hydrogen bond dissociates, and the channel closes. Mutations of the proximal Tyr-263, which also links TM α-helices 2 and 3 via a hydrogen bond, alter proton sensitivity over a 1.5 pH unit range. MD simulations show that mutations of Tyr-263 alter the hydrogen bonding capacity of His-235. The overall findings show that His-235 in the TM region of GLIC is a novel proton binding site linked to channel activation.     

Cofactor activation and substrate binding in pyruvate decarboxylase. Insights into the reaction mechanism from molecular dynamics simulations

We have performed long-term molecular dynamics simulations of pyruvate decarboxylase from Zymomonas mobilis. Nine structures were modeled to investigate mechanistic questions related to binding of the cofactor, thiamin diphosphate (ThDP), and the substrate in the active site. The simulations reveal that the proposed three ThDP-tautomers all can bind in the active site and indicate that the equilibrium is shifted toward 4'-aminopyrimidine ThDP in the absence of substrate. 4'-Aminopyrimidinium ThDP is found to be a likely intermediate in the equilibrium. Mutations of important active site residues, Glu473Ala and Glu50Ala, were modeled to further elucidate their catalytic role. Formation of the catalytic important ylide by deprotonation of ThDP(C2) is investigated. Only the less favored tautomer, 1',4'-iminopyrimidine ThDP (imino-ThDP), could be deprotonated. The two other tautomers of ThDP could not be activated at the C2-position, thus, explaining the mechanistic importance of the less stable imino-ThDP. Finally, binding of pyruvate in the active site with the cofactor modeled as the nucleophilic ylide (ylide-ThDP) is studied. The carbonyl group of the substrate forms a hydrogen bond to Tyr290(OH). No hydrogen bond could be identified between ThDP(N4') and the substrate. The geometry of the substrate binding is well-suited for a nucleophilic attack by ylide-ThDP(C2). We propose that a proton relay from His113 via Asp27 and Tyr290 to the carbonyl oxygen atom of the substrate may be involved in the mechanism.     

Hypothesis of a proton switch in QM/MM modelling of interaction of dUMP analogues with thymidylate synthase

Quantum mechanical, molecular mechanics and molecular dynamics (MD) methods were used to investigate initial steps of 2′-deoxyuridine-5′-monophosphate (dUMP) methylation catalysed by the thymidylate synthase (TS) enzyme. The amino acid residues surrounding the active site within a 10 Å radius sphere were modelled with the combined quantum mechanical (B3LYP/LANL2DZ) and molecular mechanics ONIOM double-layer method. The results indicated the initial nucleophilic attack of Cys146 on dUMP to be concerted with formation of a hydrogen bond to the oxygen O4 of dUMP. Moreover, the proton in the vicinity of the O4 atom appears to act as a ‘proton switch’: if a proton is present near O4, it stabilises the S(Cys146)–C6(dUMP) sulphur–carbon bond, but if it is absent, the sulphur–carbon bond does not form. If the O4 oxygen is replaced by sulphur atom, the ‘switch effect’ does not occur. The suggested correlation between the strength of hydrogen bond involving O4 oxygen and the ability of dUMP to form bonds at C6 corresponds well to the crystal structures of TS complexes available in the Protein Data Bank. In the vast majority of crystal structures, the presence of the S(Cys146)–C6(dUMP) bond was coupled with the presence of hydrogen bond between the dUMP O4 atom and the conserved Asn177. The ‘proton switch’ hypothesis is supported also by the results of MD studies of TS binary complexes, suggesting that average distance separating S(Cys146) and C6(dUMP) becomes distinctly shorter in the presence of hydrogen bonding between Asn177 and O4.     

Theoretical Insights into Catalytic Mechanism of Protein Arginine Methyltransferase 1

Protein arginine methyltransferase 1 (PRMT1), the major arginine asymmetric dimethylation enzyme in mammals, is emerging as a potential drug target for cancer and cardiovascular disease. Understanding the catalytic mechanism of PRMT1 will facilitate inhibitor design. However, detailed mechanisms of the methyl transfer process and substrate deprotonation of PRMT1 remain unclear. In this study, we present a theoretical study on PRMT1 catalyzed arginine dimethylation by employing molecular dynamics (MD) simulation and quantum mechanics/molecular mechanics (QM/MM) calculation. Ternary complex models, composed of PRMT1, peptide substrate, and S-adenosyl-methionine (AdoMet) as cofactor, were constructed and verified by 30-ns MD simulation. The snapshots selected from the MD trajectory were applied for the QM/MM calculation. The typical S_N2-favored transition states of the first and second methyl transfers were identified from the potential energy profile. Deprotonation of substrate arginine occurs immediately after methyl transfer, and the carboxylate group of E144 acts as proton acceptor. Furthermore, natural bond orbital analysis and electrostatic potential calculation showed that E144 facilitates the charge redistribution during the reaction and reduces the energy barrier. In this study, we propose the detailed mechanism of PRMT1-catalyzed asymmetric dimethylation, which increases insight on the small-molecule effectors design, and enables further investigations into the physiological function of this family.     

Effects of the position and manner of hydration on the stability of solvated N-methylacetamides and the strength of binding between N-methylacetamide and water clusters: a computational study

This work mainly studies the effects of the position (there are two possible hydrated sites) and the manner (i.e., whether water acts as a proton donor or acceptor) of hydration by various numbers of water molecules on the stability of 14 solvated N-methylacetamide structures, NMA-(H₂O) _n (n = 1–3), as well as the binding strength between the NMA and the water cluster, using molecular dynamics (MD) and B3LYP methods. Natural bond orbital (NBO) analysis is used to explore the origin of these effects. Some novel observations are obtained from the work. Our results show that monohydration at the carbonyl site favors stability and binding strength compared to monohydration at the amino site. Similarly, the preferred hydration at the carbonyl site is observed for dihydrated NMAs when the second water is added as a proton donor to the C=O group or the first water is H-bonded to the C=O group. However, unfavorable hydration at the C=O site occurs if the second water acts as a proton acceptor. Trihydration by a ring cluster of three water molecules at either the carbonyl site or the amino one yields relatively stable complexes, but significantly disfavors binding strength. The other trihydrated NMAs show similar behavior to dihydrated NMAs. In addition, our results show that the C=O and N–H frequencies can still be utilized to examine the H-bond effects of the water cluster.     

Reactivity of an Excess Electron with Monovalent Cations in Bulk Water by Mixed Quantum Classical Molecular Dynamics Simulations

A mixed quantum classical molecular dynamics (QCMD) simulation of the silver and sodium cations in presence of an excess electron is reported. The silver cation is shown to be reduced by the hydrated electron and to form a stable, highly polarized, neutral atom. On the contrary, the sodium cation is not reduced and a metastable contact ion pair is observed. The resulting absorption spectra of both species are compared with experiments and shown to be in good agreement. Furthermore, the free energy curve for the charge separation was calculated and rationalized in terms of a thermodynamic cycle. Finally, a direct, reactive, molecular dynamics trajectory provides some useful informations on the reduction mechanism.     

Structure and dynamics of a proton wire: a theoretical study of H+ translocation along the single-file water chain in the gramicidin A channel.

The rapid translocation of H+ along a chain of hydrogen-bonded water molecules, or proton wire, is thought to be an important mechanism for proton permeation through transmembrane channels. Computer simulations are used to study the properties of the proton wire formed by the single-file waters in the gramicidin A channel. The model includes the polypeptidic dimer, with 22 water molecules and one excess proton. The dissociation of the water molecules is taken into account by the "polarization model" of Stillinger and co-workers. The importance of quantum effects due to the light mass of the hydrogen nuclei is examined with the use of discretized Feynman path integral molecular dynamics simulations. Results show that the presence of an excess proton in the pore orients the single-file water molecules and affects the geometry of water-water hydrogen bonding interactions. Rather than a well-defined hydronium ion OH3+ in the single-file region, the protonated species is characterized by a strong hydrogen bond resembling that of O2H5+. The quantum dispersion of protons has a small but significant effect on the equilibrium structure of the hydrogen-bonded water chain. During classical trajectories, proton transfer between consecutive water molecules is a very fast spontaneous process that takes place in the subpicosecond time scale. The translocation along extended regions of the chain takes place neither via a totally concerted mechanism in which the donor-acceptor pattern would flip over the entire chain in a single step, nor via a succession of incoherent hops between well-defined intermediates. Rather, proton transfer in the wire is a semicollective process that results from the subtle interplay of rapid hydrogen-bond length fluctuations along the water chain. These rapid structural fluctuations of the protonated single file of waters around an average position and the slow movements of the average position of the excess proton along the channel axis occur on two very different time scales. Ultimately, it is the slow reorganization of hydrogen bonds between single-file water molecules and channel backbone carbonyl groups that, by affecting the connectivity and the dynamics of the single-file water chain, also limits the translocation of the proton across the pore.     

Role of methionine 71 in substrate recognition and structural integrity of bacterial peptidyl-tRNA hydrolase

BackgroundBacterial peptidyl-tRNA hydrolase (Pth) is an essential enzyme that alleviates tRNA starvation by recycling prematurely dissociated peptidyl-tRNAs. The specificity of Pth for N-blocked-aminoacyl-tRNA has been proposed to be contingent upon conserved residue N14 forming a hydrogen bond with the carbonyl of the first peptide bond in the substrate. M71 is involved in forming a conserved hydrogen bond with N14. Other interactions facilitating this recognition are not known.MethodsThe structure, dynamics, and stability of the M71A mutant of Pth from Vibrio cholerae (VcPth) were characterized by X-ray crystallography, NMR spectroscopy, MD simulations and DSC.ResultsCrystal structure of M71A mutant was determined. In the structure, the dimer interface is formed by the insertion of six C-terminal residues of one molecule into the active site of another molecule. The side-chain amide of N14 was hydrogen bonded to the carbonyl of the last peptide bond formed between residues A196 and E197, and also to A71. The CSP profile of mutation was similar to that observed for the N14D mutant. M71A mutation lowered the thermal stability of the protein.ConclusionOur results indicate that the interactions of M71 with N14 and H24 play an important role in optimal positioning of their side-chains relative to the peptidyl-tRNA substrate. Overall, these interactions of M71 are important for the activity, stability, and compactness of the protein.SignificanceThe work presented provides original and new structural and dynamics information that significantly enhances our understanding of the network of interactions that govern this enzyme's activity and selectivity.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Simulation of the substrate cavity dynamics of quercetinase

Molecular dynamics (MD) simulations have been performed on quercetin 2,3 dioxygenase (2,3QD) to study the mobility and flexibility of the substrate cavity. 2,3QD is the only firmly established Cu-containing dioxygenase known so far. It catalyses the breakage of the O-heterocycle of flavonols. The substrates occupy a shallow and overall hydrophobic cavity proximal to the metal centre of the homo-dimeric enzyme. The linker connecting the C-terminal and N-terminal domains in the monomer is partly disordered in the crystal structure and part of it forms a flexible lid at the entrance of the substrate cavity. This loop has been tentatively assigned a role in the enzyme mechanism: it helps lock the substrate into place. The dynamics of this loop has been investigated by MD simulation. The initial coordinates were taken from the crystal structure of 2,3QD in the presence of the substrate kaempferol (KMP). After equilibration and simulation over 7.2ns the substrate was removed and another equilibration and simulation of 7.2ns was performed. The results show that the structures of the free enzyme as well as of the enzyme-substrate complex are stable in MD simulation. The linker shows strongly enhanced mobility in the loop region that is close to the entrance to the substrate cavity (residues 154-169). Movement of the loop takes place on a timescale of 5-10ns. To confirm the conclusions about the loop dynamics drawn from the 7.2ns simulation, the simulation was extended with another 8ns. When substrate binds into the cavity the loop orders remarkably, although mobility is retained by residues 155-158. Some regions of the loop (residues 154-160 and 164-176) move over a considerable distance and approach the substrate closely, reinforcing the idea that they lock the substrate in the substrate cavity. The enthalpic component of the interaction of the loop with the protein and the KMP appears to favour the locking of the substrate. Two water molecules were found immobilised in the cavity, one of which exhibited rotation on the picosecond timescale. When the substrate is removed, the empty cavity fills up with water within 200ps.     

Structural insights by molecular dynamics simulations into specificity of the major human AP endonuclease toward the benzene-derived DNA adduct, pBQ-C

The benzetheno exocyclic adduct of the cytosine (C) base (pBQ-C) is a product of reaction between DNA and a stable metabolite of the human carcinogen benzene, p-benzoquinone (pBQ). We reported previously that the pBQ-C-containing duplex is a substrate for the human AP endonuclease (APE1), an enzyme that cleaves an apurinic/apyrimidinic (AP) site from double stranded DNA. In this work, using molecular dynamics simulation (MD), we provided a structural explanation for the recognition of the pBQ-C adduct by APE1. Molecular modeling of the DNA duplex containing pBQ-C revealed significant displacement of this adduct toward the major groove with pronounced kinking of the DNA at the lesion site, which could serve as a structural element recognized by the APE1 enzyme. Using 3 ns MD it was shown that the position of the pBQ-C adduct is stabilized by two hydrogen bonds formed between the adduct and the active site amino acids Asp 189 and Ala 175. The pBQ-C/APE1 complex, generated by MD, has a similar hydrogen bond network between target phosphodiester bond at the pBQ-C site and key amino acids at the active site, as in the crystallographically determined APE1 complexed with an AP site-containing DNA duplex. The position of the adduct at the enzyme active site, together with the hydrogen bond network, suggests a similar reaction mechanism for phosphodiester bond cleavage of oligonucleotide containing pBQ-C as reported for the AP site.     

Ab initio models for receptor-ligand interactions in proteins. 4. Model assembly study of the catalytic mechanism of triosephosphate isomerase

The catalytic mechanism of triosephosphate isomerase (TIM) was investigated with ab initio quantum mechanical calculations. Electrostatic interactions between the quantum mechanical active site and the protein and solvent environment were modeled using the finite difference Poission-Boltzman method. The complexes of TIM with the substrate dihydroxyacetone phosphate (DHAP), five possible intermediates and the product glyceraldehyde-3-phosphate (GAP) were optimized in the active-site model at the 3-21G^(*) level and energy profile for the proton abstraction from DHAP by the active-site Glu167 was calculated at the MP2/3-21G^(*)//3-21G^(*) level. Calculated energetics of the enzyme reaction were found to be in reasonable agreement with the experimental findings. Calculations revealed that an enediol of the substrate is a probable intermediate in the enzyme reaction. It was suggested that the proton abstracted from the substrate by the active-site glutamate goes to the carbonyl oxygen of the substrate producing enediol intermediate either directly or after it is exchanged with solvent. © 1996 Wiley-Liss, Inc.     

Solution conformation of a cyclic neurokinin antagonist: A nmr and molecular dynamics study

The solution structure of a hexapeptide, cyclo(Gln-Trp-Phe-Gly-Leu-Met), which is a selective NK-2 antagonist, has been studied by a combination of two-dimensional nmr and molecular dynamics (MD) techniques. The simulation based on nmr and MD data resulted in the convergence to a family of structures. Free molecular dynamics for 50 ps in the presence of DMSO solvent molecules shows that the structure is energetically stable. One intramolecular hydrogen bond between the amide proton of Gin and the carbonyl oxygen of Gly was revealed. This result is consistent with the results from the measurement of the temperature coefficient of the amide protons. The extent of intermolecular hydrogen bonding between the amide protons of the peptide and DMSO was also revealed by the free MD simulation. The resulting structure of the cyclic peptide contains a variation type I′ β-turn in the Gly-Leu-Met-Gln segment. Comparison of the structure of this peptide with that of other NK-2 antagonist cyclic hexapeptides was made, and the activity of cyclic antagonists appears to be inversely related to the conformational rigidity of the cyclic peptides. © 1994 John Wiley & Sons, Inc.     

Computer simulation of primary kinetic isotope effects in the proposed rate-limiting step of the glyoxalase I catalyzed reaction

The proposed rate-limiting step of the glyoxalase I catalyzed reaction is the proton abstraction from the C1 carbon of the substrate by Glu(172). Here we examine primary kinetic isotope effects and the influence of quantum dynamics on this process by computer simulations. The calculations utilize the empirical valence bond method in combination with the molecular dynamics free energy perturbation technique and path integral simulations. For the enzyme-catalyzed reaction a H/D kinetic isotope effect of 5.0 +/- 1. 3 is predicted in reasonable agreement with the experimental result of about 3. Furthermore, the magnitude of quantum mechanical effects is found to be very similar for the enzyme reaction and the corresponding uncatalyzed process in solution, in agreement with other studies. The problems associated with attaining the required accuracy in order for the present approach to be useful as a diagnostic tool for the study of enzyme reactions are also discussed.     

Approximate quantum trajectory dynamics for reactive processes in condensed phase

A method of molecular dynamics with quantum corrections, practical for studies of large molecular systems, is reviewed. The approach is based on the Bohmian formulation of the time-dependent Schrödinger equation in which a wavefunction is represented by an ensemble of interdependent trajectories. The quantum effects come from the quantum potential acting on trajectories on par with the usual classical potential. The quantum potential is determined from the evolving nuclear wavefunction, i.e. from the quantum trajectory (QT) ensemble itself. For practical and conceptual reasons the quantum potential and corresponding quantum nuclear effect are computed only for the selected light nuclei. For studies of reactive chemical processes, the classical potential is computed on-the-fly using the density functional tight binding method of electronic structure. A massively parallel implementation, based on the message passing interface allows for efficient simulations of ensembles of thousands of trajectories describing systems of up to 200 atoms. As a biochemical application, the approximate QT approach is used to model the tunnelling-dominated proton transfer in soybean-lipoxygenase-1. A materials science application is represented by a study of the nuclear quantum effect on adsorption of hydrogen and deuterium on a C₃₇H₁₅ molecule, which is a model ‘flake’ of graphene.     

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.     

Insights into Cu(I) exchange in HAH1 using quantum mechanical and molecular simulations

The human antioxidant protein, HAH1, is an important participant in a Cu(I) transport chain, delivering one Cu(I) ion to the Wilson's (WND) or Menkes disease protein (MNK). Full geometry optimizations and second-derivative calculations were performed on several binding site models using the B3LYP functional to derive parameters for the construction of a novel molecular mechanical (MM) force field for Cu(I) and its ligating residues in HAH1. MM minimization and molecular dynamics (MD) calculations were then performed using the AMBER suite to validate the newly generated force field. The X-ray crystal structure of the protein and the geometry of the Cu(I) binding site within the protein were reproduced by the MD simulations on the protein based on rmsd and visual inspection, validating the new force field parameters. The results from the quantum mechanical (QM) and MD simulations suggest that either a two- or three-coordinate exchange reaction is preferred and that it is unlikely that a four-coordinate Cu(I) species plays a role in copper exchange.