Unraveling the Catalytic Pathway of Metalloenzyme Farnesyltransferase through QM/MM Computation

Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations

The initial phases of drug discovery - in silico drug design - could benefit from first principle Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics (MD) simulations in explicit solvent, yet many applications are currently limited by the short time scales that this approach can cover. Developing scalable first principle QM/MM MD interfaces fully exploiting current exascale machines - so far an unmet and crucial goal - will help overcome this problem, opening the way to the study of the thermodynamics and kinetics of ligand binding to protein with first principle accuracy. Here, taking two relevant case studies involving the interactions of ligands with rather large enzymes, we showcase the use of our recently developed massively scalable Multiscale Modeling in Computational Chemistry (MiMiC) QM/MM framework (currently using DFT to describe the QM region) to investigate reactions and ligand binding in enzymes of pharmacological relevance. We also demonstrate for the first time strong scaling of MiMiC-QM/MM MD simulations with parallel efficiency of ∼70% up to >80,000 cores. Thus, among many others, the MiMiC interface represents a promising candidate toward exascale applications by combining machine learning with statistical mechanics based algorithms tailored for exascale supercomputers.

Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures

Drug action is inherently multiscale: it connects molecular interactions to emergent properties at cellular and larger scales. Simulation techniques at each of these different scales are already central to drug design and development, but methods capable of connecting across these scales will extend understanding of complex mechanisms and the ability to predict biological effects. Improved algorithms, ever-more-powerful computing architectures and the accelerating growth of rich datasets are driving advances in multiscale modeling methods capable of bridging chemical and biological complexity from the atom to the cell. Particularly exciting is the development of highly detailed, structure-based, physical simulations of biochemical systems, which are now able to access experimentally relevant timescales for large systems and, at the same time, achieve unprecedented accuracy. In this Perspective, we discuss how emerging data-rich, physics-based multiscale approaches are of the cusp of realizing long-promised impact in the discovery, design and development of novel therapeutics. We highlight emerging methods and applications in this growing field, and outline how different scales can be combined in practical modelling and simulation strategies.

Role of substrate dynamics in protein prenylation reactions

CONSPECTUS: The role dynamics plays in proteins is of intense contemporary interest. Fundamental insights into how dynamics affects reactivity and product distributions will facilitate the design of novel catalysts that can produce high quality compounds that can be employed, for example, as fuels and life saving drugs. We have used molecular dynamics (MD) methods and combined quantum mechanical/molecular mechanical (QM/MM) methods to study a series of proteins either whose substrates are too far away from the catalytic center or whose experimentally resolved substrate binding modes cannot explain the observed product distribution. In particular, we describe studies of farnesyl transferase (FTase) where the farnesyl pyrophosphate (FPP) substrate is ∼8 Å from the zinc-bound peptide in the active site of FTase. Using MD and QM/MM studies, we explain how the FPP substrate spans the gulf between it and the active site, and we have elucidated the nature of the transition state (TS) and offered an alternate explanation of experimentally observed kinetic isotope effects (KIEs). Our second story focuses on the nature of substrate dynamics in the aromatic prenyltransferase (APTase) protein NphB and how substrate dynamics affects the observed product distribution. Through the examples chosen we show the power of MD and QM/MM methods to provide unique insights into how protein substrate dynamics affects catalytic efficiency. We also illustrate how complex these reactions are and highlight the challenges faced when attempting to design de novo catalysts. While the methods used in our previous studies provided useful insights, several clear challenges still remain. In particular, we have utilized a semiempirical QM model (self-consistent charge density functional tight binding, SCC-DFTB) in our QM/MM studies since the problems we were addressing required extensive sampling. For the problems illustrated, this approach performed admirably (we estimate for these systems an uncertainty of ∼2 kcal/mol), but it is still a semiempirical model, and studies of this type would benefit greatly from more accurate ab initio or DFT models. However, the challenge with these methods is to reach the level of sampling needed to study systems where large conformational changes happen in the many nanoseconds to microsecond time regimes. Hence, how to couple expensive and accurate QM methods with sophisticated sampling algorithms is an important future challenge especially when large-scale studies of catalyst design become of interest. The use of MD and QM/MM models to elucidate enzyme catalytic pathways and to design novel catalytic agents is in its infancy but shows tremendous promise. While this Account summarizes where we have been, we also discuss briefly future directions that improve our fundamental ability to understand enzyme catalysis.

Protein prenylation: enzymes, therapeutics, and biotechnology applications

Protein prenylation is a ubiquitous covalent post-translational modification found in all eukaryotic cells, comprising attachment of either a farnesyl or a geranylgeranyl isoprenoid. It is essential for the proper cellular activity of numerous proteins, including Ras family GTPases and heterotrimeric G-proteins. Inhibition of prenylation has been extensively investigated to suppress the activity of oncogenic Ras proteins to achieve antitumor activity. Here, we review the biochemistry of the prenyltransferase enzymes and numerous isoprenoid analogs synthesized to investigate various aspects of prenylation and prenyltransferases. We also give an account of the current status of prenyltransferase inhibitors as potential therapeutics against several diseases including cancers, progeria, aging, parasitic diseases, and bacterial and viral infections. Finally, we discuss recent progress in utilizing protein prenylation for site-specific protein labeling for various biotechnology applications.

Farnesyltransferase and geranylgeranyltransferase I: structures, mechanism, inhibitors and molecular modeling

Farnesyltransferase (FTase) and geranylgeranyltransferase type I (GGTase-I) have crucial roles in the post-translational modifications of Ras proteins and, therefore, they are promising therapeutic targets for the treatment of various Ras-induced cancers and several other kinds of diseases. In this review, we provide an overview of the structures and biological functions of FTase and GGTase-I. Then, we summarize the typical inhibitors of FTase and GGTase-I, and highlight the drug candidates in clinical trials. In addition, we survey some recent advances in computer-aided drug design (CADD) and molecular modeling studies of FTase and GGTase-I.

Computational insights into function and inhibition of fatty acid amide hydrolase

The Fatty Acid Amide Hydrolase (FAAH) enzyme is a membrane-bound serine hydrolase responsible for the deactivating hydrolysis of a family of naturally occurring fatty acid amides. FAAH is a critical enzyme of the endocannabinoid system, being mainly responsible for regulating the level of its main cannabinoid substrate anandamide. For this reason, pharmacological inhibition of FAAH, which increases the level of endogenous anandamide, is a promising strategy to cure a variety of diseases including pain, inflammation, and cancer. Much structural, mutagenesis, and kinetic data on FAAH has been generated over the last couple of decades. This has prompted several informative computational investigations to elucidate, at the atomic-level, mechanistic details on catalysis and inhibition of this pharmaceutically relevant enzyme. Here, we review how these computational studies - based on classical molecular dynamics, full quantum mechanics, and hybrid QM/MM methods - have clarified the binding and reactivity of some relevant substrates and inhibitors of FAAH. We also discuss the experimental implications of these computational insights, which have provided a thoughtful elucidation of the complex physical and chemical steps of the enzymatic mechanism of FAAH. Finally, we discuss how computations have been helpful for building structure-activity relationships of potent FAAH inhibitors.

Mechanistic insights into Mg2+-independent prenylation by CloQ from classical molecular mechanics and hybrid quantum mechanics/molecular mechanics molecular dynamics simulations

Understanding the mechanism of prenyltransferases is important to the design of engineered proteins capable of synthesizing derivatives of naturally occurring therapeutic agents. CloQ is a Mg(2+)-independent aromatic prenyltransferase (APTase) that transfers a dimethylallyl group to 4-hydroxyphenylpyruvate in the biosynthetic pathway for clorobiocin. APTases consist of a common ABBA fold that defines a β-barrel containing the reaction cavity. Positively charged basic residues line the inside of the β-barrel of CloQ to activate the pyrophosphate leaving group to replace the function of the Mg(2+) cofactor in other APTases. Classical molecular dynamics simulations of CloQ, its E281G and F68S mutants, and the related NovQ were used to explore the binding of the 4-hydroxyphenylpyruvate (4HPP) and dimethylallyl diphosphate substrates in the reactive cavity and the role of various conserved residues. Hybrid quantum mechanics/molecular mechanics potential of mean force (PMF) calculations show that the effect of the replacement of the Mg(2+) cofactor with basic residues yields a similar activation barrier for prenylation to Mg(2+)-dependent APTases like NphB. The topology of the binding pocket for 4HPP is important for selective prenylation at the ortho position of the ring. Methylation at this position alters the conformation of the substrate for O-prenylation at the phenol group. Further, a two-dimensional PMF scan shows that a "reverse" prenylation product may be a possible target for protein engineering.

Prenyltransferase Inhibitors: Treating Human Ailments from Cancer to Parasitic Infections

The posttranslational modification of protein prenylation is a covalent lipid modification on the C-terminus of substrate proteins that serves to enhance membrane affinity. Oncogenic proteins such as Ras have this modification and significant effort has been placed into developing inhibitors of the prenyltransferase enzymes for clinical therapy. In addition to cancer therapy, prenyltransferase inhibitors have begun to find important therapeutic uses in other diseases, including progeria, hepatitis C and D, parasitic infections, and other maladies. This review will trace the evolution of prenyltransferase inhibitors from their initial use as cancer therapeutics to their expanded applications for other diseases.

Catalytic mechanism of aromatic prenylation by NphB

NphB is an aromatic prenyltransferase that catalyzes the attachment of a 10-carbon geranyl group to aromatic substrates. Importantly, NphB exhibits a rich substrate selectivity and product regioselectivity. A systematic computational study has been conducted in order to address several question associated with NphB-catalyzed geranylation. The reaction mechanism of the prenylation step has been characterized as a S(N)1 type dissociative mechanism with a weakly stable carbocation intermediate. A novel π-chamber composed of Tyr121, Tyr216, and 1,6-DHN is found to be important in stabilizing the carbocation. The observed difference in the rates of product formation from 5- and 2-prenylation arises from the differing orientations of the aromatic substrate in the resting state. 4-Prenylation shares the same resting state with 5-prenylation, but the lower free energy barrier for carbocation formation makes the latter reaction more facile. The high free energy barrier associated with 7-prenylation is caused by the unfavorable orientation of 1,6-DHN in active site pocket, along with the difficulty of proton elimination after the prenylation step. A water-mediated proton transfer facilitates the loss of hydrogen at the prenylation site to form the final prenylated product. Interestingly, the same crystallographically observed water molecule has been found to be responsible for proton loss in all three experimentally identified products. After proton transfer, the relaxation of the final product from a sp(3) carbon center to a sp(2) center triggers a "spring-loaded" product release mechanism which pushes the final product out of the binding pocket toward the edge of the active site. The hydrogen bond interactions between the two hydroxyl groups of the aromatic product and the side chains of Ser214 and Tyr288 help to "steer" the movement of the product. In addition, mutagenesis studies identify these same two side chains as being responsible for the observed regioselectivity, particularly 2-prenylation. These observations provide valuable insights into NphB chemistry, offering an opportunity to better engineer the active site and to control the reactivity in order to obtain high yields of the desired product(s). Furthermore, the S(N)1 reaction mechanism observed for NphB differs from the prenylation reaction found in, for example, the farnesyltransferase, which proceeds via an S(N)2-like reaction pathway. The spring-loaded release mechanism highlighted herein also offers novel insights into how enzymes facilitate product release.

Insights into the mechanistic dichotomy of the protein farnesyltransferase peptide substrates CVIM and CVLS

Protein farnesyltransferase (FTase) catalyzes farnesylation of a variety of peptide substrates. (3)H α-secondary kinetic isotope effect (α-SKIE) measurements of two peptide substrates, CVIM and CVLS, are significantly different and have been proposed to reflect a rate-limiting S(N)2-like transition state with dissociative characteristics for CVIM, while, due to the absence of an isotope effect, CVLS was proposed to have a rate-limiting peptide conformational change. Potential of mean force quantum mechanical/molecular mechanical studies coupled with umbrella sampling techniques were performed to further probe this mechanistic dichotomy. We observe the experimentally proposed transition state (TS) for CVIM but find that CVLS has a symmetric S(N)2 TS, which is also consistent with the absence of a (3)H α-SKIE. These calculations demonstrate facile substrate-dependent alterations in the transition state structure catalyzed by FTase.

Finding a needle in the haystack: computational modeling of Mg2+ binding in the active site of protein farnesyltransferase

Studies aimed at elucidating the unknown Mg2+ binding site in protein farnesyltransferase (FTase) are reported. FTase catalyzes the transfer of a farnesyl group to a conserved cysteine residue (Cys1p) on a target protein, an important step for proteins in the signal transduction pathways (e.g., Ras). Mg2+ ions accelerate the protein farnesylation reaction by up to 700-fold. The exact function of Mg2+ in catalysis and the structural characteristics of its binding remain unresolved to date. Molecular dynamics (MD) simulations addressing the role of magnesium ions in FTase are presented, and relevant octahedral binding motifs for Mg2+ in wild-type (WT) FTase and the Dβ352A mutant are explored. Our simulations suggest that the addition of Mg2+ ions causes a conformational change to occur in the FTase active site, breaking interactions known to keep FPP in its inactive conformation. Two relevant Mg2+ ion binding motifs were determined in WT FTase. In the first binding motif, WT1, the Mg2+ ion is coordinated to D352β, zinc-bound D297β, two water molecules, and one oxygen atom from the α- and β-phosphates of farnesyl diphosphate (FPP). The second binding motif, WT2, is identical with the exception of the zinc-bound D297β being replaced by a water molecule in the Mg2+ coordination complex. In the Dβ352A mutant Mg2+ binding motif, D297β, three water molecules, and one oxygen atom from the α- and β-phosphates of FPP complete the octahedral coordination sphere of Mg2+. Simulations of WT FTase, in which Mg2+ was replaced by water in the active site, recreated the salt bridges and hydrogen-bonding patterns around FPP, validating these simulations. In all Mg2+ binding motifs, a key hydrogen bond was identified between a magnesium-bound water and Cys1p, bridging the two metallic binding sites and, thereby, reducing the equilibrium distance between the reacting atoms of FPP Cys1p. The free energy profiles calculated for these systems provide a qualitative understanding of experimental results. They demonstrate that the two reactive atoms approach each other more readily in the presence of Mg2+ in WT FTase and mutant. The flexible WT2 model was found to possess the lowest barrier toward the conformational change, suggesting it is the preferred Mg2+ binding motif in WT FTase. In the mutant, the absence of D352β makes the transition toward a conformational change harder. Our calculations find support for the proposal that D352β performs a critical role in Mg2+ binding and Mg2+ plays an important role in the conformational transition step.

Biosynthesis via carbocations: theoretical studies on terpene formation

This review describes applications of quantum chemical calculations in the field of terpene biosynthesis, with a focus on insights into the mechanisms of terpene-forming carbocation rearrangements arising from theoretical studies.