Catalytic metal ions and enzymatic processing of DNA and RNA

Biochemical Characterization of Yeast Xrn1

Messenger RNA degradation is an important component of overall gene expression. During the final step of eukaryotic mRNA degradation, exoribonuclease 1 (Xrn1) carries out 5' → 3' processive, hydrolytic degradation of RNA molecules using divalent metal ion catalysis. To initiate studies of the 5' → 3' RNA decay machinery in our lab, we expressed a C-terminally truncated version of Saccharomyces cerevisiae Xrn1 and explored its enzymology using a second-generation, time-resolved fluorescence RNA degradation assay. Using this system, we quantitatively explored Xrn1's preference for 5'-monophosphorylated RNA substrates, its pH dependence, and the importance of active site mutations in the molecule's conserved catalytic core. Furthermore, we explore Xrn1's preference for RNAs containing a 5' single-stranded region both in an intermolecular hairpin structure and in an RNA-DNA hybrid duplex system. These results both expand and solidify our understanding of Xrn1, a centrally important enzyme whose biochemical properties have implications in numerous RNA degradation and processing pathways.

Design, Synthesis, Dynamic Docking, Biochemical Characterization, and in Vivo Pharmacokinetics Studies of Novel Topoisomerase II Poisons with Promising Antiproliferative Activity

We previously reported a first set of hybrid topoisomerase II (topoII) poisons whose chemical core merges key pharmacophoric elements of etoposide and merbarone, which are two well-known topoII blockers. Here, we report on the expansion of this hybrid molecular scaffold and present 16 more hybrid derivatives that have been designed, synthesized, and characterized for their ability to block topoII and for their overall drug-like profile. Some of these compounds act as topoII poison and exhibit good solubility, metabolic (microsomal) stability, and promising cytotoxicity in three cancer cell lines (DU145, HeLa, A549). Compound 3f (ARN24139) is the most promising drug-like candidate, with a good pharmacokinetics profile in vivo. Our results indicate that this hybrid new chemical class of topoII poisons deserves further exploration and that 3f is a favorable lead candidate as a topoII poison, meriting future studies to test its efficacy in in vivo tumor models.

Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5' Structure-Specific Nucleases

Enzymes of the 5′ structure-specific nuclease family are crucial for DNA repair, replication, and recombination. One such enzyme is the human exonuclease 1 (hExo1) metalloenzyme, which cleaves DNA strands, acting primarily as a processive 5′-3′ exonuclease and secondarily as a 5′-flap endonuclease. Recently, in crystallo reaction intermediates have elucidated how hExo1 exerts hydrolysis of DNA phosphodiester bonds. These hExo1 structures show a third metal ion intermittently bound close to the two-metal-ion active site, to which recessed ends or 5′-flap substrates bind. Evidence of this third ion has been observed in several nucleic-acid-processing metalloenzymes. However, there is still debate over what triggers the (un)binding of this transient third ion during catalysis and whether this ion has a catalytic function. Using extended molecular dynamics and enhanced sampling free-energy simulations, we observed that the carboxyl side chain of Glu89 (located along the arch motif in hExo1) flips frequently from the reactant state to the product state. The conformational flipping of Glu89 allows one metal ion to be recruited from the bulk and promptly positioned near the catalytic center. This is in line with the structural evidence. Additionally, our simulations show that the third metal ion assists the departure, through the mobile arch, of the nucleotide monophosphate product from the catalytic site. Structural comparisons of nuclease enzymes suggest that this Glu(Asp)-mediated mechanism for third ion recruitment and nucleic acid hydrolysis may be shared by other 5′ structure-specific nucleases.

DOTA-Branched Organic Frameworks as Giant and Potent Metal Chelators

Multinuclear complexes as metallo-agents for clinical use have caught extensive attention. In this paper, using 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) as both a functioning unit and a constructing junction, we build a series of DOTA-branched organic frameworks with multiple chelating holes by organizing DOTA layer by layer. These giant chelators are well characterized, which reveals their nanosized and soft structures. Further experiments demonstrate that they could efficiently hold abundant metal ions with much higher kinetic stabilities than the conventional small DOTA chelator. Their corresponding polynuclear complexes containing Gd³⁺, Tb³⁺, or both show superior imaging properties, excellent feasibility for peripheral modification, and unusual kinetic stability. This work can be easily extended to the fabrication of diverse homomultinuclear complexes and core/shell heteromultinuclear complexes with multifunctional properties. We expect that this new type of giant molecules and the ligand-branching strategy would open up a new avenue for the design and construction of next-generation polymetallic agents with high performance and stabilities for biomedical applications.

NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes

Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.

Study of the role of Mg2+ in dsRNA processing mechanism by bacterial RNase III through QM/MM simulations

The ribonuclease III (RNase III) cleaves dsRNA in specific positions generating mature RNAs. RNase III enzymes play important roles in RNA processing, post-transcriptional gene expression, and defense against viral infection. The enzyme's active site contains Mg²⁺ ions bound by a network of acidic residues and water molecules, but there is a lack of information about their specific roles. In this work, multiple steered molecular dynamics simulations at QM/MM level were performed to explore the hydrolysis reaction carried out by the enzyme. Free energy profiles modifying the features of the active site are obtained and the role of Mg²⁺ ions, the solvent molecules and the residues of the active site are discussed in detail. Our results show that Mg²⁺ ions carry out different roles in the hydrolysis process positioning the substrate for the attack from a coordinated nucleophile and activating it to perform hydrolysis reaction, cleaving the dsRNA backbone in a S_N2 substitution. In addition, water molecules present in the active site lower the energy barrier of the process. RNase III hydrolyzes dsRNA to generate mature RNAs. For this purpose, its active site contains Mg²⁺ which has an important role during the reaction. Results show that the Mg²⁺ activates the solvent molecule that produces the nucleophilic attack and the surrounding waters contribute significantly to the hydrolysis process.

Structural characterization of Class 2 OLD family nucleases supports a two-metal catalysis mechanism for cleavage

Overcoming lysogenization defect (OLD) proteins constitute a family of uncharacterized nucleases present in bacteria, archaea, and some viruses. These enzymes contain an N-terminal ATPase domain and a C-terminal Toprim domain common amongst replication, recombination, and repair proteins. The in vivo activities of OLD proteins remain poorly understood and no definitive structural information exists. Here we identify and define two classes of OLD proteins based on differences in gene neighborhood and amino acid sequence conservation and present the crystal structures of the catalytic C-terminal regions from the Burkholderia pseudomallei and Xanthamonas campestris p.v. campestris Class 2 OLD proteins at 2.24 Å and 1.86 Å resolution respectively. The structures reveal a two-domain architecture containing a Toprim domain with altered architecture and a unique helical domain. Conserved side chains contributed by both domains coordinate two bound magnesium ions in the active site of B. pseudomallei OLD in a geometry that supports a two-metal catalysis mechanism for cleavage. The spatial organization of these domains additionally suggests a novel mode of DNA binding that is distinct from other Toprim containing proteins. Together, these findings define the fundamental structural properties of the OLD family catalytic core and the underlying mechanism controlling nuclease activity.

Smoothed Potential MD Simulations for Dissociation Kinetics of Etoposide To Unravel Isoform Specificity in Targeting Human Topoisomerase II

Human type II topoisomerases (TopoII) are essential for controlling DNA topology within the cell. For this reason, there are a number of TopoII-targeted anticancer drugs that act by inducing DNA cleavage mediated by both TopoII isoforms (TopoIIα and TopoIIβ) in cells. However, recent studies suggest that specific poisoning of TopoIIα may be a safer strategy for treating cancer. This is because poisoning of TopoIIβ appears to be linked to the generation of secondary leukemia in patients. We recently reported that enzyme-mediated DNA cleavage complexes (in which TopoII is covalently linked to the cleaved DNA during catalysis) formed in the presence of the anticancer drug etoposide persisted approximately 3-fold longer with TopoIIα than TopoIIβ. Notably, enhanced drug-target residence time may reduce the adverse effects of specific TopoIIα poisons. However, it is still not clear how to design drugs that are specific for the α isoform. In this study, we report the results of classical molecular dynamics (MD) simulations to comparatively analyze the molecular interactions formed within the TopoII/DNA/etoposide complex with both isoforms. We also used smoothed potential MD to estimate etoposide dissociation kinetics from the two isoform complexes. These extensive classical and enhanced sampling simulations revealed stabilizing interactions of etoposide with two serine residues (Ser763 and Ser800) in TopoIIα. These interactions are missing in TopoIIβ, where both amino acids are alanine residues. This may explain the greater persistence of etoposide-stabilized cleavage complexes formed with Topo TopoIIα. These findings could be useful for the rational design of specific TopoIIα poisons.

Investigating mycobacterial topoisomerase I mechanism from the analysis of metal and DNA substrate interactions at the active site

We have obtained new crystal structures of Mycobacterium tuberculosis topoisomerase I, including structures with ssDNA substrate bound to the active site, with and without Mg²⁺ ion present. Significant enzyme conformational changes upon DNA binding place the catalytic tyrosine in a pre-transition state position for cleavage of a specific phosphodiester linkage. Meanwhile, the enzyme/DNA complex with bound Mg²⁺ ion may represent the post-transition state for religation in the enzyme's multiple-step DNA relaxation catalytic cycle. The first observation of Mg²⁺ ion coordinated with the TOPRIM residues and DNA phosphate in a type IA topoisomerase active site allows assignment of likely catalytic role for the metal and draws a comparison to the proposed mechanism for type IIA topoisomerases. The critical function of a strictly conserved glutamic acid in the DNA cleavage step was assessed through site-directed mutagenesis. The functions assigned to the observed Mg²⁺ ion can account for the metal requirement for DNA rejoining but not DNA cleavage by type IA topoisomerases. This work provides new structural insights into a more stringent requirement for DNA rejoining versus cleavage in the catalytic cycle of this essential enzyme, and further establishes the potential for selective interference of DNA rejoining by this validated TB drug target.

Factors Influencing the Activity of Nanozymes in the Cleavage of an RNA Model Substrate

A series of 2-nm gold nanoparticles passivated with different thiols all featuring at least one triazacyclonanone-Zn(II) complex and different flanking units (a second Zn(II) complex, a triethyleneoxymethyl derivative or a guanidinium of arginine of a peptide) were prepared and studied for their efficiency in the cleavage of the RNA-model substrate 2-hydroxypropyl-p-nitrophenyl phosphate. The source of catalysis for each of them was elucidated from the kinetic analysis (Michaelis–Menten profiles, pH dependence and kinetic isotope effect). The data indicated that two different mechanisms were operative: One involving two Zn(II) complexes and the other one involving a single Zn(II) complex and a flanking guanidinium cation. The mechanism based on a dinuclear catalytic site appeared more efficient than the one based on the cooperativity between a metal complex and a guanidinium.

Structural insights into the mechanism of double strand break formation by Hermes, a hAT family eukaryotic DNA transposase

Some DNA transposons relocate from one genomic location to another using a mechanism that involves generating double-strand breaks at their transposon ends by forming hairpins on flanking DNA. The same double-strand break mode is employed by the V(D)J recombinase at signal-end/coding-end junctions during the generation of antibody diversity. How flanking hairpins are formed during DNA transposition has remained elusive. Here, we describe several co-crystal structures of the Hermes transposase bound to DNA that mimics the reaction step immediately prior to hairpin formation. Our results reveal a large DNA conformational change between the initial cleavage step and subsequent hairpin formation that changes which strand is acted upon by a single active site. We observed that two factors affect the conformational change: the complement of divalent metal ions bound by the catalytically essential DDE residues, and the identity of the –2 flanking base pair. Our data also provides a mechanistic link between the efficiency of hairpin formation (an A:T basepair is favored at the –2 position) and Hermes' strong target site preference. Furthermore, we have established that the histidine residue within a conserved C/DxxH motif present in many transposase families interacts directly with the scissile phosphate, suggesting a crucial role in catalysis.

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.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Understanding the mechanistic basis of non-coding RNA through molecular dynamics simulations

Noncoding RNA (ncRNA) has a key role in regulating gene expression, mediating fundamental processes and diseases via a variety of yet unknown mechanisms. Here, we review recent applications of conventional and enhanced Molecular Dynamics (MD) simulations methods to address the mechanistic function of large biomolecular systems that are tightly involved in the ncRNA function and that are of key importance in life sciences. This compendium focuses of three biomolecular systems, namely the CRISPR-Cas9 genome editing machinery, group II intron ribozyme and the ribonucleoprotein complex of the spliceosome, which edit and process ncRNA. We show how the application of a novel accelerated MD simulations method has been key in disclosing the conformational transitions underlying RNA binding in the CRISPR-Cas9 complex, suggesting a mechanism for RNA recruitment and clarifying the conformational changes required for attaining genome editing. As well, we discuss the use of mixed quantum-classical MD simulations in deciphering the catalytic mechanism of RNA splicing as operated by group II intron ribozyme, one of the largest ncRNA structures crystallized so far. Finally, we debate the future challenges and opportunities in the field, discussing the recent application of MD simulations for unraveling the functional biophysics of the spliceosome, a multi-mega Dalton complex of proteins and small nuclear RNAs that performs RNA splicing in humans. This showcase of applications highlights the current talent of MD simulations to dissect atomic-level details of complex biomolecular systems instrumental for the design of finely engineered genome editing machines. As well, this review aims at inspiring future investigations of several other ncRNA regulatory systems, such as micro and small interfering RNAs, which achieve their function and specificity using RNA-based recognition and targeting strategies.

Deciphering Off-Target Effects in CRISPR-Cas9 through Accelerated Molecular Dynamics

CRISPR-Cas9 is the state-of-the-art technology for editing and manipulating nucleic acids. However, the occurrence of off-target mutations can limit its applicability. Here, all-atom enhanced molecular dynamics (MD) simulations-using Gaussian accelerated MD (GaMD)-are used to decipher the mechanism of off-target binding at the molecular level. GaMD reveals that base pair mismatches in the target DNA at distal sites with respect to the protospacer adjacent motif (PAM) can induce an extended opening of the RNA:DNA heteroduplex, which leads to newly formed interactions between the unwound DNA and the L2 loop of the catalytic HNH domain. These conserved interactions constitute a "lock" effectively decreasing the conformational freedom of the HNH domain and hampering its activation for cleavage. Remarkably, depending on their positions at PAM distal sites, DNA mismatches responsible for off-target cleavages are unable to "lock" the HNH domain, thereby leading to the unselective cleavage of DNA sequences. In consistency with the available experimental data, the ability to "lock" the catalytic HNH domain in an inactive "conformational checkpoint" is shown to be a key determinant in the onset of off-target effects. This mechanistic rationale contributes in clarifying a long lasting open issue in the CRISPR-Cas9 function and poses the foundation for designing novel and more specific Cas9 variants, which could be obtained by magnifying the "locking" interactions between HNH and the target DNA in the presence of any incorrect off-target sequence, thus preventing undesired cleavages.

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

The invisible dance of CRISPR-Cas9. Simulations unveil the molecular side of the gene-editing revolution

Since the discovery of the DNA double helix, the main molecular repository of genetic information, scientists have been struggling to find ways to efficiently manipulate genes. The ability to mark, modify, or regulate specific sequences of DNA in a controlled fashion is of key importance because of the ways that gene editing could be used to improve human life. For example, genetic therapies are being developed to permanently cure cancer and other life-threatening diseases.

Structure and Dynamics of the CRISPR-Cas9 Catalytic Complex

CRISPR-Cas9 is a bacterial immune system with exciting applications for genome editing. In spite of extensive experimental characterization, the active site chemistry of the RuvC domain-which performs DNA cleavages-has remained elusive. Its knowledge is key for structure-based engineering aimed at improving DNA cleavages. Here, we deliver an in-depth characterization by using quantum-classical (QM/MM) molecular dynamics (MD) simulations and a Gaussian accelerated MD method, coupled with bioinformatics analysis. We disclose a two-metal aided architecture in the RuvC active site, which is poised to operate DNA cleavages, in analogy with other DNA/RNA processing enzymes. The conformational dynamics of the RuvC domain further reveals that an "arginine finger" stably contacts the scissile phosphate, with the function of stabilizing the active complex. Remarkably, the formation of a catalytically competent state of the RuvC domain is only observed upon the conformational activation of the other nuclease domain of CRISPR-Cas9-i.e., the HNH domain-such allowing concerted cleavages of double stranded DNA. This structure is in agreement with the available experimental data and remarkably differs from previous models based on classical mechanics, demonstrating also that only quantum mechanical simulations can accurately describe the metal-aided active site in CRISPR-Cas9. This fully catalytic structure-in which both the HNH and RuvC domains are prone to perform DNA cleavages-constitutes a stepping-stone for understanding DNA cleavage and specificity. It calls for novel experimental verifications and offers the structural foundations for engineering efforts aimed at improving the genome editing capability of CRISPR-Cas9.

Two symmetric arginine residues play distinct roles in Thermus thermophilus Argonaute DNA guide strand-mediated DNA target cleavage

Bacterium Thermus thermophilus Argonaute (Ago; TtAgo) is a prokaryotic Ago (pAgo) that acts as the host defense against the uptake and propagation of foreign DNA by catalyzing the DNA cleavage reaction. The TtAgo active site consists of a plugged-in glutamate finger with two arginine residues (R545 and R486) located symmetrically around it. An interesting challenge is to understand how they can collaboratively facilitate enzymatic catalysis. In Kluyveromyces polysporus Ago, a eukaryotic Ago, the evolutionarily symmetrical residues are arginine and histidine, both of which function to stabilize the plugged-in catalytic tetrad conformation. Surprisingly, our simulation results indicated that, in TtAgo, only R545 is involved in the cleavage reaction by serving as a critical structural anchor to stabilize the catalytic tetrad Asp-Glu-Asp-Asp that is completed by the insertion of the glutamate finger, whereas R486 is not involved in target cleavage. The TtAgo-mediated target DNA cleavage occurs in a substrate-assisted mechanism, in which the pro-Rp (Rp, a tetrahedral phosphorus center with "R-type" chirality) oxygen of scissile phosphate acts as a general base to activate the nucleophilic water. Our unexpected theoretical findings on distinct roles played by R545 and R486 in TtAgo catalysis have been validated by single-point site-mutagenesis experiments, wherein the target cleavage is abolished for all mutants of R545. In sharp contrast, the cleavage activity is maintained for all mutants of R486. Our work provides mechanistic insights on the catalytic specificity of Ago proteins and could facilitate the design of new gene-editing tools in the long term.

Dynamic coordination of two-metal-ions orchestrates λ-exonuclease catalysis

Metal ions at the active site of an enzyme act as cofactors, and their dynamic fluctuations can potentially influence enzyme activity. Here, we use λ-exonuclease as a model enzyme with two Mg²⁺ binding sites and probe activity at various concentrations of magnesium by single-molecule-FRET. We find that while Mg_A²⁺ and Mg_B²⁺ have similar binding constants, the dissociation rate of Mg_A²⁺ is two order of magnitude lower than that of Mg_B²⁺ due to a kinetic-barrier-difference. At physiological Mg²⁺ concentration, the Mg_B²⁺ ion near the 5'-terminal side of the scissile phosphate dissociates each-round of degradation, facilitating a series of DNA cleavages via fast product-release concomitant with enzyme-translocation. At a low magnesium concentration, occasional dissociation and slow re-coordination of Mg_A²⁺ result in pauses during processive degradation. Our study highlights the importance of metal-ion-coordination dynamics in correlation with the enzymatic reaction-steps, and offers insights into the origin of dynamic heterogeneity in enzymatic catalysis.

Key role of the REC lobe during CRISPR-Cas9 activation by 'sensing', 'regulating', and 'locking' the catalytic HNH domain

Understanding the conformational dynamics of CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 is of the utmost importance for improving its genome editing capability. Here, molecular dynamics simulations performed using Anton-2 - a specialized supercomputer capturing micro-to-millisecond biophysical events in real time and at atomic-level resolution - reveal the activation process of the endonuclease Cas9 toward DNA cleavage. Over the unbiased simulation, we observe that the spontaneous approach of the catalytic domain HNH to the DNA cleavage site is accompanied by a remarkable structural remodeling of the recognition (REC) lobe, which exerts a key role for DNA cleavage. Specifically, the significant conformational changes and the collective conformational dynamics of the REC lobe indicate a mechanism by which the REC1-3 regions 'sense' nucleic acids, 'regulate' the HNH conformational transition, and ultimately 'lock' the HNH domain at the cleavage site, contributing to its catalytic competence. By integrating additional independent simulations and existing experimental data, we provide a solid validation of the activated HNH conformation, which had been so far poorly characterized, and we deliver a comprehensive understanding of the role of REC1-3 in the activation process. Considering the importance of the REC lobe in the specificity of Cas9, this study poses the basis for fully understanding how the REC components control the cleavage of off-target sequences, laying the foundation for future engineering efforts toward improved genome editing.

Phosphodiester models for cleavage of nucleic acids

Nucleic acids that store and transfer biological information are polymeric diesters of phosphoric acid. Cleavage of the phosphodiester linkages by protein enzymes, nucleases, is one of the underlying biological processes. The remarkable catalytic efficiency of nucleases, together with the ability of ribonucleic acids to serve sometimes as nucleases, has made the cleavage of phosphodiesters a subject of intensive mechanistic studies. In addition to studies of nucleases by pH-rate dependency, X-ray crystallography, amino acid/nucleotide substitution and computational approaches, experimental and theoretical studies with small molecular model compounds still play a role. With small molecules, the importance of various elementary processes, such as proton transfer and metal ion binding, for stabilization of transition states may be elucidated and systematic variation of the basicity of the entering or departing nucleophile enables determination of the position of the transition state on the reaction coordinate. Such data is important on analyzing enzyme mechanisms based on synergistic participation of several catalytic entities. Many nucleases are metalloenzymes and small molecular models offer an excellent tool to construct models for their catalytic centers. The present review tends to be an up to date summary of what has been achieved by mechanistic studies with small molecular phosphodiesters.

Crystal structures of Lymphocytic choriomeningitis virus endonuclease domain complexed with diketo-acid ligands

The Arenaviridae family, together with the Bunyaviridae and Orthomyxoviridae families, is one of the three negative-stranded RNA viral families that encode an endonuclease in their genome. The endonuclease domain is at the N-terminus of the L protein, a multifunctional protein that includes the RNA-dependent RNA polymerase. The synthesis of mRNA in arenaviruses is a process that is primed by capped nucleotides that are 'stolen' from the cellular mRNA by the endonuclease domain in cooperation with other domains of the L protein. This molecular mechanism has been demonstrated previously for the endonuclease of the prototype Lymphocytic choriomeningitis virus (LCMV). However, the mode of action of this enzyme is not fully understood as the original structure did not contain catalytic metal ions. The pivotal role played by the cap-snatching process in the life cycle of the virus and the highly conserved nature of the endonuclease domain make it a target of choice for the development of novel antiviral therapies. Here, the binding affinities of two diketo-acid (DKA) compounds (DPBA and L-742,001) for the endonuclease domain of LCMV were evaluated using biophysical methods. X-ray structures of the LCMV endonuclease domain with catalytic ions in complex with these two compounds were determined, and their efficacies were assessed in an in vitro endonuclease-activity assay. Based on these data and computational simulation, two new DKAs were synthesized. The LCMV endonuclease domain exhibits a good affinity for these DKAs, making them a good starting point for the design of arenavirus endonuclease inhibitors. In addition to providing the first example of an X-ray structure of an arenavirus endonuclease incorporating a ligand, this study provides a proof of concept that the design of optimized inhibitors against the arenavirus endonuclease is possible.

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.

Structural Origin of Metal Specificity in Isatin Hydrolase from Labrenzia aggregata Investigated by Computer Simulations

We performed quantum-chemical calculations, ab initio molecular dynamics, hybrid quantum mechanics/molecular mechanics (QM/MM) and enhanced sampling metadynamics simulations to investigate the origin of metal specificity in isatin hydrolase from Labrenzia aggregata. The peculiar octahedral binding geometry of the Mn²⁺ ion in the Michaelis complex includes both the isatin substrate and the catalytic water within the first coordination shell of the cation. Our calculations show that the same arrangement of the ligands cannot be efficiently achieved in the presence of other small divalent metal cations such as Zn²⁺ or Cu²⁺ . On the contrary, bulkier alkaline-earth cations such as Mg²⁺ , which allow octahedral coordination, are not able to activate the catalytic water into the stronger OH^- nucleophile required to attack the stable N-aryl-amide moiety of isatin.

Second-Shell Basic Residues Expand the Two-Metal-Ion Architecture of DNA and RNA Processing Enzymes

Synthesis and scission of phosphodiester bonds in DNA and RNA regulate vital processes within the cell. Enzymes that catalyze these reactions operate mostly via the recognized two-metal-ion mechanism. Our analysis reveals that basic amino acids and monovalent cations occupy structurally conserved positions nearby the active site of many two-metal-ion enzymes for which high-resolution (<3 Å) structures are known, including DNA and RNA polymerases, nucleases such as Cas9, and splicing ribozymes. Integrating multiple-sequence and structural alignments with molecular dynamics simulations, electrostatic potential maps, and mutational data, we found that these elements always interact with the substrates, suggesting that they may play an active role for catalysis, in addition to their electrostatic contribution. We discuss possible mechanistic implications of this expanded two-metal-ion architecture, including inferences on medium-resolution cryoelectron microscopy structures. Ultimately, our analysis may inspire future experiments and strategies for enzyme engineering or drug design to modulate nucleic acid processing.

•

Basic residues in the active site of two-metal-ion enzymes are structurally conserved

•

These residues are also conserved in evolution

•

Mutagenesis suggests these residues may exert an effect on DNA- and RNA processing

•

Our work offers insights into CRISPR/Cas9, spliceosome, and DNA/RNA polymerases

Structural insights into DNA cleavage activation of CRISPR-Cas9 system

CRISPR-Cas9 technology has been widely used for genome engineering. Its RNA-guided endonuclease Cas9 binds specifically to target DNA and then cleaves the two DNA strands with HNH and RuvC nuclease domains. However, structural information regarding the DNA cleavage-activating state of two nuclease domains remains sparse. Here, we report a 5.2 Å cryo-EM structure of Cas9 in complex with sgRNA and target DNA. This structure reveals a conformational state of Cas9 in which the HNH domain is closest to the DNA cleavage site. Compared with two known HNH states, our structure shows that the HNH active site moves toward the cleavage site by about 25 and 13 Å, respectively. In combination with EM-based molecular dynamics simulations, we show that residues of the nuclease domains in our structure could form cleavage-compatible conformations with the target DNA. Together, these results strongly suggest that our cryo-EM structure resembles a DNA cleavage-activating architecture of Cas9.

Novel xanthone-polyamine conjugates as catalytic inhibitors of human topoisomerase IIα

It has been proposed that xanthone derivatives with anticancer potential act as topoisomerase II inhibitors because they interfere with the ability of the enzyme to bind its ATP cofactor. In order to further characterize xanthone mechanism and generate compounds with potential as anticancer drugs, we synthesized a series of derivatives in which position 3 was substituted with different polyamine chains. As determined by DNA relaxation and decatenation assays, the resulting compounds are potent topoisomerase IIα inhibitors. Although xanthone derivatives inhibit topoisomerase IIα-catalyzed ATP hydrolysis, mechanistic studies indicate that they do not act at the ATPase site. Rather, they appear to function by blocking the ability of DNA to stimulate ATP hydrolysis. On the basis of activity, competition, and modeling studies, we propose that xanthones interact with the DNA cleavage/ligation active site of topoisomerase IIα and inhibit the catalytic activity of the enzyme by interfering with the DNA strand passage step.

Protospacer Adjacent Motif-Induced Allostery Activates CRISPR-Cas9

CRISPR-Cas9 is a genome editing technology with major impact in life sciences. In this system, the endonuclease Cas9 generates double strand breaks in DNA upon RNA-guided recognition of a complementary DNA sequence, which strictly requires the presence of a protospacer adjacent motif (PAM) next to the target site. Although PAM recognition is essential for cleavage, it is unknown whether and how PAM binding activates Cas9 for DNA cleavage at spatially distant sites. Here, we find evidence of a PAM-induced allosteric mechanism revealed by microsecond molecular dynamics simulations. PAM acts as an allosteric effector and triggers the interdependent conformational dynamics of the Cas9 catalytic domains (HNH and RuvC), responsible for concerted cleavage of the two DNA strands. Targeting such an allosteric mechanism should enable control of CRISPR-Cas9 functionality.

Activity inhibition and crystal polymorphism induced by active-site metal swapping

The Arenaviridae family is one of the two RNA viral families that encode a 3'-5' exonuclease in their genome. An exonuclease domain is found in the Arenaviridae nucleoprotein and targets dsRNA specifically. This domain is directly involved in suppression of innate immunity in the host cell. Like most phosphate-processing enzymes, it requires a divalent metal ion such as Mg²⁺ (or Mn²⁺) as a cofactor to catalyse nucleotide-cleavage and nucleotide-transfer reactions. On the other hand, calcium (Ca²⁺) inhibits this enzymatic activity, in spite of the fact that Mg²⁺ and Ca²⁺ present comparable binding affinities and biological availabilities. Here, the molecular and structural effects of the replacement of magnesium by calcium and its inhibition mechanism for phosphodiester cleavage, an essential reaction in the viral process of innate immunity suppression, are studied. Biochemical data and high-resolution structures of the Mopeia virus exonuclease domain complexed with each ion are reported for the first time. The consequences of the ion swap for the stability of the protein, the catalytic site and the functional role of a specific metal ion in enabling the catalytic cleavage of a dsRNA substrate are outlined.

CRISPR-Cas9 conformational activation as elucidated from enhanced molecular simulations

CRISPR-Cas9 has become a facile genome editing technology, yet the structural and mechanistic features underlying its function are unclear. Here, we perform extensive molecular simulations in an enhanced sampling regime, using a Gaussian-accelerated molecular dynamics (GaMD) methodology, which probes displacements over hundreds of microseconds to milliseconds, to reveal the conformational dynamics of the endonuclease Cas9 during its activation toward catalysis. We disclose the conformational transition of Cas9 from its apo form to the RNA-bound form, suggesting a mechanism for RNA recruitment in which the domain relocations cause the formation of a positively charged cavity for nucleic acid binding. GaMD also reveals the conformation of a catalytically competent Cas9, which is prone for catalysis and whose experimental characterization is still limited. We show that, upon DNA binding, the conformational dynamics of the HNH domain triggers the formation of the active state, explaining how the HNH domain exerts a conformational control domain over DNA cleavage [Sternberg SH et al. (2015) Nature, 527, 110-113]. These results provide atomic-level information on the molecular mechanism of CRISPR-Cas9 that will inspire future experimental investigations aimed at fully clarifying the biophysics of this unique genome editing machinery and at developing new tools for nucleic acid manipulation based on CRISPR-Cas9.

Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides

Ribonucleic acid (RNA) is involved in many regulatory and catalytic processes in the cell. The function of any RNA molecule is intimately related with its structure. In-line probing experiments provide valuable structural data sets for a variety of RNAs and are used to characterize conformational changes in riboswitches. However, the structural determinants that lead to differential reactivities in unpaired nucleotides have not been investigated yet. In this work, we used a combination of theoretical approaches, i.e., classical molecular dynamics simulations, multiscale quantum mechanical/molecular mechanical calculations, and enhanced sampling techniques in order to compute and interpret the differential reactivity of individual residues in several RNA motifs, including members of the most important GNRA and UNCG tetraloop families. Simulations on the multinanosecond timescale are required to converge the related free-energy landscapes. The results for uGAAAg and cUUCGg tetraloops and double helices are compared with available data from in-line probing experiments and show that the introduced technique is able to distinguish between nucleotides of the uGAAAg tetraloop based on their structural predispositions toward phosphodiester backbone cleavage. For the cUUCGg tetraloop, more advanced ab initio calculations would be required. This study is the first attempt to computationally classify chemical probing experiments and paves the way for an identification of tertiary structures based on the measured reactivity of nonreactive nucleotides.

Development of Site-Specific Mg(2+)-RNA Force Field Parameters: A Dream or Reality? Guidelines from Combined Molecular Dynamics and Quantum Mechanics Simulations

The vital contribution of Mg²⁺ ions to RNA biology is challenging to dissect at the experimental level. This calls for the integrative support of atomistic simulations, which at the classical level are plagued by limited accuracy. Indeed, force fields intrinsically neglect nontrivial electronic effects that Mg²⁺ exerts on its surrounding ligands in varying RNA coordination environments. Here, we present a combined computational study based on classical molecular dynamics (MD) and Density Functional Theory (DFT) calculations, aimed at characterizing (i) the performance of five Mg²⁺ force field (FF) models in RNA systems and (ii) how charge transfer and polarization affect the binding of Mg²⁺ ions in different coordination motifs. As a result, a total of ∼2.5 μs MD simulations (100/200 ns for each run) for two prototypical Mg²⁺-dependent ribozymes showed remarkable differences in terms of populations of inner-sphere coordination site types. Most importantly, complementary DFT calculations unveiled that differences in charge transfer and polarization among recurrent Mg²⁺-RNA coordination motifs are surprisingly small. In particular, the charge of the Mg²⁺ ions substantially remains constant through different coordination sites, suggesting that the common philosophy of developing site-specific Mg²⁺ ion parameters is not in line with the physical origin of the Mg²⁺-RNA MD simulations inaccuracies. Overall, this study constitutes a guideline for an adept use of current Mg²⁺ models and provides novel insights for the rational development of next-generation Mg²⁺ FFs to be employed for atomistic simulations of RNA.

Isotope effect analyses provide evidence for an altered transition state for RNA 2'-O-transphosphorylation catalyzed by Zn(2+)

Solvent D2O and (18)O kinetic isotope effects on RNA 2'-O-transphosphorylation catalyzed by Zn(2+) demonstrate an altered transition state relative to specific base catalysis. A recent model from DFT calculations involving inner sphere coordination to the non-bridging and leaving group oxygens is consistent with the data.