A hybrid elastic band string algorithm for studies of enzymatic reactions

Design of prodrug for stereoisomers of omapatrilat to cross the blood-brain barrier using docking, homology modeling, MD, and QM/MM methods

In this study, we designed a suitable ester prodrug for omapatrilat to penetrate the blood-brain barrier and treat CNS diseases. Based on the ADMET properties, the methyl carboxylate ester of omapatrilat was chosen from among several prodrug structures. Sixteen methyl carboxylate esters were constructed for omapatrilat. The structure of brain carboxylesterase was derived via homology modeling, and molecular docking was used to determine the most potent stereoisomers against brain carboxylesterase. The top three stereoisomer complexes, and the apo form of the protein, were then considered using molecular dynamics simulation and MM/GBSA analysis. Following the simulation, structural analysis was performed using RMSD, RMSF, Rg, and hydrogen bond analysis tools. Our data demonstrated that the prodrug of RSSR is a suitable structure for crossing the blood-brain barrier and binding to brain carboxylesterase. In addition, we found via QM/MM calculation that the catalytic reaction of the prodrug of RSSR against brain carboxylesterase occurs via two steps, including acylation and diacylation steps. Based on our findings, we propose a clinical trial of a methyl carboxylate ester prodrug of omapatrilat's RSSR for the treatment of brain diseases.Communicated by Ramaswamy H. Sarma.

Theoretical characterization of the shikimate 5-dehydrogenase reaction from Mycobacterium tuberculosis by hybrid QC/MM simulations and quantum chemical descriptors

In this study, we have investigated the enzyme shikimate 5-dehydrogenase from the causative agent of tuberculosis, Mycobacterium tuberculosis. We have employed a mixture of computational techniques, including molecular dynamics, hybrid quantum chemical/molecular mechanical potentials, relaxed surface scans, quantum chemical descriptors and free-energy simulations, to elucidate the enzyme's reaction pathway. Overall, we find a two-step mechanism, with a single transition state, that proceeds by an energetically uphill hydride transfer, followed by an energetically downhill proton transfer. Our mechanism and calculated free energy barrier for the reaction, 64.9 kJ mol^- 1, are in good agreement with those predicted from experiment. An analysis of quantum chemical descriptors along the reaction pathway indicated a possibly important, yet currently unreported, role of the active site threonine residue, Thr65.

Hierarchical Markov State Model Building to Describe Molecular Processes

Markov state models can describe ensembles of pathways via kinetic networks but are difficult to create when large free-energy barriers limit unbiased sampling. Chain-of-states simulations allow sampling over large free-energy barriers but are often constructed using a single pathway that is unlikely to thermodynamically average over orthogonal degrees of freedom in complex systems. Here, we combine the advantages of these two approaches in the form of a Markov state model of Markov state models, which we call a Hierarchical Markov state model. In this approach, independent Markov models are constructed in regions of configuration space that are locally well sampled but are separated by large free-energy barriers from other regions. A string method is used to construct an ensemble of pathways connecting the states of these different local Markov models, and the rate through each pathway is then estimated. These rates are then combined with the rate information from the local Markov models in a master equation to predict global rates, fluxes, and populations. By applying this hierarchical approach to tractable systems, a toy potential and dipeptides, we demonstrate that it is more accurate than the conventional single-pathway description. The advantages of this approach are that it (i) is more realistic than the conventional chain-of-states approach, as an ensemble of pathways rather than a single pathway is used to describe processes in high-dimensional systems, and (ii) it resolves the issue of poor sampling in Markov State model building when large free-energy barriers are present. The divide-and-conquer strategy inherent to this approach should make this procedure straightforward to apply to more complex systems.

Negative Impact of Carbapenem Methylation on the Reactivity of β-Lactams for Cysteine Acylation as Revealed by Quantum Calculations and Kinetic Analyses

The Enterococcus faecium l,d-transpeptidase (Ldt_fm) mediates resistance to most β-lactam antibiotics in this bacterium by replacing classical peptidoglycan polymerases. The catalytic Cys of Ldt_fm is rapidly acylated by β-lactams belonging to the carbapenem class but not by penams or cephems. We previously reported quantum calculations and kinetic analyses for Ldt_fm and showed that the inactivation profile is not determined by differences in drug binding (K_D [equilibrium dissociation constant] values in the 50 to 80 mM range). In this study, we analyzed the reaction of a Cys sulfhydryl with various β-lactams in the absence of the enzyme environment in order to compare the intrinsic reactivity of drugs belonging to the penam, cephem, and carbapenem classes. For this purpose, we synthesized cyclic Cys-Asn (cCys-Asn) to generate a soluble molecule with a sulfhydryl closely mimicking a cysteine in a polypeptide chain, thereby avoiding free reactive amino and carboxyl groups. Computational studies identified a thermodynamically favored pathway involving a concerted rupture of the β-lactam amide bond and formation of an amine anion. Energy barriers indicated that the drug reactivity was the highest for nonmethylated carbapenems, intermediate for methylated carbapenems and cephems, and the lowest for penams. Electron-withdrawing groups were key reactivity determinants by enabling delocalization of the negative charge of the amine anion. Acylation rates of cCys-Asn determined by spectrophotometry revealed the same order in the reactivity of β-lactams. We concluded that the rate of Ldt_fm acylation is largely determined by the β-lactam reactivity with one exception, as the enzyme catalytic pocket fully compensated for the detrimental effect of carbapenem methylation.

Cyclization Reaction Catalyzed by Cyclodipeptide Synthases Relies on a Conserved Tyrosine Residue

Cyclodipeptide synthases (CDPSs) form various cyclodipeptides from two aminoacyl tRNAs via a stepwise mechanism with the formation of a dipeptidyl enzyme intermediate. As a final step of the catalytic reaction, the dipeptidyl group undergoes intramolecular cyclization to generate the target cyclodipeptide product. In this work, we investigated the cyclization reaction in the cyclodipeptide synthase AlbC using QM/MM methods and free energy simulations. The results indicate that the catalytic Y202 residue is in its neutral protonated form, and thus, is not likely to serve as a general base during the reaction. We further demonstrate that the reaction relies on the conserved residue Y202 serving as a proton relay, and the direct proton transfer from the amino group to S37 of AlbC is unlikely. Calculations reveal that the hydroxyl group of tyrosine is more suitable for the proton transfer than hydroxyl groups of other amino acids, such as serine and threonine. Results also show that the residues E182, N40, Y178 and H203 maintain the correct conformation of the dipeptide needed for the cyclization reaction. The mechanism discovered in this work relies on the amino groups conserved among the entire CDPS family and, thus is expected to be universal among CDPSs.

Leaching from Palladium Nanoparticles in an Ionic Liquid Leads to the Formation of Ionic Monometallic Species

Molecular dynamics simulations and DFT calculations suggest that leaching of palladium species from Pd nanoparticles in ionic liquids does not involve "naked" Pd(0) atoms or neutral ArPdX species formed by oxidative addition of arylhalides. Instead, the ionic liquid contributes largely to leaching of ionic PdX^- or PdAr⁺ species.

Mechanisms and Specificity of Phenazine Biosynthesis Protein PhzF

Phenazines are bacterial virulence and survival factors with important roles in infectious disease. PhzF catalyzes a key reaction in their biosynthesis by isomerizing (2 S,3 S)-2,3-dihydro-3-hydroxy anthranilate (DHHA) in two steps, a [1,5]-hydrogen shift followed by tautomerization to an aminoketone. While the [1,5]-hydrogen shift requires the conserved glutamate E45, suggesting acid/base catalysis, it also shows hallmarks of a sigmatropic rearrangement, namely the suprafacial migration of a non-acidic proton. To discriminate these mechanistic alternatives, we employed enzyme kinetic measurements and computational methods. Quantum mechanics/molecular mechanics (QM/MM) calculations revealed that the activation barrier of a proton shuttle mechanism involving E45 is significantly lower than that of a sigmatropic [1,5]-hydrogen shift. QM/MM also predicted a large kinetic isotope effect, which was indeed observed with deuterated substrate. For the tautomerization, QM/MM calculations suggested involvement of E45 and an active site water molecule, explaining the observed stereochemistry. Because these findings imply that PhzF can act only on a limited substrate spectrum, we also investigated the turnover of DHHA derivatives, of which only O-methyl and O-ethyl DHHA were converted. Together, these data reveal how PhzF orchestrates a water-free with a water-dependent step. Its unique mechanism, specificity and essential role in phenazine biosynthesis may offer opportunities for inhibitor development.

Aminoacetylation Reaction Catalyzed by Leucyl-tRNA Synthetase Operates via a Self-Assisted Mechanism Using a Conserved Residue and the Aminoacyl Substrate

Leucyl-tRNA synthetase catalyzes attachment of leucine amino acid to its cognate tRNA. During the second, aminoacetylation, step of the reaction, the leucyl moiety is transferred from leucyl-adenylate to the terminal A76 adenosine of tRNA. In this work, we have investigated the aminoacetylation step catalyzed by leucyl-tRNA synthase, using ab initio quantum chemical/molecular mechanical hybrid potentials in conjunction with reaction-path-location algorithms and molecular dynamics free energy simulations. We have modeled reaction mechanisms arising from both crystallographic studies and computational work. We invoke various groups as potential proton acceptors-namely, the phosphate and leucyl amino groups of leucyl-adenylate, the A76 base of tRNA, and the Asp80 and Glu532 residues of the protein-and consider both metal-assisted and metal-free reactions. Free energy calculations indicate that both the phosphate group of leucyl adenylate and Glu532 are not strong bases. This agrees with the results of the quantum chemical/molecular mechanical reaction path calculations which give high free energy barriers for the studied pathways involving these groups. A self-assisted mechanism with the leucyl amino group and Asp80 as proton acceptors is the most likely. Furthermore, in this mechanism the presence of a metal ion coordinated by the phosphate group and Glu532 strongly activates the reaction.

QM/MM investigation of ATP hydrolysis in aqueous solution

Adenosine-5'-triphosphate (ATP) hydrolysis represents a most important reaction in biology. Despite extensive research efforts, the mechanism for ATP hydrolysis in aqueous solution still remains under debate. Previous theoretical studies often predefined reaction coordinates to characterize the mechanism for ATP hydrolysis in water with Mg(2+) by evaluating free energy profiles through these preassumed reaction paths. In the present work, a nudged elastic band method is applied to identify the minimum energy path calculated with a hybrid quantum mechanics and molecular mechanics approach. Along the reaction path, the free energy profile was obtained to have a single transition state and the activation energy of 32.5 kcal/mol. This transition state bears a four-centered structure that describes a concerted nature of the reaction. In the More-O'Ferrall-Jencks diagram, the results show that the reaction proceeds through a concerted path before the system reaches the transition state and along an associative path after the transition state. In addition, the calculated reaction free energy is -7.0 kcal/mol, in good agreement with experiment, capturing the exothermic feature of MgATP(2-) hydrolysis in aqueous solution, whereas the reaction was often shown to be endothermic in the previous theoretical studies. As Mg(2+) is required for ATP hydrolysis in cells, its role in the reaction is also elucidated.

Identification of a second GTP-bound magnesium ion in archaeal initiation factor 2

Eukaryotic and archaeal translation initiation processes involve a heterotrimeric GTPase e/aIF2 crucial for accuracy of start codon selection. In eukaryotes, the GTPase activity of eIF2 is assisted by a GTPase-activating protein (GAP), eIF5. In archaea, orthologs of eIF5 are not found and aIF2 GTPase activity is thought to be non-assisted. However, no in vitro GTPase activity of the archaeal factor has been reported to date. Here, we show that aIF2 significantly hydrolyses GTP in vitro. Within aIF2γ, H97, corresponding to the catalytic histidine found in other translational GTPases, and D19, from the GKT loop, both participate in this activity. Several high-resolution crystal structures were determined to get insight into GTP hydrolysis by aIF2γ. In particular, a crystal structure of the H97A mutant was obtained in the presence of non-hydrolyzed GTP. This structure reveals the presence of a second magnesium ion bound to GTP and D19. Quantum chemical/molecular mechanical simulations support the idea that the second magnesium ion may assist GTP hydrolysis by helping to neutralize the developing negative charge in the transition state. These results are discussed in light of the absence of an identified GAP in archaea to assist GTP hydrolysis on aIF2.

Enzymatic minimum free energy path calculations using swarms of trajectories

The development of approaches for simulating rare events in complex molecular systems is a central concern in chemical physics. In recent work, Roux and co-workers proposed a novel, swarms of trajectories (SoT) method for determining the transition paths of such events. It consists of the dynamical refinement on the system's free energy surface of a putative transition path that is parametrized in terms of a set of collective variables (CVs) that are identified as being important for the transition. In this work, we have implemented the SoT method and used it to investigate the catalytic mechanisms of two enzymatic reactions using hybrid QM/MM potentials. Our aim has been to test the performance of SoT for enzyme systems and to devise robust simulation protocols that can be employed in future studies of this type. We identify the conditions under which converged results can be obtained using inertial and Brownian dynamical evolutions of the CVs, show that the inclusion of several CVs can give significant additional insight into the mechanisms of the reactions, and show that the use of minimum energy paths as starting guesses can greatly accelerate path refinement.

Challenges in computational studies of enzyme structure, function and dynamics

In this review we give an overview of the field of Computational enzymology. We start by describing the birth of the field, with emphasis on the work of the 2013 chemistry Nobel Laureates. We then present key features of the state-of-the-art in the field, showing what theory, accompanied by experiments, has taught us so far about enzymes. We also briefly describe computational methods, such as quantum mechanics-molecular mechanics approaches, reaction coordinate treatment, and free energy simulation approaches. We finalize by discussing open questions and challenges.

The mechanism of citryl-coenzyme A formation catalyzed by citrate synthase

The enzyme citrate synthase is used by all living cells to catalyze the first step of the citric acid cycle. In this work, we have investigated the enolization and condensation steps catalyzed by citrate synthase, using ab initio (B3LYP/def2-TZVP and MP2/aug-cc-pVDZ) quantum chemical/molecular mechanical hybrid potentials in conjunction with reaction-path-location algorithms and molecular dynamics free energy simulations. The results of the latter indicate that the catalytic His238 residue is in its neutral form, and also argue strongly for the presence of a water molecule in the enzyme's catalytic center. Such a water is observed in some, but not all, of the experimentally resolved structures of the protein. The mechanism itself starts with an enolization that proceeds via an enolate intermediate rather than the enol form, which is much more unstable. This is in agreement with the results of other workers. For the condensation step, we investigated two mechanisms in which there is a direct nucleophilic attack of the enolate intermediate on the oxaloacetate carbonyl carbon, and found the one in which there is no proton transfer from the neighboring arginine to be preferred. Although this residue, Arg329, is not implicated directly in the reaction, it helps to stabilize the negative citryl-CoA formed during the condensation step.

Mechanism of activation of elongation factor Tu by ribosome: catalytic histidine activates GTP by protonation

Elongation factor Tu (EF-Tu) is central to prokaryotic protein synthesis as it has the role of delivering amino-acylated tRNAs to the ribosome. Release of EF-Tu, after correct binding of the EF-Tu:aa-tRNA complex to the ribosome, is initiated by GTP hydrolysis. This reaction, whose mechanism is uncertain, is catalyzed by EF-Tu, but requires activation by the ribosome. There have been a number of mechanistic proposals, including those spurred by a recent X-ray crystallographic analysis of a ribosome:EF-Tu:aa-tRNA:GTP-analog complex. In this work, we have investigated these and alternative hypotheses, using high-level quantum chemical/molecular mechanical simulations for the wild-type protein and its His85Gln mutant. For both proteins, we find previously unsuggested mechanisms as being preferred, in which residue 85, either His or Gln, directly assists in the reaction. Analysis shows that the RNA has a minor catalytic effect in the wild-type reaction, but plays a significant role in the mutant by greatly stabilizing the reaction's transition state. Given the similarity between EF-Tu and other members of the translational G-protein family, it is likely that these mechanisms of ribosome-activated GTP hydrolysis are pertinent to all of these proteins.