Dynamics Studies of a Malachite Green−RNA Complex Revealing the Origin of the Red-Shift and Energetic Contributions of Stacking Interactions

Photoinduced adsorption of oligonucleotides on polyvinyl chloride films containing malachite green derivative

DNA adsorption on the micrometer scale in a simple and cost-effective manner has received considerable interest. We prepared a film by casting an organic solvent containing polyvinyl chloride and a malachite green derivative, which can be photoionized to afford a cationic moiety for interaction with DNA. In this article, we report photoinduced oligonucleotide adsorption on a film that offers spatial and temporal control over oligonucleotide adsorption. Fluorescence microscopy was used to observe the oligonucleotide adsorption. Oligonucleotides of various sequences and lengths were also examined. UV irradiation using a photomask having 100 μm-diameter holes promoted the oligonucleotide adsorption on the film, whereas there was hardly any oligonucleotide adsorbed on the non-irradiated area. We found that the nucleobase contributed to the adsorption and part of the anchor in the oligonucleotide chain was responsible for the adsorption on the film.

Sampling Performance of Multiple Independent Molecular Dynamics Simulations of an RNA Aptamer

Using multiple independent simulations instead of one long simulation has been shown to improve the sampling performance attained with the molecular dynamics (MD) simulation method. However, it is generally not known how long each independent simulation should be, how many independent simulations should be used, or to what extent either of these factors affects the overall sampling performance achieved for a given system. The goal of the present study was to assess the sampling performance of multiple independent MD simulations, where each independent simulation begins from a different initial molecular conformation. For this purpose, we used an RNA aptamer that is 25 nucleotides long as a case study. The initial conformations of the aptamer are derived from six de novo predicted 3D structures. Each of the six de novo predicted structures is energy minimized in solution and equilibrated with MD simulations at high temperature. Ten conformations from these six high-temperature equilibration runs are selected as initial conformations for further simulations at ambient temperature. In total, we conducted 60 independent MD simulations, each with a duration of 100 ns, to study the conformation and dynamics of the aptamer. For each group of 10 independent simulations that originated from a particular de novo predicted structure, we evaluated the potential energy distribution of the RNA and used recurrence quantification analysis to examine the sampling of RNA conformational transitions. To assess the impact of starting from different de novo predicted structures, we computed the density of structure projection on principal components to compare the regions sampled by the different groups of ten independent simulations. The recurrence rate and dependence of initial conformation among the groups were also compared. We stress the necessity of using different initial configurations as simulation starting points by showing long simulations from different initial structures suffer from being trapped in different states. Finally, we summarized the sampling efficiency for the complete set of 60 independent simulations and determined regions of under-sampling on the potential energy landscape. The results suggest that conducting multiple independent simulations using a diverse set of de novo predicted structures is a promising approach to achieve sufficient sampling. This approach avoids undesirable outcomes, such as the problem of the RNA aptamer being trapped in a local minimum. For others wishing to conduct multiple independent simulations, the analysis protocol presented in this study is a guide for examining overall sampling and determining if more simulations are necessary for sufficient sampling.

Thioflavin T specifically brightening “Guanine Island” in duplex-DNA: a novel fluorescent probe for single-nucleotide mutation

Here, we found that Thioflavin T (ThT) could specifically bind with a G-GGG unit (named as “Guanine Island”) in double stranded DNA (ds-DNA).

Optimization of Experimental Parameters to Explore Small-Ligand/Aptamer Interactions through Use of (1) H NMR Spectroscopy and Molecular Modeling

Aptamers constitute an emerging class of molecules designed and selected to recognize any given target that ranges from small compounds to large biomolecules, and even cells. However, the underlying physicochemical principles that govern the ligand-binding process still have to be clarified. A major issue when dealing with short oligonucleotides is their intrinsic flexibility that renders their active conformation highly sensitive to experimental conditions. To overcome this problem and determine the best experimental parameters, an approach based on the design-of-experiments methodology has been developed. Here, the focus is on DNA aptamers that possess high specificity and affinity for small molecules, L-tyrosinamide, and adenosine monophosphate. Factors such as buffer, pH value, ionic strength, Mg(2+) -ion concentration, and ligand/aptamer ratio have been considered to find the optimal experimental conditions. It was then possible to gain new insight into the conformational features of the two ligands by using ligand-observed NMR spectroscopic techniques and molecular mechanics.

Thermodynamics and kinetics of adaptive binding in the malachite green RNA aptamer

Adaptive binding, the ability of molecules to fold themselves around the structure of a ligand and thereby incorporating it into their three-dimensional fold, is a key feature of most RNA aptamers. The malachite green aptamer (MGA) has been shown to bind several closely related triphenyl dyes with planar and nonplanar structures in this manner. Competitive binding studies using isothermal titration calorimetry and stopped flow kinetics have been conducted with the aim of understanding the adaptive nature of RNA-ligand interaction. The results of these studies reveal that binding of one ligand can reduce the ability of the aptamer pocket to adapt to another ligand, even if this second ligand has a significantly higher affinity to the free aptamer. A similar effect is observed in the presence of Mg(2+) ions which stabilize the binding pocket in a more ligand bound-like conformation.

Entropy and Mg2+ control ligand affinity and specificity in the malachite green binding RNA aptamer

The binding of small molecule targets by RNA aptamers provides an excellent model to study the versatility of RNA function. The malachite green aptamer binds and recognizes its ligand via stacking and electrostatic interactions. The binding of the aptamer to its original selection target and three related molecules was determined by isothermal titration calorimetry, equilibrium dialysis, and fluorescence titration. The results reveal that the entropy of complex formation plays a large role in determining binding affinity and ligand specificity. These data combined with previous structural studies show that metal ions are required to stabilize the complexes with non-native ligands whereas the complex with the original selection target is stable at low salt and in the absence of divalent metal ions.

Computational and experimental analyses converge to reveal a coherent yet malleable aptamer structure that controls chemical reactivity

As short nucleic acids, aptamers in solution are believed to be structurally flexible. Consistent with this view, most aptamers examined for this property have been shown to bind their target molecules by mechanisms that can be described as "induced fit". But, it is not known to what extent this structural flexibility affects the integrity of the target-aptamer interaction. Using the malachite green aptamer (MGA) as a model system, we show that the MGA can protect its bound target, malachite green (MG), from oxidation over several days. Protection is reversed by an oligonucleotide complementary to the MGA binding pocket. Computational cavity analysis of the MGA-MG structure predicted that MG oxidation is protected because a molecule as small as an OH(-) is sterically excluded from the C1 position of the bound MG. These results suggest that, while the MGA-MG interface is sufficiently coherent to prevent OH(-) penetration, the bases involved in the interaction are sufficiently mobile that they can exchange out of the MG binding interface to hybridize with a complementary oligonucleotide. The computational predictions were confirmed experimentally using variants of the MGA with single base changes in the binding pocket. This work demonstrates the successful application of molecular dynamics simulations and cavity analysis in determining the effects of sequence variations on the structure of a small single-stranded nucleic acid. It also shows that a nucleic acid aptamer can control access to specific chemical groups on its target, which suggests that aptamers might be applied for selectively protecting small molecules from modification.

Theoretical studies on the interaction of modified pyrimidines and purines with purine riboswitch

Recent experimental study [S.D. Gilbert, S.J. Mediatore, R.T. Batey, Modified pyrimidine specifically bind the purine riboswitch, J. Am. Chem. Soc. 128 (2006) 14214-14215] demonstrated that the purine riboswitch could specifically bind some ligands other than purines such as amino-pyrimidines, and the authors proposed that the five-membered ring of purine was not required for recognition. To get insight into the interaction details, we used molecular docking method to investigate the interactions of a mutant form of guanine riboswitch with a series of amino-purines, amino-pyrimidines and imidazole derivatives, and employed molecular simulation method to study the dynamic behavior of the selected complexes. The calculation results reveal that (1) all the amino-purines and amino-pyrimidines bind in a same cavity composed of four nucleobases including U22, U47, U51 and U74, which is consistent with the experimental results, while the two imidazole derivatives adopt other binding modes; (2) the purines are engulfed within three-way junction motifs, but most pyrimidines only form two-way junctions with the riboswitch; (3) the number and position of amino substituents could seriously affect the binding of pyrimidines. As riboswitches are potentially excellent candidates for antibiotic therapeutics, these findings may be useful for understanding the range of compounds that riboswitch can specifically recognize.

Fluorinated analogs of malachite green: synthesis and toxicity

A series of fluorinated analogs of malachite green (MG) have been synthesized and their toxicity to Saccharomyces cerevisiae and a human ovarian epithelial cell line examined. The toxicity profiles were found to be different for these two species. Two analogs, one with 2,4-difluoro substitution and the other with 2-fluoro substitution seem to be the most promising analogs because they showed the lowest toxicity to the human cells.

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.