1
|
Bodosa J, Pane AJ, Klauda JB. Modeling asymmetric cell membranes at all-atom resolution. Methods Enzymol 2024; 701:157-174. [PMID: 39025571 DOI: 10.1016/bs.mie.2024.03.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Molecular dynamics (MD) simulations are a useful tool when studying the properties of membranes as they allow for a molecular view of lipid interactions with proteins, nucleic acids, or small molecules. While model membranes are usually symmetric in their lipid composition between leaflets and include a small number of lipid components, physiological membranes are highly complex and vary in the level of asymmetry. Simulation studies have shown that changes in leaflet asymmetry can alter the properties of a membrane. It is therefore necessary to carefully build asymmetric membranes to accurately simulate membranes. This chapter carefully describes the different methods for building asymmetric membranes and the advantages/disadvantages of each method. The simplest methods involve building a membrane with either an equal number of lipids per leaflet or an equal initial surface area (SA) estimated by the area per lipid. More detailed methods include combining two symmetric membranes of equal SA or altering an asymmetric membrane and adjusting the number of lipids after equilibration to minimize an observable such as differential stress (0-DS). More complex methods that require specific simulation software are also briefly described. The challenges and assumptions are listed for each method which should help guide the researcher to choose the best method for their unique MD simulation of an asymmetric membrane.
Collapse
Affiliation(s)
- Jessica Bodosa
- Institute for Physical Science and Technology, Biophysics Program, University of Maryland, College Park, MD, United States
| | - Anthony J Pane
- Institute for Physical Science and Technology, Biophysics Program, University of Maryland, College Park, MD, United States
| | - Jeffery B Klauda
- Institute for Physical Science and Technology, Biophysics Program, University of Maryland, College Park, MD, United States; Department of Chemical and Biomolecular Engineering, College Park, MD, United States.
| |
Collapse
|
2
|
Wan J, Dai L, Xiao H, Zhang W, Zhang R, Xie T, Jia Y, Gao X, Huang J, Liu F. Biological characteristics of mechanosensitive channels MscS and MscL in Actinobacillus pleuropneumoniae. J Bacteriol 2024; 206:e0042923. [PMID: 38391161 PMCID: PMC10955882 DOI: 10.1128/jb.00429-23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Accepted: 01/29/2024] [Indexed: 02/24/2024] Open
Abstract
Actinobacillus pleuropneumoniae is an important respiratory pathogen that can cause porcine contagious pleuropneumonia (PCP), resulting in significant economic losses in swine industry. Microorganisms are subjected to drastic changes in environmental osmolarity. In order to alleviate the drastic rise or fall of osmolarity, cells activate mechanosensitive channels MscL and MscS through tension changes. MscL not only regulates osmotic pressure but also has been reported to secrete protein and uptake aminoglycoside antibiotic. However, MscL and MscS, as the most common mechanosensitive channels, have not been characterized in A. pleuropneumoniae. In this study, the osmotic shock assay showed that MscL increased sodium adaptation by regulating cell length. The results of MIC showed that deletion of mscL decreased the sensitivity of A. pleuropneumoniae to multiple antibiotics, while deletion of mscS rendered A. pleuropneumoniae hypersensitive to penicillin. Biofilm assay demonstrated that MscL contributed the biofilm formation but MscS did not. The results of animal assay showed that MscL and MscS did not affect virulence in vivo. In conclusion, MscL is essential for sodium hyperosmotic tolerance, biofilm formation, and resistance to chloramphenicol, erythromycin, penicillin, and oxacillin. On the other hand, MscS is only involved in oxacillin resistance.IMPORTANCEBacterial resistance to the external environment is a critical function that ensures the normal growth of bacteria. MscL and MscS play crucial roles in responding to changes in both external and internal environments. However, the function of MscL and MscS in Actinobacillus pleuropneumoniae has not yet been reported. Our study shows that MscL plays a significant role in osmotic adaptation, antibiotic resistance, and biofilm formation of A. pleuropneumoniae, while MscS only plays a role in antibiotic resistance. Our findings provide new insights into the functional characteristics of MscL and MscS in A. pleuropneumoniae. MscL and MscS play a role in antibiotic resistance and contribute to the development of antibiotics for A. pleuropneumoniae.
Collapse
Affiliation(s)
- Jiajia Wan
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Lu Dai
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Huasong Xiao
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Wendie Zhang
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Rui Zhang
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Tingting Xie
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Yizhen Jia
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Xuejun Gao
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| | - Jing Huang
- College of Arts and Sciences, Yangtze University, Jingzhou, Hubei, China
| | - Feng Liu
- College of Animal Science and Technology, Yangtze University, Jingzhou, Hubei, China
| |
Collapse
|
3
|
Immadisetty K, Polasa A, Shelton R, Moradi M. Elucidating the molecular basis of spontaneous activation in an engineered mechanosensitive channel. Comput Struct Biotechnol J 2022; 20:2539-2550. [PMID: 35685356 PMCID: PMC9156883 DOI: 10.1016/j.csbj.2022.05.022] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 05/11/2022] [Accepted: 05/12/2022] [Indexed: 12/11/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Kalyan Immadisetty
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, United States
| | - Adithya Polasa
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, United States
| | - Reid Shelton
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, United States
| | - Mahmoud Moradi
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, United States
| |
Collapse
|
4
|
Prasad S, Simmonett AC, Meana-Pañeda R, Brooks BR. The Extended Eighth-Shell method for periodic boundary conditions with rotational symmetry. J Comput Chem 2021; 42:1373-1383. [PMID: 33977553 DOI: 10.1002/jcc.26545] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Revised: 02/18/2021] [Accepted: 03/25/2021] [Indexed: 11/07/2022]
Abstract
The Eighth-Shell method for parallelization of molecular dynamics simulations has previously been shown to be the most optimal for efficiency at large process counts. However, in its current formulation only the P1 space group is supported for periodic boundary conditions (PBC) and thus reflection and/or rotational crystal symmetries are not supported. In this work, we outline the development and implementation of the Extended Eighth-Shell (EES) method that allows rotational symmetry by using an extended import region compared to the ES method. It simulates only the asymmetric unit and communicates coordinates and forces with images that correspond to P21 PBC. The P21 PBC has application in lipid bilayer simulations as it can be used to allow lipids to switch leaftlets, thus rapidly balancing the chemical potential difference between the two layers. Our results show that the EES method scales efficiently over large number of processes and can be used for simulations with P21 symmetry in an orthorhombic crystal.
Collapse
Affiliation(s)
- Samarjeet Prasad
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, USA
| | - Rubén Meana-Pañeda
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute (NHLBI), National Institutes of Health (NIH), Bethesda, Maryland, USA
| |
Collapse
|
5
|
Balleza D, Rosas ME, Romero-Romero S. Voltage vs. Ligand I: Structural basis of the intrinsic flexibility of S3 segment and its significance in ion channel activation. Channels (Austin) 2019; 13:455-476. [PMID: 31647368 PMCID: PMC6833973 DOI: 10.1080/19336950.2019.1674242] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
We systematically predict the internal flexibility of the S3 segment, one of the most mobile elements in the voltage-sensor domain. By analyzing the primary amino acid sequences of V-sensor containing proteins, including Hv1, TPC channels and the voltage-sensing phosphatases, we established correlations between the local flexibility and modes of activation for different members of the VGIC superfamily. Taking advantage of the structural information available, we also assessed structural aspects to understand the role played by the flexibility of S3 during the gating of the pore. We found that S3 flexibility is mainly determined by two specific regions: (1) a short NxxD motif in the N-half portion of the helix (S3a), and (2) a short sequence at the beginning of the so-called paddle motif where the segment has a kink that, in some cases, divide S3 into two distinct helices (S3a and S3b). A good correlation between the flexibility of S3 and the reported sensitivity to temperature and mechanical stretch was found. Thus, if the channel exhibits high sensitivity to heat or membrane stretch, local S3 flexibility is low. On the other hand, high flexibility of S3 is preferentially associated to channels showing poor heat and mechanical sensitivities. In contrast, we did not find any apparent correlation between S3 flexibility and voltage or ligand dependence. Overall, our results provide valuable insights into the dynamics of channel-gating and its modulation.
Collapse
Affiliation(s)
- Daniel Balleza
- Departamento de Química ICET, Universidad Autónoma de Guadalajara , Zapopan Jalisco , Mexico
| | - Mario E Rosas
- Departamento de Química ICET, Universidad Autónoma de Guadalajara , Zapopan Jalisco , Mexico
| | - Sergio Romero-Romero
- Facultad de Medicina, Departamento de Bioquímica, Universidad Nacional Autónoma de México, 04510 Mexico City, Mexico. Current address: Department of Biochemistry, University of Bayreuth , Bayreuth , Germany
| |
Collapse
|
6
|
Muller MP, Jiang T, Sun C, Lihan M, Pant S, Mahinthichaichan P, Trifan A, Tajkhorshid E. Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation. Chem Rev 2019; 119:6086-6161. [PMID: 30978005 PMCID: PMC6506392 DOI: 10.1021/acs.chemrev.8b00608] [Citation(s) in RCA: 153] [Impact Index Per Article: 30.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.
Collapse
Affiliation(s)
- Melanie P. Muller
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Tao Jiang
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Chang Sun
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Muyun Lihan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Shashank Pant
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Paween Mahinthichaichan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Anda Trifan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Emad Tajkhorshid
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
7
|
Zhang Y, Zheng QC. What are the effects of the serine triad on proton conduction of an influenza B M2 channel? An investigation by molecular dynamics simulations. Phys Chem Chem Phys 2019; 21:8820-8826. [PMID: 30968902 DOI: 10.1039/c9cp00612e] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The tetrameric influenza B M2 channel (BM2), an acid activated proton channel, is important in the influenza virus B lifecycle. A conserved HxxxW motif is responsible for proton conduction and channel gating. In this study, to explore the effects of the serine triad (S9, S12 and S16) on proton conduction, we performed classical molecular dynamics (CMD) simulations and adaptive steered molecular dynamics (ASMD) simulations at different protonation states of the H19 tetrad. The results of the pore radius and the C-terminal tilt angle show that the electrostatic repulsion induced by protonated H19 is the key driving force for opening the BM2 channel. The open states could be stabilized by the hydrogen bonds between S16 and protonated H19. The solvent accessible surface area and water density indicate that the polar hydrophilic environment provided by the serine triad facilitates the formation of a water wire, and then exhibits favourable effects on proton conduction. The mutant research verifies and supports these views. Our work clarifies the effects of the serine triad on proton conduction in the BM2 channel, which would help us deeply understand the proton conduction mechanism in BM2 and provides a new perspective for antiviral drug design against BM2.
Collapse
Affiliation(s)
- Yue Zhang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun 130023, People's Republic of China.
| | | |
Collapse
|
8
|
Melo M, Arnarez C, Sikkema H, Kumar N, Walko M, Berendsen HJC, Kocer A, Marrink SJ, Ingólfsson HI. High-Throughput Simulations Reveal Membrane-Mediated Effects of Alcohols on MscL Gating. J Am Chem Soc 2017; 139:2664-2671. [PMID: 28122455 PMCID: PMC5343553 DOI: 10.1021/jacs.6b11091] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2016] [Indexed: 12/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Manuel
N. Melo
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG, Groningen, The Netherlands
| | - Clément Arnarez
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG, Groningen, The Netherlands
| | - Hendrik Sikkema
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG, Groningen, The Netherlands
| | - Neeraj Kumar
- Groningen
Institute for Evolutionary Life Sciences, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
| | - Martin Walko
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
| | - Herman J. C. Berendsen
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG, Groningen, The Netherlands
| | - Armagan Kocer
- Department
of Neuroscience, University Medical Center
Groningen, University of Groningen, Antonius Deusinglaan 1, 99713 AV, Groningen, The
Netherlands
| | - Siewert J. Marrink
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG, Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Groningen
Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG, Groningen, The Netherlands
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG, Groningen, The Netherlands
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California
| |
Collapse
|
9
|
Park S, Beaven AH, Klauda JB, Im W. How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets? J Chem Theory Comput 2016; 11:3466-77. [PMID: 26575780 PMCID: PMC4685941 DOI: 10.1021/acs.jctc.5b00232] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Difficulties in estimating the correct
number of lipids in each
leaflet of complex bilayer membrane simulation systems make it inevitable
to introduce a mismatch in lipid packing (i.e., area per lipid) and
thus alter the lateral pressure of each leaflet. To investigate potential
impacts of such mismatch on simulation results, we performed molecular
dynamics simulations of saturated and monounsaturated lipid bilayers
with and without gramicidin A or WALP23 at various mismatches by adjusting
the number of lipids in the lower leaflet from no mismatch to a 25%
reduction compared to that in the upper leaflet. All simulations were
stable under the constant pressure barostat, but the mismatch induces
asymmetric lipid packing between the leaflets, so that the upper leaflet
becomes more ordered, and the lower leaflet becomes less ordered.
The mismatch impacts on various bilayer properties are mild up to
5–10% mismatch, and bilayers with fully saturated chains appear
to be more prone to these impacts than those with unsaturated tails.
The nonvanishing leaflet surface tensions and the free energy derivatives
with respect to the bilayer curvature indicate that the bilayer would
be energetically unstable in the presence of mismatch. We propose
a quantitative criterion for allowable mismatch based on the energetics
derived from a continuum elastic model, which grows as a square root
of the number of the lipids in the system. On the basis of this criterion,
we infer that the area per lipid mismatch up to 5% would be tolerable
in various membrane simulations of reasonable all-atom system sizes
(40–160 lipids per leaflet).
Collapse
Affiliation(s)
| | | | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering and the Biophysics Program, University of Maryland , College Park, Maryland 20742, United States
| | | |
Collapse
|
10
|
Cheng X, Jo S, Qi Y, Marassi FM, Im W. Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes. Biophys J 2016; 108:1954-62. [PMID: 25902435 DOI: 10.1016/j.bpj.2015.03.012] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2014] [Revised: 02/05/2015] [Accepted: 03/10/2015] [Indexed: 10/23/2022] Open
Abstract
Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys(40), residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.
Collapse
Affiliation(s)
- Xi Cheng
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | - Sunhwan Jo
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | - Yifei Qi
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | | | - Wonpil Im
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas.
| |
Collapse
|
11
|
Wu EL, Fleming PJ, Yeom MS, Widmalm G, Klauda JB, Fleming KG, Im W. E. coli outer membrane and interactions with OmpLA. Biophys J 2014; 106:2493-502. [PMID: 24896129 PMCID: PMC4052237 DOI: 10.1016/j.bpj.2014.04.024] [Citation(s) in RCA: 109] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Revised: 04/17/2014] [Accepted: 04/21/2014] [Indexed: 10/25/2022] Open
Abstract
The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.
Collapse
Affiliation(s)
- Emilia L Wu
- Department of Molecular Biosciences and Center for Bioinformatics, The University of Kansas, Lawrence, Kansas
| | - Patrick J Fleming
- T. C. Jenkins Department of Biophysics, John Hopkins University, Baltimore, Maryland
| | - Min Sun Yeom
- Korean Institute of Science and Technology Information, Daejeon, Korea
| | - Göran Widmalm
- Department of Organic Chemistry and Stockholm Center for Biomembrane Research, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, The University of Maryland, College Park, Maryland
| | - Karen G Fleming
- T. C. Jenkins Department of Biophysics, John Hopkins University, Baltimore, Maryland.
| | - Wonpil Im
- Department of Molecular Biosciences and Center for Bioinformatics, The University of Kansas, Lawrence, Kansas.
| |
Collapse
|
12
|
Ingólfsson HI, Lopez CA, Uusitalo JJ, de Jong DH, Gopal SM, Periole X, Marrink SJ. The power of coarse graining in biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2014; 4:225-248. [PMID: 25309628 PMCID: PMC4171755 DOI: 10.1002/wcms.1169] [Citation(s) in RCA: 325] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Computational modeling of biological systems is challenging because of the multitude of spatial and temporal scales involved. Replacing atomistic detail with lower resolution, coarse grained (CG), beads has opened the way to simulate large-scale biomolecular processes on time scales inaccessible to all-atom models. We provide an overview of some of the more popular CG models used in biomolecular applications to date, focusing on models that retain chemical specificity. A few state-of-the-art examples of protein folding, membrane protein gating and self-assembly, DNA hybridization, and modeling of carbohydrate fibers are used to illustrate the power and diversity of current CG modeling.
Collapse
Affiliation(s)
- Helgi I Ingólfsson
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Cesar A Lopez
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Jaakko J Uusitalo
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Djurre H de Jong
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Srinivasa M Gopal
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Xavier Periole
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| |
Collapse
|
13
|
Martinac B, Nomura T, Chi G, Petrov E, Rohde PR, Battle AR, Foo A, Constantine M, Rothnagel R, Carne S, Deplazes E, Cornell B, Cranfield CG, Hankamer B, Landsberg MJ. Bacterial mechanosensitive channels: models for studying mechanosensory transduction. Antioxid Redox Signal 2014; 20:952-69. [PMID: 23834368 DOI: 10.1089/ars.2013.5471] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Boris Martinac
- 1 Molecular Cardiology and Biophysics Division/Mechanosensory Biophysics Laboratory, Victor Chang Cardiac Research Institute , Darlinghurst, Australia
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
14
|
Sawada Y, Murase M, Sokabe M. The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations: from tension sensing to channel opening. Channels (Austin) 2013; 6:317-31. [PMID: 23146938 DOI: 10.4161/chan.21895] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Yasuyuki Sawada
- Department of Physiology, Nagoya University Graduate School of Medicine, Nagoya, Japan
| | | | | |
Collapse
|
15
|
Wassenaar TA, Ingólfsson HI, Priess M, Marrink SJ, Schäfer LV. Mixing MARTINI: electrostatic coupling in hybrid atomistic-coarse-grained biomolecular simulations. J Phys Chem B 2013; 117:3516-30. [PMID: 23406326 DOI: 10.1021/jp311533p] [Citation(s) in RCA: 119] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Hybrid molecular dynamics simulations of atomistic (AA) solutes embedded in coarse-grained (CG) environment can substantially reduce the computational cost with respect to fully atomistic simulations. However, interfacing both levels of resolution is a major challenge that includes a balanced description of the relevant interactions. This is especially the case for polar solvents such as water, which screen the electrostatic interactions and thus require explicit electrostatic coupling between AA and CG subsystems. Here, we present and critically test computationally efficient hybrid AA/CG models. We combined the Gromos atomistic force field with the MARTINI coarse-grained force field. To enact electrostatic coupling, two recently developed CG water models with explicit electrostatic interactions were used: the polarizable MARTINI water model and the BMW model. The hybrid model was found to be sensitive to the strength of the AA-CG electrostatic coupling, which was adjusted through the relative dielectric permittivity εr(AA-CG). Potentials of mean force (PMFs) between pairs of amino acid side chain analogues in water and partitioning free enthalpies of uncharged amino acid side chain analogues between apolar solvent and water show significant differences between the hybrid simulations and the fully AA or CG simulations, in particular for charged and polar molecules. For apolar molecules, the results obtained with the hybrid AA/CG models are in better agreement with the fully atomistic results. The structures of atomistic ubiquitin solvated in CG water and of a single atomistic transmembrane α-helix and the transmembrane portion of an atomistic mechanosensitive channel in CG lipid bilayers were largely maintained during 50-100 ns of AA/CG simulations, partly due to an overstabilization of intramolecular interactions. This work highlights some key challenges on the way toward hybrid AA/CG models that are both computationally efficient and sufficiently accurate for biomolecular simulations.
Collapse
Affiliation(s)
- Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | | | | | | | | |
Collapse
|
16
|
Deplazes E, Louhivuori M, Jayatilaka D, Marrink SJ, Corry B. Structural investigation of MscL gating using experimental data and coarse grained MD simulations. PLoS Comput Biol 2012; 8:e1002683. [PMID: 23028281 PMCID: PMC3447979 DOI: 10.1371/journal.pcbi.1002683] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2012] [Accepted: 07/21/2012] [Indexed: 12/02/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Evelyne Deplazes
- School of Biomedical, Biomolecular and Chemical Sciences, The University of Western Australia, Perth, Australia
| | | | | | | | | |
Collapse
|