Catalytic metal ions and enzymatic processing of DNA and RNA

Correct Nucleotide Selection Is Confined at the Binding Site of Polymerase Enzymes

DNA polymerases (Pols) add incoming nucleotides (deoxyribonucleoside triphosphate (dNTPs)) to growing DNA strands, a crucial step for DNA synthesis. The insertion of correct (vs incorrect) nucleotides relates to Pols' fidelity, which defines Pols' ability to faithfully replicate DNA strands in a template-dependent manner. We and others have demonstrated that reactant alignment and correct base pairing at the Pols catalytic site are crucial structural features to fidelity. Here, we first used equilibrium molecular simulations to demonstrate that the local dynamics at the protein-DNA interface in the proximity of the catalytic site is different when correct vs incorrect dNTPs are bound to polymerase β (Pol β). Formation and dynamic stability of specific interatomic interactions around the incoming nucleotide influence the overall binding site architecture. This explains why certain Pols' mutants can affect the local catalytic environment and influence the selection of correct vs incorrect nucleotides. In particular, this is here demonstrated by analyzing the interaction network formed by the residue R283, whose mutant R283A has an experimentally measured lower capacity of differentiating correct (G:dCTP) vs incorrect (G:dATP) base pairing in Pol β. We also used alchemical free-energy calculations to quantify the G:dCTP →G:dATP transformation in Pol β wild-type and mutant R283A. These results correlate well with the experimental trend, thus corroborating our mechanistic insights. Sequence and structural comparisons with other Pols from the same family suggest that these findings may also be valid in similar enzymes.

Honeycomb-inspired ZIF-sealed interface enhances osseointegration via anti-infection and osteoimmunomodulation

Implant-associated infections (IAIs) pose a significant threat to orthopedic surgeries. Bacteria colonizing the surface of implants disrupt bone formation-related cells and interfere with the osteoimmune system, resulting in an impaired immune microenvironment and osteogenesis disorders. Inspired by nature, a zeolitic imidazolate framework (ZIF)-sealed smart drug delivery system on Ti substrates (ZSTG) was developed for the "natural-artificial dual-enzyme intervention (NADEI)" strategy to address these challenges. The subtle sealing design of ZIF-8 on the TiO2 nanotubes ensured glucose oxidase (GOx) activity and prevented its premature leakage. In the acidic infection microenvironment, the degradation of ZIF-8 triggered the rapid release of GOx, which converted glucose into H2O2 for disinfection. The Zn2+ released from degraded ZIF-8, as a DNase mimic, can hydrolyze extracellular DNA, which further enhances H2O2-induced disinfection and prevents biofilm formation. Importantly, Zn2+-mediated M2 macrophage polarization significantly improved the impaired osteoimmune microenvironment, accelerating bone repair. Transcriptomics revealed that ZSTG effectively suppressed the inflammatory cascade induced by lipopolysaccharide while promoting cell proliferation, homeostasis maintenance, and bone repair. In vitro and in vivo results confirmed the superior anti-infective, osteoimmunomodulatory, and osteointegrative capacities of the ZSTG-mediated NADEI strategy. Overall, this smart bionic platform has significant potential for future clinical applications to treat IAIs.

RNA targeting and cleavage by the type III-Dv CRISPR effector complex

CRISPR-Cas are adaptive immune systems in bacteria and archaea that utilize CRISPR RNA-guided surveillance complexes to target complementary RNA or DNA for destruction1-5. Target RNA cleavage at regular intervals is characteristic of type III effector complexes6-8. Here, we determine the structures of the Synechocystis type III-Dv complex, an apparent evolutionary intermediate from multi-protein to single-protein type III effectors9,10, in pre- and post-cleavage states. The structures show how multi-subunit fusion proteins in the effector are tethered together in an unusual arrangement to assemble into an active and programmable RNA endonuclease and how the effector utilizes a distinct mechanism for target RNA seeding from other type III effectors. Using structural, biochemical, and quantum/classical molecular dynamics simulation, we study the structure and dynamics of the three catalytic sites, where a 2'-OH of the ribose on the target RNA acts as a nucleophile for in line self-cleavage of the upstream scissile phosphate. Strikingly, the arrangement at the catalytic residues of most type III complexes resembles the active site of ribozymes, including the hammerhead, pistol, and Varkud satellite ribozymes. Our work provides detailed molecular insight into the mechanisms of RNA targeting and cleavage by an important intermediate in the evolution of type III effector complexes.

Flexibility of Binding Site is Essential to the Ca2+ Selectivity in EF-Hand Calcium-Binding Proteins

High binding affinity and selectivity of metal ions are essential to the function of metalloproteins. Thus, understanding the factors that determine these binding characteristics is of major interest for both fundamental mechanistic investigations and guiding of the design of novel metalloproteins. In this work, we perform QM cluster model calculations and quantum mechanics/molecular mechanics (QM/MM) free energy simulations to understand the binding selectivity of Ca2+ and Mg2+ in the wild-type carp parvalbumin and its mutant. While a nonpolarizable MM model (CHARMM36) does not lead to the correct experimental trend, treatment of the metal binding site with the DFTB3 model in a QM/MM framework leads to relative binding free energies (ΔΔGbind) comparable with experimental data. For the wild-type (WT) protein, the calculated ΔΔGbind is ∼6.6 kcal/mol in comparison with the experimental value of 5.6 kcal/mol. The good agreement highlights the value of a QM description of the metal binding site and supports the role of electronic polarization and charge transfer to metal binding selectivity. For the D51A/E101D/F102W mutant, different binding site models lead to considerable variations in computed binding affinities. With a coordination number of seven for Ca2+, which is shown by QM/MM metadynamics simulations to be the dominant coordination number for the mutant, the calculated relative binding affinity is ∼4.8 kcal/mol, in fair agreement with the experimental value of 1.6 kcal/mol. The WT protein is observed to feature a flexible binding site that accommodates a range of coordination numbers for Ca2+, which is essential to the high binding selectivity for Ca2+ over Mg2+. In the mutant, the E101D mutation reduces the flexibility of the binding site and limits the dominant coordination number of Ca2+ to be seven, thereby leading to reduced binding selectivity against Mg2+. Our results highlight that the binding selectivity of metal ions depends on both the structural and dynamical properties of the protein binding site.

Elucidation of the catalytic mechanism of a single-metal dependent homing endonuclease using QM and QM/MM approaches: the case study of I-PpoI

Homing endonucleases (HEs) are highly specific DNA cleaving enzymes, with I-PpoI having been suggested to use a single metal to accelerate phosphodiester bond cleavage. Although an I-PpoI mechanism has been proposed based on experimental structural data, no consensus has been reached regarding the roles of the metal or key active site amino acids. This study uses QM cluster and QM/MM calculations to provide atomic-level details of the I-PpoI catalytic mechanism. Minimal QM cluster and large-scale QM/MM models demonstrate that the experimentally-proposed pathway involving direct Mg2+ coordination to the substrate coupled with leaving group protonation through a metal-activated water is not feasible due to an inconducive I-PpoI active site alignment. Despite QM cluster models of varying size uncovering a pathway involving leaving group protonation by a metal-activated water, indirect (water-mediated) metal coordination to the substrate is required to afford this pathway, which renders this mechanism energetically infeasible. Instead, QM cluster models reveal that the preferred pathway involves direct Mg2+-O3' coordination to stabilize the charged substrate and assist leaving group departure, while H98 activates the water nucleophile. These calculations also underscore that both catalytic residues that directly interact with the substrate and secondary amino acids that position or stabilize these residues are required for efficient catalysis. QM/MM calculations on the solvated enzyme-DNA complex verify the preferred mechanism, which is fully consistent with experimental kinetic, structural, and mutational data. The fundamental understanding of the I-PpoI mechanism of action, gained from the present work can be used to further explore potential uses of this enzyme in biotechnology and medicine, and direct future computational investigations of other members of the understudied HE family.

NaOAc-Assisted Aerobic Oxidation Protocol for the Synthesis of Pentacoordinate Chalcogenyl Spirophosphoranes with P-Se/P-S Bonds under Open Air

The NaOAc-assisted aerobic oxidation reaction of pentacoordinate hydrospirophosphoranes and dichalcogenyl compounds with open air as a green oxidant has been developed under mild conditions. A series of novel pentacoordinate spirophosphoranes with P-Se/P-S bonds were synthesized in excellent yields. The reaction mechanism was determined by 31P nuclear magnetic resonance tracing experiments, high-resolution mass spectrometry tracing experiments, and X-ray diffraction analysis. The method features a broad substrate scope, good functional group tolerance, and a high degree of atomic utilization and is meaningful for the synthesis of bioactive chalcogenphosphate compounds with chalcogen and phosphorus moieties.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

Everything OLD is new again: How structural, functional, and bioinformatic advances have redefined a neglected nuclease family

Overcoming lysogenization defect (OLD) proteins are a conserved family of ATP-powered nucleases that function in anti-phage defense. Recent bioinformatic, genetic, and crystallographic studies have yielded new insights into the structure, function, and evolution of these enzymes. Here we review these developments and propose a new classification scheme to categorize OLD homologs that relies on gene neighborhoods, biochemical properties, domain organization, and catalytic machinery. This taxonomy reveals important similarities and differences between family members and provides a blueprint to contextualize future in vivo and in vitro findings. We also detail how OLD nucleases are related to PARIS and Septu anti-phage defense systems and discuss important mechanistic questions that remain unanswered.

An unnatural enzyme with endonuclease activity towards small non-coding RNAs

Endonucleases are enzymes that cleave internal phosphodiester bonds within double-stranded DNA or RNA and are essential for biological functions. Herein, we use genetic code expansion to create an unnatural endonuclease that cleaves non-coding RNAs including short interfering RNA (siRNA) and microRNAs (miRNAs), a function that does not exist in nature. We introduce a metal-chelating unnatural amino acid, (2,2'-bipyridin-5-yl)alanine (BpyAla) to impart endonuclease activity to the viral suppressor of RNA silencing protein p19. Upon binding of copper, the mutant p19-T111BpyAla displays catalytic site-specific cleavage of siRNA and human miRNAs. Catalysis is confirmed using fluorescence polarization and fluorescence turn-on. Global miRNA profiling reveals that the engineered enzyme cleaves miRNAs in a human cell line. The therapeutic potential is demonstrated by targeting miR-122, a critical host factor for the hepatitis C virus (HCV). Unnatural endonuclease function is shown to deplete miR-122 levels with similar effects to an antagomir that reduces HCV levels therapeutically.

Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides

Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

Molecular Complementarity of Proteomimetic Materials for Target-Specific Recognition and Recognition-Mediated Complex Functions

As biomolecules essential for sustaining life, proteins are generated from long chains of 20 different α-amino acids that are folded into unique 3D structures. In particular, many proteins have molecular recognition functions owing to their binding pockets, which have complementary shapes, charges, and polarities for specific targets, making these biopolymers unique and highly valuable for biomedical and biocatalytic applications. Based on the understanding of protein structures and microenvironments, molecular complementarity can be exhibited by synthesizable and modifiable materials. This has prompted researchers to explore the proteomimetic potentials of a diverse range of materials, including biologically available peptides and oligonucleotides, synthetic supramolecules, inorganic molecules, and related coordination networks. To fully resemble a protein, proteomimetic materials perform the molecular recognition to mediate complex molecular functions, such as allosteric regulation, signal transduction, enzymatic reactions, and stimuli-responsive motions; this can also expand the landscape of their potential bio-applications. This review focuses on the recognitive aspects of proteomimetic designs derived for individual materials and their conformations. Recent progress provides insights to help guide the development of advanced protein mimicry with material heterogeneity, design modularity, and tailored functionality. The perspectives and challenges of current proteomimetic designs and tools are also discussed in relation to future applications.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

On the Role of Molecular Conformation of the 8-Oxoguanine Lesion in Damaged DNA Processing by Polymerases

A common and insidious DNA damage is 8-oxoguanine (8OG), bypassed with low catalytic efficiency and high error frequency by polymerases (Pols) during DNA replication. This is a fundamental process with far-reaching implications in cell function and diseases. However, the molecular determinants of how 8OG exactly affects the catalytic efficiency of Pols remain largely unclear. By examining ternary deoxycytidine triphosphate/DNA/Pol complexes containing the 8OG damage, we found that 8OG consistently adopts different conformations when bound to Pols, compared to when in isolated DNA. Equilibrium molecular dynamics and metadynamics free energy calculations quantified that 8OG is in the lowest energy conformation in isolated DNA. In contrast, 8OG adopts high-energy conformations often characterized by intramolecular steric repulsion when bound to Pols. We show that the 8OG conformation can be regulated by mutating Pol residues interacting with the 8OG phosphate group. These findings propose the 8OG conformation as a factor in Pol-mediated processing of damaged DNA.

The Electronic Structure of Genome Editors from the First Principles

Genome editing based on the CRISPR-Cas9 system has paved new avenues for medicine, pharmaceutics, biotechnology, and beyond. This article reports the role of first-principles (ab-initio) molecular dynamics (MD) in the CRISPR-Cas9 revolution, achieving a profound understanding of the enzymatic function and offering valuable insights for enzyme engineering. We introduce the methodologies and explain the use of ab-initio MD simulations to characterize the two-metal dependent mechanism of DNA cleavage in the RuvC domain of the Cas9 enzyme, and how a second catalytic domain, HNH, cleaves the target DNA with the aid of a single metal ion. A detailed description of how ab-initio MD is combined with free-energy methods - i.e., thermodynamic integration and metadynamics - to break and form chemical bonds is given, explaining the use of these methods to determine the chemical landscape and establish the catalytic mechanism in CRISPR-Cas9. The critical role of classical methods is also discussed, explaining theory and application of constant pH MD simulations, used to accurately predict the catalytic residues' protonation states. Overall, first-principles methods are shown to unravel the electronic structure of the Cas9 enzyme, providing valuable insights that can serve for the design of genome editing tools with improved catalytic efficiency or controllable activity.

Distal Mutations in the β-Clamp of DNA Polymerase III* Disrupt DNA Orientation and Affect Exonuclease Activity

DNA polymerases are responsible for the replication and repair of DNA found in all DNA-based organisms. DNA Polymerase III is the main replicative polymerase of E. coli and is composed of over 10 proteins. A subset of these proteins (Pol III*) includes the polymerase (α), exonuclease (ϵ), clamp (β), and accessory protein (θ). Mutations of residues in, or around the active site of the catalytic subunits (α and ϵ), can have a significant impact on catalysis. However, the effects of distal mutations in noncatalytic subunits on the activity of catalytic subunits are less well-characterized. Here, we investigate the effects of two Pol III* variants, β-L82E/L82'E and β-L82D/L82'D, on the proofreading reaction catalyzed by ϵ. MD simulations reveal major changes in the dynamics of Pol III*, which extend throughout the complex. These changes are mostly induced by a shift in the position of the DNA substrate inside the β-clamp, although no major structural changes are observed in the protein complex. Quantum mechanics/molecular mechanics (QM/MM) calculations indicate that the β-L82D/L82'D variant has reduced catalytic proficiency due to highly endoergic reaction energies resulting from structural changes in the active site and differences in the electric field at the active site arising from the protein and substrate. Conversely, the β-L82E/L82'E variant is predicted to maintain proofreading activity, exhibiting a similar reaction barrier for nucleotide excision compared with the WT system. However, significant differences in the reaction mechanism are obtained due to the changes induced by the mutations on the β-clamp.

A Practical Approach to Bicyclic Carbamoyl Pyridones with Application to the Synthesis of HIV-1 Integrase Strand Transfer Inhibitors

An efficient one-pot synthetic method has been developed for the preparation of bicyclic carbamoyl pyridones from the known common intermediate methyl 5-((2,4-difluorobenzyl)carbamoyl)-1-(2,2-dimethoxyethyl)-3-methoxy-4-oxo-1,4-dihydropyridine-2-carboxylate (8). The scalable protocol is facile and employs readily available reagents, needing only a single purification as the final step. The utility of the approach was demonstrated by preparing a library of HIV-1 integrase strand transfer inhibitors (INSTIs) that differ by the presence or absence of a double bond in the B-ring of the bicyclic carbamoyl pyridines 6 and 7. Several of the analogs show good antiviral potencies in single-round HIV-1 replication antiviral assays and show no cytotoxicity in cell culture assays. In general, the compounds with a B-ring double bond have higher antiviral potencies than their saturated congeners. Our methodology should be applicable to the synthesis of a range of new metal-chelating analogs.

Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions

Coulomb interactions play a crucial role in a wide array of processes in aqueous solutions but present conceptual and computational challenges to both theory and simulations. We review recent developments in an approach addressing these challenges─local molecular field (LMF) theory. LMF theory exploits an exact and physically suggestive separation of intermolecular Coulomb interactions into strong short-range and uniformly slowly varying long-range components. This allows us to accurately determine the averaged effects of the long-range components on the short-range structure using effective single particle fields and analytical corrections, greatly reducing the need for complex lattice summation techniques used in most standard approaches. The simplest use of these ideas in aqueous solutions leads to the short solvent (SS) model, where both solvent-solvent and solute-solvent Coulomb interactions have only short-range components. Here we use the SS model to give a simple description of pairing of nucleobases and biologically relevant ions in water.

Trypanosoma brucei RRP44: a versatile enzyme for processing structured and non-structured RNA substrates

Rrp44/Dis3 is a conserved eukaryotic ribonuclease that acts on processing and degradation of nearly all types of RNA. It contains an endo- (PIN) and an exonucleolytic (RNB) domain and, its depletion in model organisms supports its essential function for cell viability. In Trypanosoma brucei, depletion of Rrp44 (TbRRP44) blocks maturation of ribosomal RNA, leading to disruption of ribosome synthesis and inhibition of cell proliferation. We have determined the crystal structure of the exoribonucleolytic module of TbRRP44 in an active conformation, revealing novel details of the catalytic mechanism of the RNB domain. For the first time, the position of the second magnesium involved in the two-metal-ion mechanism was determined for a member of the RNase II family. In vitro, TbRRP44 acts preferentially on non-structured uridine-rich RNA substrates. However, we demonstrated for the first time that both TbRRP44 and its homologue from Saccharomyces cerevisiae can also degrade structured substrates without 3'-end overhang, suggesting that Rrp44/Dis3 ribonucleases may be involved in degradation of a wider panel of RNA than has been assumed. Interestingly, deletion of TbRRP44 PIN domain impairs RNA binding to different extents, depending on the type of substrate.

Inspection on the Mechanism of SARS-CoV-2 Inhibition by Penciclovir: A Molecular Dynamic Study

In recent years, humanity has had to face a critical pandemic due to SARS-CoV-2. In the rapid search for effective drugs against this RNA-positive virus, the repurposing of already existing nucleotide/nucleoside analogs able to stop RNA replication by inhibiting the RNA-dependent RNA polymerase enzyme has been evaluated. In this process, a valid contribution has been the use of in silico experiments, which allow for a rapid evaluation of the possible effectiveness of the proposed drugs. Here we propose a molecular dynamic study to provide insight into the inhibition mechanism of Penciclovir, a nucleotide analog on the RNA-dependent RNA polymerase enzyme. Besides the presented results, in this article, for the first time, molecular dynamic simulations have been performed considering not only the RNA-dependent RNA polymerase protein, but also its cofactors (fundamental for RNA replication) and double-strand RNA.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

On the Recognition of Natural Substrate CTP and Endogenous Inhibitor ddhCTP of SARS-CoV-2 RNA-Dependent RNA Polymerase: A Molecular Dynamics Study

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 outbreak that is affecting the entire planet. As the pandemic is still spreading worldwide, with multiple mutations of the virus, it is of interest and of help to employ computational methods for identifying potential inhibitors of the enzymes responsible for viral replication. Attractive antiviral nucleotide analogue RNA-dependent RNA polymerase (RdRp) chain terminator inhibitors are investigated with this purpose. This study, based on molecular dynamics (MD) simulations, addresses the important aspects of the incorporation of an endogenously synthesized nucleoside triphosphate, ddhCTP, in comparison with the natural nucleobase cytidine triphosphate (CTP) in RdRp. The ddhCTP species is the product of the viperin antiviral protein as part of the innate immune response. The absence of the ribose 3'-OH in ddhCTP could have important implications in its inhibitory mechanism of RdRp. We built an in silico model of the RNA strand embedded in RdRp using experimental methods, starting from the cryo-electron microscopy structure and exploiting the information obtained by spectrometry on the RNA sequence. We determined that the model was stable during the MD simulation time. The obtained results provide deeper insights into the incorporation of nucleoside triphosphates, whose molecular mechanism by the RdRp active site still remains elusive.

Principles of target DNA cleavage and the role of Mg2+ in the catalysis of CRISPR-Cas9

At the core of the CRISPR-Cas9 genome-editing technology, the endonuclease Cas9 introduces site-specific breaks in DNA. However, precise mechanistic information to ameliorating Cas9 function is still missing. Here, multi-microsecond molecular dynamics, free-energy and multiscale simulations are combined with solution NMR and DNA cleavage experiments to resolve the catalytic mechanism of target DNA cleavage. We show that the conformation of an active HNH nuclease is tightly dependent on the catalytic Mg²⁺, unveiling its cardinal structural role. This activated Mg²⁺-bound HNH is consistently described through molecular simulations, solution NMR and DNA cleavage assays, revealing also that the protonation state of the catalytic H840 is strongly affected by active site mutations. Finally, ab-initio QM(DFT)/MM simulations and metadynamics establish the catalytic mechanism, showing that the catalysis is activated by H840 and completed by K866, rationalising DNA cleavage experiments. This information is critical to enhance the enzymatic function of CRISPR-Cas9 toward improved genome-editing.

Implementing accelerated dynamics to unravel the effects of high-fidelity Cas9 mutants on target DNA and guide RNA hybrid stability

The clustered regularly interspersed short palindromic repeats (CRISPR) and its associated nuclease (Cas9) offers a unique and easily reprogrammable system for editing eukaryotic genomes. Cas9 is guided to the target by an RNA strand, and precise edits are created by introducing double-stranded breaks. However, nuclease activity of Cas9 is also triggered at other sites other than the target sit, which is a major limitation for various applications. Cas9 variants have been designed to improve the efficacy of the tool by introducing certain mutations. However, the on-target activity of such Cas9 variants is often seen as compromised. Hence, understanding the sub-molecular differences in the variants is essential to elucidate the factors that contribute to efficiency. The study reveals distortions in the PAM-distal regions of the nucleic hybrids as well as changes in the interactions between the Cas9 variants and RNA-DNA hybrid, contributing to the explanation for differences in on-target activity.Communicated by Ramaswamy H. Sarma.

Structural and Molecular Kinetic Features of Activities of DNA Polymerases

DNA polymerases catalyze DNA synthesis during the replication, repair, and recombination of DNA. Based on phylogenetic analysis and primary protein sequences, DNA polymerases have been categorized into seven families: A, B, C, D, X, Y, and RT. This review presents generalized data on the catalytic mechanism of action of DNA polymerases. The structural features of different DNA polymerase families are described in detail. The discussion highlights the kinetics and conformational dynamics of DNA polymerases from all known polymerase families during DNA synthesis.

Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations

Considering that practically all reactions that involve nucleotides also involve metal ions, it is evident that the coordination chemistry of nucleotides and their derivatives is an essential corner stone of biological inorganic chemistry. Nucleotides are either directly or indirectly involved in all processes occurring in Nature. It is therefore no surprise that the constituents of nucleotides have been chemically altered-that is, at the nucleobase residue, the sugar moiety, and also at the phosphate group, often with the aim of discovering medically useful compounds. Among such derivatives are acyclic nucleoside phosphonates (ANPs), where the sugar moiety has been replaced by an aliphatic chain (often also containing an ether oxygen atom) and the phosphate group has been replaced by a phosphonate carrying a carbon-phosphorus bond to make the compounds less hydrolysis-sensitive. Several of these ANPs show antiviral activity, and some of them are nowadays used as drugs. The antiviral activity results from the incorporation of the ANPs into the growing nucleic acid chain-i.e., polymerases accept the ANPs as substrates, leading to chain termination because of the missing 3'-hydroxyl group. We have tried in this review to describe the coordination chemistry (mainly) of the adenine nucleotides AMP and ATP and whenever possible to compare it with that of the dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2- = adenine(N9)-CH2-CH2-O-CH2-PO32) [or its diphosphate (PMEApp4-)] as a representative of the ANPs. Why is PMEApp4- a better substrate for polymerases than ATP4-? There are three reasons: (i) PMEA2- with its anti-like conformation (like AMP2-) fits well into the active site of the enzyme. (ii) The phosphonate group has an enhanced metal ion affinity because of its increased basicity. (iii) The ether oxygen forms a 5-membered chelate with the neighboring phosphonate and favors thus coordination at the Pα group. Research on ANPs containing a purine residue revealed that the kind and position of the substituent at C2 or C6 has a significant influence on the biological activity. For example, the shift of the (C6)NH2 group in PMEA to the C2 position leads to 9-[2-(phosphonomethoxy)ethyl]-2-aminopurine (PME2AP), an isomer with only a moderate antiviral activity. Removal of (C6)NH2 favors N7 coordination, e.g., of Cu2+, whereas the ether O atom binding of Cu2+ in PMEA facilitates N3 coordination via adjacent 5- and 7-membered chelates, giving rise to a Cu(PMEA)cl/O/N3 isomer. If the metal ions (M2+) are M(α,β)-M(γ)-coordinated at a triphosphate chain, transphosphorylation occurs (kinases, etc.), whereas metal ion binding in a M(α)-M(β,γ)-type fashion is relevant for polymerases. It may be noted that with diphosphorylated PMEA, (PMEApp4-), the M(α)-M(β,γ) binding is favored because of the formation of the 5-membered chelate involving the ether O atom (see above). The self-association tendency of purines leads to the formation of dimeric [M2(ATP)]2(OH)- stacks, which occur in low concentration and where one half of the molecule undergoes the dephosphorylation reaction and the other half stabilizes the structure-i.e., acts as the "enzyme" by bridging the two ATPs. In accord herewith, one may enhance the reaction rate by adding AMP2- to the [Cu2(ATP)]2(OH)- solution, as this leads to the formation of mixed stacked Cu3(ATP)(AMP)(OH)- species, in which AMP2- takes over the structuring role, while the other "half" of the molecule undergoes dephosphorylation. It may be added that Cu3(ATP)(PMEA) or better Cu3(ATP)(PMEA)(OH)- is even a more reactive species than Cu3(ATP)(AMP)(OH)-. - The matrix-assisted self-association and its significance for cell organelles with high ATP concentrations is summarized and discussed, as is, e.g., the effect of tryptophanate (Trp-), which leads to the formation of intramolecular stacks in M(ATP)(Trp)3- complexes (formation degree about 75%). Furthermore, it is well-known that in the active-site cavities of enzymes the dielectric constant, compared with bulk water, is reduced; therefore, we have summarized and discussed the effect of a change in solvent polarity on the stability and structure of binary and ternary complexes: Opposite effects on charged O sites and neutral N sites are observed, and this leads to interesting insights.

Two-Metal-Ion Catalysis: Inhibition of DNA Polymerase Activity by a Third Divalent Metal Ion

Almost all DNA polymerases (pols) exhibit bell-shaped activity curves as a function of both pH and Mg²⁺ concentration. The pol activity is reduced when the pH deviates from the optimal value. When the pH is too low the concentration of a deprotonated general base (namely, the attacking 3′-hydroxyl of the 3′ terminal residue of the primer strand) is reduced exponentially. When the pH is too high the concentration of a protonated general acid (i.e., the leaving pyrophosphate group) is reduced. Similarly, the pol activity also decreases when the concentration of the divalent metal ions deviates from its optimal value: when it is too low, the binding of the two catalytic divalent metal ions required for the full activity is incomplete, and when it is too high a third divalent metal ion binds to pyrophosphate, keeping it in the replication complex longer and serving as a substrate for pyrophosphorylysis within the complex. Currently, there is a controversy about the role of the third metal ion which we will address in this review.

Controlled Trafficking of Multiple and Diverse Cations Prompts Nucleic Acid Hydrolysis

Recent in crystallo reaction intermediates have detailed how nucleic acid hydrolysis occurs in the RNA ribonuclease H1 (RNase H1), a fundamental metalloenzyme involved in maintaining the human genome. At odds with the previous characterization, these in crystallo structures unexpectedly captured multiple metal ions (K⁺ and Mg²⁺) transiently bound in the vicinity of the two-metal-ion active site. Using multi-microsecond all-atom molecular dynamics and free-energy simulations, we investigated the functional implications of the dynamic exchange of multiple K⁺ and Mg²⁺ ions at the RNase H1 reaction center. We found that such ions are timely positioned at non-overlapping locations near the active site, at different stages of catalysis, being crucial for both reactants' alignment and leaving group departure. We also found that this cation trafficking is tightly regulated by variations of the solution's ionic strength and is aided by two conserved second-shell residues, E188 and K196, suggesting a mechanism for the cations' recruitment during catalysis. These results indicate that controlled trafficking of multi-cation dynamics, opportunely prompted by second-shell residues, is functionally essential to the complex enzymatic machinery of the RNase H1. These findings revise the current knowledge on the RNase H1 catalysis and open new catalytic possibilities for other similar metalloenzymes including, but not limited to, CRISPR-Cas9, group II intron ribozyme and the human spliceosome.

On the Metal-Aided Catalytic Mechanism for Phosphodiester Bond Cleavage Performed by Nanozymes

Recent studies have shown that gold nanoparticles (AuNPs) functionalized with Zn(II) complexes can cleave phosphate esters and nucleic acids. Remarkably, such synthetic nanonucleases appear to catalyze metal (Zn)-aided hydrolytic reactions of nucleic acids similar to metallonuclease enzymes. To clarify the reaction mechanism of these nanocatalysts, here we have comparatively analyzed two nanonucleases with a >10-fold difference in the catalytic efficiency for the hydrolysis of the 2-hydroxypropyl-4-nitrophenylphosphate (HPNP, a typical RNA model substrate). We have used microsecond-long atomistic simulations, integrated with NMR experiments, to investigate the structure and dynamics of the outer coating monolayer of these nanoparticles, either alone or in complex with HPNP, in solution. We show that the most efficient one is characterized by coating ligands that promote a well-organized monolayer structure, with the formation of solvated bimetallic catalytic sites. Importantly, we have found that these nanoparticles can mimic two-metal-ion enzymes for nucleic acid processing, with Zn ions that promote HPNP binding at the reaction center. Thus, the two-metal-ion-aided hydrolytic strategy of such nanonucleases helps in explaining their catalytic efficiency for substrate hydrolysis, in accordance with the experimental evidence. These mechanistic insights reinforce the parallelism between such functionalized AuNPs and proteins toward the rational design of more efficient catalysts.

Finding the Ion in the RNA-Stack: Can Computational Models Accurately Predict Key Functional Elements in Large Macromolecular Complexes?

This viewpoint discusses the predictive power and impact of computational analyses and simulations to gain prospective, experimentally supported mechanistic insights into complex biological systems. Remarkably, two newly resolved cryoEM structures have confirmed the previous, and independent, prediction of the precise localization and dynamics of key catalytic ions in megadalton-large spliceosomal complexes. This outstanding outcome endorses a prominent synergy of computational and experimental methods in the prospective exploration of such large multicomponent biosystems.

The Mechanism of Cleavage of RNA Phosphodiesters by a Gold Nanoparticle Nanozyme

The cleavage of uridine 3'-phosphodiesters bearing alcohols with pKa ranging from 7.14 to 14.5 catalyzed by AuNPs functionalized with 1,4,7-triazacyclononane-Zn(II) complexes has been studied to unravel the source of catalysis by these nanosystems (nanozymes). The results have been compared with those obtained with two Zn(II) dinuclear catalysts for which the mechanism is fairly understood. Binding to the Zn(II) ions by the substrate and the uracil of uridine was observed. The latter leads to inhibition of the process and formation of less productive binding complexes than in the absence of the nucleobase. The nanozyme operates with these substrates mostly via a nucleophilic mechanism with little stabilization of the pentacoordinated phosphorane and moderate assistance in leaving group departure. This is attributed to a decrease of binding strength of the substrate to the catalytic site in reaching the transition state due to an unfavorable binding mode with the uracil. The nanozyme favors substrates with better leaving groups than the less acidic ones.

Establishing the allosteric mechanism in CRISPR-Cas9

Allostery is a fundamental property of proteins, which regulates biochemical information transfer between spatially distant sites. Here, we report on the critical role of molecular dynamics (MD) simulations in discovering the mechanism of allosteric communication within CRISPR-Cas9, a leading genome editing machinery with enormous promises for medicine and biotechnology. MD revealed how allostery intervenes during at least three steps of the CRISPR-Cas9 function: affecting DNA recognition, mediating the cleavage and interfering with the off-target activity. An allosteric communication that activates concerted DNA cleavages was found to led through the L1/L2 loops, which connect the HNH and RuvC catalytic domains. The identification of these "allosteric transducers" inspired the development of novel variants of the Cas9 protein with improved specificity, opening a new avenue for controlling the CRISPR-Cas9 activity. Discussed studies also highlight the critical role of the recognition lobe in the conformational activation of the catalytic HNH domain. Specifically, the REC3 region was found to modulate the dynamics of HNH by sensing the formation of the RNA:DNA hybrid. The role of REC3 was revealed to be particularly relevant in the presence of DNA mismatches. Indeed, interference of REC3 with the RNA:DNA hybrid containing mismatched pairs at specific positions resulted in locking HNH in an inactive "conformational checkpoint" conformation, thereby hampering off-target cleavages. Overall, MD simulations established the fundamental mechanisms underlying the allosterism of CRISPR-Cas9, aiding engineering strategies to develop new CRISPR-Cas9 variants for improved genome editing.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

Active-Site Models of Streptococcus pyogenes Cas9 in DNA Cleavage State

CRISPR-Cas9 is a powerful tool for target genome editing in living cells. Significant advances have been made to understand how this system cleaves target DNA. HNH is a nuclease domain, which shares structural similarity with the HNH endonuclease characterzied by a beta-beta-alpha-metal fold. Therefore, based on one- and two-metal-ion mechanisms, homology modeling and molecular dynamics (MD) simulation are suitable tools for building an atomic model of Cas9 in the DNA cleavage state. Here, by modeling and MD, we presented an atomic model of SpCas9-sgRNA-DNA complex with the cleavage state. This model shows that the HNH and RuvC conformations resemble their DNA cleavage state where the active-sites in the complex coordinate with DNA, Mg²⁺ ions, and water. Among them, residues D10, E762, H983, and D986 locate at the first shell of the RuvC active-site and interact with the ions directly, residues H982 or/and H985 are general (Lewis) bases, and the coordinated water is located at the positions for nucleophilic attack of the scissile phosphate. Meanwhile, this catalytic model led us to engineer a new SpCas9 variant (SpCas9-H982A + H983D) with reduced off-target effects. Thus, our study provided new mechanistic insights into the CRISPR-Cas9 system in the DNA cleavage state and offered useful guidance for engineering new CRISPR-Cas9 editing systems with improved specificity.

Insights into RNA-dependent RNA Polymerase Inhibitors as Antiinfluenza Virus Agents

BACKGROUND

Influenza is a seasonal disease that affects millions of people every year and has a significant economic impact. Vaccines are the best strategy to fight this viral pathology, but they are not always available or administrable, prompting the search for antiviral drugs. RNA-dependent RNA polymerase (RdRp) recently emerged as a promising target because of its key role in viral replication and its high conservation among viral strains.

DISCUSSION

This review presents an overview of the most interesting RdRp inhibitors that have been discussed in the literature since 2000. Compounds already approved or in clinical trials and a selection of inhibitors endowed with different scaffolds are described, along with the main features responsible for their activity.

RESULTS

RdRp inhibitors are emerging as a new strategy to fight viral infections and the importance of this class of drugs has been confirmed by the FDA approval of baloxavir marboxil in 2018. Despite the complexity of the RdRp machine makes the identification of new compounds a challenging research topic, it is likely that in the coming years, this field will attract the interest of a number of academic and industrial scientists because of the potential strength of this therapeutic approach.

R, Sissi C, Pecina A, Riccardi L, De Vivo M, Mancin F, Scrimin P. A Gold Nanoparticle Nanonuclease Relying on a Zn(II) Mononuclear Complex

Similarly to enzymes, functionalized gold nanoparticles efficiently catalyze chemical reactions, hence the term nanozymes. Herein, we present our results showing how surface-passivated gold nanoparticles behave as synthetic nanonucleases, able to cleave pBR322 plasmid DNA with the highest efficiency reported so far for catalysts based on a single metal ion mechanism. Experimental and computational data indicate that we have been successful in creating a catalytic site precisely mimicking that suggested for natural metallonucleases relying on a single metal ion for their activity. It comprises one Zn(II) ion to which a phosphate diester of DNA is coordinated. Importantly, as in nucleic acids-processing enzymes, a positively charged arginine plays a key role by assisting with transition state stabilization and by reducing the pKa of the nucleophilic alcohol of a serine. Our results also show how designing a catalyst for a model substrate (bis-p-nitrophenylphosphate) may provide wrong indications as for its efficiency when it is tested against the real target (plasmid DNA).

Catalytic Mechanism of Non-Target DNA Cleavage in CRISPR-Cas9 Revealed by Ab Initio Molecular Dynamics

CRISPR-Cas9 is a cutting-edge genome editing technology, which uses the endonuclease Cas9 to introduce mutations at desired sites of the genome. This revolutionary tool is promising to treat a myriad of human genetic diseases. Nevertheless, the molecular basis of DNA cleavage, which is a fundamental step for genome editing, has not been established. Here, quantum-classical molecular dynamics (MD) and free energy methods are used to disclose the two-metal-dependent mechanism of phosphodiester bond cleavage in CRISPR-Cas9. Ab initio MD reveals a conformational rearrangement of the Mg2+-bound RuvC active site, which entails the relocation of H983 to act as a general base. Then, the DNA cleavage proceeds through a concerted associative pathway fundamentally assisted by the joint dynamics of the two Mg2+ ions. This clarifies previous controversial experimental evidence, which could not fully establish the catalytic role of the conserved H983 and the metal cluster conformation. The comparison with other two-metal-dependent enzymes supports the identified mechanism and suggests a common catalytic strategy for genome editing and recombination. Overall, the non-target DNA cleavage catalysis described here resolves a fundamental open question in the CRISPR-Cas9 biology and provides valuable insights for improving the catalytic efficiency and the metal-dependent function of the Cas9 enzyme, which are at the basis of the development of genome editing tools.

Novel, Potent, and Druglike Tetrahydroquinazoline Inhibitor That Is Highly Selective for Human Topoisomerase II α over β

We disclose a novel class of 6-amino-tetrahydroquinazoline derivatives that inhibit human topoisomerase II (topoII), a validated target of anticancer drugs. In contrast to topoII-targeted drugs currently in clinical use, these compounds do not act as topoII poisons that enhance enzyme-mediated DNA cleavage, a mechanism that is linked to the development of secondary leukemias. Instead, these tetrahydroquinazolines block the topoII function with no evidence of DNA intercalation. We identified a potent lead compound [compound 14 (ARN-21934) IC₅₀ = 2 μM for inhibition of DNA relaxation, as compared to an IC₅₀ = 120 μM for the anticancer drug etoposide] with excellent metabolic stability and solubility. This new compound also shows ~100-fold selectivity for topoIIα over topoβ, a broad antiproliferative activity toward cultured human cancer cells, a favorable in vivo pharmacokinetic profile, and the ability to penetrate the blood–brain barrier. Thus, ARN-21934 is a highly promising lead for the development of novel and potentially safer topoII-targeted anticancer drugs.