Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches

Since no single experimental or modeling technique provides data that allow a description of transport processes in clays and clay minerals at all relevant scales, several complementary approaches have to be combined to understand and explain the interplay between transport relevant phenomena. In this paper molecular dynamics simulations (MD) were used to investigate the mobility of water in the interlayer of montmorillonite (Mt), and to estimate the influence of mineral surfaces and interlayer ions on the water diffusion. Random Walk (RW) simulations based on a simplified representation of pore space in Mt were used to estimate and understand the effect of the arrangement of Mt particles on the meso- to macroscopic diffusivity of water. These theoretical calculations were complemented with quasielastic neutron scattering (QENS) measurements of aqueous diffusion in Mt with two pseudo-layers of water performed at four significantly different energy resolutions (i.e. observation times). The size of the interlayer and the size of Mt particles are two characteristic dimensions which determine the time dependent behavior of water diffusion in Mt. MD simulations show that at very short time scales water dynamics has the characteristic features of an oscillatory motion in the cage formed by neighbors in the first coordination shell. At longer time scales, the interaction of water with the surface determines the water dynamics, and the effect of confinement on the overall water mobility within the interlayer becomes evident. At time scales corresponding to an average water displacement equivalent to the average size of Mt particles, the effects of tortuosity are observed in the meso- to macroscopic pore scale simulations. Consistent with the picture obtained in the simulations, the QENS data can be described using a (local) 3D diffusion at short observation times, whereas at sufficiently long observation times a 2D diffusive motion is clearly observed. The effects of tortuosity measured in macroscopic tracer diffusion experiments are in qualitative agreement with RW simulations. By using experimental data to calibrate molecular and mesoscopic theoretical models, a consistent description of water mobility in clay minerals from the molecular to the macroscopic scale can be achieved. In turn, simulations help in choosing optimal conditions for the experimental measurements and the data interpretation. (C) 2014 Elsevier B.V. All rights reserved.

The Hydroxyl Side Chain of a Highly Conserved Serine Residue Is Required for Cation Selectivity and Substrate Transport in the Glial Glutamate Transporter GLT-1/SLC1A2

Glutamate transporters maintain synaptic concentration of the excitatory neurotransmitter below neurotoxic levels. Their transport cycle consists of cotransport of glutamate with three sodium ions and one proton, followed by countertransport of potassium. Structural studies proposed that a highly conserved serine located in the binding pocket of the homologous GltPh coordinates l-aspartate as well as the sodium ion Na1. To experimentally validate these findings, we generated and characterized several mutants of the corresponding serine residue, Ser-364, of human glutamate transporter SLC1A2 (solute carrier family 1 member 2), also known as glutamate transporter GLT-1 and excitatory amino acid transporter EAAT2. S364T, S364A, S364C, S364N, and S364D were expressed in HEK cells and Xenopus laevis oocytes to measure radioactive substrate transport and transport currents, respectively. All mutants exhibited similar plasma membrane expression when compared with WT SLC1A2, but substitutions of serine by aspartate or asparagine completely abolished substrate transport. On the other hand, the threonine mutant, which is a more conservative mutation, exhibited similar substrate selectivity, substrate and sodium affinities as WT but a lower selectivity for Na(+) over Li(+). S364A and S364C exhibited drastically reduced affinities for each substrate and enhanced selectivity for l-aspartate over d-aspartate and l-glutamate, and lost their selectivity for Na(+) over Li(+). Furthermore, we extended the analysis of our experimental observations using molecular dynamics simulations. Altogether, our findings confirm a pivotal role of the serine 364, and more precisely its hydroxyl group, in coupling sodium and substrate fluxes.

IGHD II: A Novel GH-1 Gene Mutation (GH-L76P) Severely Affects GH Folding, Stability, and Secretion

CONTEXT The autosomal dominant form of GH deficiency (IGHD II) is characterized by markedly reduced GH secretion combined with low concentrations of IGF-1 leading to short stature. OBJECTIVE Structure-function analysis of a missense mutation in the GH-1 gene converting codon 76 from leucine (L) to proline (P) yielding a mutant GH-L76P peptide. DESIGN, SETTINGS, AND PATIENTS Heterozygosity for GH-L76P/wt-GH was identified in a nonconsanguineous Spanish family. The index patients, two siblings, a boy and a girl, were referred for assessment of their short stature (-3.2 and -3.8 SD). Their grandmother, father, and aunt were also carrying the same mutation and showed severe short stature; therefore, IGHD II was diagnosed. INTERVENTIONS AND RESULTS AtT-20 cells coexpressing both wt-GH and GH-L76P showed a reduced GH secretion (P < .001) after forskolin stimulation when compared with the cells expressing only wt-GH. In silico mutagenesis and molecular dynamics simulations presented alterations of correct folding and mutant stability compared with wt-GH. Therefore, further structural analysis of the GH-L76P mutant was performed using expressed and purified proteins in Escherichia coli by thermofluor assay and fast degradation proteolysis assay. Both assays revealed that the GH-L76P mutant is unstable and misfolded compared to wt-GH confirming the bioinformatic model prediction. CONCLUSIONS This is the first report of a family suffering from short stature caused by IGHD II, which severely affects intracellular GH folding and stability as well as secretion, highlighting the necessity of functional analysis of any GH variant for defining new mechanisms as a cause for IGHD II.

A novel computer simulation method for simulating the multiscale transduction dynamics of signal proteins

Signal proteins are able to adapt their response to a change in the environment, governing in this way a broad variety of important cellular processes in living systems. While conventional molecular-dynamics (MD) techniques can be used to explore the early signaling pathway of these protein systems at atomistic resolution, the high computational costs limit their usefulness for the elucidation of the multiscale transduction dynamics of most signaling processes, occurring on experimental timescales. To cope with the problem, we present in this paper a novel multiscale-modeling method, based on a combination of the kinetic Monte-Carlo- and MD-technique, and demonstrate its suitability for investigating the signaling behavior of the photoswitch light-oxygen-voltage-2-Jα domain from Avena Sativa (AsLOV2-Jα) and an AsLOV2-Jα-regulated photoactivable Rac1-GTPase (PA-Rac1), recently employed to control the motility of cancer cells through light stimulus. More specifically, we show that their signaling pathways begin with a residual re-arrangement and subsequent H-bond formation of amino acids near to the flavin-mononucleotide chromophore, causing a coupling between β-strands and subsequent detachment of a peripheral α-helix from the AsLOV2-domain. In the case of the PA-Rac1 system we find that this latter process induces the release of the AsLOV2-inhibitor from the switchII-activation site of the GTPase, enabling signal activation through effector-protein binding. These applications demonstrate that our approach reliably reproduces the signaling pathways of complex signal proteins, ranging from nanoseconds up to seconds at affordable computational costs.

Conformational analysis of a potent SSTR3-selective somatostatin analogue by NMR in water solution

The three-dimensional structure of a potent SSTR3-selective analogue of somatostatin, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-Agl(8)(N(beta) Me, 2-naphthoyl)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-Agl(8)(N(beta) Me, 2-naphthoyl)]-SRIF) (peptide 1) has been determined by (1)H NMR in water and molecular dynamics (MD) simulations. The peptide exists in two conformational isomers differing mainly by the cis/trans isomerization of the side chain in residue 8. The structure of 1 is compared with the consensus structural motifs of other somatostatin analogues that bind predominantly to SSTR1, SSTR2/SSTR5 and SSTR4 receptors, and to the 3D structure of a non-selective SRIF analogue, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-2Nal(8)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-2Nal(8)]-SRIF) (peptide 2). The structural determinant factors that could explain selectivity of peptide 1 for SSTR3 receptors are discussed.