Principles of membrane protein interactions with annular lipids deduced from aquaporin-0 2D crystals

Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts

Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.

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.

Solid-state NMR structure determination of a membrane protein in E. coli cellular inner membrane

Structure determination of membrane proteins in native cellular membranes is critical to precisely reveal their structures in physiological conditions. However, it remains challenging for solid-state nuclear magnetic resonance (ssNMR) due to the low sensitivity and high complexity of ssNMR spectra of cellular membranes. Here, we present the structure determination of aquaporin Z (AqpZ) by ssNMR in Escherichia coli inner membranes. To enhance the signal sensitivity of AqpZ, we optimized protein overexpression and removed outer membrane components. To suppress the interference of background proteins, we used a "dual-media" expression approach and antibiotic treatment. Using 1017 distance restraints obtained from two-dimensional 13C-13C spectra based on the complete chemical shift assignments, the 1.7-Å ssNMR structure of AqpZ is determined in E. coli inner membranes. This cellular ssNMR structure determination paves the way for analyzing the atomic structural details for membrane proteins in native cellular membranes.

Structural analysis of the water channel AQP2 by single-particle cryo-EM

Water channels, which are small membrane proteins almost entirely buried in lipid membranes, are challenging research targets for single-particle cryo-electron microscopy (cryo-EM), a powerful technique routinely used to determine the structures of membrane proteins. Because the single-particle method enables structural analysis of a whole protein with flexible parts that interfere with crystallization, we have focused our efforts on analyzing water channel structures. Here, utilizing this system, we analyzed the structure of full-length aquaporin-2 (AQP2), a primary regulator of vasopressin-dependent reabsorption of water at the renal collecting ducts. The 2.9 Å resolution map revealed a cytoplasmic extension of the cryo-EM density that was presumed to be the highly flexible C-terminus at which the localization of AQP2 is regulated in the renal collecting duct cells. We also observed a continuous density along the common water pathway inside the channel pore and lipid-like molecules at the membrane interface. Observations of these constructions in the AQP2 structure analyzed without any fiducial markers (e.g., a rigidly bound antibody) indicate that single-particle cryo-EM will be useful for investigating water channels in native states as well as in complexes with chemical compounds.

Mechanosensitive aquaporins

Cellular systems must deal with mechanical forces to satisfy their physiological functions. In this context, proteins with mechanosensitive properties play a crucial role in sensing and responding to environmental changes. The discovery of aquaporins (AQPs) marked a significant breakthrough in the study of water transport. Their transport capacity and regulation features make them key players in cellular processes. To date, few AQPs have been reported to be mechanosensitive. Like mechanosensitive ion channels, AQPs respond to tension changes in the same range. However, unlike ion channels, the aquaporin's transport rate decreases as tension increases, and the molecular features of the mechanism are unknown. Nevertheless, some clues from mechanosensitive ion channels shed light on the AQP-membrane interaction. The GxxxG motif may play a critical role in the water permeation process associated with structural features in AQPs. Consequently, a possible gating mechanism triggered by membrane tension changes would involve a conformational change in the cytoplasmic extreme of the single file region of the water pathway, where glycine and histidine residues from loop B play a key role. In view of their transport capacity and their involvement in relevant processes related to mechanical forces, mechanosensitive AQPs are a fundamental piece of the puzzle for understanding cellular responses.

Aquaporin Gating: A New Twist to Unravel Permeation through Water Channels

Aquaporins (AQPs) are small transmembrane tetrameric proteins that facilitate water, solute and gas exchange. Their presence has been extensively reported in the biological membranes of almost all living organisms. Although their discovery is much more recent than ion transport systems, different biophysical approaches have contributed to confirm that permeation through each monomer is consistent with closed and open states, introducing the term gating mechanism into the field. The study of AQPs in their native membrane or overexpressed in heterologous systems have experimentally demonstrated that water membrane permeability can be reversibly modified in response to specific modulators. For some regulation mechanisms, such as pH changes, evidence for gating is also supported by high-resolution structures of the water channel in different configurations as well as molecular dynamics simulation. Both experimental and simulation approaches sustain that the rearrangement of conserved residues contributes to occlude the cavity of the channel restricting water permeation. Interestingly, specific charged and conserved residues are present in the environment of the pore and, thus, the tetrameric structure can be subjected to alter the positions of these charges to sustain gating. Thus, is it possible to explore whether the displacement of these charges (gating current) leads to conformational changes? To our knowledge, this question has not yet been addressed at all. In this review, we intend to analyze the suitability of this proposal for the first time.

Native-like environments afford novel mechanistic insights into membrane proteins

Advances in cryogenic electron microscopy (cryo-EM) enabled routine near-atomic structure determination of membrane proteins, while nanodisc technology has provided a way to provide membrane proteins with a native or native-like lipid environment. After giving a brief history of membrane mimetics, we present example structures of membrane proteins in nanodiscs that revealed information not provided by structures obtained in detergent. We describe how the lipid environment surrounding the membrane protein can be custom designed during nanodisc assembly and how it can be modified after assembly to test functional hypotheses. Because nanodiscs most closely replicate the physiologic environment of membrane proteins and often afford novel mechanistic insights, we propose that nanodiscs ought to become the standard for structural studies on membrane proteins.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Current Developments in Native Nanometric Discoidal Membrane Bilayer Formed by Amphipathic Polymers

Unlike cytosolic proteins, membrane proteins (MPs) are embedded within the plasma membrane and the lipid bilayer of intracellular organelles. MPs serve in various cellular processes and account for over 65% of the current drug targets. The development of membrane mimetic systems such as bicelles, short synthetic polymers or amphipols, and membrane scaffold proteins (MSP)-based nanodiscs has facilitated the accommodation of synthetic lipids to stabilize MPs, yet the preparation of these membrane mimetics remains detergent-dependent. Bio-inspired synthetic polymers present an invaluable tool for excision and liberation of superstructures of MPs and their surrounding annular lipid bilayer in the nanometric discoidal assemblies. In this article, we discuss the significance of self-assembling process in design of biomimetic systems, review development of multiple series of amphipathic polymers and the significance of these polymeric "belts" in biomedical research in particular in unraveling the structures, dynamics and functions of several high-value membrane protein targets.

Aquaporins - Expression, purification and characterization

Aquaporin water channels facilitate the bi-directional flow of water and small, neutral solutes down an osmotic gradient in all kingdoms of life. Over the last two decades, the availability of high-quality protein has underpinned progress in the structural and functional characterization of these water channels. In particular, recombinant protein technology has guaranteed the supply of aquaporin samples that were of sufficient quality and quantity for further study. Here we review the features of successful expression, purification and characterization strategies that have underpinned these successes and that will drive further breakthroughs in the field. Overall, Escherichia coli is a suitable host for prokaryotic isoforms, while Pichia pastoris is the most commonly-used recombinant host for eukaryotic variants. Generally, a two-step purification procedure is suitable after solubilization in glucopyranosides and most structures are determined by X-ray following crystallization.

Sample Preparation and Data Collection for Electron Crystallographic Studies on Membrane Protein Structures and Lipid-Protein Interaction

Electron crystallography is a unique tool to study membrane protein structures and lipid-protein interactions in their native-like environments. Two-dimensional (2D) protein crystallization enables the lipids immobilized by the proteins, and the generated high-resolution density map allows us to model the atomic coordinates of the surrounding lipids to study lipid-protein interaction. This protocol describes the sample preparation for electron crystallographic studies, including back-injection method and carbon sandwich method. The protocols of data collection for electron crystallography, including electron imaging and diffraction, of the 2D membrane crystal will be followed.

Cholesterol Localization around the Metabotropic Glutamate Receptor 2

The metabotropic glutamate receptor (mGluR) 2 plays a key role in the central nervous system. mGluR2 has been shown to be regulated by its surrounding lipid environment, especially by cholesterol, by an unknown mechanism. Here, using a combination of biochemical approaches, photo-cross-linking experiments, and molecular dynamics simulations we show the interaction of cholesterol with at least two, but potentially five more, preferential sites on the mGluR2 transmembrane domain. Our simulations demonstrate that surface matching, rather than electrostatic interactions with specific amino acids, is the main factor defining cholesterol localization. Moreover, the cholesterol localization observed here is similar to the sterol-binding pattern previously described in silico for other members of the mGluR family. Biochemical assays suggest little influence of cholesterol on trafficking or dimerization of mGluR2. Nevertheless, simulations revealed a significant reduction of residue-residue contacts together with an alteration in the internal mechanical stress at the cytoplasmic side of the helical bundle when cholesterol was present in the membrane. These alterations may be related to destabilization of the basal state of mGluR2. Due to the high sequence conservation of the transmembrane domains of mGluRs, the molecular interaction of cholesterol and mGluR2 described here is also likely to be relevant for other members of the mGLuR family.

Molecular dynamics of the histamine H3 membrane receptor reveals different mechanisms of GPCR signal transduction

In this work, we studied the mechanisms of classical activation and inactivation of signal transduction by the histamine H3 receptor, a 7-helix transmembrane bundle G-Protein Coupled Receptor through long-time-scale atomistic molecular dynamics simulations of the receptor embedded in a hydrated double layer of dipalmitoyl phosphatidyl choline, a zwitterionic polysaturated ordered lipid. Three systems were prepared: the apo receptor, representing the constitutively active receptor; and two holo-receptors-the receptor coupled to the antagonist/inverse agonist ciproxifan, representing the inactive state of the receptor, and the receptor coupled to the endogenous agonist histamine and representing the active state of the receptor. An extensive analysis of the simulation showed that the three states of H3R present significant structural and dynamical differences as well as a complex behavior given that the measured properties interact in multiple and interdependent ways. In addition, the simulations described an unexpected escape of histamine from the orthosteric binding site, in agreement with the experimental modest affinities and rapid off-rates of agonists.

Degeneracy in molecular scale organization of biological membranes

The scale-rich spatiotemporal organization in biological membranes has its origin in the differential inter- and intra-molecular interactions among their constituents. In this work, we explore the molecular-origin behind that variety and possible degeneracy in lateral organization in membranes. For our study, we post-process microsecond long all-atom molecular dynamics trajectories for three systems that exhibit fluid phase coexistence: (i) PSM/POPC/Chol (0.47/0.32/0.21), (ii) PSM/DOPC/Chol (0.43/0.38/0.19) and (iii) DPPC/DOPC/Chol (0.37/0.36/0.27). To distinguish the liquid ordered and disordered regions at molecular scales, we calculate the degree of non-affineness of individual lipids in their neighbourhood and track their topological rearrangements. Disconnectivity graph analysis with respect to membrane organization shows that the DPPC/DOPC/Chol and PSM/DOPC/Chol systems exhibit funnel-like energy landscapes as opposed to a highly frustrated energy landscape for the more biomimetic PSM/POPC/Chol system. We use these measurements to develop a continuous lattice Hamiltonian and evolve that using Monte Carlo simulated annealing to explore the possibility of structural degeneracy in membrane organization. Our data show that model membranes with lipid constituents that are biomimetic (PSM/POPC/Chol) have the ability to access a large range of membrane sub-structure space (higher degeneracy) as compared to the other two systems, which form only one kind of substructure even with changing composition. Since the spatiotemporal organization in biological membranes dictates the "molecular encounters" and in turn larger scale biological processes such as molecular transport, trafficking and cellular signalling, we posit that this structural degeneracy could enable access to a larger repository to functionally important molecular organization in systems with physiologically relevant compositions.

Helical Inclusions in Phospholipid Membranes: Lipid Adaptation and Chiral Order

The lipid bilayer is a flexible matrix that is able to adapt in response to the perturbation induced by inclusions, such as peptides and proteins. Here we use molecular dynamics simulations with a coarse-grained model to investigate the effect of a helical inclusion on a lipid bilayer in the liquid disordered phase. We show that the helical inclusion induces a collective tilt of acyl chains, with a small, yet unambiguous difference between a right- and a left-handed inclusion. This behavior is rationalized using the elastic continuum theory: The magnitude of the chiral (twist) deformation of the bilayer is determined by the interaction at the lipid/inclusion interface, and the decay length is controlled by the elastic properties of the bilayer. The lipid reorganization can thus be identified as a generic mechanism that, together with specific interactions, contributes to chiral recognition in phospholipid bilayers. An enhanced response is expected in highly ordered environments, such as rafts in biomembranes, with a potential impact on membrane-mediated interactions between inclusions.

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.

Experimental and Simulation Studies of Aquaporin 0 Water Permeability and Regulation

We begin with the history of aquaporin zero (AQP0), the most prevalent membrane protein in the eye lens, from the early days when AQP0 was a protein of unknown function known as Major Intrinsic Protein 26. We progress through its joining the aquaporin family as a water channel in its own right and discuss how regulation of its water permeability by pH and calcium came to be discovered experimentally and linked to lens homeostasis and development. We review the development of molecular dynamics (MD) simulations of lipid bilayers and membrane proteins, including aquaporins, with an emphasis on simulation studies that have elucidated the mechanisms of water conduction, selectivity, and proton exclusion by aquaporins in general. We also review experimental and theoretical progress toward understanding why mammalian AQP0 has a lower water permeability than other aquaporins and the evolution of our present understanding of how its water permeability is regulated by pH and calcium. Finally, we discuss how MD simulations have elucidated the nature of lipid interactions with AQP0.

Hierarchical Optimization of High-Performance Biomimetic and Bioinspired Membranes

Biomimetic and bioinspired membranes have emerged as an innovative platform for water purification and aqueous separations. They are inspired by the exceptional water permeability (∼10⁹ water molecules per second per channel) and perfect selectivity of biological water channels, aquaporins. However, only few successes have been reported for channel-based membrane fabrication due to inherent challenges of realizing coherence between channel design at the angstrom level and development of scalable membranes that maintain these molecular properties at practice-relevant scales. In this article, we feature recent progress toward practical biomimetic membranes, with the review organized along a hierarchical structural perspective that biomimetic membranes commonly share. These structures range from unitary pore shapes and tubular hydrophobic channel geometries to self-assembled bilayer structures and finally to macroscale membranes covering a size range from the angstrom, to the micrometer scale, and finally to the centimeter and larger scales. To maximize the advantage of water channel implementation into membranes, each feature needs to be optimized in an appropriate manner that provides a path to successful scale-up to achieve high performance in practical biomimetic and bioinspired membranes.

Roughness of Transmembrane Helices Reduces Lipid Membrane Dynamics

The dynamics of cellular membranes is primarily determined by lipid species forming a bilayer. Proteins are considered mainly as effector molecules of diverse cellular processes. In addition to large assemblies of proteins, which were found to influence properties of fluid membranes, biological membranes are densely populated by small, highly mobile proteins. However, little is known about the effect of such proteins on the dynamics of membranes. Using synthetic peptides, we demonstrate that transmembrane helices interfere with the mobility of membrane components by trapping lipid acyl chains on their rough surfaces. The effect is more pronounced in the presence of cholesterol, which segregates from the rough surface of helical peptides. This may contribute to the formation or stabilization of membrane heterogeneities. Since roughness is a general property of helical transmembrane segments, our results suggest that, independent of their size or cytoskeleton linkage, integral membrane proteins affect local membrane dynamics and organization.

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.

The lipid bilayer membrane and its protein constituents

Robertson recounts the historical experiments that gave rise to our current understanding of cell membranes.

In 1918, the year the Journal of General Physiology was founded, there was little understanding of the structure of the cell membrane. It was evident that cells had invisible barriers separating the cytoplasm from the external solution. However, it would take decades before lipid bilayers were identified as the essential constituent of membranes. It would take even longer before it was accepted that there existed hydrophobic proteins that were embedded within the membrane and that these proteins were responsible for selective permeability in cells. With a combination of intuitive experiments and quantitative thinking, the last century of cell membrane research has led us to a molecular understanding of the structure of the membrane, as well as many of the proteins embedded within. Now, research is turning toward a physical understanding of the reactions of membrane proteins and lipids in this unique and incredibly complex solvent environment.

Plant and Mammal Aquaporins: Same but Different

Aquaporins (AQPs) constitute an ancient and diverse protein family present in all living organisms, indicating a common ancient ancestor. However, during evolution, these organisms appear and evolve differently, leading to different cell organizations and physiological processes. Amongst the eukaryotes, an important distinction between plants and animals is evident, the most conspicuous difference being that plants are sessile organisms facing ever-changing environmental conditions. In addition, plants are mostly autotrophic, being able to synthesize carbohydrates molecules from the carbon dioxide in the air during the process of photosynthesis, using sunlight as an energy source. It is therefore interesting to analyze how, in these different contexts specific to both kingdoms of life, AQP function and regulation evolved. This review aims at highlighting similarities and differences between plant and mammal AQPs. Emphasis is given to the comparison of isoform numbers, their substrate selectivity, the regulation of the subcellular localization, and the channel activity.

Plant and animal aquaporins crosstalk: what can be revealed from distinct perspectives

Aquaporins (AQPs) can be revisited from a distinct and complementary perspective: the outcome from analyzing them from both plant and animal studies. (1) The approach in the study. Diversity found in both kingdoms contrasts with the limited number of crystal structures determined within each group. While the structure of almost half of mammal AQPs was resolved, only a few were resolved in plants. Strikingly, the animal structures resolved are mainly derived from the AQP2-lineage, due to their important roles in water homeostasis regulation in humans. The difference could be attributed to the approach: relevance in animal research is emphasized on pathology and in consequence drug screening that can lead to potential inhibitors, enhancers and/or regulators. By contrast, studies on plants have been mainly focused on the physiological role that AQPs play in growth, development and stress tolerance. (2) The transport capacity. Besides the well-described AQPs with high water transport capacity, large amount of evidence confirms that certain plant AQPs can carry a large list of small solutes. So far, animal AQP list is more restricted. In both kingdoms, there is a great amount of evidence on gas transport, although there is still an unsolved controversy around gas translocation as well as the role of the central pore of the tetramer. (3) More roles than expected. We found it remarkable that the view of AQPs as specific channels has evolved first toward simple transporters to molecules that can experience conformational changes triggered by biochemical and/or mechanical signals, turning them also into signaling components and/or behave as osmosensor molecules.

Stoichiometry determination of macromolecular membrane protein complexes

Gaining knowledge of the structural makeup of protein complexes is critical to advance our understanding of their formation and functions. This task is particularly challenging for transmembrane protein complexes, and grows ever more imposing with increasing size of these large macromolecular structures. The last 10 years have seen a steep increase in solved high-resolution membrane protein structures due to both new and improved methods in the field, but still most structures of large transmembrane complexes remain elusive. An important first step towards the structure elucidation of these difficult complexes is the determination of their stoichiometry, which we discuss in this review. Knowing the stoichiometry of complex components not only answers unresolved structural questions and is relevant for understanding the molecular mechanisms of macromolecular machines but also supports further attempts to obtain high-resolution structures by providing constraints for structure calculations.

Temperature dependence of protein-water interactions in a gated yeast aquaporin

Regulation of aquaporins is a key process of living organisms to counteract sudden osmotic changes. Aqy1, which is a water transporting aquaporin of the yeast Pichia pastoris, is suggested to be gated by chemo-mechanical stimuli as a protective regulatory-response against rapid freezing. Here, we tested the influence of temperature by determining the X-ray structure of Aqy1 at room temperature (RT) at 1.3 Å resolution, and by exploring the structural dynamics of Aqy1 during freezing through molecular dynamics simulations. At ambient temperature and in a lipid bilayer, Aqy1 adopts a closed conformation that is globally better described by the RT than by the low-temperature (LT) crystal structure. Locally, for the blocking-residue Tyr31 and the water molecules inside the pore, both LT and RT data sets are consistent with the positions observed in the simulations at room-temperature. Moreover, as the temperature was lowered, Tyr31 adopted a conformation that more effectively blocked the channel, and its motion was accompanied by a temperature-driven rearrangement of the water molecules inside the channel. We therefore speculate that temperature drives Aqy1 from a loosely- to a tightly-blocked state. This analysis provides high-resolution structural evidence of the influence of temperature on membrane-transport channels.

Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects

Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.

Sub-nanometre mapping of the aquaporin-water interface using multifrequency atomic force microscopy

Aquaporins are integral membrane proteins that regulate the transport of water and small molecules in and out of the cell. In eye lens tissue, circulation of water, ions and metabolites is ensured by a microcirculation system in which aquaporin-0 (AQP0) plays a central role. AQP0 allows water to flow beyond the diffusion limit through lens membranes. AQP0 naturally arranges in a square lattice. The malfunction of AQP0 is related to numerous diseases such as cataracts. Despite considerable research into its structure, function and dynamics, the interface between the protein and the surrounding liquid and the effect of the lattice arrangement on the behaviour of water at the interface with the membrane are still not fully understood. Here we use a multifrequency atomic force microscopy (AFM) approach to map both the liquid at the interface with AQP0 and the protein itself with sub-nanometer resolution. Imaging using the fundamental eigenmode of the AFM cantilever probes mainly the interfacial water at the surface of the membrane. The results highlight a well-defined region that surrounds AQP0 tetramers and where water exhibits a higher affinity for the protein. Imaging in the second eigenmode is dominated by the mechanical response of the protein and provides sub-molecular details of the protein surface and the sub-surface structure. The relationship between modes and harmonics is also examined.

A sliding selectivity scale for lipid binding to membrane proteins

Biological membranes form barriers that are essential for cellular integrity and compartmentalisation. Proteins in the membrane have co-evolved with their hydrophobic lipid environment, which serves as a solvent for proteins with very diverse requirements. As a result, their interactions range from non-selective to highly discriminating. Mass spectrometry enables us to monitor how lipids interact with membrane proteins and assess their effects on structure and dynamics. Recent studies illustrate the ability to differentiate specific lipid binding, preferential interactions with lipid subsets, and nonselective annular contacts. Here, we consider the biological implications of different lipid-binding scenarios and propose that binding occurs on a sliding selectivity scale, in line with the view of biological membranes as facilitators of dynamic protein and lipid organization.

Control of B Cell Responsiveness by Isotype and Structural Elements of the Antigen Receptor

Expression of a functional B cell antigen receptor (BCR) plays a central role in regulating B cell development, maturation, and effector functions. Although IgM is solely expressed in immature B cell stages, the presence of both IgM- and IgD-BCR isotypes on mature naïve B cells raises the question of whether IgD has a unique role in B cell activation and function. While earlier studies suggested a broad functional redundancy between IgM and IgD, recent data point to an important immune regulatory role of IgD. Herein, we review these findings and discuss how the structural flexibility, mode of antigen binding, and co-receptor interactions, enable the IgD-BCR to act as a 'rheostat', regulating the activation and function of mature naïve B cells.

Single-particle electron microscopy in the study of membrane protein structure

Single-particle electron microscopy (EM) provides the great advantage that protein structure can be studied without the need to grow crystals. However, due to technical limitations, this approach played only a minor role in the study of membrane protein structure. This situation has recently changed dramatically with the introduction of direct electron detection device cameras, which allow images of unprecedented quality to be recorded, also making software algorithms, such as three-dimensional classification and structure refinement, much more powerful. The enhanced potential of single-particle EM was impressively demonstrated by delivering the first long-sought atomic model of a member of the biomedically important transient receptor potential channel family. Structures of several more membrane proteins followed in short order. This review recounts the history of single-particle EM in the study of membrane proteins, describes the technical advances that now allow this approach to generate atomic models of membrane proteins and provides a brief overview of some of the membrane protein structures that have been studied by single-particle EM to date.

Evaluation of super-resolution performance of the K2 electron-counting camera using 2D crystals of aquaporin-0

The K2 Summit camera was initially the only commercially available direct electron detection camera that was optimized for high-speed counting of primary electrons and was also the only one that implemented centroiding so that the resolution of the camera can be extended beyond the Nyquist limit set by the physical pixel size. In this study, we used well-characterized two-dimensional crystals of the membrane protein aquaporin-0 to characterize the performance of the camera below and beyond the physical Nyquist limit and to measure the influence of electron dose rate on image amplitudes and phases.

Aquaporin-Based Biomimetic Polymeric Membranes: Approaches and Challenges

In recent years, aquaporin biomimetic membranes (ABMs) for water separation have gained considerable interest. Although the first ABMs are commercially available, there are still many challenges associated with further ABM development. Here, we discuss the interplay of the main components of ABMs: aquaporin proteins (AQPs), block copolymers for AQP reconstitution, and polymer-based supporting structures. First, we briefly cover challenges and review recent developments in understanding the interplay between AQP and block copolymers. Second, we review some experimental characterization methods for investigating AQP incorporation including freeze-fracture transmission electron microscopy, fluorescence correlation spectroscopy, stopped-flow light scattering, and small-angle X-ray scattering. Third, we focus on recent efforts in embedding reconstituted AQPs in membrane designs that are based on conventional thin film interfacial polymerization techniques. Finally, we describe some new developments in interfacial polymerization using polyhedral oligomeric silsesquioxane cages for increasing the physical and chemical durability of thin film composite membranes.

Highly permeable artificial water channels that can self-assemble into two-dimensional arrays

Bioinspired artificial water channels aim to combine the high permeability and selectivity of biological aquaporin (AQP) water channels with chemical stability. Here, we carefully characterized a class of artificial water channels, peptide-appended pillar[5]arenes (PAPs). The average single-channel osmotic water permeability for PAPs is 1.0(± 0.3) × 10(-14) cm(3)/s or 3.5(± 1.0) × 10(8) water molecules per s, which is in the range of AQPs (3.4 ∼ 40.3 × 10(8) water molecules per s) and their current synthetic analogs, carbon nanotubes (CNTs, 9.0 × 10(8) water molecules per s). This permeability is an order of magnitude higher than first-generation artificial water channels (20 to ∼ 10(7) water molecules per s). Furthermore, within lipid bilayers, PAP channels can self-assemble into 2D arrays. Relevant to permeable membrane design, the pore density of PAP channel arrays (∼ 2.6 × 10(5) pores per μm(2)) is two orders of magnitude higher than that of CNT membranes (0.1 ∼ 2.5 × 10(3) pores per μm(2)). PAP channels thus combine the advantages of biological channels and CNTs and improve upon them through their relatively simple synthesis, chemical stability, and propensity to form arrays.

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

CholMine: Determinants and Prediction of Cholesterol and Cholate Binding Across Nonhomologous Protein Structures

Identifying physiological ligands is necessary for annotating new protein structures, yet this presents a significant challenge to biologists and pharmaceutical chemists. Here we develop a predictor of cholesterol and cholate binding that works across diverse protein families, extending beyond sequence motif-based prediction. This approach combines SimSite3D site comparison with the detection of conserved interactions in cholesterol/cholate bound crystal structures to define three-dimensional interaction motifs. The resulting predictor identifies cholesterol sites with an ∼82% unbiased true positive rate in both membrane and soluble proteins, with a very low false positive rate relative to other predictors. The CholMine Web server can analyze users' structures, detect those likely to bind cholesterol/cholate, and predict the binding mode and key interactions. By deciphering the determinants of binding for these important steroids, CholMine may also aid in the design of selective inhibitors and detergents for targets such as G protein coupled receptors and bile acid receptors.

Molecular dynamics simulations of the bacterial UraA H+-uracil symporter in lipid bilayers reveal a closed state and a selective interaction with cardiolipin

The Escherichia coli UraA H⁺-uracil symporter is a member of the nucleobase/ascorbate transporter (NAT) family of proteins, and is responsible for the proton-driven uptake of uracil. Multiscale molecular dynamics simulations of the UraA symporter in phospholipid bilayers consisting of: 1) 1-palmitoyl 2-oleoyl-phosphatidylcholine (POPC); 2) 1-palmitoyl 2-oleoyl-phosphatidylethanolamine (POPE); and 3) a mixture of 75% POPE, 20% 1-palmitoyl 2-oleoyl-phosphatidylglycerol (POPG); and 5% 1-palmitoyl 2-oleoyl-diphosphatidylglycerol/cardiolipin (CL) to mimic the lipid composition of the bacterial inner membrane, were performed using the MARTINI coarse-grained force field to self-assemble lipids around the crystal structure of this membrane transport protein, followed by atomistic simulations. The overall fold of the protein in lipid bilayers remained similar to the crystal structure in detergent on the timescale of our simulations. Simulations were performed in the absence of uracil, and resulted in a closed state of the transporter, due to relative movement of the gate and core domains. Anionic lipids, including POPG and especially CL, were found to associate with UraA, involving interactions between specific basic residues in loop regions and phosphate oxygens of the CL head group. In particular, three CL binding sites were identified on UraA: two in the inner leaflet and a single site in the outer leaflet. Mutation of basic residues in the binding sites resulted in the loss of CL binding in the simulations. CL may play a role as a “proton trap” that channels protons to and from this transporter within CL-enriched areas of the inner bacterial membrane.

Symporters are proteins that are responsible for the co-transport of ions and small molecule solutes across cell membranes. UraA is an example of a symporter, and is responsible for the proton-driven uptake of uracil in bacteria like E. coli. Despite its importance as a member of a large family of nucleobase/ascorbate transporters (NAT) and the existence of structural and functional data, the mechanism by which UraA transports uracil across the bacterial membrane, and in particular the role of its diverse and complex lipid environment in the transport mechanism, remains elusive. In this study, we have used a multiscale computational methodology to examine the dynamics of UraA and to elucidate its interactions with lipids that resemble its native environment in the bacterial inner membrane. Our results demonstrate that negatively-charged lipids in the membrane (phosphatidylglycerol and cardiolipin) associate preferentially with UraA and may play a role in its function. Additionally, our simulations resulted in a closed state of UraA, a likely intermediate in the transport mechanism that may not be readily accessible by experimental methods.

Effect of lipid head groups on double-layered two-dimensional crystals formed by aquaporin-0

Aquaporin-0 (AQP0) is a lens-specific water channel that also forms membrane junctions. Reconstitution of AQP0 with dimyristoyl phosphatidylcholine (DMPC) and E. coli polar lipids (EPL) yielded well-ordered, double-layered two-dimensional (2D) crystals that allowed electron crystallographic structure determination of the AQP0-mediated membrane junction. The interacting tetramers in the two crystalline layers are exactly in register, resulting in crystals with p422 symmetry. The high-resolution density maps also allowed modeling of the annular lipids surrounding the tetramers. Comparison of the DMPC and EPL bilayers suggested that the lipid head groups do not play an important role in the interaction of annular lipids with AQP0. We now reconstituted AQP0 with the anionic lipid dimyristoyl phosphatidylglycerol (DMPG), which yielded a mixture of 2D crystals with different symmetries. The different crystal symmetries result from shifts between the two crystalline layers, suggesting that the negatively charged PG head group destabilizes the interaction between the extracellular AQP0 surfaces. Reconstitution of AQP0 with dimyristoyl phosphatidylserine (DMPS), another anionic lipid, yielded crystals that had the usual p422 symmetry, but the crystals showed a pH-dependent tendency to stack through their cytoplasmic surfaces. Finally, AQP0 failed to reconstitute into membranes that were composed of more than 40% dimyristoyl phosphatidic acid (DMPA). Hence, although DMPG, DMPS, and DMPA are all negatively charged lipids, they have very different effects on AQP0 2D crystals, illustrating the importance of the specific lipid head group chemistry beyond its mere charge.

The nanoscale organization of the B lymphocyte membrane

The fluid mosaic model of Singer and Nicolson correctly predicted that the plasma membrane (PM) forms a lipid bi-layer containing many integral trans-membrane proteins. This model also suggested that most of these proteins were randomly dispersed and freely diffusing moieties. Initially, this view of a dynamic and rather unorganized membrane was supported by early observations of the cell surfaces using the light microscope. However, recent studies on the PM below the diffraction limit of visible light (~250nm) revealed that, at nanoscale dimensions, membranes are highly organized and compartmentalized structures. Lymphocytes are particularly useful to study this nanoscale membrane organization because they grow as single cells and are not permanently engaged in cell:cell contacts within a tissue that can influence membrane organization. In this review, we describe the methods that can be used to better study the protein:protein interaction and nanoscale organization of lymphocyte membrane proteins, with a focus on the B cell antigen receptor (BCR). Furthermore, we discuss the factors that may generate and maintain these membrane structures.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Internal lipid architecture of the hetero-oligomeric cytochrome b6f complex

The role of lipids in the assembly, structure, and function of hetero-oligomeric membrane protein complexes is poorly understood. The dimeric photosynthetic cytochrome b6f complex, a 16-mer of eight distinct subunits and 26 transmembrane helices, catalyzes transmembrane proton-coupled electron transfer for energy storage. Using a 2.5 Å crystal structure of the dimeric complex, we identified 23 distinct lipid-binding sites per monomer. Annular lipids are proposed to provide a connection for super-complex formation with the photosystem-I reaction center and the LHCII kinase enzyme for transmembrane signaling. Internal lipids mediate crosslinking to stabilize the domain-swapped iron-sulfur protein subunit, dielectric heterogeneity within intermonomer and intramonomer electron transfer pathways, and dimer stabilization through lipid-mediated intermonomer interactions. This study provides a complete structure analysis of lipid-mediated functions in a multi-subunit membrane protein complex and reveals lipid sites at positions essential for assembly and function.

Membrane proteins bind lipids selectively to modulate their structure and function

Previous studies have established that the folding, structure and function of membrane proteins are influenced by their lipid environments and that lipids can bind to specific sites, for example, in potassium channels. Fundamental questions remain however regarding the extent of membrane protein selectivity towards lipids. Here we report a mass spectrometry approach designed to determine the selectivity of lipid binding to membrane protein complexes. We investigate the mechanosensitive channel of large conductance (MscL) from Mycobacterium tuberculosis and aquaporin Z (AqpZ) and the ammonia channel (AmtB) from Escherichia coli, using ion mobility mass spectrometry (IM-MS), which reports gas-phase collision cross-sections. We demonstrate that folded conformations of membrane protein complexes can exist in the gas phase. By resolving lipid-bound states, we then rank bound lipids on the basis of their ability to resist gas phase unfolding and thereby stabilize membrane protein structure. Lipids bind non-selectively and with high avidity to MscL, all imparting comparable stability; however, the highest-ranking lipid is phosphatidylinositol phosphate, in line with its proposed functional role in mechanosensation. AqpZ is also stabilized by many lipids, with cardiolipin imparting the most significant resistance to unfolding. Subsequently, through functional assays we show that cardiolipin modulates AqpZ function. Similar experiments identify AmtB as being highly selective for phosphatidylglycerol, prompting us to obtain an X-ray structure in this lipid membrane-like environment. The 2.3 Å resolution structure, when compared with others obtained without lipid bound, reveals distinct conformational changes that re-position AmtB residues to interact with the lipid bilayer. Our results demonstrate that resistance to unfolding correlates with specific lipid-binding events, enabling a distinction to be made between lipids that merely bind from those that modulate membrane protein structure and/or function. We anticipate that these findings will be important not only for defining the selectivity of membrane proteins towards lipids, but also for understanding the role of lipids in modulating protein function or drug binding.

Growth of large and highly ordered 2D crystals of a K⁺ channel, structural role of lipidic environment

2D crystallography has proven to be an excellent technique to determine the 3D structure of membrane proteins. Compared to 3D crystallography, it has the advantage of visualizing the protein in an environment closer to the native one. However, producing good 2D crystals is still a challenge and little statistical knowledge can be gained from literature. Here, we present a thorough screening of 2D crystallization conditions for a prokaryotic inwardly rectifying potassium channel (>130 different conditions). Key parameters leading to very large and well-organized 2D crystals are discussed. In addition, the problem of formation of multilayers during the growth of 2D crystals is also addressed. An intermediate resolution projection map of KirBac3.1 at 6 Å is presented, which sheds (to our knowledge) new light on the structure of this channel in a lipid environment.