Position of the Third Na+ Site in the Aspartate Transporter GltPh and the Human Glutamate Transporter, EAAT1

Glutamate transport via the human excitatory amino acid transporters is coupled to the co-transport of three Na+ ions, one H+ and the counter-transport of one K+ ion. Transport by an archaeal homologue of the human glutamate transporters, GltPh, whose three dimensional structure is known is also coupled to three Na+ ions but only two Na+ ion binding sites have been observed in the crystal structure of GltPh. In order to fully utilize the GltPh structure in functional studies of the human glutamate transporters, it is essential to understand the transport mechanism of GltPh and accurately determine the number and location of Na+ ions coupled to transport. Several sites have been proposed for the binding of a third Na+ ion from electrostatic calculations and molecular dynamics simulations. In this study, we have performed detailed free energy simulations for GltPh and reveal a new site for the third Na+ ion involving the side chains of Threonine 92, Serine 93, Asparagine 310, Aspartate 312, and the backbone of Tyrosine 89. We have also studied the transport properties of alanine mutants of the coordinating residues Threonine 92 and Serine 93 in GltPh, and the corresponding residues in a human glutamate transporter, EAAT1. The mutant transporters have reduced affinity for Na+ compared to their wild type counterparts. These results confirm that Threonine 92 and Serine 93 are involved in the coordination of the third Na+ ion in GltPh and EAAT1

Optimization of the CHARMM Additive Force Field for DNA: Improved Treatment of the BI/BII Conformational Equilibrium

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36

X-ray structure of a ClC chloride channel at 3.0 A reveals the molecular basis of anion selectivity

The ClC chloride channels catalyse the selective flow of Cl- ions across cell membranes, thereby regulating electrical excitation in skeletal muscle and the flow of salt and water across epithelial barriers. Genetic defects in ClC Cl- channels underlie several familial muscle and kidney diseases. Here we present the X-ray structures of two prokaryotic ClC Cl- channels from Salmonella enterica serovar typhimurium and Escherichia coli at 3.0 and 3.5 A, respectively. Both structures reveal two identical pores, each pore being formed by a separate subunit contained within a homodimeric membrane protein. Individual subunits are composed of two roughly repeated halves that span the membrane with opposite orientations. This antiparallel architecture defines a selectivity filter in which a Cl- ion is stabilized by electrostatic interactions with alpha-helix dipoles and by chemical coordination with nitrogen atoms and hydroxyl groups. These findings provide a structural basis for further understanding the function of ClC Cl- channels, and establish the physical and chemical basis of their anion selectivity

The Accuracy of Generalized Born Forces

Generalized Born (GB) models play an important role in biomolecular dynamics simulations. The issue of the accuracy of solvation forces computed according to these models has not been investigated in much detail to the best of our knowledge. Optimal atomic radii have been derived for most used force fields in order to reproduce molecular dynamics results under the Poisson\u2013Boltzmann (PB) continuum model. In principle, if optimized radii are used, the PB model reproduces well MD simulations. For this reason in this contribution, we consider the most used GB model implementations in publicly available molecular dynamics simulation softwares and assess the accuracy of GB solvation forces compared to the reference PB ones. The results show that current implementations that compute forces by pairwise summations correlate very well with PB forces if judicious choice of model parameters is performed. The latter can be obtained by linear fitting of solvation self-energies computed by the most accurate GBR6 GB model

Identification and expression of Acetylcholinesterase in Octopus vulgaris arm development and regeneration: A conservative role for ACHE?

Acetylcholinesterase (ACHE) is a glycoprotein with a key role in terminating synaptic transmission in cholinergic neurons of both vertebrates and invertebrates. ACHE is also involved in the regulation of cell growth and morphogenesis during embryogenesis and regeneration acting through its non-cholinergic sites. The mollusk Octopus vulgaris provides a powerful model for investigating the mechanisms underlying tissue morphogenesis due to its high regenerative power. Here, we performed a comparative investigation of arm morphogenesis during adult arm regeneration and embryonic arm development which may provide insights on the conserved ACHE pathways. In this study, we cloned and characterized O. vulgaris ACHE, finding a single highly conserved ACHE hydrophobic variant, characterized by prototypical catalytic sites and a putative consensus region for a glycosylphosphatidylinositol (GPI)-anchor attachment at the COOH-terminus. We then show that its expression level is correlated to the stage of morphogenesis in both adult and embryonic arm. In particular, ACHE is localized in typical neuronal sites when adult-like arm morphology is established and in differentiating cell locations during the early stages of arm morphogenesis. This possibility is also supported by the presence in the ACHE sequence and model structure of both cholinergic and non-cholinergic sites. This study provides insights into ACHE conserved roles during processes of arm morphogenesis. In addition, our modeling study offers a solid basis for predicting the interaction of the ACHE domains with pharmacological blockers for in vivo investigations. We therefore suggest ACHE as a target for the regulation of tissue morphogenesis.Versión del edito