Proton transfer function of carbonic anhydrase: Insights from QM/MM simulations

Recent QM/MM analyses of proton transfer function of human carbonic anhydrase II (CAII) are briefly reviewed. The topics include a preliminary analysis of nuclear quadrupole coupling constant calculations for the zinc ion and more detailed analyses of microscopic pK_a of the zinc-bound water and free energy profile for the proton transfer. From a methodological perspective, our results emphasize that performing sufficient sampling is essential to the calculation of all these quantities, which reflects the well solvated nature of CAII active site. From a mechanistic perspective, our analyses highlight the importance of electrostatics in shaping the energetics and kinetics of proton transfer in CAII for its function. We argue that once the pK_a for the zinc-bound water is modulated to be in the proper range (~ 7.0), proton transfer through a relatively well solvated cavity towards/from the protein surface (His64) does not require any major acceleration. Therefore, although structural details like the length of the water wire between the donor and acceptor groups still may make a non-negligible contribution, our computational results and the framework of analysis suggest that the significance of such “fine-tuning” is likely secondary to the modulation of pK_a of the zinc-bound water. We encourage further experimental analysis with mutation of (charged) residues not in the immediate neighborhood of the zinc ion to quantitatively test this electrostatics based framework; in particular, Φ analysis based on these mutations may shed further light into the relative importance of the classical Grotthus mechanism and the “proton hole” pathway that we have proposed recently for CAII.     

Carbinolamine formation and dehydration in a DNA repair enzyme active site

In order to suggest detailed mechanistic hypotheses for the formation and dehydration of a key carbinolamine intermediate in the T4 pyrimidine dimer glycosylase (T4PDG) reaction, we have investigated these reactions using steered molecular dynamics with a coupled quantum mechanics-molecular mechanics potential (QM/MM). We carried out simulations of DNA abasic site carbinolamine formation with and without a water molecule restrained to remain within the active site quantum region. We recovered potentials of mean force (PMF) from thirty replicate reaction trajectories using Jarzynski averaging. We demonstrated feasible pathways involving water, as well as those independent of water participation. The water-independent enzyme-catalyzed reaction had a bias-corrected Jarzynski-average barrier height of approximately (6.5 kcal mol(-1) (27.2 kJ mol(-1)) for the carbinolamine formation reaction and 44.5 kcal mol(-1) (186 kJ mol(-1)) for the reverse reaction at this level of representation. When the proton transfer was facilitated with an intrinsic quantum water, the barrier height was approximately 15 kcal mol(-1) (62.8 kJ mol(-1)) in the forward (formation) reaction and 19 kcal mol(-1) (79.5 kJ mol(-1)) for the reverse. In addition, two modes of unsteered (free dynamics) carbinolamine dehydration were observed: in one, the quantum water participated as an intermediate proton transfer species, and in the other, the active site protonated glutamate hydrogen was directly transferred to the carbinolamine oxygen. Water-independent unforced proton transfer from the protonated active site glutamate carboxyl to the unprotonated N-terminal amine was also observed. In summary, complex proton transfer events, some involving water intermediates, were studied in QM/MM simulations of T4PDG bound to a DNA abasic site. Imine carbinolamine formation was characterized using steered QM/MM molecular dynamics. Dehydration of the carbinolamine intermediate to form the final imine product was observed in free, unsteered, QM/MM dynamics simulations, as was unforced acid-base transfer between the active site carboxylate and the N-terminal amine.     

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.     

Role of active site histidines in the two half-reactions of the aryl-alcohol oxidase catalytic cycle

The crystal structure of aryl-alcohol oxidase (AAO), a flavoenzyme involved in lignin degradation, reveals two active-site histidines, whose role in the two enzyme half-reactions was investigated. The redox state of flavin during turnover of the variants obtained show a stronger histidine involvement in the reductive than in the oxidative half-reaction. This was confirmed by the k(cat)/K(m(Al)) and reduction constants that are 2-3 orders of magnitude decreased for the His546 variants and up to 5 orders for the His502 variants, while the corresponding O(2) constants only decreased up to 1 order of magnitude. These results confirm His502 as the catalytic base in the AAO reductive half-reaction. The solvent kinetic isotope effect (KIE) revealed that hydroxyl proton abstraction is partially limiting the reaction, while the α-deuterated alcohol KIE showed a stereoselective hydride transfer. Concerning the oxidative half-reaction, directed mutagenesis and computational simulations indicate that only His502 is involved. Quantum mechanical/molecular mechanical (QM/MM) reveals an initial partial electron transfer from the reduced FADH(-) to O(2), without formation of a flavin-hydroperoxide intermediate. Reaction follows with a nearly barrierless His502H(+) proton transfer that decreases the triplet/singlet gap. Spin inversion and second electron transfer, concomitant with a slower proton transfer from flavin N5, yields H(2)O(2). No solvent KIE was found for O(2) reduction confirming that the His502 proton transfer does not limit the oxidative half-reaction. However, the small KIE on k(cat)/K(m(Ox)), during steady-state oxidation of α-deuterated alcohol, suggests that the second proton transfer from N5H is partially limiting, as predicted by the QM/MM simulations.     

Hydrogen donor-acceptor fluctuations from kinetic isotope effects: a phenomenological model

Kinetic isotope effects (KIEs) and their temperature dependence can probe the structural and dynamic nature of enzyme-catalyzed proton or hydride transfers. The molecular interpretation of their temperature dependence requires expensive and specialized quantum mechanics/molecular mechanics (QM/MM) calculations to provide a quantitative molecular understanding. Currently available phenomenological models use a nonadiabatic assumption that is not appropriate for most hydride and proton-transfer reactions, while others require more parameters than the experimental data justify. Here we propose a phenomenological interpretation of KIEs based on a simple method to quantitatively link the size and temperature dependence of KIEs to a conformational distribution of the catalyzed reaction. This model assumes adiabatic hydrogen tunneling, and by fitting experimental KIE data, the model yields a population distribution for fluctuations of the distance between donor and acceptor atoms. Fits to data from a variety of proton and hydride transfers catalyzed by enzymes and their mutants, as well as nonenzymatic reactions, reveal that steeply temperature-dependent KIEs indicate the presence of at least two distinct conformational populations, each with different kinetic behaviors. We present the results of these calculations for several published cases and discuss how the predictions of the calculations might be experimentally tested. This analysis does not replace molecular QM/MM investigations, but it provides a fast and accessible way to quantitatively interpret KIEs in the context of a Marcus-like model.     

Catalytic mechanism and product specificity of rubisco large subunit methyltransferase: QM/MM and MD investigations

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of Rubisco large subunit methyltransferase (LSMT). It was found that the appearance of a water channel is required for the stepwise methylation by S-adenosylmethionine (AdoMet). The water channel appears in the presence of AdoMet (LSMT.Lys-NH3+.AdoMet), but is not present immediately after methyl transfer (LSMT.Lys-N(Me)H2+.AdoHcy). The water channel allows proton dissociation from both LSMT.AdoMet.Lys-NH3+ and LSMT.AdoMet.Lys-N(Me)H2+. The water channel does not appear for proton dissociation from LSMT.AdoMet.Lys-N(Me)2H+, and a third methyl transfer does not occur. By QM/MM, the calculated free energy barrier of the first methyl transfer reaction catalyzed by LSMT (Lys-NH2 + AdoMet --> Lys-N(Me)H2+ + AdoHcy) is DeltaG++ = 22.8 +/- 3.3 kcal/mol. This DeltaG++ is in remarkable agreement with the value 23.0 kcal/mol calculated from the experimental rate constant (6.2 x 10-5 s-1). The calculated DeltaG++ of the second methyl transfer reaction (AdoMet + Lys-N(Me)H --> AdoHcy + Lys-N(Me)2H+) at the QM/MM level is 20.5 +/- 3.6 kcal/mol, which is in agreement with the value 22.0 kcal/mol calculated from the experimental rate constant (2.5 x 10-4 s-1). The third methyl transfer (Lys-N(Me)2 + AdoMet --> Lys-N(Me)3+ + AdoHcy) is associated with an allowed DeltaG++ of 25.9 +/- 3.2 kcal/mol. However, this reaction does not occur because a water channel does not form to allow the proton dissociation of Lys-N(Me)2H+. Future studies will determine whether the product specificity of lysine (mono, di, and tri) methyltransferases is determined by the formation of water channels.     

Protein effects in non-heme iron enzyme catalysis: insights from multiscale models

Many non-heme iron enzymes have similar sets of ligands but still catalyze widely different reactions. A key question is, therefore, the role of the protein in controlling reactivity and selectivity. Examples from multiscale simulations, primarily QM/MM, of both mono- and binuclear non-heme iron enzymes are used to analyze the stability of these models and what they reveal about the protein effects. Consistent results from QM/MM modeling are the importance of the hydrogen bond network to control reactivity and electrostatic stabilization of electron transfer from second-sphere residues. The long-range electrostatic effects on reaction barriers are small for many systems. In the systems where large electrostatic effects have been reported, these lead to higher barriers. There is thus no evidence of any significant long-range electrostatic effects contributing to the catalytic efficiency of non-heme iron enzymes. However, the correct evaluation of electrostatic contributions is challenging, and the correlation between calculated residue contributions and the effects of mutation experiments is not very strong. The largest benefits of QM/MM models are thus the improved active-site geometries, rather than the calculation of accurate energies. Reported differences in mechanistic predictions between QM and QM/MM models can be explained by differences in hydrogen bonding patterns in and around the active site. Correctly constructed cluster models can give results with similar accuracy as those from multiscale models, but the latter reduces the risk of drawing the wrong mechanistic conclusions based on incorrect geometries and are preferable for all types of modeling, even when using very large QM parts.     

An Alternative Mechanism for the Methylation of Phosphoethanolamine Catalyzed by Plasmodium falciparum Phosphoethanolamine Methyltransferase

The phosphobase methylation pathway catalyzed by the phosphoethanolamine methyltransferase in Plasmodium falciparum (PfPMT), the malaria parasite, offers an attractive target for anti-parasitic drug development. PfPMT methylates phosphoethanolamine (pEA) to phosphocholine for use in membrane biogenesis. Quantum mechanics and molecular mechanics (QM/MM) calculations tested the proposed reaction mechanism for methylation of pEA involving the previously identified Tyr-19–His-132 dyad, which indicated an energetically unfavorable mechanism. Instead, the QM/MM calculations suggested an alternative mechanism involving Asp-128. The reaction coordinate involves the stepwise transfer of a proton to Asp-128 via a bridging water molecule followed by a typical S_n2-type methyl transfer from S-adenosylmethionine to pEA. Functional analysis of the D128A, D128E, D128Q, and D128N PfPMT mutants shows a loss of activity with pEA but not with the final substrate of the methylation pathway. X-ray crystal structures of the PfPMT-D128A mutant in complex with S-adenosylhomocysteine and either pEA or phosphocholine reveal how mutation of Asp-128 disrupts a hydrogen bond network in the active site. The combined QM/MM, biochemical, and structural studies identify a key role for Asp-128 in the initial step of the phosphobase methylation pathway in Plasmodium and provide molecular insight on the evolution of multiple activities in the active site of the PMT.     

Enzyme structure captures four cysteines aligned for disulfide relay

Thioredoxin superfamily proteins introduce disulfide bonds into substrates, catalyze the removal of disulfides, and operate in electron relays. These functions rely on one or more dithiol/disulfide exchange reactions. The flavoenzyme quiescin sulfhydryl oxidase (QSOX), a catalyst of disulfide bond formation with an interdomain electron transfer step in its catalytic cycle, provides a unique opportunity for exploring the structural environment of enzymatic dithiol/disulfide exchange. Wild‐type Rattus norvegicus QSOX1 (RnQSOX1) was crystallized in a conformation that juxtaposes the two redox‐active di‐cysteine motifs in the enzyme, presenting the entire electron‐transfer pathway and proton‐transfer participants in their native configurations. As such a state cannot generally be enriched and stabilized for analysis, RnQSOX1 gives unprecedented insight into the functional group environments of the four cysteines involved in dithiol/disulfide exchange and provides the framework for analysis of the energetics of electron transfer in the presence of the bound flavin adenine dinucleotide cofactor. Hybrid quantum mechanics/molecular mechanics (QM/MM) free energy simulations based on the X‐ray crystal structure suggest that formation of the interdomain disulfide intermediate is highly favorable and secures the flexible enzyme in a state from which further electron transfer via the flavin can occur.     

Optimization of cutting schemes for the evaluation of molecular electrostatic potentials in proteins via Moving-Domain QM/MM

This work presents new developments of the moving-domain QM/MM (MoD-QM/MM) method for modeling protein electrostatic potentials. The underlying goal of the method is to map the electronic density of a specific protein configuration into a point-charge distribution. Important modifications of the general strategy of the MoD-QM/MM method involve new partitioning and fitting schemes and the incorporation of dynamic effects via a single-step free energy perturbation approach (FEP). Selection of moderately sized QM domains partitioned between and C (from C=O), with incorporation of delocalization of electrons over neighboring domains, results in a marked improvement of the calculated molecular electrostatic potential (MEP). More importantly, we show that the evaluation of the electrostatic potential can be carried out on a dynamic framework by evaluating the free energy difference between a non-polarized MEP and a polarized MEP. A simplified form of the potassium ion channel protein Gramicidin-A from Bacillus brevis is used as the model system for the calculation of MEP. Figure Schematic representation of the Moving Domain QM/MM method     

Substrate diffusion and oxidation in GMC oxidoreductases: an experimental and computational study on fungal aryl-alcohol oxidase

AAO (aryl-alcohol oxidase) provides H?O? in fungal degradation of lignin, a process of high biotechnological interest. The crystal structure of AAO does not show open access to the active site, where different aromatic alcohols are oxidized. In the present study we investigated substrate diffusion and oxidation in AAO compared with the structurally related CHO (choline oxidase). Cavity finder and ligand diffusion simulations indicate the substrate-entrance channel, requiring side-chain displacements and involving a stacking interaction with Tyr?2. Mixed QM (quantum mechanics)/MM (molecular mechanics) studies combined with site-directed mutagenesis showed two active-site catalytic histidine residues, whose substitution strongly decreased both catalytic and transient-state reduction constants for p-anisyl alcohol in the H502A (over 1800-fold) and H546A (over 35-fold) variants. Combination of QM/MM energy profiles, protonation predictors, molecular dynamics, mutagenesis and pH profiles provide a robust answer regarding the nature of the catalytic base. The histidine residue in front of the FAD ring, AAO His??2 (and CHO His???), acts as a base. For the two substrates assayed, it was shown that proton transfer preceded hydride transfer, although both processes are highly coupled. No stable intermediate was observed in the energy profiles, in contrast with that observed for CHO. QM/MM, together with solvent KIE (kinetic isotope effect) results, suggest a non-synchronous concerted mechanism for alcohol oxidation by AAO.     

Study of angiotensin-(1-7) vasoactive peptide and its beta-cyclodextrin inclusion complexes: complete sequence-specific NMR assignments and structural studies

We report the complete sequence-specific hydrogen NMR assignments of vasoactive peptide angiotensin-(1-7) (Ang-(1-7)). Assignments of the majority of the resonances were accomplished by COSY, TOCSY, and ROESY peak coordinates at 400MHz and 600MHz. Long-side-chain amino acid spin system identification was facilitated by long-range coherence transfer experiments (TOCSY). Problems with overlapped resonance signals were solved by analysis of heteronuclear 2D experiments (HSQC and HMBC). Nuclear Overhauser effects (NOE) results were used to probe peptide conformation. We show that the inclusion of the angiotensin-(1-7) tyrosine residue is favored in inclusion complexes with beta-cyclodextrin. QM/MM simulations at the DFTB/UFF level confirm the experimental NMR findings and provide detailed structural information on these compounds in aqueous solution.     

Copper binding sites in the C-terminal domain of mouse prion protein: A hybrid (QM/MM) molecular dynamics study

We present a hybrid QM/MM Car-Parrinello molecular dynamics study of the copper-loaded C-terminal domain of the mouse prion protein. By means of a statistical analysis of copper coordination in known protein structures, we localized the protein regions with the highest propensity for copper ion binding. The identified candidate structures were subsequently refined via QM/MM simulations. Their EPR characteristics were computed to make contact with the experimental data and to probe the sensitivity to structural and chemical changes. Overall best agreement with the experimental EPR data (Van Doorslaer et al., J Phys Chem B 2001; 105: 1631-1639) and the information currently available in the literature is observed for a binding site involving H187. Moreover, a reinterpretation of the experimental proton hyperfine couplings was possible in the light of the present computational findings.     

DRF90: a polarizable force field

The direct reaction field (DRF) approach has proven to be a useful tool to investigate the influence of solvents on the quantum/classical behaviour of solute molecules. In this paper, we report the latest extension of this DRF approach, which consists of the gradient of the completely classical energy expressions of this otherwise QM/MM method. They can be used in (completely classical) molecular dynamics (MD) simulations and geometry optimizations, that can be followed by a number of single point QM/MM calculations on configurations obtained in these simulations/optimizations. We report all energy and gradient expressions, and results for a number of interesting (model) systems. They include geometry optimization of the benzene dimer as well as MD simulations of some solvents. The most stable configuration for the benzene dimer is shown to be the parallel-displaced form, which is slightly more stable (0.3 kcal/mol) than the T-shaped dimer.     

QM/MM free energy simulations: recent progress and challenges

Due to the higher computational cost relative to pure molecular mechanical (MM) simulations, hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations particularly require a careful consideration of balancing computational cost and accuracy. Here, we review several recent developments in free energy methods most relevant to QM/MM simulations and discuss several topics motivated by these developments using simple but informative examples that involve processes in water. For chemical reactions, we highlight the value of invoking enhanced sampling technique (e.g. replica-exchange) in umbrella sampling calculations and the value of including collective environmental variables (e.g. hydration level) in metadynamics simulations; we also illustrate the sensitivity of string calculations, especially free energy along the path, to various parameters in the computation. Alchemical free energy simulations with a specific thermodynamic cycle are used to probe the effect of including the first solvation shell into the QM region when computing solvation free energies. For cases where high-level QM/MM potential functions are needed, we analyse two different approaches: the QM/MM–MFEP method of Yang and co-workers and perturbative correction to low-level QM/MM free energy results. For the examples analysed here, both approaches seem productive although care needs to be exercised when analysing the perturbative corrections.     

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations

Recent computer simulations of the cysteine nucleophilic attack on propanal in human mitochondrial aldehyde dehydrogenase (ALDH2) yielded an unexpected result: the chemically reasonable formation of a dead-end cysteine-cofactor adduct when NAD+ was in the "hydride transfer" position. More recently, this adduct found independent crystallographic support in work on formyltetrahydrofolate dehydrogenase, work which further found evidence of the same adduct on re-examination of deposited electron densities of ALDH2. Although the experimental data showed that this adduct was reversible, several mechanistic questions arise from the fact that it forms at all. Here, we present results from further quantum mechanical/molecular mechanical (QM/MM) simulations toward understanding the mechanistic implications of adduct formation. These simulations revealed formation of the oxyanion thiohemiacetal intermediate only when the nicotinamide ring of NAD+ is oriented away from the active site, contrary to prior arguments. In contrast, and in seeming paradox, when NAD is oriented to receive the hydride, disassociation of the oxyanion intermediate to form the dead-end adduct is more thermodynamically favored than maintaining the oxyanion intermediate necessary for catalysis to proceed. However, this disassociation to the adduct could be avoided through proton transfer from the enzyme to the intermediate. Our results continue to indicate that the unlikely source of this proton is the Cys302 main chain amide.     

The Role of Gln61 in HRas GTP hydrolysis: a quantum mechanics/molecular mechanics study

Activation of the water molecule involved in GTP hydrolysis within the HRas·RasGAP system is analyzed using a tailored approach based on hybrid quantum mechanics/molecular mechanics (QM/MM) simulation. A new path emerges: transfer of a proton from the attacking water molecule to a second water molecule, then a different proton is transferred from this second water molecule to the GTP. Gln(61) will stabilize the transient OH(-) and H(3)O(+) molecules thus generated. This newly proposed mechanism was generated by using, for the first time to our knowledge, the entire HRas-RasGAP protein complex in a QM/MM simulation context. It also offers a rational explanation for previous experimental results regarding the decrease of GTPase rate found in the HRas Q61A mutant and the increase exhibited by the HRas Q61E mutant.     

Recent developments in QM/MM methods towards open-boundary multi-scale simulations

Combined quantum mechanics/molecular mechanics (QM/MM) methods have been widely used in multi-scale modelling and simulations of physical, chemical and biological processes in complex environments. In this review, we provide an overview of the recently developed QM/MM algorithms, with emphasis on our works, towards the ultimate goal of establishing an open boundary between the QM and MM subsystems. The open boundary is characterised by on-the-fly exchanges of partial charges and atoms between the QM and MM subsystems, allowing us to focus on the small QM subsystem of primary interest in dynamics simulations. An open-boundary scheme has the promise to the utilisations of small QM subsystems, high-levels of QM theory and long simulation times, which can potentially lead to new insights.