Calcium is one of the biologically most important ions; however, its accurate description by classical molecular dynamics simulations is complicated by strong electrostatic and polarization ...interactions with surroundings due to its divalent nature. Here, we explore the recently suggested approach for effectively accounting for polarization effects via ionic charge rescaling and develop a new and accurate parametrization of the calcium dication. Comparison to neutron scattering and viscosity measurements demonstrates that our model allows for an accurate description of concentrated aqueous calcium chloride solutions. The present model should find broad use in efficient and accurate modeling of calcium in aqueous environments, such as those encountered in biological and technological applications.
It is a textbook knowledge that charges of the same polarity repel each other. For two monovalent ions in the gas phase at a close contact this repulsive interaction amounts to hundreds of kilojoules ...per mole. In aqueous solutions, however, this Coulomb repulsion is strongly attenuated by a factor equal to the dielectric constant of the medium. The residual repulsion, which now amounts only to units of kilojoules per mole, may be in principle offset by attractive interactions. Probably the smallest cationic pair, where a combination of dispersion and cavitation forces overwhelms the Coulomb repulsion, consists of two guanidinium ions in water. Indeed, by a combination of molecular dynamics with electronic structure calculations and electrophoretic, as well as spectroscopic, experiments, we have demonstrated that aqueous guanidinium cations form (weakly) thermodynamically stable like-charge ion pairs. The importance of pairing of guanidinium cations in aqueous solutions goes beyond a mere physical curiosity, since it has significant biochemical implications. Guanidinium chloride is known to be an efficient and flexible protein denaturant. This is due to the ability of the orientationally amphiphilic guanidinium cations to disrupt various secondary structural motifs of proteins by pairing promiscuously with both hydrophobic and hydrophilic groups, including guanidinium-containing side chains of arginines. The fact that the cationic guanidinium moiety forms the dominant part of the arginine side chain implies that the like-charge ion pairing may also play a role for interactions between peptides and proteins. Indeed, arginine–arginine pairing has been frequently found in structural protein databases. In particular, when strengthened by a presence of negatively charged glutamate, aspartate, or C-terminal carboxylic groups, this binding motif helps to stabilize peptide or protein dimers and is also found in or near active sites of several enzymes. The like-charge pairing of the guanidinium side-chain groups may also hold the key to the understanding of the arginine “magic”, that is, the extraordinary ability of arginine-rich polypeptides to passively penetrate across cellular membranes. Unlike polylysines, which are also highly cationic but lack the ease in crossing membranes, polyarginines do not exhibit mutual repulsion. Instead, they accumulate at the membrane, weaken it, and might eventually cross in a concerted, “train-like” manner. This behavior of arginine-rich cell penetrating peptides can be exploited when devising smart strategies how to deliver in a targeted way molecular cargos into the cell.
Arginine-rich cell-penetrating peptides do not enter cells by directly passing through a lipid membrane; they instead passively enter vesicles and live cells by inducing membrane multilamellarity and ...fusion. The molecular picture of this penetration mode, which differs qualitatively from the previously proposed direct mechanism, is provided by molecular dynamics simulations. The kinetics of vesicle agglomeration and fusion by an iconic cellpenetrating peptide—nonaarginine—are documented via real-time fluorescence techniques, while the induction of multilamellar phases in vesicles and live cells is demonstrated by a combination of electron and fluorescence microscopies. This concert of experiments and simulations reveals that the identified passive cell penetration mechanism bears analogy to vesicle fusion induced by calcium ions, indicating that the two processes may share a common mechanistic origin.
Magnesium and zinc dications possess the same charge and have an almost identical size, yet they behave very differently in aqueous solutions and play distinct biological roles. It is thus crucial to ...identify the origins of such different behaviors and to assess to what extent they can be captured by force-field molecular dynamics simulations. In this work, we combine neutron scattering experiments in a specific mixture of H2O and D2O (the so-called null water) with ab initio molecular dynamics simulations to probe the difference in the hydration structure and ion-pairing properties of chloride solutions of the two cations. The obtained data are used as a benchmark to develop a scaled-charge force field for Mg2+ that includes electronic polarization in a mean field way. We show that using this electronic continuum correction we can describe aqueous magnesium chloride solutions well. However, in aqueous zinc chloride specific interaction terms between the ions need to be introduced to capture ion pairing quantitatively.
The transcription factor RUNX1 is frequently mutated in myelodysplastic syndrome and leukemia. RUNX1 mutations can be early events, creating preleukemic stem cells that expand in the bone marrow. ...Here we show, counterintuitively, that Runx1-deficient hematopoietic stem and progenitor cells (HSPCs) have a slow growth, low biosynthetic, small cell phenotype and markedly reduced ribosome biogenesis (Ribi). The reduced Ribi involved decreased levels of rRNA and many mRNAs encoding ribosome proteins. Runx1 appears to directly regulate Ribi; Runx1 is enriched on the promoters of genes encoding ribosome proteins and binds the rDNA repeats. Runx1-deficient HSPCs have lower p53 levels, reduced apoptosis, an attenuated unfolded protein response, and accordingly are resistant to genotoxic and ER stress. The low biosynthetic activity and corresponding stress resistance provides a selective advantage to Runx1-deficient HSPCs, allowing them to expand in the bone marrow and outcompete normal HSPCs.
Display omitted
•Runx1-deficient HSCs are resistant to endogenous and genotoxic stress•Runx1-deficient HSCs have decreased ribosome biogenesis•Runx1 directly regulates ribosome biogenesis
Loss-of-function RUNX1 mutations are common in myelodysplastic syndrome and leukemia and can be early events. Cai et al. demonstrate that Runx1 loss decreases ribosome biogenesis and translation in hematopoietic stem and progenitor cells and confers resistance to endogenous and genotoxic stress.
Molecular dynamics simulations of concentrated aqueous solutions of LiCl and Li2SO4 were conducted in order to provide molecular insight into recent neutron scattering data. The structures predicted ...from the molecular dynamics simulations using standard nonpolarizable force fields provided a very poor fit to the experiment; therefore, refinement was needed. The electronic polarizability of the medium was effectively accounted for by implementing the electronic continuum correction, which practically means rescaling the ionic charges. Consistent with previous studies, we found that this approach in each case provided a significantly improved fit to the experimental data, which was further enhanced by slightly adjusting the radius of the lithium ion. The polarization effect was particularly pronounced in the Li2SO4 solution where the ions in the nonpolarizable simulations tended to cluster unphysically. With the above alterations, the employed force field displayed an excellent fit to the neutron scattering data and provided a useful interpretative framework for the experimental measurements. At the same time, the present study underlines the importance of solvent polarization effects in hydration of ions with high charge density.
We employed density functional theory-based ab initio molecular dynamics simulations to examine the hydration structure of several common alkali and alkali earth metal cations. We found that the ...commonly used atom pairwise dispersion correction scheme D3, which assigns dispersion coefficients based on the neutral form of the atom rather than its actual oxidation state, leads to inaccuracies in the hydration structures of these cations. We evaluated this effect for lithium, sodium, potassium, and calcium and found that the inaccuracies are particularly pronounced for sodium and potassium compared to the experiment. To remedy this issue, we propose disabling the D3 correction specifically for all cation-including pairs, which leads to a much better agreement with experimental data.
Gating charges in voltage-sensing domains (VSD) of voltage-sensitive ion channels and enzymes are carried on arginine side chains rather than lysine. This arginine preference may result from the ...unique hydration properties of the side chain guanidinium group which facilitates its movement through a hydrophobic plug that seals the center of the VSD, as suggested by molecular dynamics simulations. To test for side chain interactions implicit in this model we inspected interactions of the side chains of arginine and lysine with each of the 19 non-glycine amino acids in proteins in the protein data bank. The arginine guanidinium interacts with non-polar aromatic and aliphatic side chains above and below the guanidinium plane while hydrogen bonding with polar side chains is restricted to in-plane positions. In contrast, non-polar side chains interact largely with the aliphatic part of the lysine side chain. The hydration properties of arginine and lysine are strongly reflected in their respective interactions with non-polar and polar side chains as observed in protein structures and in molecular dynamics simulations, and likely underlie the preference for arginine as a mobile charge carrier in VSD.
The carbonate ion plays a central role in the biochemical formation of the shells of aquatic life, which is an important path for carbon dioxide sequestration. Given the vital role of carbonate in ...this and other contexts, it is imperative to develop accurate models for such a high charge density ion. As a divalent ion, carbonate has a strong polarizing effect on surrounding water molecules. This raises the question whether it is possible to describe accurately such systems without including polarization. It has recently been suggested the lack of electronic polarization in nonpolarizable water models can be effectively compensated by introducing an electronic dielectric continuum, which is with respect to the forces between atoms equivalent to rescaling the ionic charges. Given how widely nonpolarizable models are used to model electrolyte solutions, establishing the experimental validity of this suggestion is imperative. Here, we examine a stringent test for such models: a comparison of the difference of the neutron scattering structure factors of K2CO3 vs KNO3 solutions and that predicted by molecular dynamics simulations for various models of the same systems. We compare standard nonpolarizable simulations in SPC/E water to analogous simulations with effective ion charges, as well as simulations in explicitly polarizable POL3 water (which, however, has only about half the experimental polarizability). It is found that the simulation with rescaled charges is in a very good agreement with the experimental data, which is significantly better than for the nonpolarizable simulation and even better than for the explicitly polarizable POL3 model.
PDF
