Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Predicting crossover generation in DNA shuffling

We introduce a quantitative framework for assessing the generation of crossovers in DNA shuffling experiments. The approach uses free energy calculations and complete sequence information to model the annealing process. Statistics obtained for the annealing events then are combined with a reassembly algorithm to infer crossover allocation in the reassembled sequences. The fraction of reassembled sequences containing zero, one, two, or more crossovers and the probability that a given nucleotide position in a reassembled sequence is the site of a crossover event are estimated. Comparisons of the predictions against experimental data for five example systems demonstrate good agreement despite the fact that no adjustable parameters are used. An in silico case study of a set of 12 subtilases examines the effect of fragmentation length, annealing temperature, sequence identity and number of shuffled sequences on the number, type, and distribution of crossovers. A computational verification of crossover aggregation in regions of near-perfect sequence identity and the presence of synergistic reassembly in family DNA shuffling is obtained.

A signature of the T → R transition in human hemoglobin

Allosteric effects in hemoglobin arise from the equilibrium between at least two energetic states of the molecule: a tense state, T, and a relaxed state, R. The two states differ from each other in the number and energy of the interactions between hemoglobin subunits. In the T state, constraints between subunits oppose the structural changes resulting from ligand binding. In the R state, these constraints are released, thus enhancing ligand-binding affinity. In the present work, we report the presence of four sites in hemoglobin that are structurally stabilized in the R relative to the T state. These sites are Hisα103(G10) and Hisα122(H5) in each α subunit of hemoglobin. They are located at the α1β1 and α2β2 interfaces of the hemoglobin tetramer, where the histidine side chains form hydrogen bonds with specific residues from the β chains. We have measured the solvent exchange rates of side chain protons of Hisα103(G10) and Hisα122(H5) in both deoxygenated and ligated hemoglobin by NMR spectroscopy. The exchange rates were found to be higher in the deoxygenated-T than in ligated-R state. Analysis of exchange rates in terms of the local unfolding model revealed that the structural stabilization free energy at each of these two histidines is larger by ≈1.5 kcal/(mol tetramer) in the R relative to the T state. The location of these histidines at the intradimeric α1β1 and α2β2 interfaces also suggests a role for these interfaces in the allosteric equilibrium of hemoglobin.

Site-specifically located 8-amino-2′-deoxyguanosine: thermodynamic stability and mutagenic properties in Escherichia coli

2-Nitropropane (2-NP), an important industrial solvent and a component of cigarette smoke, is mutagenic in bacteria and carcinogenic in rats. 8-Amino-2′-deoxyguanosine (8-amino-dG) is one of the types of DNA damage found in liver, the target organ in 2-NP-treated rats. To investigate the thermodynamic properties of 8-amino-dG opposite each of the four DNA bases, we have synthesized an 11mer, d(CCATCG*CTACC), in which G* represents the modified base. By annealing a complementary DNA strand to this modified 11mer, four sets of duplexes were generated each containing one of the four DNA bases opposite the lesion. Circular dichroism studies indicated that 8-amino-dG did not alter the global helical properties of natural right-handed B-DNA. The thermal stability of each duplex was examined by UV melting measurements and compared with its unmodified counterpart. For the unmodified 11mer, the relative stability of the complementary DNA bases opposite G was in the order C > T > G > A, as determined from their –ΔG° values. The free energy change of each modified duplex was lower than its unmodified counterpart, except for the G*:G pair that exhibited a higher melting transition and a larger –ΔG° than the G:G duplex. Nevertheless, the stability of the modified 11mer duplex also followed the order C > T > G > A when placed opposite 8-amino-dG. To explore if 8-amino-dG opposite another 8-amino-dG has any advantage in base pairing, a G*:G* duplex was evaluated, which showed that the stability of this duplex was similar to the G*:G duplex. Mutagenesis of 8-amino-dG in this sequence context was studied in Escherichia coli, which showed that the lesion is weakly mutagenic (mutation frequency ∼10–3) but still can induce a variety of targeted and semi-targeted mutations.

Deuterium isotope effect on the intramolecular electron transfer in Pseudomonas aeruginosa azurin

Intramolecular electron transfer in azurin in water and deuterium oxide has been studied over a broad temperature range. The kinetic deuterium isotope effect, kH/kD, is smaller than unity (0.7 at 298 K), primarily caused by the different activation entropies in water (−56.5 J K−1 mol−1) and in deuterium oxide (−35.7 J K−1 mol−1). This difference suggests a role for distinct protein solvation in the two media, which is supported by the results of voltammetric measurements: the reduction potential (E0′) of Cu2+/+ at 298 K is 10 mV more positive in D2O than in H2O. The temperature dependence of E0′ is also different, yielding entropy changes of −57 J K−1 mol−1 in water and −84 J K−1 mol−1 in deuterium oxide. The driving force difference of 10 mV is in keeping with the kinetic isotope effect, but the contribution to ΔS‡ from the temperature dependence of E0′ is positive rather than negative. Isotope effects are, however, also inherent in the nuclear reorganization Gibbs free energy and in the tunneling factor for the electron transfer process. A slightly larger thermal protein expansion in H2O than in D2O (0.001 nm K−1) is sufficient both to account for the activation entropy difference and to compensate for the different temperature dependencies of E0′. Thus, differences in driving force and thermal expansion appear as the most straightforward rationale for the observed isotope effect.

Balanced branching in transcription termination

The theory of stochastic transcription termination based on free-energy competition [von Hippel, P. H. & Yager, T. D. (1992) Science 255, 809–812 and von Hippel, P. H. & Yager, T. D. (1991) Proc. Natl. Acad. Sci. USA 88, 2307–2311] requires two or more reaction rates to be delicately balanced over a wide range of physical conditions. A large body of work on glasses and large molecules suggests that this balancing should be impossible in such a large system in the absence of a new organizing principle of matter. We review the experimental literature of termination and find no evidence for such a principle, but do find many troubling inconsistencies, most notably, anomalous memory effects. These effects suggest that termination has a deterministic component and may conceivably not be stochastic at all. We find that a key experiment by Wilson and von Hippel [Wilson, K. S. & von Hippel, P. H. (1994) J. Mol. Biol. 244, 36–51] thought to demonstrate stochastic termination was an incorrectly analyzed regulatory effect of Mg2+ binding.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

Mechanism for nucleic acid chaperone activity of HIV-1 nucleocapsid protein revealed by single molecule stretching

The nucleocapsid protein (NC) of HIV type 1 is a nucleic acid chaperone that facilitates the rearrangement of nucleic acids into conformations containing the maximum number of complementary base pairs. We use an optical tweezers instrument to stretch single DNA molecules from the helix to coil state at room temperature in the presence of NC and a mutant form (SSHS NC) that lacks the two zinc finger structures present in NC. Although both NC and SSHS NC facilitate annealing of complementary strands through electrostatic attraction, only NC destabilizes the helical form of DNA and reduces the cooperativity of the helix-coil transition. In particular, we find that the helix-coil transition free energy at room temperature is significantly reduced in the presence of NC. Thus, upon NC binding, it is likely that thermodynamic fluctuations cause continuous melting and reannealing of base pairs so that DNA strands are able to rapidly sample configurations to find the lowest energy state. The reduced cooperativity allows these fluctuations to occur in the middle of complex double-stranded structures. The reduced stability and cooperativity, coupled with the electrostatic attraction generated by the high charge density of NC, is responsible for the nucleic acid chaperone activity of this protein.

Monitoring the GAP catalyzed H-Ras GTPase reaction at atomic resolution in real time

The molecular reaction mechanism of the GTPase-activating protein (GAP)-catalyzed GTP hydrolysis by Ras was investigated by time resolved Fourier transform infrared (FTIR) difference spectroscopy using caged GTP (P3-1-(2-nitro)phenylethyl guanosine 5′-O-triphosphate) as photolabile trigger. This approach provides the complete GTPase reaction pathway with time resolution of milliseconds at the atomic level. Up to now, one structural model of the GAP⋅Ras⋅GDP⋅AlFx transition state analog is known, which represents a “snap shot” along the reaction-pathway. As now revealed, binding of GAP to Ras⋅GTP shifts negative charge from the γ to β phosphate. Such a shift was already identified by FTIR in GTP because of Ras binding and is now shown to be enhanced by GAP binding. Because the charge distribution of the GAP⋅Ras⋅GTP complex thus resembles a more dissociative-like transition state and is more like that in GDP, the activation free energy is reduced. An intermediate is observed on the reaction pathway that appears when the bond between β and γ phosphate is cleaved. In the intermediate, the released Pi is strongly bound to the protein and surprisingly shows bands typical of those seen for phosphorylated enzyme intermediates. All these results provide a mechanistic picture that is different from the intrinsic GTPase reaction of Ras. FTIR analysis reveals the release of Pi from the protein complex as the rate-limiting step for the GAP-catalyzed reaction. The approach presented allows the study not only of single proteins but of protein–protein interactions without intrinsic chromophores, in the non-crystalline state, in real time at the atomic level.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.

Polarity and permeation profiles in lipid membranes

The isotropic 14N-hyperfine coupling constant, a\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{o}^{N}}}\end{equation*}\end{document}, of nitroxide spin labels is dependent on the local environmental polarity. The dependence of a\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{o}^{N}}}\end{equation*}\end{document} in fluid phospholipid bilayer membranes on the C-atom position, n, of the nitroxide in the sn-2 chain of a spin-labeled diacyl glycerophospholipid therefore determines the transmembrane polarity profile. The polarity variation in phospholipid membranes, with and without equimolar cholesterol, is characterized by a sigmoidal, trough-like profile of the form {1 + exp [(n − no)/λ]}−1, where n = no is the point of maximum gradient, or polarity midpoint, beyond which the free energy of permeation decreases linearly with n, on a characteristic length-scale, λ. Integration over this profile yields a corresponding expression for the permeability barrier to polar solutes. For fluid membranes without cholesterol, no ≈ 8 and λ ≈ 0.5–1 CH2 units, and the permeability barrier introduces an additional diffusive resistance that is equivalent to increasing the effective membrane thickness by 35–80%, depending on the lipid. For membranes containing equimolar cholesterol, no ≈ 9–10, and the total change in polarity is greater than for membranes without cholesterol, increasing the permeability barrier by a factor of 2, whereas the decay length remains similar. The permeation of oxygen into fluid lipid membranes (determined by spin-label relaxation enhancements) displays a profile similar to that of the transmembrane polarity but of opposite sense. For fluid membranes without cholesterol no ≈ 8 and λ ≈ 1 CH2 units, also for oxygen. The permeation profile for polar paramagnetic ion complexes is closer to a single exponential decay, i.e., no lies outside the acyl-chain region of the membrane. These results are relevant not only to the permeation of water and polar solutes into membranes and their permeabilities, but also to depth determinations of site-specifically spin-labeled protein residues by using paramagnetic relaxation agents.

Differences in water release for the binding of EcoRI to specific and nonspecific DNA sequences.

The free energy difference between complexes of the restriction nuclease EcoRI with nonspecific DNA and with the enzyme's recognition sequence is linearly dependent on the water chemical potential of the solution, set using several very different solutes, ranging from glycine and glycerol to triethylene glycol and sucrose. This osmotic dependence indicates that the nonspecific complex sequesters some 110 waters more than the specific complex with the recognition sequence. The insensitivity of the difference in number of waters released to the solute identity further indicates that this water is sequestered in a space that is sterically inaccessible to solutes, most likely at the protein-DNA interface of the nonspecific complex. Calculations based on the structure of the specific complex suggest that the apposing DNA and protein surfaces in the nonspecific complex retain approximately a full hydration layer of water.

Mimicry of the calcium-induced conformational state of troponin C by low temperature under pressure.

Calcium binding to the N-domain of troponin C initiates a series of conformational changes that lead to muscle contraction. Calcium binding provides the free energy for a hydrophobic region in the core of N-domain to assume a more open configuration. Fluorescence measurements on a tryptophan mutant (F29W) show that a similar conformational change occurs in the absence of Ca2+ when the temperature is lowered under pressure. The conformation induced by subzero temperatures binds the hydrophobic probe bis-aminonaphthalene sulfonate, and the tryptophan has the same fluorescence lifetime (7 ns) as in the Ca2+-bound form. The decrease in volume (delta V = -25.4 ml/mol) corresponds to an increase in surface area. Thermodynamic measurements suggest an enthalpy-driven conformational change that leads to an intermediate with an exposed N-domain core and a high affinity for Ca2+.

On the mechanism of action of ribonuclease A: relevance of enzymatic studies with a p-nitrophenylphosphate ester and a thiophosphate ester.

It has been reported that His-119 of ribonuclease A plays a major role as an imidazolium ion acid catalyst in the cyclization/cleavage of normal dinucleotides but that it is not needed for the cyclization/cleavage of 3'-uridyl p-nitrophenyl phosphate. We see that this is also true for simple buffer catalysis, where imidazole (as in His-12 of the enzyme), but not imidazolium ion, plays a significant catalytic role with the nitrophenyl substrate, but both are catalytic for normal dinucleotides such as uridyluridine. Rate studies show that the enzyme catalyzes the cyclization of the nitrophenylphosphate derivative 47,000,000 times less effectively (kcat/kuncat) than it does uridyladenosine, indicating that approximately 50% of the catalytic free energy change is lost with this substrate. This suggests that the nitrophenyl substrate is not correctly bound to take full advantage of the catalytic groups of the enzyme and is thus not a good guide to the mechanism used by normal nucleotides. The published data on kinetic effects with ribonuclease A of substituting thiophosphate groups for the phosphate groups of normal substrates has been discussed elsewhere, and it was argued that these effects are suggestive of the classical mechanism for ribonuclease action, not the novel mechanism we have recently proposed. The details of these rate effects, including stereochemical preferences in the thiophosphate series, can be invoked as support for our newer mechanism.