The cost of living in the membrane: a case study of hydrophobic mismatch for the multi-segment protein LeuT

Functional regulation of aquaporin dynamics by lipid bilayer composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Mechanisms underlying drug-mediated regulation of membrane protein function

The hydrophobic coupling between membrane proteins and their host lipid bilayer provides a mechanism by which bilayer-modifying drugs may alter protein function. Drug regulation of membrane protein function thus may be mediated by both direct interactions with the protein and drug-induced alterations of bilayer properties, in which the latter will alter the energetics of protein conformational changes. To tease apart these mechanisms, we examine how the prototypical, proton-gated bacterial potassium channel KcsA is regulated by bilayer-modifying drugs using a fluorescence-based approach to quantify changes in both KcsA function and lipid bilayer properties (using gramicidin channels as probes). All tested drugs inhibited KcsA activity, and the changes in the different gating steps varied with bilayer thickness, suggesting a coupling to the bilayer. Examining the correlations between changes in KcsA gating steps and bilayer properties reveals that drug-induced regulation of membrane protein function indeed involves bilayer-mediated mechanisms. Both direct, either specific or nonspecific, binding and bilayer-mediated mechanisms therefore are likely to be important whenever there is overlap between the concentration ranges at which a drug alters membrane protein function and bilayer properties. Because changes in bilayer properties will impact many diverse membrane proteins, they may cause indiscriminate changes in protein function.

Unusual mode of dimerization of retinitis pigmentosa-associated F220C rhodopsin

Mutations in the G protein-coupled receptor (GPCR) rhodopsin are a common cause of autosomal dominant retinitis pigmentosa, a blinding disease. Rhodopsin self-associates in the membrane, and the purified monomeric apo-protein opsin dimerizes in vitro as it transitions from detergent micelles to reconstitute into a lipid bilayer. We previously reported that the retinitis pigmentosa-linked F220C opsin mutant fails to dimerize in vitro, reconstituting as a monomer. Using fluorescence-based assays and molecular dynamics simulations we now report that whereas wild-type and F220C opsin display distinct dimerization propensities in vitro as previously shown, they both dimerize in the plasma membrane of HEK293 cells. Unexpectedly, molecular dynamics simulations show that F220C opsin forms an energetically favored dimer in the membrane when compared with the wild-type protein. The conformation of the F220C dimer is unique, with transmembrane helices 5 and 6 splayed apart, promoting widening of the intracellular vestibule of each protomer and influx of water into the protein interior. FRET experiments with SNAP-tagged wild-type and F220C opsin expressed in HEK293 cells are consistent with this conformational difference. We speculate that the unusual mode of dimerization of F220C opsin in the membrane may have physiological consequences.

Confinement in Nanodiscs Anisotropically Modifies Lipid Bilayer Elastic Properties

Lipid nanodiscs are small synthetic lipid bilayer structures that are stabilized in solution by special circumscribing (or scaffolding) proteins or polymers. Because they create native-like environments for transmembrane proteins, lipid nanodiscs have become a powerful tool for structural determination of this class of systems when combined with cryo-electron microscopy or nuclear magnetic resonance. The elastic properties of lipid bilayers determine how the lipid environment responds to membrane protein perturbations, and how the lipid in turn modifies the conformational state of the embedded protein. However, despite the abundant use of nanodiscs in determining membrane protein structure, the elastic material properties of even pure lipid nanodiscs (i.e., without embedded proteins) have not yet been quantitatively investigated. A major hurdle is due to the inherently nonlocal treatment of the elastic properties of lipid systems implemented by most existing methods, both experimental and computational. In addition, these methods are best suited for very large “infinite” size lipidic assemblies, or ones that contain periodicity, in the case of simulations. We have previously described a computational analysis of molecular dynamics simulations designed to overcome these limitations, so it allows quantification of the bending rigidity (K_C) and tilt modulus (κ_t) on a local scale even for finite, nonperiodic systems, such as lipid nanodiscs. Here we use this computational approach to extract values of K_C and κ_t for a set of lipid nanodisc systems that vary in size and lipid composition. We find that the material properties of lipid nanodiscs are different from those of infinite bilayers of corresponding lipid composition, highlighting the effect of nanodisc confinement. Nanodiscs tend to show higher stiffness than their corresponding macroscopic bilayers, and moreover, their material properties vary spatially within them. For small-size MSP1 nanodiscs, the stiffness decreases radially, from a value that is larger in their center than the moduli of the corresponding bilayers by a factor of ∼2–3. The larger nanodiscs (MSP1E3D1 and MSP2N2) show milder spatial changes of moduli that are composition dependent and can be maximal in the center or at some distance from it. These trends in moduli correlate with spatially varying structural properties, including the area per lipid and the nanodisc thickness. Finally, as has previously been reported, nanodiscs tend to show deformations from perfectly flat circular geometries to varying degrees, depending on size and lipid composition. The modulations of lipid elastic properties that we find should be carefully considered when making structural and functional inferences concerning embedded proteins.

Membrane lipids are both the substrates and a mechanistically responsive environment of TMEM16 scramblase proteins

Recent discoveries about functional mechanisms of proteins in the TMEM16 family of phospholipid scramblases have illuminated the dual role of the membrane as both the substrate and a mechanistically responsive environment in the wide range of physiological processes and genetic disorders in which they are implicated. This is highlighted in the review of recent findings from our collaborative investigations of molecular mechanisms of TMEM16 scramblases that emerged from iterative functional, structural, and computational experimentation. In the context of this review, we present new MD simulations and trajectory analyses motivated by the fact that new structural information about the TMEM16 scramblases is emerging from cryo-EM determinations in lipid nanodiscs. Because the functional environment of these proteins in in vivo and in in vitro is closer to flat membranes, we studied comparatively the responses of the membrane to the TMEM16 proteins in flat membranes and nanodiscs. We find that bilayer shapes in the nanodiscs are very different from those observed in the flat membrane systems, but the function-related slanting of the membrane observed at the nhTMEM16 boundary with the protein is similar in the nanodiscs and in the flat bilayers. This changes, however, in the bilayer composed of longer-tail lipids, which is thicker near the phospholipid translocation pathway, which may reflect an enhanced tendency of the long tails to penetrate the pathway and create, as shown previously, a nonconductive environment. These findings support the correspondence between the mechanistic involvement of the lipid environment in the flat membranes, and the nanodiscs. © 2019 Wiley Periodicals, Inc.

The fine art of preparing membrane transport proteins for biomolecular simulations: Concepts and practical considerations

Molecular dynamics (MD) simulations have developed into an invaluable tool in bimolecular research, due to the capability of the method in capturing molecular events and structural transitions that describe the function as well as the physiochemical properties of biomolecular systems. Due to the progressive development of more efficient algorithms, expansion of the available computational resources, as well as the emergence of more advanced methodologies, the scope of computational studies has increased vastly over time. We now have access to a multitude of online databases, software packages, larger molecular systems and novel ligands due to the phenomenon of emerging novel psychoactive substances (NPS). With so many advances in the field, it is understandable that novices will no doubt find it challenging setting up a protein-ligand system even before they run their first MD simulation. These initial steps, such as homology modelling, ligand docking, parameterization, protein preparation and membrane setup have become a fundamental part of the drug discovery pipeline, and many areas of biomolecular sciences benefit from the applications provided by these technologies. However, there still remains no standard on their usage. Therefore, our aim within this review is to provide a clear overview of a variety of concepts and methodologies to consider, providing a workflow for a case study of a membrane transport protein, the full-length human dopamine transporter (hDAT) in complex with different stimulants, where MD simulations have recently been applied successfully.

Molecular Basis of the Glucose Transport Mechanism in Plants

The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane (SWEET, Sugars Will Eventually Be Exported Transporter). However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in the OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as for the glucose transport mechanism. SWEETs undergo a structural transition to outward-facing (OF), occluded (OC), and inward-facing (IF) and strongly support an alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly the same for OF, OC, and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and acts as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass-like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

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.

The Role of Lipid Interactions in Simulations of the α-Hemolysin Ion-Channel-Forming Toxin

Molecular dynamics simulations were performed to describe the function of the ion-channel-forming toxin α-hemolysin (αHL) in lipid membranes that were composed of either 1,2-diphytanoyl-sn-glycero-3-phospho-choline or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline. The simulations highlight the importance of lipid type in maintaining αHL structure and function, enabling direct comparison to experiments for biosensing applications. We determined that although the two lipids studied are similar in structure, 1,2-diphytanoyl-sn-glycero-3-phospho-choline membranes better match the hydrophobic thickness of αHL compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline membranes. This hydrophobic match is essential to maintaining proper alignment of β-sheet loops at the trans entrance of αHL, which, when disrupted, creates an additional constriction to ion flow that decreases the channel current below experimental values and creates greater variability in channel conductance. Agreement with experiments was further improved with sufficient lipid membrane equilibration and allowed the discrimination of subtle αHL conduction states with lipid type. Finally, we explore the effects of truncating the extramembrane cap of αHL and its role in maintaining proper alignment of αHL in the membrane and channel conductance. Our results demonstrate the essential role of lipid type and lipid-protein interactions in simulations of αHL and will considerably improve the interpretation of experimental data.

A partially-open inward-facing intermediate conformation of LeuT is associated with Na+ release and substrate transport

Neurotransmitter:sodium symporters (NSS), targets of antidepressants and psychostimulants, clear neurotransmitters from the synaptic cleft through sodium (Na+)-coupled transport. Substrate and Na+ are thought to be transported from the extracellular to intracellular space through an alternating access mechanism by coordinated conformational rearrangements in the symporter that alternately expose the binding sites to each side of the membrane. However, the mechanism by which the binding of ligands coordinates conformational changes occurring on opposite sides of the membrane is not well understood. Here, we report the use of single-molecule fluorescence resonance energy transfer (smFRET) techniques to image transitions between distinct conformational states on both the extracellular and intracellular sides of the prokaryotic NSS LeuT, including partially open intermediates associated with transport activity. The nature and functional context of these hitherto unidentified intermediate states shed new light on the allosteric mechanism that couples substrate and Na+ symport by the NSS family through conformational dynamics.

New Continuum Approaches for Determining Protein-Induced Membrane Deformations

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

Grafting Charged Species to Membrane-Embedded Scaffolds Dramatically Increases the Rate of Bilayer Flipping

The cell membrane is a barrier to the passive diffusion of charged molecules due to the chemical properties of the lipid bilayer. Surprisingly, recent experiments have identified processes in which synthetic and biological charged species directly transfer across lipid bilayers on biologically relevant time scales. In particular, amphiphilic nanoparticles have been shown to insert into lipid bilayers, requiring the transport of charged species across the bilayer. The molecular factors facilitating this rapid insertion process remain unknown. In this work, we use atomistic molecular dynamics simulations to calculate the free energy barrier associated with "flipping" charged species across a lipid bilayer for species that are grafted to a membrane-embedded scaffold, such as a membrane-embedded nanoparticle. We find that the free energy barrier for flipping a grafted ligand can be over 7 kcal/mol lower than the barrier for translocating an isolated, equivalent ion, yielding a 5 order of magnitude decrease in the corresponding flipping time scale. Similar results are found for flipping charged species grafted to either nanoparticle or protein scaffolds. These results reveal new mechanistic insight into the flipping of charged macromolecular components that might play an important, yet overlooked, role in signaling and charge transport in biological settings. Furthermore, our results suggest guidelines for the design of synthetic materials capable of rapidly flipping charged moieties across the cell membrane.

What Can and Cannot Be Learned from Molecular Dynamics Simulations of Bacterial Proton-Coupled Oligopeptide Transporter GkPOT?

We have performed an extensive set of all-atom molecular dynamics (MD) simulations of a bacterial proton-coupled oligopeptide transporter (POT) in an explicit membrane environment. We have characterized both the local and global conformational dynamics of the transporter upon the proton and/or substrate binding, within a statistical framework. Our results reveal a clearly distinct behavior for local conformational dynamics in the absence and presence of the proton at the putative proton binding residue E310. Particularly, we find that the substrate binding conformation is drastically different in the two conditions, where the substrate binds to the protein in a lateral/vertical manner, in the presence/absence of the proton. We do not observe any statistically significant distinctive behavior in terms of the global conformational changes in different simulation conditions, within the time scales of our simulations. Our extensive simulations and analyses call into question the implicit assumption of many MD studies that local conformational changes observed in short simulations could provide clues to the global conformational changes that occur on much longer time scales. The linear regression analysis of quantities associated with the global conformational fluctuations, however, provides an indication of a mechanism involving the concerted motion of the transmembrane helices, consistent with the rocker-switch mechanism.

The Environment Shapes the Inner Vestibule of LeuT

Human neurotransmitter transporters are found in the nervous system terminating synaptic signals by rapid removal of neurotransmitter molecules from the synaptic cleft. The homologous transporter LeuT, found in Aquifex aeolicus, was crystallized in different conformations. Here, we investigated the inward-open state of LeuT. We compared LeuT in membranes and micelles using molecular dynamics simulations and lanthanide-based resonance energy transfer (LRET). Simulations of micelle-solubilized LeuT revealed a stable and widely open inward-facing conformation. However, this conformation was unstable in a membrane environment. The helix dipole and the charged amino acid of the first transmembrane helix (TM1A) partitioned out of the hydrophobic membrane core. Free energy calculations showed that movement of TM1A by 0.30 nm was driven by a free energy difference of ~15 kJ/mol. Distance measurements by LRET showed TM1A movements, consistent with the simulations, confirming a substantially different inward-open conformation in lipid bilayer from that inferred from the crystal structure.

Crystal structures of the bacterial small amino acid transporter LeuT provided structural evidence for the alternating access model. Thereby, these structures shaped our understanding of the mechanisms underlying substrate translocation by neurotransmitter transporters. However, it has been questioned, if the crystallized inward-open conformation of LeuT can exist in the membrane environment. Here we show that, while stable in detergent micelles, the inward-open conformation of LeuT is of high energy and undergoes structural readjustments. We use a multi-faceted approach including molecular dynamics simulations, scintillation proximity assays, free energy calculations and apply for the first time lanthanide resonance energy transfer measurements to verify the in silico predictions. In silico and in vitro approaches using the same conditions allowed us to combine the macroscopic experimental data with microscopic all atom results from simulations to identify the underlying driving forces: partitioning of charged and polar groups from the hydrophobic membrane interior to the hydrophilic environment. We propose that the inward-facing state shows a much smaller movement of TM1A, but large enough to create an access path to the S1 substrate binding site from the vestibule.

Conformational Dynamics on the Extracellular Side of LeuT Controlled by Na+ and K+ Ions and the Protonation State of Glu290

Ions play key mechanistic roles in the gating dynamics of neurotransmitter:sodium symporters (NSSs). In recent microsecond scale molecular dynamics simulations of a complete model of the dopamine transporter, a NSS protein, we observed a partitioning of K(+) ions from the intracellular side toward the unoccupied Na2 site of dopamine transporter following the release of the Na2-bound Na(+) Here we evaluate with computational simulations and experimental measurements of ion affinities under corresponding conditions, the consequences of K(+) binding in the Na2 site of LeuT, a bacterial homolog of NSS, when both Na(+) ions and substrate have left, and the transporter prepares for a new cycle. We compare the results with the consequences of binding Na(+) in the same apo system. Analysis of >50-μs atomistic molecular dynamics and enhanced sampling trajectories of constructs with Glu(290), either charged or neutral, point to the Glu(290) protonation state as a main determinant in the structural reconfiguration of the extracellular vestibule of LeuT in which a "water gate" opens through coordinated motions of residues Leu(25), Tyr(108), and Phe(253) The resulting water channel enables the binding/dissociation of the Na(+) and K(+) ions that are prevalent, respectively, in the extracellular and intracellular environments.

Continuum descriptions of membranes and their interaction with proteins: Towards chemically accurate models

Biological membranes deform in response to resident proteins leading to a coupling between membrane shape and protein localization. Additionally, the membrane influences the function of membrane proteins. Here we review contributions to this field from continuum elastic membrane models focusing on the class of models that couple the protein to the membrane. While it has been argued that continuum models cannot reproduce the distortions observed in fully-atomistic molecular dynamics simulations, we suggest that this failure can be overcome by using chemically accurate representations of the protein. We outline our recent advances along these lines with our hybrid continuum-atomistic model, and we show the model is in excellent agreement with fully-atomistic simulations of the nhTMEM16 lipid scramblase. We believe that the speed and accuracy of continuum-atomistic methodologies will make it possible to simulate large scale, slow biological processes, such as membrane morphological changes, that are currently beyond the scope of other computational approaches. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Properties of an inward-facing state of LeuT: conformational stability and substrate release

The leucine transporter (LeuT) is a bacterial homolog of the human monoamine transporters, which are important pharmaceutical targets. There are no high-resolution structures of the human transporters available; however, LeuT has been crystallized in several different conformational states. Recently, an inward-facing conformation of LeuT was solved revealing an unexpectedly large movement of transmembrane helix 1a (TM1a). We have performed molecular dynamics simulations of the mutated and wild-type transporter, with and without the cocrystallized Fab antibody fragment, to investigate the properties of this inward-facing conformation in relation to transport by LeuT within the membrane environment. In all of the simulations, local conformational changes with respect to the crystal structure are consistently observed, especially in TM1a. Umbrella sampling revealed a soft potential for TM1a tilting. Furthermore, simulations of inward-facing LeuT with Na(+) ions and substrate bound suggest that one of the Na(+) ion binding sites is fully disrupted. Release of alanine and the second Na(+) ion is also observed, giving insight into the final stage of the translocation process in atomistic detail.

Allosteric Mechanisms of Molecular Machines at the Membrane: Transport by Sodium-Coupled Symporters

Solute transport across cell membranes is ubiquitous in biology as an essential physiological process. Secondary active transporters couple the unfavorable process of solute transport against its concentration gradient to the energetically favorable transport of one or several ions. The study of such transporters over several decades indicates that their function involves complex allosteric mechanisms that are progressively being revealed in atomistic detail. We focus on two well-characterized sodium-coupled symporters: the bacterial amino acid transporter LeuT, which is the prototype for the "gated pore" mechanism in the mammalian synaptic monoamine transporters, and the archaeal GltPh, which is the prototype for the "elevator" mechanism in the mammalian excitatory amino acid transporters. We present the evidence for the role of allostery in the context of a quantitative formalism that can reconcile biochemical and biophysical data and thereby connects directly to recent insights into the molecular structure and dynamics of these proteins. We demonstrate that, while the structures and mechanisms of these transporters are very different, the available data suggest a common role of specific models of allostery in their functions. We argue that such allosteric mechanisms appear essential not only for sodium-coupled symport in general but also for the function of other types of molecular machines in the membrane.

Role of Annular Lipids in the Functional Properties of Leucine Transporter LeuT Proteomicelles

Recent work has shown that the choice of the type and concentration of detergent used for the solubilization of membrane proteins can strongly influence the results of functional experiments. In particular, the amino acid transporter LeuT can bind two substrate molecules in low concentrations of n-dodecyl β-d-maltopyranoside (DDM), whereas high concentrations reduce the molar binding stoichiometry to 1:1. Subsequent molecular dynamics (MD) simulations of LeuT in DDM proteomicelles revealed that DDM can penetrate to the extracellular vestibule and make stable contacts in the functionally important secondary substrate binding site (S2), suggesting a potential competitive mechanism for the reduction in binding stoichiometry. Because annular lipids can be retained during solubilization, we performed MD simulations of LeuT proteomicelles at various stages of the solubilization process. We find that at low DDM concentrations, lipids are retained around the protein and penetration of detergent into the S2 site does not occur, whereas at high concentrations, lipids are displaced and the probability of DDM binding in the S2 site is increased. This behavior is dependent on the type of detergent, however, as we find in the simulations that the detergent lauryl maltose-neopentyl glycol, which is approximately twice the size of DDM and structurally more closely resembles lipids, does not penetrate the protein even at very high concentrations. We present functional studies that confirm the computational findings, emphasizing the need for careful consideration of experimental conditions, and for cautious interpretation of data in gathering mechanistic information about membrane proteins.

Membrane Protein Properties Revealed through Data-Rich Electrostatics Calculations

The electrostatic properties of membrane proteins often reveal many of their key biophysical characteristics, such as ion channel selectivity and the stability of charged membrane-spanning segments. The Poisson-Boltzmann (PB) equation is the gold standard for calculating protein electrostatics, and the software APBSmem enables the solution of the PB equation in the presence of a membrane. Here, we describe significant advances to APBSmem, including full automation of system setup, per-residue energy decomposition, incorporation of PDB2PQR, calculation of membrane-induced pKa shifts, calculation of non-polar energies, and command-line scripting for large-scale calculations. We highlight these new features with calculations carried out on a number of membrane proteins, including the recently solved structure of the ion channel TRPV1 and a large survey of 1,614 membrane proteins of known structure. This survey provides a comprehensive list of residues with large electrostatic penalties for being embedded in the membrane, potentially revealing interesting functional information.

MemProtMD: Automated Insertion of Membrane Protein Structures into Explicit Lipid Membranes

There has been exponential growth in the number of membrane protein structures determined. Nevertheless, these structures are usually resolved in the absence of their lipid environment. Coarse-grained molecular dynamics (CGMD) simulations enable insertion of membrane proteins into explicit models of lipid bilayers. We have automated the CGMD methodology, enabling membrane protein structures to be identified upon their release into the PDB and embedded into a membrane. The simulations are analyzed for protein-lipid interactions, identifying lipid binding sites, and revealing local bilayer deformations plus molecular access pathways within the membrane. The coarse-grained models of membrane protein/bilayer complexes are transformed to atomistic resolution for further analysis and simulation. Using this automated simulation pipeline, we have analyzed a number of recently determined membrane protein structures to predict their locations within a membrane, their lipid/protein interactions, and the functional implications of an enhanced understanding of the local membrane environment of each protein.

•

A simulation pipeline for predicting the location of a membrane protein in a bilayer

•

A protocol for identifying novel membrane protein structures in the PDB

•

Analysis of lipid binding sites and local bilayer deformation by membrane proteins

•

Functional implications from enhanced understanding of local membrane environments

Functional mechanisms of neurotransmitter transporters regulated by lipid-protein interactions of their terminal loops

The physiological functions of neurotransmitter:sodium symporters (NSS) in reuptake of neurotransmitters from the synapse into the presynaptic nerve have been shown to be complemented by their involvement, together with non-plasma membrane neurotransmitter transporters, in the reverse transport of substrate (efflux) in response to psychostimulants. Recent experimental evidence implicates highly anionic phosphatidylinositol 4,5-biphosphate (PIP(2)) lipids in such functions of the serotonin (SERT) and dopamine (DAT) transporters. Thus, for both SERT and DAT, neurotransmitter efflux has been shown to be strongly regulated by the presence of PIP(2) lipids in the plasma membrane, and the electrostatic interaction of the N-terminal region of DAT with the negatively charged PIP(2) lipids. We examine the experimentally established phenotypes in a structural context obtained from computational modeling based on recent crystallographic data. The results are shown to set the stage for a mechanistic understanding of physiological actions of neurotransmitter transporters in the NSS family of membrane proteins. This article is part of a Special Issue entitled: Lipid-protein interactions.

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.

Computational modeling of the N-terminus of the human dopamine transporter and its interaction with PIP2 -containing membranes

The dopamine transporter (DAT) is a transmembrane protein belonging to the family of neurotransmitter:sodium symporters (NSS). Members of the NSS are responsible for the clearance of neurotransmitters from the synaptic cleft, and for their translocation back into the presynaptic nerve terminal. The DAT contains long intracellular N- and C-terminal domains that are strongly implicated in the transporter function. The N-terminus (N-term), in particular, regulates the reverse transport (efflux) of the substrate through DAT. Currently, the molecular mechanisms of the efflux remain elusive in large part due to lack of structural information on the N-terminal segment. Here we report a computational model of the N-term of the human DAT (hDAT), obtained through an ab initio structure prediction, in combination with extensive atomistic molecular dynamics (MD) simulations in the context of a lipid membrane. Our analysis reveals that whereas the N-term is a highly dynamic domain, it contains secondary structure elements that remain stable in the long MD trajectories of interactions with the bilayer (totaling >2.2 μs). Combining MD simulations with continuum mean-field modeling we found that the N-term engages with lipid membranes through electrostatic interactions with the charged lipids PIP2 (phosphatidylinositol 4,5-Biphosphate) or PS (phosphatidylserine) that are present in these bilayers. We identify specific motifs along the N-term implicated in such interactions and show that differential modes of N-term/membrane association result in differential positioning of the structured segments on the membrane surface. These results will inform future structure-based studies that will elucidate the mechanistic role of the N-term in DAT function.

Transport domain unlocking sets the uptake rate of an aspartate transporter

Glutamate transporters terminate neurotransmission by clearing synaptically released glutamate from the extracellular space, allowing repeated rounds of signaling and preventing glutamate-mediated excitotoxicity. Crystallographic studies on an archaeal homologue, Glt_Ph, showed that distinct transport domains translocate substrates into the cytoplasm by moving across the membrane within a central trimerization scaffold. Here, we report direct observations of these 'elevator-like' transport domain motions in the context of reconstituted proteoliposomes and physiological ion gradients using single-molecule fluorescence resonance energy transfer (smFRET) imaging. We show that Glt_Ph bearing two “humanizing” mutations exhibits markedly increased transport domain dynamics, which parallels an increased rate of substrate transport, thereby establishing a direct temporal relationship between transport domain motions and substrate uptake. Crystallographic and computational investigations reveal that these mutations favor structurally “unlocked” states with increased solvent occupancy at the interface between the transport domain and the trimeric scaffold.

A structural model of the active ribosome-bound membrane protein insertase YidC

The integration of most membrane proteins into the cytoplasmic membrane of bacteria occurs co-translationally. The universally conserved YidC protein mediates this process either individually as a membrane protein insertase, or in concert with the SecY complex. Here, we present a structural model of YidC based on evolutionary co-variation analysis, lipid-versus-protein-exposure and molecular dynamics simulations. The model suggests a distinctive arrangement of the conserved five transmembrane domains and a helical hairpin between transmembrane segment 2 (TM2) and TM3 on the cytoplasmic membrane surface. The model was used for docking into a cryo-electron microscopy reconstruction of a translating YidC-ribosome complex carrying the YidC substrate FOc. This structure reveals how a single copy of YidC interacts with the ribosome at the ribosomal tunnel exit and identifies a site for membrane protein insertion at the YidC protein-lipid interface. Together, these data suggest a mechanism for the co-translational mode of YidC-mediated membrane protein insertion.

Simulating G protein-coupled receptors in native-like membranes: from monomers to oligomers

G protein-coupled receptors (GPCRs) are one of the most relevant superfamilies of transmembrane proteins as they participate in an important variety of biological events. Recently, the scientific community is witnessing an advent of a GPCR crystallization age along with impressive improvements achieved in the field of computer simulations during the last two decades. Computer simulation techniques such as molecular dynamics (MD) simulations are now frequent tools to study the dynamic behavior of GPCRs and, more importantly, to model the complex membrane environment where these proteins spend their lifetime. Thanks to these tools, GPCRs can be simulated not only longer but also in a more "physiological" fashion. In this scenario, scientists are taking advantage of such advances to approach certain phenomena such as GPCR oligomerization occurring only at timescales not reachable until now. Thus, despite current MD simulations having important limitations today, they have become an essential tool to study key biophysical properties of GPCRs and GPCR oligomers.

Not just an oil slick: how the energetics of protein-membrane interactions impacts the function and organization of transmembrane proteins

The membrane environment, its composition, dynamics, and remodeling, have been shown to participate in the function and organization of a wide variety of transmembrane (TM) proteins, making it necessary to study the molecular mechanisms of such proteins in the context of their membrane settings. We review some recent conceptual advances enabling such studies, and corresponding computational models and tools designed to facilitate the concerted experimental and computational investigation of protein-membrane interactions. To connect productively with the high resolution achieved by cognate experimental approaches, the computational methods must offer quantitative data at an atomistically detailed level. We show how such a quantitative method illuminated the mechanistic importance of a structural characteristic of multihelical TM proteins, that is, the likely presence of adjacent polar and hydrophobic residues at the protein-membrane interface. Such adjacency can preclude the complete alleviation of the well-known hydrophobic mismatch between TM proteins and the surrounding membrane, giving rise to an energy cost of residual hydrophobic mismatch. The energy cost and biophysical formulation of hydrophobic mismatch and residual hydrophobic mismatch are reviewed in the context of their mechanistic role in the function of prototypical members of multihelical TM protein families: 1), LeuT, a bacterial homolog of mammalian neurotransmitter sodium symporters; and 2), rhodopsin and the β1- and β2-adrenergic receptors from the G-protein coupled receptor family. The type of computational analysis provided by these examples is poised to translate the rapidly growing structural data for the many TM protein families that are of great importance to cell function into ever more incisive insights into mechanisms driven by protein-ligand and protein-protein interactions in the membrane environment.

Structural insights into functional lipid-protein interactions in secondary transporters

BACKGROUND

Structural evidences with functional corroborations have revealed distinct features of lipid-protein interactions especially in channels and receptors. Many membrane embedded transporters are also known to require specific lipids for their functions and for some of them cellular and biochemical data suggest tight regulation by the lipid bilayer. However, molecular details on lipid-protein interactions in transporters are sparse since lipids are either depleted from the detergent solubilized transporters in three-dimensional crystals or not readily resolved in crystal structures. Nevertheless the steady increase in the progress of transporter structure determination contributed more examples of structures with resolved lipids.

SCOPE OF REVIEW

This review gives an overview on transporter structures in complex with lipids reported to date and discusses commonly encountered difficulties in the identification of functionally significant lipid-protein interactions based on those structures and functional in vitro data. Recent structures provided molecular details into regulation mechanism of transporters by specific lipids. The review highlights common findings and conserved patterns for distantly related transporter families to draw a more general picture on the regulatory role of lipid-protein interactions.

MAJOR CONCLUSIONS

Several common themes of the manner in which lipids directly influence membrane-mediated folding, oligomerization and structure stability can be found. Especially for LeuT-like fold transporters similarities in structurally resolved lipid-protein interactions suggest a common way in which transporter conformations are affected by lipids even in evolutionarily distinct transporters. Lipids appear to play an additional role as joints mechanically reinforcing the inverted repeat topology, which is a major determinant in the alternating access mechanism of secondary transporters.

GENERAL SIGNIFICANCE

This review brings together and adds to the repertoire of knowledge on lipid-protein interactions of functional significance presented in structures of membrane transporters. Knowledge of specific lipid-binding sites and modes of lipid influence on these proteins not only accomplishes the molecular description of transport cycle further, but also sheds light into localization dependent differences of transporter function. This article is part of a Special Issue entitled Structural biochemistry and biophysics of membrane proteins.

NbIT--a new information theory-based analysis of allosteric mechanisms reveals residues that underlie function in the leucine transporter LeuT

Complex networks of interacting residues and microdomains in the structures of biomolecular systems underlie the reliable propagation of information from an input signal, such as the concentration of a ligand, to sites that generate the appropriate output signal, such as enzymatic activity. This information transduction often carries the signal across relatively large distances at the molecular scale in a form of allostery that is essential for the physiological functions performed by biomolecules. While allosteric behaviors have been documented from experiments and computation, the mechanism of this form of allostery proved difficult to identify at the molecular level. Here, we introduce a novel analysis framework, called N-body Information Theory (NbIT) analysis, which is based on information theory and uses measures of configurational entropy in a biomolecular system to identify microdomains and individual residues that act as (i)-channels for long-distance information sharing between functional sites, and (ii)-coordinators that organize dynamics within functional sites. Application of the new method to molecular dynamics (MD) trajectories of the occluded state of the bacterial leucine transporter LeuT identifies a channel of allosteric coupling between the functionally important intracellular gate and the substrate binding sites known to modulate it. NbIT analysis is shown also to differentiate residues involved primarily in stabilizing the functional sites, from those that contribute to allosteric couplings between sites. NbIT analysis of MD data thus reveals rigorous mechanistic elements of allostery underlying the dynamics of biomolecular systems.

We developed the new information theory-based analysis framework presented here, NbIT analysis, for the study of allosteric mechanisms in biomolecular systems from Molecular Dynamics trajectories. The illustrative application of NbIT to the analysis of the occluded state in the bacterial transporter LeuT, produced a quantitative representation of the allosteric behavior, and identified intramolecular channels that enable the long-distance information transmission. Our findings, identifying the roles of specific residues in the communication of the allosteric information, were validated by the recognition of residues that have been previously shown to play functional roles in this very well studied system. In addition, we show that application of NbIT analysis leads to the discrimination of functional roles by differentiating between residues that are essential to the dynamics within functional sites (e.g., the substrate binding sites), and residues whose role is to communicate between such functional sites. These results demonstrate that the information theoretical analysis presented here is a powerful tool for quantifying complex allosteric behavior in biomolecular systems and for identifying the crucial components underlying those behaviors.

Protein and lipid interactions driving molecular mechanisms of in meso crystallization

The recent advances in the in meso crystallization technique for the structural characterization of G-protein coupled receptor (GPCR) proteins have established the usefulness of the lipidic-cubic phases (LCPs) in the field of crystallography of membrane proteins. It is surprising that despite the success of the approach, the molecular mechanisms of the in meso method are still not well understood. Therefore, the approach must rely on extensive screening for a suitable protein construct, for host and additive lipids, and for the appropriate precipitants and temperature. To shed light on the in meso crystallization mechanisms, we used extensive coarse-grained molecular dynamics simulations to study, in molecular detail, LCPs under different conditions (compositions and temperatures relevant to crystallogenesis) and their interactions with different types of GPCR constructs. The results presented show how the modulation of the lattice constant of the LCP (triggered by the addition of precipitant during the in meso assay), or of the host lipid type, can destabilize monomeric proteins in the bilayer of the LCP and thus drive their aggregation into the stacked lamellae, where the residual hydrophobic mismatch between the protein and the membrane can drive the formation of lateral contacts leading to nucleation and crystal growth. Moreover, we demonstrate how particular protein designs (such as transmembrane proteins engineered to contain large polar regions) can promote protein stacking interactions in the third, out-of-plane, dimension. The insights provided by the new aspects of the specific molecular mechanisms responsible for protein-protein interactions inside the cubic phase presented here should be helpful in guiding the rational design of future in meso trials with successful outcomes.

Structural evidence for functional lipid interactions in the betaine transporter BetP

Bilayer lipids contribute to the stability of membrane transporters and are crucially involved in their proper functioning. However, the molecular knowledge of how surrounding lipids affect membrane transport is surprisingly limited and despite its general importance is rarely considered in the molecular description of a transport mechanism. One reason is that only few atomic resolution structures of channels or transporters reveal a functional interaction with lipids, which are difficult to detect in X-ray structures per se. Overcoming these difficulties, we report here on a new structure of the osmotic stress-regulated betaine transporter BetP in complex with anionic lipids. This lipid-associated BetP structure is important in the molecular understanding of osmoregulation due to the strong dependence of activity regulation in BetP on the presence of negatively charged lipids. We detected eight resolved palmitoyl-oleoyl phosphatidyl glycerol (PG) lipids mimicking parts of the membrane leaflets and interacting with key residues in transport and regulation. The lipid-protein interactions observed here in structural detail in BetP provide molecular insights into the role of lipids in osmoregulated secondary transport.

Membrane driven spatial organization of GPCRs

Spatial organization of G-protein coupled receptors (GPCRs) into dimers and higher order oligomers has been demonstrated in vitro and in vivo. The pharmacological readout was shown to depend on the specific interfaces, but why particular regions of the GPCR structure are involved, and how ligand-determined states change them remains unknown. Here we show why protein-membrane hydrophobic matching is attained upon oligomerization at specific interfaces from an analysis of coarse-grained molecular dynamics simulations of the spontaneous diffusion-interaction of the prototypical beta2-adrenergic (β₂AR) receptors in a POPC lipid bilayer. The energy penalty from mismatch is significantly reduced in the spontaneously emerging oligomeric arrays, making the spatial organization of the GPCRs dependent on the pattern of mismatch in the monomer. This mismatch pattern is very different for β₂AR compared to the highly homologous and structurally similar β₁AR, consonant with experimentally observed oligomerization patterns of β₂AR and β₁AR. The results provide a mechanistic understanding of the structural context of oligomerization.