Exploring nucleoside hydrolase catalysis in silico: molecular dynamics study of enzyme-bound substrate and transition state

The mechanism of action of inosine-uridine nucleoside hydrolase has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Five MD studies have been performed with enzyme substrate complex (E.S), enzyme substrate complex with protonated His241 (EH.S), enzyme transition state complex (E.TS), enzyme transition state complex with protonated His241 (EH.TS), and His241Ala transition state complex E(H241A).TS. Special attention has been given to the role of His241, which has been considered as the general acid catalyst to assist departure of the leaving nucleobase on the basis of its location in the active site in the X-ray crystal structure (). Yet on the basis of the location in the active site, Tyr229 is closer to the aniline ring of pAPIR as compared to His241. On initiation of MD simulations, His241 does not approach the nucleobase in the structures of EH.S, E.S, EH.TS, and E.TS. In the solvated enzyme, Tyr229, which is a member of the hydrogen bonding network inosine O2'.Asp14.His241.Tyr229.inosine N7, serves as a proton source to the leaving nucleobase. The loss of significant activity of His241Ala mutant is shown to be related to the disruption of the above hydrogen bonded network and the distancing of Tyr229 from inosine N7. The structures of the enzyme complexes with substrate or TS are not visibly altered on protonation of His241, a most unusual outcome. The bell-shaped pH dependence upon pK(app)'s of 7.1 and 9.1 may be attributed to the necessity of the dissociation of Asp10 or Asp15 and the acid form of Tyr229, respectively. In TS, the residue Ile81 migrated closer, whereas Arg233 moved away from the nucleobase. The probability of ribooxocarbenium ion stabilization by Asn168 and Asp14 is discussed. The Asp14-CO(2)(-) is hydrogen bonded to the ribose 2'-OH for 96% of the MD simulation time. Nucleophilic addition of water138 to ribooxocarbenium ion is suggested to be assisted by the proton shuttle from water138 --> Asp10 --> Asp15 --> water pool. An anticorrelation motion between Tyr229-OH and Asn168-OD1 in EH.S and E.S is observed. The relationship of this anticorrelated motion to mechanism, if any, deserves further exploration, perhaps the formation of a near attack conformation.     

Effect of monovalent ion binding on molecular dynamics of the S100-family calcium-binding protein calbindin D9k

Calbindin D_9k is a member of the S100 subfamily of EF-hand calcium binding proteins, and has served as an important model system for biophysical studies. The fast timescale dynamics of the calcium-free (apo) state is characterized using molecular dynamics simulations. Order parameters for the backbone NH bond vectors are determined from the simulations and compared with experimentally derived values, with a focus on the dynamics of calcium-binding site I. There is a significant discrepancy between simulated and experimental order parameters for site I residues in the case of no ion bound in site I. However, it was found in the simulations that a Na⁺ ion can bind in site I, and the resulting order parameters determined from the simulations are in excellent agreement with experiment. Comparisons are made to X-ray structures of other S100 family members in which Na⁺ ions were observed or suggested to be bound in site I. © 2019 Wiley Periodicals, Inc.     

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.     

Characterization of the Structure and Dynamics of the HDV Ribozyme at Different Stages Along the Reaction Path

The structure and dynamics of the hepatitis delta virus ribozyme (HDVr) are studies using molecular dynamics simulations at several stages along its catalytic reaction path, including reactant, activated precursor, transition state mimic and product states, departing from an initial structure based on the C75U mutant crystal structure (PDB: 1VC7). Results of five 350 ns molecular dynamics simulations reveal a spontaneous rotation of U-1 that leads to an in-line conformation and support the role of protonated C75 as the general acid in the transition state. Our results provide rationale for the interpretation of several important experimental results, and make experimentally testable predictions regarding the roles of key active site residues that are not obvious from any available crystal structures.     

Calculations on the structure and spectral properties of cytochrome c551 using DFT and ONIOM methods

The active site geometry of cytochrome (Cyt) c(551) and its mutated form containing Fe(II) and Fe(III) ions have been calculated using density functional theory (DFT)-based Becke's three-parameter hybrid exchange and Lee-Yang-Parr correlation (B3LYP) method. In addition, calculations have also been carried out using hybrid meta DFT-based M06 functional. The effect of the protein milieu on the active site geometry has also been probed using two-layer via our own N-layered integrated molecular orbital + molecular mechanics (ONIOM) method. Evidence from the calculations reveal that the active site geometry is not significantly affected by the oxidation state of metal ion. The difference in the geometry of the active site and that of the same with the entire protein environment is only minimal, which shows that the protein milieu does not influence the structure of the active site. The calculated electronic transition energies from the time-dependent DFT (TDDFT) calculations are in close agreement with the experimental values. Although there are no significant variations in the active site geometry upon oxidation, the changes in the electronic transition energies have been attributed to the reduction in the overlap of metal ion with the ligand orbitals. In addition, it is found that mutation does not influence the active site geometry and the electronic transition energies. Nevertheless, mutation leads to the formation of more compact structure than the native Cyt c(551).     

Coordination number of zinc ions in the phosphotriesterase active site by molecular dynamics and quantum mechanics

We have run several molecular dynamics (MD) simulations on zinc-containing phosphotriesterase (PTE) with two bound substrates, sarin and paraoxon, and with the substrate analog diethyl 4-methylbenzylphosphonate. A standard nonbonded model was employed to treat the zinc ions with the commonly used charge of +2. In all the trajectories, we observed a tightly bound water (TBW) molecule in the active site that was coordinated to the less buried zinc ion. The phosphoryl oxygen of the substrate/inhibitor was found to be coordinated to the same zinc ion so that, considering all ligands, the less buried zinc was hexa-coordinated. The hexa-coordination of this zinc ion was not seen in the deposited X-ray pdb files for PTE. Several additional MD simulations were then performed using different charges (+1, +1.5) on the zinc ions, along with ab initio and density functional theory (DFT) calculations, to evaluate the following possibilities: the crystal diffraction data were not correctly interpreted; the hexa-coordinated zinc ion in PTE is only present in solution and not in the crystal; and the hexa-coordinated zinc ion in PTE is an artifact of the force field used. A charge of +1.5 leads to a coordination number (CN) of 5 on both zinc ions, which is consistent with the results from ab initio and DFT calculations and with the latest high resolution X-ray crystal structure. The commonly used charge of +2 produces a CN of 6 on the less buried zinc. The CN on the more buried zinc ion is 5 when the substrate/inhibitor is present in the simulation, and increases to 6 when the substrate/inhibitor is removed prior to the simulation. The results of both of the MD and quantum mechanical calculations lead to the conclusion that the zinc ions in the PTE active site are both penta-coordinated, and that the MD simulations performed with the charge of +2 overestimate the CN of the zinc ions in the PTE active site. The overall protein structures in the simulations remain unaffected by the change in zinc charge from +2 to +1.5. The results also suggest that the charge +1.5 is the most appropriate for the molecular dynamics simulations on zinc-containing PTE when a nonbonded model is used and no global thermodynamic conclusion is sought. We also show that the standard nonbonded model is not able to properly treat the CN and energy at the same time. A preliminary, promising charge-transfer model is discussed with the use of the zinc charge of +1.5.     

MM-PBSA free energy analysis of endo-1,4-xylanase II (XynII)-substrate complexes: binding of the reactive sugar in a skew boat and chair conformation

The binding of xylotetraose in different conformations to the active site of endo-1,4-beta-xylanase II (XynII) from Trichoderma reesei was studied using molecular dynamics (MD) simulations and free energy analyses employing the MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Surface Area) method. MD simulations of 1 ns were done for the substrate xylotetraose having the reactive sugar, which is bound in the -1 subsite of XynII in the 4C1 (chair) and 2So (skew boat) ground state conformations, and for the transition state of the XynII catalysed hydrolysis of the beta-glycosidic linkage. According to the simulations and free energy analysis, XynII binds the substrate with the -1 sugar in the 2So conformation 59.8 kJ mol(-1) tighter than the substrate with the sugar in the 4C1 conformation. The reactive 2So conformation resembles closely the reaction transition state and has the breaking glycosidic bond in a pseudo-axial orientation ready for facile bond cleavage. The transition state was calculated to be bound 77.1 kJ mol(-1) tighter than the 4C1 ground state conformation. The molecular mechanical interaction energy between the enzyme and the reactive pyranoside unit at the -1 subsite was 75.7 kJ mol(-1) more favorable for the binding of the 2So conformation than the 2C1 conformation, explaining the clearly tighter binding of the reactive structure The results of this study indicate that in the Michaelis complex XynII, a member of the family 11 xylanase, the substrate is bound in a skew boat conformation and in the catalytic reaction, the -1 sugar proceeds from the 4C1 conformation through 2So to the transition state with the -1 sugar in the 2,5B conformation.     

Microsecond timescale MD simulations at the transition state of PmHMGR predict remote allosteric residues

Understanding the mechanisms of enzymatic catalysis requires a detailed understanding of the complex interplay of structure and dynamics of large systems that is a challenge for both experimental and computational approaches. More importantly, the computational demands of QM/MM simulations mean that the dynamics of the reaction can only be considered on a timescale of nanoseconds even though the conformational changes needed to reach the catalytically active state happen on a much slower timescale. Here we demonstrate an alternative approach that uses transition state force fields (TSFFs) derived by the quantum-guided molecular mechanics (Q2MM) method that provides a consistent treatment of the entire system at the classical molecular mechanics level and allows simulations at the microsecond timescale. Application of this approach to the second hydride transfer transition state of HMG-CoA reductase from Pseudomonas mevalonii (PmHMGR) identified three remote residues, R396, E399 and L407, (15–27 Å away from the active site) that have a remote dynamic effect on enzyme activity. The predictions were subsequently validated experimentally via site-directed mutagenesis. These results show that microsecond timescale MD simulations of transition states are possible and can predict rather than just rationalize remote allosteric residues.

Transition state force fields enable MD simulations at the transition state of HMGCoA reductase that sample the transition state ensemble on the μs timescale to identify remote residues that affect the reaction rate.     

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.     

Ion size influence on the Ar solvation shells of M(+)-C6F6 clusters (M = Na, K, Rb, Cs)

The size-specific influence of alkali metal ions in the gradual transition from cluster rearrangement to solvation dynamics is investigated by means of molecular dynamics simulations for alkali metal cation-hexafluorobenzene systems, M(+)-C(6)F(6) (M = Na, K, Rb and Cs), surrounded by Ar atoms. To analyze such transition, different small aggregates of the M(+)-C(6)F(6)-Ar(n) (n = 1, ..., 30) type and M(+)-C(6)F(6) clusters solvated by about 500 Ar atoms are considered. The Ar-C(6)F(6) interaction contribution has been described using two different formalisms, based on the interaction decomposition in atom-bond and in atom-effective atom terms, which have been applied to study the small aggregates and to investigate the Ar solvated M(+)-C(6)F(6) clusters, respectively. The selectivity of the promoted phenomena from the M(+) ion size and their dependence from the number of Ar atoms is characterized.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

Reversibly controllable guest binding in precisely defined cavities: selectivity, induced fit, and switching in novel resorcin[4]arene-based container molecules

Two molecular baskets are presented, which were constructed based on a resorcin[4]arene platform. The molecules completely surround suitable guests, such as cyclo- or oxacycloalkanes, and bind them with high strength. The thermodynamic parameters for inclusion complexation were determined as well as the influence of encapsulation on the ring inversion barrier of bound cyclohexane. Two-dimensional NMR spectroscopy clearly shows the existence of a directed attractive interaction between oxacyclohexane and one of the hosts, which constrains the rotation of the bound molecule. Both containers exhibit remarkable binding selectivity as a consequence of their precisely defined structures. They both differentiate between homologous cycloalkanes, and whereas cyclohexane binds best within the larger of the two interior cavities, cyclopentane fits best in the smaller one. The selectivity is governed by ideal filling of space. We have conducted molecular dynamics experiments to understand the thermal fluctuations in the cavity sizes when a guest is bound. The simulations show that within a very narrow range the hosts adapt their binding site to different guests in order to optimize the fraction of occupied space. Once a binding geometry is established, it is characterized by a very low degree of flexibility. The guest-hosting properties of both molecules can be suspended by an external stimulus: addition of acid induces an opening of portals in the structures and immediately releases all bound cargo. Neutralization of the solution completely restores the initial state.     

Link between allosteric signal transduction and functional dynamics in a multisubunit enzyme: S-adenosylhomocysteine hydrolase

S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key role in methylation reactions including those required for maturation of viral mRNA, is an important drug target in the discovery of antiviral agents. While targeting the active site is a straightforward strategy of enzyme inhibition, evidence of allosteric modulation of active site in many enzymes underscores the molecular origin of signal transduction. Information of co-evolving sequences in SAHH family and the key residues for functional dynamics that can be identified using native topology of the enzyme provide glimpses into how the allosteric signaling network, dispersed over the molecular structure, coordinates intra- and intersubunit conformational dynamics. To study the link between the allosteric communication and functional dynamics of SAHHs, we performed Brownian dynamics simulations by building a coarse-grained model based on the holo and ligand-bound structures. The simulations of ligand-induced transition revealed that the signal of intrasubunit closure dynamics is transmitted to form intersubunit contacts, which in turn invoke a precise alignment of active site, followed by the dimer-dimer rotation that compacts the whole tetrameric structure. Further analyses of SAHH dynamics associated with ligand binding provided evidence of both induced fit and population shift mechanisms and also showed that the transition-state ensemble is akin to the ligand-bound state. Besides the formation of enzyme-ligand contacts at the active site, the allosteric couplings from the residues distal to the active site are vital to the enzymatic function.     

Dissection of the difference between the group I metal ions in inhibiting GSK3β: a computational study

Glycogen synthase kinase 3β (GSK3β) is a serine/threonine kinase that requires two cofactor Mg(2+) ions for catalysis in regulating many important cellular signals. Experimentally, Li(+) is a competitive inhibitor of GSK3β relative to Mg(2+), while this mechanism is not experienced with other group I metal ions. Herein, we use native Mg(2)(2+)-Mg(1)(2+) GSK3β and its Mg(2)(2+)-M(1)(+) (M = Li, Na, K, and Rb) derivatives to investigate the effect of metal ion substitution on the mechanism of inhibition through two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations and molecular dynamics (MD) simulations. The results of ONIOM calculations elucidate that the interaction of Na(+), K(+), and Rb(+) with ATP is weaker compared to that of Mg(2+) and Li(+) with ATP, and the critical triphosphate moiety of ATP undergoes a large conformational change in the Na(+), K(+), and Rb(+) substituted systems. As a result, the three metal ions (Na(+), K(+), and Rb(+)) are not stable and depart from the active site, while Mg(2+) and Li(+) can stabilize in the active site, evident in MD simulations. Comparisons of Mg(2)(2+)-Mg(1)(2+) and Mg(2)(2+)-Li(1)(+) systems reveal that the inline phosphor-transfer of ATP and the two conserved hydrogen bonds between Lys85 and ATP, together with the electrostatic potential at the Li(1)(+) site, are disrupted in the Mg(2)(2+)-Li(1)(+) system. These computational results highlight the possible mechanism why Li(+) inhibits GSK3β.     

Mechanism of OMP decarboxylation in orotidine 5'-monophosphate decarboxylase

Despite extensive experimental and theoretical studies, the detailed catalytic mechanism of orotidine 5'-monophosphate decarboxylase (ODCase) remains controversial. In particular simulation studies using high level quantum mechanics have failed to reproduce experimental activation free energy. One common feature of many previous simulations is that there is a water molecule in the vicinity of the leaving CO2 group whose presence was only observed in the inhibitor bound complex of ODCase/BMP. Various roles have even been proposed for this water molecule from the perspective of stabilizing the transition state and/or intermediate state. We hypothesize that this water molecule is not present in the active ODCase/OMP complex. Based on QM/MM minimum free energy path simulations with accurate density functional methods, we show here that in the absence of this water molecule the enzyme functions through a simple direct decarboxylation mechanism. Analysis of the interactions in the active site indicates multiple factors contributing to the catalysis, including the fine-tuned electrostatic environment of the active site and multiple hydrogen-bonding interactions. To understand better the interactions between the enzyme and the inhibitor BMP molecule, simulations were also carried out to determine the binding free energy of this special water molecule in the ODCase/BMP complex. The results indicate that the water molecule in the active site plays a significant role in the binding of BMP by contributing approximately -3 kcal/mol to the binding free energy of the complex. Therefore, the complex of BMP plus a water molecule, instead of the BMP molecule alone, better represents the tight binding transition state analogue of ODCase. Our simulation results support the direct decarboxylation mechanism and highlight the importance of proper recognition of protein bound water molecules in the protein-ligand binding and the enzyme catalysis.     

Hydration structure of water confined between mica surfaces

We report further molecular dynamics simulations on the structure of bound hydration layers under extreme confinement between mica surfaces. We find that the liquid phase of water is maintained down to 2 monolayer (ML) thick, whereas the structure of the K(+) ion hydration shell is close to the bulk structure even under D = 0.92 nm confinement. Unexpectedly, the density of confined water remains approximately the bulk value or less, whereas the diffusion of water molecules decreases dramatically. Further increase in confinement leads to a transition to a bilayer ice, whose density is much less than that of ice Ih due to the formation of a specific hydrogen-bonding network.     

Binding conformation of 2-oxoamide inhibitors to group IVA cytosolic phospholipase A2 determined by molecular docking combined with molecular dynamics

The group IVA cytosolic phospholipase A(2) (GIVA cPLA(2)) plays a central role in inflammation. Long chain 2-oxoamides constitute a class of potent GIVA cPLA(2) inhibitors that exhibit potent in vivo anti-inflammatory and analgesic activity. We have now gained insight into the binding of 2-oxoamide inhibitors in the GIVA cPLA(2) active site through a combination of molecular docking calculations and molecular dynamics simulations. Recently, the location of the 2-oxoamide inhibitor AX007 within the active site of the GIVA cPLA(2) was determined using a combination of deuterium exchange mass spectrometry followed by molecular dynamics simulations. After the optimization of the AX007-GIVA cPLA(2) complex using the docking algorithm Surflex-Dock, a series of additional 2-oxoamide inhibitors have been docked in the enzyme active site. The calculated binding affinity presents a good statistical correlation with the experimental inhibitory activity (r(2) = 0.76, N = 11). A molecular dynamics simulation of the docking complex of the most active compound has revealed persistent interactions of the inhibitor with the enzyme active site and proves the stability of the docking complex and the validity of the binding suggested by the docking calculations. The combination of molecular docking calculations and molecular dynamics simulations is useful in defining the binding of small-molecule inhibitors and provides a valuable tool for the design of new compounds with improved inhibitory activity against GIVA cPLA(2).     

HIV-1 protease flaps spontaneously close to the correct structure in simulations following manual placement of an inhibitor into the open state

We report unrestrained, all-atom molecular dynamics simulations of HIV-1 protease (HIV-PR) with a continuum solvent model that reproducibly sample closing of the active site flaps following manual placement of a cyclic urea inhibitor into the substrate binding site of the open protease. The open form was obtained from the unbound, semi-open HIV-PR crystal structure, which we recently reported (Hornak, V.; et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 915-920.) to have spontaneously opened during unrestrained dynamics. In those simulations, the transiently open flaps always returned to the semi-open form that is observed in all crystal structures of the free protease. Here, we show that manual docking of the inhibitor reproducibly induces spontaneous conversion to the closed form as seen in all inhibitor-bound HIV-PR crystal structures. These simulations reproduced not only the greater degree of flap closure, but also the striking difference in flap "handedness" between bound and free enzyme. In most of the simulations, the final structures were highly accurate. Root-mean-square deviations (RMSD) from the crystal structure of the complex were approximately 1.5 A (averaged over the last 100 ps) for the inhibitor and each flap despite initial RMSD of 2-5 A for the inhibitors and 6-11 A for the flaps. Key hydrogen bonds were formed between the flap tips and between flaps and inhibitor that match those seen in the crystal structure. The results demonstrate that all-atom simulations have the ability to significantly improve poorly docked ligand conformations and reproduce large-scale receptor conformational changes that occur upon binding.     

Quantum Mechanical/Molecular Mechanical Study of the HDV Ribozyme: Impact of the Catalytic Metal Ion on the Mechanism

A recent crystal structure of the precleaved HDV ribozyme along with biochemical data support a mechanism for phosphodiester bond self-cleavage in which C75 acts as a general acid and bound Mg(2+) ion acts as a Lewis acid. Herein this precleaved crystal structure is used as the basis for quantum mechanical/molecular mechanical calculations. These calculations indicate that the self-cleavage reaction is concerted with a phosphorane-like transition state when a divalent ion, Mg(2+) or Ca(2+), is bound at the catalytic site but is sequential with a phosphorane intermediate when a monovalent ion, such as Na(+), is at this site. Electrostatic potential calculations suggest that the divalent metal ion at the catalytic site lowers the pK(a) of C75, leading to the concerted mechanism in which the proton is partially transferred to the leaving group in the phosphorane-like transition state. These observations are consistent with experimental data, including pK(a) measurements, reaction kinetics, and proton inventories with divalent and monovalent ions.     

The catalytic mechanism of an aspartic proteinase explored with neutron and X-ray diffraction

Hydrogen atoms play key roles in enzyme mechanism, but as this study shows, even high-quality X-ray data to a resolution of 1 A cannot directly visualize them. Neutron diffraction, however, can locate deuterium atoms even at resolutions around 2 A. Both neutron and X-ray diffraction data have been used to investigate the transition state of the aspartic proteinase endothiapepsin. The different techniques reveal a different part of the story, revealing the clearest picture yet of the catalytic mechanism by which the enzyme operates. Room temperature neutron and X-ray diffraction data were used in a newly developed joint refinement software package to visualize deuterium atoms within the active site of the enzyme when a gem-diol transition state analogue inhibitor is bound at the active site. These data were also used to estimate their individual occupancy, while analysis of the differences between the bond lengths of the catalytic aspartates was performed using atomic resolution X-ray data. The two methods are in agreement on the protonation state of the active site with a transition state analogue inhibitor bound confirming the catalytic mechanism at which the enzyme operates.