To Switch or Not To Switch: The Effects of Potassium and Sodium Ions on α-Poly-l-glutamate Conformations in Aqueous Solutions

K+ promotes the favorable effect of polyamine on gene expression better than Na

BACKGROUND

Polyamines are involved in a wide variety of biological processes including a marked effect on the structure and function of DNA. During our study on the interaction of polyamines with DNA, we found that K+ enhanced in vitro gene expression in the presence of polyamine more strongly than Na+. Thus, we sought to clarify the physico-chemical mechanism underlying this marked difference between the effects of K+ and Na+.

PRINCIPAL FINDINGS

It was found that K+ enhanced gene expression in the presence of spermidine, SPD(3+), much more strongly than Na+, through in vitro experiments with a Luciferase assay on cell extracts. Single-DNA observation by fluorescence microscopy showed that Na+ prevents the folding transition of DNA into a compact state more strongly than K+. 1H NMR measurement revealed that Na+ inhibits the binding of SPD to DNA more strongly than K+. Thus, SPD binds to DNA more favorably in K+-rich medium than in Na+-rich medium, which leads to favorable conditions for RNA polymerase to access DNA by decreasing the negative charge.

CONCLUSION AND SIGNIFICANCE

We found that Na+ and K+ exhibit markedly different effects through competitive binding with a cationic polyamine, SPD, to DNA, which causes a large difference in the higher-order structure of genomic DNA. It is concluded that the larger favorable effect of Na+ than K+ on in vitro gene expression observed in this study is well attributable to the significant difference between Na+ and K+ on the competitive binding inducing conformational transition of DNA.

EDTA aggregates induce SYPRO orange-based fluorescence in thermal shift assay

Ethylenediaminetetraacetic acid (EDTA) is widely used in the life sciences as chelating ligand of metal ions. However, formation of supramolecular EDTA aggregates at pH > 8 has been reported, which may lead to artifactual assay results. When applied as a buffer component at pH ≈ 10 in differential scanning fluorimetry (TSA) using SYPRO Orange as fluorescent dye, we observed a sharp change in fluorescence intensity about 20°C lower than expected for the investigated protein. We hypothesized that this change results from SYPRO Orange/EDTA interactions. TSA experiments in the presence of SYPRO Orange using solutions that contain EDTA-Na+ but no protein were performed. The TSA experiments provide evidence that suggests that at pH > 9, EDTA4- interacts with SYPRO Orange in a temperature-dependent manner, leading to a fluorescence signal yielding a "denaturation temperature" of ~68°C. Titrating Ca2+ to SYPRO Orange and EDTA solutions quenched fluorescence. Ethylene glycol tetraacetic acid (EGTA) behaved similarly to EDTA. Analytical ultracentrifugation corroborated the formation of EDTA aggregates. Molecular dynamics simulations of free diffusion of EDTA-Na+ and SYPRO Orange of in total 27 μs suggested the first structural model of EDTA aggregates in which U-shaped EDTA4- arrange in an inverse bilayer-like manner, exposing ethylene moieties to the solvent, with which SYPRO Orange interacts. We conclude that EDTA aggregates induce a SYPRO Orange-based fluorescence in TSA. These results make it relevant to ascertain that future TSA results are not influenced by interference between EDTA, or EDTA-related molecules, and the fluorescent dye.

Beyond the Hofmeister Series: Ion-Specific Effects on Proteins and Their Biological Functions

Ions differ in their ability to salt out proteins from solution as expressed in the lyotropic or Hofmeister series of cations and anions. Since its first formulation in 1888, this series has been invoked in a plethora of effects, going beyond the original salting out/salting in idea to include enzyme activities and the crystallization of proteins, as well as to processes not involving proteins like ion exchange, the surface tension of electrolytes, or bubble coalescence. Although it has been clear that the Hofmeister series is intimately connected to ion hydration in homogeneous and heterogeneous environments and to ion pairing, its molecular origin has not been fully understood. This situation could have been summarized as follows: Many chemists used the Hofmeister series as a mantra to put a label on ion-specific behavior in various environments, rather than to reach a molecular level understanding and, consequently, an ability to predict a particular effect of a given salt ion on proteins in solutions. In this Feature Article we show that the cationic and anionic Hofmeister series can now be rationalized primarily in terms of specific interactions of salt ions with the backbone and charged side chain groups at the protein surface in solution. At the same time, we demonstrate the limitations of separating Hofmeister effects into independent cationic and anionic contributions due to the electroneutrality condition, as well as specific ion pairing, leading to interactions of ions of opposite polarity. Finally, we outline the route beyond Hofmeister chemistry in the direction of understanding specific roles of ions in various biological functionalities, where generic Hofmeister-type interactions can be complemented or even overruled by particular steric arrangements in various ion binding sites.

Quantification of cation–anion interactions in crystalline monopotassium and monosodium glutamate salts

Specific anion–cation orbital interactions lead to the large structural and spectral differences observed in crystalline monosodium and monopotassium glutamates.

Improving the description of interactions between Ca2+ and protein carboxylate groups, including γ-carboxyglutamic acid: revised CHARMM22* parameters

We show that the CHARMM22* force-field over-binds the interaction between aqueous carboxylates and Ca²⁺, and introduce a modification that can recover experimentally-determined binding free energies for these systems.

Structural conservation of the short α-helix in modified higher and lower polarity water solutions

Which conformation is preferred when the polarity of water molecules is scaled byE_P=E_LJ+S²E_C?

Strong non-classical induction forces in ion-surface interactions: general origin of Hofmeister effects

Hofmeister effects continue to defy all-encompassing theories and their origin is still a matter of debate. We observed strong Hofmeister effects in Ca²⁺/Na⁺ exchange on a permanently charged surface over a wide range of ionic strengths. They could not be attributed to dispersion forces, classical induction forces, ionic size, or hydration effects. We demonstrated that another stronger force was active in the ion-surface interactions and which would create Hofmeister effects in general. The strength of this force was up to 10⁴ times that of the classical induction force and could be comparable to the Coulomb force. Coulomb, dispersion and hydration effects appeared to be interwined to affect the force. The presence of the observed strong non-classical induction force implied that energies of non-valence electrons of ions/atoms at the interface might be heavily underestimated in current theories and possibly just those underestimated energies of non-valence electrons determined Hofmeister effects.

Sodium versus potassium effects on the glutamic acid side-chains interaction on a heptapeptide

Equilibrium peptide conformations in solution, especially in the presence of salts, has been of interest for several decades. The fundamental interactions that determine the dominant peptide conformations in solution have been experimentally and computationally probed; however, a unified understanding has not yet emerged. In a previous study, we performed metadynamics simulations on the heptapeptide AEAAAEA in Sodium Chloride ( NaCl ) and Potassium Chloride ( KCl ) solutions at concentrations ranging from 0.5–2.0 M. Using a three-dimensional collective variable coordinate system, we computed the free energy landscapes in each saline environment as well as in pure water. We found that the presence of Na + and K + ions induces some changes in the stability of the conformers that define the state space, but does not alter the overall energetics between conformers and does not favor helical conformations. We investigate here, how the presence of salts ( NaCl and KCl ) affects the glutamic–glutamic interaction and its consequences on the stability of each equilibrium conformation. We perform this study through fixed backbone simulations for the most populated conformations identified in our previous work: the α-helix, 310-helix, π-helix, the extended polyproline II (PPII) and 2.51-helix conformations. It was found that for each conformation, there exists stable substates determined by the glutamic acid side-chains distance and orientation, and that Na + and K + cations (de)stabilize preferentially each conformation. It was also found that intramolecular single water mediated hydrogen bonds play a crucial role in the observed (de) stabilization of each equilibrium conformation.

Mixed-salt effects on the conformation of a short salt-bridge-forming α helix: a simulation study

The structure of a single alanine-based ACE-AEAAAKEAAAKA-NH2 peptide in explicit aqueous solutions with mixed inorganic salts (NaCl and KCl) is investigated by using molecular simulations. The concentration of Na(+), c(Na(+)), varies from 0.0M to 1.0M, whereas the concentration of K(+) is 1-c(Na(+)). The simulated peptide is very sensitive to the change of concentration ratio between Na(+) and K(+). When the concentration ratio between Na^{+} and K^{+} is changed from 0.5/0.5, the structure of the peptide becomes loose or disordered. This specific phenomenon is confirmed via checking the changes of helix parameters and mapping the free energy along different coordinates. The higher normalized probability of forming direct and indirect salt bridges between residues Glu7(+) and Lys11(+) and the smallest probability of forming ringlike structures should be responsible for the stabilized helix structure in the 0.5 Na(+)/0.5 K(+) solution. Furthermore, a noticeable conformational transition from an extended helix to an α helix is found in the 0.5 Na(+)/0.5 K(+) solution, where a local ion cloud shows that some Na(+) ions in the inner shells are still directly binding with the peptide, while K(+) in the outer shells are moving into the inner shells, keeping the peptide in the collapsed state.

Salt effects on hydrophobic-core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica-exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt-bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near-native conformations in pure water. Our analyses also show that such salt-bridge formation hinders the fully native hydrophobic-core packing at the final stages of folding.

A systematic investigation and insight into the formation mechanism of bilayers of fatty acid/soap mixtures in aqueous solutions

Vesicles are the most common form of bilayer structures in fatty acid/soap mixtures in aqueous solutions; however, a peculiar bilayer structure called a "planar sheet" was found for the first time in the mixtures. In the past few decades, considerable research has focused on the formation theory of bilayers in fatty acid/soap mixtures. The hydrogen bond theory has been widely accepted by scientists to explain the formation of bilayers. However, except for the hydrogen bond, no other driving forces were proposed systematically. In this work, three kinds of weak interactions were investigated in detail, which could perfectly demonstrate the formation mechanism of bilayer structures in the fatty acid/soap mixtures in aqueous solutions. (i) The influence of hydrophobic interaction was detected by changing the chain length of fatty acid (C(n)H(2n+1)COOH), in which n = 10 to 18, the phase behavior was investigated, and the phase region was presented. With the help of cryogenic transmission electron microscopy (cryo-TEM) observations, deuterium nuclear magnetic resonance ((2)H NMR), and X-ray diffraction (XRD) measurements, the vesicles and planar sheets were determined. The chain length of C(n)H(2n+1)COOH has an important effect on the physical state of the hydrophobic chain, resulting in an obvious difference in the viscoelasticity of the solution samples. (ii) The existence of hydrogen bonds between fatty acids and their soaps in aqueous solutions was demonstrated by Fourier transform infrared (FT-IR) spectroscopy and molecule dynamical simulation. From the pH measurements, the pH ranges of the bilayer formation were at the pKa values of fatty acids, respectively. (iii) Counterions can be embedded in the stern layer of the bilayers and screen the electrostatic repulsion between the COO(-) anionic headgroups. FT-IR characterization demonstrated a bidentate bridging coordination mode between counterions and carboxylates. The conductivity measurements provided the degree of counterion binding (β = 0.854), indicating the importance of the counterions.

Salting-in with a salting-out agent: explaining the cation specific effects on the aqueous solubility of amino acids

Although the understanding of ion specific effects on the aqueous solubilities of biomolecules is crucial for the development of many areas of biochemistry and life sciences, a consensual and well-supported molecular picture of the phenomena has not yet been established. Mostly, the influence of cations and the nature of the molecular interactions responsible for the reversal of the Hofmeister trend in aqueous solutions of amino acids and proteins are still defectively understood. Aiming at contributing to the understanding of the molecular-level mechanisms governing the cation specific effects on the aqueous solubilities of biocompounds, experimental solubility measurements and classical molecular dynamics simulations were performed for aqueous solutions of three amino acids (alanine, valine, and isoleucine), in the presence of a series of inorganic salts. The evidence gathered suggests that the mechanism by which salting-in inducing cations operate in aqueous solutions of amino acids is different from that of anions, and allows for a novel and consistent molecular description of the effect of the cation on the solubility based on specific interactions of the cations with the negatively charged moieties of the biomolecules.

Conformational free-energy landscapes for a peptide in saline environments

The conformations that proteins adopt in solution are a function of both their primary structure and surrounding aqueous environment. Recent experimental and computational work on small peptides, e.g., polyK, polyE, and polyR, have highlighted an interesting and unusual behavior in the presence of aqueous ions such as ClO₄⁻, Na⁺, and K⁺. Notwithstanding the aforementioned studies, as of this writing, the nature of the driving force induced by the presence of ions and its role on the conformational stability of peptides remains only partially understood. Molecular-dynamics simulations have been performed on the heptapeptide AEAAAEA in NaCl and KCl solutions at concentrations of 0.5, 1.0, and 2.0 M. Metadynamics in conjunction with a three-dimensional model reaction coordinate was used to sample the conformational space of the peptide. All simulations were run for 2 μs. Free-energy landscapes were computed over the model reaction coordinate for the peptide in each saline assay as well as in the absence of ions. Circular dichroism spectra were also calculated from each trajectory. In the presence of Na⁺ and K⁺ ions, no increase in helicity is observed with respect to the conformation in pure water.

Salt-induced control of supramolecular order in biocatalytic hydrogelation

Biocatalytic action and specific ion effects are both known to have dramatic effects on molecular self-assembly and hydrogelation. In this paper, we demonstrate that these effects are highly cooperative. Biocatalytic hydrogelation of Fmoc peptides in the presence of salts combines kinetic (through enzymatic catalysis) and thermodynamic (specific ion and protein templating) contributions when applied in combination. Spectroscopic data (obtained by fluorescence spectroscopy and circular dichroism) revealed that hydrophobic interactions are greatly affected, giving rise to differential chiral organization and supramolecular structure formation. The kinetic effects of catalytic action could be removed from the system by applying a heat/cool cycle, giving insight into the thermodynamic influence of both protein and salt on these systems and showing that the effects of catalysis, templating, and salts are cooperative. The variable molecular interactions are expressed as variable material properties, such as thermal stability and mechanical strength of the final gel-phase material. To gain more insight into the role of the enzyme, beyond catalysis, in the underlying mechanism, static light scattering is performed, which indicates the different mode of aggregation of the enzyme molecules in the presence of different salts in aqueous solution that may play a role to direct the assembly via templating. Overall, the results show that the combination of specific salts and enzymatic hydrogelation can give rise to complex self-assembly behaviors that may be exploited to tune hydrogel properties.

Hydrogen-bond stabilization in oxyanion holes: grand jeté to three dimensions

We recently reported crystallographic evidence that the hydrogen bonds which can stabilize oxygen-centered negative charge within enzyme oxyanion holes are rarely found in the place they should be expected on the basis of the analysis of small-molecule crystal structures. We investigated this phenomenon using calculations on simplified active site models. A recent paper suggested that several aspects of the analysis required further exploration. In this paper we: (i) review the results of our crystallographic study; (ii) report molecular dynamics studies which investigate the effect of protein movement; (iii) report ONIOM calculations which trace the reaction coordinate for an oxyanion hole reaction in the presence of a complete enzyme active site. These results show that the limitations of gas phase calculations on simplified models do not invalidate our comparison of competing active site geometries. These new results reaffirm the conclusion that oxyanion holes are not usually stabilized by planar arrangements of H-bonds, and that this sub-optimal transition state stabilization leads to better overall catalysis.

Conformation of poly-L-glutamate is independent of ionic strength

CD and UV resonance Raman measurements surprisingly find that the charge screening of even 2 M concentrations of NaCl and KCl does not alter the unfolded PPII and 2.5(1)-helix conformations of poly-L-glutamate. These salts appear to be excluded from the region between the side chain charges and the peptide backbone. Furthermore, no direct ion pairing occurs between these salts and the side chain carboxylates.

Differential geometry based solvation model II: Lagrangian formulation

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation models. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The optimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and PB equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of computation, thanks to the equivalence of the Laplace-Beltrami operator in the two representations. The coupled partial differential equations (PDEs) are solved with an iterative procedure to reach a steady state, which delivers desired solvent-solute interface and electrostatic potential for problems of interest. These quantities are utilized to evaluate the solvation free energies and protein-protein binding affinities. A number of computational methods and algorithms are described for the interconversion of Lagrangian and Eulerian representations, and for the solution of the coupled PDE system. The proposed approaches have been extensively validated. We also verify that the mean curvature flow indeed gives rise to the minimal molecular surface and the proposed variational procedure indeed offers minimal total free energy. Solvation analysis and applications are considered for a set of 17 small compounds and a set of 23 proteins. The salt effect on protein-protein binding affinity is investigated with two protein complexes by using the present model. Numerical results are compared to the experimental measurements and to those obtained by using other theoretical methods in the literature.

Energy landscape analyses of disordered histone tails reveal special organization of their conformational dynamics

Histone tails are highly flexible N- or C-terminal protrusions of histone proteins which facilitate the compaction of DNA into dense superstructures known as chromatin. On a molecular scale histone tails are polyelectrolytes with high degree of conformational disorder which allows them to function as biomolecular "switches", regulating various genetic processes. Unfortunately, their intrinsically disordered nature creates obstacles for comprehensive experimental investigation of both the structural and dynamical aspects of histone tails, because of which their conformational behaviors are still not well understood. In this work we have carried out ∼3 microsecond long all atom replica exchange molecular dynamics (REMD) simulations for each of four histone tails, H4, H3, H2B, and H2A, and probed their intrinsic conformational preferences. Our subsequent free energy landscape analysis demonstrated that most tails are not fully disordered, but show distinct conformational organization, containing specific flickering secondary structural elements. In particular, H4 forms β-hairpins, H3 and H2B adopt α-helical elements, while H2A is fully disordered. We rationalized observed patterns of conformational dynamics of various histone tails using ideas from physics of polyelectrolytes and disordered systems. We also discovered an intriguing re-entrant contraction-expansion of the tails upon heating, which is caused by subtle interplay between ionic screening and chain entropy.