Investigating lipid composition effects on the mechanosensitive channel of large conductance (MscL) using molecular dynamics simulations

Dependence of protein-induced lipid bilayer deformations on protein shape

Membrane proteins typically deform the surrounding lipid bilayer membrane, which can play an important role in the function, regulation, and organization of membrane proteins. Membrane elasticity theory provides a beautiful description of protein-induced lipid bilayer deformations, in which all physical parameters can be directly determined from experiments. While analytic solutions of protein-induced elastic bilayer deformations are most easily developed for proteins with approximately circular cross sections, structural biology has shown that membrane proteins come in a variety of distinct shapes, with often considerable deviations from a circular cross section. We develop here a boundary value method (BVM) that permits the construction of analytic solutions of protein-induced elastic bilayer deformations for protein shapes with arbitrarily large deviations from a circular cross section, for constant as well as variable boundary conditions along the bilayer-protein interface. We apply this BVM to protein-induced lipid bilayer thickness deformations. Our BVM reproduces available analytic solutions for proteins with circular cross section and yields, for proteins with noncircular cross section, excellent agreement with numerical, finite element solutions. On this basis, we formulate a simple analytic approximation of the bilayer thickness deformation energy associated with general protein shapes and show that, for modest deviations from rotational symmetry, this analytic approximation is in good agreement with BVM solutions. Using the BVM, we survey the dependence of protein-induced elastic bilayer thickness deformations on protein shape, and thus explore how the coupling of protein shape and bilayer thickness deformations affects protein oligomerization and transitions in protein conformational state.

Decisive structural elements in water and ion permeation through mechanosensitive channels of large conductance: insights from molecular dynamics simulation

In this paper, a series of equilibrium molecular dynamics simulations (EMD), steered molecular dynamics (SMD), and computational electrophysiology methods are carried out to explore water and ion permeation through mechanosensitive channels of large conductance (MscL). This research aims to identify the pore-lining side chains of the channel in different conformations of MscL homologs by analyzing the pore size. The distribution of permeating water dipole angles through the pore domains enclosed by VAL21 and GLU104 demonstrated that water molecules are oriented toward the charged oxygen headgroups of GLU104 from their hydrogen atoms to retain this interaction in a stabilized fashion. Although, this behavior was not perceived for VAL21. Numerical assessments of the secondary structure clarified that, during the ion permeation, in addition to the secondary structure alterations, the structure of Tb-MscL would also undergo significant conformational changes. It was elucidated that VAL21, GLU104, and water molecules accomplish a fundamental task in ion permeation. The mentioned residues hinder ion permeation so that the pulling SMD force is increased remarkably when the ions permeate through the domains enclosed by VAL21 and GLU102. The hydration level and potassium diffusivity in the hydrophobic gate of the transmembrane domain were promoted by applying the external electric field. Furthermore, the implementation of an external electric field altered the distribution pattern for potassium ions in the system while intensifying the accumulation of Cl⁻ in the vicinity of ARG11 and ARG98.

Graphical representation of the most determinant pore-lining side chains of Tb-MscL along with the solid surfaces depicting the spatial shape of the interior pore.

Synthetic cells in biomedical applications

Synthetic cells are engineered vesicles that can mimic one or more salient features of life. These features include directed localization, sense-and-respond behavior, gene expression, metabolism, and high stability. In nanomedicine, many of these features are desirable capabilities of drug delivery vehicles but are difficult to engineer. In this focus article, we discuss where synthetic cells offer unique advantages over nanoparticle and living cell therapies. We review progress in the engineering of the above life-like behaviors and how they are deployed in nanomedicine. Finally, we assess key challenges synthetic cells face before being deployed as drugs and suggest ways to overcome these challenges. This article is categorized under: Therapeutic Approaches and Drug Discovery > Emerging Technologies Biology-Inspired Nanomaterials > Lipid-Based Structures.

A Single Mutation in the Outer Lipid-Facing Helix of a Pentameric Ligand-Gated Ion Channel Affects Channel Function Through a Radially-Propagating Mechanism

Pentameric ligand-gated ion channels (pLGICs) mediate fast synaptic transmission and are crucial drug targets. Their gating mechanism is triggered by ligand binding in the extracellular domain that culminates in the opening of a hydrophobic gate in the transmembrane domain. This domain is made of four α-helices (M1 to M4). Recently the outer lipid-facing helix (M4) has been shown to be key to receptor function, however its role in channel opening is still poorly understood. It could act through its neighboring helices (M1/M3), or via the M4 tip interacting with the pivotal Cys-loop in the extracellular domain. Mutation of a single M4 tyrosine (Y441) to alanine renders one pLGIC-the 5-HT3A receptor-unable to function despite robust ligand binding. Using Y441A as a proxy for M4 function, we here predict likely paths of Y441 action using molecular dynamics, and test these predictions with functional assays of mutant receptors in HEK cells and Xenopus oocytes using fluorescent membrane potential sensitive dye and two-electrode voltage clamp respectively. We show that Y441 does not act via the M4 tip or Cys-loop, but instead connects radially through M1 to a residue near the ion channel hydrophobic gate on the pore-lining helix M2. This demonstrates the active role of the M4 helix in channel opening.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

When Stiffness Matters: Mechanosensing in Heart Development and Disease

During embryonic morphogenesis, the heart undergoes a complex series of cellular phenotypic maturations (e.g., transition of myocytes from proliferative to quiescent or maturation of the contractile apparatus), and this involves stiffening of the extracellular matrix (ECM) acting in concert with morphogenetic signals. The maladaptive remodeling of the myocardium, one of the processes involved in determination of heart failure, also involves mechanical cues, with a progressive stiffening of the tissue that produces cellular mechanical damage, inflammation, and ultimately myocardial fibrosis. The assessment of the biomechanical dependence of the molecular machinery (in myocardial and non-myocardial cells) is therefore essential to contextualize the maturation of the cardiac tissue at early stages and understand its pathologic evolution in aging. Because systems to perform multiscale modeling of cellular and tissue mechanics have been developed, it appears particularly novel to design integrated mechano-molecular models of heart development and disease to be tested in ex vivo reconstituted cells/tissue-mimicking conditions. In the present contribution, we will discuss the latest implication of mechanosensing in heart development and pathology, describe the most recent models of cell/tissue mechanics, and delineate novel strategies to target the consequences of heart failure with personalized approaches based on tissue engineering and induced pluripotent stem cell (iPSC) technologies.

Influence of the Headgroup on the Interaction of Poly(ethylene oxide)-Poly(propylene oxide) Block Copolymers with Lipid Bilayers

The lipid headgroup plays an important role in the association of polymers with lipid bilayer membranes. Herein, we report how a glycerol headgroup versus a choline headgroup affects the interaction of poly(ethylene oxide)-b-poly(propylene oxide) (PEO-PPO) block copolymers with lipid bilayer vesicles. Unilamellar vesicles composed of phosphatidylcholine and phosphatidylglycerol at various molar ratios were used as model membranes. The interactions between the block copolymers and lipid bilayers were quantified by pulsed-field gradient nuclear magnetic resonance (PFG-NMR) based on the distinctly different mobilities of free and bound polymers. All the investigated polymer species showed significantly higher binding with 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) sodium salt (POPG) liposomes than with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) liposomes, indicating stronger association with the glycerol headgroup compared to the choline headgroup. This effect did not become significant until the composition of mixed POPC/POPG liposomes contained more than 20 mol % POPG. A plausible explanation for the enhanced polymer binding with POPG invokes the role of hydrogen bonding between the glycerol headgroup and the ether moieties of the polymers.

Applications of Molecular Simulation in the Discovery of Antituberculosis Drugs: A Review

After decades of efforts, tuberculosis has been well controlled in most places. The existing drugs are no longer sufficient for the treatment of drug-resistant Mycobacterium tuberculosis due to significant toxicity and selective pressure, especially for XDR-TB. In order to accelerate the development of high-efficiency, low-toxic antituberculosis drugs, it is particularly important to use Computer Aided Drug Design (CADD) for rational drug design. Here, we systematically reviewed the specific role of molecular simulation in the discovery of new antituberculosis drugs. The purpose of this review is to overview current applications of molecular simulation methods in the discovery of antituberculosis drugs. Furthermore, the unique advantages of molecular simulation was discussed in revealing the mechanism of drug resistance. The comprehensive use of different molecular simulation methods will help reveal the mechanism of drug resistance and improve the efficiency of rational drug design. With the help of molecular simulation methods such as QM/MM method, the mechanisms of biochemical reactions catalyzed by enzymes at atomic level in Mycobacterium tuberculosis has been deeply analyzed. QSAR and virtual screening both accelerate the development of highefficiency, low-toxic potential antituberculosis drugs. Improving the accuracy of existing algorithms and developing more efficient new methods for CADD will always be a hot topic in the future. It is of great value to utilize molecular dynamics simulation to investigate complex systems that cannot be studied in experiments, especially for drug resistance of Mycobacterium tuberculosis.

Biophysical Mechanisms of Membrane-Thickness-Dependent MscL Gating: An All-Atom Molecular Dynamics Study

The bacterial mechanosensitive channel, MscL, is activated by membrane tension, acting as a safety valve to prevent cell lysis against hypotonic challenge. It has been established that its activation threshold decreases with membrane thickness, while the underlying mechanism remains to be solved. We performed all-atom molecular dynamics (MD) simulations for the initial opening process of MscL embedded in four different types of lipid bilayers with different thicknesses: 1,2-dilauroyl- sn-glycero-3-phosphocholine (DLPC)), 1,2-dimyristoyl-glycero-3-phosphorylcholine (DMPC), 1,2-dipalmitoyl- sn-glycero-3-phosphocholine (DPPC), and 1,2-distearoyl- sn-glycero-3-phosphocholine (DSPC). In response to membrane stretching, channel opening occurred only in the thinner membranes (DLPC and DMPC) in a thickness-dependent way. We found that the MscL opening was governed by the rate and degree of membrane thinning and that the channel opening was tightly associated with the tilting of transmembrane (TM) helices of MscL toward the membrane plane. Upon membrane stretching, the order parameter of acyl chains of thinner membranes (DLPC and DMPC) became smaller, whereas other thicker membranes (DPPC and DSPC) showed interdigitation with little changes in the order parameter. The decreased order parameter contributed much more to membrane thinning than did interdigitation. We conclude that the membrane-thickness-dependent MscL opening mainly arises from structural changes in MscL to match the altered membrane thickness by stretching.

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.

Membrane Thickness as a Key Factor Contributing to the Activation of Osmosensors and Essential Ras Signaling Pathways

The cell membrane provides a functional link between the external environment and the replicating DNA genome by using ligand-gated receptors and chemical signals to activate signaling transduction pathways. However, increasing evidence has also indicated that the phospholipid bilayer itself by altering various physical parameters serves as a sensor that regulate membrane proteins in a specific manner. Changes in thickness and/or curvature of the membrane have been shown to be induced by mechanical forces and transmitted through the transmembrane helices of several types of mechanosensitive (MS) ion channels underlying functions such as osmoregulation in bacteria and sensory processing in mammalian cells. This review focus on recent protein functional and structural data indicating that the activation of bacterial and yeast osmosensors is consistent with thickness-induced tilting changes of the transmembrane domains of these proteins. Membrane thinning in combination with curvature changes may also lead to the lateral transfer of the small lipid-anchored GTPases Ras1 and H-Ras out of lipid rafts for clustering and signaling. The modulation of signaling pathways by amphiphilic peptides and the membrane-active antibiotics colistin and Amphotericin B is also discussed.

Bilayer-Mediated Structural Transitions Control Mechanosensitivity of the TREK-2 K2P Channel

The mechanosensitive two-pore domain (K2P) K⁺ channels (TREK-1, TREK-2, and TRAAK) are important for mechanical and thermal nociception. However, the mechanisms underlying their gating by membrane stretch remain controversial. Here we use molecular dynamics simulations to examine their behavior in a lipid bilayer. We show that TREK-2 moves from the “down” to “up” conformation in direct response to membrane stretch, and examine the role of the transmembrane pressure profile in this process. Furthermore, we show how state-dependent interactions with lipids affect the movement of TREK-2, and how stretch influences both the inner pore and selectivity filter. Finally, we present functional studies that demonstrate why direct pore block by lipid tails does not represent the principal mechanism of mechanogating. Overall, this study provides a dynamic structural insight into K2P channel mechanosensitivity and illustrates how the structure of a eukaryotic mechanosensitive ion channel responds to changes in forces within the bilayer.

•

Mechanogating of TREK-2 involves movement from the down to up conformation

•

Simulations sample a wide range of mechanosensitive K2P channel structures

•

Changes in the pressure profile and state-dependent lipid interactions play a key role

•

Lipid block of the inner pore does not mediate stretch activation

Dihydrostreptomycin Directly Binds to, Modulates, and Passes through the MscL Channel Pore

The primary mechanism of action of the antibiotic dihydrostreptomycin is binding to and modifying the function of the bacterial ribosome, thus leading to decreased and aberrant translation of proteins; however, the routes by which it enters the bacterial cell are largely unknown. The mechanosensitive channel of large conductance, MscL, is found in the vast majority of bacterial species, where it serves as an emergency release valve rescuing the cell from sudden decreases in external osmolarity. While it is known that MscL expression increases the potency of dihydrostreptomycin, it has remained unclear if this effect is due to a direct interaction. Here, we use a combination of genetic screening, MD simulations, and biochemical and mutational approaches to determine if dihydrostreptomycin directly interacts with MscL. Our data strongly suggest that dihydrostreptomycin binds to a specific site on MscL and modifies its conformation, thus allowing the passage of K+ and glutamate out of, and dihydrostreptomycin into, the cell.

Bilayer-thickness-mediated interactions between integral membrane proteins

Hydrophobic thickness mismatch between integral membrane proteins and the surrounding lipid bilayer can produce lipid bilayer thickness deformations. Experiment and theory have shown that protein-induced lipid bilayer thickness deformations can yield energetically favorable bilayer-mediated interactions between integral membrane proteins, and large-scale organization of integral membrane proteins into protein clusters in cell membranes. Within the continuum elasticity theory of membranes, the energy cost of protein-induced bilayer thickness deformations can be captured by considering compression and expansion of the bilayer hydrophobic core, membrane tension, and bilayer bending, resulting in biharmonic equilibrium equations describing the shape of lipid bilayers for a given set of bilayer-protein boundary conditions. Here we develop a combined analytic and numerical methodology for the solution of the equilibrium elastic equations associated with protein-induced lipid bilayer deformations. Our methodology allows accurate prediction of thickness-mediated protein interactions for arbitrary protein symmetries at arbitrary protein separations and relative orientations. We provide exact analytic solutions for cylindrical integral membrane proteins with constant and varying hydrophobic thickness, and develop perturbative analytic solutions for noncylindrical protein shapes. We complement these analytic solutions, and assess their accuracy, by developing both finite element and finite difference numerical solution schemes. We provide error estimates of our numerical solution schemes and systematically assess their convergence properties. Taken together, the work presented here puts into place an analytic and numerical framework which allows calculation of bilayer-mediated elastic interactions between integral membrane proteins for the complicated protein shapes suggested by structural biology and at the small protein separations most relevant for the crowded membrane environments provided by living cells.

Computational Investigation of the Effect of Lipid Membranes on Ion Permeation in Gramicidin A

Membrane proteins are embedded in a lipid bilayer and interact with the lipid molecules in subtle ways. This can be studied experimentally by examining the effect of different lipid bilayers on the function of membrane proteins. Understanding the causes of the functional effects of lipids is difficult to dissect experimentally but more amenable to a computational approach. Here we perform molecular dynamics simulations and free energy calculations to study the effect of two lipid types (POPC and NODS) on the conductance of the gramicidin A (gA) channel. A larger energy barrier is found for the K⁺ potential of mean force in gA embedded in POPC compared to that in NODS, which is consistent with the enhanced experimental conductance of cations in gA embedded in NODS. Further analysis of the contributions to the potential energy of K⁺ reveals that gA and water molecules in gA make similar contributions in both bilayers but there are significant differences between the two bilayers when the lipid molecules and interfacial waters are considered. It is shown that the stronger dipole moments of the POPC head groups create a thicker layer of interfacial waters with better orientation, which ultimately is responsible for the larger energy barrier in the K⁺ PMF in POPC.

The Combined Effect of Hydrophobic Mismatch and Bilayer Local Bending on the Regulation of Mechanosensitive Ion Channels

The hydrophobic mismatch between the lipid bilayer and integral membrane proteins has well-defined effect on mechanosensitive (MS) ion channels. Also, membrane local bending is suggested to modulate MS channel activity. Although a number of studies have already shown the significance of each individual factor, the combined effect of these physical factors on MS channel activity have not been investigated. Here using finite element simulation, we study the combined effect of hydrophobic mismatch and local bending on the archetypal mechanosensitive channel MscL. First we show how the local curvature direction impacts on MS channel modulation. In the case of MscL, we show inward (cytoplasmic) bending can more effectively gate the channel compared to outward bending. Then we indicate that in response to a specific local curvature, MscL inserted in a bilayer with the same hydrophobic length is more expanded in the constriction pore region compared to when there is a protein-lipid hydrophobic mismatch. Interestingly in the presence of a negative mismatch (thicker lipids), MscL constriction pore is more expanded than in the presence of positive mismatch (thinner lipids) in response to an identical membrane curvature. These results were confirmed by a parametric energetic calculation provided for MscL gating. These findings have several biophysical consequences for understanding the function of MS channels in response to two major physical stimuli in mechanobiology, namely hydrophobic mismatch and local membrane curvature.

The efficacy of trivalent cyclic hexapeptides to induce lipid clustering in PG/PE membranes correlates with their antimicrobial activity

Various models have been proposed for the sequence of events occurring after binding of specific antimicrobial peptides to lipid membranes. The lipid clustering model arose by the finding that antimicrobial peptides can induce a segregation of certain negatively charged lipids in lipid model membranes. Anionic lipid segregation by cationic peptides is initially an effect of charge interaction where the ratio of peptide and lipid charges is thought to be the decisive parameter in the peptide induced lipid demixing. However, the sequence of events following this initial lipid clustering is more complex and can lead to deactivation of membrane proteins involved in cell division or perturbation of lipid reorganization essential for cell division. In this study we used DSC and ITC techniques to investigate the effect of binding different cyclic hexapeptides with varying antimicrobial efficacy, to phosphatidylglycerol (PG)/phosphatidylethanolamine (PE) lipid membranes and their ability to induce lipid segregation in these mixtures. We found that these cyclic hexapeptides consisting of three charged and three aromatic amino acids showed indeed different abilities to induce lipid demixing depending on their amino acid composition and their sequence. The results clearly showed that the cationic amino acids are essential for electrostatic binding but that the three hydrophobic amino acids in the peptides and their position in the sequence also contribute to binding affinity and to the extent of induction of lipid clustering. The efficacy of these different hexapeptides to induce PG clusters in PG/PE membranes was found to be correlated with their antimicrobial activity.

Lipid-protein interactions: Lessons learned from stress

Biological membranes are essential for normal function and regulation of cells, forming a physical barrier between extracellular and intracellular space and cellular compartments. These physical barriers are subject to mechanical stresses. As a consequence, nature has developed proteins that are able to transpose mechanical stimuli into meaningful intracellular signals. These proteins, termed Mechanosensitive (MS) proteins provide a variety of roles in response to these stimuli. In prokaryotes these proteins form transmembrane spanning channels that function as osmotically activated nanovalves to prevent cell lysis by hypoosmotic shock. In eukaryotes, the function of MS proteins is more diverse and includes physiological processes such as touch, pain and hearing. The transmembrane portion of these channels is influenced by the physical properties such as charge, shape, thickness and stiffness of the lipid bilayer surrounding it, as well as the bilayer pressure profile. In this review we provide an overview of the progress to date on advances in our understanding of the intimate biophysical and chemical interactions between the lipid bilayer and mechanosensitive membrane channels, focusing on current progress in both eukaryotic and prokaryotic systems. These advances are of importance due to the increasing evidence of the role the MS channels play in disease, such as xerocytosis, muscular dystrophy and cardiac hypertrophy. Moreover, insights gained from lipid-protein interactions of MS channels are likely relevant not only to this class of membrane proteins, but other bilayer embedded proteins as well. This article is part of a Special Issue entitled: Lipid-protein interactions.

Breaking the hydrophobicity of the MscL pore: insights into a charge-induced gating mechanism

The mechanosensitive channel of large conductance (MscL) is a protein that responds to membrane tension by opening a transient pore during osmotic downshock. Due to its large pore size and functional reconstitution into lipid membranes, MscL has been proposed as a promising artificial nanovalve suitable for biotechnological applications. For example, site-specific mutations and tailored chemical modifications have shown how MscL channel gating can be triggered in the absence of tension by introducing charged residues at the hydrophobic pore level. Recently, engineered MscL proteins responsive to stimuli like pH or light have been reported. Inspired by experiments, we present a thorough computational study aiming at describing, with atomistic detail, the artificial gating mechanism and the molecular transport properties of a light-actuated bacterial MscL channel, in which a charge-induced gating mechanism has been enabled through the selective cleavage of photo-sensitive alkylating agents. Properties such as structural transitions, pore dimension, ion flux and selectivity have been carefully analyzed. Besides, the effects of charge on alternative sites of the channel with respect to those already reported have been addressed. Overall, our results provide useful molecular insights into the structural events accompanying the engineered MscL channel gating and the interplay of electrostatic effects, channel opening and permeation properties. In addition, we describe how the experimentally observed ionic current in a single-subunit charged MscL mutant is obtained through a hydrophobicity breaking mechanism involving an asymmetric inter-subunit motion.

The mechanosensitive channel of small conductance (MscS) functions as a Jack-in-the box

Phenotypical analysis of the lipid interacting residues in the closed state of the mechanosensitive channel of small conductance (MscS) from Escherichia coli (E. coli) has previously shown that these residues are critical for channel function. In the closed state, mutation of individual hydrophobic lipid lining residues to alanine, thus reducing the hydrophobicity, resulted in phenotypic changes that were observable using in vivo assays. Here, in an analogous set of experiments, we identify eleven residues in the first transmembrane domain of the open state of MscS that interact with the lipid bilayer. Each of these residues was mutated to alanine and leucine to modulate their hydrophobic interaction with the lipid tail-groups in the open state. The effects of these changes on channel function were analyzed using in vivo bacterial assays and patch clamp electrophysiology. Mutant channels were found to be functionally indistinguishable from wildtype MscS. Thus, mutation of open-state lipid interacting residues does not differentially stabilize or destabilize the open, closed, intermediate, or transition states of MscS. Based on these results and other data from the literature, we propose a new gating paradigm for MscS where MscS acts as a "Jack-In-The-Box" with the intrinsic bilayer lateral pressure holding the channel in the closed state. In this model, upon application of extrinsic tension the channel springs into the open state due to relief of the intrinsic lipid bilayer pressure.

Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE

In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein's β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.

Bacterial mechanosensitive channels: models for studying mechanosensory transduction

SIGNIFICANCE

Sensations of touch and hearing are manifestations of mechanical contact and air pressure acting on touch receptors and hair cells of the inner ear, respectively. In bacteria, osmotic pressure exerts a significant mechanical force on their cellular membrane. Bacteria have evolved mechanosensitive (MS) channels to cope with excessive turgor pressure resulting from a hypo-osmotic shock. MS channel opening allows the expulsion of osmolytes and water, thereby restoring normal cellular turgor and preventing cell lysis.

RECENT ADVANCES

As biological force-sensing systems, MS channels have been identified as the best examples of membrane proteins coupling molecular dynamics to cellular mechanics. The bacterial MS channel of large conductance (MscL) and MS channel of small conductance (MscS) have been subjected to extensive biophysical, biochemical, genetic, and structural analyses. These studies have established MscL and MscS as model systems for mechanosensory transduction.

CRITICAL ISSUES

In recent years, MS ion channels in mammalian cells have moved into focus of mechanotransduction research, accompanied by an increased awareness of the role they may play in the pathophysiology of diseases, including cardiac hypertrophy, muscular dystrophy, or Xerocytosis.

FUTURE DIRECTIONS

A recent exciting development includes the molecular identification of Piezo proteins, which function as nonselective cation channels in mechanosensory transduction associated with senses of touch and pain. Since research on Piezo channels is very young, applying lessons learned from studies of bacterial MS channels to establishing the mechanism by which the Piezo channels are mechanically activated remains one of the future challenges toward a better understanding of the role that MS channels play in mechanobiology.

Connection between oligomeric state and gating characteristics of mechanosensitive ion channels

The mechanosensitive channel of large conductance (MscL) is capable of transducing mechanical stimuli such as membrane tension into an electrochemical response. MscL provides a widely-studied model system for mechanotransduction and, more generally, for how bilayer mechanical properties regulate protein conformational changes. Much effort has been expended on the detailed experimental characterization of the molecular structure and biological function of MscL. However, despite its central significance, even basic issues such as the physiologically relevant oligomeric states and molecular structures of MscL remain a matter of debate. In particular, tetrameric, pentameric, and hexameric oligomeric states of MscL have been proposed, together with a range of detailed molecular structures of MscL in the closed and open channel states. Previous theoretical work has shown that the basic phenomenology of MscL gating can be understood using an elastic model describing the energetic cost of the thickness deformations induced by MscL in the surrounding lipid bilayer. Here, we generalize this elastic model to account for the proposed oligomeric states and hydrophobic shapes of MscL. We find that the oligomeric state and hydrophobic shape of MscL are reflected in the energetic cost of lipid bilayer deformations. We make quantitative predictions pertaining to the gating characteristics associated with various structural models of MscL and, in particular, show that different oligomeric states and hydrophobic shapes of MscL yield distinct membrane contributions to the gating energy and gating tension. Thus, the functional properties of MscL provide a signature of the oligomeric state and hydrophobic shape of MscL. Our results suggest that, in addition to the hydrophobic mismatch between membrane proteins and the surrounding lipid bilayer, the symmetry and shape of the hydrophobic surfaces of membrane proteins play an important role in the regulation of protein function by bilayer membranes.

Controlled delivery of bioactive molecules into live cells using the bacterial mechanosensitive channel MscL

Bacterial mechanosensitive channels are some of the largest pores in nature. In particular, MscL, with a pore diameter >25 Å, allows passage of large organic ions and small proteins. Functional MscL reconstitution into lipids has been proposed for applications in vesicular-based drug release. Here we show that these channels can be functionally expressed in mammalian cells to afford rapid controlled uptake of membrane-impermeable molecules. We first demonstrate that MscL gating in response to increased membrane tension is preserved in mammalian cell membranes. Molecular delivery is controlled by adopting an established method of MscL charge-induced activation. We then determine pore size limitations using fluorescently labelled model cargoes. Finally, we activate MscL to introduce the cell-impermeable bi-cyclic peptide phalloidin, a specific marker for actin filaments, into cells. We propose that MscL will be a useful tool for gated and controlled delivery of bioactive molecules into cells.

The determinants of hydrophobic mismatch response for transmembrane helices

Hydrophobic mismatch arises from a difference in the hydrophobic thickness of a lipid membrane and a transmembrane protein segment, and is thought to play an important role in the folding, stability and function of membrane proteins. We have investigated the possible adaptations that lipid bilayers and transmembrane α-helices undergo in response to mismatch, using fully-atomistic molecular dynamics simulations totaling 1.4 μs. We have created 25 different tryptophan-alanine-leucine transmembrane α-helical peptide systems, each composed of a hydrophobic alanine-leucine stretch, flanked by 1-4 tryptophan side chains, as well as the β-helical peptide dimer, gramicidin A. Membrane responses to mismatch include changes in local bilayer thickness and lipid order, varying systematically with peptide length. Adding more flanking tryptophan side chains led to an increase in bilayer thinning for negatively mismatched peptides, though it was also associated with a spreading of the bilayer interface. Peptide tilting, bending and stretching were systematic, with tilting dominating the responses, with values of up to ~45° for the most positively mismatched peptides. Peptide responses were modulated by the number of tryptophan side chains due to their anchoring roles and distributions around the helices. Potential of mean force calculations for local membrane thickness changes, helix tilting, bending and stretching revealed that membrane deformation is the least energetically costly of all mismatch responses, except for positively mismatched peptides where helix tilting also contributes substantially. This comparison of energetic driving forces of mismatch responses allows for deeper insight into protein stability and conformational changes in lipid membranes.

Defining the role of the tension sensor in the mechanosensitive channel of small conductance

Mutations that alter the phenotypic behavior of the Escherichia coli mechanosensitive channel of small conductance (MscS) have been identified; however, most of these residues play critical roles in the transition between the closed and open states of the channel and are not directly involved in lipid interactions that transduce the tension response. In this study, we use molecular dynamic simulations to predict critical lipid interacting residues in the closed state of MscS. The physiological role of these residues was then investigated by performing osmotic downshock assays on MscS mutants where the lipid interacting residues were mutated to alanine. These experiments identified seven residues in the first and second transmembrane helices as lipid-sensing residues. The majority of these residues are hydrophobic amino acids located near the extracellular interface of the membrane. All of these residues interact strongly with the lipid bilayer in the closed state of MscS, but do not face the bilayer directly in structures associated with the open and desensitized states of the channel. Thus, the position of these residues relative to the lipid membrane appears related to the ability of the channel to sense tension in its different physiological states.

An in vivo screen reveals protein-lipid interactions crucial for gating a mechanosensitive channel

The bacterial mechanosensitive channel MscL is the best-studied mechanosensor, thus serving as a paradigm of how a protein senses and responds to mechanical force. Models for the transition of Escherichia coli MscL from closed to open states propose a tilting of the transmembrane domains in the plane of the membrane, suggesting dynamic protein-lipid interactions. Here, we used a rapid in vivo assay to assess the function of channels that were post-translationally modified at several different sites in a region just distal to the cytoplasmic end of the second transmembrane helix. We utilized multiple probes with various affinities for the membrane environment. The in vivo functional data, combined with site-directed mutagenesis, single-channel analyses, and tryptophan fluorescence measurements, confirmed that lipid interactions within this region are critical for MscL gating. The data suggest a model in which this region acts as an anchor for the transmembrane domain tilting during gating. Furthermore, the conservation of analogous motifs among many other channels suggests a conserved protein-lipid dynamic mechanism.

Mechanosensitivity of ion channels based on protein-lipid interactions

Ion channels form a group of membrane proteins that pass ions through a pore beyond the energy barrier of the lipid bilayer. The structure of the transmembrane segment of membrane proteins is influenced by the charges and the hydrophobicity of the surrounding lipids and the pressure on its surface. A mechanosensitive channel is specifically designed to change its conformation in response to changes in the membrane pressure (tension). However, mechanosensitive channels are not the only group that is sensitive to the physical environment of the membrane: voltage-gated channels are also amenable to the lipid environment. In this article, we review the structure and gating mechanisms of the mechanosensitive channels and voltage-gated channels and discuss how their functions are affected by the physical properties of the lipid bilayer.

Mechanosensitivity of ion channels based on protein-lipid interactions

Ion channels form a group of membrane proteins that pass ions through a pore beyond the energy barrier of the lipid bilayer. The structure of the transmembrane segment of membrane proteins is influenced by the charges and the hydrophobicity of the surrounding lipids and the pressure on its surface. A mechanosensitive channel is specifically designed to change its conformation in response to changes in the membrane pressure (tension). However, mechanosensitive channels are not the only group that is sensitive to the physical environment of the membrane: voltage-gated channels are also amenable to the lipid environment. In this article, we review the structure and gating mechanisms of the mechanosensitive channels and voltage-gated channels and discuss how their functions are affected by the physical properties of the lipid bilayer.

Screening for the drug-phospholipid interaction: correlation to phospholipidosis

Phospholipid bilayers represent a complex, anisotropic environment fundamentally different from bulk oil or octanol, for instance. Even "simple" drug association to phospholipid bilayers can only be fully understood if the slab-of-hydrocarbon approach is abandoned and the complex, anisotropic properties of lipid bilayers reflecting the chemical structures and organization of the constituent phospholipids are considered. The interactions of drugs with phospholipids are important in various processes, such as drug absorption, tissue distribution, and subcellular distribution. In addition, drug-lipid interactions may lead to changes in lipid-dependent protein activities, and further, to functional and morphological changes in cells, a prominent example being the phospholipidosis (PLD) induced by cationic amphiphilic drugs. Herein we briefly review drug-lipid interactions in general and the significance of these interactions in PLD in particular. We also focus on a potential causal connection between drug-induced PLD and steatohepatitis, which is induced by some cationic amphiphilic drugs.

Emerging roles for lipids in shaping membrane-protein function

Studies of membrane proteins have revealed a direct link between the lipid environment and the structure and function of some of these proteins. Although some of these effects involve specific chemical interactions between lipids and protein residues, many can be understood in terms of protein-induced perturbations to the membrane shape. The free-energy cost of such perturbations can be estimated quantitatively, and measurements of channel gating in model systems of membrane proteins with their lipid partners are now confirming predictions of simple models.

Effect of lipid composition on buforin II structure and membrane entry

Buforin II is a 21-amino acid polycationic antimicrobial peptide derived from a peptide originally isolated from the stomach tissue of the Asian toad Bufo bufo gargarizans. It is hypothesized to target a wide range of bacteria by translocating into cells without membrane permeabilization and binding to nucleic acids. Previous research found that the structure and membrane interactions of buforin II are related to lipid composition. In this study, we used molecular dynamics (MD) simulations along with lipid vesicle experiments to gain insight into how buforin II interacts differently with phosphatidylcholine (PC), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE) lipids. Fluorescent spectroscopic measurements agreed with the previous assertion that buforin II does not interact with pure PC vesicles. Nonetheless, the reduced entry of the peptide into anionic PG membranes versus neutral PC membranes during simulations correlates with the experimentally observed reduction in BF2 translocation through pure PG membranes. Simulations showing membrane entry into PC also provide insight into how buforin II may initially penetrate cell membranes. Our MD simulations also allowed us to consider how neutral PE lipids affect the peptide differently than PC. In particular, the peptide had a more helical secondary structure in simulations with PE lipids. A change in structure was also apparent in circular dichroism measurements. PE also reduced membrane entry in simulations, which correlates with decreased translocation in the presence of PE observed in previous studies. Together, these results provide molecular-level insight into how lipid composition can affect buforin II structure and function and will be useful in efforts to design peptides with desired antimicrobial and cell-penetrating properties.

Polyunsaturated fatty acid modulation of voltage-gated ion channels

Arachidonic acid (AA) was found to inhibit the function of whole-cell voltage-gated (VG) calcium currents nearly 16 years ago. There are now numerous examples demonstrating that AA and other polyunsaturated fatty acids (PUFAs) modulate the function of VG ion channels, primarily in neurons and muscle cells. We will review and extract some common features about the modulation by PUFAs of VG calcium, sodium, and potassium channels and discuss the impact of this modulation on the excitability of neurons and cardiac myocytes. We will describe the fatty acid nature of the membrane, how fatty acids become available to function as modulators of VG channels, and the physiologic importance of this type of modulation. We will review the evidence for molecular mechanisms and assess our current understanding of the structural basis for modulation. With guidance from research on the structure of fatty acid binding proteins, the role of lipids in gating mechanosensitive (MS) channels, and the impact of membrane lipid composition on membrane-embedded proteins, we will highlight some avenues for future investigations.

Simulations of lipid adsorption on TiO2 surfaces in solution

Molecular dynamics simulations are carried out to study the adsorption of three lipids, namely, DOPC, DOPS, and DMTAP, on TiO2(110) rutile surfaces and the influence of the interface on their conformational properties. Three types of rutile (110) surfaces, characterized by a different degree of hydroxylation (the neutral nonhydroxylated and hydroxylated surfaces and a partially hydroxylated surface with charge density corresponding to physiological pH) are investigated using force fields derived from ab initio calculations and experimental data. It is found that the stability of the adsorbate and the strength of the attachment are strictly connected with the nature of both the lipid and the surface. Direct coordination of the phosphate or carbonyl oxygens of the lipids with available titanium sites, observed in the case of partially or nonhydroxylated layers, determines stronger adsorption and, as a consequence, reduced dynamics. For a given hydration state of the surface, the adsorption strengths are in the order DOPS > DOPC >> DMTAP, in agreement with experimental data according to which the presence of DOPS units inside lipid bilayers favors stronger adsorption and lower mobility. The adsorption geometry, the hydration state of the lipid headgroups, and the dynamical processes (detachment, diffusion, etc.) occurring at the lipid/oxide interface are analyzed in detail, putting on a roughly quantitative basis time scales and energy barriers of the latter processes.

Mechanosensitive channels: insights from continuum-based simulations

Mechanotransduction plays an important role in regulating cell functions and it is an active topic of research in biophysics. Despite recent advances in experimental and numerical techniques, the intrinsic multiscale nature imposes tremendous challenges for revealing the working mechanisms of mechanosensitive channels. Recently, a continuum-mechanics-based hierarchical modeling and simulation framework has been established and applied to study the mechanical responses and gating behaviors of a prototypical mechanosensitive channel, the mechanosensitive channel of large conductance (MscL) in bacteria Escherichia coli (E. coli), from which several putative gating mechanisms have been tested and new insights are deduced. This article reviews these latest findings using the continuum mechanics framework and suggests possible improvements for future simulation studies. This computationally efficient and versatile continuum-mechanics-based protocol is poised to make contributions to the study of a variety of mechanobiology problems.

Bacterial mechanosensitive channels: Experiment and theory

Since their discovery in Escherichia coli some 20 years ago, studies of bacterial mechanosensitive (MS) ion channels have been at the forefront of the MS channel research field. Two major events greatly advanced the research on bacterial MS channels: (i) cloning of MscL and MscS, the MS channels of Large and Small conductance, and (ii) solving their 3D crystal structure. These events enabled further experimental studies employing EPR and FRET spectroscopy in addition to patch clamp and molecular biological techniques that have successfully been used in characterization of the structure and function of bacterial MS channels. In parallel with the experimental studies computational modelling has been applied to elucidate the molecular dynamics of MscL and MscS, which has significantly contributed to our understanding of basic physical principles of the mechanosensory transduction in living organisms.

Mechanosensitive membrane channels in action

The tension-driven gating process of MscL from Mycobacterium tuberculosis, Tb-MscL, has been addressed at near-atomic detail using coarse-grained molecular dynamics simulations. To perform the simulations, a novel coarse-grained peptide model based on a thermodynamic parameterization of the amino-acid side chains has been applied. Both the wild-type Tb-MscL and its gain-of-function mutant V21D embedded in a solvated lipid bilayer have been studied. To mimic hypoosmotic shock conditions, simulations were performed at increasing levels of membrane tension approaching the rupture threshold of the lipid bilayer. Both the wild-type and the mutant channel are found to undergo significant conformational changes in accordance with an irislike expansion mechanism, reaching a conducting state on a microsecond timescale. The most pronounced expansion of the pore has been observed for the V21D mutant, which is consistent with the experimentally shown gain-of-function phenotype of the V21D mutant.

Anionic phospholipids affect the rate and extent of flux through the mechanosensitive channel of large conductance MscL

The mechanosensitive channel of large conductance MscL from Escherichia coli has been reconstituted into sealed vesicles, and the effects of lipid structure on the flux of the fluorescent molecule calcein through the open channel have been studied. The channel was opened by reaction of the G22C mutant of MscL with the reagent [2-(triethylammonium)ethyl]methanethiosulfonate (MTSET) which introduces five positive charges within the pore constriction. Flux through the channel was small when the lipid was phosphatidylcholine, but addition of the anionic lipids phosphatidylglycerol, phosphatidic acid, or cardiolipin up to 50 mol % resulted in increases in the amplitudes and rates of release of calcein. Similar effects were seen when either wild-type MscL or the G22C mutant was opened by osmotic pressure difference; rates of release of calcein were very slow in the absence of anionic lipid but increased with increasing concentrations of phosphatidylglycerol to 50 mol %. The observed partial release of trapped calcein following activation of MscL was attributed to the formation of a long-lived subconductance state of MscL following channel opening. Effects of anionic lipid were attributed to an increase in the rate of the transition from closed to fully open state and to a decrease in the rate of the transition from the fully open state to the subconductance state. Higher concentrations of anionic lipid led to a decrease in the rate and amplitude of release of calcein, possibly due to a decreased rate of flux through the open channel. In mixtures with anionic lipids, phosphatidylethanolamine resulted in lower rates and amplitude of release than phosphatidylcholine.

Gating of the mechanosensitive channel protein MscL: the interplay of membrane and protein

The mechanosensitive channel of large conductance (MscL) belongs to a family of transmembrane channel proteins in bacteria and functions as a safety valve that relieves the turgor pressure produced by osmotic downshock. MscL gating can be triggered solely by stretching of the membrane. This work reports an effort to understand this mechanotransduction by means of molecular dynamics (MD) simulation on the MscL of mycobacterium tuberculosis embedded in a palmitoyloleoylphosphatidylethanolamine membrane. Equilibrium MD under zero membrane tension produced a more compact protein structure, as measured by its radii of gyration, compared to the crystal structure, in agreement with previous experimental findings. Even under a large applied tension up to 1000 dyn/cm, the MscL lateral dimension largely remained unchanged after up to 20 ns of simulation. A nonequilibrium MD simulation of 3% membrane expansion showed a significant increase in membrane rigidity upon MscL inclusion, which can contribute to efficient mechanotransduction. Direct observation of channel opening was possible only when an explicit lateral bias force was applied to each of the five subunits of MscL in the radially outward direction. Using this force, open structures with a large pore of radius 10 A could be obtained. The channel opening takes place in a stepwise manner and concurrently with the water chain formation across the channel, which occurs without direct involvement of protein hydrophilic residues. The N-terminal S1 helices stabilize the open structure, and the membrane asymmetry (different lipid density on the two leaflets of membrane) promotes channel opening.

Different effects of lipid chain length on the two sides of a membrane and the lipid annulus of MscL

Quenching of the fluorescence of Trp residues in a membrane protein by lipids with bromine-containing fatty acyl chains provides a powerful technique for measuring lipid-protein binding constants. Single Trp residues have been placed on the periplasmic and cytoplasmic sides of the mechanosensitive channel of large conductance MscL from Mycobacterium tuberculosis to measure, separately, lipid binding constants on the two faces of MscL. The chain-length dependence of lipid binding was found to be different on the two sides of MscL, the chain-length dependence being more marked on the cytoplasmic than on the periplasmic side. To determine if lipid binding constants are affected by the properties of the lipid molecules not in direct contact with MscL (the bulk lipid), the amount of bulk lipid present in the system was varied. The binding constant of the short-chain phospholipid didodecylphosphatidylcholine was found to be independent of the molar ratio of lipid/MscL pentamer over the range 500:1-50:1, suggesting that lipid binding constants are determined largely by the properties of the lipid molecules interacting directly with MscL. These results point to a model in which lipid molecules located on the transmembrane surface of a membrane protein (the annular lipid molecules), by playing a dominant role in the interaction between a membrane protein and the surrounding lipid bilayer, could effectively buffer the membrane protein from changes in the properties of the bulk lipid bilayer.

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.