Molecular driving forces defining lipid positions around aquaporin-0

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.

Ionic Liquid Gating Induces Anomalous Permeation through Membrane Channel Proteins

Membrane channel proteins (MCPs) play key roles in matter transport through cell membranes and act as major targets for vaccines and drugs. For emerging ionic liquid (IL) drugs, a rational understanding of how ILs affect the structure and transport function of MCP is crucial to their design. In this work, GPU-accelerated microsecond-long molecular dynamics simulations were employed to investigate the modulating mechanism of ILs on MCP. Interestingly, ILs prefer to insert into the lipid bilayer and channel of aquaporin-2 (AQP2) but adsorb on the entrance of voltage-gated sodium channels (Nav). Molecular trajectory and free energy analysis reflect that ILs have a minimal impact on the structure of MCPs but significantly influence MCP functions. It demonstrates that ILs can decrease the overall energy barrier for water through AQP2 by 1.88 kcal/mol, whereas that for Na+ through Nav is increased by 1.70 kcal/mol. Consequently, the permeation rates of water and Na+ can be enhanced and reduced by at least 1 order of magnitude, respectively. Furthermore, an abnormal IL gating mechanism was proposed by combining the hydrophobic nature of MCP and confined water/ion coordination effects. More importantly, we performed experiments to confirm the influence of ILs on AQP2 in human cells and found that treatment with ILs significantly accelerated the changes in cell volume in response to altered external osmotic pressure. Overall, these quantitative results will not only deepen the understanding of IL-cell interactions but may also shed light on the rational design of drugs and disease diagnosis.

Non-affine deformation analysis and 3D packing defects: A new way to probe membrane heterogeneity in molecular simulations

Here, we discuss a new framework developed over the last 5 years in our group to probe nanoscale membrane heterogeneity. The framework is based on the idea of characterizing lateral heterogeneity through non-affine deformation (NAD) measurements, transverse heterogeneity through three dimensional (3D) lipid packing defects, and using these approaches to formalize the seemingly trivial correlation between lateral organization and lipid packing in biological membranes. We find that measurements from NAD analysis, a prescription which is borrowed from Physics of glasses and granular material, can faithfully distinguish between liquid-ordered and disordered phases in membranes at molecular length scales and, can also be used to identify phase boundaries with high precision. Concomitantly, 3D-packing defects can not only distinguish between the two co-existing fluid phases based on their molecular scale packing (or membrane free volume), but also provide a route to connect the membrane domains to their functionality, such as exploring the molecular origins of inter-leaflet domain registration and peptide partitioning. The correlation between lateral membrane order and transverse packing presents novel molecular design-level features that can explain functions such as protein/peptide partitioning and small-molecule permeation dynamics in complex and heterogeneous membranes with high-fidelity. The framework allows us to explore the nature of lateral organization and molecular packing as a manifestation of intricate molecular interactions among a chemically rich variety of lipids and other molecules in a membrane with complex membrane composition and asymmetry across leaflets.

Protein Interface Regulating the Inserting Process of Imidazole Ionic Liquids into the Cell Membrane

Ionic liquids (ILs) have shown promising potential in membrane protein extraction; however, the underlying mechanism remains unclear. Herein, we employed GPU-accelerated molecular dynamics (MD) simulations to investigate the dynamic insertion process of ILs into cell membranes containing membrane proteins. Our findings reveal that ILs spontaneously insert into the membrane, and the presence of membrane proteins significantly decelerates the rate of IL insertion into the membrane. Specifically, the relationship between the insertion rate and inserting free energy exhibits non-monotonic changes, which can be attributed to interfacial effects. The protein-water interface acts as trap for free ions and ionic clusters, while free ions preferentially insert into the membrane from the protein-lipid interface, which limits the insertion rate due to its narrowness. Thus, the insertion rate is governed by a combination of the free energy and interfacial effects. These findings provide valuable insights into the interfacial effects of protein-lipid bilayers and have implications for various biochemical-related applications.

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.

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.

Translocating Peptides of Biomedical Interest Obtained from the Spike (S) Glycoprotein of the SARS-CoV-2

At the beginning of 2020, the pandemic caused by the SARS-CoV-2 virus led to the fast sequencing of its genome to facilitate molecular engineering strategies to control the pathogen’s spread. The spike (S) glycoprotein has been identified as the leading therapeutic agent due to its role in localizing the ACE2 receptor in the host’s pulmonary cell membrane, binding, and eventually infecting the cells. Due to the difficulty of delivering bioactive molecules to the intracellular space, we hypothesized that the S protein could serve as a source of membrane translocating peptides. AHB-1, AHB-2, and AHB-3 peptides were identified and analyzed on a membrane model of DPPC (dipalmitoylphosphatidylcholine) using molecular dynamics (MD) simulations. An umbrella sampling approach was used to quantify the energy barrier necessary to cross the boundary (13.2 to 34.9 kcal/mol), and a flat-bottom pulling helped to gain a deeper understanding of the membrane’s permeation dynamics. Our studies revealed that the novel peptide AHB-1 exhibited comparable penetration potential of already known potent cell-penetrating peptides (CPPs) such as TP2, Buforin II, and Frenatin 2.3s. Results were confirmed by in vitro analysis of the peptides conjugated to chitosan nanoparticles, demonstrating its ability to reach the cytosol and escape endosomes, while maintaining high biocompatibility levels according to standardized assays.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Aquaporin ion conductance properties defined by membrane environment, protein structure, and cell physiology

Aquaporins (AQPs) are multifunctional transmembrane channel proteins permeable to water and an expanding array of solutes. AQP-mediated ion channel activity was first observed when purified AQP0 from bovine lens was incorporated into lipid bilayers. Electrophysiological properties of ion-conducting AQPs since discovered in plants, invertebrates, and mammals have been assessed using native, reconstituted, and heterologously expressed channels. Accumulating evidence is defining amino acid residues that govern differential solute permeability through intrasubunit and central pores of AQP tetramers. Rings of charged and hydrophobic residues around pores influence AQP selectivity, and are candidates for further work to define motifs that distinguish ion conduction capability, versus strict water and glycerol permeability. Similarities between AQP ion channels thus far include large single channel conductances and long open times, but differences in ionic selectivity, permeability to divalent cations, and mechanisms of gating (e.g., by voltage, pH, and cyclic nucleotides) are unique to subtypes. Effects of lipid environments in modulating parameters such as single channel amplitude could explain in part the variations in AQP ion channel properties observed across preparations. Physiological roles of the ion-conducting AQP classes span diverse processes including regulation of cell motility, organellar pH, neural development, signaling, and nutrient acquisition. Advances in computational methods can generate testable predictions of AQP structure-function relationships, which combined with innovative high-throughput assays could revolutionize the field in defining essential properties of ion-conducting AQPs, discovering new AQP ion channels, and understanding the effects of AQP interactions with proteins, signaling cascades, and membrane lipids.

The Transporter-Mediated Cellular Uptake and Efflux of Pharmaceutical Drugs and Biotechnology Products: How and Why Phospholipid Bilayer Transport Is Negligible in Real Biomembranes

Over the years, my colleagues and I have come to realise that the likelihood of pharmaceutical drugs being able to diffuse through whatever unhindered phospholipid bilayer may exist in intact biological membranes in vivo is vanishingly low. This is because (i) most real biomembranes are mostly protein, not lipid, (ii) unlike purely lipid bilayers that can form transient aqueous channels, the high concentrations of proteins serve to stop such activity, (iii) natural evolution long ago selected against transport methods that just let any undesirable products enter a cell, (iv) transporters have now been identified for all kinds of molecules (even water) that were once thought not to require them, (v) many experiments show a massive variation in the uptake of drugs between different cells, tissues, and organisms, that cannot be explained if lipid bilayer transport is significant or if efflux were the only differentiator, and (vi) many experiments that manipulate the expression level of individual transporters as an independent variable demonstrate their role in drug and nutrient uptake (including in cytotoxicity or adverse drug reactions). This makes such transporters valuable both as a means of targeting drugs (not least anti-infectives) to selected cells or tissues and also as drug targets. The same considerations apply to the exploitation of substrate uptake and product efflux transporters in biotechnology. We are also beginning to recognise that transporters are more promiscuous, and antiporter activity is much more widespread, than had been realised, and that such processes are adaptive (i.e., were selected by natural evolution). The purpose of the present review is to summarise the above, and to rehearse and update readers on recent developments. These developments lead us to retain and indeed to strengthen our contention that for transmembrane pharmaceutical drug transport "phospholipid bilayer transport is negligible".

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.

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.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

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.

Understanding the Role of Lipids in Signaling Through Atomistic and Multiscale Simulations of Cell Membranes

Cell signaling controls essentially all cellular processes. While it is often assumed that proteins are the key architects coordinating cell signaling, recent studies have shown more and more clearly that lipids are also involved in signaling processes in a number of ways. Lipids do, for instance, act as messengers, modulate membrane receptor conformation and dynamics, and control membrane receptor partitioning. Further, through structural modifications such as oxidation, the functions of lipids as part of signaling processes can be modified. In this context, in this article we discuss the understanding recently revealed by atomistic and coarse-grained computer simulations of nanoscale processes and underlying physicochemical principles related to lipids’ functions in cellular signaling.

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.

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.

Aquaporins: More Than Functional Monomers in a Tetrameric Arrangement

Aquaporins (AQPs) function as tetrameric structures in which each monomer has its own permeable pathway. The combination of structural biology, molecular dynamics simulations, and experimental approaches has contributed to improve our knowledge of how protein conformational changes can challenge its transport capacity, rapidly altering the membrane permeability. This review is focused on evidence that highlights the functional relationship between the monomers and the tetramer. In this sense, we address AQP permeation capacity as well as regulatory mechanisms that affect the monomer, the tetramer, or tetramers combined in complex structures. We therefore explore: (i) water permeation and recent evidence on ion permeation, including the permeation pathway controversy-each monomer versus the central pore of the tetramer-and (ii) regulatory mechanisms that cannot be attributed to independent monomers. In particular, we discuss channel gating and AQPs that sense membrane tension. For the latter we propose a possible mechanism that includes the monomer (slight changes of pore shape, the number of possible H-bonds between water molecules and pore-lining residues) and the tetramer (interactions among monomers and a positive cooperative effect).

Prediction of Thylakoid Lipid Binding Sites on Photosystem II

The thylakoid membrane has a unique lipid composition, consisting mostly of galactolipids. These thylakoid lipids have important roles in photosynthesis. Here, we investigate to what extent these lipids bind specifically to the Photosystem II complex. To this end, we performed coarse-grain MD simulations of the Photosystem II complex embedded in a thylakoid membrane with realistic composition. Based on >85 μs simulation time, we find that monogalactosyldiacylglycerol and sulfoquinovosyldiacylglycerol lipids are enriched in the annular shell around the protein, and form distinct binding sites. From the analysis of residue contacts, we conclude that electrostatic interactions play an important role in stabilizing these binding sites. Furthermore, we find that chlorophyll a has a prevalent role in the coordination of the lipids. In addition, we observe lipids to diffuse in and out of the plastoquinone exchange cavities, allowing exchange of cocrystallized lipids with the bulk membrane and suggesting a more open nature of the plastoquinone exchange cavity. Together, our data provide a wealth of information on protein-lipid interactions for a key protein in photosynthesis.

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.

Lipid perturbation by membrane proteins and the lipophobic effect

Understanding how membrane proteins interact with their environment is fundamental to the understanding of their structure, function and interactions. We have performed coarse-grained molecular dynamics simulations on a series of membrane proteins in a membrane environment to examine the perturbations of the lipids by the presence of protein. We analyze these perturbations in terms of elastic membrane deformations and local lipid protein interactions. However these two factors are insufficient to describe the variety of effects that we observe and the changes caused by membranes proteins to the structure and dynamics of their lipid environment. Other factors that change the conformation available to lipid molecules are evident and are able to modify lipid structure far from the protein surface, and thus mediate long-range interactions between membrane proteins. We suggest that these multiple modifications to lipid behavior are responsible, at the molecular level, for the lipophobic effect we have proposed to account for our observations of membrane protein organization.

Molecular dynamics simulations of membrane proteins and their interactions: from nanoscale to mesoscale

Molecular dynamics simulations provide a computational tool to probe membrane proteins and systems at length scales ranging from nanometers to close to a micrometer, and on microsecond timescales. All atom and coarse-grained simulations may be used to explore in detail the interactions of membrane proteins and specific lipids, yielding predictions of lipid binding sites in good agreement with available structural data. Building on the success of protein-lipid interaction simulations, larger scale simulations reveal crowding and clustering of proteins, resulting in slow and anomalous diffusional dynamics, within realistic models of cell membranes. Current methods allow near atomic resolution simulations of small membrane organelles, and of enveloped viruses to be performed, revealing key aspects of their structure and functionally important dynamics.

Insights into the molecular basis for substrate binding and specificity of the wild-type L-arginine/agmatine antiporter AdiC

Pathogenic enterobacteria need to survive the extreme acidity of the stomach to successfully colonize the human gut. Enteric bacteria circumvent the gastric acid barrier by activating extreme acid-resistance responses, such as the arginine-dependent acid resistance system. In this response, l-arginine is decarboxylated to agmatine, thereby consuming one proton from the cytoplasm. In Escherichia coli, the l-arginine/agmatine antiporter AdiC facilitates the export of agmatine in exchange of l-arginine, thus providing substrates for further removal of protons from the cytoplasm and balancing the intracellular pH. We have solved the crystal structures of wild-type AdiC in the presence and absence of the substrate agmatine at 2.6-Å and 2.2-Å resolution, respectively. The high-resolution structures made possible the identification of crucial water molecules in the substrate-binding sites, unveiling their functional roles for agmatine release and structure stabilization, which was further corroborated by molecular dynamics simulations. Structural analysis combined with site-directed mutagenesis and the scintillation proximity radioligand binding assay improved our understanding of substrate binding and specificity of the wild-type l-arginine/agmatine antiporter AdiC. Finally, we present a potential mechanism for conformational changes of the AdiC transport cycle involved in the release of agmatine into the periplasmic space of E. coli.

Mutual Interactions between Aquaporins and Membrane Components

In recent years, a number of studies have been focused on the structural evaluation of protein complexes in order to get mechanistic insights into how proteins communicate at the molecular level within the cell. Specific sites of protein-aquaporin interaction have been evaluated and new forms of regulation of aquaporins described, based on these associations. Heterotetramerizations of aquaporin isoforms are considered as novel regulatory mechanisms for plasma membrane (PIPs) and tonoplast (TIPs) proteins, influencing their intrinsic permeability and trafficking dynamics in the adaptive response to changing environmental conditions. However, protein-protein interaction is an extensive theme that is difficult to tackle and new methodologies are being used to study the physical interactions involved. Bimolecular fluorescence complementation and the identification of cross-linked peptides based on tandem mass spectra, that are complementary to other methodologies such as heterologous expression, co-precipitation assays or confocal fluorescence microscopy, are discussed in this review. The chemical composition and the physical characteristics of the lipid bilayer also influence many aspects of membrane aquaporins, including their functionality. The molecular driving forces stabilizing the positions of the lipids around aquaporins could define their activity, thereby altering the conformational properties. Therefore, an integrative approach to the relevance of the membrane-aquaporin interaction to different processes related to plant cell physiology is provided. Finally, it is described how the interactions between aquaporins and copolymer matrixes or biological compounds offer an opportunity for the functional incorporation of aquaporins into new biotechnological advances.

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.

Human Aquaporin-4 and Molecular Modeling: Historical Perspective and View to the Future

Among the different aquaporins (AQPs), human aquaporin-4 (hAQP4) has attracted the greatest interest in recent years as a new promising therapeutic target. Such a membrane protein is, in fact, involved in a multiple sclerosis-like immunopathology called Neuromyelitis Optica (NMO) and in several disorders resulting from imbalanced water homeostasis such as deafness and cerebral edema. The gap of knowledge in its functioning and dynamics at the atomistic level of detail has hindered the development of rational strategies for designing hAQP4 modulators. The application, lately, of molecular modeling has proved able to fill this gap providing a breeding ground to rationally address compounds targeting hAQP4. In this review, we give an overview of the important advances obtained in this field through the application of Molecular Dynamics (MD) and other complementary modeling techniques. The case studies presented herein are discussed with the aim of providing important clues for computational chemists and biophysicists interested in this field and looking for new challenges.

Lipid interaction sites on channels, transporters and receptors: Recent insights from molecular dynamics simulations

Lipid molecules are able to selectively interact with specific sites on integral membrane proteins, and modulate their structure and function. Identification and characterization of these sites are of importance for our understanding of the molecular basis of membrane protein function and stability, and may facilitate the design of lipid-like drug molecules. Molecular dynamics simulations provide a powerful tool for the identification of these sites, complementing advances in membrane protein structural biology and biophysics. We describe recent notable biomolecular simulation studies which have identified lipid interaction sites on a range of different membrane proteins. The sites identified in these simulation studies agree well with those identified by complementary experimental techniques. This demonstrates the power of the molecular dynamics approach in the prediction and characterization of lipid interaction sites on integral membrane proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Structural Properties of Protein-Detergent Complexes from SAXS and MD Simulations

In experimental studies of solubilized membrane proteins, the detergent corona influences the protein behavior and the resulting measurement. Thus, combinations of experimental techniques with atomistic modeling have been used to resolve corona structural parameters and distributions. Here, we used small-angle X-ray scattering (SAXS) data and molecular dynamics simulations to study a model protein-detergent complex (PDC) consisting of aquaporin-0 and dodecyl-β-maltoside molecules (βDDM). The corona morphology of single snapshots was found to be rough, but it is smooth and compacted in 100-ns-scale ensemble averages. Individual snapshots therefore were unable to accurately represent the ensemble information as captured by experimental SAXS. Mimicking of annular lipids by detergent was also observed. SAXS prediction using different published methods was used to identify optimal βDDM numbers. Explicit-solvent methods predicted best agreement using 290-βDDM PDCs, but implicit-solvent methods gave unclear predictions due to overcompensation by free solvation-layer density parameters. Thus, ensemble-based approaches and physically motivated constraints will help to extract structural information from SAXS data.

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.

Molecular dynamics simulations of the transport of reactive oxygen species by mammalian and plant aquaporins

BACKGROUND

Aquaporins are responsible for water transport across lipid membranes. They are also able to transport reactive oxygen species, playing an important role in redox signaling. Certain plant aquaporins have even the ability to be regulated by oxidative stress. However, the underlying mechanisms are still not fully understood.

METHODS

Here, molecular dynamics simulations were employed to determine the activation free energies related to the transport of reactive oxygen species through both mammalian and plant aquaporin models.

RESULTS AND CONCLUSIONS

Both aquaporins may transport hydrogen peroxide (H2O2) and the protonated form of superoxide radicals (HO2). The solution-to-pore transfer free energies were low for small oxy-radicals, suggesting that even highly reactive hydroxyl radicals (HO) might have access to the pore interior and oxidize amino acids responsible for channel selectivity. In the plant aquaporin, no significant change in water permeability was observed upon oxidation of the solvent-exposed disulfide bonds at the extracellular region. During the simulated time scale, the existence of a direct oxidative gating mechanism involving these disulfide bonds could not be demonstrated.

GENERAL SIGNIFICANCE

Simulation results may improve the understanding of redox signaling mechanisms and help in the interpretation of protein oxidative labeling experiments.

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.

Characterization of 3D Voronoi tessellation nearest neighbor lipid shells provides atomistic lipid disruption profile of protein containing lipid membranes

Quantifying protein-induced lipid disruptions at the atomistic level is a challenging problem in membrane biophysics. Here we propose a novel 3D Voronoi tessellation nearest-atom-neighbor shell method to classify and characterize lipid domains into discrete concentric lipid shells surrounding membrane proteins in structurally heterogeneous lipid membranes. This method needs only the coordinates of the system and is independent of force fields and simulation conditions. As a proof-of-principle, we use this multiple lipid shell method to analyze the lipid disruption profiles of three simulated membrane systems: phosphatidylcholine, phosphatidylcholine/cholesterol, and beta-amyloid/phosphatidylcholine/cholesterol. We observed different atomic volume disruption mechanisms due to cholesterol and beta-amyloid. Additionally, several lipid fractional groups and lipid-interfacial water did not converge to their control values with increasing distance or shell order from the protein. This volume divergent behavior was confirmed by bilayer thickness and chain orientational order calculations. Our method can also be used to analyze high-resolution structural experimental data.

Folding and stability of integral membrane proteins in amphipols

Amphipols (APols) are a family of amphipathic polymers designed to keep transmembrane proteins (TMPs) soluble in aqueous solutions in the absence of detergent. APols have proven remarkably efficient at (i) stabilizing TMPs, as compared to detergent solutions, and (ii) folding them from a denatured state to a native, functional one. The underlying physical-chemical mechanisms are discussed.

All-atom molecular dynamics simulation of a photosystem i/detergent complex

All-atom molecular dynamics (MD) simulation was used to investigate the solution structure and dynamics of the photosynthetic pigment-protein complex photosystem I (PSI) from Thermosynechococcus elongatus embedded in a toroidal belt of n-dodecyl-β-d-maltoside (DDM) detergent. Evaluation of root-mean-square deviations (RMSDs) relative to the known crystal structure show that the protein complex surrounded by DDM molecules is stable during the 200 ns simulation time, and root-mean-square fluctuation (RMSF) analysis indicates that regions of high local mobility correspond to solvent-exposed regions such as turns in the transmembrane α-helices and flexible loops on the stromal and lumenal faces. Comparing the protein-detergent complex to a pure detergent micelle, the detergent surrounding the PSI trimer is found to be less densely packed but with more ordered detergent tails, contrary to what is seen in most lipid bilayer models. We also investigated any functional implications for the observed conformational dynamics and protein-detergent interactions, discovering interesting structural changes in the psaL subunits associated with maintaining the trimeric structure of the protein. Importantly, we find that the docking of soluble electron mediators such as cytochrome c6 and ferredoxin to PSI is not significantly impacted by the solubilization of PSI in detergent.

Intact and N- or C-terminal end truncated AQP0 function as open water channels and cell-to-cell adhesion proteins: end truncation could be a prelude for adjusting the refractive index of the lens to prevent spherical aberration

BACKGROUND

Investigate the impact of natural N- or C-terminal post-translational truncations of lens mature fiber cell Aquaporin 0 (AQP0) on water permeability (Pw) and cell-to-cell adhesion (CTCA) functions.

METHODS

The following deletions/truncations were created by site-directed mutagenesis (designations in parentheses): Amino acid residues (AA) 2-6 (AQP0-N-del-2-6), AA235-263 (AQP0-1-234), AA239-263 (AQP0-1-238), AA244-263 (AQP0-1-243), AA247-263 (AQP0-1-246), AA250-263 (AQP0-1-249) and AA260-263 (AQP0-1-259). Protein expression was studied using immunostaining, fluorescent tags and organelle-specific markers. Pw was tested by expressing the respective complementary ribonucleic acid (cRNA) in Xenopus oocytes and conducting osmotic swelling assay. CTCA was assessed by transfecting intact or mutant AQP0 into adhesion-deficient L-cells and performing cell aggregation and adhesion assays.

RESULTS

AQP0-1-234 and AQP0-1-238 did not traffic to the plasma membrane. Trafficking of AQP0-N-del-2-6 and AQP0-1-243 was reduced causing decreased membrane Pw and CTCA. AQP0-1-246, AQP0-1-249 and AQP0-1-259 mutants trafficked properly and functioned normally. Pw and CTCA functions of the mutants were directly proportional to the respective amount of AQP0 expressed at the plasma membrane and remained comparable to those of intact AQP0 (AQP0-1-263).

CONCLUSIONS

Post-translational truncation of N- or C-terminal end amino acids does not alter the basal water permeability of AQP0 or its adhesive functions. AQP0 may play a role in adjusting the refractive index to prevent spherical aberration in the constantly growing lens.

GENERAL SIGNIFICANCE

Similar studies can be extended to other lens proteins which undergo post-translational truncations to find out how they assist the lens to maintain transparency and homeostasis for proper focusing of objects on to the retina.

A theoretical view of protein dynamics

Proteins are fascinating supramolecular structures, which are able to recognize ligands transforming binding information into chemical signals. They can transfer information across the cell, can catalyse complex chemical reactions, and are able to transform energy into work with much more efficiency than any human engine. The unique abilities of proteins are tightly coupled with their dynamic properties, which are coded in a complex way in the sequence and carefully refined by evolution. Despite its importance, our experimental knowledge of protein dynamics is still rather limited, and mostly derived from theoretical calculations. I will review here, in a systematic way, the current state-of-the-art theoretical approaches to the study of protein dynamics, emphasizing the most recent advances, examples of use and the expected lines of development in the near future.

Water permeability of aquaporin-4 channel depends on bilayer composition, thickness, and elasticity

Aquaporin-4 (AQP4) is the primary water channel in the mammalian brain, particularly abundant in astrocytes, whose plasma membranes normally contain high concentrations of cholesterol. Here we test the hypothesis that the water permeabilities of two naturally occurring isoforms (AQP4-M1 and AQP4-M23) depend on bilayer mechanical/structural properties modulated by cholesterol and phospholipid composition. Osmotic stress measurements were performed with proteoliposomes containing AQP4 and three different lipid mixtures: 1), phosphatidylcholine (PC) and phosphatidylglycerol (PG); 2), PC, PG, with 40 mol % cholesterol; and 3), sphingomyelin (SM), PG, with 40 mol % cholesterol. The unit permeabilities of AQP4-M1 were 3.3 ± 0.4 × 10(-13) cm(3)/s (mean ± SE), 1.2 ± 0.1 × 10(-13) cm(3)/s, and 0.4 ± 0.1 × 10(-13) cm(3)/s in PC:PG, PC:PG:cholesterol, and SM:PG:cholesterol, respectively. The unit permeabilities of AQP4-M23 were 2.1 ± 0.2 × 10(-13) cm(3)/s, 0.8 ± 0.1 × 10(-13) cm(3)/s, and 0.3 ± 0.1 × 10(-13) cm(3)/s in PC:PG, PC:PG:cholesterol, and SM:PG:cholesterol, respectively. Thus, for each isoform the unit permeabilities strongly depended on bilayer composition and systematically decreased with increasing bilayer compressibility modulus and bilayer thickness. These observations suggest that altering lipid environment provides a means of regulating water channel permeability. Such permeability changes could have physiological consequences, because AQP4 water permeability would be reduced by its sequestration into SM:cholesterol-enriched raft microdomains. Conversely, under ischemic conditions astrocyte membrane cholesterol content decreases, which could increase AQP4 permeability.

Multiscale simulations reveal conserved patterns of lipid interactions with aquaporins

Interactions of membrane proteins with lipid molecules are central to their stability and function. We have used multiscale molecular dynamics simulations to determine the extent to which interactions with lipids are conserved across the aquaporin (Aqp) family of membrane proteins. Simulation-based assessment of the lipid interactions made by Aqps when embedded within a simple phospholipid bilayer agrees well with the protein-lipid contacts determined by electron diffraction from 2D crystals. Extending this simulation-based analysis to all Aqps of known structure reveals a degree of conservation of such interactions across the Aqp structural proteome. Despite similarities in the binding orientations and interactions of the lipids, there do not appear to be distinct, high-specificity lipid binding sites on the surface of Aqps. Rather Aqps exhibit a more broadly conserved protein/lipid interface, suggestive of interchange between annular and bulk lipids, instead of a fixed annular "shell" of lipids.