A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.

Dynamic Origin of Chirality Transfer between Chiral Surface and Achiral Ligand in Au38 Clusters

The transfer of chirality between nanomolecules is at the core of several applications in chiral technology such as sensing and catalysis. However, the origin of this phenomenon and how exactly nanoscale objects transfer chirality to molecules in their vicinity remain largely obscure. Here, we show that the transfer of chirality for the intrinsically chiral gold cluster Au₃₈(SR)₂₄ is site dependent; that is, it differs depending on the ligand-binding sites. This is closely related to the dynamic nature of the ligands on the cluster surface. Using a combination of NMR techniques and molecular dynamics simulations, we could assign the four symmetry-unique ligands on the cluster. The study reveals largely different conformational dynamics of the bound ligands, explaining the diverse diastereotopicities observed for the CH₂ protons of the ligands. Although chirality is a structural property, our study reveals the importance of dynamics for the transfer of chirality.

Molecular-Dynamics-Simulation-Directed Rational Design of Nanoreceptors with Targeted Affinity

Here, we demonstrate the possibility of rationally designing nanoparticle receptors with targeted affinity and selectivity for specific small molecules. We used atomistic molecular-dynamics (MD) simulations to gradually mutate and optimize the chemical structure of the molecules forming the coating monolayer of gold nanoparticles (1.7 nm gold-core size). The MD-directed design resulted in nanoreceptors with a 10-fold improvement in affinity for the target analyte (salicylate) and a 100-fold decrease of the detection limit by NMR-chemosensing from the millimolar to the micromolar range. We could define the exact binding mode, which features prolonged contacts and deep penetration of the guest into the monolayer, as well as a distinct shape of the effective binding pockets characterized by exposed interacting points.

NanoModeler: A Webserver for Molecular Simulations and Engineering of Nanoparticles

Functionalized nanoparticles (NPs) are at the frontier of nanoscience. They hold the promise of innovative applications for human health and technology. In this context, molecular dynamics (MD) simulations of NPs are increasingly employed to understand the fundamental structural and dynamical features of NPs. While informative, such simulations demand a laborious two-step process for their setup. In-house scripts are required to (i) construct complex 3D models of the inner metal core and outer layer of organic ligands, and (ii) correctly assign force-field parameters to these composite systems. Here, we present NanoModeler ( www.nanomodeler.it ), the first Webserver designed to automatically generate and parametrize model systems of monolayer-protected gold NPs and gold nanoclusters. The only required input is a structure file of one or two ligand(s) to be grafted onto the gold core, with the option of specifying homogeneous or heterogeneous NP morphologies. NanoModeler then generates 3D models of the nanosystem and the associated topology files. These files are ready for use with the Gromacs MD engine, and they are compatible with the AMBER family of force fields. We illustrate NanoModeler's capabilities with MD simulations of selected representative NP model systems. NanoModeler is the first platform to automate and standardize the construction and parametrization of realistic models for atomistic simulations of gold NPs and gold nanoclusters.

Abietane-Type Diterpenoids Inhibit Protein Tyrosine Phosphatases by Stabilizing an Inactive Enzyme Conformation

Protein tyrosine phosphatases (PTPs) contribute to a striking variety of human diseases, yet they remain vexingly difficult to inhibit with uncharged, cell-permeable molecules; no inhibitors of PTPs have been approved for clinical use. This study uses a broad set of biophysical analyses to evaluate the use of abietane-type diterpenoids, a biologically active class of phytometabolites with largely nonpolar structures, for the development of pharmaceutically relevant PTP inhibitors. Results of nuclear magnetic resonance analyses, mutational studies, and molecular dynamics simulations indicate that abietic acid can inhibit protein tyrosine phosphatase 1B, a negative regulator of insulin signaling and an elusive drug target, by binding to its active site in a non-substrate-like manner that stabilizes the catalytically essential WPD loop in an inactive conformation; detailed kinetic studies, in turn, show that minor changes in the structures of abietane-type diterpenoids (e.g., the addition of hydrogens) can improve potency (i.e., lower IC₅₀) by 7-fold. These findings elucidate a previously uncharacterized mechanism of diterpenoid-mediated inhibition and suggest, more broadly, that abietane-type diterpenoids are a promising source of structurally diverse-and, intriguingly, microbially synthesizable-molecules on which to base the design of new PTP-inhibiting therapeutics.

Polyamine-containing etoposide derivatives as poisons of human type II topoisomerases: Differential effects on topoisomerase IIα and IIβ

Etoposide is an anticancer drug that acts by inducing topoisomerase II-mediated DNA cleavage. Despite its wide use, etoposide is associated with some very serious side-effects including the development of treatment-related acute myelogenous leukemias. Etoposide targets both human topoisomerase IIα and IIβ. However, the contributions of the two enzyme isoforms to the therapeutic vs. leukemogenic properties of the drug are unclear. In order to develop an etoposide-based drug with specificity for cancer cells that express an active polyamine transport system, the sugar moiety of the drug has been replaced with a polyamine tail. To analyze the effects of this substitution on the specificity of hybrid molecules toward the two enzyme isoforms, we analyzed the activity of a series of etoposide-polyamine hybrids toward human topoisomerase IIα and IIβ. All of the compounds displayed an ability to induce enzyme-mediated DNA cleavage that was comparable to or higher than that of etoposide. Relative to the parent drug, the hybrid compounds displayed substantially higher activity toward topoisomerase IIβ than IIα. Modeling studies suggest that the enhanced specificity may result from interactions with Gln778 in topoisomerase IIβ. The corresponding residue in the α isoform is a methionine.

A Strategically Located Arg/Lys Residue Promotes Correct Base Paring During Nucleic Acid Biosynthesis in Polymerases

Polymerases (Pols) synthesize the double-stranded nucleic acids in the Watson-Crick (W-C) conformation, which is critical for DNA and RNA functioning. Yet, the molecular basis to catalyze the W-C base pairing during Pol-mediated nucleic acids biosynthesis remains unclear. Here, through bioinformatics analyses on a large data set of Pol/DNA structures, we first describe the conserved presence of one positively charged residue (Lys or Arg), which is similarly located near the enzymatic two-metal active site, always interacting directly with the incoming substrate (d)NTP. Incidentally, we noted that some Pol/DNA structures showing the alternative Hoogsteen base pairing were often solved with this specific residue either mutated, displaced, or missing. We then used quantum and classical simulations coupled to free-energy calculations to illustrate how, in human DNA Pol-η, the conserved Arg61 favors W-C base pairing through defined interactions with the incoming nucleotide. Taken together, these structural observations and computational results suggest a structural framework in which this specific residue is critical for stabilizing the incoming (d)NTP nucleotide and base pairing during Pol-mediated nucleic acid biosynthesis. These results may benefit enzyme engineering for nucleic acid processing and encourage new drug discovery strategies to modulate Pols function.

Novel Bacterial Topoisomerase Inhibitors Exploit Asp83 and the Intrinsic Flexibility of the DNA Gyrase Binding Site

DNA gyrases are enzymes that control the topology of DNA in bacteria cells. This is a vital function for bacteria. For this reason, DNA gyrases are targeted by widely used antibiotics such as quinolones. Recently, structural and biochemical investigations identified a new class of DNA gyrase inhibitors called NBTIs (i.e., novel bacterial topoisomerase inhibitors). NBTIs are particularly promising because they are active against multi-drug resistant bacteria, an alarming clinical issue. Structural data recently demonstrated that these NBTIs bind tightly to a newly identified pocket at the dimer interface of the DNA–protein complex. In the present study, we used molecular dynamics (MD) simulations and docking calculations to shed new light on the binding of NBTIs to this site. Interestingly, our MD simulations demonstrate the intrinsic flexibility of this binding site, which allows the pocket to adapt its conformation and form optimal interactions with the ligand. In particular, we examined two ligands, AM8085 and AM8191, which induced a repositioning of a key aspartate (Asp83B), whose side chain can rotate within the binding site. The conformational rearrangement of Asp83B allows the formation of a newly identified H-bond interaction with an NH on the bound NBTI, which seems important for the binding of NBTIs having such functionality. We validated these findings through docking calculations using an extended set of cognate oxabicyclooctane-linked NBTIs derivatives (~150, in total), screened against multiple target conformations. The newly identified H-bond interaction significantly improves the docking enrichment. These insights could be helpful for future virtual screening campaigns against DNA gyrase.

Allosteric Communication Networks in Proteins Revealed through Pocket Crosstalk Analysis

The detection and characterization of binding pockets and allosteric communication in proteins is crucial for studying biological regulation and performing drug design. Nowadays, ever-longer molecular dynamics (MD) simulations are routinely used to investigate the spatiotemporal evolution of proteins. Yet, there is no computational tool that can automatically detect all the pockets and potential allosteric communication networks along these extended MD simulations. Here, we use a novel and fully automated algorithm that examines pocket formation, dynamics, and allosteric communication embedded in microsecond-long MD simulations of three pharmaceutically relevant proteins, namely, PNP, A2A, and Abl kinase. This dynamic analysis uses pocket crosstalk, defined as the temporal exchange of atoms between adjacent pockets, along the MD trajectories as a fingerprint of hidden allosteric communication networks. Importantly, this study indicates that dynamic pocket crosstalk analysis provides new mechanistic understandings on allosteric communication networks, enriching the available experimental data. Thus, our results suggest the prospective use of this unprecedented dynamic analysis to characterize transient binding pockets for structure-based drug design.

Nanoparticle-Based Receptors Mimic Protein-Ligand Recognition

The self-assembly of a monolayer of ligands on the surface of noble-metal nanoparticles dictates the fundamental nanoparticle's behavior and its functionality. In this combined computational-experimental study, we analyze the structure, organization, and dynamics of functionalized coating thiols in monolayer-protected gold nanoparticles (AuNPs). We explain how functionalized coating thiols self-organize through a delicate and somehow counterintuitive balance of interactions within the monolayer itself and with the solvent. We further describe how the nature and plasticity of these interactions modulate nanoparticle-based chemosensing. Importantly, we found that self-organization of coating thiols can induce the formation of binding pockets in AuNPs. These transient cavities can accommodate small molecules, mimicking protein-ligand recognition, which could explain the selectivity and sensitivity observed for different organic analytes in NMR chemosensing experiments. Thus, our findings advocate for the rational design of tailored coating groups to form specific recognition binding sites on monolayer-protected AuNPs.

•

Synthesis and molecular simulations of AuNPs for chemosensing

•

A rationale for the molecular recognition ability of functionalized AuNPs

•

Functionalized coating ligands form transient protein-like binding pockets

•

Toward the computational nanodesign of intelligent nanoreceptors for chemosensing

The functionalization of monolayer-protected nanoparticles is at the frontier of nanotechnology, such that innovative applications are emerging in fields such as nanomedicine, chemosensing, and even catalysis. Importantly, the nanoparticle's functionality is mainly defined by the nature of the ligands forming the coating monolayer. Here, we show how the self-organization of functionalized coating ligands in monolayer-protected gold nanoparticles (AuNPs) affects their solubility and molecular recognition abilities. We found that coating ligands form transient, protein-like binding pockets in functionalized AuNPs. Thus, we reveal that nanoparticle-based chemosensing operates through a recognition process that is similar to that for protein-ligand complex formation. These findings could now herald the arrival of the computational nanodesign of intelligent nanodevices with recognition abilities toward small molecules such as drugs, metabolites, illegal drugs, and small molecular markers for cancer.

Pharmacophore Identification and Scaffold Exploration to Discover Novel, Potent, and Chemically Stable Inhibitors of Acid Ceramidase in Melanoma Cells

Acid ceramidase (AC) hydrolyzes ceramides, which are central lipid messengers for metabolism and signaling of sphingolipids. A growing body of evidence links deregulation of sphingolipids to several diseases, including cancer. Indeed, AC expression is abnormally high in melanoma cells. AC inhibition may thus be key to treating malignant melanoma. Here, we have used a systematic scaffold exploration to design a general pharmacophore for AC inhibition. This pharmacophore comprises a 6 + 5 fused ring heterocycle linked to an aliphatic substituent via a urea moiety. We have thus identified the novel benzimidazole derivatives 10, 21, 27, and 30, which are highly potent AC inhibitors. Their chemical and metabolic stabilities are comparable or superior to those of previously reported AC inhibitors. Moreover, they are potent against endogenous AC in intact melanoma cells. These novel inhibitors merit further characterization and can serve as a promising starting point for the discovery of new antimelanoma therapeutics.

Lid domain plasticity and lipid flexibility modulate enzyme specificity in human monoacylglycerol lipase

Human monoacylglycerol lipase (MAGL) is a membrane-interacting enzyme that generates pro-inflammatory signaling molecules. For this reason, MAGL inhibition is a promising strategy to treat pain, cancer, and neuroinflammatory diseases. MAGL can hydrolyze monoacylglycerols bearing an acyl chain of different lengths and degrees of unsaturation, cleaving primarily the endocannabinoid 2-arachidonoylglycerol. Importantly, the enzymatic binding site of MAGL is confined by a 75-amino-acid-long, flexible cap domain, named 'lid domain', which is structurally similar to that found in several other lipases. However, it is unclear how lid domain plasticity affects catalysis in MAGL. By integrating extensive molecular dynamics simulations and free-energy calculations with mutagenesis and kinetic experiments, we here define a lid-domain-mediated mechanism for substrate selection and binding in MAGL catalysis. In particular, we clarify the key role of Phe159 and Ile179, two conserved residues within the lid domain, in regulating substrate specificity in MAGL. We conclude by proposing that other structurally related lipases may share this lid-domain-mediated mechanism for substrate specificity.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Carnosine and Homocarnosine Degradation Mechanisms by the Human Carnosinase Enzyme CN1: Insights from Multiscale Simulations

The endogenous dipeptide l-carnosine, and its derivative homocarnosine, prevent and reduce several pathologies like amytrophic lateral sclerosis (ALS), Alzheimer's disease, and Parkinson's disease. Their beneficial action is severely hampered because of the hydrolysis by carnosinase enzymes, in particular the human carnosinase, hCN1. This belongs to the metallopeptidase M20 family, where a cocatalytic active site is formed by two Zn(2+) ions, bridged by a hydroxide anion. The protein may exist as a monomer and as a dimer in vivo. Here we used hybrid quantum mechanics/molecular mechanics simulations based on the dimeric apoenzyme's structural information to predict the Michaelis complexes with l-carnosine and its derivative homocarnosine. On the basis of our calculations, we suggest that (i) l-carnosine degradation occurs through a nucleophilic attack of a Zn(2+)-coordinated bridging moiety for both monomer and dimer. This mechanistic hypothesis for hCN1 catalysis differs from previous proposals, while it is in agreement with available experimental data. (ii) The experimentally measured higher affinity of homocarnosine for the enzyme relative to l-carnosine might be explained, at least in part, by more extensive interactions inside the monomeric and dimeric hCN1's active site. (iii) Hydrogen bonds at the binding site, present in the dimer but absent in the monomer, might play a role in the experimentally observed higher activity of the dimeric form. Investigations of the enzymatic reaction are required to establish or disprove this hypothesis. Our results may serve as a basis for the design of potent hCN1 inhibitors.

Cooperative motion of a key positively charged residue and metal ions for DNA replication catalyzed by human DNA Polymerase-η

Trans-lesion synthesis polymerases, like DNA Polymerase-η (Pol-η), are essential for cell survival. Pol-η bypasses ultraviolet-induced DNA damages via a two-metal-ion mechanism that assures DNA strand elongation, with formation of the leaving group pyrophosphate (PPi). Recent structural and kinetics studies have shown that Pol-η function depends on the highly flexible and conserved Arg61 and, intriguingly, on a transient third ion resolved at the catalytic site, as lately observed in other nucleic acid-processing metalloenzymes. How these conserved structural features facilitate DNA replication, however, is still poorly understood. Through extended molecular dynamics and free energy simulations, we unravel a highly cooperative and dynamic mechanism for DNA elongation and repair, which is here described by an equilibrium ensemble of structures that connect the reactants to the products in Pol-η catalysis. We reveal that specific conformations of Arg61 help facilitate the recruitment of the incoming base and favor the proper formation of a pre-reactive complex in Pol-η for efficient DNA editing. Also, we show that a third transient metal ion, which acts concertedly with Arg61, serves as an exit shuttle for the leaving PPi. Finally, we discuss how this effective and cooperative mechanism for DNA repair may be shared by other DNA-repairing polymerases.

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

The protein farnesyltransferase (FTase) is a Zn(2+)-metalloenzyme that catalyzes the farnesylation reaction, i.e., the transfer of the 15-carbon atom farnesyl group from farnesyl diphosphate (FPP) to a specific cysteine of protein substrates. Oncogenic Ras proteins, which are among the FTase substrates, are observed in about 20-30% of human cancer cells. Thus, FTase represents a target for anticancer drug design. Herein, we present a classical force-field-based and quantum mechanics/molecular mechanics (QM/MM) computational study of the FTase reaction mechanism. Our findings offer a detailed picture of the FTase catalytic pathway, describing structural features and the energetics of its saddle points. A moderate dissociation of the diphosphate group from the FPP is observed during the nucleophilic attack of the zinc-bound thiolate. At the transition state, a resonance structure is observed, which indicates the formation of a metastable carbocation. However, no stable intermediate is found along the reaction pathway. Thus, the reaction occurs via an associative mechanism with dissociative character, in agreement with the mechanism proposed by Fierke et al. ( Biochemistry 2000, 39, 2593-2602 and Biochemistry 2003, 42, 9741-9748 ). Moreover, a fluorine-substituted FPP analogue (CF3-FPP) is used to investigate the inhibitory effect of fluorine, which in turn provides additional agreement with experimental data.

The Molecular Basis for Dual Fatty Acid Amide Hydrolase (FAAH)/Cyclooxygenase (COX) Inhibition

The design of multitarget‐directed ligands is a promising strategy for discovering innovative drugs. Here, we report a mechanistic study that clarifies key aspects of the dual inhibition of the fatty acid amide hydrolase (FAAH) and the cyclooxygenase (COX) enzymes by a new multitarget‐directed ligand named ARN2508 (2‐[3‐fluoro‐4‐[3‐(hexylcarbamoyloxy)phenyl]phenyl]propanoic acid). This potent dual inhibitor combines, in a single scaffold, the pharmacophoric elements often needed to block FAAH and COX, that is, a carbamate moiety and the 2‐arylpropionic acid functionality, respectively. Molecular modeling and molecular dynamics simulations suggest that ARN2508 uses a noncovalent mechanism of inhibition to block COXs, while inhibiting FAAH via the acetylation of the catalytic Ser241, in line with previous experimental evidence for covalent FAAH inhibition. This study proposes the molecular basis for the dual FAAH/COX inhibition by this novel hybrid scaffold, stimulating further experimental studies and offering new insights for the rational design of novel anti‐inflammatory agents that simultaneously act on FAAH and COX.

An optimized polyamine moiety boosts the potency of human type II topoisomerase poisons as quantified by comparative analysis centered on the clinical candidate F14512

Combined computational-experimental analyses explain and quantify the spermine-vectorized F14512's boosted potency as a topoII poison. We found that an optimized polyamine moiety boosts drug binding to the topoII/DNA cleavage complex, rather than to the DNA alone. These results provide new structural bases and key reference data for designing new human topoII poisons.

Keys to Lipid Selection in Fatty Acid Amide Hydrolase Catalysis: Structural Flexibility, Gating Residues and Multiple Binding Pockets

The fatty acid amide hydrolase (FAAH) regulates the endocannabinoid system cleaving primarily the lipid messenger anandamide. FAAH has been well characterized over the years and, importantly, it represents a promising drug target to treat several diseases, including inflammatory-related diseases and cancer. But its enzymatic mechanism for lipid selection to specifically hydrolyze anandamide, rather than similar bioactive lipids, remains elusive. Here, we clarify this mechanism in FAAH, examining the role of the dynamic paddle, which is formed by the gating residues Phe432 and Trp531 at the boundary between two cavities that form the FAAH catalytic site (the “membrane-access” and the “acyl chain-binding” pockets). We integrate microsecond-long MD simulations of wild type and double mutant model systems (Phe432Ala and Trp531Ala) of FAAH, embedded in a realistic membrane/water environment, with mutagenesis and kinetic experiments. We comparatively analyze three fatty acid substrates with different hydrolysis rates (anandamide > oleamide > palmitoylethanolamide). Our findings identify FAAH’s mechanism to selectively accommodate anandamide into a multi-pocket binding site, and to properly orient the substrate in pre-reactive conformations for efficient hydrolysis that is interceded by the dynamic paddle. Our findings therefore endorse a structural framework for a lipid selection mechanism mediated by structural flexibility and gating residues between multiple binding cavities, as found in FAAH. Based on the available structural data, this exquisite catalytic strategy for substrate specificity seems to be shared by other lipid-degrading enzymes with similar enzymatic architecture. The mechanistic insights for lipid selection might assist de-novo enzyme design or drug discovery efforts.

We describe a new structural enzymatic framework to regulate substrate specificity in lipid-degrading enzymes such as fatty acid amide hydrolase (FAAH), a key enzyme for the endocannabinoid lipid signaling that hydrolyzes a variety of lipids, however with different catalytic rates. The identified novel mechanism and key features for lipid selection in FAAH are then analysed in the context of other relevant lipid-degrading enzymes. Through the integration of microsecond-long molecular dynamics simulations with mutagenesis and kinetic experiments, our study suggests that structural flexibility, gating residues and multiple cavities in one catalytic site are keys to lipid selection in the endocannabinoid system. Our results suggest that the structural framework proposed here could likely be a general enzymatic strategy of other lipid-degrading enzymes to select the preferred lipid substrate within a broad spectrum of biologically active lipids. This new, and likely general, structural framework for lipid selection in FAAH could therefore now encourage additional experimental verifications of the role of ligand and structural flexibility, as regulated by key gating residues at the boundaries of multiple cavities forming a single catalytic site, as observed in several other lipid-degrading enzymes.

Catalytic metal ions and enzymatic processing of DNA and RNA

CONSPECTUS: Two-metal-ion-dependent nucleases cleave the phosphodiester bonds of nucleic acids via the two-metal-ion (2M) mechanism. Several high-resolution X-ray structures portraying the two-metal-aided catalytic site, together with mutagenesis and kinetics studies, have demonstrated a functional role of the ions for catalysis in numerous metallonucleases. Overall, the experimental data confirm the general mechanistic hypothesis for 2M-aided phosphoryl transfer originally reported by Steitz and Steitz ( Proc. Natl. Acad. Sci. U.S.A. 1993 , 90 ( 14 ), 6498 - 6502 ). This seminal paper proposed that one metal ion favors the formation of the nucleophile, while the nearby second metal ion facilitates leaving group departure during RNA hydrolysis. Both metals were suggested to stabilize the enzymatic transition state. Nevertheless, static X-ray structures alone cannot exhaustively unravel how the two ions execute their functional role along the enzymatic reaction during processing of DNA or RNA strands when moving from reactants to products, passing through metastable intermediates and high-energy transition states. In this Account, we discuss the role of multiscale molecular simulations in further disclosing mechanistic insights of 2M-aided catalysis for two prototypical enzymatic targets for drug discovery, namely, ribonuclease H (RNase H) and type II topoisomerase (topoII). In both examples, first-principles molecular simulations, integrated with structural data, emphasize a cooperative motion of the bimetal motif during catalysis. The coordinated motion of both ions is crucial for maintaining a flexible metal-centered structural architecture exquisitely tailored to accommodate the DNA or RNA sugar-phosphate backbone during phosphodiester bond cleavage. Furthermore, our analysis of RNase H and the N-terminal domain (PAN) of influenza polymerase shows that classical molecular dynamics simulations coupled with enhanced sampling techniques have contributed to describe the modulatory effect of metal ion concentration and metal uptake on the 2M mechanism and efficiency. These aspects all point to the emerging and intriguing role of additional adjacent ions potentially involved in the modulation of phosphoryl transfer reactions and enzymatic turnover in 2M-catalysis, as recently observed experimentally in polymerase η and homing endonuclease I-DmoI. These computational results, integrated with experimental findings, describe and reinforce the nascent concept of a functional and cooperative dynamics of the catalytic metal ions during the 2M-dependent enzymatic processing of DNA and RNA. Encouraged by the insights provided by computational approaches, we foresee further experiments that will feature the functional and joint dynamics of the catalytic metal ions for nucleic acid processing. This could impact the de novo design of artificial metallonucleases and the rational design of potent metal-chelating inhibitors of pharmaceutically relevant enzymes.

Benzoxazolone carboxamides: potent and systemically active inhibitors of intracellular acid ceramidase

The ceramides are a family of bioactive lipid-derived messengers involved in the control of cellular senescence, inflammation, and apoptosis. Ceramide hydrolysis by acid ceramidase (AC) stops the biological activity of these substances and influences survival and function of normal and neoplastic cells. Because of its central role in the ceramide metabolism, AC may offer a novel molecular target in disorders with dysfunctional ceramide-mediated signaling. Here, a class of benzoxazolone carboxamides is identified as the first potent and systemically active inhibitors of AC. Prototype members of this class inhibit AC with low nanomolar potency by covalent binding to the catalytic cysteine. Their metabolic stability and high in vivo efficacy suggest that these compounds may be used as probes to investigate the roles of ceramide in health and disease, and that this scaffold may represent a promising starting point for the development of novel therapeutic agents.

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.

Wagging the Tail: Essential Role of Substrate Flexibility in FAAH Catalysis

The serine hydrolase, fatty acid amide hydrolase (FAAH), is responsible for the intracellular degradation of anandamide and other bioactive fatty acid ethanolamides involved in the regulation of pain, inflammation, and other pathophysiological processes. The catalytic site of FAAH is composed of multiple cavities with mixed hydrophobic and hydrophilic properties, the role of which remains incompletely understood. Anandamide is thought to enter the active site through a "membrane-access" (MA) channel and position its flexible fatty acyl chain in a highly hydrophobic "acyl chain-binding" (AB) cavity to allow for hydrolysis to occur. Using microsecond molecular dynamics (MD) simulations of FAAH embedded in a realistic membrane/water environment, we show now that anandamide may not lock itself into the AB cavity but may rather assume catalytically significant conformations required for hydrolysis by moving its flexible arachidonoyl tail between the MA and AB cavities. This process is regulated by a phenylalanine residue (Phe432) located at the boundary between the two cavities, which may act as a "dynamic paddle." The results identify structural flexibility as a key determinant by which FAAH recognizes its primary lipid substrate.

A binding site for nonsteroidal anti-inflammatory drugs in fatty acid amide hydrolase

In addition to inhibiting the cyclooxygenase (COX)-mediated biosynthesis of prostanoids, various widely used nonsteroidal anti-inflammatory drugs (NSAIDs) enhance endocannabinoid signaling by blocking the anandamide-degrading membrane enzyme fatty acid amide hydrolase (FAAH). The X-ray structure of FAAH in complex with the NSAID carprofen, along with site-directed mutagenesis, enzyme activity assays, and NMR analysis, has revealed the molecular details of this interaction, providing information that may guide the design of dual FAAH-COX inhibitors with superior analgesic efficacy.

Identification and characterization of carprofen as a multitarget fatty acid amide hydrolase/cyclooxygenase inhibitor

Pain and inflammation are major therapeutic areas for drug discovery. Current drugs for these pathologies have limited efficacy, however, and often cause a number of unwanted side effects. In the present study, we identify the nonsteroidal anti-inflammatory drug carprofen as a multitarget-directed ligand that simultaneously inhibits cyclooxygenase-1 (COX-1), COX-2, and fatty acid amide hydrolase (FAAH). Additionally, we synthesized and tested several derivatives of carprofen, sharing this multitarget activity. This may result in improved analgesic efficacy and reduced side effects (Naidu et al. J. Pharmacol. Exp. Ther.2009, 329, 48-56; Fowler, C. J.; et al. J. Enzyme Inhib. Med. Chem.2012, in press; Sasso et al. Pharmacol. Res.2012, 65, 553). The new compounds are among the most potent multitarget FAAH/COX inhibitors reported so far in the literature and thus may represent promising starting points for the discovery of new analgesic and anti-inflammatory drugs.

A catalytic mechanism for cysteine N-terminal nucleophile hydrolases, as revealed by free energy simulations

The N-terminal nucleophile (Ntn) hydrolases are a superfamily of enzymes specialized in the hydrolytic cleavage of amide bonds. Even though several members of this family are emerging as innovative drug targets for cancer, inflammation, and pain, the processes through which they catalyze amide hydrolysis remains poorly understood. In particular, the catalytic reactions of cysteine Ntn-hydrolases have never been investigated from a mechanistic point of view. In the present study, we used free energy simulations in the quantum mechanics/molecular mechanics framework to determine the reaction mechanism of amide hydrolysis catalyzed by the prototypical cysteine Ntn-hydrolase, conjugated bile acid hydrolase (CBAH). The computational analyses, which were confirmed in water and using different CBAH mutants, revealed the existence of a chair-like transition state, which might be one of the specific features of the catalytic cycle of Ntn-hydrolases. Our results offer new insights on Ntn-mediated hydrolysis and suggest possible strategies for the creation of therapeutically useful inhibitors.

Covalent inhibitors of fatty acid amide hydrolase: a rationale for the activity of piperidine and piperazine aryl ureas

Recently, covalent drugs have attracted great interest in the drug discovery community, with successful examples that have demonstrated their therapeutic effects. Here, we focus on the covalent inhibition of the fatty acid amide hydrolase (FAAH), which is a promising strategy in the treatment of pain and inflammation. Among the most recent and potent FAAH inhibitors (FAAHi), there are the cyclic piperidine and piperazine aryl ureas. FAAH hydrolyzes efficiently the amide bond of these compounds, forming a covalent enzyme-inhibitor adduct. To rationalize this experimental evidence, we performed an extensive computational analysis centered on piperidine-based PF750 (1) and piperazine-based JNJ1661010 (2), two potent lead compounds used to generate covalent inhibitors as clinical candidates. We found that FAAH induces a distortion of the amide bond of the piperidine and piperazine aryl ureas. Quantum mechanics/molecular mechanics ΔE(LUMO-HOMO) energies indicate that the observed enzyme-induced distortion of the amide bond favors the formation of a covalent FAAH-inhibitor adduct. These findings could help in the rational structure-based design of novel covalent FAAHi.

Understanding the effect of magnesium ion concentration on the catalytic activity of ribonuclease H through computation: does a third metal binding site modulate endonuclease catalysis?

Ribonuclease H (RNase H) belongs to the nucleotidyl-transferase superfamily and hydrolyzes the phosphodiester linkage on the RNA strand of a DNA/RNA hybrid duplex. Due to its activity in HIV reverse transcription, it represents a promising target for anti-HIV drug design. While crystallographic data have located two ions in the catalytic site, there is ongoing debate concerning just how many metal ions bound at the active site are optimal for catalysis. Indeed, experiments have shown a dependency of the catalytic activity on the Mg(2+) concentration. Moreover, in RNase H, the glutamate residue E188 has been shown to be essential for full enzymatic activation, regardless of the Mg(2+) concentration. The catalytic center is known to contain two Mg(2+) ions, and E188 is not one of the primary metal ligands. Herein, classical molecular dynamics (MD) simulations are employed to study the metal-ligand coordination in RNase H at different concentration of Mg(2+). Importantly, the presence of a third Mg(2+) ion, bound to the second-shell ligand E188, is a persistent feature of the MD simulations. Free energy calculations have identified two distinct conformations, depending on the concentration of Mg(2+). At standard concentration, a third Mg(2+) is found in the catalytic pocket, but it does not perturb the optimal RNase H active conformation. However, at higher concentration, the third Mg(2+) ion heavily perturbs the nucleophilic water and thereby influences the catalytic efficiency of RNase H. In addition, the E188A mutant shows no ability to engage additional Mg(2+) ions near the catalytic pocket. This finding likely explains the decrease in catalytic activity of E188A and also supports the key role of E188 in localizing the third Mg(2+) ion at the active site. Glutamate residues are commonly found surrounding the metal center in the endonuclease family, which suggests that this structural motif may be an important feature to enhance catalytic activity. The present MD calculations support the hypothesis that RNase H can accommodate three divalent metal ions in its catalytic pocket and provide an in-depth understanding of their dynamic role for catalysis.

Bridging quantum mechanics and structure-based drug design

The last decade has seen great advances in the use of quantum mechanics (QM) to solve biological problems of pharmaceutical relevance. For instance, enzymatic catalysis is often investigated by means of the so-called QM/MM approach, which uses QM and molecular mechanics (MM) methods to determine the (free) energy landscape of the enzymatic reaction mechanism. Here, I will discuss a few representative examples of QM and QM/MM studies of important metalloenzymes of pharmaceutical interest (i.e. metallophosphatases and metallo-beta-lactamases). This review article aims to show how QM-based methods can be used to elucidate ligand-receptor interactions. The challenge is then to exploit this knowledge for the structure-based design of new and potent inhibitors, such as transition state (TS) analogues that resemble the structure and physicochemical properties of the enzymatic TS. Given the results and potential expressed to date by QM-based methods in studying biological problems, the application of QM in structure-based drug design will likely increase, making of these once-prohibitive computations a routinely used tool for drug design.

Phosphodiester cleavage in ribonuclease H occurs via an associative two-metal-aided catalytic mechanism

Ribonuclease H (RNase H) belongs to the nucleotidyl-transferase (NT) superfamily and hydrolyzes the phosphodiester linkages that form the backbone of the RNA strand in RNA x DNA hybrids. This enzyme is implicated in replication initiation and DNA topology restoration and represents a very promising target for anti-HIV drug design. Structural information has been provided by high-resolution crystal structures of the complex RNase H/RNA x DNA from Bacillus halodurans (Bh), which reveals that two metal ions are required for formation of a catalytic active complex. Here, we use classical force field-based and quantum mechanics/molecular mechanics calculations for modeling the nucleotidyl transfer reaction in RNase H, clarifying the role of the metal ions and the nature of the nucleophile (water versus hydroxide ion). During the catalysis, the two metal ions act cooperatively, facilitating nucleophile formation and stabilizing both transition state and leaving group. Importantly, the two Mg(2+) metals also support the formation of a meta-stable phosphorane intermediate along the reaction, which resembles the phosphorane intermediate structure obtained only in the debated beta-phosphoglucomutase crystal (Lahiri, S. D.; et al. Science 2003, 299 (5615), 2067-2071). The nucleophile formation (i.e., water deprotonation) can be achieved in situ, after migration of one proton from the water to the scissile phosphate in the transition state. This proton transfer is actually mediated by solvation water molecules. Due to the highly conserved nature of the enzymatic bimetal motif, these results might also be relevant for structurally similar enzymes belonging to the NT superfamily.

Computational Study of the Phosphoryl Transfer Catalyzed by a Cyclin-Dependent Kinase

A cyclin-dependent kinase, Cdk2, catalyzes the transfer of the gamma-phosphate from ATP to a threonine or serine residue of its polypeptide substrates. Here, we investigate aspects of the reaction mechanism of Cdk2 by gas-phase density functional calculations, classical molecular dynamics, and Car-Parrinello QM/MM simulations. We focus on the role of the conserved Asp127 and on the nature of the phosphoryl transfer reaction mechanism catalyzed by Cdk2. Our findings suggest that Asp127 is active in its deprotonated form by assisting the formation of the near-attack orientation of the substrate serine or threonine. Therefore, the residue does not act as a general base during the catalysis. The mechanism for the phosphoryl transfer is a single SN2-like concerted step, which shows a phosphorane-like transition state geometry. Although the resulting reaction mechanism is in agreement with a previous density functional study of the same catalytic reaction mechanism (Cavalli et al., Chem. Comm. 2003, 1308-1309), the reaction barrier is considerably lower when QM/MM calculations are performed, as in this study ( approximately 42 kcal mol(-1) QM vs. approximately 24 kcal mol(-1) QM/MM); this indicates that important roles for the catalysis are played by the protein environment and solvent waters. Because of the high amino acid sequence conservation among the whole family of cyclin-dependent kinases (CDKs), these results could be general for the CDK family.

Proton shuttles and phosphatase activity in soluble epoxide hydrolase

Recently, a novel metal Mg2+-dependent phosphatase activity has been discovered in the N-terminal domain of the soluble epoxide hydrolase (sEH), opening a new branch of fatty acid metabolism and providing an additional site for drug targeting. Importantly, the sEH N-terminal fold belongs to the haloacid dehalogenase (HAD) superfamily, which comprises a vast majority of phosphotransferases. Herein, we present the results of a computational study of the sEH phosphatase activity, which includes classical molecular dynamics (MD) simulations and mixed quantum mechanical/molecular mechanics (QM/MM) calculations. On the basis of experimental results, a two-step mechanism has been proposed and herein investigated: (1) phosphoenzyme intermediate formation and (2) phosphoenzyme intermediate hydrolysis. Building on our earlier work, we now provide a detailed description of the reaction mechanism for the whole catalytic cycle along with its free energy profile. The present computations suggest metaphosphate-like transition states for these phosphoryl transfers. They also reveal that the enzyme promotes water deprotonation and facilitates shuttling of protons via a metal-ligand connecting water bridge (WB). These WB-mediated proton shuttles are crucial for the activation of the solvent nucleophile and for the stabilization of the leaving group. Moreover, due to the conservation of structural features in the N-terminal catalytic site of sEH and other members of the HAD superfamily, we suggest a generalization of our findings to these other metal-dependent phosphatases.

Modeling the charge distribution at metal sites in proteins for molecular dynamics simulations

Almost half of the proteome of living organisms is constituted of metalloproteins. Unfortunately, the ability of the current generation of molecular dynamics pairwise-additive forcefields to properly describe metal pockets is severely lacking due to the intrinsic difficulty of handling polarization and charge transfer contributions. In order to improve the description of metalloproteins, a simple reparameterization strategy is proposed herein that does not involve artificial constraints. Specifically, a non-bonded quantum mechanical-based model is used to capture the mean polarization and charge transfer contributions to the interatomic forces within the metal site. The present approach is demonstrated to provide enough accuracy to maintain the integrity of the metal pocket for a variety of metalloproteins during extended (multi-nanosecond) molecular dynamics simulations. The method enables the sampling of small conformational changes and the relaxation of local frustrations in NMR structures.

Metadynamics as a tool for exploring free energy landscapes of chemical reactions

The metadynamics or hills method is a relatively new molecular dynamics technique aimed to enhance the sampling of separated regions in phase space and map out the underlying free energy landscape as a function of a small number of order parameters or collective variables. The high efficiency allows for the application of metadynamics in combination with first principles dynamics methods, in particular with Car-Parrinello molecular dynamics, to study processes in which changes in the electronic structure play a dominant role, such as chemical reactions. The option to choose several independent collective variables is important to tackle complex and concerted transformations that lack an obvious a priori choice for a single reaction coordinate. In this Account, we discuss the role of metadynamics in the search of transition states, local minima, reaction paths, free energy profiles, and reaction coordinates among a growing list of alternative methods.

Role of phosphorylated Thr160 for the activation of the CDK2/Cyclin A complex

The enzymatic activity of the CDK2/Cyclin A complex increases upon the specific phosphorylation of Thr160@CDK2. In the present study, we have performed a comparative molecular dynamics (MD) study of models of the complex CDK2/Cyclin A/Substrate, which differ for the presence or absence of the phosphate group bound to Thr160. The models are based on two X-ray structures available for CDK2/CyclinA and pCDK2/CyclinA/Substrate complexes. In this way, we analyze the influence of the phosphorylated Thr160 (pThr160) on both the flexibility of CDK2 activation loop (AL) and substrate binding in CDK2. Our calculations point to a decreased flexibility of the AL in the phosphorylated model, in fairly good agreement with experimental data, and to a key role of pThr160 for substrate recognition and stability. Multiple alignments of the CDKs sequences point to the very high conservation of the AL sequence among the CDKs, thus extending our results to all CDKs.