Influence of sodium ions on the dynamics and structure of single-stranded DNA oligomers: a molecular dynamics study

Construction of Tough Hydrogel Cross-Linked via Ionic Interaction by Protection Effect of Hydrophobic Domains

Most hydrogels have poor mechanical properties, severely limiting their potential applications, and numerous approaches have been introduced to fabricate more robust and durable examples. However, these systems consist of nonbiodegradable polymers which limit their application in tissue engineering. Herein, we focus on the fabrication and investigate the influence of hydrophobic segments on ionic cross-linking properties for the construction of a tough, biodegradable hydrogel. A biodegradable, poly(γ-glutamic acid) polymer conjugated with a hydrophobic amino acid, l-phenylalanine ethyl ester (Phe), together with an ionic cross-linking group, alendronic acid (Aln) resulting in γ-PGA-Aln-Phe, was initially synthesized. Rheological assessments through time sweep oscillation testing revealed that the presence of hydrophobic domains accelerated gelation. Comparing gels with and without hydrophobic domains, the compressive strength of γ-PGA-Aln-Phe was found to be six times higher and exhibited longer stability properties in ethylenediaminetetraacetic acid solution, lasting for up to a month. Significantly, the contribution of the hydrophobic domains to the mechanical strength and stability of ionic cross-linking properties of the gel was found to be the dominant factor for the fabrication of a tough hydrogel. As a result, this study provides a new strategy for mechanical enhancement and preserves ionic cross-linked sites by the addition of hydrophobic domains. The development of tough, biodegradable hydrogels reported herein will open up new possibilities for applications in the field of biomaterials.

Single Stranded DNA Immune Modulators with Unmethylated CpG Motifs: Structure and Molecular Recognition by Toll-Like Receptor 9

Single stranded microbial DNA fragments with unmethylated deoxycytidylyldeoxyguanosine dinucleotide (CpG) motifs are interpreted as danger signals by the innate immune system via recognition by the Toll-like Receptor 9 (TLR9). Their synthetic analogues, Oligodeoxynucleotides (ODN) comprise a promising class of immune modulators with potential applications in the treatment of multiple diseases, such as cancer, autoimmune diseases or allergy. ODN molecules contain a core hexamer sequence, which is species specific consisting of GACGTT and AACGT for mouse and GTCGTT in humans. Assessment of structural features of different type of ODNs is highly challenging. NMR spectroscopic insights were gained for a short, single CpG motif containing ODN 1668. The structural basis of ODN recognition by TLR9 recently started to unravel as crystal structures of TLR9 orthologues in complex with ODN 1668 were solved. Systematic investigations of ODN sequences revealed that ODNs with a single CpG motif are capable of activating mouse TLR9, but two closely positioned CpG motifs are necessary for activation of human TLR9. Furthermore, longer ODNs with TCC and TCG sequences at the 5' end were shown to activate TLR9 with higher efficiency. It was revealed that 5'-xCx motif containing short ODNs (sODN) are able to augment the immune response of short, single CpG containing ODNs, which are incapable of activating of TLR9 alone. All these observations pointed to the existence of a second binding site on TLR9, which was characterized in crystal structures that delivered further insights of the nucleic acid recognition of the innate immune system by TLR9.

Structural Properties of Small Single-Stranded Circular Nucleic Acids

One strategy to avoid rapid degradation of small nucleic acids in biomedical applications is to covalently link their 3'- and 5'-ends, turning them into circular nucleic acids (circNAs). Here, we examine structural properties of flexible non-base-paired circNAs, containing 6-48 nucleotides, in aqueous solution, using microsecond long molecular dynamics simulations. Analyses of conformational ensembles of circular DNA (circDNA) and RNA (circRNA) molecules of different lengths and sequences reveal how their structures and dynamics are affected by the constraints of their geometries. The circDNAs are more bent and flexible than circRNAs, with distinctly different arrangements of phosphate backbones and bases. Small circNAs can sequester counterions in conformations that resemble crown ethers for the smallest (6-8 nucleotide long) molecules examined, in contrast to their linear counterparts. At millimolar concentrations (7.9 mM), circNA molecules were observed to aggregate, adopting linear chain shapes at physiological ionic strengths.

Dynamics of physiologically relevant noncanonical DNA structures: an overview from experimental and theoretical studies

DNA is a complex molecule with phenomenal inherent plasticity and the ability to form different hydrogen bonding patterns of varying stabilities. These properties enable DNA to attain a variety of structural and conformational polymorphic forms. Structurally, DNA can exist in single-stranded form or as higher-order structures, which include the canonical double helix as well as the noncanonical duplex, triplex and quadruplex species. Each of these structural forms in turn encompasses an ensemble of dynamically heterogeneous conformers depending on the sequence composition and environmental context. In vivo, the widely populated canonical B-DNA attains these noncanonical polymorphs during important cellular processes. While several investigations have focused on the structure of these noncanonical DNA, studying their dynamics has remained nontrivial. Here, we outline findings from some recent advanced experimental and molecular simulation techniques that have significantly contributed toward understanding the complex dynamics of physiologically relevant noncanonical forms of DNA.

DNA quadruplexes as molecular scaffolds for controlled assembly of fluorogens with aggregation-induced emission

Aggregation-induced emission (AIE) can be generated due to the restriction of intramolecular motions. The controllable assembly of fluorogens with AIE properties (AIEgens) is able to provide a new opportunity for precise manipulation of fluorescent signal transduction. Here, a tetrapod DNA quadruplex (TP-G4) was designed as a molecular scaffold for assembly and precise modulation of light emission of an oligonucleotide-grafted fluorogen with aggregation-induced emission (Oligo-AIEgen). The Oligo-AIEgen was synthesized by attaching the AIEgen to the 3'-terminus of the oligonucleotide through a dibenzylcyclooctyne mediated coupling reaction. The AIEgen emitted no detectable fluorescence in the context of a double-stranded structure. When hybridized to the parallel-stranded TP-G4, several AIEgens were located in close proximity to generate fluorescence. The fluorescence intensity has been precisely regulated by manipulation of the spacer length between the core structure of the scaffold and AIEgen, as well as by altering the quartet number of the G-quadruplex. Similar control of fluorescence was also demonstrated using tetramolecular and bimolecular i-motif quadruplex structures as the scaffolds. These scaffolds provide a proof of concept on the manipulation of molecular interactions, which forms a universal molecular tool for the design of new biosensing strategies.

Surface-Mediated DNA Hybridization: Effects of DNA Conformation, Surface Chemistry, and Electrostatics

Single-molecule Förster Resonance Energy Transfer (FRET) was used to study the dynamic association of mobile donor-labeled ssDNA oligonucleotides ("target") with covalently immobilized complementary acceptor-labeled ssDNA oligonucleotides ("probe"). While probe-target association events were resolved for all experiments, such FRET events were far more likely to occur in systems with complementarity and on hydrophobic, as compared to hydrophilic, surfaces. The distribution of donor-acceptor association-time intervals did not exhibit simple first-order kinetics, and when decomposed into a superposition of first-order processes, only a small fraction of events corresponded to a long-lived state that was presumed to represent true DNA hybridization, while the majority of association events were transient, representing nonspecific associations or partial hybridization. The structure of the DNA target and probe affected both the stability of the hybridized state, as well as the likelihood that an association between the two led to hybridization. In particular, the likelihood of hybridization decreased for longer target strands and for targets with stem-loop secondary structure. The presence of oligonucleotide secondary structure reduced the stability of hybridization, while greater complementarity increased stability of the hybridized state. Interestingly, increased ionic strength (i.e., greater electrostatic screening) increased the probability of hybridization but did not influence the lifetime of the hybridized state. Combined, these observations provide a nuanced view of surface-mediated DNA hybridization, where various factors independently influence the probability and stability of hybridization.

Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations

Given the ubiquitous nature of protein-DNA interactions, it is important to understand the interaction thermodynamics of individual amino acid side chains for DNA. One way to assess these preferences is to perform molecular dynamics (MD) simulations. Here we report MD simulations of 20 amino acid side chain analogs interacting simultaneously with both a 70-base-pair double-stranded DNA and with a 70-nucleotide single-stranded DNA. The relative preferences of the amino acid side chains for dsDNA and ssDNA match well with values deduced from crystallographic analyses of protein-DNA complexes. The estimated apparent free energies of interaction for ssDNA, on the other hand, correlate well with previous simulation values reported for interactions with isolated nucleobases, and with experimental values reported for interactions with guanosine. Comparisons of the interactions with dsDNA and ssDNA indicate that, with the exception of the positively charged side chains, all types of amino acid side chain interact more favorably with ssDNA, with intercalation of aromatic and aliphatic side chains being especially notable. Analysis of the data on a base-by-base basis indicates that positively charged side chains, as well as sodium ions, preferentially bind to cytosine in ssDNA, and that negatively charged side chains, and chloride ions, preferentially bind to guanine in ssDNA. These latter observations provide a novel explanation for the lower salt dependence of DNA duplex stability in GC-rich sequences relative to AT-rich sequences.

Exploring ion induced folding of a single-stranded DNA oligomer from molecular simulation studies

One crucial issue in DNA hydration is the effect of salts on its conformational features. This has relevance in biology as cations present in the cellular environment shield the negative charges on the DNA backbone, thereby reducing the repulsive force between them. By screening the negative charges along the backbone, cations stabilize the folded structure of DNA. To study the effect of the added salt on single-stranded DNA (ss-DNA) conformations, we have performed room temperature molecular dynamics simulations of an aqueous solution containing the ss-DNA dodecamer with the 5'-CGCGAATTCGCG-3' sequence in the presence of 0.2, 0.5, and 0.8 M NaCl. Our calculations reveal that in the presence of the salt, the DNA molecule forms more collapsed coil-like conformations due to the screening of negative charges along the backbone. Additionally, we demonstrated that the formation of an octahedral inner-sphere complex by the strongly bound ion plays an important role in the stabilization of such folded conformation of DNA. Importantly, it is found that ion-DNA interactions can also explain the formation of non-sequential base stackings with longer lifetimes. Such non-sequential base stackings further stabilize the collapsed coil-like folded form of the DNA oligomer.

Effect of C5-Methylation of Cytosine on the UV-Induced Reactivity of Duplex DNA: Conformational and Electronic Factors

C5-methylation of cytosines is strongly correlated with UV-induced mutations detected in skin cancers. Mutational hot-spots appearing at TCG sites are due to the formation of pyrimidine cyclobutane dimers (CPDs). The present study, performed for the model DNA duplex (TCGTA)3·(TACGA)3 and the constitutive single strands, examines the factors underlying the effect of C5-methylation on pyrimidine dimerization at TCG sites. This effect is quantified for the first time by quantum yields ϕ. They were determined following irradiation at 255, 267, and 282 nm and subsequent photoproduct analysis using HPLC coupled to mass spectrometry. C5-methylation leads to an increase of the CPD quantum yield up to 80% with concomitant decrease of that of pyrimidine(6-4) pyrimidone adducts (64PPs) by at least a factor of 3. The obtained ϕ values cannot be explained only by the change of the cytosine absorption spectrum upon C5-methylation. The conformational and electronic factors that may affect the dimerization reaction are discussed in light of results obtained by fluorescence spectroscopy, molecular dynamics simulations, and quantum mechanical calculations. Thus, it appears that the presence of an extra methyl on cytosine affects the sugar puckering, thereby enhancing conformations of the TC step that are prone to CPD formation but less favorable to 64PPs. In addition, C5-methylation diminishes the amplitude of conformational motions in duplexes; in the resulting stiffer structure, ππ* excitations may be transferred from initially populated exciton states to reactive pyrimidines giving rise to CPDs.

Effect of temperature on the low-frequency vibrational spectrum and relative structuring of hydration water around a single-stranded DNA

Molecular dynamics simulations of the single-stranded DNA oligomer (5'-CGCGAAT TCGCG-3') in aqueous solution have been carried out at different temperatures between 160 K and 300 K. The effects of temperature on the low-frequency vibrational spectrum and local structural arrangements of water molecules hydrating the DNA strand have been explored in detail. The low-frequency density of states distributions reveal that increasingly trapped transverse water motions play a dominant role in controlling the band corresponding to O⋯O⋯O bending or transverse oscillations of hydration water at supercooled temperatures. In addition, presence of a broad band around 260 (±20) cm(-1) under supercooled conditions indicates transformation from high density liquid-like structuring of hydration water at higher temperatures to that of a low density liquid at lower temperatures. It is found that long-range correlations between the supercooled hydration water molecules arise due to such local structural transition around the DNA oligomer.

Molecular dynamics simulation of a single-stranded DNA with heterogeneous distribution of nucleobases in aqueous medium

The DNA metabolic processes often involve single-stranded DNA (ss-DNA) molecules as important intermediates. In the absence of base complementarity, ss-DNAs are more flexible and interact strongly with water in aqueous media. Ss-DNA-water interactions are expected to control the conformational flexibility of the DNA strand, which in turn should influence the properties of the surrounding water molecules. We have performed room temperature molecular dynamics simulation of an aqueous solution containing the ss-DNA dodecamer, 5'-CGCGAATTCGCG-3'. The conformational flexibility of the DNA strand and the microscopic structure and ordering of water molecules around it have been explored. The simulation reveals transformation of the initial base-stacked form of the ss-DNA to a fluctuating collapsed coil-like conformation with the formation of a few non-sequentially stacked base pairs. A preliminary analysis shows further collapse of the DNA conformation in presence of additional salt (NaCl) due to screening of negative charges along the backbone by excess cations. Additionally, higher packing of water molecules within a short distance from the DNA strand is found to be associated with realignment of water molecules by breaking their regular tetrahedral ordering.

UV activation of polymeric high aspect ratio microstructures: ramifications in antibody surface loading for circulating tumor cell selection

The need to activate thermoplastic surfaces using robust and efficient methods has been driven by the fact that replication techniques can be used to produce microfluidic devices in a high production mode and at low cost, making polymer microfluidics invaluable for in vitro diagnostics, such as circulating tumor cell (CTC) analysis, where device disposability is critical to mitigate artifacts associated with sample carryover. Modifying the surface chemistry of thermoplastic devices through activation techniques can be used to increase the wettability of the surface or to produce functional scaffolds to allow for the covalent attachment of biologics, such as antibodies for CTC recognition. Extensive surface characterization tools were used to investigate UV activation of various surfaces to produce uniform and high surface coverage of functional groups, such as carboxylic acids in microchannels of different aspect ratios. We found that the efficiency of the UV activation process is highly dependent on the microchannel aspect ratio and the identity of the thermoplastic substrate. Colorimetric assays and fluorescence imaging of UV-activated microchannels following EDC/NHS coupling of Cy3-labeled oligonucleotides indicated that UV-activation of a PMMA microchannel with an aspect ratio of ~3 was significantly less efficient toward the bottom of the channel compared to the upper sections. This effect was a consequence of the bulk polymer's damping of the modifying UV radiation due to absorption artifacts. In contrast, this effect was less pronounced for COC. Moreover, we observed that after thermal fusion bonding of the device's cover plate to the substrate, many of the generated functional groups buried into the bulk rendering them inaccessible. The propensity of this surface reorganization was found to be higher for PMMA compared to COC. As an example of the effects of material and microchannel aspect ratios on device functionality, thermoplastic devices for the selection of CTCs from whole blood were evaluated, which required the immobilization of monoclonal antibodies to channel walls. From our results, we concluded the CTC yield and purity of isolated CTCs were dependent on the substrate material with COC producing the highest clinical yields for CTCs as well as better purities compared to PMMA.

Effects of molecular size and surface hydrophobicity on oligonucleotide interfacial dynamics

Single-molecule total internal reflection fluorescence microscopy was used to observe the dynamic behavior of polycytosine single-stranded DNA (ssDNA) (1-50 nucleotides long) at the interface between aqueous solution and hydrophilic (oligoethylene glycol-modified fused silica, OEG) and hydrophobic (octadecyltriethoxysilane-modified fused silica, OTES) solid surfaces. High throughput molecular tracking was used to determine >75,000 molecular trajectories for each molecular length, which were then used to calculate surface residence time and squared displacement (i.e., "step-size") distributions. On hydrophilic OEG surfaces, the surface residence time increased systematically with ssDNA chain length, as expected due to increasing molecule-surface interactions. Interestingly, the residence time decreased with increasing ssDNA length on the hydrophobic OTES surface, particularly for longer chains. Similarly, the interfacial mobility of polynucleotides slowed with increasing chain length on OEG, but became faster on OTES. On OTES surfaces, the rates associated with desorption and surface diffusion exhibited the distinctive anomalous temperature dependence that is characteristic of hydrophobic interactions for short-chain species but not for longer chains. These combined observations suggest that long oligonucleotides adopt conformations minimizing hydrophobic interactions, e.g., by internal sequestration of hydrophobic nucleobases.

Mass spectrometric study on sodium ion induced central nucleotide deletion in the gas phase

We report a mass spectrometric study on sodium ion induced central nucleotide deletion from protonated oligonucleotides (ONTs) and the concurrent recombination of the terminal nucleotides. To shed some light on the mechanism behind this intriguing fragmentation channel, we have studied the metastable decay of a number of different protonated hexameric and octameric oligonucleotides with 0-6 and 0-8 of their exchangeable protons replaced with sodium ions, respectively. In selected cases, we have also studied the further fragmentation of the parent ions after initial base loss. Our findings are concurrent with a reaction mechanism where the initial step is the elimination of a protonated, high proton affinity (PA) base from the center of the ONTs. This is followed by an elimination of a (next neighbour) nucleotide that contains a second high PA base and the concurrent recombination of the terminal nucleotides. To our knowledge, such central nucleotide deletion in the gas phase has only been reported in one previous study (Flosadóttir et al., J. Am. Soc. Mass Spectrom 20:689-696, 2009), and this is the first systematic approach to understand the mechanism behind this channel.

Effect of oligonucleotide length on the assembly of DNA materials: molecular dynamics simulations of layer-by-layer DNA films

DNA strand length has been found to be an important factor in many DNA-based nanoscale systems. Here, we apply molecular dynamics simulations in a synergistic effort with layer-by-layer experimental data to understand the effect of DNA strand length on the assembly of DNA films. The results indicate that short (less than 10 bases) and long (more than 30 bases) single-stranded DNAs do not exhibit optimal film growth, and this can be associated with the limited accessibility of the bases on the surface due to formation of self-protected interactions that prevent efficient hybridization. Interestingly, the presence of a duplex attached to a single strand significantly alters the persistence length of the polyT strands. Our study suggests that restrained polyT, compared to labile suspensions of free polyT, are more capable of hybridization and hence DNA-based assembly.

Simulating DNA hybridization on an amine-functionalized silicon substrate

The ability of DNA biosensors to capture oligonucleotide molecules in solution is of great importance in genetics, medical diagnostics, and drug discovery. The DNA hybridization event in which the probe, which is usually a single-stranded DNA segment covalently immobilized on a functionalized surface via a cross-linker molecule, recognizes the complementary target and forms a stable duplex structure is the basis of highly specific biorecognizing devices. The rate of hybridization depends on the solvent, length of the strands, complexity of the system, and other factors and could be considerably altered by the type of attachment and by the density of the probe on the substrates. Recent experimental investigations have shown that some probes can hybridize directly from bulk solutions. In this computational study, we provide a model for the behavior of these systems choosing cross-linker, probe, and target on the basis of experimental data. MD simulations of the single-stranded DNA fragment 5'-d(TGGC)-3' attached to an allylamine-functionalized Si(111) surface through an oxanine cross-linker in aqueous solution containing the complementary sequence, i.e. 5'-d(CGCCA)-3', are presented. A possible probe-target capture mechanism obtained using explicit solvent and state-of-the-art classical molecular dynamics simulation protocols is described. The hybridization process of the tethered DNA single strand, the intermediate structures appeared during the formation of the double helix, their internal dynamics and their behavior with respect to the substrate are characterized in detail.

Thymine dimerization in DNA model systems: cyclobutane photolesion is predominantly formed via the singlet channel

UV-induced formation of cylcobutane pyrimidine dimers (CPD) in all thymine DNA models have been studied by femtosecond IR spectroscopy. CPDs are shown to form within approximately 1 ps during the decay of the initially excited (1)pi pi * state. The quantum yields phi(D)(ps) determined after the (1)pi pi * decay equal the final yield phi(D)(cw). This gives evidence for a predominance of the singlet channel in CPD formation.

DNA excited-state dynamics: from single bases to the double helix

Ultraviolet light is strongly absorbed by DNA, producing excited electronic states that sometimes initiate damaging photochemical reactions. Fully mapping the reactive and nonreactive decay pathways available to excited electronic states in DNA is a decades-old quest. Progress toward this goal has accelerated rapidly in recent years, in large measure because of ultrafast laser experiments. Here we review recent discoveries and controversies concerning the nature and dynamics of excited states in DNA model systems in solution. Nonradiative decay by single, solvated nucleotides occurs primarily on the subpicosecond timescale. Surprisingly, excess electronic energy relaxes one or two orders of magnitude more slowly in DNA oligo- and polynucleotides. Highly efficient nonradiative decay pathways guarantee that most excited states do not lead to deleterious reactions but instead relax back to the electronic ground state. Understanding how the spatial organization of the bases controls the relaxation of excess electronic energy in the double helix and in alternative structures is currently one of the most exciting challenges in the field.

UV excitation of single DNA and RNA strands produces high yields of exciplex states between two stacked bases

Excited electronic states created by UV excitation of the diribonucleoside monophosphates ApA, ApG, ApC, ApU, and CpG were studied by the femtosecond transient-absorption technique. Bleach recovery signals recorded at 252 nm show that long-lived excited states are formed in all five dinucleosides. The lifetimes of these states exceed those measured in equimolar mixtures of the constituent mononucleotides by one to two orders of magnitude, indicating that electronic coupling between proximal nucleobases dramatically slows the relaxation of excess electronic energy. The decay rates of the long-lived states decrease with increasing energy of the charge-transfer state produced by transferring an electron from one base to another. The charge-transfer character of the long-lived states revealed by this analysis supports their assignment to excimer or exciplex states. Identical bleach recovery signals were seen for ApA, (A)(4), and poly(A) at delay times >10 ps after photoexcitation. This indicates that excited states localized on a stack of just two bases are the common trap states independent of the number of stacked nucleotides. The fraction of initial excitations that decay to long-lived exciplex states is approximately equal to the fraction of stacked bases determined by NMR measurements. This supports a model in which excitations associated with two stacked bases decay to exciplex states, whereas excitations in unstacked bases decay via ultrafast internal conversion. These results establish the importance of charge transfer-quenching pathways for UV-irradiated RNA and DNA in room-temperature solution.

Electronic excitations and spectra in single-stranded DNA

Using density functional theory and molecular dynamics simulations, we show that delocalized states extending over three bases can be directly excited in single-stranded poly(A) DNA. The results are in semiquantitative agreement with recent experimental results for the delocalization length of these states in single- and double-stranded DNA. The structures used in these molecular dynamics calculations are validated by comparing calculated circular dichroic spectra for d(A)2 and d(A)4 with experiment. These spectra, which arise from highly stacked structures, are in good agreement with experiment, suggesting that the short delocalization in ssDNA arises in spite of strong stacking.

A doorway state leads to photostability or triplet photodamage in thymine DNA

Ultraviolet irradiation of DNA produces electronic excited states that predominantly eliminate the excitation energy by returning to the ground state (photostability) or following minor pathways into mutagenic photoproducts (photodamage). The cyclobutane pyrimidine dimer (CPD) formed from photodimerization of thymines in DNA is the most common form of photodamage. The underlying molecular processes governing photostability and photodamage of thymine-constituted DNA remain unclear. Here, a combined femtosecond broadband time-resolved fluorescence and transient absorption spectroscopies were employed to study a monomer thymidine and a single-stranded thymine oligonucleotide. We show that the protecting deactivation of a thymine multimer is due to an ultrafast single-base localized stepwise mechanism where the initial excited state decays via a doorway state to the ground state or proceeds via the doorway state to a triplet state identified as a major precursor for CPD photodamage. These results provide new mechanistic characterization of and a dynamic link between the photoexcitation of DNA and DNA photostability and photodamage.

Structure and dynamics of poly(T) single-strand DNA: implications toward CPD formation

The formation of cyclobutane pyrimidine dimers between adjacent thymines by UV radiation is thought to be the first event in a cascade leading to skin cancer. Recent studies showed that thymine dimers are fully formed within 1 ps of UV irradiation, suggesting that the conformation at the moment of excitation is the determining factor in whether a given base pair dimerizes. MD simulations on the 50 ns time scale are used to study the populations of reactive conformers that exist at any given time in T18 single-strand DNA. Trajectory analysis shows that only a small percentage of the conformations fulfill distance and dihedral requirements for thymine dimerization, in line with the experimentally observed quantum yield of 3%. Plots of the pairwise interactions in the structures predict hot spots of DNA damage where dimerization in the ssT18 is predicted to be most favored. The importance of hairpin formation by intra-strand base pairing for distinguishing reactive and unreactive base pairs is discussed in detail. The data presented thus explain the structural origin of the results from the ultrafast studies of thymine dimer formation.

The polyamines spermidine and spermine retard the platination rate of single-stranded DNA oligomers and plasmids

The presence of polyamines in living cells is crucial for survival. Due to their high net charge at physiological pH, polyamines effectively charge neutralize the phosphodiester backbone of DNA in an interaction that also may protect the DNA from external damage. We here present a study illustrating the influence of spermidine and spermine on the platination reactions of the model oligonucleotides d(T(6)GT(6)), d(T(12)GT(12)), and d(T(24)GT(24)), and the pUC18 DNA plasmid. The aquated forms of the anticancer active compounds cisplatin (cis-[Pt(NH(3))(2)Cl(2)]) and the major Pt(II) metabolite of JM216 (cis-[PtCl(2)(NH(3))(c-C(6)H(11)NH(2))], JM118) were used as platination reagents. The study shows that the kinetics for formation of the coordinative Pt-DNA adduct are strongly influenced by the presence of sub-millimolar polyamine concentrations. At polyamine concentrations in the muM-range, the reactions remain salt-dependent. In contrast, platination of pUC18 is effectively prevented at mM concentrations of both spermidine and spermine with the latter as the more potent inhibitor. The results suggest that variations of intracellular polyamine concentrations may have a profound influence on the efficacy by which cationically charged reagents interfere with DNA function in vivo by modulation of the preassociation conditions.

Thymine dimerization in DNA is an ultrafast photoreaction

Femtosecond time-resolved infrared spectroscopy was used to study the formation of cyclobutane dimers in the all-thymine oligodeoxynucleotide (dT)18 by ultraviolet light at 272 nanometers. The appearance of marker bands in the time-resolved spectra indicates that the dimers are fully formed approximately 1 picosecond after ultraviolet excitation. The ultrafast appearance of this mutagenic photolesion points to an excited-state reaction that is approximately barrierless for bases that are properly oriented at the instant of light absorption. The low quantum yield of this photoreaction is proposed to result from infrequent conformational states in the unexcited polymer, revealing a strong link between conformation before light absorption and photodamage.

Steered Molecular Dynamics Studies of the Potential of Mean Force of a Na+ or K+ Ion in a Cyclic Peptide Nanotube

Potential of mean force (PMF) profiles of a single Na+ or K+ ion passing through a cyclic peptide nanotube, cyclo[-(D-Ala-Glu-D-Ala-Gln)2-], in water are calculated to provide insight into ion transport and to understand the conductance difference between these two ions. The PMF profiles are obtained by performing steered molecular dynamics (SMD) simulations that are based on the Jarzynski equality. The computed PMF profiles for both ions show barriers of around 2.4 kcal/mol at the channel entrances and exits and energy wells in the middle of the tube. The energy barriers, so-called dielectric energy barriers, arise due to the desolvation of water molecules when ions move across the nanotube, and the energy wells appear as a result of attractive interactions between the cations and negatively charged carbonyl oxygens on the backbone of the tube. We find more and deeper energy wells in the PMF profile for Na+ than for K+, which suggests that Na+ ions have a longer residence time inside the nanotube and that permeation of Na+ ions is reduced compared to K+ ions. Calculations of the radial distribution functions (RDF) between the ions and oxygens in the water molecules and in carbonyl groups on the tube and an investigation of the orientations of the carbonyl groups show that, in contrast with the dynamic carbonyl groups observed in the selectivity filter of the KcsA ion channel, the carbonyl groups in the cyclic peptide nanotube are relatively rigid, with only slight reorientation of the carbonyl groups as the cations pass through. The rigidity of the carbonyl groups in the cyclic peptide nanotube can be attributed to their role in hydrogen bonding, which is responsible for the tube structure. Comparison of the PMF profiles with the electrostatic energy profiles calculated from the Poisson-Boltzmann (PB) equation, a dielectric continuum model, reveals that the dielectric continuum model breaks down in the confined region within the tube that governs ion transport.

Single-strand conformation polymorphism (SSCP) of oligodeoxyribonucleotides: an insight into solution structural dynamics of DNAs provided by gel electrophoresis and molecular dynamics simulations

Studies on the solution structure dynamics of RNA/DNA are becoming crucially important. The phenomena of SSCP (single-strand conformation polymorphism), small RNA dynamics in a cell, and others can be related to the conformational changes of single-stranded (ss) RNAs/DNAs in solution. However, little is known about those dynamics. Only the intra-structural transition of ssDNAs in solution has been reported based on Watson-Crick (W-C) base-pairing. Here, we found a general feature of the SSCP phenomenon by studying the simpler molecules of ss-oligodeoxyribonucleotides. A single base substitution or a positional exchange of nucleotide in a highly homologous series of ss-dodecanucleotides led to a change in the mobility-in-gel. This was unexpected, since most of these nucleotides [such as d(A(11)G) or d(A(11)C)] have no possibility of forming W-C base-pairing. MD (molecular dynamics) experiments revealed differences in shape and size between the dynamic structures of these molecules which could affect their mobility-in-gel. In addition, a high correlation was observed between the electrophoretic mobility and the size-related parameters such as end-to-end distance obtained from MD simulations. Because the simulation was considerably shorter (nanosecond) than the experimental time-scale (second), the result must be considered conservatively; but it is nevertheless encouraging for utilizing MD simulation for structural analysis of oligonucleotides.

Orientational dynamics and dye-DNA interactions in a dye-labeled DNA aptamer

We report the picosecond and nanosecond timescale rotational dynamics of a dye-labeled DNA oligonucleotide or "aptamer" designed to bind specifically to immunoglobulin E. Rotational dynamics in combination with fluorescence lifetime measurements provide information about dye-DNA interactions. Comparison of Texas Red (TR), fluorescein, and tetramethylrhodamine (TAMRA)-labeled aptamers reveals surprising differences with significant implications for biophysical studies employing such conjugates. Time-resolved anisotropy studies demonstrate that the TR- and TAMRA-aptamer anisotropy decays are dominated by the overall rotation of the aptamer, whereas the fluorescein-aptamer anisotropy decay displays a subnanosecond rotational correlation time much shorter than that expected for the overall rotation of the aptamer. Docking and molecular dynamics simulations suggest that the low mobility of TR is a result of binding in the groove of the DNA helix. Additionally, associated anisotropy analysis of the TAMRA-aptamer reveals both quenched and unquenched states that experience significant coupling to the DNA motion. Therefore, quenching of TAMRA by guanosine must depend on the configuration of the dye bound to the DNA. The strong coupling of TR to the rotational dynamics of the DNA aptamer, together with the absence of quenching of its fluorescence by DNA, makes it a good probe of DNA orientational dynamics. The understanding of the nature of dye-DNA interactions provides the basis for the development of bioconjugates optimized for specific biophysical measurements and is important for the sensitivity of anisotropy-based DNA-protein interaction studies employing such conjugates.

Direct measurement of the association constant of HER2/neu antisense oligonucleotide to its target RNA sequence using a molecular beacon

A molecular beacon approach was developed to directly determine the association constant of RNA-DNA hybrid formation. The molecular beacon was composed of a 15-nt loop structure containing the antisense sequence that can hybridize with the AUG translational start site of the HER2/neu gene, which is overexpressed in a significant proportion of breast, ovarian, and lung tumors. The equilibrium association constant (Ka) of DNA binding to the RNA oligonucleotide was 6.4 +/- 0.14 x 10(7) M(-1) in the presence of 150 mM NaCl at 22 degrees C. The free energy change (AG) associated with RNA-DNA hybrid formation was -10.7 kcal/mole. The melting temperature (Tm) of RNA-DNA hybrid was 64.4 degrees C +/- 1 degree C in the presence of 150 mM NaCl. The RNA-DNA hybrid was more stable than the corresponding DNA-DNA duplex in 150 mM NaCl, as judged by both Ka and Tm data. We also determined the Ka, deltaG, and Tm values of RNA-DNA and DNA-DNA duplex formation in the presence of three monovalent cations, Li+, K+, and Cs+. The feasibility of this method was also investigated using a phosphorothioate molecular beacon. The information generated through this new approach for thermodynamic measurements might be useful for the design of oligonucleotides for antisense therapeutics.