1
|
Madsen JJ, Ohkubo YZ. Elucidating the complex membrane binding of a protein with multiple anchoring domains using extHMMM. PLoS Comput Biol 2024; 20:e1011421. [PMID: 38976709 DOI: 10.1371/journal.pcbi.1011421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2023] [Accepted: 06/19/2024] [Indexed: 07/10/2024] Open
Abstract
Membrane binding is a crucial mechanism for many proteins, but understanding the specific interactions between proteins and membranes remains a challenging endeavor. Coagulation factor Va (FVa) is a large protein whose membrane interactions are complicated due to the presence of multiple anchoring domains that individually can bind to lipid membranes. Using molecular dynamics simulations, we investigate the membrane binding of FVa and identify the key mechanisms that govern its interaction with membranes. Our results reveal that FVa can either adopt an upright or a tilted molecular orientation upon membrane binding. We further find that the domain organization of FVa deviates (sometimes significantly) from its crystallographic reference structure, and that the molecular orientation of the protein matches with domain reorganization to align the C2 domain toward its favored membrane-normal orientation. We identify specific amino acid residues that exhibit contact preference with phosphatidylserine lipids over phosphatidylcholine lipids, and we observe that mostly electrostatic effects contribute to this preference. The observed lipid-binding process and characteristics, specific to FVa or common among other membrane proteins, in concert with domain reorganization and molecular tilt, elucidate the complex membrane binding dynamics of FVa and provide important insights into the molecular mechanisms of protein-membrane interactions. An updated version of the HMMM model, termed extHMMM, is successfully employed for efficiently observing membrane bindings of systems containing the whole FVa molecule.
Collapse
Affiliation(s)
- Jesper J Madsen
- Department of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, Florida, United States of America
- Center for Global Health and Infectious Diseases Research, Global and Planetary Health, College of Public Health, University of South Florida, Tampa, Florida, United States of America
| | - Y Zenmei Ohkubo
- Department of Bioinformatics, School of Life and Natural Sciences, Abdullah Gül University, Kayseri, Turkey
| |
Collapse
|
2
|
Lihan M, Tajkhorshid E. Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions. J Chem Inf Model 2024; 64:4822-4834. [PMID: 38844760 DOI: 10.1021/acs.jcim.4c00619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.
Collapse
Affiliation(s)
- Muyun Lihan
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Emad Tajkhorshid
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Biochemistry, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
3
|
Ohkubo YZ, Radulovic PW, Kahira AN, Madsen JJ. Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V. Curr Res Struct Biol 2024; 7:100149. [PMID: 38766652 PMCID: PMC11098723 DOI: 10.1016/j.crstbi.2024.100149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Revised: 03/28/2024] [Accepted: 04/29/2024] [Indexed: 05/22/2024] Open
Abstract
Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.
Collapse
Affiliation(s)
- Y. Zenmei Ohkubo
- Department of Bioinformatics, School of Life and Natural Sciences, Abdullah Gül University, Kayseri, Turkey
| | - Peter W. Radulovic
- Graduate Programs, Taneja College of Pharmacy, University of South Florida, Tampa, FL 33612, USA
| | - Albert N. Kahira
- Graduate Programs, School of Engineering, Abdullah Gül University, Kayseri, Turkey
| | - Jesper J. Madsen
- Department of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, FL 33612, USA
- Center for Global Health and Infectious Diseases Research, Global and Planetary Health, College of Public Health, University of South Florida, Tampa, FL 33612, USA
| |
Collapse
|
4
|
Hasdemir HS, Pozzi N, Tajkhorshid E. Atomistic Characterization of Beta-2-Glycoprotein I Domain V Interaction with Anionic Membranes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.19.585743. [PMID: 38562685 PMCID: PMC10983932 DOI: 10.1101/2024.03.19.585743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Background Interaction of beta-2-glycoprotein I ( β 2 GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of it's membrane bind-ing domain, Domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood. Objectives This study aims to elucidate the mechanism by which DV of β 2 GPI binds to anionic membranes. Methods We utilized molecular dynamics (MD) simulations to investigate the struc-tural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the HMMM simulations, we conducted additional simulations using a full mem-brane model. Results The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into β 2 GPI's membrane interaction. Simulations revealed two distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of β 2 GPI by antiphospholipid antibodies. Conclusion The study advances our understanding of the molecular interactions be-tween β 2 GPI's DV and anionic membranes, crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant bind-ing mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.
Collapse
|
5
|
Radka CD, Grace CR, Hasdemir HS, Li Y, Rodriguez CC, Rodrigues P, Oldham ML, Qayyum MZ, Pitre A, MacCain WJ, Kalathur RC, Tajkhorshid E, Rock CO. The carboxy terminus causes interfacial assembly of oleate hydratase on a membrane bilayer. J Biol Chem 2024; 300:105627. [PMID: 38211817 PMCID: PMC10847778 DOI: 10.1016/j.jbc.2024.105627] [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: 11/15/2023] [Revised: 12/19/2023] [Accepted: 12/22/2023] [Indexed: 01/13/2024] Open
Abstract
The soluble flavoprotein oleate hydratase (OhyA) hydrates the 9-cis double bond of unsaturated fatty acids. OhyA substrates are embedded in membrane bilayers; OhyA must remove the fatty acid from the bilayer and enclose it in the active site. Here, we show that the positively charged helix-turn-helix motif in the carboxy terminus (CTD) is responsible for interacting with the negatively charged phosphatidylglycerol (PG) bilayer. Super-resolution microscopy of Staphylococcus aureus cells expressing green fluorescent protein fused to OhyA or the CTD sequence shows subcellular localization along the cellular boundary, indicating OhyA is membrane-associated and the CTD sequence is sufficient for membrane recruitment. Using cryo-electron microscopy, we solved the OhyA dimer structure and conducted 3D variability analysis of the reconstructions to assess CTD flexibility. Our surface plasmon resonance experiments corroborated that OhyA binds the PG bilayer with nanomolar affinity and we found the CTD sequence has intrinsic PG binding properties. We determined that the nuclear magnetic resonance structure of a peptide containing the CTD sequence resembles the OhyA crystal structure. We observed intermolecular NOE from PG liposome protons next to the phosphate group to the CTD peptide. The addition of paramagnetic MnCl2 indicated the CTD peptide binds the PG surface but does not insert into the bilayer. Molecular dynamics simulations, supported by site-directed mutagenesis experiments, identify key residues in the helix-turn-helix that drive membrane association. The data show that the OhyA CTD binds the phosphate layer of the PG surface to obtain bilayer-embedded unsaturated fatty acids.
Collapse
Affiliation(s)
- Christopher D Radka
- Department of Microbiology, Immunology, and Molecular Genetics, University of Kentucky, Lexington, Kentucky, USA; Department of Host Microbe Interactions, St Jude Children's Research Hospital, Memphis, Tennessee, USA.
| | - Christy R Grace
- Department of Structural Biology, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - Hale S Hasdemir
- Theoretical and Computational Biophysics Group, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | - Yupeng Li
- Theoretical and Computational Biophysics Group, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | - Carlos C Rodriguez
- Theoretical and Computational Biophysics Group, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | - Patrick Rodrigues
- Hartwell Center of Biotechnology, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - Michael L Oldham
- Department of Structural Biology, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - M Zuhaib Qayyum
- Department of Structural Biology, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - Aaron Pitre
- Cell and Tissue Imaging Center, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - William J MacCain
- Department of Host Microbe Interactions, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - Ravi C Kalathur
- Department of Structural Biology, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | - Charles O Rock
- Department of Host Microbe Interactions, St Jude Children's Research Hospital, Memphis, Tennessee, USA
| |
Collapse
|
6
|
Vo A, Murphy M, Phan P, Stone T, Prabhu R. Molecular Dynamics Simulation of Membrane Systems in the context of Traumatic Brain Injury. CURRENT OPINION IN BIOMEDICAL ENGINEERING 2023. [DOI: 10.1016/j.cobme.2023.100453] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
|
7
|
Muller MP, Morrissey JH, Tajkhorshid E. Molecular View into Preferential Binding of the Factor VII Gla Domain to Phosphatidic Acid. Biochemistry 2022; 61:1694-1703. [PMID: 35853076 PMCID: PMC9637449 DOI: 10.1021/acs.biochem.2c00266] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Factor VII (FVII) is a serine protease with a key role in initiating the coagulation cascade. It is part of a family of vitamin K-dependent clotting proteins, which require vitamin K for formation of their specialized membrane-binding domains (Gla domains). Membrane binding of the FVII Gla domain is critical to the activity of FVII, mediating the formation of its complex with other clotting factors. While Gla domains among coagulation factors are highly conserved in terms of amino acid sequence and structure, they demonstrate differential binding specificity toward anionic lipids. Although most Gla domain-containing clotting proteins display a strong preference for phosphatidylserine (PS), it has been demonstrated that FVII and protein C instead bind preferentially to phosphatidic acid (PA). We have developed the first model of the FVII Gla domain bound to PA lipids in membranes containing PA, the highly mobile membrane mimetic model, which accelerates slow diffusion of lipids in molecular dynamics simulations and therefore facilitates the membrane binding process and enhances sampling of lipid interactions. Simulations were performed using atomic level molecular dynamics, requiring a fixed charge to all atoms. The overall charge assigned to each PA lipid for this study was -1. We also developed an additional model of the FVII Gla domain bound to a 1:1 PS/PC membrane and compared the modes of binding of PS and PA lipids to FVII, allowing us to identify potential PA-specific binding sites.
Collapse
Affiliation(s)
- Melanie P Muller
- 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 61801, United States
| | - James H Morrissey
- Department of Biological Chemistry, University of Michigan Medical School, Ann Arbor, Michigan 48109, United States
| | - 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 61801, United States
| |
Collapse
|
8
|
Wu W, Zhou Q, Xu B, Zhong Y, Xu H, Zhang L, Mao Z. Study on the dispersion mechanism of the polycarboxylic acid dispersant for disperse dyes. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
9
|
Braun L, Schoen I, Vogel V. PIP 2-induced membrane binding of the vinculin tail competes with its other binding partners. Biophys J 2021; 120:4608-4622. [PMID: 34411575 DOI: 10.1016/j.bpj.2021.08.018] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Revised: 06/14/2021] [Accepted: 08/11/2021] [Indexed: 01/09/2023] Open
Abstract
Vinculin plays a key role during the first phase of focal adhesion formation and interacts with the plasma membrane through specific binding of its tail domain to the lipid phosphatidylinositol 4,5-bisphosphate (PIP2). Our understanding of the PIP2-vinculin interaction has been hampered by contradictory biochemical and structural data. Here, we used a multiscale molecular dynamics simulation approach, in which unbiased coarse-grained molecular dynamics were used to generate starting structures for subsequent microsecond-long all-atom simulations. This allowed us to map the interaction of the vinculin tail with PIP2-enriched membranes in atomistic detail. In agreement with experimental data, we have shown that membrane binding is sterically incompatible with the intramolecular interaction between vinculin's head and tail domain. Our simulations further confirmed biochemical and structural results, which identified two positively charged surfaces, the basic collar and the basic ladder, as the main PIP2 interaction sites. By introducing a valency-disaggregated binding network analysis, we were able to map the protein-lipid interactions in unprecedented detail. In contrast to the basic collar, in which PIP2 is specifically recognized by an up to hexavalent binding pocket, the basic ladder forms a series of low-valency binding sites. Importantly, many of these PIP2 binding residues are also involved in maintaining vinculin in a closed, autoinhibited conformation. These findings led us to propose a molecular mechanism for the coupling between vinculin activation and membrane binding. Finally, our refined binding site suggests an allosteric relationship between PIP2 and F-actin binding that disfavors simultaneous interaction with both ligands, despite nonoverlapping binding sites.
Collapse
Affiliation(s)
- Lukas Braun
- Laboratory of Applied Mechanobiology, Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland.
| | - Ingmar Schoen
- School of Pharmacy and Biomolecular Sciences, Irish Centre for Vascular Biology, Royal College of Surgeons in Ireland, Dublin 2, Ireland
| | - Viola Vogel
- Laboratory of Applied Mechanobiology, Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland
| |
Collapse
|
10
|
Principles and Methods in Computational Membrane Protein Design. J Mol Biol 2021; 433:167154. [PMID: 34271008 DOI: 10.1016/j.jmb.2021.167154] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2021] [Revised: 07/03/2021] [Accepted: 07/06/2021] [Indexed: 01/13/2023]
Abstract
After decades of progress in computational protein design, the design of proteins folding and functioning in lipid membranes appears today as the next frontier. Some notable successes in the de novo design of simplified model membrane protein systems have helped articulate fundamental principles of protein folding, architecture and interaction in the hydrophobic lipid environment. These principles are reviewed here, together with the computational methods and approaches that were used to identify them. We provide an overview of the methodological innovations in the generation of new protein structures and functions and in the development of membrane-specific energy functions. We highlight the opportunities offered by new machine learning approaches applied to protein design, and by new experimental characterization techniques applied to membrane proteins. Although membrane protein design is in its infancy, it appears more reachable than previously thought.
Collapse
|
11
|
Gorgun D, Lihan M, Kapoor K, Tajkhorshid E. Binding mode of SARS-CoV-2 fusion peptide to human cellular membrane. Biophys J 2021; 120:2914-2926. [PMID: 33675757 PMCID: PMC7929786 DOI: 10.1016/j.bpj.2021.02.041] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2020] [Revised: 01/29/2021] [Accepted: 02/15/2021] [Indexed: 01/08/2023] Open
Abstract
Infection of human cells by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics simulations that take advantage of the highly mobile membrane mimetic model to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas, each for 300 ns. In 73% of the simulations, the FP reaches a stable, membrane-bound configuration, in which the peptide deeply penetrated into the membrane. Clustering of the results reveals three major membrane-binding modes (binding modes 1-3), in which binding mode 1 populates over half of the data points. Taking into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, the significant depth of penetration of the whole peptide, and the dense population of the respective cluster, we propose that the most deeply inserted membrane-bound form (binding mode 1) represents more closely the biologically relevant form. Analysis of FP-lipid interactions shows the involvement of specific residues, previously described as the "fusion-active core residues," in membrane binding. Taken together, the results shed light on a key step involved in SARS-CoV2 infection, with potential implications in designing novel inhibitors.
Collapse
Affiliation(s)
- Defne Gorgun
- 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
| | - Muyun Lihan
- 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
| | - Karan Kapoor
- 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
| | - 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.
| |
Collapse
|
12
|
Protein Structure Prediction and Design in a Biologically Realistic Implicit Membrane. Biophys J 2020; 118:2042-2055. [PMID: 32224301 DOI: 10.1016/j.bpj.2020.03.006] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 02/04/2020] [Accepted: 03/09/2020] [Indexed: 11/19/2022] Open
Abstract
Protein design is a powerful tool for elucidating mechanisms of function and engineering new therapeutics and nanotechnologies. Although soluble protein design has advanced, membrane protein design remains challenging because of difficulties in modeling the lipid bilayer. In this work, we developed an implicit approach that captures the anisotropic structure, shape of water-filled pores, and nanoscale dimensions of membranes with different lipid compositions. The model improves performance in computational benchmarks against experimental targets, including prediction of protein orientations in the bilayer, ΔΔG calculations, native structure discrimination, and native sequence recovery. When applied to de novo protein design, this approach designs sequences with an amino acid distribution near the native amino acid distribution in membrane proteins, overcoming a critical flaw in previous membrane models that were prone to generating leucine-rich designs. Furthermore, the proteins designed in the new membrane model exhibit native-like features including interfacial aromatic side chains, hydrophobic lengths compatible with bilayer thickness, and polar pores. Our method advances high-resolution membrane protein structure prediction and design toward tackling key biological questions and engineering challenges.
Collapse
|
13
|
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
|
14
|
A network of phosphatidylinositol 4,5-bisphosphate binding sites regulates gating of the Ca 2+-activated Cl - channel ANO1 (TMEM16A). Proc Natl Acad Sci U S A 2019; 116:19952-19962. [PMID: 31515451 DOI: 10.1073/pnas.1904012116] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
ANO1 (TMEM16A) is a Ca2+-activated Cl- channel that regulates diverse cellular functions including fluid secretion, neuronal excitability, and smooth muscle contraction. ANO1 is activated by elevation of cytosolic Ca2+ and modulated by phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2]. Here, we describe a closely concerted experimental and computational study, including electrophysiology, mutagenesis, functional assays, and extended sampling of lipid-protein interactions with molecular dynamics (MD) to characterize PI(4,5)P2 binding modes and sites on ANO1. ANO1 currents in excised inside-out patches activated by 270 nM Ca2+ at +100 mV are increased by exogenous PI(4,5)P2 with an EC50 = 1.24 µM. The effect of PI(4,5)P2 is dependent on membrane voltage and Ca2+ and is explained by a stabilization of the ANO1 Ca2+-bound open state. Unbiased atomistic MD simulations with 1.4 mol% PI(4,5)P2 in a phosphatidylcholine bilayer identified 8 binding sites with significant probability of binding PI(4,5)P2 Three of these sites captured 85% of all ANO1-PI(4,5)P2 interactions. Mutagenesis of basic amino acids near the membrane-cytosol interface found 3 regions of ANO1 critical for PI(4,5)P2 regulation that correspond to the same 3 sites identified by MD. PI(4,5)P2 is stabilized by hydrogen bonding between amino acid side chains and phosphate/hydroxyl groups on PI(4,5)P2 Binding of PI(4,5)P2 alters the position of the cytoplasmic extension of TM6, which plays a crucial role in ANO1 channel gating, and increases the accessibility of the inner vestibule to Cl- ions. We propose a model consisting of a network of 3 PI(4,5)P2 binding sites at the cytoplasmic face of the membrane allosterically regulating ANO1 channel gating.
Collapse
|
15
|
Muller MP, Jiang T, Sun C, Lihan M, Pant S, Mahinthichaichan P, Trifan A, Tajkhorshid E. Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation. Chem Rev 2019; 119:6086-6161. [PMID: 30978005 PMCID: PMC6506392 DOI: 10.1021/acs.chemrev.8b00608] [Citation(s) in RCA: 133] [Impact Index Per Article: 26.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.
Collapse
Affiliation(s)
- Melanie P. Muller
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Tao Jiang
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Chang Sun
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Muyun Lihan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Shashank Pant
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Paween Mahinthichaichan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Anda Trifan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Emad Tajkhorshid
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
16
|
Vermaas JV, Bentley GJ, Beckham GT, Crowley MF. Membrane Permeability of Terpenoids Explored with Molecular Simulation. J Phys Chem B 2018; 122:10349-10361. [DOI: 10.1021/acs.jpcb.8b08688] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Josh V. Vermaas
- Biosciences Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Gayle J. Bentley
- National Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Gregg T. Beckham
- National Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Michael F. Crowley
- Biosciences Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| |
Collapse
|
17
|
Vermaas JV, Beckham GT, Crowley MF. Membrane Permeability of Fatty Acyl Compounds Studied via Molecular Simulation. J Phys Chem B 2017; 121:11311-11324. [PMID: 29040809 DOI: 10.1021/acs.jpcb.7b08233] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Interest in fatty acid-derived products as fuel and chemical precursors has grown substantially. Microbes can be genetically engineered to produce fatty acid-derived products that are able to cross host membranes and can be extracted into an applied organic overlay. This process is thought to be passive, with a rate dependent on the chemistry of the crossing compound. The relationship between the chemical composition and the energetics and kinetics of product accumulation within the overlay is not well understood. Through biased and unbiased molecular simulation, we compute the membrane permeability coefficients from production to extraction for different fatty acyl products, including fatty acids, fatty alcohols, fatty aldehydes, alkanes, and alkenes. These simulations identify specific interactions that accelerate the transit of aldehydes across the membrane bilayer relative to other oxidized products, specifically the lack of hydrogen bonds to the surrounding membrane environment. However, since extraction from the outer membrane leaflet into the organic phase is found to be rate limiting for the entire process, we find that fatty alcohols and fatty aldehydes would both manifest similar fluxes into a dodecane overlay under equivalent conditions, outpacing the accumulation of acids or alkanes into the organic phase. Since aldehydes are known to be highly reactive as well as toxic in high quantities, the findings suggest that indeed fatty alcohols are the optimal long-tail fatty acyl product for extraction.
Collapse
Affiliation(s)
- Josh V Vermaas
- Biosciences Center, National Renewable Energy Laboratory , Golden, Colorado 80401, United States
| | - Gregg T Beckham
- National Bioenergy Center, National Renewable Energy Laboratory , Golden, Colorado 80401, United States
| | - Michael F Crowley
- Biosciences Center, National Renewable Energy Laboratory , Golden, Colorado 80401, United States
| |
Collapse
|