Analysis of the stability and flexibility of RNA complexes containing bulge loops of different sizes

Further Probing of Cu2+-Dependent PNAzymes Acting as Artificial RNA Restriction Enzymes

Peptide nucleic acid (PNA)-neocuproine conjugates have been shown to efficiently catalyse the cleavage of RNA target sequences in the presence of Cu²⁺ ions in a site-specific manner. These artificial enzymes are designed to force the formation of a bulge in the RNA target, the sequence of which has been shown to be key to the catalytic activity. Here, we present a further investigation into the action of Cu²⁺-dependent PNAzymes with respect to the dependence on bulge composition in 3- and 4-nucleotide bulge systems. Cu²⁺-dependent PNAzymes were shown to have a clear preference for 4-nucleotide bulges, as the cleavage of 3-nucleotide bulge-forming RNA sequences was significantly slower, which is illustrated by a shift in the half-lives from approximately 30 min to 24 h. Nonetheless, the nucleotide preferences at different positions in the bulge displayed similar trends in both systems. Moreover, the cleavage site was probed by introducing critical chemical modifications to one of the cleavage site nucleotides of the fastest cleaved 4-nucleotide RNA bulge. Namely, the exclusion of the exocyclic amine of the central adenine and the replacement of the 2′-hydroxyl nucleophile with 2′-H or 2′-OMe substituents in the RNA severely diminished the rate of RNA cleavage by the Cu²⁺-dependent PNAzyme, giving insight into the mechanism of cleavage. Moreover, the shorter recognition arm of the RNA/PNAzyme complex was modified by extending the PNAzyme by two additional nucleobases. The new PNAzyme was able to efficiently promote the cleavage of RNA when fully hybridised to a longer RNA target and even outperform the previous fastest PNAzyme. The improvement was demonstrated in cleavage studies with stoichiometric amounts of either PNAzyme present, and the extended PNAzyme was also shown to give turnover with a 10-fold excess of the RNA target.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Influence of conjugation and other structural changes on the activity of Cu²⁺ based PNAzymes

We have previously shown that PNA-neocuproine conjugates can act as artificial RNA restriction enzymes. In the present study we have additionally conjugated the PNA with different entities, such as oligoethers, peptides etc. and also constructed systems where the PNA is designed to clamp the target RNA forming a triplex. Some conjugations are detrimental for the activity while most are silent which means that conjugation can be done to alter physical properties without losing activity. Conjugation with a single oligoether close to the neocuproine does enhance the rate almost twofold compared to the system without the oligoether. The systems designed to clamp the RNA target by forming a triplex retain the activity if the added oligoT sequence is 5 PNA units or shorter and extends the arsenal of artificial RNA restriction enzymes. Changing the direction of a closing base pair, where the target RNA forms a bulge, from a GC to a CG pair enhances the rate of cleavage somewhat without compromising the selectivity, leading to the so far most efficient artificial nuclease reported.

Characterizing the bending and flexibility induced by bulges in DNA duplexes

Advances in DNA nanotechnology have stimulated the search for simple motifs that can be used to control the properties of DNA nanostructures. One such motif, which has been used extensively in structures such as polyhedral cages, two-dimensional arrays, and ribbons, is a bulged duplex, that is, two helical segments that connect at a bulge loop. We use a coarse-grained model of DNA to characterize such bulged duplexes. We find that this motif can adopt structures belonging to two main classes: one where the stacking of the helices at the center of the system is preserved, the geometry is roughly straight, and the bulge is on one side of the duplex and the other where the stacking at the center is broken, thus allowing this junction to act as a hinge and increasing flexibility. Small loops favor states where stacking at the center of the duplex is preserved, with loop bases either flipped out or incorporated into the duplex. Duplexes with longer loops show more of a tendency to unstack at the bulge and adopt an open structure. The unstacking probability, however, is highest for loops of intermediate lengths, when the rigidity of single-stranded DNA is significant and the loop resists compression. The properties of this basic structural motif clearly correlate with the structural behavior of certain nano-scale objects, where the enhanced flexibility associated with larger bulges has been used to tune the self-assembly product as well as the detailed geometry of the resulting nanostructures. We further demonstrate the role of bulges in determining the structure of a "Z-tile," a basic building block for nanostructures.

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.

Impact of 2'-hydroxyl sampling on the conformational properties of RNA: update of the CHARMM all-atom additive force field for RNA

Here, we present an update of the CHARMM27 all-atom additive force field for nucleic acids that improves the treatment of RNA molecules. The original CHARMM27 force field parameters exhibit enhanced Watson-Crick base pair opening which is not consistent with experiment, whereas analysis of molecular dynamics (MD) simulations show the 2'-hydroxyl moiety to almost exclusively sample the O3' orientation. Quantum mechanical (QM) studies of RNA related model compounds indicate the energy minimum associated with the O3' orientation to be too favorable, consistent with the MD results. Optimization of the dihedral parameters dictating the energy of the 2'-hydroxyl proton targeting the QM data yielded several parameter sets, which sample both the base and O3' orientations of the 2'-hydroxyl to varying degrees. Selection of the final dihedral parameters was based on reproduction of hydration behavior as related to a survey of crystallographic data and better agreement with experimental NMR J-coupling values. Application of the model, designated CHARMM36, to a collection of canonical and noncanonical RNA molecules reveals overall improved agreement with a range of experimental observables as compared to CHARMM27. The results also indicate the sensitivity of the conformational heterogeneity of RNA to the orientation of the 2'-hydroxyl moiety and support a model whereby the 2'-hydroxyl can enhance the probability of conformational transitions in RNA.

Conformational analysis of genotoxic benzo[a]pyrene-7,8-dione-duplex DNA adducts using a molecular dynamics method (II)

The conformations of the benzo[a]pyrene-7,8-quinone (BPQ) modified oligonucleotide were investigated using molecular dynamic simulation. In the initial structures, the central guanine base was modified with BPQ resulting in the formation of four structurally distinguishable 10-(N2-deoxyguanosyl)-9,10-dihydro-9-hydroxy benzo[a]pyrene-7,8-dione adducts (BPQ-G3,4). Each of the oligonucleotide adduct consisted of two conformers, namely syn and anti conformations, depending on the rotation around the glycosidic bond between BPQ and the guanine base. The results revealed that the BPQ moiety was located in the major groove for all four syn conformers. The relative energies of these conformers were high, and the backbone largely deviated from the B-form. On the other hand, BPQ was located in the minor groove with relatively low energies, and backbone was retained in all of the anti conformer cases. The most conceivable BPQ-modified double stranded oligonucleotide structure was proposed from the energy calculation and the structural analysis.

RNA stability under different combinations of amber force fields and solvation models

The proper matching of force field and solvent is critical to obtain correct result in molecular dynamics simulation of bio-molecules. This problem has been intensively investigated for protein but not for RNA yet. In this paper, we use standard molecular dynamics and replica exchange molecular dynamics to take a series of tests on the RNA stability under different combinations of Amber force field parameters (ff98, ff99 and ff99bsc0) and the general Born implicit solvent models (igb1, igb2 and igb5). It is found that only ff98 and ff99bsc0 with igb1 can keep the native conformations of RNA hairpin and duplex. Our results suggest that ff98 plus igb1 may be reasonable choice for molecular dynamics simulation of RNA dynamics.

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.

Molecular dynamics of potential rRNA binders: single-stranded nucleic acids and some analogues

By hindering or "silencing" protein translation in vivo, antisense nucleic acid analogues that hybridize to bacterial rRNA could serve as a promising class of antibacterial compounds. Thus, we performed a comparative analysis of the dynamical properties of modified oligonucleotides based upon a sequence (5')r(UGUUACGACU)(3') that is complementary to bacterial ribosomal A-site RNA. In particular, 25 ns explicit solvent molecular dynamics simulations were computed for the following six single-stranded decamers: (1) the above RNA in unmodified form; (2) the 2'-O-methyl-modified RNA; (3) peptide nucleic acid (PNA) analogues of the above sequence, containing either (a) T or (b) U; and (4) two serine-substituted PNAs. Our results show that 2'-O-methylation attenuates RNA backbone dynamics, thereby preventing interconversion between stacked and unstacked conformations. The PNA analogue is rendered less flexible by replacing uracil with thymine; in addition, we found that derivatizing the PNA backbone with serine leads to enhanced base-stacking interactions. Consistent with known solubility properties of these classes of molecules, both RNAs exhibited greater localization of water molecules than did PNA. In terms of counterions, the initially helical conformation of the 2'-O-methyl RNA exhibits the highest Na(+) density among all the simulated decamers, while Na(+) build-up was most negligible for the neutral PNA systems. Further studies of the conformational and physicochemical properties of such modified single-stranded oligomers may facilitate better design of nucleic acid analogues, particularly those capable of serving as specific, high-affinity ribosomal A-site binders.

A theoretical elucidation of glucose interaction with HSA's domains

The interaction of different domains belonging to Human Serum Albumin (HSA) with open form of glucose have been investigated using molecular dynamics simulation methods. Applying docking, primary structures involving interaction of some residues with glucose have been obtained. Subsequently, equilibrium geometries at 300 K and minimum geometries have been determined for each of aforementioned structures by employing MD simulation and simulated annealing. The stability of species has been evaluated using a SAWSA v2.0 model. Ultimately, NBO analysis has been carried out to specify possible hydrogen bonding regarding the HSA interaction with glucose. Results obtained show that glucose can interact with Lys195, Lys199, and Glu153. In these interactions, each lysine forms an H-bonding with glucose. The H-bonding is obtained by stretching of N-H bond belonging to NH(3)(+) group of lysine along an oxygen atom of glucose. In addition, the above mentioned lysines are protonated, and there is an electrostatic interaction between glucose with Lys195 or Lys199. In addition, an H-bonding is formed between O atom of -COO group belonging to Glu153 and H atom of OH group belonging to glucose. Because, the N-H group of Lys195 interacts with the O atom of latter OH group, reaction of Lys195 is more desirable than that of Lys199. In fact, glucose is placed in the vicinity of Lys195 along with electrostatic interaction and H-bonding to Lys195 and Lys199 as well as H-bonding with Glu153, which subsequently reacts with Lys195. Thus, Lys195 is the primary site in reaction of glucose with HSA.

Insight derived from molecular dynamics simulation into substrate-induced changes in protein motions of proteinase K

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.

Molecular dynamics studies on troponin (TnI-TnT-TnC) complexes: insight into the regulation of muscle contraction

Mutations of any subunit of the troponin complex may lead to serious disorders. Rational approaches to managing these disorders require knowledge of the complex interactions among the three subunits that are required for proper function. Molecular dynamics (MD) simulations were performed for both skeletal (sTn) and cardiac (cTn) troponin. The interactions and correlated motions among the three components of the troponin complex were analyzed using both Molecular Mechanics-Generalized Born Surface Area (MMGBSA) and cross-correlation techniques. The TnTH2 helix was strongly positively correlated with the two long helices of TnI. The C domain of TnC was positively correlated with TnI and TnT. The N domain of TnC was negatively correlated with TnI and TnT in cTn, but not in sTn. The two C-domain calcium-binding sites of TnC were dynamically correlated. The two regulatory N-domain calcium-binding sites of TnC were dynamically correlated, even though the calcium-binding site I is dysfunctional. The strong interaction residue pairs and the strong dynamically correlated residues pairs among the three components of troponin complexes were identified. These correlated motions are consistent with the idea that there is a high degree of cooperativity among the components of the regulatory complex in response to Ca(2+) and other effectors. This approach may give insight into the mechanism by which mutations of troponin cause disease. It is interesting that some observed disease causing mutations fall within regions of troponin that are strongly correlated or interacted.

PNAzymes that are artificial RNA restriction enzymes

DNA-cleaving restriction enzymes are well-known tools in biomedical and biotechnological research. There are, however, no corresponding enzymes known for RNA cleavage. There has been an ongoing development of artificial ribonucleases, including some attempts at sequence selectivity. However, so far these systems have displayed modest rates of cleavage, and in most cases, the cleaver has been used in excess or in stoichiometric amounts. In the current work, we present PNA-based systems (PNAzymes) that carry a Cu(II)-2,9-dimethylphenanthroline group and that act as site and sequence specific RNases. The general basis for the systems is that the target is cleaved at a nonbase paired region (RNA bulge) which is formed in the substrate upon binding of the PNAzyme. With this copper based system, cleavage takes place at virtually only one site and with a half-life of down to 30 min under stoichiometric conditions. Efficient turnover of RNA-substrate is shown with a 100-fold excess of substrate, thus, demonstrating true enzyme behavior. In addition, alteration of the sequence in the RNA bulge or a mismatch in the base-pairing region leads to substantial decreases in rate showing both kinetic resolution and binding discrimination in the substrate selectivity. The selectivity is further demonstrated by the substrates, with two potential cleavage sites differing in only one base, are cleaved only at the site that either does not have a mismatch or is kinetically preferred. We suggest that these systems can serve as a basis for construction of RNA restriction enzymes for in vitro manipulations.

The therapeutically anti-prion active antibody-fragment scFv-W226: paramagnetic relaxation-enhanced NMR spectroscopy aided structure elucidation of the paratope-epitope interface

Antibodies have become indispensable reagents with numerous applications in biological and biotechnical analysis, in diagnostics as well as in therapy. In all cases, selective interaction with an epitope is crucial and depends on the conformation of the paratope. While epitopes are routinely mapped at high throughput, methods revealing structural insights on a rather short timescale are rare. We here demonstrate paramagnetic relaxation-enhanced (PRE) NMR spectroscopy to be a powerful tool unraveling structural information about epitope-orientation in a groove spanned by the complementary determining regions. In particular, we utilize the spin label TOAC, which is fused to the peptidic epitope using standard solid-phase chemistry and which is characterized by a reduced mobility compared to, e.g., spin labels attached to the side-chain functionalities of cysteine or lysine residues. We apply the method to determine the orientation of helix 1 of the prion protein, which is the epitope for the therapeutically anti-prion active scF(v) fragment W226.

Conformational and oligomeric effects on the cysteine pK(a) of tryparedoxin peroxidase

Typical 2-Cys peroxiredoxins (Prxs) are peroxidases which regulate cell signaling pathways, apoptosis, and differentiation. These enzymes are obligate homodimers, and can form decamers in solution. During catalysis, Prxs exhibit cysteine-dependent reactivity which requires the deprotonation of the peroxidatic cysteine (C(p)) supported by a lowered pK(a) in the initial step. We present the results of molecular dynamics simulations combined with pKa calculations on the monomeric, dimeric and decameric forms of one typical 2-Cys Prx, the tryparedoxin peroxidase from Trypanosoma cruzi (PDB id, 1uul). The calculations indicate that C(p) (C52) pK(a) values are highly affected by oligomeric state; an unshifted C(p) pK(a) (approximately 8.3, comparable to the pK(a) of isolated cysteine) is calculated for the monomer. In the dimers, starting with essentially identical structures, the C(p)s evolve dynamically asymmetric pK(a)s during the simulations; one subunit's C(p) pK(a) is shifted downward at a time, while the other is unshifted. However, when averaged over time, or multiple simulations, the two subunits within a dimer exhibit the same C(p), showing no preference for a lowered pK(a) in either subunit. Two conserved pathways that communicate the asymmetric pK(a)s between C(p)s of different subunits can be identified. In the decamer, all the C(p) pK(a)s are shifted downward, with slight asymmetry in the dimers which form the decamers. Structural analyses implicate oligomerization effects as responsible for these oligomeric state-dependent C(p) pK(a) shifts. The intra-dimer and the inter-dimer subunit contacts in the decamer restrict the conformations of the side chains of several residues (T49, T54 and E55) calculated to be key in shifting the C(p) pK(a). In addition, the backbone fluctuations of a few residues (M46, D47 and F48) result in a different electrostatic environment for the C(p) in dimers relative to the monomers. These side chain and backbone interactions which contribute to pK(a) modulation indicate the importance of oligomerization to the function of the typical 2-Cys Prxs.

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.