Molecular dynamics study of MscL interactions with a curved lipid bilayer

Molecular Dynamics Simulations of Curved Lipid Membranes

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.

Elucidating the molecular basis of spontaneous activation in an engineered mechanosensitive channel

Mechanosensitive channel of large conductance (MscL) detects and responds to changes in the pressure profile of cellular membranes and transduces the mechanical energy into electrical and/or chemical signals. MscL can be activated using ultrasonic or chemical activation methods to improve the absorption of medicines and bioactive compounds into cells. However, re-engineering chemical signals such as pH change can trigger channel activation in MscL. This study elucidates the activation mechanism of an engineered MscL at an atomic level through a combination of equilibrium and non-equilibrium (NE) molecular dynamics (MD) simulations. Comparing the wild-type (WT) and engineered MscL activation processes suggests that the two systems are likely associated with different active states and different transition pathways. These findings indicate that (1) periplasmic loops play a key role in the activation process of MscL, (2) the loss of various backbone-backbone hydrogen bonds and salt bridge interactions in the engineered MscL channel causes the spontaneous opening of the channel, and (3) the most significant interactions lost during the activation process are between the transmembrane helices 1 and 2 in engineered MscL channel. The orientation-based biasing approach for producing and optimizing an open MscL model used in this work is a promising way to characterize unknown protein functional states and investigate the activation processes in ion channels and transmembrane proteins in general. This work paves the way for a computational framework for engineering more efficient pH-sensing mechanosensitive channels.

Biophysical Principles of Ion-Channel-Mediated Mechanosensory Transduction

Recent rapid progress in the field of mechanobiology has been driven by novel emerging tools and methodologies and growing interest from different scientific disciplines. Specific progress has been made toward understanding how cell mechanics is linked to intracellular signaling and the regulation of gene expression in response to a variety of mechanical stimuli. There is a direct link between the mechanoreceptors at the cell surface and intracellular biochemical signaling, which in turn controls downstream effector molecules. Among the mechanoreceptors in the cell membrane, mechanosensitive (MS) ion channels are essential for the ultra-rapid (millisecond) transduction of mechanical stimuli into biologically relevant signals. The three decades of research on mechanosensitive channels resulted in the formulation of two basic principles of mechanosensitive channel gating: force-from-lipids and force-from-filament. In this review, we revisit the biophysical principles that underlie the innate force-sensing ability of mechanosensitive channels as contributors to the force-dependent evolution of life forms.

Remembering Mechanosensitivity of NMDA Receptors

An increase in post-synaptic Ca²⁺ conductance through activation of the ionotropic N-methyl-D-aspartate receptor (NMDAR) and concomitant structural changes are essential for the initiation of long-term potentiation (LTP) and memory formation. Memories can be initiated by coincident events, as occurs in classical conditioning, where the NMDAR can act as a molecular coincidence detector. Binding of glutamate and glycine, together with depolarization of the postsynaptic cell membrane to remove the Mg²⁺ channel pore block, results in NMDAR opening for Ca²⁺ conductance. Accumulating evidence has implicated both force-from-lipids and protein tethering mechanisms for mechanosensory transduction in NMDAR, which has been demonstrated by both, membrane stretch and application of amphipathic molecules such as arachidonic acid (AA). The contribution of mechanosensitivity to memory formation and consolidation may be to increase activity of the NMDAR leading to facilitated memory formation. In this review we look back at the progress made toward understanding the physiological and pathological role of NMDA receptor channels in mechanobiology of the nervous system and consider these findings in like of their potential functional implications for memory formation. We examine recent studies identifying mechanisms of both NMDAR and other mechanosensitive channels and discuss functional implications including gain control of NMDA opening probability. Mechanobiology is a rapidly growing area of biology with many important implications for understanding form, function and pathology in the nervous system.

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.

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.

Computational Design of Membrane Curvature-Sensing Peptides

Computer simulations have become an indispensable tool in studying molecular biological systems. The unmatched spatial and temporal resolution that it offers enables for microscopic-level views into the dynamics and mechanics of biological systems. Recent advances in hardware resources have also opened up to computer simulations the investigation of longer timescale biological processes and larger systems. The study of membrane proteins or peptides especially benefits from simulations due to difficulties related to crystallization of such proteins in a membrane environment. In this chapter, we outline the method of molecular dynamics and how it is applied to simulations that involve a peptide and lipid bilayers. In particular, the simulation of a membrane-curvature sensing peptide is examined, and ways of employing computational simulations to design such peptides are discussed.

Nanomechanical properties of MscL α helices: A steered molecular dynamics study

Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a 'lubricant' (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.

Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study

Interaction of fullerenes with asymmetric and curved DOPC/DOPS bicelles is studied by means of coarse-grained molecular dynamics simulations. The effects caused by asymmetric lipid composition of the membrane leaflets and the curvature of the membrane are analyzed. It is shown that the aggregates of fullerenes prefer to penetrate into the membrane in the regions of the moderately positive mean curvature. Upon penetration into the hydrophobic core of the membrane fullerenes avoid the regions of the extreme positive or the negative curvature. Fullerenes increase the ordering of lipid tails, which are in direct contact with them, but do not influence other lipids significantly. Our data suggest that the effects of the membrane curvature should be taken into account in the studies concerning permeability of the membranes to fullerenes and fullerene-based drug delivery systems.

Functional similarities between heterogeneously and homogenously expressed MscL constructs

The mechanosensitive channel of large conductance MscL is a well-characterized mechanically gated non-selective ion channel, which often serves as a prototype mechanosensitive channel for mechanotransduction studies. However, there are some discrepancies between MscL constructs used in these studies, most notably unintended heterogeneous expression from some MscL expression constructs. In this study we investigate the possible cause of this expression pattern, and compare the original non-homogenously expressing constructs with our new homogeneously expressing one to confirm that there is little functional difference between them. In addition, a new MscL construct has been developed with an improved molar extinction coefficient at 280 nm, enabling more accurate protein quantification.

Molecular dynamics study on protein-water interplay in the mechanogating of the bacterial mechanosensitive channel MscL

One of the goals of mechanosensitive channel (MSC) studies is to understand the underlying molecular and biophysical mechanisms of the mechano-gating process from force sensing to gate opening. We focus on the latter process and investigate the role of water in the bacterial MSC MscL, which is activated by membrane tension. We analyze the interplay between water and the gate-constituting amino acids, Leu19-Gly26, through molecular dynamics simulations. To highlight the role of water, specifically hydration of the gate, in MscL gating, we restrain lateral movements of the water molecules along the water-vapor interfaces at the top and bottom of the vapor bubble, plugging the closed gate. The gating behaviors in this model and the normal MscL model, in which water movements are unrestrained, are compared. In the normal model, increased membrane tension breaks the hydrogen bond between Leu19 and Val 23 of the inner helix, exposing the backbone carbonyl oxygen of Leu19 to the water-accessible lumen side of the gate. Associated with this activity, water comes to access the vapor region and stably interacts with the carbonyl oxygen to induce a dewetting to wetting transition that facilitates gate expansion toward channel opening. By contrast, in the water-restrained model, carbonyl oxygen is also exposed, but no further conformational changes occur at the gate. This suggests that gate opening relies on a conformational change initiated by wetting. The penetrated water weakens the hydrophobic interaction between neighboring transmembrane inner helices called the "hydrophobic lock" by wedging into the space between their interacting portions.

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.

Dynamics of protein-protein interactions at the MscL periplasmic-lipid interface

MscL, the highly conserved bacterial mechanosensitive channel of large conductance, is one of the best studied mechanosensors. It is a homopentameric channel that serves as a biological emergency release valve that prevents cell lysis from acute osmotic stress. We previously showed that the periplasmic region of the protein, particularly a single residue located at the TM1/periplasmic loop interface, F47 of Staphylococcus aureus and I49 of Escherichia coli MscL, plays a major role in both the open dwell time and mechanosensitivity of the channel. Here, we introduced cysteine mutations at these sites and found they formed disulfide bridges that decreased the channel open dwell time. By scanning a likely interacting domain, we also found that these sites could be disulfide trapped by addition of cysteine mutations in other locations within the periplasmic loop of MscL, and this also led to rapid channel kinetics. Together, the data suggest structural rearrangements and protein-protein interactions that occur within this region upon normal gating, and further suggest that locking portions of the channel into a transition state decreases the stability of the open state.

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.

The effect of local bending on gating of MscL using a representative volume element and finite element simulation

Many physiological processes such as cell division, endocytosis and exocytosis cause severe local curvature of the cell membrane. Local curvature has been shown experimentally to modulate numerous mechanosensitive (MS) ion channels. In order to quantify the effects of local curvature we introduced a coarse grain representative volume element for the bacterial mechanosensitive ion channel of large conductance (MscL) using continuum elasticity. Our model is designed to be consistent with the channel conformation in the closed and open states to capture its major continuum rheological behavior in response to the local membrane curvature. Herein we show that change in the local curvature of the lipid bilayer can modulate MscL activity considerably by changing both bilayer thickness and lateral pressure profile. Intriguingly, although bending in any direction results in almost the same free-energy cost, inward (cytoplasmic) bending favors channel opening, whereas outward (periplasmic) bending facilitates closing of the narrowest part of the MscL pore. This quantitative study using MscL as a model channel may have wide reaching consequences for the effect of local curvature on the physiological function of other types of prokaryotic and eukaryotic membrane proteins.

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations: from tension sensing to channel opening

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.

Three-dimensional stress field around a membrane protein: atomistic and coarse-grained simulation analysis of gramicidin A

Using both atomistic and coarse-grained (CG) models, we compute the three-dimensional stress field around a gramicidin A (gA) dimer in lipid bilayers that feature different degrees of negative hydrophobic mismatch. The general trends in the computed stress field are similar at the atomistic and CG levels, supporting the use of the CG model for analyzing the mechanical features of protein/lipid/water interfaces. The calculations reveal that the stress field near the protein-lipid interface exhibits a layered structure with both significant repulsive and attractive regions, with the magnitude of the stress reaching 1000 bar in certain regions. Analysis of density profiles and stress field distributions helps highlight the Trp residues at the protein/membrane/water interface as mechanical anchors, suggesting that similar analysis is useful for identifying tension sensors in other membrane proteins, especially membrane proteins involved in mechanosensation. This work fosters a connection between microscopic and continuum mechanics models for proteins in complex environments and makes it possible to test the validity of assumptions commonly made in continuum mechanics models for membrane mediated processes. For example, using the calculated stress field, we estimate the free energy of membrane deformation induced by the hydrophobic mismatch, and the results for regions beyond the annular lipids are in general consistent with relevant experimental data and previous theoretical estimates using elasticity theory. On the other hand, the assumptions of homogeneous material properties for the membrane and a bilayer thickness at the protein/lipid interface being independent of lipid type (e.g., tail length) appear to be oversimplified, highlighting the importance of annular lipids of membrane proteins. Finally, the stress field analysis makes it clear that the effect of even rather severe hydrophobic mismatch propagates to only about two to three lipid layers, thus putting a limit on the range of cooperativity between membrane proteins in crowded cellular membranes.

Changes in lipid density induce membrane curvature

Highly curved bilayer lipid membranes make up the shell of many intra- and extracellular compartments, including organelles and vesicles. Using all-atom molecular dynamics simulations, we show that increasing the density of lipids in the bilayer membrane can induce the membrane to form a curved shape.

Structural investigation of MscL gating using experimental data and coarse grained MD simulations

The mechanosensitive channel of large conductance (MscL) has become a model system in which to understand mechanosensation, a process involved in osmoregulation and many other physiological functions. While a high resolution closed state structure is available, details of the open structure and the gating mechanism remain unknown. In this study we combine coarse grained simulations with restraints from EPR and FRET experiments to study the structural changes involved in gating with much greater level of conformational sampling than has previously been possible. We generated a set of plausible open pore structures that agree well with existing open pore structures and gating models. Most interestingly, we found that membrane thinning induces a kink in the upper part of TM1 that causes an outward motion of the periplasmic loop away from the pore centre. This previously unobserved structural change might present a new mechanism of tension sensing and might be related to a functional role in osmoregulation.

Cells in biological organisms have to be able to respond to mechanical forces during processes such as touch, hearing, pain sensation and tissue growth. One way this is achieved is through mechanosensitive ion channels, membrane embedded proteins that initiate electrical signalling upon tension within the cell or cell membrane. The malfunction of such channels is also associated with a range of diseases including muscular dystrophy and cardiac arrhythmia. In this manuscript, we study in detail the mechanosensitive channel of large conductance (MscL) from bacteria, a model system in which to understand the principles of mechanosensation. Despite many years of investigative work the details of how the protein senses tension in the surrounding membrane remain unknown. By combining structural data from experiments with computer simulation we are able to model the open channel structure of the protein and report previously unobserved structural changes that might present a new mechanism of sensing tension. The methods developed in this paper are not limited to the study of mechanosensitive ion channels and may be useful in understanding the structure and function of other membrane proteins.

Bacterial mechanosensitive channels as a paradigm for mechanosensory transduction

Research on bacterial mechanosensitive (MS) channels has since their discovery been at the forefront of the MS channel field due to extensive studies of the structure and function of MscL and MscS, two of the several different types of MS channels found in bacteria. Just a few years after these two MS channels were cloned their 3D structure was solved by X-ray crystallography. Today, the repertoire of multidisciplinary approaches used in experimental and theoretical studies following the cloning and crystallographic determination of the MscL and MscS structure has expanded by including electronparamagnetic resonance (EPR) and Förster resonance energy transfer (FRET) spectroscopy aided by computational modelling employing molecular dynamics as well as Brownian dynamics simulations, which significantly advanced the understanding of structural determinants of the gating and conduction properties of these two MS channels. These extensive multidisciplinary studies of MscL and MscS have greatly contributed to elucidation of the basic physical principles of MS channel gating by mechanical force. This review summarizes briefly the major experimental and conceptual advancements, which helped in establishing MscL and MscS as a major paradigm of mechanosensory transduction in living cells.

Membrane tension, lipid adaptation, conformational changes, and energetics in MscL gating

This study aims to explore gating mechanisms of mechanosensitive channels in terms of membrane tension, membrane adaptation, protein conformation, and energetics. The large conductance mechanosensitive channel from Mycobacterium tuberculosis (Tb-MscL) is used as a model system; Tb-MscL acts as a safety valve by releasing small osmolytes through the channel opening under extreme hypoosmotic conditions. Based on the assumption that the channel gating involves tilting of the transmembrane (TM) helices, we have performed free energy simulations of Tb-MscL as a function of TM helix tilt angle in a dimyristoylphosphatidylcholine bilayer. Based on the change in system dimensions, TM helix tilting is shown to be essentially equivalent to applying an excess surface tension to the membrane, causing channel expansion, lipid adaptation, and membrane thinning. Such equivalence is further corroborated by the observation that the free energy cost of Tb-MscL channel expansion is comparable to the work done by the excess surface tension. Tb-MscL TM helix tilting results in an expanded water-conducting channel of an outer dimension similar to the proposed fully open MscL structure. The free energy decomposition indicates a possible expansion mechanism in which tilting and expanding of TM2 facilitates the iris-like motion of TM1, producing an expanded Tb-MscL.

Mechanism of membrane curvature sensing by amphipathic helix containing proteins

There are several examples of membrane-associated protein domains that target curved membranes. This behavior is believed to have functional significance in a number of essential pathways, such as clathrin-mediated endocytosis, which involve dramatic membrane remodeling and require the recruitment of various cofactors at different stages of the process. This work is motivated in part by recent experiments that demonstrated that the amphipathic N-terminal helix of endophilin (H0) targets curved membranes by binding to hydrophobic lipid bilayer packing defects which increase in number with increasing membrane curvature. Here we use state-of-the-art atomistic simulation to explore the packing defect structure of curved membranes, and the effect of this structure on the folding of H0. We find that not only are packing defects increased in number with increasing membrane curvature, but also that their size distribution depends nontrivially on the curvature, falling off exponentially with a decay constant that depends on the curvature, and crucially that even on highly curved membranes defects large enough to accommodate the hydrophobic face of H0 are never observed. We furthermore find that a percolation model for the defects explains the defect size distribution, which implies that larger defects are formed by coalescence of noninteracting smaller defects. We also use the recently developed metadynamics algorithm to study in detail the effect of such defects on H0 folding. It is found that the comparatively larger defects found on a convex membrane promote H0 folding by several kcal/mol, while the smaller defects found on flat and concave membrane surfaces inhibit folding by kinetically trapping the peptide. Together, these observations suggest H0 folding is a cooperative process in which the folding peptide changes the defect structure relative to an unperturbed membrane.

Nanomechanics of multiple units in the erythrocyte membrane skeletal network

Erythrocytes undergo deformations when they transport O₂ and CO₂ across the membrane, yet the 3D nanomechanics of the skeletal network remains poorly understood. Expanding from our previous single isolated unit, we now simulate networks consisting of 1–10 concentric rings of repeating units in equibiaxial deformation. The networks are organized with (1) a 3D model for a single unit, (2) a wrap-around mode between Sp and actin protofilament in the intra-unit interaction, and (3) a random inter-unit connectivity. These assumptions permit efficient five-degrees-of-freedom (5DOF) simulations when up to 30 pN of radial forces are applied to the boundary spectrin (Sp) and the center and other units are analyzed. As 6 Sp balance their tensions, hexagonal units become irregular. While actin protofilaments remain tangent to the network, their yaw (Φ) angles change drastically with addition of neighboring units or an Sp unfolding. It is anticipated that during deformation, transmembrane complexes associated with the network move laterally through the lipid bilayer and increase the diffusion of molecules across the membrane. When protofilament/Sp sweeps under the lipid bilayer, they mix up the submembrane concentration gradient. Thus, the nanomechanics of actin protofilaments and Sp may enhance the transport of molecules during erythrocyte deformation.

Curvature generation and pressure profile modulation in membrane by lysolipids: insights from coarse-grained simulations

Although many membrane additives are known to modulate the activities of membrane proteins via perturbing the properties of lipid membrane, the underlying mechanism is often not precisely understood. In this study, we investigate the impact of asymmetrically incorporating single-tailed lysophosphatidylcholine (LPC) into a membrane bilayer using coarse-grained molecular dynamics simulations. Using a simple computational protocol designed to approximately mimic a micropipette setting, we show that asymmetric incorporation of LPC can lead to significant curvature in a bilayer. Detailed analysis of geometrical and mechanical properties (pressure profile) of the resulting mound structure indicates that the degree of pressure profile perturbation is determined not by the local curvature per se but by the packing of lipid headgroups (i.e., area-per-lipid). The findings help provide a concrete basis for understanding the activation mechanism of mechanosensitive channels by asymmetric incorporation of LPC into membrane patches in patch-clamp experiments. The calculated local pressure profiles are valuable to the construction of realistic membrane models for the analysis of mechanosensation in a continuum mechanics framework.

Thermo-sensitive response based on the membrane fluidity adaptation in Paramecium multimicronucleatum

Relationships between the thermo-sensitive response and membrane lipid fluidity were studied using a ciliated protozoan, Paramecium multimicronucleatum. Paramecium elicits a transient membrane depolarization in response to a cooling stimulus (temperature drop). The depolarization amplitude was largest when the cooling stimulus was started from the culture temperature, whilst when cooling started at a temperature more than 5°C higher or lower than the culture temperature, only a small depolarization was induced. Therefore, the cooling-induced response was dependent on the culture temperature and its sensitivity to the cooling stimulus was highest at the culture temperature. Membrane fluidity measurements of living cells using the fluorescent dye 6-lauroyl-2-dimethylaminonaphthalene (laurdan) showed that the fluidity measured at the culture temperature was almost constant irrespective of the temperature at which the cells had been cultured and adapted, which is consistent with homeoviscous adaptation. The constant fluidity at the culture temperature quickly decreased within a few seconds of application of the cooling stimulus, and the decreased fluidity gradually readapted to a constant level at the decreased temperature within 1 h. When the constant fluidity at culture temperature was modified by the addition of procaine or benzyl alcohol, the cooling-induced depolarization was completely abolished. These results suggest the possibility that the adaptation of fluidity to a constant level and its quick decrease below the constant level activate cooling-sensitive channels to elicit the transient depolarization.

Lipids on the move: simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.

Ion channels in microbes

Studies of ion channels have for long been dominated by the animalcentric, if not anthropocentric, view of physiology. The structures and activities of ion channels had, however, evolved long before the appearance of complex multicellular organisms on earth. The diversity of ion channels existing in cellular membranes of prokaryotes is a good example. Although at first it may appear as a paradox that most of what we know about the structure of eukaryotic ion channels is based on the structure of bacterial channels, this should not be surprising given the evolutionary relatedness of all living organisms and suitability of microbial cells for structural studies of biological macromolecules in a laboratory environment. Genome sequences of the human as well as various microbial, plant, and animal organisms unambiguously established the evolutionary links, whereas crystallographic studies of the structures of major types of ion channels published over the last decade clearly demonstrated the advantage of using microbes as experimental organisms. The purpose of this review is not only to provide an account of acquired knowledge on microbial ion channels but also to show that the study of microbes and their ion channels may also hold a key to solving unresolved molecular mysteries in the future.

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.

Gating mechanisms of mechanosensitive channels of large conductance, II: systematic study of conformational transitions

Part II of this study is based on the continuum mechanics-based molecular dynamics-decorated finite element method (MDeFEM) framework established in Part I. In Part II, the gating pathways of Escherichia coli-MscL channels under various basic deformation modes are simulated. Upon equibiaxial tension (which is verified to be the most effective mode for gating), the MDeFEM results agree well with both experiments and all-atom simulations in literature, as well as the analytical continuum models and elastic network models developed in Part I. Different levels of model sophistication and effects of structural motifs are explored in detail, where the importance of mechanical roles of transmembrane helices, cytoplasmic helices, and loops are discussed. The conformation transitions under complex membrane deformations are predicted, including bending, torsion, cooperativity, patch clamp, and indentation. Compared to atom-based molecular dynamics simulations and elastic network models, the MDeFEM framework is unusually well-suited for simulating complex deformations at large length scales. The versatile hierarchical framework can be further applied to simulate the gating transition of other mechanosensitive channels and other biological processes where mechanical perturbation is important.

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.

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.

Mechanosensitive channel of large conductance

Microbial cells constitutively express the Large Conductance Mechanosensitive Channel which opens in response to stretch forces in the lipid bilayer. The channel protein forms a homopentamer with each subunit containing two transmembrane regions and gates via the bilayer mechanism evoked by hydrophobic mismatch and changes in the membrane curvature and/or transbilayer pressure profile. During the stationary phase and during osmotic shock the channel protein is up-regulated to prevent cell lysis. Pharmacological potential of MscL may involve discovery of new age antibiotics to combat multiple drug-resistant bacterial strains.

Molecular dynamics simulations of potassium channels

Despite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.

Lipid stress at play: mechanosensitivity of voltage-gated channels

Membrane stretch modulates the activity of voltage-gated channels (VGCs). These channels are nearly ubiquitous among eukaryotes and they are present, too, in prokaryotes, so the potential ramifications of VGC mechanosensitivity are diverse. In situ traumatic stretch can irreversibly alter VGC activity with lethal results but that is pathology. This chapter discusses the reversible responses of VGCs to stretch, with the general relation of stretch stimuli to other forms of lipid stress, and briefly, with some irreversible stretch effects (=stretch trauma). A working assumption throughout is that mechanosensitive (MS) VGC motions-that is, motions that respond reversibly to bilayer stretch-are susceptible to other forms of lipid stress, such as the stresses produced when amphiphilic molecules (anesthetics, lipids, alcohols, and lipophilic drugs) are inserted into the bilayer. Insofar as these molecules change the bilayer's lateral pressure profile, they can be termed bilayer mechanical reagents (BMRs). The chapter also discusses the MS VGC behavior against the backdrop of eukaryotic channels more widely accepted as "MS channels"--namely, the transient receptor potential (TRP)-based MS cation channels.

Assessing the nature of lipid raft membranes

The paradigm of biological membranes has recently gone through a major update. Instead of being fluid and homogeneous, recent studies suggest that membranes are characterized by transient domains with varying fluidity. In particular, a number of experimental studies have revealed the existence of highly ordered lateral domains rich in sphingomyelin and cholesterol (CHOL). These domains, called functional lipid rafts, have been suggested to take part in a variety of dynamic cellular processes such as membrane trafficking, signal transduction, and regulation of the activity of membrane proteins. However, despite the proposed importance of these domains, their properties, and even the precise nature of the lipid phases, have remained open issues mainly because the associated short time and length scales have posed a major challenge to experiments. In this work, we employ extensive atom-scale simulations to elucidate the properties of ternary raft mixtures with CHOL, palmitoylsphingomyelin (PSM), and palmitoyloleoylphosphatidylcholine. We simulate two bilayers of 1,024 lipids for 100 ns in the liquid-ordered phase and one system of the same size in the liquid-disordered phase. The studies provide evidence that the presence of PSM and CHOL in raft-like membranes leads to strongly packed and rigid bilayers. We also find that the simulated raft bilayers are characterized by nanoscale lateral heterogeneity, though the slow lateral diffusion renders the interpretation of the observed lateral heterogeneity more difficult. The findings reveal aspects of the role of favored (specific) lipid-lipid interactions within rafts and clarify the prominent role of CHOL in altering the properties of the membrane locally in its neighborhood. Also, we show that the presence of PSM and CHOL in rafts leads to intriguing lateral pressure profiles that are distinctly different from corresponding profiles in nonraft-like membranes. The results propose that the functioning of certain classes of membrane proteins is regulated by changes in the lateral pressure profile, which can be altered by a change in lipid content.