Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Origin of DNA helical structure and its sequence dependence

Conformational analysis of DNA shows that the origin of the B-form double helix can be attributed in large part to the atomic charge pattern in the base pairs. The charge patterns favor specific helical stacking of the base pairs. Base pairs alone--without backbones--have a strong tendency to form helix, indicating that the backbones play a rather passive role in determining the basic helical structure of DNA. It is mainly the electrostatic interactions determined by the charge pattern on base pairs that stabilize a particular helical conformation. The charge pattern in the base pairs appears to be responsible for much of the sequence dependence of DNA conformation, rather than steric clashes.     

Analysis of Four-Way Junctions in RNA Structures

RNA secondary structures can be divided into helical regions composed of canonical Watson-Crick and related base pairs, as well as single-stranded regions such as hairpin loops, internal loops, and junctions. These elements function as building blocks in the design of diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze existing RNA four-way junctions in terms of base-pair interactions and 3D configurations. Specifically, we identify nine broad junction families according to coaxial stacking patterns and helical configurations. We find that helices within junctions tend to arrange in roughly parallel and perpendicular patterns and stabilize their conformations using common tertiary motifs such as coaxial stacking, loop-helix interaction, and helix packing interaction. Our analysis also reveals a number of highly conserved base-pair interaction patterns and novel tertiary motifs such as A-minor-coaxial stacking combinations and sarcin/ricin motif variants. Such analyses of RNA building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

De-intercalation of ethidium bromide and propidium iodine from DNA in the presence of caffeine

Caffeine (CAF) is capable of interacting directly with several genotoxic aromatic ligands by stacking aggregation. Formation of such hetero-complexes may diminish pharmacological activity of these ligands, which is often related to its direct interaction with DNA. To check these interactions we performed three independent series of spectroscopic titrations for each ligand (ethidium bromide, EB, and propidium iodine, PI) according to the following setup: DNA with ligand, ligand with CAF and DNA-ligand mixture with CAF. We analyzed DNA-ligand and ligand-CAF mixtures numerically using well known models: McGhee-von Hippel model for ligand-DNA interactions and thermodynamic-statistical model of mixed association of caffeine with aromatic ligands developed by Zdunek et al. (2000). Based on these models we calculated association constants and concentrations of mixture components using a novel method developed here. Results are in good agreement with parameters calculated in separate experiments and demonstrate de-intercalation of EB and PI molecules from DNA caused by CAF.     

Transmembrane domains in the functions of Fc receptors

In the present study, we use a novel method, PHDhtm, to predict the exact locations and extents of the transmembrane (TM) domains of multisubunit immunoglobulin Fc-receptors. Whereas most previous studies have used single residue hydrophobicity plots for characterizing of these domains, PHDhtm utilizes a system of neural networks and the evolutionary information contained in multiple alignments of related sequences to predict the above. Present PHDhtm application predicts TM domains of immunoglobulin Fc-receptors that in many cases differ significantly from those derived by using earlier methods. Comparisons of helical wheel projections of the presently derived TM domains from PHDhtm with those produced earlier reveal different hydrophobic moments as well as hydrophobic and hydrophilic surfaces. These differences probably alter the character of subunit association within the receptor complexes. This new algorithm can also be used for other membrane protein complexes and may advance both understanding the principles underlying such complexes formation and design of peptides that can interfere with such TM domain association so as to modulate specific cellular responses.     

The spatial configuration of ordered polynucleotide chains. I. Helix formation and base stacking.

A single virtual bond scheme set forth previously for the treatment of average properties of randomly coiling polynucleotides is here applied to the calculation of helical parameters which characterize a regularly repeating polynucleotide molecule. Only a fraction of the enormous number of conformationally feasible helixes fulfill the geometric criteria of vertical base stacking usually associated with ordered polynucleotide chains. Detailed examination of the nature and mode of base stacking feasible in a single helical backbone structure indicates that the handedness of a base stacking arrangement does not correlate either quantitatively or qualitatively with the handedness of the polymer backbone. A number of polynucleotide chains which exhibit lefthanded base stacking patterns in nmr and CD studies may, in fact, be righthanded helixes.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

Tertiary Motifs Revealed in Analyses of Higher-Order RNA Junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Cleavage of a four-way DNA junction by a restriction enzyme spanning the point of strand exchange.

The four-way DNA junction is believed to fold in the presence of metal ions into an X-shaped structure, in which there is pairwise coaxial stacking of helical arms. A restriction enzyme MboII has been used to probe this structure. A junction was constructed containing a recognition site for MboII in one helical arm, positioned such that stacking of arms would result in cleavage in a neighbouring arm. Strong cleavage was observed, at the sites expected on the basis of coaxial stacking. An additional cleavage was seen corresponding to the formation of an alternative stacking isomer, suggesting that the two isomeric forms are in dynamic equilibrium in solution.     

Fluorescent stains for quantification of DNA by confocal laser scanning microscopy in 3-D

Confocal microscopy requires the use of fluorophores to visualize structures of interest within a specimen. To perform reliable measurements of the intensity of fluorescence, the stain should be specific, penetrate well into tissue sections, and bind stoichiometrically. Furthermore, emission must be linear with respect to DNA content and brightness, and fluorescence should be stable. Confocal microscopy is used to determine DNA ploidy and to analyze texture of nuclei, which is accomplished in three dimensions, because nuclei can be measured within the original tissue context. For this purpose the sample must be stained with a DNA binding fluorophore with the properties described above. Stains with different properties have been developed for different applications. We review here the advantages and disadvantages of these different stains for analyzing DNA ploidy and nuclear texture using three-dimensional microscopy. We conclude that SYBR green I and TO-PRO-3 are the most suitable stains for this purpose at present.     

What drives the binding of minor groove-directed ligands to DNA hairpins?

Understanding the molecular basis of ligand-DNA-binding events, and its application to the rational design of novel drugs, requires knowledge of the structural features and forces that drive the corresponding recognition processes. Existing structural evidence on DNA complexation with classical minor groove-directed ligands and the corresponding studies of binding energetics have suggested that this type of binding can be described as a rigid-body association. In contrast, we show here that the binding-coupled conformational changes may be crucial for the interpretation of DNA (hairpin) association with a classical minor groove binder (netropsin). We found that, although the hairpin form is the only accessible state of ligand-free DNA, its association with the ligand may lead to its transition into a duplex conformation. It appears that formation of the fully ligated duplex from the ligand-free hairpin, occurring via two pathways, is enthalpically driven and accompanied by a significant contribution of the hydrophobic effect. Our thermodynamic and structure-based analysis, together with corresponding theoretical studies, shows that none of the predicted binding steps can be considered as a rigid-body association. In this light we anticipate our thermodynamic approach to be the basis of more sophisticated nucleic acid recognition mechanisms, which take into account the dynamic nature of both the nucleic acid and the ligand molecule.     

Novel N-(3-carboxyl-9-benzyl-carboline-1-yl)ethylamino acids: synthesis, anti-proliferation activity and two-step-course of intercalation with calf thymus DNA

To explore the intercalating mechanism of -carbolines, four novel N-(3-carboxyl-9-benzyl-carboline-1-yl)ethylamino acids [-phenylalanine (6a), -alanine (6b), -isoleucine (6c) and -glycine (6d)] were prepared here. Their in vitro anticancer activities were examined by their anti-proliferation for 5 human carcinoma cell lines. The average IC50s against 5 human carcinoma cell lines are 53.54 microM, 118.77 microM, 147.34 microM and greater than 291.63 microM for 6a, 6b, 6c and 6d, respectively. The DNA intercalating mechanism of 6a-d was approved by the comparison of the parameters and signals of UV, CD and fluorescence spectra of calf thymus DNA (CT DNA) alone and the CT DNA/6a-d system. Using fluorescence titration based kinetic analysis a two-step-course consisting of stacking and intercalating was described and the stacking was considered as the key step to the CT DNA intercalating mechanism of 6a-d. Using fluorescence titration based thermomechanical analysis, the stacking complexes of 6a-d with CT DNA were described to be formed spontaneously and to be stabilized predominantly by their hydrophobic interactions. The intercalation itself goes very fast and only has limited contribution to their anticancer activities.     

The thermodynamics of polyamide-DNA recognition: hairpin polyamide binding in the minor groove of duplex DNA

Crescent-shaped synthetic ligands containing aromatic amino acids have been designed for specific recognition of predetermined DNA sequences in the minor groove of DNA. Simple rules have been developed that relate the side-by-side pairings of Imidazole (Im) and Pyrrole (Py) amino acids to their predicted target DNA sequences. We report here thermodynamic characterization of the DNA-binding properties of the six-ring hairpin polyamide, ImImPy-gamma-PyPyPy-beta-Dp (where gamma = gamma-aminobutyric acid, beta = beta-alanine, and Dp = dimethylaminopropylamide). Our data reveal that, at 20 degrees C, this ligand binds with a relatively modest 1.8-fold preference for the designated match site, 5'-TGGTA-3', over the single base pair mismatch site, 5'-TGTTA-3'. By contrast, we find that the ligand exhibits a 102-fold greater affinity for its designated match site relative to the double base pair mismatch site, 5'-TATTA-3'. These results demonstrate that the energetic cost of binding to a double mismatch site is not necessarily equal to twice the energetic cost of binding to a single mismatch site. Our calorimetrically measured binding enthalpies and calculated entropy data at 20 degrees C reveal the ligand sequence specificity to be enthalpic in origin. We have compared the DNA-binding properties of ImImPy-gamma-PyPyPy-beta-Dp with the hairpin polyamide, ImPyPy-gamma-PyPyPy-beta-Dp (an Im --> Py "mutant"). Our data reveal that both ligands exhibit high affinities for their designated match sites, consistent with the Dervan pairing rules. Our data also reveal that, relative to their corresponding single mismatch sites, ImImPy-gamma-PyPyPy-beta-Dp is less selective than ImPyPy-gamma-PyPyPy-beta-Dp for its designated match site. This result suggests, at least in this case, that enhanced binding affinity can be accompanied by some loss in sequence specificity. Such systematic comparative studies allow us to begin to establish the thermodynamic database required for the rational design of synthetic polyamides with predictable DNA-binding affinities and specificities.