Influence of Spermine on DNA Conformation in a Molecular Dynamics Trajectory of d(CGCGAATTCGCG)2: Major Groove Binding by One Spermine Molecule Delays the A→B Transition

Molecular Dynamics Simulations of Dna and Its Complexes

This article describes how classical molecular simulation methods are being used to gain a molecular-level understanding of the interaction mechanisms responsible for DNA–ligand recognition, and that govern the response of DNA to ligand binding. Case studies using a variety of different ligands—including small pharmaceutical drugs, proteins and lipids—are used to illustrate the power of modern molecular dynamics simulation methods for understanding how we may control the function and structure of DNA.

Single-molecule portrait of DNA and RNA double helices

The composition and geometry of the genetic information carriers were described as double-stranded right helices sixty years ago. The flexibility of their sugar-phosphate backbones and the chemistry of their nucleotide subunits, which give rise to the RNA and DNA polymers, were soon reported to generate two main structural duplex states with biological relevance: the so-called A and B forms. Double-stranded (ds) RNA adopts the former whereas dsDNA is stable in the latter. The presence of flexural and torsional stresses in combination with environmental conditions in the cell or in the event of specific sequences in the genome can, however, stabilize other conformations. Single-molecule manipulation, besides affording the investigation of the elastic response of these polymers, can test the stability of their structural states and transition models. This approach is uniquely suited to understanding the basic features of protein binding molecules, the dynamics of molecular motors and to shedding more light on the biological relevance of the information blocks of life. Here, we provide a comprehensive single-molecule analysis of DNA and RNA double helices in the context of their structural polymorphism to set a rigorous interpretation of their material response both inside and outside the cell. From early knowledge of static structures to current dynamic investigations, we review their phase transitions and mechanochemical behaviour and harness this fundamental knowledge not only through biological sciences, but also for Nanotechnology and Nanomedicine.

Prediction of hammerhead ribozyme intracellular activity with the catalytic core fingerprint

Hammerhead ribozyme is a versatile tool for down-regulation of gene expression in vivo. Owing to its small size and high activity, it is used as a model for RNA structure–function relationship studies. In the present paper we describe a new extended hammerhead ribozyme HH-2 with a tertiary stabilizing motif constructed on the basis of the tetraloop receptor sequence. This ribozyme is very active in living cells, but shows low activity in vitro. To understand it, we analysed tertiary structure models of substrate–ribozyme complexes. We calculated six unique catalytic core geometry parameters as distances and angles between particular atoms that we call the ribozyme fingerprint. A flanking sequence and tertiary motif change the geometry of the general base, general acid, nucleophile and leaving group. We found almost complete correlation between these parameters and the decrease of target gene expression in the cells. The tertiary structure model calculations allow us to predict ribozyme intracellular activity. Our approach could be widely adapted to characterize catalytic properties of other RNAs.

Molecular dynamics simulations of DNA-polycation complex formation

Complexes formed from DNA and polycations are of interest because of their potential use in gene therapy; however, there remains a lack of understanding of the structure and formation of DNA-polycation complexes at atomic scale. In this work, molecular dynamics simulations of the DNA duplex d(CGCGAATTCGCG) in the presence of polycation chains are carried out to shed light on the specific atomic interaction that result in complex formation. The structures of complexes formed from DNA with polyethylenimine, which is considered one of the most promising DNA vector candidates, and a second polycation, poly-L-lysine, are compared. After an initial separation of approximately 50 A, the DNA and polycation come together and form a stable complex within 10 ns. The DNA does not undergo any major structural changes on complexation and remains in the B-form. In the formed complex, the charged amine groups of the polycation mainly interact with DNA phosphate groups, with polycation intrusion into the major and minor grooves dependent on the identity and charge state of the polycation. The ability of the polycation to effectively neutralize the charge of the DNA phosphate groups and the resulting influence on the DNA helix interaction are discussed.

Structural Analysis of DNA Interactions with Biogenic Polyamines and Cobalt(III)hexamine Studied by Fourier Transform Infrared and Capillary Electrophoresis

Biogenic polyamines, such as putrescine, spermidine, and spermine are small organic polycations involved in numerous diverse biological processes. These compounds play an important role in nucleic acid function due to their binding to DNA and RNA. It has been shown that biogenic polyamines cause DNA condensation and aggregation similar to that of inorganic cobalt(III)hexamine cation, which has the ability to induce DNA conformational changes. However, the nature of the polyamine.DNA binding at the molecular level is not clearly established and is the subject of much controversy. In the present study the effects of spermine, spermidine, putrescine, and cobalt(III)hexamine on the solution structure of calf-thymus DNA were investigated using affinity capillary electrophoresis, Fourier transform infrared, and circular dichroism spectroscopic methods. At low polycation concentrations, putrescine binds preferentially through the minor and major grooves of double strand DNA, whereas spermine, spermidine, and cobalt(III)hexamine bind to the major groove. At high polycation concentrations, putrescine interaction with the bases is weak, whereas strong base binding occurred for spermidine in the major and minor grooves of DNA duplex. However, major groove binding is preferred by spermine and cobalt(III)hexamine cations. Electrostatic attractions between polycation and the backbone phosphate group were also observed. No major alterations of B-DNA were observed for biogenic polyamines, whereas cobalt(III)hexamine induced a partial B --> A transition. DNA condensation was also observed for cobalt(III)hexamine cation, whereas organic polyamines induced duplex stabilization. The binding constants calculated for biogenic polyamines are K(Spm) = 2.3 x 10(5) M(-1), K(Spd) = 1.4 x 10(5) M(-1), and K(Put) = 1.02 x 10(5) M(-1). Two binding constants have been found for cobalt(III)hexamine with K(1) = 1.8 x 10(5) M(-1) and K(2) = 9.2 x 10(4) M(-1). The Hill coefficients indicate a positive cooperativity binding for biogenic polyamines and a negative cooperativity for cobalt(III)hexamine.

Simulation and modeling of nucleic acid structure, dynamics and interactions

In moving towards the simulation of larger nucleic acid assemblies over longer timescales that include more accurate representations of the environment, we are nearing the end of an era characterized by single nanosecond molecular dynamics simulation of nucleic acids. We are excited by the promise and predictability of the modeling methods, yet remain prudently cautious of sampling and force field limitations. Highlights include the accurate representation of subtle drug-DNA interactions, the detailed study of modified and unusual nucleic acid structures, insight into the influence of dynamics on the structure of DNA, and exploration of the interaction of solvent and ions with nucleic acids.