Nucleic acids in modern molecular therapies: A realm of opportunities for strategic drug design.

RNA vaccines have made evident to society what was already known by the scientific community: nucleic acids will be the "drugs of the future." By modifying the genome, interfering in transcription or translation, and by introducing new catalysts into the cell or by mimicking antibody effects, nucleic acids can generate therapeutic activities that are not accessible by any other therapeutic agents. There are, however, challenges that need to be solved in the next few years to make nucleic acids usable in a wide range of therapeutic scenarios. This review illustrates how simulation methods can help achieve this goal.

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.

Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case.

Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.

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) IC50 = 2 µM for inhibition of DNA relaxation, as compared to an IC50 = 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.

Visualizing group II intron dynamics between the first and second steps of splicing.

Group II introns are ubiquitous self-splicing ribozymes and retrotransposable elements evolutionarily and chemically related to the eukaryotic spliceosome, with potential applications as gene-editing tools. Recent biochemical and structural data have captured the intron in multiple conformations at different stages of catalysis. Here, we employ enzymatic assays, X-ray crystallography, and molecular simulations to resolve the spatiotemporal location and function of conformational changes occurring between the first and the second step of splicing. We show that the first residue of the highly-conserved catalytic triad is protonated upon 5'-splice-site scission, promoting a reversible structural rearrangement of the active site (toggling). Protonation and active site dynamics induced by the first step of splicing facilitate the progression to the second step. Our insights into the mechanism of group II intron splicing parallels functional data on the spliceosome, thus reinforcing the notion that these evolutionarily-related molecular machines share the same enzymatic strategy.

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.

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.

Chemoselective Synthesis of N-Terminal Cysteinyl Thioesters via ß,γ-C,S Thiol-Michael Addition.

Dehydroalanine (ΔAla) is a highly electrophilic residue that can react efficiently with sulfur nucleophiles to furnish cysteinyl analogues. Herein, we report an efficient synthesis of N-terminal cysteinyl thioesters, suitable for S, N-acyl transfer, based on ß,γ-C,S thiol-Michael addition. Both ionic and radical-based methodologies were found to be efficient for this process.

A synthetic genetic polymer with an uncharged backbone chemistry based on alkyl phosphonate nucleic acids.

The physicochemical properties of nucleic acids are dominated by their highly charged phosphodiester backbone chemistry. This polyelectrolyte structure decouples information content (base sequence) from bulk properties, such as solubility, and has been proposed as a defining trait of all informational polymers. However, this conjecture has not been tested experimentally. Here, we describe the encoded synthesis of a genetic polymer with an uncharged backbone chemistry: alkyl phosphonate nucleic acids (phNAs) in which the canonical, negatively charged phosphodiester is replaced by an uncharged P-alkyl phosphonodiester backbone. Using synthetic chemistry and polymerase engineering, we describe the enzymatic, DNA-templated synthesis of P-methyl and P-ethyl phNAs, and the directed evolution of specific streptavidin-binding phNA aptamer ligands directly from random-sequence mixed P-methyl/P-ethyl phNA repertoires. Our results establish an example of the DNA-templated enzymatic synthesis and evolution of an uncharged genetic polymer and provide a foundational methodology for their exploration as a source of novel functional molecules.

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.

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.

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

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

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

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

Antigenic peptide molecular recognition by the DRB1-DQB1 haplotype modulates multiple sclerosis susceptibility.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system that has a notably high incidence in Sardinia. Our study focuses on two HLA class II haplotypes associated with the disease in Sardinia, the rare predisposing DRB1*15:01-DQB1*06:02 and the widespread protective DRB1*16:01-DQB1*05:02. This framework enabled the highlighting of HLA binding pocket specificity and peptide recognition mechanisms by employing molecular dynamics simulations of the whole DRB1-DQB1 haplotype interacting with MBP- and EBV-derived peptides. We analyzed peptide-protein interaction networks and temporal evolution of the original complexes and after key amino acid mutations. The mutation G86V of the protective DRB1 allele exerted its effect mainly in the presence of the EBV viral peptide, with local and long range outcomes. However, the V38A mutation of the protective DQB1 showed a long range effect only in the case of the MBP myelin peptide. Our findings also demonstrate a DRB1/DQB1 complementary molecular recognition of peptides. This mechanism could provide a robust synergistic action and a differential role of DRB1 and DQB1 in tissues and in the time-steps towards autoimmunity. In addition, we demonstrate that negatively charged residues in pockets 4 and 9 play a role in MS susceptibility. Our findings are supported by recent experiments using a closely related MS animal model. Overall, our analysis confirms the role of the DRB1-DQB1 haplotype in conferring disease predisposition and could provide a valuable aid in designing optimal therapeutic peptides for MS therapy.