Evaluation of weak interactions of proteins and organic cations with DNA duplex structures

Molecular interactions and reactions in living cells occur with high background concentrations of organic compounds including proteins. Uncharged water-soluble polymers are commonly used cosolutes in studies on molecular crowding, and most studies argue about the effects of intracellular crowding based on results obtained using polymer cosolutes. Further investigations using protein crowders and organic cations are important in understanding the effects of cellular environments on nucleic acids with negatively charged surfaces. We assessed the effects of using model globular proteins, serum proteins, histone proteins, structurally flexible polypeptides, di- and polyamines, and uncharged polymers. Thermal stability analysis of DNA oligonucleotide structures revealed that unlike conventional polymer cosolutes, basic globular proteins (lysozyme and cytochrome c) at high concentrations stabilized long internal and bulge loop structures but not fully matched duplexes. The selective stabilization of long loop structures suggests preferential binding to unpaired nucleotides in loops through weak electrostatic interactions. Furthermore, the ability of the proteins to stabilize the loop structures was enhanced under macromolecular crowding conditions. Remarkably, the effects of basic proteins on the stability of fully matched duplexes were dissimilar to those of basic amino-acid-rich polypeptides and polyamines. This study provides new insights into the interaction of nucleic acid structures with organic cations.     

Thermodynamic modulation of folding and aggregation energy landscape by DNA binding of functional domains of TDP-43

TDP-43 is a vital nucleic acid binding protein which forms stress-induced aberrant aggregates in around 97% cases of ALS, a fatal neurodegenerative disease. The functional tandem RRM domain of the protein (TDP-43^tRRM) has been shown to undergo amyloid-like aggregation under stress in a pH-dependent fashion. However, the underlying thermodynamic and molecular basis of aggregation and how the energy landscape of folding, stability, and aggregation are coupled and modulated by nucleic acid binding is poorly understood. Here, we show that the pH stress thermodynamically destabilizes the native protein and systematically populates the unfolded-like aggregation-prone molecules which leads to amyloid-like aggregation. We observed that specific DNA binding inhibits aggregation and populates native-like compact monomeric state even under low-pH stress as measured by circular dichroism, ANS binding, size exclusion chromatography, and transmission electron microscopy. We show that DNA-binding thermodynamically stabilizes and populates the native state even under stress and reduces the population of unfolded-like aggregation-prone molecules which leads to systematic aggregation inhibition. Our results suggest that thermodynamic modulation of the folding and aggregation energy landscape by nucleic-acid-like molecules could be a promising approach for effective therapeutic intervention in TDP-43-associated proteinopathies.     

A uniform mechanism correlating dangling-end stabilization and stacking geometry

The geometry of the dangling base in 105 published structures (from X-ray/NMR) containing single-stranded overhangs has been analyzed and correlated to the thermodynamic stabilization found (UV) for the corresponding dangling base/closing basepair combination in short oligonucleotides. The study considers most combinations of closing basepairs, sequence and dangling base residue type, attached in both the 3'- and 5'-ends of both DNA and RNA. Linear regression analysis showed a straightforward correlation (R = 0.873) between the degree of screening for the hydrogen bonds of the closing basepair provided by the dangling base and the resulting thermodynamic stabilization in both DNA and RNA series with dangling ends either at the 3'- or at the 5'-terminus. Regression analysis of only the datasets from RNA gives an improved correlation, R = 0.934, showing that dangling ends on RNA are more ordered than the dangling ends on DNA, R = 0.376. This study highlights the gain in the free energy of stabilization owing to the favorable stacking between the dangling nucleobase and the neighboring basepair and the resulting strengthening of the hydrogen bond of the closing basepair. By acting as a hydrophobic cap on the terminal of the DNA or RNA duplex, the dangling-end residue restricts the bulk water access to the terminal basepair, thereby providing it with a microenvironment devoid of water, which consequently enhances its thermodynamic stability, making it energetically comparable to the corresponding internal basepair. Thus, one single structural model consisting of the interplay of the above electrostatic interactions can be used to explain the molecular basis of the observed thermodynamic effects for dangling-end attachment to the 3'- and 5'-ends of both DNA and RNA duplexes, which is a key step toward accurate dangling-end effect prediction.     

The stability of intramolecular DNA G-quadruplexes compared with other macromolecules

DNA quadruplexes are often conceived as very stable structures. However, most of the free energy of stabilization derives from specific ion binding via inner sphere coordination of the GO6 of the guanine residues comprising the basic quartet. When compared with other nucleic acid structures such as DNA or RNA duplexes and hairpins, or proteins of the same number of atoms, metal-coordinated intramolecular quadruplexes are found to be of comparable or lower thermodynamic stability under similar solution conditions. Furthermore, intramolecular quadruplexes are actually less stable kinetically, than DNA duplexes or hairpins of the same size.     

Unfolding of DNA quadruplexes induced by HIV-1 nucleocapsid protein

The human immunodeficiency virus type 1 nucleocapsid protein (NC) is a nucleic acid chaperone that catalyzes the rearrangement of nucleic acids into their thermodynamically most stable structures. In the present study, a combination of optical and thermodynamic techniques were used to characterize the influence of NC on the secondary structure, thermal stability and energetics of monomolecular DNA quadruplexes formed by the sequence d(GGTTGGTGTGGTTGG) in the presence of K⁺ or Sr²⁺. Circular dichroism studies demonstrate that NC effectively unfolds the quadruplexes. Studies carried out with NC variants suggest that destabilization is mediated by the zinc fingers of NC. Calorimetric studies reveal that NC destabilization is enthalpic in origin, probably owing to unstacking of the G-quartets upon protein binding. In contrast, parallel studies performed on a related DNA duplex reveal that under conditions where NC readily destabilizes and unfolds the quadruplexes, its effect on the DNA duplex is much less pronounced. The differences in NC's ability to destabilize quadruplex versus duplex is in accordance with the higher ΔG of melting for the latter, and with the inverse correlation between nucleic acid stability and the destabilizing activity of NC.     

Conformation and the sodium ion condensation on DNA and RNA structures in the presence of a neutral cosolute as a mimic of the intracellular media

Water-soluble neutral cosolutes can be used to quantify biomolecular properties in the particular molecular environment occurring in a cell. We studied the conformation and the thermal stability of DNA and RNA structures in the presence of PEG [poly(ethylene glycol)] and smaller cosolutes of glycerol, ethylene glycol, 1,3-propanediol, 2-methoxyethanol, and 1,2-dimethoxyethane. Although the neutral cosolutes destabilized the oligonucleotide duplex and the hairpin structures, the left-handed Z-form duplex was more energetically favored in the cosolute-containing solutions. These observations were due to the contribution of water molecule on the nucleotide structure formations because the cosolutes act as an osmolyte to reduce the water activity of a solution. Moreover, the sodium ion condensation for the duplex and the hairpin formations was reduced in the presence of PEG, while that for the transition from the B-form to the Z-form was unaltered. The CD (circular dichroism) and EPR (electron paramagnetic resonance) spectra demonstrated that the cosolutes changed the helical conformation of the unstructured oligonucleotides, but not those of the ordered structures. The results of the favorable formation of the noncanonical nucleotide structures, and minimized conformational and thermal perturbations of the ordered nucleotide structures in the cosolute-containing solutions implicate the significance of the intracellular environment on DNA and RNA structures in a cell.     

Predicting the effects of basepair mutations in DNA-protein complexes by thermodynamic integration

Thermodynamically rigorous free energy methods in principle allow the exact computation of binding free energies in biological systems. Here, we use thermodynamic integration together with molecular dynamics simulations of a DNA-protein complex to compute relative binding free energies of a series of mutants of a protein-binding DNA operator sequence. A guanine-cytosine basepair that interacts strongly with the DNA-binding protein is mutated into adenine-thymine, cytosine-guanine, and thymine-adenine. It is shown that basepair mutations can be performed using a conservative protocol that gives error estimates of ∼10% of the change in free energy of binding. Despite the high CPU-time requirements, this work opens the exciting opportunity of being able to perform basepair scans to investigate protein-DNA binding specificity in great detail computationally.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Unique properties of purine/pyrimidine asymmetric PNA.DNA duplexes: differential stabilization of PNA.DNA duplexes by purines in the PNA strand

PNA.DNA duplexes are significantly stabilized by purine nucleobases in the PNA strand. To elucidate and understand the effect of switching the backbone in a nucleic acid duplex, we now report a thermodynamics study along with a solution conformations study of two purine/pyrimidine strand asymmetric duplexes and a strand symmetrical control by comparing the behavior of all four possible PNA/DNA combinations. In essence, we are comparing an identical basepair stack connected by either an aminoethyl glycine PNA or a deoxyribose DNA backbone. We show that the PNA.DNA duplexes containing purine-rich PNA strands are stabilized with regard to the thermal melting temperature and free energy as well as enthalpy (and concomitantly relatively less entropically disfavored). Based on our data, we find it unlikely that differences in counterion binding (identical ionic-strength dependence was observed), hydration (identical and insignificant water release was observed), or single-strand conformation can be responsible for the difference in duplex stability. The only consistent difference observed between the purine-rich PNA versus the pyrimidine-rich PNA in isosequential PNA.DNA duplexes is the significant increase in both binding enthalpy and entropy for the PNA.DNA duplexes containing pyrimidine-rich PNA in organic solvent, which would indicate that these duplexes are relatively enthalpically disfavored in water. Although our results so far do not allow us to identify the origin of the different stabilities of homopurine/homopyrimidine PNA.DNA duplexes, the evidence does point to a significant structural component, which involves enthalpic contributions both within the duplex structure and also from bound water molecules.     

A sequence-independent analysis of the loop length dependence of intramolecular RNA G-quadruplex stability and topology

G-Quadruplexes are noncanonical nucleic acid secondary structures based on guanine association that are readily adopted by G-rich RNA and DNA sequences. Naturally occurring genomic G-quadruplex-forming sequences have functional roles in biology that are mediated through structure. To appreciate how this is achieved, an understanding of the likelihood of G-quadruplex formation and the structural features of the folded species under a defined set of conditions is informative. We previously systematically investigated the thermodynamic stability and folding topology of DNA G-quadruplexes and found a strong dependence of these properties on loop length and loop arrangement [Bugaut, A., and Balasubramanian, S. (2008) Biochemistry 47, 689-697]. Here we report on a complementary analysis of RNA G-quadruplexes using UV melting and circular dichroism spectroscopy that also serves as a comparison to the equivalent DNA G-quadruplex-forming sequences. We found that the thermodynamic stability of G-quadruplex RNA can be modulated by loop length while the overall structure is largely unaffected. The systematic design of our study also revealed subtle loop length dependencies in RNA G-quadruplex structure.     

Differences between EcoRI nonspecific and "star" sequence complexes revealed by osmotic stress

The binding of the restriction endonuclease EcoRI to DNA is exceptionally specific. Even a single basepair change ("star" sequence) from the recognition sequence, GAATTC, decreases the binding free energy of EcoRI to values nearly indistinguishable from nonspecific binding. The difference in the number of waters sequestered by the protein-DNA complexes of the "star" sequences TAATTC and CAATTC and by the specific sequence complex determined from the dependence of binding free energy on water activity is also practically indistinguishable at low osmotic pressures from the 110 water molecules sequestered by nonspecific sequence complexes. Novel measurements of the dissociation rates of noncognate sequence complexes and competition equilibrium show that sequestered water can be removed from "star" sequence complexes by high osmotic pressure, but not from a nonspecific complex. By 5 Osm, the TAATTC "star" sequence complex has lost almost 90 of the approximately 110 waters initially present. It is more difficult to remove water from the CAATTC "star" sequence complex. The sequence dependence of water loss correlates with the known sequence dependence of "star" cleavage activity.     

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

The Relationship between the Effects of Cytokinin on the Expension and Metabolism of Excised Cucumber Cotyledons and Water Stress

Excised cucumber (Cucumis sativus, var. Jin-ian No. 4) cotyledons were incubated with 10 ppm of BA (benzyladenine) or water for 1 h, then thouroughly rinsed with water and grown in darkness on filter paper saturated with different concentrations of mannitol solution. Up to 24 h, the fresh weight, carotenoid, RNA, DNA and lipid contents of cotyledons were determined. Although mannitol solution reduced the effectiveness of BA treatment, in the same condition of osmotic potential, the increases of fresh weight, carotenoid, RNA and DNA contents, as well as the decrease of lipid per cotyledon were always much higher in BA treated tissues. BA enhanced the rate of water uptake by the cotyledons. The fresh weight of BA and 0.2 M mannitol treated cotyledons was equal to that of water control, but the increases of carotenoid and nucleic acid contents and the decrease of lipid were much higher in tho former than the latter. 0.3 M and 0.5 M mannitol solu- tions almost interrupted the water uptake of water and BA treated cotyledons respectively. However, the increases of carotenoid and nucleic acid contents as well as decrease of lipid were still occurred in these conditions. The different osmotic potential did nearly not affect the ratio of the increases of carotenoid and nucleic acid contents between BA treatment and control. It means the effectiveness of BA was almost the same under different osmotic potential It is evident that BA stimulated simultaneously the water uptake and metabolism of the cotyledons. They are probably different processes but closely related to each other.     

Effects of a protecting osmolyte on the ion atmosphere surrounding DNA duplexes

Osmolytes are small, chemically diverse, organic solutes that function as an essential component of cellular stress response. Protecting osmolytes enhance protein stability via preferential exclusion, and nonprotecting osmolytes, such as urea, destabilize protein structures. Although much is known about osmolyte effects on proteins, less is understood about osmolyte effects on nucleic acids and their counterion atmospheres. Nonprotecting osmolytes destabilize nucleic acid structures, but effects of protecting osmolytes depend on numerous factors including the type of nucleic acid and the complexity of the functional fold. To begin quantifying protecting osmolyte effects on nucleic acid interactions, we used small-angle X-ray scattering (SAXS) techniques to monitor DNA duplexes in the presence of sucrose. This protecting osmolyte is a commonly used contrast matching agent in SAXS studies of protein-nucleic acid complexes; thus, it is important to characterize interaction changes induced by sucrose. Measurements of interactions between duplexes showed no dependence on the presence of up to 30% sucrose, except under high Mg(2+) conditions where stacking interactions were disfavored. The number of excess ions associated with DNA duplexes, reported by anomalous small-angle X-ray scattering (ASAXS) experiments, was sucrose independent. Although protecting osmolytes can destabilize secondary structures, our results suggest that ion atmospheres of individual duplexes remain unperturbed by sucrose.     

The effect of water on the rate of conformational change in protein allostery

The influence of solvation on the rate of quaternary structural change is investigated in human hemoglobin, an allosteric protein in which reduced water activity destabilizes the R state relative to T. Nanosecond absorption spectroscopy of the heme Soret band was used to monitor protein relaxation after photodissociation of aqueous HbCO complex under osmotic stress induced by the nonbinding cosolute poly(ethylene glycol) (PEG). Photolysis data were analyzed globally for six exponential time constants and amplitudes as a function of osmotic stress and viscosity. Increases in time constants associated with geminate rebinding, tertiary relaxation, and quaternary relaxation were observed in the presence of PEG, along with a decrease in the fraction of hemes rebinding CO with the slow rate constant characteristic of the T state. An analysis of these results along with those obtained by others for small cosolutes showed that both osmotic stress and solvent viscosity are important determinants of the microscopic R --> T rate constant. The size and direction of the osmotic stress effect suggests that at least nine additional water molecules are required to solvate the allosteric transition state relative to the R-state hydration, implying that the transition state has a greater solvent-exposed area than either end state.     

Influence of 5-N-carboxamide modifications on the thermodynamic stability of oligonucleotides

We have recently shown that the incorporation of modified nucleotides such as 5-N-carboxamide-deoxyuridines into random nucleic acid libraries improves success rates in SELEX experiments and facilitates the identification of ligands with slow off-rates. Here we report the impact of these modifications on the thermodynamic stability of both duplexes and intramolecular ‘single-stranded’ structures. Within duplexes, large, hydrophobic naphthyl groups were destabilizing relative to the all natural DNA duplex, while the hydrophilic groups exhibited somewhat improved duplex stability. All of the significant changes in stability were driven by opposing contributions from the enthalpic and entropic terms. In contrast, both benzyl and naphthyl modifications stabilized intramolecular single-stranded structures relative to their natural DNA analogs, consistent with the notion that intramolecular folding allows formation of novel, stabilizing hydrophobic interactions. Imino proton NMR data provided evidence that elements of the folded structure form at temperatures well below the T_m, with a melting transition that is distinctly less cooperative when compared to duplex DNA. Although there are no data to suggest that the unmodified DNA sequences fold into structures similar to their modified analogs, this still represents clear evidence that these modifications impart thermodynamic stability to the folded structure not achievable with unmodified DNA.