1
|
Khalid S, Schroeder C, Bond PJ, Duncan AL. What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes? MICROBIOLOGY (READING, ENGLAND) 2022; 168. [PMID: 35294337 PMCID: PMC9558347 DOI: 10.1099/mic.0.001165] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.
Collapse
Affiliation(s)
- Syma Khalid
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | - Cyril Schroeder
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | - Peter J Bond
- Bioinformatics Institute (A*STAR), Singapore 138671, Singapore.,Department of Biological Sciences, National University of Singapore, Singapore 117543, Singapore
| | - Anna L Duncan
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| |
Collapse
|
2
|
Microbial cell surface display of oxidoreductases: Concepts and applications. Int J Biol Macromol 2020; 165:835-841. [DOI: 10.1016/j.ijbiomac.2020.09.237] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2020] [Revised: 09/21/2020] [Accepted: 09/27/2020] [Indexed: 12/17/2022]
|
3
|
Rahmat Ullah S, Majid M, Rashid MI, Mehmood K, Andleeb S. Immunoinformatics Driven Prediction of Multiepitopic Vaccine Against Klebsiella pneumoniae and Mycobacterium tuberculosis Coinfection and Its Validation via In Silico Expression. Int J Pept Res Ther 2020; 27:987-999. [PMID: 33281529 PMCID: PMC7703501 DOI: 10.1007/s10989-020-10144-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/16/2020] [Indexed: 01/20/2023]
Abstract
Klebsiella pneumoniae and Mycobacterium tuberculosis coinfection is one of the most lethal combinations that has been becoming frequent yet, not diagnosed and reported properly. Due to the simultaneous occurrence of both infections, diagnosis is delayed leading to inadequate treatments and mortality. With the rise of MDR Klebsiella and Mycobacterium, a prophylactic and an immunotherapeutic vaccine has to be entailed for preemptive and adroit therapeutic approach. In this study, we aim to implement reverse vaccinology approach that encompasses a comprehensive evaluation of vital aspects of the pathogens to explore immunogenic epitopes against Omp A of Klebsiella and Rv1698, Rv1973 of Mtb that may help in vaccine development. The designed multi-epitopic vaccine was assessed for antigenicity, allergenicity and various physiochemical parameters. Molecular docking and simulations were executed to assess the immunogenicity and complex stability of the vaccine. The final multi-epitopic vaccine is validated to be highly immunogenic and can serve as a valuable proactive remedy for subject pathogens.
Collapse
Affiliation(s)
- Sidra Rahmat Ullah
- Department of Industrial Biotechnology, Atta ur Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology (NUST), Islamabad, Pakistan
| | - Mahnoor Majid
- Department of Industrial Biotechnology, Atta ur Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology (NUST), Islamabad, Pakistan
| | - Muhammad Ibrahim Rashid
- Department of Industrial Biotechnology, Atta ur Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology (NUST), Islamabad, Pakistan
- Institute of Basic Medical Sciences, Khyber Medical University (KMU), Peshawar, Pakistan
| | - Khalid Mehmood
- Department of Pharmaceutics, College of Pharmacy, Kingdom of Saudi Arabia, University of Hail, Hail, Kingdom of Saudi Arabia
| | - Saadia Andleeb
- Department of Industrial Biotechnology, Atta ur Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology (NUST), Islamabad, Pakistan
| |
Collapse
|
4
|
Lopez-Barbosa N, Suárez-Arnedo A, Cifuentes J, Gonzalez Barrios AF, Silvera Batista CA, Osma JF, Muñoz-Camargo C, Cruz JC. Magnetite-OmpA Nanobioconjugates as Cell-Penetrating Vehicles with Endosomal Escape Abilities. ACS Biomater Sci Eng 2019; 6:415-424. [PMID: 33463215 DOI: 10.1021/acsbiomaterials.9b01214] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Outer membrane protein A (OmpA) has been extensively studied in Gram-negative bacteria due to its relevance in the adhesion of pathogens to host cells and its surfactant capabilities. It consists of a hydrophobic β-barrel domain and a hydrophilic periplasmic domain, that confers OmpA an amphiphilic structure. This study aims to elucidate the capacity of Escherichia coli OmpA to translocate liposomal membranes and serve as a potential cell-penetrating vehicle. We immobilized OmpA on magnetite nanoparticles and investigated the possible functional changes exhibited by OmpA after immobilization. Liposomal intake was addressed using egg lecithin liposomes as a model, where magnetite-OmpA nanobioconjugates were able to translocate the liposomal membrane and caused a disruptive effect when subjected to a magnetic field. Nanobioconjugates showed both low cytotoxicity and hemolytic tendency. Additional interactions within the intracellular space led to altered viability results via 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT). Confocal microscopy images revealed that immobilized nanoparticles effectively enter the cytoplasm of THP-1 and Vero cells by different routes, and, subsequently, some escape endosomes, lysosomes, and other intracellular compartments with relatively high efficiencies. This was demonstrated by co-localization analyses with LysoTracker green that showed Pearson correlations of about 80 and 28%.
Collapse
Affiliation(s)
| | | | | | | | - Carlos A Silvera Batista
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37235, United States
| | | | | | | |
Collapse
|
5
|
Huber RG, Carpenter TS, Dube N, Holdbrook DA, Ingólfsson HI, Irvine WA, Marzinek JK, Samsudin F, Allison JR, Khalid S, Bond PJ. Multiscale Modeling and Simulation Approaches to Lipid-Protein Interactions. Methods Mol Biol 2019; 2003:1-30. [PMID: 31218611 DOI: 10.1007/978-1-4939-9512-7_1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Lipid membranes play a crucial role in living systems by compartmentalizing biological processes and forming a barrier between these processes and the environment. Naturally, a large apparatus of biomolecules is responsible for construction, maintenance, transport, and degradation of these lipid barriers. Additional classes of biomolecules are tasked with transport of specific substances or transduction of signals from the environment across lipid membranes. In this article, we intend to describe a set of techniques that enable one to build accurate models of lipid systems and their associated proteins, and to simulate their dynamics over a variety of time and length scales. We discuss the methods and challenges that allow us to derive structural, mechanistic, and thermodynamic information from these models. We also show how these models have recently been applied in research to study some of the most complex lipid-protein systems to date, including bacterial and viral envelopes, neuronal membranes, and mammalian signaling systems.
Collapse
Affiliation(s)
- Roland G Huber
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), Singapore, Singapore
| | - Timothy S Carpenter
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Namita Dube
- Department of Chemistry, University of Cambridge, Cambridge, UK
| | - Daniel A Holdbrook
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), Singapore, Singapore
| | - Helgi I Ingólfsson
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - William A Irvine
- Centre for Theoretical Chemistry and Physics, Institute of Natural and Mathematical Sciences, Massey University, Auckland, New Zealand
| | - Jan K Marzinek
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), Singapore, Singapore
| | | | - Jane R Allison
- School of Biological Sciences and Maurice Wilkins Centre for Molecular Biodiscovery, The University of Auckland, Auckland, New Zealand
- Biomolecular Interaction Centre, University of Canterbury, Christchurch, New Zealand
| | - Syma Khalid
- School of Chemistry, University of Southampton, Southampton, UK
| | - Peter J Bond
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), Singapore, Singapore.
- Department of Biological Sciences, National University of Singapore, Singapore, Singapore.
| |
Collapse
|
6
|
Biagini T, Chillemi G, Mazzoccoli G, Grottesi A, Fusilli C, Capocefalo D, Castellana S, Vescovi AL, Mazza T. Molecular dynamics recipes for genome research. Brief Bioinform 2017; 19:853-862. [DOI: 10.1093/bib/bbx006] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2016] [Indexed: 01/17/2023] Open
Affiliation(s)
| | | | | | | | | | | | | | | | - Tommaso Mazza
- IRCCS Casa Sollievo della Sofferenza, Bioinformatics unit, viale Regina Margherita, Rome, Italy
| |
Collapse
|
7
|
Hung HM, Nguyen VP, Ngo ST, Nguyen MT. Theoretical study of the interactions between the first transmembrane segment of NS2 protein and a POPC lipid bilayer. Biophys Chem 2016; 217:1-7. [DOI: 10.1016/j.bpc.2016.07.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2016] [Revised: 07/16/2016] [Accepted: 07/16/2016] [Indexed: 01/22/2023]
|
8
|
Ramakrishna S, Padhi S, Priyakumar UD. Modeling the structure of SARS 3a transmembrane protein using a minimum unfavorable contact approach. J CHEM SCI 2015; 127:2159-2169. [PMID: 32218650 PMCID: PMC7090505 DOI: 10.1007/s12039-015-0982-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2015] [Revised: 09/02/2015] [Accepted: 09/04/2015] [Indexed: 12/11/2022]
Abstract
3a is an accessory protein from SARS coronavirus that is known to play a significant role in the proliferation of the virus by forming tetrameric ion channels. Although the monomeric units are known to consist of three transmembrane (TM) domains, there are no solved structures available for the complete monomer. The present study proposes a structural model for the transmembrane region of the monomer by employing our previously tested approach, which predicts potential orientations of TM α-helices by minimizing the unfavorable contact surfaces between the different TM domains. The best model structure comprising all three α-helices has been subjected to MD simulations to examine its quality. The TM bundle was found to form a compact and stable structure with significant intermolecular interactions. The structural features of the proposed model of 3a account for observations from previous experimental investigations on the activity of the protein. Further analysis indicates that residues from the TM2 and TM3 domains are likely to line the pore of the ion channel, which is in good agreement with a recent experimental study. In the absence of an experimental structure for the protein, the proposed structure can serve as a useful model for inferring structure-function relationships about the protein. The structure of the membrane protein 3a from SARS coronavirus is modeled using an approach that minimizes unfavorable contacts between transmembrane domains. A structure for a complete monomeric form of the protein thereby proposed is able to account for the behavior of the protein reported in previous experimental studies. ![]()
Collapse
Affiliation(s)
- S Ramakrishna
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, 500 032 India
| | - Siladitya Padhi
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, 500 032 India
| | - U Deva Priyakumar
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, 500 032 India
| |
Collapse
|
9
|
Gao J, Chen J. Effects of flanking loops on membrane insertion of transmembrane helices: a role for peptide conformational equilibrium. J Phys Chem B 2013; 117:8330-9. [PMID: 23786317 DOI: 10.1021/jp402356c] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The ability of a transmembrane helix (TMH) to insert into a lipid bilayer has been mainly understood based on the total hydrophobicity of the peptide sequence. Recently, Hedin et al. investigated the influence of flanking loops on membrane insertion of a set of marginally hydrophobic TMHs using translocon-based membrane integration assays. While the flanking loops were found to facilitate the insertion in most cases, counter examples also emerged where the flanking loops hinder membrane insertion and contradict the hydrophobicity and charge distribution analyses. Here, coarse-grained free energy calculations and atomistic simulations were performed to investigate the energetics and conformational details of the membrane insertion of two representative marginally hydrophobic TMHs with (NhaL and EmrL) and without (NhaA and EmrD) the flanking loops. The simulations fail to directly recapitulate the contrasting effects of the flanking loops for these two TMHs, due to systematic overprediction of the stabilities of the transmembrane states that has also been consistently observed in previous studies. Nonetheless, detailed force decomposition and peptide conformation analyses suggest a novel mechanism on how the peptide conformational equilibrium in the aqueous phase may modulate the effects of flanking loops on membrane insertion. Specifically, the flanking loops in peptide EmrL interact strongly with the TMH segment and form stable compact conformations in the aqueous phase, which can hinder membrane absorption and insertion as these processes require extended conformations with minimal interactions between the flanking loops and TMH segment. This work also emphasizes the general importance of considering the peptide conformational equilibrium for understanding the mechanism and energetics of membrane insertion, an aspect that has not yet been sufficiently addressed in the literature.
Collapse
Affiliation(s)
- Jian Gao
- Department of Biochemistry and Molecular Biophysics, Kansas State University, Manhattan, Kansas 66506, United States
| | | |
Collapse
|
10
|
The role of short-chain conjugated poly-(R)-3-hydroxybutyrate (cPHB) in protein folding. Int J Mol Sci 2013; 14:10727-48. [PMID: 23702844 PMCID: PMC3709699 DOI: 10.3390/ijms140610727] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2013] [Revised: 05/09/2013] [Accepted: 05/09/2013] [Indexed: 01/25/2023] Open
Abstract
Poly-(R)-3-hydroxybutyrate (PHB), a linear polymer of R-3-hydroxybutyrate (R-3HB), is a fundamental constituent of biological cells. Certain prokaryotes accumulate PHB of very high molecular weight (10,000 to >1,000,000 residues), which is segregated within granular deposits in the cytoplasm; however, all prokaryotes and all eukaryotes synthesize PHB of medium-chain length (~100-200 residues) which resides within lipid bilayers or lipid vesicles, and PHB of short-chain length (<12 residues) which is conjugated to proteins (cPHB), primarily proteins in membranes and organelles. The physical properties of cPHB indicate it plays important roles in the targeting and folding of cPHB-proteins. Here we review the occurrence, physical properties and molecular characteristics of cPHB, and discuss its influence on the folding and structure of outer membrane protein A (OmpA) of Escherichia coli.
Collapse
|
11
|
Paramo T, Garzón D, Holdbrook DA, Khalid S, Bond PJ. The simulation approach to lipid-protein interactions. Methods Mol Biol 2013; 974:435-455. [PMID: 23404287 DOI: 10.1007/978-1-62703-275-9_19] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.
Collapse
Affiliation(s)
- Teresa Paramo
- Department of Chemistry, Unilever Centre for Molecular Informatics, University of Cambridge, Cambridge, UK
| | | | | | | | | |
Collapse
|
12
|
Abstract
The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.
Collapse
Affiliation(s)
- Syma Khalid
- School of Chemistry, University of Southampton, Southampton, UK
| | | |
Collapse
|
13
|
Maffeo C, Bhattacharya S, Yoo J, Wells D, Aksimentiev A. Modeling and simulation of ion channels. Chem Rev 2012; 112:6250-84. [PMID: 23035940 PMCID: PMC3633640 DOI: 10.1021/cr3002609] [Citation(s) in RCA: 148] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Christopher Maffeo
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - Swati Bhattacharya
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - Jejoong Yoo
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - David Wells
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - Aleksei Aksimentiev
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| |
Collapse
|
14
|
Kim S, Lee Y, Tak HM, Park HJ, Sohn YS, Hwang S, Han J, Kang D, Lee KW. Identification of blocker binding site in mouse TRESK by molecular modeling and mutational studies. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2012. [PMID: 23200789 DOI: 10.1016/j.bbamem.2012.11.021] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
TWIK (tandem-pore domain weak inward rectifying K(+))-related spinal cord K(+) channel, TRESK, a member of the tandem-pore domain K(+) channel family, is the most recently cloned K(2P) channel. TRESK is highly expressed in dorsal root ganglion neuron, a pain sensing neuron, which is a target for analgesics. In this study, a reliable 3D structure for transmembrane (TM) region of mouse TRESK (mTRESK) was constructed, and then the reasonable blocker binding mode of the protein was investigated. The 3D structure of the mTRESK built by homology modeling method was validated with recommend value of stereochemical quality. Based on the validated structure, K(+) channel blocker-bound conformation was obtained by molecular docking and 5ns MD simulation with DPPC lipid bilayer. Our docking study provides the plausible binding mode of known blockers with key interacting residues, especially, F156 and F364. Finally, these modeling results were verified by experimental study with mutation from phenylalanine to alanine (F156A, F364A and F156A/F364A) at the TM2 and TM4. This is the first modeling study for TRESK that can provide structural information of the protein including ligand binding information. These results can be useful in structure based drug design for finding new blockers of the TRESK as potential therapeutic target of pain treatment.
Collapse
Affiliation(s)
- Songmi Kim
- Gyeongsang National University, Gazha-dong, Jinju, Republic of Korea
| | | | | | | | | | | | | | | | | |
Collapse
|
15
|
Monajjemi M, Farahani N, Mollaamin F. Thermodynamic study of solvent effects on nanostructures: phosphatidylserine and phosphatidylinositol membranes. PHYSICS AND CHEMISTRY OF LIQUIDS 2012. [DOI: 10.1080/00319104.2010.527842] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
16
|
Choutko A, Glättli A, Fernández C, Hilty C, Wüthrich K, van Gunsteren WF. Membrane protein dynamics in different environments: simulation study of the outer membrane protein X in a lipid bilayer and in a micelle. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2010; 40:39-58. [PMID: 20922370 DOI: 10.1007/s00249-010-0626-7] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/18/2010] [Revised: 08/23/2010] [Accepted: 08/24/2010] [Indexed: 11/24/2022]
Abstract
The bacterial outer membrane protein OmpX from Escherichia coli has been investigated by molecular dynamics simulations when embedded in a phospholipid bilayer and as a protein-micelle aggregate. The resulting simulation trajectories were analysed in terms of structural and dynamic properties of the membrane protein. In agreement with experimental observations, highest relative stability was found for the β-barrel region that is embedded in the lipophilic phase, whereas an extracellular protruding β-sheet, which is a unique structural feature of OmpX that supposedly plays an important role in cell adhesion and invasion, shows larger structure fluctuations. Additionally, we investigated water permeation into the core of the β-barrel protein, which contains a tight salt-bridge and hydrogen-bond network, so that extensive water flux is unlikely. Differences between the bilayer and the micellar system were observed in the length of the barrel and its position inside the lipid environment, and in the protein interactions with the hydrophilic part of the lipids near the lipid/water interface. Those variations suggest that micelles and other detergent environments might not offer a wholly membrane-like milieu to promote adoption of the physiological conformational state by OmpX.
Collapse
Affiliation(s)
- Alexandra Choutko
- Institute for Physical Chemistry, Swiss Federal Institute of Technology, ETH Zurich, 8093, Zurich, Switzerland
| | | | | | | | | | | |
Collapse
|
17
|
Scheuring S, Dufrêne YF. Atomic force microscopy: probing the spatial organization, interactions and elasticity of microbial cell envelopes at molecular resolution. Mol Microbiol 2010; 75:1327-36. [DOI: 10.1111/j.1365-2958.2010.07064.x] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
18
|
JarosÅawski S, Duquesne K, Sturgis JN, Scheuring S. High-resolution architecture of the outer membrane of the Gram-negative bacteriaRoseobacter denitrificans. Mol Microbiol 2009; 74:1211-22. [DOI: 10.1111/j.1365-2958.2009.06926.x] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
19
|
Ulmschneider MB, Ulmschneider JP. Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides. Mol Membr Biol 2009; 25:245-57. [DOI: 10.1080/09687680802020313] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
20
|
Salt bridges in the miniature viral channel Kcv are important for function. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2009; 39:1057-68. [PMID: 19390850 DOI: 10.1007/s00249-009-0451-z] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2009] [Revised: 03/14/2009] [Accepted: 03/31/2009] [Indexed: 10/20/2022]
Abstract
The viral potassium channel Kcv comprises only 94 amino acids, which represent the pore module of more complex K(+) channels. As for Kir-type channels, Kcv also has a short N-terminal helix exposed to the cytoplasm, upstream of the first transmembrane domain. Here we show that this helix is relevant for Kcv function. The presence of charged amino acids, which form dynamic inter- and intra-subunit salt bridges is crucial. Electrophysiological measurements, yeast rescue experiments and molecular dynamics simulations show that mutants in which the critical salt bridge formation is impaired have no or reduced channel activity. We conclude that these salt bridges destabilise the complexation of K(+) ions by negative charges on the inner transmembrane domain at the entrance into the cavity. This feature facilitates a continuous and coordinated transfer of ions between the cavity and the cytoplasm for channels without the canonical bundle crossing.
Collapse
|
21
|
Ulmschneider JP, Ulmschneider MB. Sampling efficiency in explicit and implicit membrane environments studied by peptide folding simulations. Proteins 2008; 75:586-97. [DOI: 10.1002/prot.22270] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
22
|
Wells DB, Abramkina V, Aksimentiev A. Exploring transmembrane transport through alpha-hemolysin with grid-steered molecular dynamics. J Chem Phys 2007; 127:125101. [PMID: 17902937 PMCID: PMC2888542 DOI: 10.1063/1.2770738] [Citation(s) in RCA: 110] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The transport of biomolecules across cell boundaries is central to cellular function. While structures of many membrane channels are known, the permeation mechanism is known only for a select few. Molecular dynamics (MD) is a computational method that can provide an accurate description of permeation events at the atomic level, which is required for understanding the transport mechanism. However, due to the relatively short time scales accessible to this method, it is of limited utility. Here, we present a method for all-atom simulation of electric field-driven transport of large solutes through membrane channels, which in tens of nanoseconds can provide a realistic account of a permeation event that would require a millisecond simulation using conventional MD. In this method, the average distribution of the electrostatic potential in a membrane channel under a transmembrane bias of interest is determined first from an all-atom MD simulation. This electrostatic potential, defined on a grid, is subsequently applied to a charged solute to steer its permeation through the membrane channel. We apply this method to investigate permeation of DNA strands, DNA hairpins, and alpha-helical peptides through alpha-hemolysin. To test the accuracy of the method, we computed the relative permeation rates of DNA strands having different sequences and global orientations. The results of the G-SMD simulations were found to be in good agreement in experiment.
Collapse
Affiliation(s)
- David B Wells
- Department of Physics, University of Illinois at Urbana--Champaign, Urbana, Illinois 61801, USA.
| | | | | |
Collapse
|
23
|
Abstract
AbstractDespite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.
Collapse
|
24
|
Patargias G, Bond PJ, Deol SS, Sansom MSP. Molecular dynamics simulations of GlpF in a micelle vs in a bilayer: conformational dynamics of a membrane protein as a function of environment. J Phys Chem B 2007; 109:575-82. [PMID: 16851049 DOI: 10.1021/jp046727h] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Octyl glucoside (OG) is a detergent widely employed in structural and functional studies of membrane proteins. To better understand the nature of protein-OG interactions, molecular dynamics simulations (duration 10 ns) have been used to explore an alpha-helical membrane protein, GlpF, in OG micelles and in DMPC bilayers. Greater conformational drift of the extramembraneous protein loops, from the initial X-ray structure, is seen for the GlpF-OG simulations than for the GlpF-DMPC simulation. The mobility of the transmembrane alpha-helices is approximately 1.3x higher in the GlpF-OG than the GlpF-DMPC simulations. The detergent is seen to form an irregular torus around the protein. The presence of the protein leads to a small perturbation in the behavior of the alkyl chains in the OG micelle, namely an approximately 15% increase in the trans-gauche(-)-gauche(+) transition time. Aromatic side chains (Trp, Tyr) and basic side chains (Arg, Lys) play an important role in both protein-detergent (OG) and protein-lipid (DMPC) interactions.
Collapse
Affiliation(s)
- George Patargias
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | | | | | | |
Collapse
|
25
|
Bucher D, Guidoni L, Rothlisberger U. The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study. Biophys J 2007; 93:2315-24. [PMID: 17526559 PMCID: PMC1965436 DOI: 10.1529/biophysj.106.102509] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.
Collapse
Affiliation(s)
- Denis Bucher
- Federal Institute of Technology EPFL, Institute of Chemical Sciences and Engineering, Lausanne, Switzerland
| | | | | |
Collapse
|
26
|
Kóňa J, Minozzi M, Torre V, Carloni P. A gate mechanism indicated in the selectivity filter of the potassium channel KscA. Theor Chem Acc 2007. [DOI: 10.1007/s00214-006-0226-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
27
|
Ulmschneider MB, Ulmschneider JP, Sansom MSP, Di Nola A. A generalized born implicit-membrane representation compared to experimental insertion free energies. Biophys J 2007; 92:2338-49. [PMID: 17218457 PMCID: PMC1864825 DOI: 10.1529/biophysj.106.081810] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.
Collapse
|
28
|
Beevers AJ, Kukol A. Systematic molecular dynamics searching in a lipid bilayer: Application to the glycophorin A and oncogenic ErbB-2 transmembrane domains. J Mol Graph Model 2006; 25:226-33. [PMID: 16434222 DOI: 10.1016/j.jmgm.2005.12.008] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2005] [Revised: 12/15/2005] [Accepted: 12/19/2005] [Indexed: 11/30/2022]
Abstract
Molecular dynamics (MD) simulations of proteins in a lipid bilayer environment are usually undertaken with one or a few starting structures. Here we report a search protocol for systematically exploring the possible interactions in helical bundle transmembrane proteins, a frequently occurring structural motif. The search protocol correctly identifies the experimentally known structure of the dimeric human glycophorin A transmembrane domain as the lowest energy structure among five different models without any prior assumptions, whilst an identical in vacuo search fails to identify the correct structure. The lowest energy structure from the search in a lipid bilayer has a root mean square deviation of 1.1A to the experimental structure. We have applied the same search protocol to the unknown transmembrane structure of the oncogenic mutant ErbB-2 protein, a member of the family of epidermal growth factor receptors. Resulting structures show the role of glutamic acid hydrogen bonding and close helical packing. Water molecules may also play a key role in stabilisation of the transmembrane helix association.
Collapse
Affiliation(s)
- Andrew J Beevers
- Department of Biological Sciences, University of Warwick, Coventry CV4 7AL, UK
| | | |
Collapse
|
29
|
Bond PJ, Faraldo-Gómez JD, Deol SS, Sansom MSP. Membrane protein dynamics and detergent interactions within a crystal: a simulation study of OmpA. Proc Natl Acad Sci U S A 2006; 103:9518-23. [PMID: 16766663 PMCID: PMC1480439 DOI: 10.1073/pnas.0600398103] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Molecular dynamics (MD) simulations are used to explore the dynamics of a membrane protein in its crystal environment. A 50-ns-duration simulation (at a temperature of 300 K) is performed for the crystallographic unit cell of the bacterial outer membrane protein OmpA. The unit cell contains four protein molecules, plus detergent molecules and water. An excellent correlation between simulated and experimental values of crystallographic B factors is observed. Effectively, 0.2 micros of protein trajectories are obtained, allowing a critical assessment of simulation quality. Some deficiency in conformational sampling is demonstrated, but averaging over multiple trajectories improves this limitation. The previously undescribed structure and dynamics of detergent molecules in a unit cell are reported here, providing insight into the interactions important in the formation and stabilization of the crystalline environment at room temperature. In particular, we show that at room temperature the detergent molecules form a dynamic, extended micellar structure spreading over adjacent OmpA monomers within the crystal.
Collapse
Affiliation(s)
- Peter J. Bond
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - José D. Faraldo-Gómez
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Sundeep S. Deol
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
- To whom correspondence should be addressed. E-mail:
| |
Collapse
|
30
|
Nishizawa M, Nishizawa K. Interaction between K+ channel gate modifier hanatoxin and lipid bilayer membranes analyzed by molecular dynamics simulation. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2006; 35:373-81. [PMID: 16453153 DOI: 10.1007/s00249-006-0044-z] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2005] [Revised: 01/03/2006] [Accepted: 01/06/2006] [Indexed: 11/28/2022]
Abstract
Hanatoxin (HaTx) is an ellipsoidal-shaped peptide that binds to the voltage sensor of voltage-dependent channels. Of physicochemical interest, HaTx has a "ring" of charged residues around its periphery and a hydrophobic protrusion. It has previously been postulated that HaTx binds to and functions on the surface of membranes, but a recent fluorescent-quenching study has implied a fairly deep positioning of HaTx in the lipid bilayer membrane. We carried out numerous molecular dynamic simulations of HaTx1, a well-studied variant of HaTx, in fully hydrated phospholipid bilayers. The system reproduced the surface-binding mode of HaTx1, in which HaTx1 resided in the extracellular side (outer) of the water/membrane interface with the hydrophobic patch of HaTx1 facing the membrane interior. On the other hand, analyses with various parameter settings suggested that the surface-binding mode was unstable because of the substantial attractive electrostatic force between HaTx1 and the lipid head groups of the inner (opposite) leaflet. Compared with this electrostatic force, the energetic cost for membrane deformation involving meniscus formation appeared to be small. In an attempt to interpret the quenching data, we consider the possibility of dimpling (meniscus formation) that brings HaTx1 inward (only ~0.7-0.8 nm above the bilayer center), while accounting for the flexibility of both leaflets of the membrane and the long-range interaction between positively charged residues of the membrane-bound peptide and the polar head groups of the opposite leaflet of the membrane. It is suggested that molecular dynamics simulations taking into account the flexibility of the membrane surface is potentially useful in interpreting the fluorescence-quenching data.
Collapse
Affiliation(s)
- Manami Nishizawa
- Department of Biochemistry, Teikyo University School of Medicine, Kaga, Itabashi, 173-8605, Tokyo, Japan.
| | | |
Collapse
|
31
|
Ulmschneider MB, Sansom MSP, Di Nola A. Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques. Biophys J 2005; 90:1650-60. [PMID: 16339877 PMCID: PMC1367316 DOI: 10.1529/biophysj.105.065367] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.
Collapse
|
32
|
Domene C, Grottesi A, Sansom MSP. Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1. Biophys J 2005; 87:256-67. [PMID: 15240462 PMCID: PMC1304348 DOI: 10.1529/biophysj.104.039917] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.
Collapse
Affiliation(s)
- Carmen Domene
- Laboratory of Molecular Biophysics, Department of Biochemistry, University of Oxford, Oxford, OX1 3QU United Kingdom
| | | | | |
Collapse
|
33
|
Ash WL, Zlomislic MR, Oloo EO, Tieleman DP. Computer simulations of membrane proteins. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2005; 1666:158-89. [PMID: 15519314 DOI: 10.1016/j.bbamem.2004.04.012] [Citation(s) in RCA: 204] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2004] [Accepted: 04/29/2004] [Indexed: 11/30/2022]
Abstract
Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.
Collapse
Affiliation(s)
- Walter L Ash
- Department of Biological Sciences, University of Calgary, 2500 University Dr. NW, Calgary AB, Canada T2N 1N4
| | | | | | | |
Collapse
|
34
|
Deol SS, Bond PJ, Domene C, Sansom MSP. Lipid-protein interactions of integral membrane proteins: a comparative simulation study. Biophys J 2004; 87:3737-49. [PMID: 15465855 PMCID: PMC1304887 DOI: 10.1529/biophysj.104.048397] [Citation(s) in RCA: 85] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.
Collapse
Affiliation(s)
- Sundeep S Deol
- Department of Biochemistry, and Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford, United Kingdom
| | | | | | | |
Collapse
|
35
|
Abstract
Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.
Collapse
Affiliation(s)
- Marc Baaden
- Laboratory of Molecular Biophysics, Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | | |
Collapse
|