Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.

• MeSH
• fosfáty * MeSH
• molekulární konformace MeSH
• nukleotidové motivy MeSH
• RNA * chemie   MeSH
• simulace molekulární dynamiky MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• fosfáty * MeSH
• RNA * MeSH

Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.

• MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• RNA MeSH

Neoscytalidium (or N.) dimidiatum and N. novaehollandiae are two aggressive plant pathogenic species that affect several agricultural crops. Early detection and identification of these fungi are of critical importance to bring about the effective minimization to the threat they pose to the infected plants. Herein, two species of Neoscytalidium were rapidly discriminated by utilizing the rRNA internal transcribed (ITS4-5.8S-ITS5) PCR primers. A total of 100 isolates of Neoscytalidium species, which were isolated from Iraqi canker-infected fig trees, were included in this study. Two discrete electrophoretic PCR bands were observed in Neoscytalidium isolates-A-variants were about 546 bp, while B-variants were about 993 bp in length. The comprehensive phylogenetic analysis of both DNA variants revealed that A-variants resided between N. novaehollandiae and N. hyalinum, while B-variants were closely related to N. dimidiatum. Furthermore, the highly specific re-constructed tree of both electrophoretic variants demonstrated that B-variants share a high similarity with N. novaehollandiae. Additionally, the secondary structures for both variants were predicted computationally to reveal the structural patterns that each variant follows. In conclusion, a small rRNA locus comprising 22 nucleotides that differs in the two variants is potentially responsible for this species-specific classification. The main divergence in the amplified loci led to the classification of these fungal variants into two main species, namely N. dimidiatum and N. novaehollandiae, demonstrating that the amplification by ITS4-ITS5 rRNA fragment is a beneficial strategy that can be employed for the assessment of Neoscytalidium diversity in the natural ecosystems.

• MeSH
• Ascomycota klasifikace   genetika   izolace a purifikace   MeSH
• DNA fungální genetika   MeSH
• druhová specificita MeSH
• elektroforéza v agarovém gelu * MeSH
• Ficus mikrobiologie   MeSH
• fylogeneze MeSH
• genetická variace * MeSH
• mezerníky ribozomální DNA chemie   genetika   MeSH
• nemoci rostlin mikrobiologie   MeSH
• polymerázová řetězová reakce MeSH
• RNA ribozomální 5.8S chemie   genetika   MeSH
• sekvenční analýza DNA MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• DNA fungální MeSH
• mezerníky ribozomální DNA MeSH
• RNA ribozomální 5.8S MeSH

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.

• MeSH
• DNA chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny * MeSH
• počítačová simulace MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• DNA MeSH
• RNA MeSH

The computer-aided folding of biomolecules, particularly RNAs, is one of the most difficult challenges in computational structural biology. RNA tetraloops are fundamental RNA motifs playing key roles in RNA folding and RNA-RNA and RNA-protein interactions. Although state-of-the-art Molecular Dynamics (MD) force fields correctly describe the native state of these tetraloops as a stable free-energy basin on the microsecond time scale, enhanced sampling techniques reveal that the native state is not the global free energy minimum, suggesting yet unidentified significant imbalances in the force fields. Here, we tested our ability to fold the RNA tetraloops in various force fields and simulation settings. We employed three different enhanced sampling techniques, namely, temperature replica exchange MD (T-REMD), replica exchange with solute tempering (REST2), and well-tempered metadynamics (WT-MetaD). We aimed to separate problems caused by limited sampling from those due to force-field inaccuracies. We found that none of the contemporary force fields is able to correctly describe folding of the 5'-GAGA-3' tetraloop over a range of simulation conditions. We thus aimed to identify which terms of the force field are responsible for this poor description of TL folding. We showed that at least two different imbalances contribute to this behavior, namely, overstabilization of base-phosphate and/or sugar-phosphate interactions and underestimated stability of the hydrogen bonding interaction in base pairing. The first artifact stabilizes the unfolded ensemble, while the second one destabilizes the folded state. The former problem might be partially alleviated by reparametrization of the van der Waals parameters of the phosphate oxygens suggested by Case et al., while in order to overcome the latter effect we suggest local potentials to better capture hydrogen bonding interactions.

• MeSH
• konformace nukleové kyseliny MeSH
• RNA chemie   metabolismus   MeSH
• sbalování RNA MeSH
• simulace molekulární dynamiky * MeSH
• stabilita RNA MeSH
• statická elektřina MeSH
• teplota MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• RNA MeSH

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

• MeSH
• konformace nukleové kyseliny MeSH
• kvantová teorie * MeSH
• nukleové kyseliny chemie   MeSH
• RNA chemie   MeSH
• výpočetní biologie * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• nukleové kyseliny MeSH
• RNA MeSH

Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair and finally by a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and base-phosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 µs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by MM-PBSA free energy calculations, locally enhanced sampling (LES) runs, targeted MD (TMD) and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and net-neutralizing Na(+) vs. excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the major groove, loss of base-phosphate interaction and other instabilities. The alternative X-ray geometry even undergoes conformational transition towards the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.

• Publikační typ
• časopisecké články MeSH

BACKGROUND: The first systematic study of small non-coding RNAs (sRNA, ncRNA) in Streptomyces is presented. Except for a few exceptions, the Streptomyces sRNAs, as well as the sRNAs in other genera of the Actinomyces group, have remained unstudied. This study was based on sequence conservation in intergenic regions of Streptomyces, localization of transcription termination factors, and genomic arrangement of genes flanking the predicted sRNAs. RESULTS: Thirty-two potential sRNAs in Streptomyces were predicted. Of these, expression of 20 was detected by microarrays and RT-PCR. The prediction was validated by a structure based computational approach. Two predicted sRNAs were found to be terminated by transcription termination factors different from the Rho-independent terminators. One predicted sRNA was identified computationally with high probability as a Streptomyces 6S RNA. Out of the 32 predicted sRNAs, 24 were found to be structurally dissimilar from known sRNAs. CONCLUSION: Streptomyces is the largest genus of Actinomyces, whose sRNAs have not been studied. The Actinomyces is a group of bacterial species with unique genomes and phenotypes. Therefore, in Actinomyces, new unique bacterial sRNAs may be identified. The sequence and structural dissimilarity of the predicted Streptomyces sRNAs demonstrated by this study serve as the first evidence of the uniqueness of Actinomyces sRNAs.

• MeSH
• algoritmy MeSH
• bakteriální RNA chemie   genetika   MeSH
• druhová specificita MeSH
• genom bakteriální MeSH
• intergenová DNA MeSH
• konformace nukleové kyseliny MeSH
• molekulární modely MeSH
• nekódující RNA chemie   genetika   MeSH
• polymerázová řetězová reakce s reverzní transkripcí MeSH
• sekvence nukleotidů MeSH
• sekvenční analýza hybridizací s uspořádaným souborem oligonukleotidů MeSH
• Streptomyces coelicolor genetika   MeSH
• Streptomyces genetika   MeSH
• terminátorové oblasti (genetika) MeSH
• výpočetní biologie MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• validační studie MeSH
• Názvy látek
• bakteriální RNA MeSH
• intergenová DNA MeSH
• nekódující RNA MeSH