Proton-coupled dynamics in lactose permease

Ion and lipid orchestration of secondary active transport

Transporting small molecules across cell membranes is an essential process in cell physiology. Many structurally diverse, secondary active transporters harness transmembrane electrochemical gradients of ions to power the uptake or efflux of nutrients, signalling molecules, drugs and other ions across cell membranes. Transporters reside in lipid bilayers on the interface between two aqueous compartments, where they are energized and regulated by symported, antiported and allosteric ions on both sides of the membrane and the membrane bilayer itself. Here we outline the mechanisms by which transporters couple ion and solute fluxes and discuss how structural and mechanistic variations enable them to meet specific physiological needs and adapt to environmental conditions. We then consider how general bilayer properties and specific lipid binding modulate transporter activity. Together, ion gradients and lipid properties ensure the effective transport, regulation and distribution of small molecules across cell membranes.

Aminoglycoside uptake, stress, and potentiation in Gram-negative bacteria: new therapies with old molecules

SUMMARYAminoglycosides (AGs) are long-known molecules successfully used against Gram-negative pathogens. While their use declined with the discovery of new antibiotics, they are now classified as critically important molecules because of their effectiveness against multidrug-resistant bacteria. While they can efficiently cross the Gram-negative envelope, the mechanism of AG entry is still incompletely understood, although this comprehension is essential for the development of new therapies in the face of the alarming increase in antibiotic resistance. Increasing antibiotic uptake in bacteria is one strategy to enhance effective treatments. This review aims, first, to consolidate old and recent knowledge about AG uptake; second, to explore the connection between AG-dependent bacterial stress and drug uptake; and finally, to present new strategies of potentiation of AG uptake for more efficient antibiotic therapies. In particular, we emphasize on the connection between sugar transport and AG potentiation.

Dynamic lipid interactions in the plasma membrane Na+,K+-ATPase

The function of ion-transporting Na+,K+-ATPases depends on the surrounding lipid environment in biological membranes. Two established lipid-interaction sites A and B within the transmembrane domain have been observed to induce protein activation and stabilization, respectively. In addition, lipid-mediated inhibition has been assigned to a site C, but with the exact location not experimentally confirmed. Also, possible effects on lipid interactions by disease mutants dwelling in the membrane-protein interface remain relatively uncharacterized. We simulated human Na+,K+-ATPase α1β1FXYD homology models in E1 and E2 states in an asymmetric, multicomponent plasma membrane to determine both wild-type and disease mutant lipid-protein interactions. The simulated wild-type lipid interactions at the established sites A and B were in agreement with experimental results thereby confirming the membrane-protein model system. The less well-characterized, proposed inhibitory site C was dominated by lipids lacking inhibitory properties. Instead, two sites hosting inhibitory lipids were identified at the extracellular side and also a cytoplasmic CHL-binding site that provide putative alternative locations of Na+,K+-ATPase inhibition. Three disease mutations, Leu302Arg, Glu840Arg and Met859Arg resided in the lipid-protein interface and caused drastic changes in the lipid interactions. The simulation results show that lipid interactions to the human Na+,K+-ATPase α1β1FXYD protein in the plasma membrane are highly state-dependent and can be disturbed by disease mutations located in the lipid interface, which can open up for new venues to understand genetic disorders.

Structural basis of purine nucleotide inhibition of human uncoupling protein 1

Mitochondrial uncoupling protein 1 (UCP1) gives brown adipose tissue of mammals its specialized ability to burn calories as heat for thermoregulation. When activated by fatty acids, UCP1 catalyzes the leak of protons across the mitochondrial inner membrane, short-circuiting the mitochondrion to generate heat, bypassing ATP synthesis. In contrast, purine nucleotides bind and inhibit UCP1, regulating proton leak by a molecular mechanism that is unclear. We present the cryo-electron microscopy structure of the GTP-inhibited state of UCP1, which is consistent with its nonconducting state. The purine nucleotide cross-links the transmembrane helices of UCP1 with an extensive interaction network. Our results provide a structural basis for understanding the specificity and pH dependency of the regulatory mechanism. UCP1 has retained all of the key functional and structural features required for a mitochondrial carrier-like transport mechanism. The analysis shows that inhibitor binding prevents the conformational changes that UCP1 uses to facilitate proton leak.

Insights into the structure and function of the human organic anion transporter 1 in lipid bilayer membranes

The human SLC22A6/OAT1 plays an important role in the elimination of a broad range of endogenous substances and xenobiotics thus attracting attention from the pharmacological community. Furthermore, OAT1 is also involved in key physiological events such as the remote inter-organ communication. Despite its significance, the knowledge about hOAT1 structure and the transport mechanism at the atomic level remains fragmented owing to the lack of resolved structures. By means of protein-threading modeling refined by μs-scaled Molecular Dynamics simulations, the present study provides the first robust model of hOAT1 in outward-facing conformation. Taking advantage of the AlphaFold 2 predicted structure of hOAT1 in inward-facing conformation, we here provide the essential structural and functional features comparing both states. The intracellular motifs conserved among Major Facilitator Superfamily members create a so-called “charge-relay system” that works as molecular switches modulating the conformation. The principal element of the event points at interactions of charged residues that appear crucial for the transporter dynamics and function. Moreover, hOAT1 model was embedded in different lipid bilayer membranes highlighting the crucial structural dependence on lipid-protein interactions. MD simulations supported the pivotal role of phosphatidylethanolamine components to the protein conformation stability. The present model is made available to decipher the impact of any observed polymorphism and mutation on drug transport as well as to understand substrate binding modes.

Mechanisms of long-distance allosteric couplings in proton-binding membrane transporters

Membrane transporters that use proton binding and proton transfer for function couple local protonation change with changes in protein conformation and water dynamics. Changes of protein conformation might be required to allow transient formation of hydrogen-bond networks that bridge proton donor and acceptor pairs separated by long distances. Inter-helical hydrogen-bond networks adjust rapidly to protonation change, and ensure rapid response of the protein structure and dynamics. Membrane transporters with known three-dimensional structures and proton-binding groups inform on general principles of protonation-coupled protein conformational dynamics. Inter-helical hydrogen bond motifs between proton-binding carboxylate groups and a polar sidechain are observed in unrelated membrane transporters, suggesting common principles of coupling protonation change with protein conformational dynamics.

Tracking the ATP-binding response in adenylate kinase in real time

The biological function of proteins is critically dependent on dynamics inherent to the native structure. Such structural dynamics obey a predefined order and temporal timing to execute the specific reaction. Determination of the cooperativity of key structural rearrangements requires monitoring protein reactions in real time. In this work, we used time-resolved x-ray solution scattering (TR-XSS) to visualize structural changes in the Escherichia coli adenylate kinase (AdK) enzyme upon laser-induced activation of a protected ATP substrate. A 4.3-ms transient intermediate showed partial closing of both the ATP- and AMP-binding domains, which indicates a cooperative closing mechanism. The ATP-binding domain also showed local unfolding and breaking of an Arg¹³¹-Asp¹⁴⁶ salt bridge. Nuclear magnetic resonance spectroscopy data identified similar unfolding in an Arg¹³¹Ala AdK mutant, which refolded in a closed, substrate-binding conformation. The observed structural dynamics agree with a “cracking mechanism” proposed to underlie global structural transformation, such as allostery, in proteins.

Proton-Binding Motifs of Membrane-Bound Proteins: From Bacteriorhodopsin to Spike Protein S

Membrane-bound proteins that change protonation during function use specific protein groups to bind and transfer protons. Knowledge of the identity of the proton-binding groups is of paramount importance to decipher the reaction mechanism of the protein, and protonation states of prominent are studied extensively using experimental and computational approaches. Analyses of model transporters and receptors from different organisms, and with widely different biological functions, indicate common structure-sequence motifs at internal proton-binding sites. Proton-binding dynamic hydrogen-bond networks that are exposed to the bulk might provide alternative proton-binding sites and proton-binding pathways. In this perspective article I discuss protonation coupling and proton binding at internal and external carboxylate sites of proteins that use proton transfer for function. An inter-helical carboxylate-hydroxyl hydrogen-bond motif is present at functionally important sites of membrane proteins from archaea to the brain. External carboxylate-containing H-bond clusters are observed at putative proton-binding sites of protonation-coupled model proteins, raising the question of similar functionality in spike protein S.

The Whole Is Bigger than the Sum of Its Parts: Drug Transport in the Context of Two Membranes with Active Efflux

Cell envelope plays a dual role in the life of bacteria by simultaneously protecting it from a hostile environment and facilitating access to beneficial molecules. At the heart of this ability lie the restrictive properties of the cellular membrane augmented by efflux transporters, which preclude intracellular penetration of most molecules except with the help of specialized uptake mediators. Recently, kinetic properties of the cell envelope came into focus driven on one hand by the urgent need in new antibiotics and, on the other hand, by experimental and theoretical advances in studies of transmembrane transport. A notable result from these studies is the development of a kinetic formalism that integrates the Michaelis-Menten behavior of individual transporters with transmembrane diffusion and offers a quantitative basis for the analysis of intracellular penetration of bioactive compounds. This review surveys key experimental and computational approaches to the investigation of transport by individual translocators and in whole cells, summarizes key findings from these studies and outlines implications for antibiotic discovery. Special emphasis is placed on Gram-negative bacteria, whose envelope contains two separate membranes. This feature sets these organisms apart from Gram-positive bacteria and eukaryotic cells by providing them with full benefits of the synergy between slow transmembrane diffusion and active efflux.

Structures and General Transport Mechanisms by the Major Facilitator Superfamily (MFS)

The major facilitator superfamily (MFS) is the largest known superfamily of secondary active transporters. MFS transporters are responsible for transporting a broad spectrum of substrates, either down their concentration gradient or uphill using the energy stored in the electrochemical gradients. Over the last 10 years, more than a hundred different MFS transporter structures covering close to 40 members have provided an atomic framework for piecing together the molecular basis of their transport cycles. Here, we summarize the remarkable promiscuity of MFS members in terms of substrate recognition and proton coupling as well as the intricate gating mechanisms undergone in achieving substrate translocation. We outline studies that show how residues far from the substrate binding site can be just as important for fine-tuning substrate recognition and specificity as those residues directly coordinating the substrate, and how a number of MFS transporters have evolved to form unique complexes with chaperone and signaling functions. Through a deeper mechanistic description of glucose (GLUT) transporters and multidrug resistance (MDR) antiporters, we outline novel refinements to the rocker-switch alternating-access model, such as a latch mechanism for proton-coupled monosaccharide transport. We emphasize that a full understanding of transport requires an elucidation of MFS transporter dynamics, energy landscapes, and the determination of how rate transitions are modulated by lipids.

Tracking Membrane Protein Dynamics in Real Time

Abstract

Membrane proteins govern critical cellular processes and are central to human health and associated disease. Understanding of membrane protein function is obscured by the vast ranges of structural dynamics—both in the spatial and time regime—displayed in the protein and surrounding membrane. The membrane lipids have emerged as allosteric modulators of membrane protein function, which further adds to the complexity. In this review, we discuss several examples of membrane dependency. A particular focus is on how molecular dynamics (MD) simulation have aided to map membrane protein dynamics and how enhanced sampling methods can enable observing the otherwise inaccessible biological time scale. Also, time-resolved X-ray scattering in solution is highlighted as a powerful tool to track membrane protein dynamics, in particular when combined with MD simulation to identify transient intermediate states. Finally, we discuss future directions of how to further develop this promising approach to determine structural dynamics of both the protein and the surrounding lipids.

Graphic Abstract

Structural and Functional Adaptability of Sucrose and Lactose Permeases from Escherichia coli to the Membrane Lipid Composition

The lipid environment in which membrane proteins are embedded can influence their structure and function. Lipid-protein interactions and lipid-induced conformational changes necessary for protein function remain intractable in vivo using high-resolution techniques. Using Escherichia coli strains in which the normal phospholipid composition can be altered or foreign lipids can be introduced, we established the importance of membrane lipid composition for the proper folding, assembly, and function of E. coli lactose (LacY) and sucrose (CscB) permeases. However, the molecular mechanism underlying the lipid dependence for active transport remains unknown. Herein, we demonstrate that the structure and function of CscB and LacY can be modulated by the composition of the lipid environment. Using a combination of assays (transport activity of the substrate, protein topology, folding, and assembly into the membrane), we found that alterations in the membrane lipid composition lead to lipid-dependent structural changes in CscB and LacY. These changes affect the orientation of residues involved in LacY proton translocation and impact the rates of protonation and deprotonation of E325 by affecting the arrangement of transmembrane domains in the vicinity of the R302-E325 charge pair. Furthermore, the structural changes caused by changes in membrane lipid composition can be altered by a single-point mutation, highlighting the adaptability of these transporters to their environment. Altogether, our results demonstrate that direct interactions between a protein and its lipid environment uniquely contribute to membrane protein organization and function. Because members of the major facilitator superfamily present with well-conserved functional architecture, we anticipate that our findings can be extrapolated to other membrane protein transporters.

The proton electrochemical gradient induces a kinetic asymmetry in the symport cycle of LacY

Protonation and deprotonation of Glu325 with a pK_a of 10.5 is required for symport. Moreover, the H⁺ electrochemical gradient (∆μ∼H+) accelerates deprotonation on the intracellular side with a 50- to 100-fold decrease in the K_m. To probe the pK on the cytoplasmic side of the membrane, rates of lactose/H⁺ efflux were determined from pH 5.0 to 9.0 without or with a membrane potential (ΔΨ, interior positive) in right-side-out membrane vesicles. WT lactose efflux has an apparent pK of ∼7.2 that is unaffected by ΔΨ, mutant E325A is defective, and pH or ΔΨ (interior positive) has no effect. The effect of ΔΨ (interior positive) on the K_m for efflux with WT LacY is insignificant relative to the marked effect on influx.

LacY catalyzes accumulation of galactosides against a concentration gradient by coupling galactoside and H⁺ transport (i.e., symport). While alternating access of sugar- and H⁺-binding sites to either side of the membrane is driven by binding and dissociation of sugar, the electrochemical H⁺ gradient (∆μ∼H+) functions kinetically by decreasing the K_m for influx 50- to 100-fold with no change in K_d. The affinity of protonated LacY for sugar has an apparent pK (pK^app) of ∼10.5, due specifically to the pK_a of Glu325, a residue that plays an irreplaceable role in coupling. In this study, rates of lactose/H⁺ efflux were measured from pH 5.0 to 9.0 in the absence or presence of a membrane potential (ΔΨ, interior positive), and the effect of the imposed ΔΨ on the kinetics of efflux was also studied in right-side-out membrane vesicles. The findings reveal that ∆μ∼H+ induces an asymmetry in the transport cycle based on the following observations: 1) the efflux rate of WT LacY exhibits a pK^app of ∼7.2 that is unaffected by the imposed ΔΨ; 2) ΔΨ increases the rate of efflux at all tested pH values, but enhancement is almost 2 orders of magnitude less than observed for influx; 3) mutant Glu325 ˗ Ala does little or no efflux in the absence or presence of ΔΨ, and ambient pH has no effect; and 4) the effect of ΔΨ (interior positive) on the K_m for efflux is almost insignificant relative to the 50- to 100-fold decrease in the K_m for influx driven by ΔΨ (interior negative).

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Secondary Active Transporters

Transport of solutes across biological membranes is essential for cellular life. This process is mediated by membrane transport proteins which move nutrients, waste products, certain drugs and ions into and out of cells. Secondary active transporters couple the transport of substrates against their concentration gradients with the transport of other solutes down their concentration gradients. The alternating access model of membrane transporters and the coupling mechanism of secondary active transporters are introduced in this book chapter. Structural studies have identified typical protein folds for transporters that we exemplify by the major facilitator superfamily (MFS) and LeuT folds. Finally, substrate binding and substrate translocation of the transporters LacY of the MFS and AdiC of the amino acid-polyamine-organocation (APC) superfamily are described.

Phylogenetic and conservation analyses of MFS transporters

Major facilitator superfamily is one of the largest superfamily of secondary transporters present across the kingdom of life. Considering the physiological and clinical importance of MFS proteins, we attempted to explore the phylogenetic and structural aspects of the superfamily. To achieve the objectives, we performed global sequence-based analyses of MFS proteins encompassing multiple taxa. Notably, phylogenetic analysis of MFS proteins resulted in the clustering of MFS proteins based on their function, rather than lineage of the respective organisms. Additionally, we employed information theoretic measures, Relative entropy (RE) and Cumulative relative entropy (CRE) to decipher fold-specific and function-specific residues, respectively, in the MFS proteins. The residues with high RE score when mapped on to the 3D-structure of MFS transporter LacY, were found to be distributed throughout the tertiary structure of the protein. On the other hand, CRE calculation was employed to contrast two subfamilies Drug H⁺ antiporter 1 and 3 (DHA1 and DHA3). The particular analysis unveiled certain differentially conserved residues in DHA1 as compared to DHA3 highlighting family-specific importance of them. Remarkably, a number of high scoring CRE residues have already established functional roles, for instance, the arginine residue present in TMH4. Altogether, the current study apart from providing an insight into the functional clustering of MFS proteins also identifies residues with established or plausible roles in the transport mechanism. Thus, the study lays a platform for future structure-function studies of these proteins.

Direct protein-lipid interactions shape the conformational landscape of secondary transporters

Secondary transporters undergo structural rearrangements to catalyze substrate translocation across the cell membrane – yet how such conformational changes happen within a lipid environment remains poorly understood. Here, we combine hydrogen-deuterium exchange mass spectrometry (HDX-MS) with molecular dynamics (MD) simulations to understand how lipids regulate the conformational dynamics of secondary transporters at the molecular level. Using the homologous transporters XylE, LacY and GlpT from Escherichia coli as model systems, we discover that conserved networks of charged residues act as molecular switches that drive the conformational transition between different states. We reveal that these molecular switches are regulated by interactions with surrounding phospholipids and show that phosphatidylethanolamine interferes with the formation of the conserved networks and favors an inward-facing state. Overall, this work provides insights into the importance of lipids in shaping the conformational landscape of an important class of transporters.

Secondary transporters catalyse substrate translocation across the cell membrane but the role of lipids during the transport cycle remains unclear. Here authors used hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations to understand how lipids regulate the conformational dynamics of secondary transporters.

Uptake dynamics in the Lactose permease (LacY) membrane protein transporter

The sugar transporter Lactose permease (LacY) of Escherichia coli has become a prototype to understand the underlying molecular details of membrane transport. Crystal structures have trapped the protein in sugar-bound states facing the periplasm, but with narrow openings unable to accommodate sugar. Therefore, the molecular details of sugar uptake remain elusive. In this work, we have used extended simulations and metadynamics sampling to explore a putative sugar-uptake pathway and associated free energy landscape. We found an entrance at helix-pair 2 and 11, which involved lipid head groups and residues Gln 241 and Gln 359. Furthermore, the protein displayed high flexibility on the periplasmic side of Phe 27, which is located at the narrowest section of the pathway. Interactions to Phe 27 enabled passage into the binding site, which was associated with a 24 ± 4 kJ/mol binding free energy in excellent agreement with an independent binding free energy calculation and experimental data. Two free energy minima corresponding to the two possible binding poses of the lactose analog β-D-galactopyranosyl-1-thio-β-D-galactopyranoside (TDG) were aligned with the crystal structure-binding pocket. This work outlines the chemical environment of a putative periplasmic sugar pathway and paves way for understanding substrate affinity and specificity in LacY.

Dynamics of Long-Distance Hydrogen-Bond Networks in Photosystem II

Photosystem II uses the energy of absorbed light to split water molecules, generating molecular oxygen, electrons, and protons. The four protons generated during each reaction cycle are released to the lumen via mechanisms that are poorly understood. Given the complexity of photosystem II, which consists of multiple protein subunits and cofactor molecules and hosts numerous waters, a fundamental issue is finding transient networks of hydrogen bonds that bridge potential proton donor and acceptor groups. Here, we address this issue by performing all-atom molecular dynamics simulations of wild-type and mutant photosystem II monomers, which we analyze using a new protocol designed to facilitate efficient analysis of hydrogen-bond networks. Our computations reveal that local protein/water hydrogen-bond networks can assemble transiently in photosystem II such that the reaction center connects to the lumen. The dynamics of the hydrogen-bond networks couple to the protonation state of specific carboxylate groups and are altered in a mutant with defective proton transfer. Simulations on photosystem II without its extrinsic PsbO subunit provide a molecular interpretation of the elusive functional role of this subunit.

Proton movement and coupling in the POT family of peptide transporters

POT transporters represent an evolutionarily well-conserved family of proton-coupled transport systems in biology. An unusual feature of the family is their ability to couple the transport of chemically diverse ligands to an inwardly directed proton electrochemical gradient. For example, in mammals, fungi, and bacteria they are predominantly peptide transporters, whereas in plants the family has diverged to recognize nitrate, plant defense compounds, and hormones. Although recent structural and biochemical studies have identified conserved sites of proton binding, the mechanism through which transport is coupled to proton movement remains enigmatic. Here we show that different POT transporters operate through distinct proton-coupled mechanisms through changes in the extracellular gate. A high-resolution crystal structure reveals the presence of ordered water molecules within the peptide binding site. Multiscale molecular dynamics simulations confirm proton transport occurs through these waters via Grotthuss shuttling and reveal that proton binding to the extracellular side of the transporter facilitates a reorientation from an inward- to outward-facing state. Together these results demonstrate that within the POT family multiple mechanisms of proton coupling have likely evolved in conjunction with variation of the extracellular gate.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Proton-solute coupling mechanism of the maltose transporter from Saccharomyces cerevisiae

Mal11 catalyzes proton-coupled maltose transport across the plasma membrane of Saccharomyces cerevisiae. We used structure-based design of mutants and a kinetic analysis of maltose transport to determine the energy coupling mechanism of transport. We find that wildtype Mal11 is extremely well coupled and allows yeast to rapidly accumulate maltose to dangerous levels, resulting under some conditions in self-lysis. Three protonatable residues lining the central membrane-embedded cavity of Mal11 were identified as having potential roles in proton translocation. We probed the mechanistic basis for proton coupling with uphill and downhill transport assays and found that single mutants can still accumulate maltose but with a lower coupling efficiency than the wildtype. Next, we combined the individual mutations and created double and triple mutants. We found some redundancy in the functions of the acidic residues in proton coupling and that no single residue is most critical for proton coupling to maltose uptake, unlike what is usually observed in related transporters. Importantly, the triple mutants were completely uncoupled but still fully active in downhill efflux and equilibrium exchange. Together, these results depict a concerted mechanism of proton transport in Mal11 involving multiple charged residues.

Revised CHARMM force field parameters for iron-containing cofactors of photosystem II

Photosystem II is a complex protein-cofactor machinery that splits water molecules into molecular oxygen, protons, and electrons. All-atom molecular dynamics simulations have the potential to contribute to our general understanding of how photosystem II works. To perform reliable all-atom simulations, we need accurate force field parameters for the cofactor molecules. We present here CHARMM bonded and non-bonded parameters for the iron-containing cofactors of photosystem II that include a six-coordinated heme moiety coordinated by two histidine groups, and a non-heme iron complex coordinated by bicarbonate and four histidines. The force field parameters presented here give water interaction energies and geometries in good agreement with the quantum mechanical target data. © 2017 Wiley Periodicals, Inc.

Critical Temperature of 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine Monolayers and Its Possible Biological Relevance

Because transmembrane proteins (TMPs) can be obtained with sufficient purity for X-ray diffraction studies more frequently than decades ago, their mechanisms of action may now be elucidated. One of the pending issues is the actual interplay between transmembrane proteins and membrane lipids. There is strong evidence of the involvement of specific lipids with some membrane proteins, such as the potassium crystallographically sited activation channel (KcsA) of Streptomyces lividans and the secondary transporter of lactose LacY of Escherichia coli, the activities of which are associated with the presence of anionic phospholipids such as the phosphatidylglycerol (PG) and phosphatidyethanolamine (PE), respectively. Other proteins such as the large conductance mechanosensitive channel (MscL) of E. coli seem to depend on the adaptation of specific phospholipids to the irregular surface of the integral membrane protein. In this work we investigated the lateral compressibility of two homoacid phosphatidylethanolamines (one with both acyl chains unsaturated (DOPE), the other with the acyl chains saturated (DPPE)) and the heteroacid phosphatidyletanolamine (POPE) and their mixtures with POPG. The liquid expanded (LE) to liquid condensed (LC) transition was observed in POPE at a temperature below its critical temperature (T_c = 36 °C). Because T_c lies below the physiological temperature, the occurrence of this phase transition may have something to do with the functioning of LacY. This magnitude is discussed within the context of the experiments carried out at temperatures below the T_c of POPE at which the activity of Lac Y and other TMPs are frequently studied.

pKa of Glu325 in LacY

Lactose permease (LacY), a paradigm for the largest family of membrane transport proteins, catalyzes the coupled translocation of a galactoside and a H⁺ across the cytoplasmic membrane of Escherichia coli (galactoside/H⁺ symport). One of the most important aspects of the mechanism is the relationship between protonation and binding of the cargo galactopyranoside. In this regard, it has been shown that protonation is required for binding. Furthermore when galactoside affinity is measured as a function of pH, an apparent pK (pK^app) of ∼10.5 is obtained. Strikingly, when Glu325, a residue long known to be involved in coupling between H⁺ and sugar translocation, is replaced with a neutral side chain, the pH effect is abolished, and high-affinity binding is observed until LacY is destabilized at alkaline pH. In this paper, infrared spectroscopy is used to identify Glu325 in situ. Moreover, it is demonstrated that this residue exhibits a pK_a of 10.5 ± 0.1 that is insensitive to the presence of galactopyranoside. Thus, it is apparent that protonation of Glu325 specifically is required for effective sugar binding to LacY.

Membrane Anchoring and Ion-Entry Dynamics in P-type ATPase Copper Transport

Cu+-specific P-type ATPase membrane protein transporters regulate cellular copper levels. The lack of crystal structures in Cu+-binding states has limited our understanding of how ion entry and binding are achieved. Here, we characterize the molecular basis of Cu+ entry using molecular-dynamics simulations, structural modeling, and in vitro and in vivo functional assays. Protein structural rearrangements resulting in the exposure of positive charges to bulk solvent rather than to lipid phosphates indicate a direct molecular role of the putative docking platform in Cu+ delivery. Mutational analyses and simulations in the presence and absence of Cu+ predict that the ion-entry path involves two ion-binding sites: one transient Met148-Cys382 site and one intramembranous site formed by trigonal coordination to Cys384, Asn689, and Met717. The results reconcile earlier biochemical and x-ray absorption data and provide a molecular understanding of ion entry in Cu+-transporting P-type ATPases.

Key Residues and Phosphate Release Routes in the Saccharomyces cerevisiae Pho84 Transceptor: THE ROLE OF TYR179 IN FUNCTIONAL REGULATION

Pho84, a major facilitator superfamily (MFS) protein, is the main high-affinity P_i transceptor in Saccharomyces cerevisiae Although transport mechanisms have been suggested for other MFS members, the key residues and molecular events driving transport by P_i:H⁺ symporters are unclear. The current Pho84 transport model is based on the inward-facing occluded crystal structure of the Pho84 homologue PiPT in the fungus Piriformospora indica However, this model is limited by the lack of experimental data on the regulatory residues for each stage of the transport cycle. In this study, an open, inward-facing conformation of Pho84 was used to study the release of P_i A comparison of this conformation with the model for P_i release in PiPT revealed that Tyr¹⁷⁹ in Pho84 (Tyr¹⁵⁰ in PiPT) is not part of the P_i binding site. This difference may be due to a lack of detailed information on the P_i release step in PiPT. Molecular dynamics simulations of Pho84 in which a residue adjacent to Tyr¹⁷⁹, Asp¹⁷⁸, is protonated revealed a conformational change in Pho84 from an open, inward-facing state to an occluded state. Tyr¹⁷⁹ then became part of the binding site as was observed in the PiPT crystal structure. The importance of Tyr¹⁷⁹ in regulating P_i release was supported by site-directed mutagenesis and transport assays. Using trehalase activity measurements, we demonstrated that the release of P_i is a critical step for transceptor signaling. Our results add to previous studies on PiPT, creating a more complete picture of the proton-coupled P_i transport cycle of a transceptor.

Emulating proton-induced conformational changes in the vesicular monoamine transporter VMAT2 by mutagenesis

Neurotransporters located in synaptic vesicles are essential for communication between nerve cells in a process mediated by neurotransmitters. Vesicular monoamine transporter (VMAT), a member of the largest superfamily of transporters, mediates transport of monoamines to synaptic vesicles and storage organelles in a process that involves exchange of two H⁺ per substrate. VMAT transport is inhibited by the competitive inhibitor reserpine, a second-line agent to treat hypertension, and by the noncompetitive inhibitor tetrabenazine, presently in use for symptomatic treatment of hyperkinetic disorders. During the transport cycle, VMAT is expected to occupy at least three different conformations: cytoplasm-facing, occluded, and lumen-facing. The lumen- to cytoplasm-facing transition, facilitated by protonation of at least one of the essential membrane-embedded carboxyls, generates a binding site for reserpine. Here we have identified residues in the cytoplasmic gate and show that mutations that disrupt the interactions in this gate also shift the equilibrium toward the cytoplasm-facing conformation, emulating the effect of protonation. These experiments provide significant insight into the role of proton translocation in the conformational dynamics of a mammalian H⁺-coupled antiporter, and also identify key aspects of the mode of action and binding of two potent inhibitors of VMAT2: reserpine binds the cytoplasm-facing conformation, and tetrabenazine binds the lumen-facing conformation.

Thermodynamics of Nanobody Binding to Lactose Permease

Camelid nanobodies (Nbs) raised against the outward-facing conformer of a double-Trp mutant of the lactose permease of Escherichia coli (LacY) stabilize the permease in outward-facing conformations. Isothermal titration calorimetry is applied herein to dissect the binding thermodynamics of two Nbs, one that markedly improves access to the sugar-binding site and another that dramatically increases the affinity for galactoside. The findings presented here show that both enthalpy and entropy contribute favorably to binding of the Nbs to wild-type (WT) LacY and that binding of Nb to double-Trp mutant G46W/G262W is driven by a greater enthalpy at an entropic penalty. Thermodynamic analyses support the interpretation that WT LacY is stabilized in outward-facing conformations like the double-Trp mutant with closure of the water-filled cytoplasmic cavity through conformational selection. The LacY conformational transition required for ligand binding is reflected by a favorable entropy increase. Molecular dynamics simulations further suggest that the entropy increase likely stems from release of immobilized water molecules primarily from the cytoplasmic cavity upon closure.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

Coarse-grained simulations of proton-dependent conformational changes in lactose permease

During lactose/H(+) symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward-facing (open to cytoplasmic side) and outward-facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large-scale structural changes are not well understood. All-atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united-atom protein models with the coarse-grained MARTINI water/lipid, to investigate the proton-dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt-bridge dynamics agreed with all-atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward-facing to outward-facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter-helical distances measured in our simulations were consistent with the double electron-electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross-domain salt-bridge (Glu130-Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067-1074. © 2016 Wiley Periodicals, Inc.

Structure and mechanism of the mammalian fructose transporter GLUT5

The altered activity of the fructose transporter GLUT5, an isoform of the facilitated-diffusion glucose transporter family, has been linked to disorders such as type 2 diabetes and obesity. GLUT5 is also overexpressed in certain tumor cells and inhibitors are potential drugs for these conditions. Here, we describe the crystal structure of GLUT5 from Rattus norvegicus and Bos taurus in open outward- and open inward-facing conformations, respectively. GLUT5 has a major facilitator superfamily fold like other homologous monosaccharide transporters. Based on a comparison of the inward-facing structures of GLUT5 and human GLUT1, a ubiquitous glucose transporter, we show that a single point mutation is enough to switch the substrate binding preference of GLUT5 from fructose to glucose. A comparison of the substrate-free structures of GLUT5 with occluded substrate-bound structures of XylE suggests that, besides global rocker-switch like re-orientation of the bundles, local asymmetric rearrangements of C-terminal bundle helices TMs 7 and 10 underlie a “gated-pore” transport mechanism in such monosaccharide transporters.

Conformational gating dynamics in the GluCl anion-selective chloride channel

Cys-loop receptors are central to propagation of signals in the nervous system. The gating of the membrane-spanning pore is triggered by structural rearrangements in the agonist-binding site, located some 50 Å away from the pore. A sequential conformational change, propagating from the ligand-binding site to the pore, has been proposed to govern gating in all Cys-loop receptors. Here, we identify structural and dynamic components of the conformational gating in the eukaryotic glutamate-gated chloride channel (GluCl) by means of molecular dynamics (MD) simulations with and without the l-glutamate agonist bound. A significant increase in pore opening and accompanying hydration is observed in the presence of glutamate. Potential of mean force calculations reveal that the barrier for ion passage drops from 15 kcal/mol to 5-10 kcal/mol with the agonist bound. This appears to be explained by agonist binding that leads to significant changes in the intersubunit hydrogen-bonding pattern, which induce a slight tilt of the extracellular domain relative to the transmembrane domain in the simulations. This rearrangement is subtle, but correspond to the direction of the quaternary twist observed as a key difference between open and closed X-ray structures. While the full reversible gating is still a much slower process, the observed structural dynamics sheds new light on the early stages of how the agonist influences the extracellular domain, how the extracellular domain interacts with the transmembrane domain, and how changes in the transmembrane domain alter the free energy of ion passage.

Probing the importance of lipid diversity in cell membranes via molecular simulation

Lipid membranes in prokaryotes and eukaryotes have a wide array of lipids that are necessary for proper membrane structure and function. In this paper, an introduction to lipid diversity in biology and a mini-review on how molecular simulations have been used to model biological membranes (primarily limited to one to three lipid types in most simulation-based models) is provided, which motivates the use of all-atom molecular dynamics (MD) simulations to study the effect of lipid diversity on properties of realistic membrane models of prokaryotes and eukaryotes. As an example, cytoplasmic membrane models of Escherichia coli were developed at different stages of the colony growth cycle (early-log, mid-log, stationary and overnight). The main difference between lipid compositions at each stage was the concentration of a cyclopropane-containing moiety on the sn-2 lipid acyl chain (cyC17:0). Triplicate MD simulations for each stage were run for 300 ns to study the influence of lipid diversity on the surface area per lipid, area compressibility modulus, deuterium order parameters, and electron density profiles. The overnight stage (also known as the death stage) had the highest average surface area per lipid, highest rigidity, and lowest bilayer thickness compare to other stages of E. coli cytoplasmic membrane. Although bilayer thickness did depend on the growth stage, the changes between these were small suggesting that the hydrophobic core of transmembrane proteins fit well with the membrane in all growth stages. Although it is still common practise in MD simulations of membrane proteins to use simple one- or two-component membranes, it can be important to use diverse lipid model membranes when membrane protein structure and function are influenced by changes in lipid membrane composition.

Computational characterization of structural dynamics underlying function in active membrane transporters

Active transport of materials across the cellular membrane is one the most fundamental processes in biology. In order to accomplish this task, membrane transporters rely on a wide range of conformational changes spanning multiple time and size scales. These molecular events govern key functional aspects in membrane transporters, namely, coordinated gating motions underlying the alternating access mode of operation, and coupling of uphill transport of substrate to various sources of energy, for example, transmembrane electrochemical gradients and ATP binding and hydrolysis. Computational techniques such as molecular dynamics simulations and free energy calculations have equipped us with a powerful repertoire of biophysical tools offering unparalleled spatial and temporal resolutions that can effectively complement experimental methodologies, and therefore help fill the gap of knowledge in understanding the molecular basis of function in membrane transporters.

A chemiosmotic mechanism of symport

Lactose permease (LacY), a paradigm for the largest family of membrane transport proteins, catalyzes the coupled translocation of a galactoside and an H(+) across the Escherichia coli membrane (galactoside/H(+) symport). Initial X-ray structures reveal N- and C-terminal domains, each with six largely irregular transmembrane helices surrounding an aqueous cavity open to the cytoplasm. Recently, a structure with a narrow periplasmic opening and an occluded galactoside was obtained, confirming many observations and indicating that sugar binding involves induced fit. LacY catalyzes symport by an alternating access mechanism. Experimental findings garnered over 45 y indicate the following: (i) The limiting step for lactose/H(+) symport in the absence of the H(+) electrochemical gradient (∆µ̃H+) is deprotonation, whereas in the presence of ∆µ̃H+, the limiting step is opening of apo LacY on the other side of the membrane; (ii) LacY must be protonated to bind galactoside (the pK for binding is ∼10.5); (iii) galactoside binding and dissociation, not ∆µ̃H+, are the driving forces for alternating access; (iv) galactoside binding involves induced fit, causing transition to an occluded intermediate that undergoes alternating access; (v) galactoside dissociates, releasing the energy of binding; and (vi) Arg302 comes into proximity with protonated Glu325, causing deprotonation. Accumulation of galactoside against a concentration gradient does not involve a change in Kd for sugar on either side of the membrane, but the pKa (the affinity for H(+)) decreases markedly. Thus, transport is driven chemiosmotically but, contrary to expectation, ∆µ̃H+ acts kinetically to control the rate of the process.

Function, Structure, and Evolution of the Major Facilitator Superfamily: The LacY Manifesto

The major facilitator superfamily (MFS) is a diverse group of secondary transporters with members found in all kingdoms of life. A paradigm for MFS is the lactose permease (LacY) of Escherichia coli, which couples the stoichiometric translocation of a galactopyranoside and an H+ across the cytoplasmic membrane. LacY has been the test bed for the development of many methods applied for the analysis of transport proteins. X-ray structures of an inward-facing conformation and the most recent structure of an almost occluded conformation confirm many conclusions from previous studies. Although structure models are critical, they are insufficient to explain the catalysis of transport. The clues to understanding transport are based on the principles of enzyme kinetics. Secondary transport is a dynamic process—static snapshots of X-ray crystallography describe it only partially. However, without structural information, the underlying chemistry is virtually impossible to conclude. A large body of biochemical/biophysical data derived from systematic studies of site-directed mutants in LacY suggests residues critically involved in the catalysis, and a working model for the symport mechanism that involves alternating access of the binding site is presented. The general concepts derived from the bacterial LacY are examined for their relevance to other MFS transporters.

Flexible gates generate occluded intermediates in the transport cycle of LacY

The major facilitator superfamily (MFS) transporter lactose permease (LacY) alternates between cytoplasmic and periplasmic open conformations to co-transport a sugar molecule together with a proton across the plasma membrane. Indirect experimental evidence suggested the existence of an occluded transition intermediate of LacY, which would prevent leaking of the proton gradient. As no experimental structure is known, the conformational transition is not fully understood in atomic detail. We simulated transition events from a cytoplasmic open conformation to a periplasmic open conformation with the dynamic importance sampling molecular dynamics method and observed occluded intermediates. Analysis of water permeation pathways and the electrostatic free-energy landscape of a solvated proton indicated that the occluded state contains a solvated central cavity inaccessible from either side of the membrane. We propose a pair of geometric order parameters that capture the state of the pathway through the MFS transporters as shown by a survey of available crystal structures and models. We present a model for the occluded state of apo-LacY, which is similar to the occluded crystal structures of the MFS transporters EmrD, PepT_So, NarU, PiPT and XylE. Our simulations are consistent with experimental double electron spin–spin distance measurements that have been interpreted to show occluded conformations. During the simulations, a salt bridge that has been postulated to be involved in driving the conformational transition formed. Our results argue against a simple rigid-body domain motion as implied by a strict “rocker-switch mechanism” and instead hint at an intricate coupling between two flexible gates.

•

The transport mechanism of LacY is hypothesized to involve an intermediate “occluded” state.

•

Such a state is observed in computer simulations of the conformational transitions.

•

Simulation data are validated with experimental double electron–electron spin resonance measurements.

•

The structural gating elements of LacY are identified.

•

Occluded LacY is similar to known occluded structures of homologous proteins.

Evolutionary mix-and-match with MFS transporters II

One fundamentally important problem for understanding the mechanism of coupling between substrate and H(+) translocation with secondary active transport proteins is the identification and physical localization of residues involved in substrate and H(+) binding. This information is exceptionally difficult to obtain with the Major Facilitator Superfamily (MFS) because of the broad sequence diversity of the members. The MFS is the largest and most diverse group of transporters, many of which are clinically important, and includes members from all kingdoms of life. A wide range of substrates is transported, in many instances against a concentration gradient by transduction of the energy stored in an H(+) electrochemical gradient using symport mechanisms, which are discussed herein. Crystallographic structures of MFS members indicate that a deep central hydrophilic cavity surrounded by 12 mostly irregular transmembrane helices represents a common structural feature. An inverted triple-helix structural symmetry motif within the N- and C-terminal six-helix bundles suggests that the proteins may have arisen by intragenic multiplication. In the work presented here, the triple-helix motifs are aligned in combinatorial fashion so as to detect functionally homologous positions with known atomic structures of MFS members. Substrate and H(+)-binding sites in symporters that transport substrates, ranging from simple ions like phosphate to more complex peptides or disaccharides, are found to be in similar locations. It also appears likely that there is a homologous ordered kinetic mechanism for the H(+)-coupled MFS symporters.

Channelrhodopsins: a bioinformatics perspective

Channelrhodopsins are microbial-type rhodopsins that function as light-gated cation channels. Understanding how the detailed architecture of the protein governs its dynamics and specificity for ions is important, because it has the potential to assist in designing site-directed channelrhodopsin mutants for specific neurobiology applications. Here we use bioinformatics methods to derive accurate alignments of channelrhodopsin sequences, assess the sequence conservation patterns and find conserved motifs in channelrhodopsins, and use homology modeling to construct three-dimensional structural models of channelrhodopsins. The analyses reveal that helices C and D of channelrhodopsins contain Cys, Ser, and Thr groups that can engage in both intra- and inter-helical hydrogen bonds. We propose that these polar groups participate in inter-helical hydrogen-bonding clusters important for the protein conformational dynamics and for the local water interactions. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

Structural advances for the major facilitator superfamily (MFS) transporters

The major facilitator superfamily (MFS) is one of the largest groups of secondary active transporters conserved from bacteria to humans. MFS proteins selectively transport a wide spectrum of substrates across biomembranes and play a pivotal role in multiple physiological processes. Despite intense investigation, only seven MFS proteins from six subfamilies have been structurally elucidated. These structures were captured in distinct states during a transport cycle involving alternating access to binding sites from either side of the membrane. This review discusses recent progress in MFS structure analysis and focuses on the molecular basis for substrate binding, co-transport coupling, and alternating access.

Proper fatty acid composition rather than an ionizable lipid amine is required for full transport function of lactose permease from Escherichia coli

Energy-dependent uphill transport but not energy-independent downhill transport by lactose permease (LacY) is impaired when expressed in Escherichia coli cells or reconstituted in liposomes lacking phosphatidylethanolamine (PE) and containing only anionic phospholipids. The absence of PE results in inversion of the N-terminal half and misfolding of periplasmic domain P7, which are required for uphill transport of substrates. Replacement of PE in vitro by lipids with no net charge (phosphatidylcholine (PC), monoglucosyl diacylglycerol (GlcDAG), or diglucosyl diacylglycerol (GlcGlcDAG)) supported wild type transmembrane topology of the N-terminal half of LacY. The restoration of uphill transport in vitro was dependent on LacY native topology and proper folding of P7. Support of uphill transport by net neutral lipids in vitro (PE > PC ≫ GlcDAG ≠ GlcGlcDAG provided that PE or PC contained one saturated fatty acid) paralleled the results observed previously in vivo (PE = PC > GlcDAG ≠ GlcGlcDAG). Therefore, a free amino group is not required for uphill transport as previously concluded based on the lack of in vitro uphill transport when fully unsaturated PC replaced E. coli-derived PE. A close correlation was observed in vivo and in vitro between the ability of LacY to carry out uphill transport, the native conformation of P7, and the lipid headgroup and fatty acid composition. Therefore, the headgroup and the fatty acid composition of lipids are important for defining LacY topological organization and catalytically important structural features, further illustrating the direct role of lipids, independent of other cellular factors, in defining membrane protein structure/function.