Catalytic metal ions and enzymatic processing of DNA and RNA

Striking Plasticity of CRISPR-Cas9 and Key Role of Non-target DNA, as Revealed by Molecular Simulations

The CRISPR (clustered regularly interspaced short palindromic repeats)-Cas9 system recently emerged as a transformative genome-editing technology that is innovating basic bioscience and applied medicine and biotechnology. The endonuclease Cas9 associates with a guide RNA to match and cleave complementary sequences in double stranded DNA, forming an RNA:DNA hybrid and a displaced non-target DNA strand. Although extensive structural studies are ongoing, the conformational dynamics of Cas9 and its interplay with the nucleic acids during association and DNA cleavage are largely unclear. Here, by employing multi-microsecond time scale molecular dynamics, we reveal the conformational plasticity of Cas9 and identify key determinants that allow its large-scale conformational changes during nucleic acid binding and processing. We show how the "closure" of the protein, which accompanies nucleic acid binding, fundamentally relies on highly coupled and specific motions of the protein domains, collectively initiating the prominent conformational changes needed for nucleic acid association. We further reveal a key role of the non-target DNA during the process of activation of the nuclease HNH domain, showing how the nontarget DNA positioning triggers local conformational changes that favor the formation of a catalytically competent Cas9. Finally, a remarkable conformational plasticity is identified as an intrinsic property of the HNH domain, constituting a necessary element that allows for the HNH repositioning. These novel findings constitute a reference for future experimental studies aimed at a full characterization of the dynamic features of the CRISPR-Cas9 system, and-more importantly-call for novel structure engineering efforts that are of fundamental importance for the rational design of new genome-engineering applications.

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.

Exploring the mechanism of DNA polymerases by analyzing the effect of mutations of active site acidic groups in Polymerase β

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase β as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc.

Conformational States of HIV-1 Reverse Transcriptase for Nucleotide Incorporation vs Pyrophosphorolysis-Binding of Foscarnet

HIV-1 reverse transcriptase (RT) catalytically incorporates individual nucleotides into a viral DNA strand complementing an RNA or DNA template strand; the polymerase active site of RT adopts multiple conformational and structural states while performing this task. The states associated are dNTP binding at the N site, catalytic incorporation of a nucleotide, release of a pyrophosphate, and translocation of the primer 3'-end to the P site. Structural characterization of each of these states may help in understanding the molecular mechanisms of drug activity and resistance and in developing new RT inhibitors. Using a 38-mer DNA template-primer aptamer as the substrate mimic, we crystallized an RT/dsDNA complex that is catalytically active, yet translocation-incompetent in crystals. The ability of RT to perform dNTP binding and incorporation in crystals permitted obtaining a series of structures: (I) RT/DNA (P-site), (II) RT/DNA/AZTTP ternary, (III) RT/AZT-terminated DNA (N-site), and (IV) RT/AZT-terminated DNA (N-site)/foscarnet complexes. The stable N-site complex permitted the binding of foscarnet as a pyrophosphate mimic. The Mg(2+) ions dissociated after catalytic addition of AZTMP in the pretranslocated structure III, whereas ions A and B had re-entered the active site to bind foscarnet in structure IV. The binding of foscarnet involves chelation with the Mg(2+) (B) ion and interactions with K65 and R72. The analysis of interactions of foscarnet and the recently discovered nucleotide-competing RT inhibitor (NcRTI) α-T-CNP in two different conformational states of the enzyme provides insights for developing new classes of polymerase active site RT inhibitors.

Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand

A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. It was predicted that Tm1631 binds with the DNA and that the Tm1631–DNA complex is an endonuclease repair system with a DNA repair function (Konc et al. PLoS Comput Biol 9(11): e1003341, 2013). We observed that the severely bent, strained DNA binds to the protein for the entire 90 ns of classical molecular dynamics (MD) performed; we could observe no significant changes in the most distorted region of the DNA, where the cleavage of phosphodiester bond occurs. In this article, we modeled the reaction mechanism at the interface between Tm1631 and its proposed ligand, the DNA molecule, focusing on cleavage of the phosphodiester bond. After addition of two Mg²⁺ ions to the reaction center and extension of classical MD by 50 ns (totaling 140 ns), the DNA ligand stayed bolted to the protein. Results from density functional theory quantum mechanics/molecular mechanics (QM/MM) calculations suggest that the reaction is analogous to known endonuclease mechanisms: an enzyme reaction mechanism with two Mg²⁺ ions in the reaction center and a pentacovalent intermediate. The minimum energy pathway profile shows that the phosphodiester bond cleavage step of the reaction is kinetically controlled and not thermodynamically because of a lack of any energy barrier above the accuracy of the energy profile calculation. The role of ions is shown by comparing the results with the reaction mechanisms in the absence of the Mg²⁺ ions where there is a significantly higher reaction barrier than in the presence of the Mg²⁺ ions.Graphical abstract

New Insight into Metal Ion-Driven Catalysis of Nucleic Acids by Influenza PA-Nter

PA subunit of influenza RNA-dependent RNA polymerase deserves constantly increasing attention due to its essential role in influenza life cycle. N-terminal domain of PA (PA-Nter) harbors endonuclease activity, which is indispensable in viral transcription and replication. Interestingly, existing literature reports on in vitro ion preferences of the enzyme are contradictory. Some show PA-Nter activity exclusively with Mn²⁺, whereas others report Mg²⁺ as a natural cofactor. To clarify it, we performed a series of experiments with varied ion concentrations and substrate type. We observed cleavage in the presence of both ions, with a slight preference for manganese, however PA-Nter activity highly depended on the amount of residual, co-purified ions. Furthermore, to quantify cleavage reaction rate, we applied fluorescence cross-correlation spectroscopy (FCCS), providing highly sensitive and real-time monitoring of single molecules. Using nanomolar ssDNA in the regime of enzyme excess, we estimated the maximum reaction rate at 0.81± 0.38 and 1.38± 0.34 nM/min for Mg²⁺ and Mn²⁺, respectively. However, our calculations of PA-Nter ion occupancy, based on thermodynamic data, suggest Mg²⁺ to be a canonical metal in PA-Nter processing of RNA in vivo. Presented studies constitute a step toward better understanding of PA-Nter ion-dependent activity, which will possibly contribute to new successful inhibitor design in the future.

Mechanism of Ribonuclease III Catalytic Regulation by Serine Phosphorylation

Ribonuclease III (RNase III) is a conserved, gene-regulatory bacterial endonuclease that cleaves double-helical structures in diverse coding and noncoding RNAs. RNase III is subject to multiple levels of control, reflective of its global regulatory functions. Escherichia coli (Ec) RNase III catalytic activity is known to increase during bacteriophage T7 infection, reflecting the expression of the phage-encoded protein kinase, T7PK. However, the mechanism of catalytic enhancement is unknown. This study shows that Ec-RNase III is phosphorylated on serine in vitro by purified T7PK, and identifies the targets as Ser33 and Ser34 in the N-terminal catalytic domain. Kinetic experiments reveal a 5-fold increase in k_cat and a 1.4-fold decrease in K_m following phosphorylation, providing a 7.4–fold increase in catalytic efficiency. Phosphorylation does not change the rate of substrate cleavage under single-turnover conditions, indicating that phosphorylation enhances product release, which also is the rate-limiting step in the steady-state. Molecular dynamics simulations provide a mechanism for facilitated product release, in which the Ser33 phosphomonoester forms a salt bridge with the Arg95 guanidinium group, thereby weakening RNase III engagement of product. The simulations also show why glutamic acid substitution at either serine does not confer enhancement, thus underscoring the specific requirement for a phosphomonoester.

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.

Fighting Cancer with Transition Metal Complexes: From Naked DNA to Protein and Chromatin Targeting Strategies

Many transition metal complexes have unique physicochemical properties that can be efficiently exploited in medicinal chemistry for cancer treatment. Traditionally, double-stranded DNA has been assumed to be the main binding target; however, recent studies have shown that nucleosomal DNA as well as proteins can act as dominant molecular binding partners. This has raised new questions about the molecular determinants that govern DNA versus protein binding selectivity, and has offered new ways to rationalize their biological activity and possible side effects. To address these questions, molecular simulations at an atomistic level of detail have been used to complement, support, and rationalize experimental data. Herein we review some relevant studies-focused on platinum and ruthenium compounds-to illustrate the power of state-of-the-art molecular simulation techniques and to demonstrate how the interplay between molecular simulations and experiments can make important contributions to elucidating the target preferences of some promising transition metal anticancer agents. This contribution aims at providing relevant information that may help in the rational design of novel drug-discovery strategies.

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.