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

Mechanosensitive channel MscL gating transitions coupling with constriction point shift

The mechanosensitive channel of large conductance (MscL) acts as an "emergency release valve" that protects bacterial cells from acute hypoosmotic stress, and it serves as a paradigm for studying the mechanism underlying the transduction of mechanical forces. MscL gating is proposed to initiate with an expansion without opening, followed by subsequent pore opening via a number of intermediate substates, and ends in a full opening. However, the details of gating process are still largely unknown. Using in vivo viability assay, single channel patch clamp recording, cysteine cross-linking, and tryptophan fluorescence quenching approach, we identified and characterized MscL mutants with different occupancies of constriction region in the pore domain. The results demonstrated the shifts of constriction point along the gating pathway towards cytoplasic side from residue G26, though G22, to L19 upon gating, indicating the closed-expanded transitions coupling of the expansion of tightly packed hydrophobic constriction region to conduct the initial ion permeation in response to the membrane tension. Furthermore, these transitions were regulated by the hydrophobic and lipidic interaction with the constricting "hot spots". Our data reveal a new resolution of the transitions from the closed to the opening substate of MscL, providing insights into the gating mechanisms of MscL.

Tight hydrophobic core and flexible helices yield MscL with a high tension gating threshold and a membrane area mechanical strain buffer

The mechanosensitive (MS) channel of large conductance, MscL, is the high-tension threshold osmolyte release valve that limits turgor pressure in bacterial cells in the event of drastic hypoosmotic shock. Despite MscL from Mycobacterium tuberculosis (TbMscL) being the first structurally characterized MS channel, its protective mechanism of activation at nearly-lytic tensions has not been fully understood. Here, we describe atomistic simulations of expansion and opening of wild-type (WT) TbMscL in comparison with five of its gain-of-function (GOF) mutants. We show that under far-field membrane tension applied to the edge of the periodic simulation cell, WT TbMscL expands into a funnel-like structure with trans-membrane helices bent by nearly 70°, but does not break its 'hydrophobic seal' within extended 20 μs simulations. GOF mutants carrying hydrophilic substitutions in the hydrophobic gate of increasing severity (A20N, V21A, V21N, V21T and V21D) also quickly transition into funnel-shaped conformations but subsequently fully open within 1-8 μs. This shows that solvation of the de-wetted (vapor-locked) constriction is the rate-limiting step in the gating of TbMscL preceded by area-buffering silent expansion. Pre-solvated gates in these GOF mutants reduce this transition barrier according to hydrophilicity and the most severe V21D eliminates it. We predict that the asymmetric shape-change of the periplasmic side of the channel during the silent expansion provides strain-buffering to the outer leaflet thus re-distributing the tension to the inner leaflet, where the gate resides.

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.

Initiating Coarse-Grained MD Simulations for Membrane-Bound Proteins

Molecular dynamics (MD) simulations have become a widely used tool in the scientific community for understanding molecular scale phenomena that are challenging to address with wet-lab techniques. Coarse-grained simulations, in which multiple atoms are combined into single beads, allow for larger systems and longer time scales to be explored than atomistic techniques. Here, we describe the procedures and equipment required to set up coarse-grained simulations of membrane-bound proteins in a lipid bilayer that can mimic many membrane environments.

Molecular Paradigms for Biological Mechanosensing

Many proteins in living cells are subject to mechanical forces, which can be generated internally by molecular machines, or externally, e.g., by pressure gradients. In general, these forces fall in the piconewton range, which is similar in magnitude to forces experienced by a molecule due to thermal fluctuations. While we would naively expect such moderate forces to produce only minimal changes, a wide variety of "mechanosensing" proteins have evolved with functions that are responsive to forces in this regime. The goal of this article is to provide a physical chemistry perspective on protein-based molecular mechanosensing paradigms used in living systems, and how these paradigms can be explored using novel computational methods.

Mechanical Activation of MscL Revealed by a Locally Distributed Tension Molecular Dynamics Approach

Membrane tension perceived by mechanosensitive (MS) proteins mediates cellular responses to mechanical stimuli and osmotic stresses, and it also guides multiple biological functions including cardiovascular control and development. In bacteria, MS channels function as tension-activated pores limiting excessive turgor pressure, with MS channel of large conductance (MscL) acting as an emergency release valve preventing cell lysis. Previous attempts to simulate gating transitions in MscL by either directly applying steering forces to the protein or by increasing the whole-system tension were not fully successful and often disrupted the integrity of the system. We present a novel, to our knowledge, locally distributed tension molecular dynamics (LDT-MD) simulation method that allows application of forces continuously distributed among lipids surrounding the channel using a specially constructed collective variable. We report reproducible and reversible transitions of MscL to the open state with measured parameters of lateral expansion and conductivity that exactly satisfy experimental values. The LDT-MD method enables exploration of the MscL-gating process with different pulling velocities and variable tension asymmetry between the inner and outer membrane leaflets. We use LDT-MD in combination with well-tempered metadynamics to reconstruct the tension-dependent free-energy landscape for the opening transition in MscL. The flexible definition of the LDT collective variable allows general application of our method to study mechanical activation of any membrane-embedded protein.

Interaction between mechanosensitive channels embedded in lipid membrane

The study of the gating mechanism of mechanosensitive channels opens a window to the exploration of how different mechanical stimuli induce adaptive cellular behaviors of both the protein and the lipid, across different time and length scales. In this work, through a molecular dynamics-decorated finite element method (MDeFEM), the gating behavior of mechanosensitive channels of small conductance (MscS) in Escherichia coli (E. coli) is studied upon membrane stretch or global bending. The local membrane curvature around MscS is incorporated, as well as multiple MscL (mechanosensitive channels of large conductance) molecules in proximity to MscS. The local membrane curvature is found to delay MscS opening and diminishes moderately upon membrane stretching. Mimicking the insertion of lysophosphatidylcholine (LPC) molecules into the lipid, both downward and upward bending can active MscS, as long as the global membrane curvature radius reaches 34 nm. Based on the different MscS pore evolutions observed with the presence of one or more MscLs nearby, we propose that when coreconstituted, multiple MscL molecules tend to be located at the local membrane curvature zone around MscS. In another word, as MscL "swims around" in the lipid bilayer, it can be trapped by the membrane's local curvature. Collectively, the current study provides valuable insights into the interplay between mechanosensitive channels and lipid membrane at structural and physical levels, and specific predictions are proposed for further experimental investigations.

Emerging Diversity in Lipid-Protein Interactions

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

Gating and inactivation of mechanosensitive channels of small conductance: A continuum mechanics study

Mechanosensitive channels of small conductance (MscS) in Escherichia coli (E. coli) serve as a paradigm for understanding the gating behaviors of the MscS family of ion channels. In this work, we develop a continuum mechanics framework to explore the conformational states of MscS during the gating transition. A complete gating transition trajectory from the closed to the open state along with partially open intermediates is obtained, and the open structure is close to the available structural model from crystallographic studies. The computational efficiency of the modeling framework makes it possible to explore the roles of various structural elements (e.g., loops that connect transmembrane helices) and specific interactions in the gating transition. It is observed that removing either the Asp62-Arg131 salt bridge or the Phe68-Leu111 non-polar interaction leads to essentially non-conducting structures even with a membrane tension close to the lysis limit. The loop connecting TM2 (the second transmembrane helix) and TM3 is found to be essential for force transmission during gating, while the loop connecting TM1 and TM2 does not make any major contribution. Based on the different structural evolutions observed when the TM3 kink is treated as a loop or a helical segment, we propose that the helical propensity of the kink plays a central role in inactivation; i.e., under prolonged sub-threshold membrane tension, transition of the initially flexible loop to a helical segment in TM3 may lead to MscS inactivation. Finally, the gating transition of MscS under different transmembrane voltages is explored and found to be essentially voltage independent. Collectively, results from the current continuum mechanics analysis provide further insights into the gating transition of MscS at structural and physical levels, and specific predictions are proposed for further experimental investigations.

Elucidating a role for the cytoplasmic domain in the Mycobacterium tuberculosis mechanosensitive channel of large conductance

Microbial survival in dynamic environments requires the ability to successfully respond to abrupt changes in osmolarity. The mechanosensitive channel of large conductance (MscL) is a ubiquitous channel that facilitates the survival of bacteria and archaea under severe osmotic downshock conditions by relieving excess turgor pressure in response to increased membrane tension. A prominent structural feature of MscL, the cytoplasmic C-terminal domain, has been suggested to influence channel assembly and function. In this report, we describe the X-ray crystal structure and electrophysiological properties of a C-terminal domain truncation of the Mycobacterium tuberculosis MscL (MtMscLΔC). A crystal structure of MtMscLΔC solubilized in the detergent n-dodecyl-β-D-maltopyranoside reveals the pentameric, closed state-like architecture for the membrane spanning region observed in the previously solved full-length MtMscL. Electrophysiological characterization demonstrates that MtMscLΔC retains mechanosensitivity, but with conductance and tension sensitivity more closely resembling full length EcMscL than MtMscL. This study establishes that the C-terminal domain of MtMscL is not required for oligomerization of the full-length channel, but rather influences the tension sensitivity and conductance properties of the channel. The collective picture that emerges from these data is that each MscL channel structure has characteristic features, highlighting the importance of studying multiple homologs.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

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.

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.

Gating mechanism of mechanosensitive channel of large conductance: a coupled continuum mechanical-continuum solvation approach

Gating transition of the mechanosensitive channel of large conductance (MscL) represents a good example of important biological processes that are difficult to describe using atomistic simulations due to the large (submicron) length scale and long (millisecond) time scale. Here we develop a novel computational framework that tightly couples continuum mechanics with continuum solvation models to study the detailed gating behavior of E. coli-MscL. The components of protein molecules are modeled by continuum elements that properly describe their shape, material properties and physicochemical features (e.g., charge distribution). The lipid membrane is modeled as a three-layer material in which the lipid head group and tail regions are treated separately, taking into account the fact that fluidic lipid bilayers do not bear shear stress. Coupling between mechanical and chemical responses of the channel is realized by an iterative integration of continuum mechanics (CM) modeling and continuum solvation (CS) computation. Compared to previous continuum mechanics studies, the present model is capable of capturing the most essential features of the gating process in a much more realistic fashion: due mainly to the apolar solvation contribution, the membrane tension for full opening of MscL is reduced substantially to the experimental measured range. Moreover, the pore size stabilizes constantly during gating because of the intricate interactions of the multiple components of the system, implying the mechanism for sub-conducting states of MscL gating. A significant fraction ([Formula: see text]2/3) of the gating membrane strain is required to reach the first sub-conducting state of our model, which is featured with a relative conductance of 0.115 to the fully opened state. These trends agree well with experimental observations. We anticipate that the coupled CM/CS modeling framework is uniquely suited for the analysis of many biomolecules and their assemblies under external mechanical stimuli.

Multifunctional, Micropipette-based Method for Incorporation And Stimulation of Bacterial Mechanosensitive Ion Channels in Droplet Interface Bilayers

MscL, a large conductance mechanosensitive channel (MSC), is a ubiquitous osmolyte release valve that helps bacteria survive abrupt hypo-osmotic shocks. It has been discovered and rigorously studied using the patch-clamp technique for almost three decades. Its basic role of translating tension applied to the cell membrane into permeability response makes it a strong candidate to function as a mechanoelectrical transducer in artificial membrane-based biomolecular devices. Serving as building blocks to such devices, droplet interface bilayers (DIBs) can be used as a new platform for the incorporation and stimulation of MscL channels. Here, we describe a micropipette-based method to form DIBs and measure the activity of the incorporated MscL channels. This method consists of lipid-encased aqueous droplets anchored to the tips of two opposing (coaxially positioned) borosilicate glass micropipettes. When droplets are brought into contact, a lipid bilayer interface is formed. This technique offers control over the chemical composition and the size of each droplet, as well as the dimensions of the bilayer interface. Having one of the micropipettes attached to a harmonic piezoelectric actuator provides the ability to deliver a desired oscillatory stimulus. Through analysis of the shapes of the droplets during deformation, the tension created at the interface can be estimated. Using this technique, the first activity of MscL channels in a DIB system is reported. Besides MS channels, activities of other types of channels can be studied using this method, proving the multi-functionality of this platform. The method presented here enables the measurement of fundamental membrane properties, provides a greater control over the formation of symmetric and asymmetric membranes, and is an alternative way to stimulate and study mechanosensitive channels.

Activation of bacterial channel MscL in mechanically stimulated droplet interface bilayers

MscL, a stretch-activated channel, saves bacteria experiencing hypo-osmotic shocks from lysis. Its high conductance and controllable activation makes it a strong candidate to serve as a transducer in stimuli-responsive biomolecular materials. Droplet interface bilayers (DIBs), flexible insulating scaffolds for such materials, can be used as a new platform for incorporation and activation of MscL. Here, we report the first reconstitution and activation of the low-threshold V23T mutant of MscL in a DIB as a response to axial compressions of the droplets. Gating occurs near maximum compression of both droplets where tension in the membrane is maximal. The observed 0.1-3 nS conductance levels correspond to the V23T-MscL sub-conductive and fully open states recorded in native bacterial membranes or liposomes. Geometrical analysis of droplets during compression indicates that both contact angle and total area of the water-oil interfaces contribute to the generation of tension in the bilayer. The measured expansion of the interfaces by 2.5% is predicted to generate a 4-6 mN/m tension in the bilayer, just sufficient for gating. This work clarifies the principles of interconversion between bulk and surface forces in the DIB, facilitates the measurements of fundamental membrane properties, and improves our understanding of MscL response to membrane tension.

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.

Coarse-grained force fields for molecular simulations

Molecular dynamics (MD) simulations at the atomic scale are a powerful tool to study the structure and dynamics of model biological systems. However, because of their high computational cost, the time and length scales of atomistic simulations are limited. Biologically important processes, such as protein folding, ion channel gating, signal transduction, and membrane remodeling, are difficult to investigate using atomistic simulations. Coarse-graining reduces the computational cost of calculations by reducing the number of degrees of freedom in the model, allowing simulations of larger systems for longer times. In the first part of this chapter we review briefly some of the coarse-grained models available for proteins, focusing on the specific scope of each model. Then we describe in more detail the MARTINI coarse-grained force field, and we illustrate how to set up and run a simulation of a membrane protein using the Gromacs software package. We explain step-by-step the preparation of the protein and the membrane, the insertion of the protein in the membrane, the equilibration of the system, the simulation itself, and the analysis of the trajectory.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) due to difficulties associated with expression and crystallization. As such, it is one area where computational studies, particularly Molecular Dynamics (MD) simulations, can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, though not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

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.

Coarse-grained molecular dynamics simulations of protein-ligand binding

Coarse-grained molecular dynamics (CGMD) simulations with the MARTINI force field were performed to reproduce the protein-ligand binding processes. We chose two protein-ligand systems, the levansucrase-sugar (glucose or sucrose), and LinB-1,2-dichloroethane systems, as target systems that differ in terms of the size and shape of the ligand-binding pocket and the physicochemical properties of the pocket and the ligand. Spatial distributions of the Coarse-grained (CG) ligand molecules revealed potential ligand-binding sites on the protein surfaces other than the real ligand-binding sites. The ligands bound most strongly to the real ligand-binding sites. The binding and unbinding rate constants obtained from the CGMD simulation of the levansucrase-sucrose system were approximately 10 times greater than the experimental values; this is mainly due to faster diffusion of the CG ligand in the CG water model. We could obtain dissociation constants close to the experimental values for both systems. Analysis of the ligand fluxes demonstrated that the CG ligand molecules entered the ligand-binding pockets through specific pathways. The ligands tended to move through grooves on the protein surface. Thus, the CGMD simulations produced reasonable results for the two different systems overall and are useful for studying the protein-ligand binding processes.

Membrane remodeling capacity of a vesicle-inducing glycosyltransferase

Intracellular vesicles are abundant in eukaryotic cells but absent in the Gram-negative bacterium Escherichia coli. However, strong overexpression of a monotopic glycolipid-synthesizing enzyme, monoglucosyldiacylglycerol synthase from Acholeplasma laidlawii (alMGS), leads to massive formation of vesicles in the cytoplasm of E. coli. More importantly, alMGS provides a model system for the regulation of membrane properties by membrane-bound enzymes, which is critical for maintaining cellular integrity. Both phenomena depend on how alMGS binds to cell membranes, which is not well understood. Here, we carry out a comprehensive investigation of the membrane binding of alMGS by combining bioinformatics methods with extensive biochemical studies, structural modeling and molecular dynamics simulations. We find that alMGS binds to the membrane in a fairly upright manner, mainly by residues in the N-terminal domain, and in a way that induces local enrichment of anionic lipids and a local curvature deformation. Furthermore, several alMGS variants resulting from substitution of residues in the membrane anchoring segment are still able to generate vesicles, regardless of enzymatic activity. These results clarify earlier theories about the driving forces for vesicle formation, and shed new light on the membrane binding properties and enzymatic mechanism of alMGS and related monotopic GT-B fold glycosyltransferases.

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.

Single molecule FRET reveals pore size and opening mechanism of a mechano-sensitive ion channel

The mechanosensitive channel of large conductance, which serves as a model system for mechanosensitive channels, has previously been crystallized in the closed form, but not in the open form. Ensemble measurements and electrophysiological sieving experiments show that the open-diameter of the channel pore is >25 Å, but the exact size and whether the conformational change follows a helix-tilt or barrel-stave model are unclear. Here we report measurements of the distance changes on liposome-reconstituted MscL transmembrane α-helices, using a ‘virtual sorting’ single-molecule fluorescence energy transfer. We observed directly that the channel opens via the helix-tilt model and the open pore reaches 2.8 nm in diameter. In addition, based on the measurements, we developed a molecular dynamics model of the channel structure in the open state which confirms our direct observations.

DOI:http://dx.doi.org/10.7554/eLife.01834.001

Bacterial cells are full of fluid, and they will burst if they are not able to respond to a build up of pressure. Fortunately, the membrane of a bacterial cell contains channels that can detect the increased mechanical stress on the cell membrane and then open to relieve the pressure.

In many bacterial cells, the last defence against the cell exploding is called the mechanosensitive channel of large conductance (MscL). This is made of five proteins, each of which consists of TM1 and TM2 helixes, which are responsible for opening and closing the channel. Two models have been proposed to explain how the channels are opened. In the barrel-stave model, the TM1 helix moves, while the TM2 helix remains stationary. This results in an open pore that is lined with TM1 and TM2 helixes in the same way that wooden staves line a barrel. In the helix-tilt model, both helixes tilt towards the membrane to open the channel.

Wang et al. have now used a technique called single-molecule fluorescence resonance energy transfer (FRET) to explore the structure of the open channel in E. coli in order to determine which model is correct. In this technique an individual channel is labeled with two different fluorescent molecules. By illuminating the channel with light of a wavelength that excites the first fluorescent molecule, and measuring the strength of the fluorescence from the second molecule, it is possible to work out the distance between the two molecules. From this, the structure of the channel and how it opens and closes can be explored.

Previous attempts to measure the diameters of open channels using fluorescence techniques have suffered from issues caused by the use of large numbers of fluorescent molecules. This has made it necessary to use computational modeling to extract the required data. By looking at a series of individual proteins, Wang et al. overcame these problems and found that the diameter of the fully open pore is 2.8 nm. The result provides strong support for the helix-tilt model.

DOI:http://dx.doi.org/10.7554/eLife.01834.002

Chimeras reveal a single lipid-interface residue that controls MscL channel kinetics as well as mechanosensitivity

MscL, the highly conserved bacterial mechanosensitive channel of large conductance, serves as an osmotic "emergency release valve," is among the best-studied mechanosensors, and is a paradigm of how a channel senses and responds to membrane tension. Although all homologs tested thus far encode channel activity, many show functional differences. We tested Escherichia coli and Staphylococcus aureus chimeras and found that the periplasmic region of the protein, particularly E. coli I49 and the equivalent S. aureus F47 at the periplasmic lipid-aqueous interface of the first transmembrane domain, drastically influences both the open dwell time and the threshold of channel opening. One mutant shows a severe hysteresis, confirming the importance of this residue in determining the energy barriers for channel gating. We propose that this site acts similarly to a spring for a clasp knife, adjusting the resistance for obtaining and stabilizing an open or closed channel structure.