Structure of Met-enkephalin in explicit aqueous solution using replica exchange molecular dynamics

Replica exchange molecular dynamics (MD) simulations of Met-enkephalin in explicit solvent reveal helical and nonhelical structures. Four predominant structures of Met-enkephalin are sampled with comparable probabilities (two helical and two nonhelical). The energy barriers between these configurations are low, suggesting that Met-enkephalin switches easily between configurations. This is consistent with the requirement that Met-enkephalin be sufficiently flexible to bind to several different receptors. Replica exchange simulations of 32 ns are shown to sample approximately five times more configurational space than constant temperature MD simulations of the same duration. The energy landscape for the replica exchange simulation is presented. A detailed study of replica trajectories demonstrates that the significant increases in temperature provided by the replica exchange technique enable transitions from nonhelical to helical structures that would otherwise be prevented by kinetic trapping. Met-enkephalin (Type Entrez Proteins; Value A61445; Service Entrez Proteins).

Molecular dynamics simulations of pressure effects on hydrophobic interactions

We report results on the pressure effects on hydrophobic interactions obtained from molecular dynamics simulations of aqueous solutions of methanes in water. A wide range of pressures that is relevant to pressure denaturation of proteins is investigated. The characteristic features of water-mediated interactions between hydrophobic solutes are found to be pressure-dependent. In particular, with increasing pressure we find that (1) the solvent-separated configurations in the solute-solute potential of mean force (PMF) are stabilized with respect to the contact configurations; (2) the desolvation barrier increases monotonically with respect to both contact and solvent-separated configurations; (3) the locations of the minima and the barrier move toward shorter separations; and (4) pressure effects are considerably amplified for larger hydrophobic solutes. Together, these observations lend strong support to the picture of the pressure denaturation process proposed previously by Hummer et al. (Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 1552): with increasing pressure, the transfer of water into protein interior becomes key to the pressure denaturation process, leading to the dissociation of close hydrophobic contacts and subsequent swelling of the hydrophobic protein interior through insertions of water molecules. The pressure dependence of the PMF between larger hydrophobic solutes shows that pressure effects on the interaction between hydrophobic amino acids may be considerably amplified compared to those on the methane-methane PMF.

Amino acid substitution at position 99 affects the rate of CRP subunit exchange

We investigated the characteristics of CRP having amino acid substitutions at position 99. Analysis of amino acid residue proximity to cAMP in molecular dynamics (MD) simulations of the CRP:(cAMP)(2) complex [García, A. E., and Harman, J. G. (1996) Protein Sci. 5, 62-71] showed repositioning of tyrosine 99 (Y99) to interact with the equatorial exocyclic oxygen atom of cAMP. To test the role of Y99 in cAMP-mediated CRP activation, Y99 was substituted with alanine (A) or phenylalanine (F). Cells that contained the WT or mutant forms of CRP induced beta-galactosidase in the presence of cAMP. Purified WT, Y99A, and Y99F CRP showed only a 3- to 4-fold difference in cAMP affinity. There were no apparent differences between the three forms of CRP in cAMP binding cooperativity, in CRP:(cAMP)(1) complex binding to lacP DNA, in the formation of CRP:cAMP:RNAP complexes at lacP, or in CRP efficacy in mediating lacP activity in vitro. The apo-form of Y99A CRP was more sensitive to protease than the apo-form of either WT CRP or Y99F CRP. Whereas the WT or Y99F CRP:(cAMP)(1) complexes were cleaved by protease at hinge-region peptide bonds, the Y99A CRP:(cAMP)(1) complex was cleaved at peptide bonds located at the subunit interface. The rates of subunit exchange for Y99A CRP, both in the apo-form and in a 1:1 complex with cAMP, were significantly greater than that measured for WT CRP. The results of this study show that tyrosine 99 contributes significant structural stability to the CRP dimer, specifically in stabilizing subunit association.

Functional dynamics of the hydrophobic cleft in the N-domain of calmodulin

Molecular dynamics studies of the N-domain (amino acids 1-77; CaM(1-77)) of Ca2+-loaded calmodulin (CaM) show that a solvent exposed hydrophobic cleft in the crystal structure of CaM exhibits transitions from an exposed (open) to a buried (closed) state over a time scale of nanoseconds. As a consequence of burying the hydrophobic cleft, the R(g) of the protein is reduced by 1.5 A. Based on this prediction, x-ray scattering experiments were conducted on this domain over a range of concentrations. Models built from the scattering data show that the R(g) and general shape is consistent with the simulation studies of CaM(1-77). Based on these observations we postulate a model in which the conformation of CaM fluctuates between two different states that expose and bury this hydrophobic cleft. In aqueous solution the closed state dominates the population, while in the presence of peptides, the open state dominates. This inherent flexibility of CaM may be the key to its versatility in recognizing structurally distinct peptide sequences. This model conflicts with the currently accepted hypothesis based on observations in the crystal structure, where upon Ca2+ binding the hydrophobic cleft is exposed to solvent. We postulate that crystal packing forces stabilize the protein conformation toward the open configuration.

Exploring the energy landscape of a beta hairpin in explicit solvent

We studied the energy landscape of the peptide Ace-GEWTYDDATKTFTVTE-Nme, taken from the C-terminal fragment (41-56) of protein G, in explicit aqueous solution by a highly parallel replica-exchange approach that combines molecular dynamics trajectories with a temperature exchange Monte Carlo process. The combined trajectories in T and configurational space allow a replica to overcome a free energy barrier present at one temperature by increasing T, changing configurations, and cooling in a self-regulated manner, thus allowing sampling of broad regions of configurational space in short (nanoseconds) time scales. The free energy landscape of this system over a wide range of temperatures shows that the system preferentially adopts a beta hairpin structure. However, the peptide also samples other stable ensembles where the peptide adopts helices and helix-turn-helix states, among others. The helical states become increasingly stable at low temperatures, but are slightly less stable than the beta turn ensemble. The energy landscape is rugged at low T, where substates are separated by large energy barriers. These barriers disappear at higher T (approximately 330 K), where the system preferentially adopts a "molten globule" state with structures similar to the beta hairpin.

Helix nucleation kinetics from molecular simulations in explicit solvent

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.

Conformational diffusion and helix formation kinetics

The time, temperature, and sequence dependences of helix formation kinetics of fully atomistic peptide models in explicit solvent are described quantitatively by a diffusive search within the coil state with barrierless transitions into the helical state. Conformational diffusion leads to nonexponential kinetics and jump-width dependences in temperature jump experiments.

Water penetration and escape in proteins

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.

Pressure-induced protein-folding/unfolding kinetics

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

Conformational dynamics of cytochrome c: correlation to hydrogen exchange

We study the dynamical fluctuations of horse heart cytochrome c by molecular dynamics (MD) simulations in aqueous solution, at four temperatures: 300 K, 360 K, 430 K, and 550 K. Each simulation covers a production time of at least 1.5 nanoseconds (ns). The conformational dynamics of the system is analyzed in terms of collective motions that involve the whole protein, and local motions that involve the formation and breaking of intramolecular hydrogen bonds. The character of the MD trajectories can be described within the framework of rugged energy landscape dynamics. The MD trajectories sample multiple conformational minima, with basins in protein conformational space being sampled for a few hundred picoseconds. The trajectories of the system in configurational space can be described in terms of diffusion of a particle in real space with a waiting time distribution due to partial trapping in shallow minima. As a consequence of the hierarchical nature of the dynamics, the mean square displacement autocorrelation function, <|x(t) - x(0)|2>, exhibits a power law dependence on time, with an exponent of around 0.5 for times shorter than 100 ps, and an exponent of 1.75 for longer times. This power law behavior indicates that the system exhibits suppressed diffusion (sub-diffusion) in sampling of configurational space at time scales shorter than 100 ps, and enhanced (super-diffusion) at longer time scales. The multi-basin feature of the trajectories is present at all temperatures simulated. Structural changes associated with inter-basin displacements correspond to collective motions of the Omega loops and coiled regions and relative motions of the alpha-helices as rigid bodies. Similar motions may be involved in experimentally observed amide hydrogen exchange. However, some groups showing large correlated motions do not expose the amino hydrogens to the solvent. We show that large fluctuations are not necessarily correlated to hydrogen exchange. For example, regions of the proteins forming alpha helices and turns show significant fluctuations, but as rigid bodies, and the hydrogen bonds involved in the formation of these structures do not break in proportion to these fluctuations. Proteins 1999;36:175-191. Published 1999 Wiley-Liss, Inc.

Temperature dependence of the solubility of non-polar gases in water

An explanation is provided for the experimentally observed temperature dependence of the solubility and the solubility minimum of non-polar gases in water. The influence of solute size and solute-water attractive interactions on the solubility minimum temperature is investigated. The transfer of a non-polar solute from the ideal gas into water is divided into two steps: formation of a cavity in water with the size and shape of the solute and insertion of the solute in this cavity which is equivalent to 'turning on' solute-water attractive interactions. This two step process divides the excess chemical potential of the non-polar solute in water into repulsive and attractive contributions, respectively. The reversible work for cavity formation is modeled using an information theory model of hydrophobic hydration. Attractive contributions are calculated by modeling the water structure in the vicinity of non-polar solutes. These models make a direct connection between microscopic quantities and macroscopic observables. Moreover, they provide an understanding of the peculiar temperature dependences of the hydration thermodynamics from properties of pure water; specifically, bulk water density and the water oxygen-oxygen radial distribution function.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a beta-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

The pressure dependence of hydrophobic interactions is consistent with the observed pressure denaturation of proteins

Proteins can be denatured by pressures of a few hundred MPa. This finding apparently contradicts the most widely used model of protein stability, where the formation of a hydrophobic core drives protein folding. The pressure denaturation puzzle is resolved by focusing on the pressure-dependent transfer of water into the protein interior, in contrast to the transfer of nonpolar residues into water, the approach commonly taken in models of protein unfolding. Pressure denaturation of proteins can then be explained by the pressure destabilization of hydrophobic aggregates by using an information theory model of hydrophobic interactions. Pressure-denatured proteins, unlike heat-denatured proteins, retain a compact structure with water molecules penetrating their core. Activation volumes for hydrophobic contributions to protein folding and unfolding kinetics are positive. Clathrate hydrates are predicted to form by virtually the same mechanism that drives pressure denaturation of proteins.

Variations on a theme by Debye and Waller: from simple crystals to proteins

Debye and Waller showed how to adjust scattering intensities in diffraction experiments for harmonic motions of atoms about an average, static reference configuration. However, many motions, particularly in biological molecules as compared to simple crystals, are far from harmonic. We show how, using a variety of simple anharmonic, multiconformational models, it is possible to construct a variety of Generalized Debye-Waller Factors, and understand their meaning. A central result for these cases is that, in principle, intensity factors cannot be obtained from true total mean square displacements of the atoms. We make the distinction between the intensity factors for unimodal quasiharmonic motions and those for the anharmonic, multimodal (valley hopping) motions. Only the former affect the conventional B factors.

Hydration of an alpha-helical peptide: comparison of theory and molecular dynamics simulation

We present a statistical mechanical description of biomolecular hydration that accurately describes the hydrophobic and hydrophilic hydration of a model alpha-helical peptide. The local density of water molecules around a biomolecule is obtained by means of a potential-of-mean-force (PMF) expansion in terms of pair- and triplet-correlation functions of bulk water and dilute solutions of nonpolar atoms. The accuracy of the method is verified by comparing PMF results with the local density and site-site correlation functions obtained by molecular dynamics simulations of a model alpha-helix in solution. The PMF approach quantitatively reproduces all features of the peptide hydration determined from the molecular dynamics simulation. Regions of hydrophobic hydration near the C alpha and C beta atoms along the helix are well reproduced. The hydration of exposed polar groups at the N- and C-termini of the helix are also well described by the theory. A detailed comparison of the local hydration by means of site-site radial distribution functions evaluated with the PMF theory shows agreement with the molecular dynamics simulations. The formulation of this theory is general and can be applied to any biomolecular system. The accuracy, speed of computation, and local character of this theory make it especially suitable for studying large biomolecular systems.

An information theory model of hydrophobic interactions

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

Three-dimensional model of a selective theophylline-binding RNA molecule

A three-dimensional (3D) model for an RNA molecule that selectively binds theophylline but not caffeine is proposed. This RNA, which was found using SELEX (Jenison et al., 1994), is 10,000 times more specific for theophylline (Kn = 320 nM) than for caffeine (KD = 3.5 mM), although the two ligands are identical except for a methyl group substituted at N7 (present only in caffeine). The binding affinity for ten xanthine-based ligands was used to derive a comparative molecular field analysis model (R2 = 0.93 for three components, with cross-validated R2 of 0.73), using the SYBYL and GOLPE programs. A pharmacophoric map was generated to locate steric and electrostatic interactions between theophylline and the RNA binding site. This information was used to identify putative functional groups of the binding pocket and to generate distance constraints. On the basis of a model for the secondary structure (Jenison et al., 1994), the 3D structure of this RNA was then generated using the following method: each helical region of the RNA molecule was treated as a rigid body; single-stranded loops with specific end-to-end distances were generated. The structures of RNA-xanthine complexes were studied using a modified Monte Carlo algorithm. The detailed structure of an RNA-ligand complex model, as well as possible explanations for the theophylline selectivity are discussed.

Three-dimensional quantitative structure-activity relationships of steroid aromatase inhibitors

Inhibition of aromatase, a cytochrome P450 that converts androgens to estrogens, is relevant in the therapeutic control of breast cancer. We investigate this inhibition using a three-dimensional quantitative structure-activity relationship (3D QSAR) method known as Comparative Molecular Field Analysis, CoMFA [Cramer III, R.D. et al., J. Am. Chem. Soc., 110 (1988) 5959]. We analyzed the data for 50 steroid inhibitors [Numazawa, M. et al., J. Med. Chem., 37 (1994) 2198, and references cited therein] assayed against androstenedione on human placental microsomes. An initial CoMFA resulted in a three-component model for log(l/Ki), with an explained variance r2 of 0.885, and a cross-validated q2 of 0.673. Chemometric studies were performed using GOLPE [Baroni, M. et al., Quant. Struct.-Act. Relatsh., 12 (1993) 9]. The CoMFA/GOLPE model is discussed in terms of robustness, predictivity, explanatory power and simplicity. After randomized exclusion of 25 or 10 compounds (repeated 25 times), the q2 for one component was 0.62 and 0.61, respectively, while r2 was 0.674. We demonstrate that the predictive r2 based on the mean activity (Ym) of the training set is misleading, while the test set Ym-based predictive r2 index gives a more accurate estimate of external predictivity. Using CoMFA, the observed differences in aromatase inhibition among C6-substituted steroids are rationalized at the atomic level. The CoMFA fields are consistent with known, potent inhibitors of aromatase, not included in the model. When positioned in the same alignment, these compounds have distinct features that overlap with the steric and electrostatic fields obtained in the CoMFA model. The presence of two hydrophobic binding pockets near the aromatase active site is discussed: a steric bulk tolerant one, common for C4, C6-alpha and C7-alpha substituents, and a smaller one at the C6-beta region.

A statistical mechanical description of biomolecular hydration

An efficient and accurate theoretical description of the structural hydration of biological macromolecules is presented. The hydration of molecules of almost arbitrary size (tRNA, antibody-antigen complexes, photosynthetic reaction centre) can be studied in solution and in the crystalline environment. The biomolecular structure obtained from X-ray crystallography, NMR or modelling is required as input information. The structural arrangement of water molecules near a biomolecular surface is represented by the local water density, analogous to the corresponding electron density in an X-ray diffraction experiment. The water-density distribution is approximated in terms of two- and three-particle correlation functions of solute atoms with water using a potentials-of-mean-force expansion.

Simulations of CRP:(cAMP)2 in noncrystalline environments show a subunit transition from the open to the closed conformation

The CRP:cAMP complex functions as a transcription factor that facilitates RNA polymerase recognition of several bacterial promoters. Detailed crystal structure information is available for CRP:(cAMP)2 and for CRP:(cAMP)2 complexed with DNA. In the crystalline environment, CRP:(cAMP)2 subunits are asymmetrically related; one subunit has a closed conformation and the other has an open conformation. The CRP:(cAMP)2 complexed with DNA shows both subunits in a closed conformation. We have studied the molecular dynamics of CRP:(cAMP)2 in noncrystalline environments. CRP:(cAMP)2 was simulated for 625 ps in vacuo and for 140 ps in solution. The crystal structure of CRP:(cAMP)2 in the absence of DNA was used as the initial conformation. Molecule optimal dynamic coordinates (MODCs) (García A, 1992, Phys Rev Lett 68:2696) were used to analyze protein conformations sampled during the course of the simulations. Two MODCs define a transition of the open subunit to a closed subunit conformation during the first 125 ps of simulation in vacuo; the resulting subunit conformation is similar to that observed in CRP:(cAMP)2:DNA crystals. Simulation of CRP:(cAMP)2 in solution showed that a transition from the open to the closed state also occurs when water is explicitly included in the calculations. These calculations suggest that the asymmetric conformation of CRP:(cAMP)2 is stabilized by crystal lattice interactions. The predicted solution conformation is more symmetric, with both subunits in a closed conformation.

Hydration of nucleic acid fragments: comparison of theory and experiment for high-resolution crystal structures of RNA, DNA, and DNA-drug complexes

A computationally efficient method to describe the organization of water around solvated biomolecules is presented. It is based on a statistical mechanical expression for the water-density distribution in terms of particle correlation functions. The method is applied to analyze the hydration of small nucleic acid molecules in the crystal environment, for which high-resolution x-ray crystal structures have been reported. Results for RNA [r(ApU).r(ApU)] and DNA [d(CpG).d(CpG) in Z form and with parallel strand orientation] and for DNA-drug complexes [d(CpG).d(CpG) with the drug proflavine intercalated] are described. A detailed comparison of theoretical and experimental data shows positional agreement for the experimentally observed water sites. The presented method can be used for refinement of the water structure in x-ray crystallography, hydration analysis of nuclear magnetic resonance structures, and theoretical modeling of biological macromolecules such as molecular docking studies. The speed of the computations allows hydration analyses of molecules of almost arbitrary size (tRNA, protein-nucleic acid complexes, etc.) in the crystal environment and in aqueous solution.

Sampling of the conformations of the d(CGCTGCGGC) hairpin in solution by two-dimensional nuclear magnetic resonance and theoretical methods

Most NMR studies of DNA oligomers have focused on rigid structures that show a strong preference for one or a small set of ground-state conformations. There is an increasing interest in extending NMR methods to investigate DNA systems in which this preference does not exist. A DNA hairpin is one such system where a large number of low-energy structures coexist in solution. In this article we show how 1D/2D NMR data of the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin are used to map the conformational space of this molecule. First, we characterize the gross morphology of the hairpin by monitoring the exchangeable imino signals in the molecule. Second, we extract a set of inter-proton distances (i.e., the average values and the associated dispersions) for various pairwise interactions by performing full-matrix NOESY simulation with respect to the observed NOESY data for mixing times of 250 and 100 ms. Third, we use these distances as structural constraints to perform a 300-ps molecular dynamics simulation at 500 K. Fourth, we extract 600 snapshots (one after every 0.5 ps) from the MD trajectory and perform constrained energy minimization to map local minima on the sampled energy surface (we call this the rapid temperature quenching step). Fifth, we assign 600 structures to 14 disjoint clusters such that conformationally similar hairpins belong to the same cluster while conformationally distinct hairpins belong to different clusters. Finally, we interpret the NOESY data in terms of conformationally distinct structures by recalculating NOESY contributions taken from representative structures of different clusters. Our analyses clearly demonstrate that the NMR data correspond to an ensemble of distinct structures, i.e., a set of energetically stable but conformationally distinct structures that satisfies the constraints of loop folding in the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin. Two types of loop folding consistent with NMR data are obtained: (i) a hairpin with two G.C pairs in the stem and four residues in the loop and (ii) a hairpin with two G.C pairs and a reverse wobble G.T pair in the stem plus two residues in the loop.

Harmonic vibrations and thermodynamic stability of a DNA oligomer in monovalent salt solution

We compute the full harmonic vibrational spectrum and eigenmodes of a DNA oligomer, d(C-G)3, in optimized B and Z conformations in various ionic environments (0.01-5.0 M NaCl). The statistical interactions of DNA with the diffuse ionic cloud surrounding it in solution are approximately represented within the potential of mean force framework. The lowest eigenfrequency of the B conformation is found to drastically decrease with increased NaCl concentration. This suggests that a soft mode mechanism may be a precursor for the B-to-Z conversion. The free energy balance governing the B-Z isomerization of d(C-G)3 is dominated by the solvent-averaged effective phosphate-phosphate interactions due to substantial cancelations between the much larger intramolecular energy contributions.