1
|
Hardy BJ, Curnow P. Computational design of de novo bioenergetic membrane proteins. Biochem Soc Trans 2024; 52:1737-1745. [PMID: 38958574 DOI: 10.1042/bst20231347] [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] [Received: 03/14/2024] [Revised: 06/11/2024] [Accepted: 06/17/2024] [Indexed: 07/04/2024]
Abstract
The major energy-producing reactions of biochemistry occur at biological membranes. Computational protein design now provides the opportunity to elucidate the underlying principles of these processes and to construct bioenergetic pathways on our own terms. Here, we review recent achievements in this endeavour of 'synthetic bioenergetics', with a particular focus on new enabling tools that facilitate the computational design of biocompatible de novo integral membrane proteins. We use recent examples to showcase some of the key computational approaches in current use and highlight that the overall philosophy of 'surface-swapping' - the replacement of solvent-facing residues with amino acids bearing lipid-soluble hydrophobic sidechains - is a promising avenue in membrane protein design. We conclude by highlighting outstanding design challenges and the emerging role of AI in sequence design and structure ideation.
Collapse
Affiliation(s)
| | - Paul Curnow
- School of Biochemistry, University of Bristol, Bristol, U.K
| |
Collapse
|
2
|
Dhiman R, Perera RS, Poojari CS, Wiedemann HTA, Kappl R, Kay CWM, Hub JS, Schrul B. Hairpin protein partitioning from the ER to lipid droplets involves major structural rearrangements. Nat Commun 2024; 15:4504. [PMID: 38802378 PMCID: PMC11130287 DOI: 10.1038/s41467-024-48843-8] [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: 01/28/2024] [Accepted: 05/15/2024] [Indexed: 05/29/2024] Open
Abstract
Lipid droplet (LD) function relies on proteins partitioning between the endoplasmic reticulum (ER) phospholipid bilayer and the LD monolayer membrane to control cellular adaptation to metabolic changes. It has been proposed that these hairpin proteins integrate into both membranes in a similar monotopic topology, enabling their passive lateral diffusion during LD emergence at the ER. Here, we combine biochemical solvent-accessibility assays, electron paramagnetic resonance spectroscopy and intra-molecular crosslinking experiments with molecular dynamics simulations, and determine distinct intramembrane positionings of the ER/LD protein UBXD8 in ER bilayer and LD monolayer membranes. UBXD8 is deeply inserted into the ER bilayer with a V-shaped topology and adopts an open-shallow conformation in the LD monolayer. Major structural rearrangements are required to enable ER-to-LD partitioning. Free energy calculations suggest that such structural transition is unlikely spontaneous, indicating that ER-to-LD protein partitioning relies on more complex mechanisms than anticipated and providing regulatory means for this trans-organelle protein trafficking.
Collapse
Affiliation(s)
- Ravi Dhiman
- Medical Biochemistry and Molecular Biology, Center for Molecular Signaling (PZMS), Faculty of Medicine, Saarland University, 66421, Homburg/Saar, Germany
| | - Rehani S Perera
- Medical Biochemistry and Molecular Biology, Center for Molecular Signaling (PZMS), Faculty of Medicine, Saarland University, 66421, Homburg/Saar, Germany
| | - Chetan S Poojari
- Theoretical Physics and Center for Biophysics, Saarland University, 66123, Saarbrücken, Germany
| | - Haakon T A Wiedemann
- Physical Chemistry and Chemistry Education, Saarland University, 66123, Saarbrücken, Germany
| | - Reinhard Kappl
- Department of Biophysics, Center for Integrative Physiology and Molecular Medicine (CIPMM), Faculty of Medicine, Saarland University, 66421, Homburg/Saar, Germany
| | - Christopher W M Kay
- Physical Chemistry and Chemistry Education, Saarland University, 66123, Saarbrücken, Germany
- London Centre for Nanotechnology, University College London, WC1H 0AH, London, UK
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, 66123, Saarbrücken, Germany
| | - Bianca Schrul
- Medical Biochemistry and Molecular Biology, Center for Molecular Signaling (PZMS), Faculty of Medicine, Saarland University, 66421, Homburg/Saar, Germany.
| |
Collapse
|
3
|
Ozturk TN, König M, Carpenter TS, Pedersen KB, Wassenaar TA, Ingólfsson HI, Marrink SJ. Building complex membranes with Martini 3. Methods Enzymol 2024; 701:237-285. [PMID: 39025573 DOI: 10.1016/bs.mie.2024.03.010] [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
The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.
Collapse
Affiliation(s)
- Tugba Nur Ozturk
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, United States
| | - Melanie König
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, The Netherlands
| | - Timothy S Carpenter
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, United States
| | | | - Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, The Netherlands; Institute for Life Science and Technology, Hanze University of Applied Sciences, Groningen, The Netherlands
| | - Helgi I Ingólfsson
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, United States.
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, The Netherlands.
| |
Collapse
|
4
|
Blazhynska M, Gumbart JC, Chen H, Tajkhorshid E, Roux B, Chipot C. A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes. J Chem Theory Comput 2023; 19:9077-9092. [PMID: 38091976 PMCID: PMC11145395 DOI: 10.1021/acs.jctc.3c00941] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.
Collapse
Affiliation(s)
- Marharyta Blazhynska
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex 54506, France
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology, 837 State Street, Atlanta, Georgia 30332, United States
| | - Haochuan Chen
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex 54506, France
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 N. Mathews Avenue, Urbana, Illinois 61801, United States
- Department of Biochemistry, University of Illinois at Urbana-Champaign, 600 S. Mathews Avenue, Urbana, Illinois 61801, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 E. 57th Street W225, Chicago, Illinois 60637, United States
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex 54506, France
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 N. Mathews Avenue, Urbana, Illinois 61801, United States
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 E. 57th Street W225, Chicago, Illinois 60637, United States
- Department of Chemistry, The University of Hawai'i at Ma̅noa, 2545 McCarthy Mall, Honolulu, Hawaii 96822, United States
| |
Collapse
|
5
|
Ljubič M, Prašnikar E, Perdih A, Borišek J. All-Atom Simulations Reveal the Intricacies of Signal Transduction upon Binding of the HLA-E Ligand to the Transmembrane Inhibitory CD94/NKG2A Receptor. J Chem Inf Model 2023. [PMID: 37207294 DOI: 10.1021/acs.jcim.3c00249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Natural killer (NK) cells play an important role in the innate immune response against tumors and various pathogens such as viruses and bacteria. Their function is controlled by a wide array of activating and inhibitory receptors, which are expressed on their cell surface. Among them is a dimeric NKG2A/CD94 inhibitory transmembrane (TM) receptor which specifically binds to the non-classical MHC I molecule HLA-E, which is often overexpressed on the surface of senescent and tumor cells. Using the Alphafold 2 artificial intelligence system, we constructed the missing segments of the NKG2A/CD94 receptor and generated its complete 3D structure comprising extracellular (EC), TM, and intracellular regions, which served as a starting point for the multi-microsecond all-atom molecular dynamics simulations of the receptor with and without the bound HLA-E ligand and its nonameric peptide. The simulated models revealed that an intricate interplay of events is taking place between the EC and TM regions ultimately affecting the intracellular immunoreceptor tyrosine-based inhibition motif (ITIM) regions that host the point at which the signal is transmitted further down the inhibitory signaling cascade. Signal transduction through the lipid bilayer was also coupled with the changes in the relative orientation of the NKG2A/CD94 TM helices in response to linker reorganization, mediated by fine-tuned interactions in the EC region of the receptor, taking place after HLA-E binding. This research provides atomistic details of the cells' protection mechanism against NK cells and broadens the knowledge regarding the TM signaling of ITIM-bearing receptors.
Collapse
Affiliation(s)
- Martin Ljubič
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
- Faculty of Pharmacy, University of Ljubljana, Aškerčeva 7, 1000 Ljubljana, Slovenia
| | - Eva Prašnikar
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| | - Andrej Perdih
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
- Faculty of Pharmacy, University of Ljubljana, Aškerčeva 7, 1000 Ljubljana, Slovenia
| | - Jure Borišek
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| |
Collapse
|
6
|
Lanrezac A, Baaden M. UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape. MEMBRANES 2023; 13:362. [PMID: 36984749 PMCID: PMC10054542 DOI: 10.3390/membranes13030362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/08/2023] [Revised: 03/13/2023] [Accepted: 03/17/2023] [Indexed: 06/18/2023]
Abstract
The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.
Collapse
|
7
|
Watkins SL. Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry. MEMBRANES 2023; 13:148. [PMID: 36837651 PMCID: PMC9961006 DOI: 10.3390/membranes13020148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Revised: 01/14/2023] [Accepted: 01/18/2023] [Indexed: 06/18/2023]
Abstract
There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.
Collapse
Affiliation(s)
- Stephan L Watkins
- Plant Pathology and CRGB, Oregon State University, 2701 SW Campus Way, Corvallis, OR 97331, USA
| |
Collapse
|
8
|
Corey RA, Baaden M, Chavent M. A brief history of visualizing membrane systems in molecular dynamics simulations. FRONTIERS IN BIOINFORMATICS 2023; 3:1149744. [PMID: 37213533 PMCID: PMC10196259 DOI: 10.3389/fbinf.2023.1149744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2023] [Accepted: 03/13/2023] [Indexed: 05/23/2023] Open
Abstract
Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.
Collapse
Affiliation(s)
- R. A. Corey
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - M. Baaden
- Centre Nationale de la Recherche Scientifique, Laboratoire de Biochimie Théorique, Université Paris Cité, Paris, France
| | - M. Chavent
- Institut de Pharmacologie et Biologie Structurale, CNRS, Université de Toulouse, Toulouse, France
- *Correspondence: M. Chavent,
| |
Collapse
|
9
|
Wang H, Li M, St Onge CM, Fuss B, Zhang Y. Elucidating the binding mechanism of LPA species and analogs in an LPA 4 receptor homology model. J Mol Graph Model 2022; 116:108274. [PMID: 35868118 DOI: 10.1016/j.jmgm.2022.108274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2022] [Revised: 07/04/2022] [Accepted: 07/10/2022] [Indexed: 12/15/2022]
Abstract
Lysophosphatidic acid receptor 4 (LPA4) has emerged as a potential therapeutic target for the treatment of a variety of diseases, including cancer and obesity-induced diabetes, but its structure remains to be revealed. In the present work, a homology model of LPA4 was built for studying the binding mechanism of LPA species and analogs. Then five selected LPA species and analogs with structural variations in their phosphate groups, substitutions on the glycerol backbone, and fatty acyl chains were docked into the LPA4 model, followed by molecular dynamics simulations and energy analyses. The computational results revealed that the aliphatic residues located at the vertical cleft of LPA4 may form a hydrophobic environment for the fatty acyl moiety of LPA species and their analogs. Meanwhile, the positively charged residues in the central cavity of LPA4 may form ionic interactions with the negatively charged hydrophilic head group of LPA species and their analogs. In addition, it was noted that a different binding mode of the hydrophilic head group in each species with the central cavity of the LPA4 might lead to a special rearrangement of the fatty acyl moiety. Taken together, these results may facilitate understanding of the activation mechanism of LPA4 and help design selective ligands to modulate its function for therapeutic purposes.
Collapse
Affiliation(s)
- Huiqun Wang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA, 23298, United States
| | - Mengchu Li
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA, 23298, United States
| | - Celsey M St Onge
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA, 23298, United States
| | - Babette Fuss
- Department of Anatomy and Neurobiology, Virginia Commonwealth University Medical Center, Richmond, VA, 23298, United States
| | - Yan Zhang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA, 23298, United States.
| |
Collapse
|
10
|
Licari G, Dehghani-Ghahnaviyeh S, Tajkhorshid E. Membrane Mixer: A Toolkit for Efficient Shuffling of Lipids in Heterogeneous Biological Membranes. J Chem Inf Model 2022; 62:986-996. [PMID: 35104125 PMCID: PMC8892574 DOI: 10.1021/acs.jcim.1c01388] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.
Collapse
Affiliation(s)
- Giuseppe Licari
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States,Current address: Pharmaceutical Development Biologicals, Boehringer Ingelheim Pharmaceuticals, Inc., Biberach An Der Riß, Germany,Contributed equally to this work
| | - Sepehr Dehghani-Ghahnaviyeh
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States,Contributed equally to this work
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States
| |
Collapse
|
11
|
Róg T, Girych M, Bunker A. Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design. Pharmaceuticals (Basel) 2021; 14:1062. [PMID: 34681286 PMCID: PMC8537670 DOI: 10.3390/ph14101062] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/14/2021] [Accepted: 10/15/2021] [Indexed: 11/17/2022] Open
Abstract
We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.
Collapse
Affiliation(s)
- Tomasz Róg
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Mykhailo Girych
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Alex Bunker
- Drug Research Program, Division of Pharmaceutical Biosciences, Faculty of Pharmacy, University of Helsinki, 00014 Helsinki, Finland;
| |
Collapse
|
12
|
Weigle AT, Carr M, Shukla D. Impact of Increased Membrane Realism on Conformational Sampling of Proteins. J Chem Theory Comput 2021; 17:5342-5357. [PMID: 34339605 DOI: 10.1021/acs.jctc.1c00276] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.
Collapse
Affiliation(s)
- Austin T Weigle
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Matthew Carr
- Independent Software Development Provider310 East Marlette Avenue, Phoenix, Arizona 85012, United States
| | - Diwakar Shukla
- Department of Chemical & Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Cancer Center at Illinois, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Center for Digital Agriculture, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,NIH Center for Macromolecular Modeling and Bioinformatics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
13
|
Wang H, Cao D, Gillespie JC, Mendez RE, Selley DE, Liu-Chen LY, Zhang Y. Exploring the putative mechanism of allosteric modulations by mixed-action kappa/mu opioid receptor bitopic modulators. Future Med Chem 2021; 13:551-573. [PMID: 33590767 PMCID: PMC8027703 DOI: 10.4155/fmc-2020-0308] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 01/14/2021] [Indexed: 12/26/2022] Open
Abstract
The modulation and selectivity mechanisms of seven mixed-action kappa opioid receptor (KOR)/mu opioid receptor (MOR) bitopic modulators were explored. Molecular modeling results indicated that the 'message' moiety of seven bitopic modulators shared the same binding mode with the orthosteric site of the KOR and MOR, whereas the 'address' moiety bound with different subdomains of the allosteric site of the KOR and MOR. The 'address' moiety of seven bitopic modulators bound to different subdomains of the allosteric site of the KOR and MOR may exhibit distinguishable allosteric modulations to the binding affinity and/or efficacy of the 'message' moiety. Moreover, the 3-hydroxy group on the phenolic moiety of the seven bitopic modulators induced selectivity to the KOR over the MOR.
Collapse
Affiliation(s)
- Huiqun Wang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Danni Cao
- Center for Substance Abuse Research, Temple University Lewis Katz School of Medicine, Philadelphia, PA 19140, USA
| | - James C Gillespie
- Department of Pharmacology & Toxicology, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Rolando E Mendez
- Department of Pharmacology & Toxicology, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Dana E Selley
- Department of Pharmacology & Toxicology, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Lee-Yuan Liu-Chen
- Center for Substance Abuse Research, Temple University Lewis Katz School of Medicine, Philadelphia, PA 19140, USA
| | - Yan Zhang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USA
| |
Collapse
|
14
|
Chebon-Bore L, Sanyanga TA, Manyumwa CV, Khairallah A, Tastan Bishop Ö. Decoding the Molecular Effects of Atovaquone Linked Resistant Mutations on Plasmodium falciparum Cytb-ISP Complex in the Phospholipid Bilayer Membrane. Int J Mol Sci 2021; 22:2138. [PMID: 33670016 PMCID: PMC7926518 DOI: 10.3390/ijms22042138] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2020] [Revised: 02/02/2021] [Accepted: 02/07/2021] [Indexed: 12/19/2022] Open
Abstract
Atovaquone (ATQ) is a drug used to prevent and treat malaria that functions by targeting the Plasmodium falciparum cytochrome b (PfCytb) protein. PfCytb catalyzes the transmembrane electron transfer (ET) pathway which maintains the mitochondrial membrane potential. The ubiquinol substrate binding site of the protein has heme bL, heme bH and iron-sulphur [2FE-2S] cluster cofactors that act as redox centers to aid in ET. Recent studies investigating ATQ resistance mechanisms have shown that point mutations of PfCytb confer resistance. Thus, understanding the resistance mechanisms at the molecular level via computational approaches incorporating phospholipid bilayer would help in the design of new efficacious drugs that are also capable of bypassing parasite resistance. With this knowledge gap, this article seeks to explore the effect of three drug resistant mutations Y268C, Y268N and Y268S on the PfCytb structure and function in the presence and absence of ATQ. To draw reliable conclusions, 350 ns all-atom membrane (POPC:POPE phospholipid bilayer) molecular dynamics (MD) simulations with derived metal parameters for the holo and ATQ-bound -proteins were performed. Thereafter, simulation outputs were analyzed using dynamic residue network (DRN) analysis. Across the triplicate MD runs, hydrophobic interactions, reported to be crucial in protein function were assessed. In both, the presence and absence of ATQ and a loss of key active site residue interactions were observed as a result of mutations. These active site residues included: Met 133, Trp136, Val140, Thr142, Ile258, Val259, Pro260 and Phe264. These changes to residue interactions are likely to destabilize the overall intra-protein residue communication network where the proteins' function could be implicated. Protein dynamics of the ATQ-bound mutant complexes showed that they assumed a different pose to the wild-type, resulting in diminished residue interactions in the mutant proteins. In summary, this study presents insights on the possible effect of the mutations on ATQ drug activity causing resistance and describes accurate MD simulations in the presence of the lipid bilayer prior to conducting inhibitory drug discovery for the PfCytb-iron sulphur protein (Cytb-ISP) complex.
Collapse
Affiliation(s)
| | | | | | | | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Makhanda 6140, South Africa; (L.C.-B.); (T.A.S.); (C.V.M.); (A.K.)
| |
Collapse
|
15
|
Glatzová D, Mavila H, Saija MC, Chum T, Cwiklik L, Brdička T, Cebecauer M. The role of prolines and glycine in the transmembrane domain of LAT. FEBS J 2021; 288:4039-4052. [PMID: 33458942 DOI: 10.1111/febs.15713] [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/18/2020] [Revised: 12/15/2020] [Accepted: 01/13/2021] [Indexed: 12/30/2022]
Abstract
Linker for activation in T cells (LAT) is a critical regulator of T-cell development and function. It organises signalling events at the plasma membrane. However, the mechanism, which controls LAT localisation at the plasma membrane, is not fully understood. Here, we studied the impact of helix-breaking amino acids, two prolines and one glycine, in the transmembrane segment on localisation and function of LAT. Using in silico analysis, confocal and super-resolution imaging and flow cytometry, we demonstrate that central proline residue destabilises transmembrane helix by inducing a kink. The helical structure and dynamics are further regulated by glycine and another proline residue in the luminal part of LAT transmembrane domain. Replacement of these residues with aliphatic amino acids reduces LAT dependence on palmitoylation for sorting to the plasma membrane. However, surface expression of these mutants is not sufficient to recover function of nonpalmitoylated LAT in stimulated T cells. These data indicate that geometry and dynamics of LAT transmembrane segment regulate its localisation and function in immune cells.
Collapse
Affiliation(s)
- Daniela Glatzová
- Department of Biophysical Chemistry, J. Heyrovsky Institute of Physical Chemistry, Czech Academy of Sciences, Prague, Czech Republic.,Laboratory of Leukocyte Signaling, Institute of Molecule Genetics, Czech Academy of Sciences, Prague, Czech Republic.,Faculty of Science, Charles University, Prague, Czech Republic
| | - Harsha Mavila
- Department of Biophysical Chemistry, J. Heyrovsky Institute of Physical Chemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Maria Chiara Saija
- Department of Computational Chemistry, J. Heyrovsky Institute of Physical Chemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Tomáš Chum
- Department of Biophysical Chemistry, J. Heyrovsky Institute of Physical Chemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Lukasz Cwiklik
- Department of Computational Chemistry, J. Heyrovsky Institute of Physical Chemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Tomáš Brdička
- Laboratory of Leukocyte Signaling, Institute of Molecule Genetics, Czech Academy of Sciences, Prague, Czech Republic
| | - Marek Cebecauer
- Department of Biophysical Chemistry, J. Heyrovsky Institute of Physical Chemistry, Czech Academy of Sciences, Prague, Czech Republic
| |
Collapse
|
16
|
Federizon J, Feugmo CGT, Huang WC, He X, Miura K, Razi A, Ortega J, Karttunen M, Lovell JF. Experimental and Computational Observations of Immunogenic Cobalt Porphyrin Lipid Bilayers: Nanodomain-Enhanced Antigen Association. Pharmaceutics 2021; 13:pharmaceutics13010098. [PMID: 33466686 PMCID: PMC7828809 DOI: 10.3390/pharmaceutics13010098] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Revised: 01/08/2021] [Accepted: 01/11/2021] [Indexed: 12/18/2022] Open
Abstract
Cobalt porphyrin phospholipid (CoPoP) can incorporate within bilayers to enable non-covalent surface-display of antigens on liposomes by mixing with proteins bearing a polyhistidine tag (his-tag); however, the mechanisms for how this occurs are poorly understood. These were investigated using the his-tagged model antigen Pfs25, a protein antigen candidate for malaria transmission-blocking vaccines. Pfs25 was found to associate with the small molecule aquocobalamin, a form of vitamin B12 and a cobalt-containing corrin macrocycle, but without particle formation, enabling comparative assessment. Relative to CoPoP liposomes, binding and serum stability studies indicated a weaker association of Pfs25 to aquocobalamin or cobalt nitrilotriacetic acid (Co-NTA) liposomes, which have cobalt displayed in the aqueous phase on lipid headgroups. Antigen internalization by macrophages was enhanced with Pfs25 bound to CoPoP liposomes. Immunization in mice with Pfs25 bound to CoPoP liposomes elicited antibodies that recognized ookinetes and showed transmission-reducing activity. To explore the physical mechanisms involved, we employed molecular dynamics (MD) simulations of bilayers containing phospholipid, cholesterol, as well as either CoPoP or NTA-functionalized lipids. The results show that the CoPoP-containing bilayer creates nanodomains that allow access for a limited but sufficient amount of water molecules that could be replaced by his-tags due to their favorable free energy properties allowing for stabilization. The position of the metal center within the NTA liposomes was much more exposed to the aqueous environment, which could explain its limited capacity for stabilizing Pfs25. This study illustrates the impact of CoPoP-induced antigen particleization in enhancing vaccine efficacy, and provides molecular insights into the CoPoP bilayer properties that enable this.
Collapse
Affiliation(s)
- Jasmin Federizon
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, NY 14260, USA; (J.F.); (W.-C.H.); (X.H.)
| | | | - Wei-Chiao Huang
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, NY 14260, USA; (J.F.); (W.-C.H.); (X.H.)
| | - Xuedan He
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, NY 14260, USA; (J.F.); (W.-C.H.); (X.H.)
| | - Kazutoyo Miura
- Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, MD 20852, USA;
| | - Aida Razi
- Department of Anatomy and Cell Biology, McGill University, Montreal, QC H3A 0C7, Canada; (A.R.); (J.O.)
| | - Joaquin Ortega
- Department of Anatomy and Cell Biology, McGill University, Montreal, QC H3A 0C7, Canada; (A.R.); (J.O.)
| | - Mikko Karttunen
- Department of Chemistry, the University of Western Ontario, London, ON N6A 3K7, Canada;
- Centre for Advanced Materials and Biomaterials Research, the University of Western Ontario, London, ON N6A 3K7, Canada
- Department of Applied Mathematics, the University of Western Ontario, London, ON N6A 5B7, Canada
- Correspondence: (M.K.); (J.F.L.)
| | - Jonathan F. Lovell
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, NY 14260, USA; (J.F.); (W.-C.H.); (X.H.)
- Correspondence: (M.K.); (J.F.L.)
| |
Collapse
|
17
|
Determinants of Endoplasmic Reticulum-to-Lipid Droplet Protein Targeting. Dev Cell 2020; 54:471-487.e7. [PMID: 32730754 DOI: 10.1016/j.devcel.2020.07.001] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2019] [Revised: 05/11/2020] [Accepted: 07/07/2020] [Indexed: 02/06/2023]
Abstract
Lipid droplet (LD) formation from the endoplasmic reticulum (ER) is accompanied by the targeting and accumulation of specific hydrophobic, membrane-embedded proteins on LDs. The determinants of this process are unknown. Here, we study the hydrophobic membrane motifs of two Drosophila melanogaster proteins, GPAT4 and ALG14, that utilize this pathway, and we identify crucial sequence features that mediate LD accumulation. Molecular dynamics simulations and studies in cells reveal that LD targeting of these motifs requires deeply inserted tryptophans that have lower free energy in the LD oil phase and positively charged residues near predicted hairpin hinges that become less constrained in the LD environment. Analyzing hydrophobic motifs from similar LD-targeting proteins, it appears that the distribution of tryptophan and positively charged residues distinguishes them from non-LD-targeting membrane motifs. Our studies identify specific sequence features and principles of hydrophobic membrane motifs that mediate their accumulation on LDs.
Collapse
|
18
|
Loschwitz J, Olubiyi OO, Hub JS, Strodel B, Poojari CS. Computer simulations of protein-membrane systems. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:273-403. [PMID: 32145948 PMCID: PMC7109768 DOI: 10.1016/bs.pmbts.2020.01.001] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.
Collapse
Affiliation(s)
- Jennifer Loschwitz
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Olujide O Olubiyi
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Obafemi Awolowo University, Ile-Ife, Nigeria
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany
| | - Birgit Strodel
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Chetan S Poojari
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany.
| |
Collapse
|
19
|
Novitskaia O, Buslaev P, Gushchin I. Assembly of Spinach Chloroplast ATP Synthase Rotor Ring Protein-Lipid Complex. Front Mol Biosci 2019; 6:135. [PMID: 31850368 PMCID: PMC6896225 DOI: 10.3389/fmolb.2019.00135] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Accepted: 11/13/2019] [Indexed: 11/19/2022] Open
Abstract
Rotor ATPases are large multisubunit membrane protein complexes found in all kingdoms of life. The membrane parts of these ATPases include a ring-like assembly, so-called c-ring, consisting of several subunits c, plugged by a patch of phospholipids. In this report, we use a nature-inspired approach to model the assembly of the spinach (Spinacia oleracea) c14 ring protein-lipid complex, where partially assembled oligomers are pulled toward each other using a biasing potential. The resulting assemblies contain 23 to 26 encapsulated plug lipids, general position of which corresponds well to experimental maps. However, best fit to experimental data is achieved with 15 to 17 lipids inside the c-ring. In all of the simulations, the lipids from one leaflet (loop side of the c subunit) are ordered and static, whereas the lipids from the other leaflet are disordered and dynamic. Spontaneous permeation of water molecules toward Glu61 at the active site is also observed. The presented assembly approach is expected to be generalizable to other protein complexes with encapsulated lipid patches.
Collapse
Affiliation(s)
- Olga Novitskaia
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Pavel Buslaev
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny, Russia.,Nanoscience Center, Department of Chemistry, University of Jyväskylä, Jyväskylä, Finland
| | - Ivan Gushchin
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| |
Collapse
|
20
|
Schott-Verdugo S, Gohlke H. PACKMOL-Memgen: A Simple-To-Use, Generalized Workflow for Membrane-Protein-Lipid-Bilayer System Building. J Chem Inf Model 2019; 59:2522-2528. [PMID: 31120747 DOI: 10.1021/acs.jcim.9b00269] [Citation(s) in RCA: 102] [Impact Index Per Article: 20.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present PACKMOL-Memgen, a simple-to-use, generalized workflow for automated building of membrane-protein-lipid-bilayer systems based on open-source tools including Packmol, memembed, pdbremix, and AmberTools. Compared with web-interface-based related tools, PACKMOL-Memgen allows setup of multiple configurations of a system in a user-friendly and efficient manner within minutes. The generated systems are well-packed and thus well-suited as starting configurations in MD simulations under periodic boundary conditions, requiring only moderate equilibration times. PACKMOL-Memgen is distributed with AmberTools and runs on most computing platforms, and its output can also be used for CHARMM or adapted to other molecular-simulation packages.
Collapse
Affiliation(s)
- Stephan Schott-Verdugo
- Institute for Pharmaceutical and Medicinal Chemistry , Heinrich Heine University Düsseldorf , 40225 Düsseldorf , Germany.,Centro de Bioinformática y Simulación Molecular (CBSM), Faculty of Engineering , Universidad de Talca , 1 Poniente 1141 , Casilla 721 , Talca , Chile
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry , Heinrich Heine University Düsseldorf , 40225 Düsseldorf , Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC) & Institute for Complex Systems-Structural Biochemistry (ICS 6) , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany
| |
Collapse
|
21
|
Kirsch SA, Böckmann RA. Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition. Biophys J 2019; 116:2131-2148. [PMID: 31103234 PMCID: PMC6554532 DOI: 10.1016/j.bpj.2019.04.024] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Revised: 04/17/2019] [Accepted: 04/18/2019] [Indexed: 12/29/2022] Open
Abstract
Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (Tm). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to Tm. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.
Collapse
Affiliation(s)
- Sonja A Kirsch
- Computational Biology, Friedrich-Alexander University of Erlangen-Nürnberg, Erlangen, Germany
| | - Rainer A Böckmann
- Computational Biology, Friedrich-Alexander University of Erlangen-Nürnberg, Erlangen, Germany.
| |
Collapse
|
22
|
Resolving the conformational dynamics of ErbB growth factor receptor dimers. J Struct Biol 2019; 207:225-233. [PMID: 31163211 DOI: 10.1016/j.jsb.2019.05.013] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Revised: 05/20/2019] [Accepted: 05/31/2019] [Indexed: 12/30/2022]
Abstract
The combinatorial dimerization of the ErbB growth factor receptors (ErbB1- ErbB4) are critical for their function. Here, we have characterized the conformational dynamics of ErbB transmembrane homo-dimers and hetero-dimers by using a coarse-grain simulation framework. All dimers, except ErbB4-4 and ErbB1-4, exhibit at least two conformations. The reported NMR structures correspond to one of these conformations, representing the N-terminal active state in ErbB1-1 (RH2), ErbB2-2 (RH1) and ErbB4-4 (RH) homo-dimers and the LH dimer in ErbB3-3 homo-dimer, validating the computational approach. Further, we predict a right-handed ErbB3-3 dimer conformer that warrants experimental testing. The five hetero-dimers that have not yet been experimentally resolved display prominent right-handed dimers associating by the SmXXXSm motif. Our results provide insights into the constitutive signaling of ErbB4 after cleavage of the extracellular region. The presence of the inactive-like dimer conformers leading to symmetric kinase domains gives clues on the autoinhibition of the receptor dimers. The dimer states characterized here represent an important step towards understanding the combinatorial cross associations in the ErbB family.
Collapse
|
23
|
Enkavi G, Javanainen M, Kulig W, Róg T, Vattulainen I. Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance. Chem Rev 2019; 119:5607-5774. [PMID: 30859819 PMCID: PMC6727218 DOI: 10.1021/acs.chemrev.8b00538] [Citation(s) in RCA: 175] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Indexed: 12/23/2022]
Abstract
Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.
Collapse
Affiliation(s)
- Giray Enkavi
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - Matti Javanainen
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Institute
of Organic Chemistry and Biochemistry of the Czech Academy
of Sciences, Flemingovo naḿesti 542/2, 16610 Prague, Czech Republic
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Waldemar Kulig
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - Tomasz Róg
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Ilpo Vattulainen
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
- MEMPHYS-Center
for Biomembrane Physics
| |
Collapse
|
24
|
Visualizing Biological Membrane Organization and Dynamics. J Mol Biol 2019; 431:1889-1919. [DOI: 10.1016/j.jmb.2019.02.018] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2018] [Revised: 02/02/2019] [Accepted: 02/13/2019] [Indexed: 11/22/2022]
|
25
|
Melnyk MI, Ivanova IV, Dryn DO, Prylutskyy YI, Hurmach VV, Platonov M, Al Kury LT, Ritter U, Soloviev AI, Zholos AV. C 60 fullerenes selectively inhibit BK Ca but not K v channels in pulmonary artery smooth muscle cells. NANOMEDICINE-NANOTECHNOLOGY BIOLOGY AND MEDICINE 2019; 19:1-11. [PMID: 30981819 DOI: 10.1016/j.nano.2019.03.018] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2018] [Revised: 02/20/2019] [Accepted: 03/19/2019] [Indexed: 01/12/2023]
Abstract
Possessing unique physical and chemical properties, C60 fullerenes are arising as a potential nanotechnological tool that can strongly affect various biological processes. Recent molecular modeling studies have shown that C60 fullerenes can interact with ion channels, but there is lack of data about possible effects of C60 molecule on ion channels expressed in smooth muscle cells (SMC). Here we show both computationally and experimentally that water-soluble pristine C60 fullerene strongly inhibits the large conductance Ca2+-dependent K+ (BKCa), but not voltage-gated K+ (Kv) channels in pulmonary artery SMC. Both molecular docking simulations and analysis of single channel activity indicate that C60 fullerene blocks BKCa channel pore in its open state. In functional tests, C60 fullerene enhanced phenylephrine-induced contraction of pulmonary artery rings by about 25% and reduced endothelium-dependent acetylcholine-induced relaxation by up to 40%. These findings suggest a novel strategy for biomedical application of water-soluble pristine C60 fullerene in vascular dysfunction.
Collapse
Affiliation(s)
- Mariia I Melnyk
- O.O. Bogomoletz Institute of Physiology, National Academy of Sciences of Ukraine, Kyiv, Ukraine; Institute of Pharmacology and Toxicology, National Academy of Medical Sciences of Ukraine, Kyiv, Ukraine
| | - Irina V Ivanova
- Institute of Pharmacology and Toxicology, National Academy of Medical Sciences of Ukraine, Kyiv, Ukraine
| | - Dariia O Dryn
- O.O. Bogomoletz Institute of Physiology, National Academy of Sciences of Ukraine, Kyiv, Ukraine; Institute of Pharmacology and Toxicology, National Academy of Medical Sciences of Ukraine, Kyiv, Ukraine; ESC "Institute of Biology and Medicine", Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Yuriy I Prylutskyy
- ESC "Institute of Biology and Medicine", Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Vasyl V Hurmach
- ESC "Institute of Biology and Medicine", Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Maxim Platonov
- Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Lina T Al Kury
- College of Natural and Health Sciences, Zayed University, Abu Dhabi, United Arab Emirates
| | - Uwe Ritter
- Technical University of Ilmenau, Institute of Chemistry and Biotechnology, Ilmenau, Germany
| | - Anatoly I Soloviev
- Institute of Pharmacology and Toxicology, National Academy of Medical Sciences of Ukraine, Kyiv, Ukraine
| | - Alexander V Zholos
- O.O. Bogomoletz Institute of Physiology, National Academy of Sciences of Ukraine, Kyiv, Ukraine; ESC "Institute of Biology and Medicine", Taras Shevchenko National University of Kyiv, Kyiv, Ukraine.
| |
Collapse
|
26
|
Smith DJ, Klauda JB, Sodt AJ. Simulation Best Practices for Lipid Membranes [Article v1.0]. LIVING JOURNAL OF COMPUTATIONAL MOLECULAR SCIENCE 2019; 1:5966. [PMID: 36204133 PMCID: PMC9534443 DOI: 10.33011/livecoms.1.1.5966] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
We establish a reliable and robust standardization of settings for practical molecular dynamics (MD) simulations of pure and mixed (single- and multi-component) lipid bilayer membranes. In lipid membranes research, particle-based molecular simulations are a powerful tool alongside continuum theory, lipidomics, and model, in vitro, and in vivo experiments. Molecular simulations can provide precise and reproducible spatiotemporal (atomic- and femtosecond-level) information about membrane structure, mechanics, thermodynamics, kinetics, and dynamics. Yet the simulation of lipid membranes can be a daunting task, given the uniqueness of lipid membranes relative to conventional liquid-liquid and solid-liquid interfaces, the immense and complex thermodynamic and statistical mechanical theory, the diversity of multiscale lipid models, limitations of modern computing power, the difficulty and ambiguity of simulation controls, finite size effects, competitive continuum simulation alternatives, and the desired application, including vesicle experiments and biological membranes. These issues can complicate an essential understanding of the field of lipid membranes, and create major bottlenecks to simulation advancement. In this article, we clarify these issues and present a consistent, thorough, and user-friendly framework for the design of state-of-the-art lipid membrane MD simulations. We hope to allow early-career researchers to quickly overcome common obstacles in the field of lipid membranes and reach maximal impact in their simulations.
Collapse
Affiliation(s)
- David J. Smith
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, CA, USA
| | - Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering and Biophysics Program, University of Maryland, College Park, MD, USA
| | - Alexander J. Sodt
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| |
Collapse
|
27
|
Cipcigan F, Carrieri AP, Pyzer-Knapp EO, Krishna R, Hsiao YW, Winn M, Ryadnov MG, Edge C, Martyna G, Crain J. Accelerating molecular discovery through data and physical sciences: Applications to peptide-membrane interactions. J Chem Phys 2018; 148:241744. [DOI: 10.1063/1.5027261] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Flaviu Cipcigan
- IBM Research UK, Hartree Centre, Daresbury WA4 4AD, United Kingdom
| | | | | | - Ritesh Krishna
- IBM Research UK, Hartree Centre, Daresbury WA4 4AD, United Kingdom
| | - Ya-Wen Hsiao
- STFC Daresbury Laboratories, Daresbury WA4 4AD, United Kingdom
| | - Martyn Winn
- STFC Daresbury Laboratories, Daresbury WA4 4AD, United Kingdom
| | - Maxim G. Ryadnov
- National Physical Laboratory, Hampton Road, Teddington, United Kingdom
| | - Colin Edge
- GSK Medicines Research Centre, Stevenage SG1 2NY, United Kingdom
| | - Glenn Martyna
- IBM T. J. Watson Research Center, Yorktown Heights, New York 10598, USA
| | - Jason Crain
- IBM Research UK, Hartree Centre, Daresbury WA4 4AD, United Kingdom
- Maxwell Centre, University of Cambridge, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
28
|
Excessive aggregation of membrane proteins in the Martini model. PLoS One 2017; 12:e0187936. [PMID: 29131844 PMCID: PMC5683612 DOI: 10.1371/journal.pone.0187936] [Citation(s) in RCA: 121] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2017] [Accepted: 10/27/2017] [Indexed: 11/19/2022] Open
Abstract
The coarse-grained Martini model is employed extensively to study membrane protein oligomerization. While this approach is exceptionally promising given its computational efficiency, it is alarming that a significant fraction of these studies demonstrate unrealistic protein clusters, whose formation is essentially an irreversible process. This suggests that the protein-protein interactions are exaggerated in the Martini model. If this held true, then it would limit the applicability of Martini to study multi-protein complexes, as the rapidly clustering proteins would not be able to properly sample the correct dimerization conformations. In this work we first demonstrate the excessive protein aggregation by comparing the dimerization free energies of helical transmembrane peptides obtained with the Martini model to those determined from FRET experiments. Second, we show that the predictions provided by the Martini model for the structures of transmembrane domain dimers are in poor agreement with the corresponding structures resolved using NMR. Next, we demonstrate that the first issue can be overcome by slightly scaling down the Martini protein-protein interactions in a manner, which does not interfere with the other Martini interaction parameters. By preventing excessive, irreversible, and non-selective aggregation of membrane proteins, this approach renders the consideration of lateral dynamics and protein-lipid interactions in crowded membranes by the Martini model more realistic. However, this adjusted model does not lead to an improvement in the predicted dimer structures. This implicates that the poor agreement between the Martini model and NMR structures cannot be cured by simply uniformly reducing the interactions between all protein beads. Instead, a careful amino-acid specific adjustment of the protein-protein interactions is likely required.
Collapse
|
29
|
Lopes D, Jakobtorweihen S, Nunes C, Sarmento B, Reis S. Shedding light on the puzzle of drug-membrane interactions: Experimental techniques and molecular dynamics simulations. Prog Lipid Res 2017; 65:24-44. [DOI: 10.1016/j.plipres.2016.12.001] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Revised: 11/30/2016] [Accepted: 12/03/2016] [Indexed: 12/20/2022]
|
30
|
Vattulainen I, Róg T. Lipid membranes: Theory and simulations bridged to experiments. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2251-2253. [PMID: 27316372 DOI: 10.1016/j.bbamem.2016.06.007] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Ilpo Vattulainen
- Department of Physics, Tampere University of Technology, FI-33101 Tampere, Finland; Department of Physics, University of Helsinki, P. O. Box 64, FI-00014, Finland; MEMPHYS-Center for Biomembrane Physics, University of Southern Denmark, Odense, Denmark.
| | - Tomasz Róg
- Department of Physics, Tampere University of Technology, FI-33101 Tampere, Finland; Department of Physics, University of Helsinki, P. O. Box 64, FI-00014, Finland
| |
Collapse
|