A two-domain elevator mechanism for sodium/proton antiport

Sodium/proton (Na+/H+) antiporters, located at the plasma membrane in every cell, are vital for cell homeostasis1. In humans, their dysfunction has been linked to diseases, such as hypertension, heart failure and epilepsy, and they are well-established drug targets2. The best understood model system for Na+/H+ antiport is NhaA from Escherichia coli1, 3, for which both electron microscopy and crystal structures are available4, 5, 6. NhaA is made up of two distinct domains: a core domain and a dimerization domain. In the NhaA crystal structure a cavity is located between the two domains, providing access to the ion-binding site from the inward-facing surface of the protein1, 4. Like many Na+/H+ antiporters, the activity of NhaA is regulated by pH, only becoming active above pH 6.5, at which point a conformational change is thought to occur7. The only reported NhaA crystal structure so far is of the low pH inactivated form4. Here we describe the active-state structure of a Na+/H+ antiporter, NapA from Thermus thermophilus, at 3 Å resolution, solved from crystals grown at pH 7.8. In the NapA structure, the core and dimerization domains are in different positions to those seen in NhaA, and a negatively charged cavity has now opened to the outside. The extracellular cavity allows access to a strictly conserved aspartate residue thought to coordinate ion binding1, 8, 9 directly, a role supported here by molecular dynamics simulations. To alternate access to this ion-binding site, however, requires a surprisingly large rotation of the core domain, some 20° against the dimerization interface. We conclude that despite their fast transport rates of up to 1,500 ions per second3, Na+/H+ antiporters operate by a two-domain rocking bundle model, revealing themes relevant to secondary-active transporters in general

JGromacs: A Java Package for Analyzing Protein Simulations

UNLABELLED: In this paper, we introduce JGromacs, a Java API (Application Programming Interface) that facilitates the development of cross-platform data analysis applications for Molecular Dynamics (MD) simulations. The API supports parsing and writing file formats applied by GROMACS (GROningen MAchine for Chemical Simulations), one of the most widely used MD simulation packages. JGromacs builds on the strengths of object-oriented programming in Java by providing a multilevel object-oriented representation of simulation data to integrate and interconvert sequence, structure, and dynamics information. The easy-to-learn, easy-to-use, and easy-to-extend framework is intended to simplify and accelerate the implementation and development of complex data analysis algorithms. Furthermore, a basic analysis toolkit is included in the package. The programmer is also provided with simple tools (e.g., XML-based configuration) to create applications with a user interface resembling the command-line interface of GROMACS applications. AVAILABILITY: JGromacs and detailed documentation is freely available from http://sbcb.bioch.ox.ac.uk/jgromacs under a GPLv3 license

Allosteric effects in cyclophilin mutants may be explained by changes in nano-microsecond time scale motions

The relationship between molecular motion and catalysis in enzymes is debated. Here, simulations of cyclophilin A and three catalytically-impaired mutants reveal a nanosecond-scale interconversion between active and inactive conformations, orders of magnitude faster than previously suggested

Structural basis of outer membrane protein insertion by the BAM complex

All Gram-negative bacteria, mitochondria and chloroplasts have outer membrane proteins (OMPs) that perform many fundamental biological processes. The OMPs in Gram-negative bacteria are inserted and folded into the outer membrane by the β-barrel assembly machinery (BAM). The mechanism involved is poorly understood, owing to the absence of a structure of the entire BAM complex. Here we report two crystal structures of the Escherichia coli BAM complex in two distinct states: an inward-open state and a lateral-open state. Our structures reveal that the five polypeptide transport-associated domains of BamA form a ring architecture with four associated lipoproteins, BamB–BamE, in the periplasm. Our structural, functional studies and molecular dynamics simulations indicate that these subunits rotate with respect to the integral membrane β-barrel of BamA to induce movement of the β-strands of the barrel and promote insertion of the nascent OMP

Clinically relevant mutations in the ABCG2 transporter uncovered by genetic analysis linked to erythrocyte membrane protein expression

The ABCG2 membrane protein is a key xeno- and endobiotic transporter, modulating the absorption and metabolism of pharmacological agents and causing multidrug resistance in cancer. ABCG2 is also involved in uric acid elimination and its impaired function is causative in gout. Analysis of ABCG2 expression in the erythrocyte membranes of healthy volunteers and gout patients showed an enrichment of lower expression levels in the patients. By genetic screening based on protein expression, we found a relatively frequent, novel ABCG2 mutation (ABCG2-M71V), which, according to cellular expression studies, causes reduced protein expression, although with preserved transporter capability. Molecular dynamics simulations indicated a stumbled dynamics of the mutant protein, while ABCG2-M71V expression in vitro could be corrected by therapeutically relevant small molecules. These results suggest that personalized medicine should consider this newly discovered ABCG2 mutation, and genetic analysis linked to protein expression provides a new tool to uncover clinically important mutations in membrane proteins. © 2018 The Author(s)

CHAMPION: Chalmers Hierarchical Atomic, Molecular, Polymeric & Ionic Analysis Toolkit

We present CHAMPION: a software developed to automatically detect time-dependent bonds between atoms based on their dynamics, classify the local graph topology around them, and analyze the physicochemical properties of these topologies by statistical physics. In stark contrast to methodologies where bonds are detected based on static conditions such as cut-off distances, CHAMPION considers pairs of atoms to be bound only if they move together and act as a bound pair over time. Furthermore, the time-dependent global bond graph is possible to split into dynamically shifting connected components or subgraphs around a certain chemical motif and thereby allow the physicochemical properties of each such topology to be analyzed by statistical physics. Applicable to condensed matter and liquids in general, and electrolytes in particular, this allows both quantitative and qualitative descriptions of local structure, as well as dynamical processes such as speciation and diffusion. We present here a detailed overview of CHAMPION, including its underlying methodology, implementation and capabilities.Comment: 11 pages, 8 figure

Ceramides bind VDAC2 to trigger mitochondrial apoptosis

Ceramides draw wide attention as tumor suppressor lipids that act directly on mitochondria to trigger apoptotic cell death. However, molecular details of the underlying mechanism are largely unknown. Using a photoactivatable ceramide probe, we here identify the voltage-dependent anion channels VDAC1 and VDAC2 as mitochondrial ceramide binding proteins. Coarse-grain molecular dynamics simulations reveal that both channels harbor a ceramide binding site on one side of the barrel wall. This site includes a membrane-buried glutamate that mediates direct contact with the ceramide head group. Substitution or chemical modification of this residue abolishes photolabeling of both channels with the ceramide probe. Unlike VDAC1 removal, loss of VDAC2 or replacing its membrane-facing glutamate with glutamine renders human colon cancer cells largely resistant to ceramide-induced apoptosis. Collectively, our data support a role of VDAC2 as direct effector of ceramide-mediated cell death, providing a molecular framework for how ceramides exert their anti-neoplastic activity

The role of interfacial lipids in stabilizing membrane protein oligomers

Oligomerization of membrane proteins in response to lipid binding has a critical role in many cell-signalling pathways1 but is often difficult to define2 or predict3. Here we report the development of a mass spectrometry platform to determine simultaneously the presence of interfacial lipids and oligomeric stability and to uncover how lipids act as key regulators of membrane-protein association. Evaluation of oligomeric strength for a dataset of 125 α-helical oligomeric membrane proteins reveals an absence of interfacial lipids in the mass spectra of 12 membrane proteins with high oligomeric stability. For the bacterial homologue of the eukaryotic biogenic transporters (LeuT4, one of the proteins with the lowest oligomeric stability), we found a precise cohort of lipids within the dimer interface. Delipidation, mutation of lipid-binding sites or expression in cardiolipin-deficient Escherichia coli abrogated dimer formation. Molecular dynamics simulation revealed that cardiolipin acts as a bidentate ligand, bridging across subunits. Subsequently, we show that for the Vibrio splendidus sugar transporter SemiSWEET5, another protein with low oligomeric stability, cardiolipin shifts the equilibrium from monomer to functional dimer. We hypothesized that lipids are essential for dimerization of the Na+/H+ antiporter NhaA from E. coli, which has the lowest oligomeric strength, but not for the substantially more stable homologous Thermus thermophilus protein NapA. We found that lipid binding is obligatory for dimerization of NhaA, whereas NapA has adapted to form an interface that is stable without lipids. Overall, by correlating interfacial strength with the presence of interfacial lipids, we provide a rationale for understanding the role of lipids in both transient and stable interactions within a range of α-helical membrane proteins, including G-protein-coupled receptors

The Integrin Receptor in Biologically Relevant Bilayers: Insights from Molecular Dynamics Simulations

Integrins are heterodimeric (αβ) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2–F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2–F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes