Gating of MscL studied by steered molecular dynamics

pSPICA Force Field Extended for Proteins and Peptides

Many coarse-grained (CG) molecular dynamics (MD) studies have been performed to investigate biological processes involving proteins and lipids. CG force fields (FFs) in these MD studies often use implicit or nonpolar water models to reduce computational costs. CG-MD using water models cannot properly describe electrostatic screening effects owing to the hydration of ionic segments and thus cannot appropriately describe molecular events involving water channels and pores through lipid membranes. To overcome this issue, we developed a protein model in the pSPICA FF, in which a polar CG water model showing the proper dielectric response was adopted. The developed CG model greatly improved the transfer free energy profiles of charged side chain analogues across the lipid membrane. Application studies on melittin-induced membrane pores and mechanosensitive channels in lipid membranes demonstrated that CG-MDs using the pSPICA FF correctly reproduced the structure and stability of the pores and channels. Furthermore, the adsorption behavior of the highly charged nona-arginine peptides on lipid membranes changed with salt concentration, indicating the pSPICA FF is also useful for simulating protein adsorption on membrane surfaces.

Challenges of Protein-Protein Docking of the Membrane Proteins

Despite the recent advances in the determination of high-resolution membrane protein (MP) structures, the structural and functional characterization of MPs remains extremely challenging, mainly due to the hydrophobic nature, low abundance, poor expression, purification, and crystallization difficulties associated with MPs. Whereby the major challenges/hurdles for MP structure determination are associated with the expression, purification, and crystallization procedures. Although there have been significant advances in the experimental determination of MP structures, only a limited number of MP structures (approximately less than 1% of all) are available in the Protein Data Bank (PDB). Therefore, the structures of a large number of MPs still remain unresolved, which leads to the availability of widely unplumbed structural and functional information related to MPs. As a result, recent developments in the drug discovery realm and the significant biological contemplation have led to the development of several novel, low-cost, and time-efficient computational methods that overcome the limitations of experimental approaches, supplement experiments, and provide alternatives for the characterization of MPs. Whereby the fine tuning and optimizations of these computational approaches remains an ongoing endeavor.Computational methods offer a potential way for the elucidation of structural features and the augmentation of currently available MP information. However, the use of computational modeling can be extremely challenging for MPs mainly due to insufficient knowledge of (or gaps in) atomic structures of MPs. Despite the availability of numerous in silico methods for 3D structure determination the applicability of these methods to MPs remains relatively low since all methods are not well-suited or adequate for MPs. However, sophisticated methods for MP structure predictions are constantly being developed and updated to integrate the modifications required for MPs. Currently, different computational methods for (1) MP structure prediction, (2) stability analysis of MPs through molecular dynamics simulations, (3) modeling of MP complexes through docking, (4) prediction of interactions between MPs, and (5) MP interactions with its soluble partner are extensively used. Towards this end, MP docking is widely used. It is notable that the MP docking methods yet few in number might show greater potential in terms of filling the knowledge gap. In this chapter, MP docking methods and associated challenges have been reviewed to improve the applicability, accuracy, and the ability to model macromolecular complexes.

A novel force transduction pathway from a tension sensor to the gate in the mechano-gating of MscL channel

The bacterial mechanosensitive channel of large conductance MscL is activated exclusively by increased tension in the membrane bilayer. Despite many proposed models for MscL opening, its precise mechano-gating mechanism, particularly how the received force at the tension sensor transmits to the gate remains incomplete. Previous studies have shown that along with amphipathic N-terminus located near the cytoplasmic surface of the membrane, Phe78 residue near the outer surface also acts as a "tension sensor," while Gly22 is a central constituent of the "hydrophobic gate." Present study focused on elucidating the force transmission mechanism from the sensor Phe78 in the outer transmembrane helix (TM2) to the gate in the inner transmembrane helix (TM1) of MscL by applying the patch clamp and molecular dynamics (MD) simulations to the wild type MscL channel and its single mutants at the sensor (F78N), the gate (G22N) and their combination (G22N/F78N) double mutant. F78N MscL resulted in a severe loss-of-function, while G22N MscL caused a gain-of-function channel exhibiting spontaneous openings at the resting membrane tension. We initially speculated that the spontaneous opening in G22N mutant might occur without tension acting on Phe78 residue. To test this hypothesis, we examined the (G22N/F78N) double mutant, which unexpectedly exhibited neither spontaneous activity nor activity by a relatively high membrane tension. To understand the underlying mechanism, we conducted MD simulations and analyzed the force transduction pathway. Results showed that the mutation at the tension sensor (F78N) in TM2 caused decreased interaction of this residue not only with lipids, but also with a group of amino acids (Ile32-Leu36-Ile40) in the neighboring TM1 helix, which resulted in an inefficient force transmission to the gate-constituting amino acids on TM1. This change also induced a slight tilting of TM1 towards the membrane plane and decreased the size of the channel pore at the gate, which seems to be the major mechanism for the inhibition of spontaneous opening of the double mutant channel. More importantly, the newly identified interaction between the TM2 (Phe78) and adjacent TM1 (Ile32-Leu36-Ile40) helices seems to be an essential force transmitting mechanism for the stretch-dependent activation of MscL given that substitution of any one of these four amino acids with Asn resulted in severe loss-of-function MscL as reported in our previous work.

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.

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

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

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

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.

Pocket delipidation induced by membrane tension or modification leads to a structurally analogous mechanosensitive channel state

The mechanosensitive ion channel of large conductance MscL gates in response to membrane tension changes. Lipid removal from transmembrane pockets leads to a concerted structural and functional MscL response, but it remains unknown whether there is a correlation between the tension-mediated state and the state derived by pocket delipidation in the absence of tension. Here, we combined pulsed electron paramagnetic resonance spectroscopy and hydrogen-deuterium exchange mass spectrometry, coupled with molecular dynamics simulations under membrane tension, to investigate the structural changes associated with the distinctively derived states. Whether it is tension- or modification-mediated pocket delipidation, we find that MscL samples a similar expanded subconducting state. This is the final step of the delipidation pathway, but only an intermediate stop on the tension-mediated path, with additional tension triggering further channel opening. Our findings hint at synergistic modes of regulation by lipid molecules in membrane tension-activated mechanosensitive channels.

Dependence of Protein Structure on Environment: FOD Model Applied to Membrane Proteins

The natural environment of proteins is the polar aquatic environment and the hydrophobic (amphipathic) environment of the membrane. The fuzzy oil drop model (FOD) used to characterize water-soluble proteins, as well as its modified version FOD-M, enables a mathematical description of the presence and influence of diverse environments on protein structure. The present work characterized the structures of membrane proteins, including those that act as channels, and a water-soluble protein for contrast. The purpose of the analysis was to verify the possibility that an external force field can be used in the simulation of the protein-folding process, taking into account the diverse nature of the environment that guarantees a structure showing biological activity.

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.

Twisting and tilting of a mechanosensitive molecular probe detects order in membranes

Lateral forces in biological membranes affect a variety of dynamic cellular processes. Recent synthetic efforts have introduced fluorescent "flippers" as environment-sensitive planarizable push-pull probes that can detect lipid packing and membrane tension, and respond to lipid-induced mechanical forces by a shift in their spectroscopic properties. Herein, we investigate the molecular origin of the mechanosensitivity of the best known flipper, Flipper-TR, by an extended set of molecular dynamics (MD) simulations in membranes of increasing complexity and under different physicochemical conditions, revealing unprecedented details of the sensing process. Simulations enabled by accurate refinement of Flipper-TR force field using quantum mechanical calculations allowed us to unambiguously correlate the planarization of the two fluorescent flippers to spectroscopic response. In particular, Flipper-TR conformation exhibits bimodal distribution in disordered membranes and a unimodal distribution in highly ordered membranes. Such dramatic change was associated with a shift in Flipper-TR excitation spectra, as supported both by our simulated and experimentally-measured spectra. Flipper-TR sensitivity to phase-transition is confirmed by a temperature-jump protocol that alters the lipid phase of an ordered membrane, triggering an instantaneous mechanical twisting of the probe. Simulations show that the probe is also sensitive to surface tension, since even in a naturally disordered membrane, the unimodal distribution of coplanar flippers can be achieved if a sufficiently negative surface tension is applied to the membrane. MD simulations in ternary mixtures containing raft-like nanodomains show that the probe can discriminate lipid domains in phase-separated complex bilayers. A histogram-based approach, called DOB-phase classification, is introduced that can differentiate regions of disordered and ordered lipid phases by comparing dihedral distributions of Flipper-TR. Moreover, a new sensing mechanism involving the orientation of Flipper-TR is elucidated, corroborating experimental evidence that the probe tilt angle is strongly dependent on lipid ordering. The obtained atomic-resolution description of Flipper-TR mechanosensitivity is key to the interpretation of experimental data and to the design of novel mechanosensors with improved spectroscopic properties.

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.

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

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

Pulling MscL open via N-terminal and TM1 helices: A computational study towards engineering an MscL nanovalve

There are great opportunities in the manipulation of bacterial mechanosensitive (MS) ion channels for specific and targeted drug delivery purposes. Recent research has shown that these ion channels have the potential to be converted into nanovalves through clever use of magnetic nanoparticles and magnetic fields. Using a combination of molecular dynamics (MD) simulations and the finite element (FE) modelling, this study investigates the theoretical feasibility of opening the MscL channel (MS channel of large conductance of E. coli) by applying mechanical force directly to its N-terminus. This region has already been reported to function as a major mechanosensor in this channel. The stress-strain behaviour of each MscL helix was obtained using all atom MD simulations. Using the same method, we simulated two models, the wild-type (WT) MscL and the G22N mutant MscL, both embedded in a POPE lipid bilayer. In addition to indicating the main interacting residues at the hydrophobic pore, their pairwise interaction energies were monitored during the channel gating. We implemented these inputs into our FE model of MscL using curve-fitting codes and continuum mechanics equations. In the FE model, the channel could be fully opened via pulling directly on the N-terminus and bottom of TM1 by mutating dominant van der Waals interactions in the channel pore; otherwise the stress generated on the channel protein can irreversibly unravel the N-secondary structure. This is a significant finding suggesting that applying force in this manner is sufficient to open an MscL nanovalve delivering various drugs used, for example, in cancer chemotherapy. More importantly, the FE model indicates that to fully operate an MscL nanovalve by pulling directly on the N-terminus and bottom of TM1, gain-of-function (GOF) mutants (e.g., G22N MscL) would have to be employed rather than the WT MscL channel.

Ion-pulling simulations provide insights into the mechanisms of channel opening of the skeletal muscle ryanodine receptor

The type 1 ryanodine receptor (RyR1) mediates Ca²⁺ release from the sarcoplasmic reticulum to initiate skeletal muscle contraction and is associated with muscle diseases, malignant hyperthermia, and central core disease. To better understand RyR1 channel function, we investigated the molecular mechanisms of channel gating and ion permeation. An adequate model of channel gating requires accurate, high-resolution models of both open and closed states of the channel. To this end, we generated an open-channel RyR1 model using molecular simulations to pull Ca²⁺ through the pore constriction site of a closed-channel RyR1 structure determined at 3.8-Å resolution. Importantly, we find that our open-channel model is consistent with the RyR1 and cardiac RyR (RyR2) open-channel structures reported while this paper was in preparation. Both our model and the published structures show similar rotation of the upper portion of the pore-lining S6 helix away from the 4-fold channel axis and twisting of Ile-4937 at the channel constriction site out of the channel pore. These motions result in a minimum open-channel pore radius of ∼3 Å formed by Gln-4933, rather than Ile-4937 in the closed-channel structure. We also present functional support for our model by mutations around the closed- and open-channel constriction sites (Gln-4933 and Ile-4937). Our results indicate that use of ion-pulling simulations produces a RyR1 open-channel model, which can provide insights into the mechanisms of channel opening complementing those from the structural data.

High-Throughput Simulations Reveal Membrane-Mediated Effects of Alcohols on MscL Gating

The mechanosensitive channels of large conductance (MscL) are bacterial membrane proteins that serve as last resort emergency release valves in case of severe osmotic downshock. Sensing bilayer tension, MscL channels are sensitive to changes in the bilayer environment and are, therefore, an ideal test case for exploring membrane protein coupling. Here, we use high-throughput coarse-grained molecular dynamics simulations to characterize MscL gating kinetics in different bilayer environments under the influence of alcohols. We performed over five hundred simulations to obtain sufficient statistics to reveal the subtle effects of changes in the membrane environment on MscL gating. MscL opening times were found to increase with the addition of the straight-chain alcohols ethanol, octanol, and to some extent dodecanol but not with hexadecanol. Increasing concentration of octanol increased the impeding effect, but only up to 10-20 mol %. Our in silico predictions were experimentally confirmed using reconstituted MscL in a liposomal fluorescent efflux assay. Our combined data reveal that the effect of alcohols on MscL gating arises not through specific binding sites but through a combination of the alcohol-induced changes to a number of bilayer properties and their alteration of the MscL-bilayer interface. Our work provides a key example of how extensive molecular simulations can be used to predict the functional modification of membrane proteins by subtle changes in their bilayer environment.

Analytical calculation of the lipid bilayer bending modulus

Bending and Gaussian moduli of a homogenious single-component lipid bilayer are calculated analytically using microscopic model of the lipid hydrocarbon chains. The approach allows for thermodynamic averaging over different chains conformations. Each chain is modeled as a flexible string with finite bending rigidity and an incompressible cross-section area. The interchain steric repulsion is accounted for self-consistently determined single-chain confining parabolic potential. The model provides a simple analytical expression for the membrane bending modulus, which falls within a range of experimental values. An observed dependence of the modulus on hydrocarbon chain length is also reproduced. Correspondence between our microscopic model and the membrane theory of elasticity is established.

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.

The role of MscL amphipathic N terminus indicates a blueprint for bilayer-mediated gating of mechanosensitive channels

The bacterial mechanosensitive channel MscL gates in response to membrane tension as a result of mechanical force transmitted directly to the channel from the lipid bilayer. MscL represents an excellent model system to study the basic biophysical principles of mechanosensory transduction. However, understanding of the essential structural components that transduce bilayer tension into channel gating remains incomplete. Here using multiple experimental and computational approaches, we demonstrate that the amphipathic N-terminal helix of MscL acts as a crucial structural element during tension-induced gating, both stabilizing the closed state and coupling the channel to the membrane. We propose that this may also represent a common principle in the gating cycle of unrelated mechanosensitive ion channels, allowing the coupling of channel conformation to membrane dynamics.

Revealing Atomic-Level Mechanisms of Protein Allostery with Molecular Dynamics Simulations

Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.

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.

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.

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.

Electromechanical Signatures for DNA Sequencing through a Mechanosensitive Nanopore

Biological nanopores have been extensively used for DNA base detection since these pores are widely available and tunable through mutations. Distinguishing bases of nucleic acids by passing them through nanopores has so far primarily relied on electrical signals-specifically, ionic currents through the nanopores. However, the low signal-to-noise ratio makes detection of ionic currents difficult. In this study, we show that the initially closed mechanosensitive channel of large conductance (MscL) protein pore opens for single-stranded DNA (ssDNA) translocation under an applied electric field. As each nucleotide translocates through the pore, a unique mechanical signal is observed-specifically, the tension in the membrane containing the MscL pore is different for each nucleotide. In addition to the membrane tension, we found that the ionic current is also different for the four nucleotide types. The initially closed MscL adapts its opening for nucleotide translocation due to the flexibility of the pore. This unique operation of MscL provides single nucleotide resolution in both electrical and mechanical signals. Finally, we also show that the speed of DNA translocation is roughly 1 order of magnitude slower in MscL compared to Mycobacterium smegmatis porin A (MspA), suggesting MscL to be an attractive protein pore for DNA sequencing.

Force transduction and lipid binding in MscL: a continuum-molecular approach

The bacterial mechanosensitive channel MscL, a small protein mainly activated by membrane tension, is a central model system to study the transduction of mechanical stimuli into chemical signals. Mutagenic studies suggest that MscL gating strongly depends on both intra-protein and interfacial lipid-protein interactions. However, there is a gap between this detailed chemical information and current mechanical models of MscL gating. Here, we investigate the MscL bilayer-protein interface through molecular dynamics simulations, and take a combined continuum-molecular approach to connect chemistry and mechanics. We quantify the effect of membrane tension on the forces acting on the surface of the channel, and identify interactions that may be critical in the force transduction between the membrane and MscL. We find that the local stress distribution on the protein surface is largely asymmetric, particularly under tension, with the cytoplasmic side showing significantly larger and more localized forces, which pull the protein radially outward. The molecular interactions that mediate this behavior arise from hydrogen bonds between the electronegative oxygens in the lipid headgroup and a cluster of positively charged lysine residues on the amphipathic S1 domain and the C-terminal end of the second trans-membrane helix. We take advantage of this strong interaction (estimated to be 10–13 kT per lipid) to actuate the channel (by applying forces on protein-bound lipids) and explore its sensitivity to the pulling magnitude and direction. We conclude by highlighting the simple motif that confers MscL with strong anchoring to the bilayer, and its presence in various integral membrane proteins including the human mechanosensitive channel K2P1 and bovine rhodopsin.

Global structural changes of an ion channel during its gating are followed by ion mobility mass spectrometry

Mechanosensitive ion channels are sensors probing membrane tension in all species; despite their importance and vital role in many cell functions, their gating mechanism remains to be elucidated. Here, we determined the conditions for releasing intact mechanosensitive channel of large conductance (MscL) proteins from their detergents in the gas phase using native ion mobility-mass spectrometry (IM-MS). By using IM-MS, we could detect the native mass of MscL from Escherichia coli, determine various global structural changes during its gating by measuring the rotationally averaged collision cross-sections, and show that it can function in the absence of a lipid bilayer. We could detect global conformational changes during MscL gating as small as 3%. Our findings will allow studying native structure of many other membrane proteins.

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 mechanoelectrical response of the cytoplasmic membrane of Vibrio cholerae

Persistence of Vibrio cholerae in waters of fluctuating salinity relies on the capacity of this facultative enteric pathogen to adapt to varying osmotic conditions. In an event of osmotic downshift, osmolytes accumulated inside the bacterium can be quickly released through tension-activated channels. With the newly established procedure of giant spheroplast preparation from V. cholerae, we performed the first patch-clamp characterization of its cytoplasmic membrane and compared tension-activated currents with those in Esherichia coli. Saturating pressure ramps revealed two waves of activation belonging to the ∼1-nS mechanosensitive channel of small conductance (MscS)-like channels and ∼3-nS mechanosensitive channel of large conductance (MscL)-like channels, with a pressure midpoint ratio p_0.5MscS/p_0.5MscL of 0.48. We found that MscL-like channels in V. cholerae present at a density three times higher than in E. coli, and yet, these vibrios were less tolerant to large osmotic downshocks. The Vibrio MscS-like channels exhibit characteristic inward rectification and subconductive states at depolarizing voltages; they also adapt and inactivate at subsaturating tensions and recover within 2 s upon tension release, just like E. coli MscS. Trehalose, a compatible internal osmolyte accumulated under hypertonic conditions, significantly shifts activation curves of both MscL- and MscS-like channels toward higher tensions, yet does not freely partition into the channel pore. Direct electrophysiology of V. cholerae offers new avenues for the in situ analysis of membrane components critical for osmotic survival and electrogenic transport in this pathogen.

Mechanosensitive channels are activated by stress in the actin stress fibres, and could be involved in gravity sensing in plants

Mechanosensitive (MS) channels are expressed in a variety of cells. The molecular and biophysical mechanism involved in the regulation of MS channel activities is a central interest in basic biology. MS channels are thought to play crucial roles in gravity sensing in plant cells. To date, two mechanisms have been proposed for MS channel activation. One is that tension development in the lipid bilayer directly activates MS channels. The second mechanism proposes that the cytoskeleton is involved in the channel activation, because MS channel activities are modulated by pharmacological treatments that affect the cytoskeleton. We tested whether tension in the cytoskeleton activates MS channels. Mammalian endothelial cells were microinjected with phalloidin-conjugated beads, which bound to stress fibres, and a traction force to the actin cytoskeleton was applied by dragging the beads with optical tweezers. MS channels were activated when the force was applied, demonstrating that a sub-pN force to the actin filaments activates a single MS channel. Plants may use a similar molecular mechanism in gravity sensing, since the cytoplasmic Ca(2+) concentration increase induced by changes in the gravity vector was attenuated by potential MS channel inhibitors, and by actin-disrupting drugs. These results support the idea that the tension increase in actin filaments by gravity-dependent sedimentation of amyloplasts activates MS Ca(2+) -permeable channels, which can be the molecular mechanism of a Ca(2+) concentration increase through gravistimulation. We review recent progress in the study of tension sensing by actin filaments and MS channels using advanced biophysical methods, and discuss their possible roles in gravisensing.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

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.

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.

Voltage-dependent conformational changes in connexin channels

Channels formed by connexins display two distinct types of voltage-dependent gating, termed V(j)- or fast-gating and loop- or slow-gating. Recent studies, using metal bridge formation and chemical cross-linking have identified a region within the channel pore that contributes to the formation of the loop-gate permeability barrier. The conformational changes are remarkably large, reducing the channel pore diameter from 15 to 20Å to less than 4Å. Surprisingly, the largest conformational change occurs in the most stable region of the channel pore, the 3(10) or parahelix formed by amino acids in the 42-51 segment. The data provide a set of positional constraints that can be used to model the structure of the loop-gate closed state. Less is known about the conformation of the V(j)-gate closed state. There appear to be two different mechanisms; one in which conformational changes in channel structure are linked to a voltage sensor contained in the N-terminus of Cx26 and Cx32 and a second in which the C-terminus of Cx43 and Cx40 may act either as a gating particle to block the channel pore or alternatively to stabilize the closed state. The later mechanism utilizes the same domains as implicated in effecting pH gating of Cx43 channels. It is unclear if the two V(j)-gating mechanisms are related or if they represent different gating mechanisms that operate separately in different subsets of connexin channels. A model of the V(j)-closed state of Cx26 hemichannel that is based on the X-ray structure of Cx26 and electron crystallographic structures of a Cx26 mutation suggests that the permeability barrier for V(j)-gating is formed exclusively by the N-terminus, but recent information suggests that this conformation may not represent a voltage-closed state. Closed state models are considered from a thermodynamic perspective based on information from the 3.5Å Cx26 crystal structure and molecular dynamics (MD) simulations. The applications of computational and experimental methods to define the path of allosteric molecular transitions that link the open and closed states are discussed. This article is part of a Special Issue entitled: The Communicating junctions, composition, structure and characteristics.

Mechanosensitive channels: what can they do and how do they do it?

While mechanobiological processes employ diverse mechanisms, at their heart are force-induced perturbations in the structure and dynamics of molecules capable of triggering subsequent events. Among the best characterized force-sensing systems are bacterial mechanosensitive channels. These channels reflect an intimate coupling of protein conformation with the mechanics of the surrounding membrane; the membrane serves as an adaptable sensor that responds to an input of applied force and converts it into an output signal, interpreted for the cell by mechanosensitive channels. The cell can exploit this information in a number of ways: ensuring cellular viability in the presence of osmotic stress and perhaps also serving as a signal transducer for membrane tension or other functions. This review focuses on the bacterial mechanosensitive channels of large (MscL) and small (MscS) conductance and their eukaryotic homologs, with an emphasis on the outstanding issues surrounding the function and mechanism of this fascinating class of molecules.

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.

Elucidation of the glucose transport pathway in glucose transporter 4 via steered molecular dynamics simulations

BACKGROUND

GLUT4 is a predominant insulin regulated glucose transporter expressed in major glucose disposal tissues such as adipocytes and muscles. Under the unstimulated state, GLUT4 resides within intracellular vesicles. Various stimuli such as insulin translocate this protein to the plasma membrane for glucose transport. In the absence of a crystal structure for GLUT4, very little is known about the mechanism of glucose transport by this protein. Earlier we proposed a homology model for GLUT4 and performed a conventional molecular dynamics study revealing the conformational rearrangements during glucose and ATP binding. However, this study could not explain the transport of glucose through the permeation tunnel.

METHODOLOGY/PRINCIPAL FINDINGS

To elucidate the molecular mechanism of glucose transport and its energetic, a steered molecular dynamics study (SMD) was used. Glucose was pulled from the extracellular end of GLUT4 to the cytoplasm along the pathway using constant velocity pulling method. We identified several key residues within the tunnel that interact directly with either the backbone ring or the hydroxyl groups of glucose. A rotation of glucose molecule was seen near the sugar binding site facilitating the sugar recognition process at the QLS binding site.

CONCLUSIONS/SIGNIFICANCE

This study proposes a possible glucose transport pathway and aids the identification of several residues that make direct interactions with glucose during glucose transport. Mutational studies are required to further validate the observation made in this study.

Direct manipulation of a single potassium channel gate with an atomic force microscope probe

Ion channels are membrane proteins that regulate cell functions by controlling the ion permeability of cell membranes. An ion channel contains an ion-selective pore that permeates ions and a sensor that senses a specific stimulus such as ligand binding to regulate the permeability. The detailed molecular mechanisms of this regulation, or gating, are unknown. Gating is thought to occur from conformational changes in the sensor domain in response to the stimulus, which results in opening the gate to permit ion conduction. Using an atomic force microscope and artificial bilayer system, a mechanical stimulus is applied to a potassium channel, and its gating is monitored in real time. The channel-open probability increases greatly when pushing the cytoplasmic domain toward the membrane. This result shows that a mechanical stimulus at the cytoplasmic domain causes changes in the gating and is the first to show direct evidence of coupling between conformational changes in the cytoplasmic domain and channel gating. This novel technology has the potential to be a powerful tool for investigating the activation dynamics in channel proteins.

Cooperative nature of gating transitions in K(+) channels as seen from dynamic importance sampling calculations

The growing dataset of K(+) channel x-ray structures provides an excellent opportunity to begin a detailed molecular understanding of voltage-dependent gating. These structures, while differing in sequence, represent either a stable open or closed state. However, an understanding of the molecular details of gating will require models for the transitions and experimentally testable predictions for the gating transition. To explore these ideas, we apply dynamic importance sampling to a set of homology models for the molecular conformations of K(+) channels for four different sets of sequences and eight different states. In our results, we highlight the importance of particular residues upstream from the Pro-Val-Pro (PVP) region to the gating transition. This supports growing evidence that the PVP region is important for influencing the flexibility of the S6 helix and thus the opening of the gating domain. The results further suggest how gating on the molecular level depends on intra-subunit motions to influence the cooperative behavior of all four subunits of the K(+) channel. We hypothesize that the gating process occurs in steps: first sidechain movement, then inter-S5-S6 subunit motions, and lastly the large-scale domain rearrangements.

Symmetry-restrained molecular dynamics simulations improve homology models of potassium channels

Most crystallized homo-oligomeric ion channels are highly symmetric, which dramatically decreases conformational space and facilitates building homology models (HMs). However, in molecular dynamics (MD) simulations channels deviate from ideal symmetry and accumulate thermal defects, which complicate the refinement of HMs using MD. In this work we evaluate the ability of symmetry constrained MD simulations to improve HMs accuracy, using an approach conceptually similar to Critical Assessment of techniques for protein Structure Prediction (CASP) competition: build HMs of channels with known structure and evaluate the efficiency of proposed methods in improving HMs accuracy (measured as deviation from experimental structure). Results indicate that unrestrained MD does not improve the accuracy of HMs, instantaneous symmetrization improves accuracy but not stability of HMs during subsequent unrestrained MD, while gradually imposing symmetry constraints improves both accuracy (by 5-50%) and stability of HMs. Moreover, accuracy and stability are strongly correlated, making stability a reliable criterion in predicting the accuracy of new HMs. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Molecular Biomechanics: The Molecular Basis of How Forces Regulate Cellular Function

Recent advances have led to the emergence of molecular biomechanics as an essential element of modern biology. These efforts focus on theoretical and experimental studies of the mechanics of proteins and nucleic acids, and the understanding of the molecular mechanisms of stress transmission, mechanosensing and mechanotransduction in living cells. In particular, single-molecule biomechanics studies of proteins and DNA, and mechanochemical coupling in biomolecular motors have demonstrated the critical importance of molecular mechanics as a new frontier in bioengineering and life sciences. To stimulate a more systematic study of the basic issues in molecular biomechanics, and attract a broader range of researchers to enter this emerging field, here we discuss its significance and relevance, describe the important issues to be addressed and the most critical questions to be answered, summarize both experimental and theoretical/computational challenges, and identify some short-term and long-term goals for the field. The needs to train young researchers in molecular biomechanics with a broader knowledge base, and to bridge and integrate molecular, subcellular and cellular level studies of biomechanics are articulated.