Design of LNA Analogues Using a Combined Density Functional Theory and Molecular Dynamics Approach for RNA Therapeutics.

Antisense therapeutics treat a wide spectrum of diseases, many of which cannot be addressed with the current drug technologies. In the quest to design better antisense oligonucleotide drugs, we propose five novel LNA analogues (A1-A5) for modifying antisense oligonucleotides and establishing each with the five standard nucleic acids: adenine (A), guanine (G), cytosine (C), thymine (T), and uracil (U). Monomer nucleotides of these modifications were considered for a detailed Density Functional Theory (DFT)-based quantum chemical analysis to determine their molecular-level structural and electronic properties. A detailed MD simulation study was done on a 14-mer ASO (5'-CTTAGCACTGGCCT-3') containing these modifications targeting PTEN mRNA. Results from both molecular- and oligomer-level analysis clearly depicted LNA-level stability of the modifications, the ASO/RNA duplexes maintaining stable Watson-Crick base pairing preferring RNA-mimicking A-form duplexes. Notably, monomer MO isosurfaces for both purines and pyrimidines were majorly distributed on the nucleobase region in modifications A1 and A2 and in the bridging unit in modifications A3, A4, and A5, suggesting that A3/RNA, A4/RNA, and A5/RNA duplexes interact more with the RNase H and solvent environment. Accordingly, solvation of A3/RNA, A4/RNA, and A5/RNA duplexes was higher compared to that of LNA/RNA, A1/RNA, and A2/RNA duplexes. This study has resulted in a successful archetype for creating advantageous nucleic acid modifications tailored for particular needs, fulfilling a useful purpose of designing novel antisense modifications, which may overcome the drawbacks and improve the pharmacokinetics of existing LNA antisense modifications.

Structural insight into locked nucleic acid based novel antisense modifications: A DFT calculations at monomer and MD simulations at oligomer level.

In the present study, five novel LNA based antisense modifications have been proposed. A conformational search was carried out using TANGO, followed by geometry optimization using MOPAC. Based on their electronic energies the most stable conformation for each modification was identified. Further, DFT based full geometry optimization on the most stable conformations at the gas phase B3LYP/6-31G(d,p) using a Gaussian03 and single point energy calculations on the optimized structures at the solvent phase B3LYP/6-311G(d,p) level of theory were done to derive their quantum chemical descriptors using the Gaussian09. A comparison of global reactivity descriptors confirmed that the LNA based modifications were the most reactive. Base-pair stability was recorded by observing the binding energies and base-pairing conformations of modified GC base pairs at the B3LYP/6-311G(d,p) level of theory. Molecular dynamics simulations have been performed at the oligomer duplex level by incorporating individual modifications on 20-mer RNA-RNA duplexes using AMBER16. Free energy calculations of duplex structures suggested that incorporation of A2 modification into the RNA-RNA duplex increased the duplex binding affinity similar to LNA. Whereas, the A3 modification showed less binding compared to LNA but improved binding compared to MOE. This computational approach using quantum chemical methods may be very useful to propose better modifications than the existing ones before performing the experiments in the area of antisense technology.

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations.

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6-41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA. Communicated by Ramaswamy H. Sarma.

Acetylcholinesterase and Aß Aggregation Inhibition by Heterometallic Ruthenium(II)-Platinum(II) Polypyridyl Complexes.

Two heteronuclear ruthenium(II)-platinum(II) complexes [Ru(bpy)2(BPIMBp)PtCl2]2+ (3) and [Ru(phen)2(BPIMBp)PtCl2]2+ (4), where bpy = 2,2'-bipyridine, phen = 1,10-phenanthroline, and BPIMBp = 1,4'-bis[(2-pyridin-2-yl)-1H-imidazol-1-ylmethyl]-1,1'-biphenyl, have been designed and synthesized from their mononuclear precursors [Ru(bpy)2(BPIMBp)]2+ (1) and [Ru(phen)2(BPIMBp)]2+ (2) as multitarget molecules for Alzheimer's disease (AD). The inclusion of the cis-PtCl2 moiety facilitates the covalent interaction of Ru(II) polypyridyl complexes with amyloid ß (Aß) peptide. These multifunctional complexes act as inhibitors of acetylcholinesterase (AChE), Aß aggregation, and Cu-induced oxidative stress and protect neuronal cells against Aß-toxicity. The study highlights the design of metal based anti-Alzheimer's disease (AD) systems.

Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study.

Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.

REMD and umbrella sampling simulations to probe the energy barrier of the folding pathways of engrailed homeodomain.

Proteins fold by diverse pathways which depend on the energy barriers involved in reaching different intermediates. There has been a lot of development in the theoretical aspects of protein folding, from force-field to simulation techniques. One such simulation approach is replica exchange molecular dynamics simulation (REMD), which provides an efficient conformational sampling method to understand the events involved in protein folding. In this study, an attempt is made to explore the folding funnel of engrailed homeodomain protein (EnHD) using REMD simulations. EnHD is a 54 residue long helix bundle protein which has a folding time of about 15 µs. The protein was represented using the Amber United atom model in order to reduce the system size which helped to speed up the simulation. Individual replicas were simulated for 1.4-2 µs making cumulative time of more than 100 µs of REMD simulations. Free energy analysis was carried out to understand the folding behavior of EnHD protein. Effects of temperature range and exchange frequency in REMD simulations have been explored. In addition to this, multiple umbrella sampling (US) simulations of a total of 320 ns were also carried out, followed by weighted histogram analysis method (WHAM) to investigate the energy barriers involved during the folding of various intermediates. US studies were also carried on mutational variants of EnHD protein to see effect of the mutations on the folding pathway of the protein. The use of US technique may be helpful for predicting fast folding mutants or protein engineering. The combination of REMD with US may help in understanding the energetics between multiple pathways of fast folding proteins and their mutant counterparts.

MD simulations of HIV-1 RT primer-template complex: effect of modified nucleosides and antisense PNA oligomer.

Human immunodeficiency virus type 1 (HIV-1) requires the human tRNA(3)(Lys) as a reverse transcriptase (RT) primer. The annealing of 3' terminal 18 nucleotides of tRNA(3)(Lys) with the primer binding site (PBS) of viral RNA (vRNA) is crucial for reverse transcription. Additional contacts between the A rich (A-loop) region of vRNA and the anticodon domain of tRNA(3)(Lys) are necessary, which show the specific requirement of tRNA(3)(Lys). The importance of modified nucleosides, present in tRNA(3)(Lys), in giving stability to the primer-template complex has been determined in earlier experiments. It has been observed that the PNA oligomer targeted to PBS of vRNA destabilized the crucial interactions between primer and template due to which the reverse transcription is inhibited. Molecular dynamics simulations have been carried out to study the effect of modified nucleosides on the vRNA-tRNA(3)(Lys) complex stability and the destabilization effect of PNA oligomer on the vRNA-tRNA(3)(Lys)-PNA complex. The root-mean-square deviation, hydrogen bonding, tertiary interactions, and free energy calculations of the simulation data support the experimental results. The analyses have revealed the structural changes in PBS region of vRNA which might be another strong reason for the inability of RT binding to 7F helix for its normal functioning of reverse transcription.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape.

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

Molecular dynamics simulations of cyclohexyl modified peptide nucleic acids (PNA).

Peptide Nucleic Acids (PNA) that bind sequence specifically to DNA/RNA are of major interest in the field of molecular biology and could form the basis for gene-targeted drugs. Molecular dynamics simulations are aimed to characterize the structural and dynamical features to understand the effect of backbone modification on the structure and dynamics along with the stability of the resulting 10mer complexes of PNA with DNA/RNA. Twelve Molecular Dynamics (MD) simulations of duplexes and triplexes with and without cyclohexyl modification were carried out for 10ns each. The simulations indicate that the cyclohexyl modification with different stereoisomers has influenced all the PNA-DNA/RNA complexes. Modification has added rigidity to backbone by restricting beta to +60 in case of (1R,2S) cyclohexyl PNA and to -60 in case of (1S,2R) cyclohexyl PNA. The results of MD simulations were able to show the backbone rigidification and preference for RNA complexes over DNA due to presence of cyclohexyl ring in the PNA backbone.

Study of early events in the protein folding of villin headpiece using molecular dynamics simulation.

Protein folding is scientifically and computationally challenging problem. The early phases of protein folding are interesting due to various events like nascent secondary structure formation, hydrophobic collapse leading to formation of non-native or meta-stable conformations. These events occur within a very short time span of 100 ns as compared to total folding time of few microseconds. It is highly difficult to observe these events experimentally due to very short lifetime. Molecular dynamics simulation technique can efficiently probe the detailed atomic level understanding about these events. In the present paper, all atom molecular dynamics simulation trajectory of nearly 200 ns was carried out for fully solvated villin headpiece with PME treatment using AMBER 7 package. Initial hydrophobic collapse along with secondary structure formation resulted into formation of partially stable non-native conformations. The formation of secondary structural elements and hydrophobic collapse takes place simultaneously in the folding process.

Conformational preferences of the base substituent in hypermodified nucleotide queuosine 5'-monophosphate 'pQ' and protonated variant 'pQH+'.

Conformational preferences of the base substituent in hypermodified nucleotide queuosine 5'-monophosphate 'pQ' and its protonated form 'pQH+' have been studied using quantum chemical Perturbative Configuration Interaction with Localized Orbitals PCILO method. The salient points have also been examined using molecular mechanics force field MMFF, parameterized modified neglect of differential overlap PM3 and Hartree Fock-Density Functional Theory HF DFT (pBP/DN*) approaches. Aqueous solvation of pQ and pQH+ has also been studied using molecular dynamics simulations. Consistent with the observed crystal structure, in isolated protonated form pQH+, the quaternary amine HN(13)(+)H, of the sidechain having 7-aminomethyl linkage, hydrogen bonds with the carbonyl oxygen O(10) of the base. However, N(13)H-O(10) hydrogen bonding is not preferred for unprotonated pQ, whether isolated or hydrated. Interaction between the 5'-phosphate and the 7-aminomethyl group is more likely for isolated pQ. The cyclopentenediol hydroxyl group O4"H may hydrogen bond with the O(10) in isolated pQ as well as in pQH+. The O4"H may hydrogen bond with the 5'-phosphate as well. The presence of -CH2-NH- and O"H groups in pQ and pQH+ allows interesting possibilities for intranucleotide hydrogen bonds and interactions across the anticodon loop. Simultaneous hydrogen bonds O2P-HN(13)+H-O(10) are indicated for hydrated pQH+. Unlike weak involvement of O4"H, these interactions also persist in hydrated pQH+ and may much reduce backbone flexibility. Resulting sub-optimal Q:C base pairing leads to unbiased reading of U or C as the third codon letter. Cyclopentenediol hydroxyl groups may interact with other biomolecules, allowing specific recognition. Prospective pQ(34) and pQ(34)H+ sites for codon-anticodon base pairing remain unhindered, but non canonical Q:G base pairing (amber-suppression) is ruled out.

Conformational preferences of anticodon 3'-adjacent hypermodified nucleic acid base cis-or trans-zeatin and its 2-methylthio derivative, cis-or trans- ms(2)zeatin.

Conformational preferences of the hypermodified nucleic acid bases N6-(Delta(2)-cis-hydroxyisopentenyl)adenine, cis-io(6)Ade also known as cis-zeatin, and N(6)-(Delta(2)-trans-hydroxyisopentenyl)adenine, trans-io(6)ade or trans-zeatin, and 2-methylthio derivatives of these cis-ms(2)io(6)Ade or cis-ms(2)zeatin, and trans-ms(2)io6Ade or trans-ms(2)zeatin have been investigated theoretically by the quantum chemical Perturbative Configuration Interaction with Localized Orbitals (PCILO) method. Automated geometry optimization using quantum chemical MNDO, AM1 and PM3 methods has also been made to compare the salient features. The predicted most stable conformation of cis-io(6)Ade, trans-io(6)Ade, cis-ms(2)io(6)Ade and trans-ms(2)io(6)Ade are such that in each of these molecules the isopentenyl substituent spreads away (has "dista" conformation) from the five membered ring imidazole moiety of the adenine. The atoms N(6), C(10) and C(11) remain coplanar with the adenine ring in the predicted preferred conformation for each of these molecules. In cis-io(6)Ade as well as cis-ms(2)io(6)Ade the hydroxyl oxygen may participate in intramolecular hydrogen bonding with the H-C(10)-H group. In trans-io(6)Ade the hydroxyl group is oriented towards the H-C(2) instead. This orientation is retained in trans-ms(2)io(6)Ade, possible O-H...S hydrogen bonding may be a stabilizing factor. In all these four modified adenines C(11)-H is favourably placed to participate in intramolecular hydrogen bonding with N(1). In cis-ms(2)io(6)Ade as well as trans-ms(2)io(6)Ade the 2-methylthio group preferentially orients on the same side as C(2)-N(3) bond, due to this non-obstrusive placing, orientation of the hydroxyisopentenyl substituent remains unaffected by 2-methylthiolation. Thus the N(1) site remains shielded irrespective of the 2-methylthiolation status in these various cis-and trans-zeatin analogs alike. Firmly held orientation of hydroxyisopentenyl substituent in zeatin isomers and derivatives, in contrast to adaptable orientation of isopentenyl substituent in i(6)Ade and ms(2)i(6)Ade, may account for the increased efficiency of suppressor tRNA and reduced codon context sensitivity accompanied with the occurrence of ms(2)-zeatin (ms(2)io(6)Ade) modification.