Complex between Human RNase HI and the phosphonate-DNA/RNA duplex: Molecular dynamics study

Potential repurposing of four FDA approved compounds with antiplasmodial activity identified through proteome scale computational drug discovery and in vitro assay

Malaria elimination can benefit from time and cost-efficient approaches for antimalarials such as drug repurposing. In this work, 796 DrugBank compounds were screened against 36 Plasmodium falciparum targets using QuickVina-W. Hits were selected after rescoring using GRaph Interaction Matching (GRIM) and ligand efficiency metrics: surface efficiency index (SEI), binding efficiency index (BEI) and lipophilic efficiency (LipE). They were further evaluated in Molecular dynamics (MD). Twenty-five protein-ligand complexes were finally retained from the 28,656 (36 × 796) dockings. Hit GRIM scores (0.58 to 0.78) showed their molecular interaction similarity to co-crystallized ligands. Minimum LipE (3), SEI (23) and BEI (7) were in at least acceptable thresholds for hits. Binding energies ranged from -6 to -11 kcal/mol. Ligands showed stability in MD simulation with good hydrogen bonding and favorable protein-ligand interactions energy (the poorest being -140.12 kcal/mol). In vitro testing showed 4 active compounds with two having IC50 values in the single-digit μM range.

Tuning the hybridization properties of modified oligonucleotides: from flexible to conformationally constrained phosphonate internucleotide linkages

The concept of conformational restriction leading to the preorganization of modified strands has proven to be successful and has afforded nucleic acid analogues with many interesting properties suitable for various biochemical applications. We utilized this concept to prepare a set of constrained oligonucleotides derived from 1,4-dioxane and 1,3-dioxolane-locked nucleoside phosphonates and evaluated their hybridization affinities towards their complementary RNA strands. With an increase of ΔT_m per modification up to +5.2 °C, the hybridization experiments revealed the (S)-2',3'-O-phosphonomethylidene internucleotide linkage as one of the most T_m-increasing modifications reported to date. Moreover, we introduced a novel prediction tool for the pre-selection of potentially interesting chemical modifications of oligonucleotides.

Stacking in RNA: NMR of Four Tetramers Benchmark Molecular Dynamics

Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against 1H-1H NOESY distances and 3J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.

Evidence from molecular dynamics simulations of conformational preorganization in the ribonuclease H active site

Ribonuclease H1 (RNase H) enzymes are well-conserved endonucleases that are present in all domains of life and are particularly important in the life cycle of retroviruses as domains within reverse transcriptase. Despite extensive study, especially of the E. coli homolog, the interaction of the highly negatively charged active site with catalytically required magnesium ions remains poorly understood. In this work, we describe molecular dynamics simulations of the E. coli homolog in complex with magnesium ions, as well as simulations of other homologs in their apo states. Collectively, these results suggest that the active site is highly rigid in the apo state of all homologs studied and is conformationally preorganized to favor the binding of a magnesium ion. Notably, representatives of bacterial, eukaryotic, and retroviral RNases H all exhibit similar active-site rigidity, suggesting that this dynamic feature is only subtly modulated by amino acid sequence and is primarily imposed by the distinctive RNase H protein fold.

5'-O-Methylphosphonate nucleic acids--new modified DNAs that increase the Escherichia coli RNase H cleavage rate of hybrid duplexes

Several oligothymidylates containing various ratios of phosphodiester and isopolar 5'-hydroxyphosphonate, 5'-O-methylphosphonate and 3'-O-methylphosphonate internucleotide linkages were examined with respect to their hybridization properties with oligoriboadenylates and their ability to induce RNA cleavage by ribonuclease H (RNase H). The results demonstrated that the increasing number of 5'-hydroxyphosphonate or 5'-O-methylphosphonate units in antisense oligonucleotides (AOs) significantly stabilizes the heteroduplexes, whereas 3'-O-methylphosphonate AOs cause strong destabilization of the heteroduplexes. Only the heteroduplexes with 5'-O-methylphosphonate units in the antisense strand exhibited a significant increase in Escherichia coli RNase H cleavage activity by up to 3-fold (depending on the ratio of phosphodiester and phosphonate linkages) in comparison with the natural heteroduplex. A similar increase in RNase H cleavage activity was also observed for heteroduplexes composed of miRNA191 and complementary AOs containing 5'-O-methylphosphonate units. We propose for this type of AOs, working via the RNase H mechanism, the abbreviation MEPNA (MEthylPhosphonate Nucleic Acid).