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

Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system

Docking and molecular dynamics simulations of the ternary complex nisin2:lipid II

Lanthionine antibiotics are an important class of naturally-occurring antimicrobial peptides. The best-known, nisin, is a commercial food preservative. However, structural and mechanistic details on nisin/lipid II membrane complexes are currently lacking. Recently, we have developed empirical force-field parameters to model lantibiotics. Docking and molecular dynamics (MD) simulations have been used to study the nisin2:lipid II complex in bacterial membranes, which has been put forward as the building block of nisin/lipid II binary membrane pores. A Ile1Trp mutation of the N-terminus of nisin has been modelled and docked onto lipid II models; the computed binding affinity increased compared to wildtype. Wild-type nisin was also docked onto three different lipid II structures and a stable 2:1 nisin:lipid II complex formed. This complex was inserted into a membrane. Six independent MD simulations revealed key interactions in the complex, specifically the N terminal engagement of nisin with lipid II at the pyrophosphate and C-terminus of the pentapeptide chain. Nisin2 inserts into the membrane and we propose this is the first step in pore formation, mediated by the nisin N-terminus–lipid II pentapeptide hydrogen bond. The lipid II undecaprenyl chain adopted different conformations in the presence of nisin, which may also have implications for pore formation

Evaluating the Effects of Cutoffs and Treatment of Long-range Electrostatics in Protein Folding Simulations

The use of molecular dynamics simulations to provide atomic-level descriptions of biological processes tends to be computationally demanding, and a number of approximations are thus commonly employed to improve computational efficiency. In the past, the effect of these approximations on macromolecular structure and stability has been evaluated mostly through quantitative studies of small-molecule systems or qualitative observations of short-timescale simulations of biological macromolecules. Here we present a quantitative evaluation of two commonly employed approximations, using a test system that has been the subject of a number of previous protein folding studies–the villin headpiece. In particular, we examined the effect of (i) the use of a cutoff-based force-shifting technique rather than an Ewald summation for the treatment of electrostatic interactions, and (ii) the length of the cutoff used to determine how many pairwise interactions are included in the calculation of both electrostatic and van der Waals forces. Our results show that the free energy of folding is relatively insensitive to the choice of cutoff beyond 9 Å, and to whether an Ewald method is used to account for long-range electrostatic interactions. In contrast, we find that the structural properties of the unfolded state depend more strongly on the two approximations examined here

Outer membrane protein size and LPS O-antigen define protective antibody targeting to the Salmonella surface

Lipopolysaccharide (LPS) O-antigen (O-Ag) is known to limit antibody binding to surface antigens, although the relationship between antibody, O-Ag and other outer-membrane antigens is poorly understood. Here we report, immunization with the trimeric porin OmpD from Salmonella Typhimurium (STmOmpD) protects against infection. Atomistic molecular dynamics simulations indicate this is because OmpD trimers generate footprints within the O-Ag layer sufficiently sized for a single IgG Fab to access. While STmOmpD differs from its orthologue in S. Enteritidis (SEn) by a single amino-acid residue, immunization with STmOmpD confers minimal protection to SEn. This is due to the OmpD-O-Ag interplay restricting IgG binding, with the pairing of OmpD with its native O-Ag being essential for optimal protection after immunization. Thus, both the chemical and physical structure of O-Ag are key for the presentation of specific epitopes within proteinaceous surface-antigens. This enhances combinatorial antigenic diversity in Gram-negative bacteria, while reducing associated fitness costs

Parameterization of a coarse-grained model of cholesterol with point-dipole electrostatics

© 2018, Springer Nature Switzerland AG. We present a new coarse-grained (CG) model of cholesterol (CHOL) for the electrostatic-based ELBA force field. A distinguishing feature of our CHOL model is that the electrostatics is modeled by an explicit point dipole which interacts through an ideal vacuum permittivity. The CHOL model parameters were optimized in a systematic fashion, reproducing the electrostatic and nonpolar partitioning free energies of CHOL in lipid/water mixtures predicted by full-detailed atomistic molecular dynamics simulations. The CHOL model has been validated by comparison to structural, dynamic and thermodynamic properties with experimental and atomistic simulation reference data. The simulation of binary DPPC/cholesterol mixtures covering the relevant biological content of CHOL in mammalian membranes is shown to correctly predict the main lipid behavior as observed experimentally

Membrainy: a ‘smart’, unified membrane analysis tool

BACKGROUND: The study of biological membranes using Molecular Dynamics has become an increasingly popular means by which to investigate the interactions of proteins, peptides and potentials with lipid bilayers. These interactions often result in changes to the properties of the lipids which can modify the behaviour of the membrane. Membrainy is a unified membrane analysis tool that contains a broad spectrum of analytical techniques to enable: measurement of acyl chain order parameters; presentation of 2D surface and thickness maps; determination of lateral and axial headgroup orientations; measurement of bilayer and leaflet thickness; analysis of the annular shell surrounding membrane-embedded objects; quantification of gel percentage; time evolution of the transmembrane voltage; area per lipid calculations; and quantification of lipid mixing/demixing entropy. RESULTS: Each analytical component within Membrainy has been tested on a variety of lipid bilayer systems and was found to be either comparable to or an improvement upon existing software. For the analytical techniques that have no direct comparable software, our results were confirmed with experimental data. CONCLUSIONS: Membrainy is a user-friendly, intelligent membrane analysis tool that automatically interprets a variety of input formats and force fields, is compatible with both single and double bilayers, and capable of handling asymmetric bilayers and lipid flip-flopping. Membrainy has been designed for ease of use, requiring no installation or configuration and minimal user-input to operate

Structure and mechanism of Zn^(2+)- transporting P-type ATPases

Zinc is an essential micronutrient for all living organisms. It is required for signalling and proper functioning of a range of proteins involved in, for example, DNA binding and enzymatic catalysis. In prokaryotes and photosynthetic eukaryotes, Zn2+-transporting P-type ATPases of class IB (ZntA) are crucial for cellular redistribution and detoxification of Zn2+ and related elements. Here we present crystal structures representing the phosphoenzyme ground state (E2P) and a dephosphorylation intermediate (E2·P_i) of ZntA from Shigella sonnei, determined at 3.2 Å and 2.7 Å resolution, respectively. The structures reveal a similar fold to Cu^+-ATPases, with an amphipathic helix at the membrane interface. A conserved electronegative funnel connects this region to the intramembranous high-affinity ion-binding site and may promote specific uptake of cellular Zn^(2+) ions by the transporter. The E2P structure displays a wide extracellular release pathway reaching the invariant residues at the high-affinity site, including C392, C394 and D714. The pathway closes in the E2·P_i state, in which D714 interacts with the conserved residue K693, which possibly stimulates Zn^(2+) release as a built-in counter ion, as has been proposed for H^+-ATPases. Indeed, transport studies in liposomes provide experimental support for ZntA activity without counter transport. These findings suggest a mechanistic link between P_(IB)-type Zn^(2+)-ATPases and P_(III)-type H^+-ATPases and at the same time show structural features of the extracellular release pathway that resemble P_(II)-type ATPases such as the sarcoplasmic/endoplasmic reticulum Ca^(2+)-ATPase (SERCA) and Na^+, K^+-ATPase. These findings considerably increase our understanding of zinc transport in cells and represent new possibilities for biotechnology and biomedicine

Structural and functional insights into oligopeptide acquisition by the RagAB transporter from Porphyromonas gingivalis

Porphyromonas gingivalis, an asaccharolytic member of the Bacteroidetes, is a keystone pathogen in human periodontitis that may also contribute to the development of other chronic inflammatory diseases. P. gingivalis utilizes protease-generated peptides derived from extracellular proteins for growth, but how these peptides enter the cell is not clear. Here, we identify RagAB as the outer-membrane importer for these peptides. X-ray crystal structures show that the transporter forms a dimeric RagA2B2 complex, with the RagB substrate-binding surface-anchored lipoprotein forming a closed lid on the RagA TonB-dependent transporter. Cryo-electron microscopy structures reveal the opening of the RagB lid and thus provide direct evidence for a ‘pedal bin’ mechanism of nutrient uptake. Together with mutagenesis, peptide-binding studies and RagAB peptidomics, our work identifies RagAB as a dynamic, selective outer-membrane oligopeptide-acquisition machine that is essential for the efficient utilization of proteinaceous nutrients by P. gingivalis

Allosteric activation of an ion channel triggered by modification of mechanosensitive nano-pockets

Lipid availability within transmembrane nano-pockets of ion channels is linked with mechanosensation. However, the effect of hindering lipid-chain penetration into nano-pockets on channel structure has not been demonstrated. Here we identify nano-pockets on the large conductance mechanosensitive channel MscL, the high-pressure threshold channel. We restrict lipid-chain access to the nano-pockets by mutagenesis and sulfhydryl modification, and monitor channel conformation by PELDOR/DEER spectroscopy. For a single site located at the entrance of the nano-pockets and distal to the channel pore we generate an allosteric response in the absence of tension. Single-channel recordings reveal a significant decrease in the pressure activation threshold of the modified channel and a sub-conducting state in the absence of applied tension. Threshold is restored to wild-type levels upon reduction of the sulfhydryl modification. The modification associated with the conformational change restricts lipid access to the nano-pocket, interrupting the contact between the membrane and the channel that mediates mechanosensitivity