Protein environment affects the water-tryptophan binding mode. MD, QM/MM, and NMR studies of engrailed homeodomain mutants.

Water molecules can interact with aromatic moieties using either their O-H bonds or their lone-pairs of electrons. In proteins, water-π interactions have been reported to occur with tryptophan and histidine residues, and dynamic exchange between O-Hπ hydrogen bonding and lone-pairπ interactions was suggested to take place, based on ab initio calculations. Here we used classical and QM/MM molecular dynamics simulations, complemented with an NMR study, to examine a specific water-indole interaction observed in the engrailed homeodomain and in its mutants.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations.

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA.

Single Stranded Loops of Quadruplex DNA As Key Benchmark for Testing Nucleic Acids Force Fields.

We have carried out a set of explicit solvent molecular dynamics (MD) simulations on two DNA quadruplex (G-DNA) molecules, namely the antiparallel d(G4T4G4)2 dimeric quadruplex with diagonal loops and the parallel-stranded human telomeric monomolecular quadruplex d[AGGG(TTAGGG)3] with three propeller loops. The main purpose of the paper was testing of the capability of the MD simulation technique to describe single-stranded topologies of G-DNA loops, which represent a very challenging task for computational methods. The total amount of conventional and locally enhanced sampling (LES) simulations analyzed in this study exceeds 1.

Molecular dynamics simulations and their application to four-stranded DNA.

This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields.

Impact of an extruded nucleotide on cleavage activity and dynamic catalytic core conformation of the hepatitis delta virus ribozyme.

The self-cleaving hepatitis delta virus (HDV) ribozyme is essential for the replication of HDV, a liver disease causing pathogen in humans. The catalytically critical nucleotide C75 of the ribozyme is buttressed by a trefoil turn pivoting around an extruded G76. In all available crystal structures, the conformation of G76 is restricted by stacking with G76 of a neighboring molecule.

Calculation of structural behavior of indirect NMR spin-spin couplings in the backbone of nucleic acids.

Calculated indirect NMR spin-spin coupling constants (J-couplings) between (31)P, (13)C, and (1)H nuclei were related to the backbone torsion angles of nucleic acids (NAs), and it was shown that J-couplings can facilitate accurate and reliable structural interpretation of NMR measurements and help to discriminate between their distinct conformational classes. A proposed stepwise procedure suggests assignment of the J-couplings to torsion angles from the sugar part to the phosphodiester link. Some J-couplings show multidimensional dependence on torsion angles, the most prominent of which is the effect of the sugar pucker.

Molecular dynamics simulations of RNA: an in silico single molecule approach.

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application.

Principles of RNA base pairing: structures and energies of the trans Watson-Crick/sugar edge base pairs.

Due to the presence of the 2'-OH hydroxyl group of ribose, RNA molecules utilize an astonishing variability of base pairing patterns to build up their structures and perform the biological functions. Many of the key RNA base pairing families have no counterparts in DNA. In this study, the trans Watson-Crick/sugar edge (trans WC/SE) RNA base pair family has been characterized using quantum chemical and molecular mechanics calculations.

Non-Watson-Crick base pairing in RNA. quantum chemical analysis of the cis Watson-Crick/sugar edge base pair family.

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/6-31G and RIMP2/aug-cc-pVDZ calculations.

Molecular dynamics simulations of sarcin-ricin rRNA motif.

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented.

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations.

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis.

Theoretical study of the guanine --> 6-thioguanine substitution in duplexes, triplexes, and tetraplexes.

Molecular dynamics and thermodynamic integration calculations have been carried out on a set of G-rich single-strand, duplex, triplex, and quadruplex DNAs to study the structural and stability changes connected with the guanine --> 6-thioguanine (G --> S) mutation. The presence of 6-thioguanine leads to a shift of the geometry from the B/A intermediate to the pure B-form in duplex DNA. The G --> S mutation does not largely affect the structure of the antiparallel triplex when it is located at the reverse-Hoogsteen position, but leads to a non-negligible local distortion in the structure when it is located at the Watson-Crick position.

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex.

Molecular dynamics simulations of RNA-protein complex between Escherichia coli loop E/helix IV (LE/HeIV) rRNA and L25 protein reveal a qualitative agreement between the experimental and simulated structures. The major groove of LE is a prominent rRNA cation-binding site. Divalent cations rigidify the LE major groove geometry whereas in the absence of divalent cations LE extensively interacts with monovalent cations via inner-shell binding.

Ribosomal RNA kink-turn motif--a flexible molecular hinge.

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries.

Molecular dynamics simulations of Guanine quadruplex loops: advances and force field limitations.

A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures.

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets.

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs.

Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates.

The formation of a cation-stabilized guanine quadruplex (G-DNA) stem is an exceptionally slow process involving complex kinetics that has not yet been characterized at atomic resolution. Here, we investigate the formation of a parallel stranded G-DNA stem consisting of four strands of d(GGGG) using molecular dynamics simulations with explicit inclusion of counterions and solvent. Due to the limitations imposed by the nanosecond timescale of the simulations, rather than watching for the spontaneous formation of G-DNA, our approach probes the stability of possible supramolecular intermediates (including two-, three-, and four-stranded assemblies with out-of-register base pairing between guanines) on the formation pathway.

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs.

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out.

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E.

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence.

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution.

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively.

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts.

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove.