From gas phase via microhydration to bulk water simulations

The salt bridge formation and stability in the terminated lysine-glutamate dipeptide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desolvated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hydrate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins

Calcium Sensing by Recoverin : Effect of Protein Conformation on Ion Affinity

The detailed functional mechanism of recoverin, which acts as a myristoyl switch at the rod outer-segment disk membrane, is elucidated by direct and replica-exchange molecular dynamics. In accord with NMR structural evidence and calcium binding assays, simulations point to the key role of enhanced calcium binding to the EF3 loop of the semiopen state of recoverin as compared to the closed state. This 2-4-order decrease in calcium dissociation constant stabilizes the semiopen state in response to the increase of cytosolic calcium concentration in the vicinity of recoverin. A second calcium ion then binds to the EF2 loop and, consequently, the structure of the protein changes from the semiopen to the open state. The latter has the myristoyl chain extruded to the cytosol, ready to act as a membrane anchor of recoverin.Peer reviewe

Scanning electron microscopy and molecular dynamics of surfaces of growing and ablating hexagonal ice crystals

We present the first clearly resolved observations of surfaces of growing and ablating hexagonal ice crystals using variable-pressure scanning electron microscopy. The ice surface develops trans-prismatic strands, separated from one another by distances of 5–10 ?m. The strands are present at a wide range of supersaturations, but are most pronounced at temperatures near the frost point. Pyramidal facets consistent with Miller-Bravais indices of 1011, and possibly also 2021, are associated with ice growth under these conditions. A molecular-dynamics model of a free-standing ice Ih nanocolumn containing 8400 water molecules does not develop trans-prismatic strands, suggesting these features originate at larger spatial or temporal scales. The possible relevance of these surface features to cirrus ice is discussed

Single Lipid Extraction: The Anchoring Strength of Cholesterol in Liquid-Ordered and Liquid-Disordered Phases

AbstractCholesterol is important for the formation of microdomains in supported lipid bilayers and is enriched in the liquid-ordered phase. To understand the interactions leading to this enrichment, we developed an AFM-based single-lipid-extraction (SLX) approach that enables us to determine the anchoring strength of cholesterol in the two phases of a phase-separated lipid membrane. As expected, the forces necessary for extracting a single cholesterol molecule from liquid-ordered phases are significantly higher than for extracting it from the liquid-disordered phases. Interestingly, application of the Bell model shows two energy barriers that correlate with the head and full length of the cholesterol molecule. The resulting lifetimes for complete extraction are 90 s and 11 s in the liquid-ordered and liquid-disordered phases, respectively. Molecular dynamics simulations of the very same experiment show similar force profiles and indicate that the stabilization of cholesterol in the liquid-ordered phase is mainly due to nonpolar contacts

Changes in the hydration structure of imidazole upon protonation: Neutron scattering and molecular simulations

International audienceThe imidazole motif is widely encountered in biomolecules, and its biological role, for instance as a proton relay, is often linked to its ability to form hydrogen bonds with water molecules. Detailed characterization of the hydration pattern of imidazole and its changes upon protonation is thus of high interest. Here, we combine neutron scattering experiments with force field simulations to provide an unprecedented characterization of the neutral and protonated imidazole solvation at the atomistic level. We show that neutron diffraction data can be used to assess the quality of the imidazole force field in molecular simulations. Simulations using the CHARMM general force field for imidazole are in excellent agreement with the experimental neutron scattering data and we use them to provide an atomic scale interpretation of the neutron scattering patterns. Upon protonation, we clearly identify the signature of the reorganization in the hydration pattern caused by the change from one H-bond donor and one H-bond acceptor group for imidazole to two H-bond donor groups for imidazolium. We also point the limits of the experiment, which is rather insensitive to details of the H-bond geometry at the deprotonated nitrogen of imidazole and further complement the description of the hydration structure with ab initio molecular dynamics simulations

Hydration and Ion Pairing in Aqueous Mg 2+ and Zn 2+ Solutions: Force-Field Description Aided by Neutron Scattering Experiments and Ab Initio Molecular Dynamics Simulations

International audienc

Membrane Binding of Recoverin : From Mechanistic Understanding to Biological Functionality

Recoverin is a neuronal calcium sensor involved in vision adaptation that reversibly associates with cellular membranes via its calcium-activated myristoyl switch. While experimental evidence shows that the myristoyl group significantly enhances membrane affinity of this protein, molecular details of the binding process are still under debate. Here, we present results of extensive molecular dynamics simulations of recoverin in the proximity of a phospholipid bilayer. We capture multiple events of spontaneous membrane insertion of the myristoyl moiety and confirm its critical role in the membrane binding. Moreover, we observe that the binding strongly depends on the conformation of the N-terminal domain. We propose that a suitable conformation of the N-terminal domain can be stabilized by the disordered C-terminal segment or by binding of the target enzyme, i.e., rhodopsin kinase. Finally, we find that the presence of negatively charged lipids in the bilayer stabilizes a physiologically functional orientation of the membrane-bound recoverin.Peer reviewe

Accurate Binding of Sodium and Calcium to a POPC Bilayer by Effective Inclusion of Electronic Polarization

Binding affinities and stoichiometries of Na+ and Ca2+ ions to phospholipid bilayers are of paramount significance in the properties and functionality of cellular membranes. Current estimates of binding affinities and stoichiometries of cations are, however, inconsistent due to limitations in the available experimental and computational methods. In this work, we improve the description of the binding details of Na+ and Ca2+ ions to a 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer by implicitly including electronic polarization as a mean field correction, known as the electronic continuum correction (ECC). This is applied by scaling the partial charges of a selected state-of-the-art POPC lipid model for molecular dynamics simulations. Our improved ECC-POPC model reproduces not only the experimentally measured structural parameters for the ion-free membrane, but also the response of lipid headgroup to a strongly bound cationic amphiphile, as well as the binding affinities of Na+ and Ca2+ ions. With our new model, we observe on the one side negligible binding of Na+ ions to POPC bilayer, while on the other side stronger interactions of Ca2+ primarily with phosphate oxygens, which is in agreement with the previous interpretations of the experimental spectroscopic data. The present model results in Ca2+ ions forming complexes with one to three POPC molecules with almost equal probabilities, suggesting more complex binding stoichiometries than those from simple models used to interpret the NMR data previously. The results of this work pave the way to quantitative molecular simulations with realistic electrostatic interactions of complex biochemical systems at cellular membranes.Peer reviewe