1
|
Mohole M, Naglekar A, Sengupta D, Chattopadhyay A. Probing the energy landscape of the lipid interactions of the serotonin 1A receptor. Biophys Chem 2024; 313:107289. [PMID: 39002247 DOI: 10.1016/j.bpc.2024.107289] [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: 02/16/2024] [Revised: 07/02/2024] [Accepted: 07/02/2024] [Indexed: 07/15/2024]
Abstract
G protein-coupled receptors (GPCRs) are lipid-regulated transmembrane proteins that play a central role in cell signaling and pharmacology. Although the role of membrane lipids in GPCR function is well established, the underlying GPCR-lipid interactions have not been thermodynamically characterized due to the complexity of these interactions. In this work, we estimate the energetics and dynamics of lipid association from coarse-grain simulations of the serotonin1A receptor embedded in a complex membrane. We show that lipids bind to the receptor with varying energetics of 1-4 kT, and timescales of 1-10 μs. The most favorable energetics and longest residence times are observed for cholesterol, glycosphingolipid GM1, phosphatidylethanolamine (PE) and phosphatidylserine (PS) lipids. Multi-exponential fitting of the contact probability suggests distinct dynamic regimes, corresponding to ps, ns and μs timescales, that we correlate with the annular, intermediate and non-annular lipid sites. The timescales of lipid binding correspond to high barrier heights, despite their relatively weaker energetics. Our results highlight that GPCR-lipid interactions are driven by both thermodynamic interactions and the dynamical features of lipid binding.
Collapse
Affiliation(s)
- Madhura Mohole
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411 008, India; Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201 002, India
| | - Amit Naglekar
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411 008, India; Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201 002, India
| | - Durba Sengupta
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411 008, India; Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201 002, India.
| | - Amitabha Chattopadhyay
- CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad 500 007, India; Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201 002, India.
| |
Collapse
|
2
|
Zoccante A, Cara E, Ferrarese Lupi F, Hönicke P, Kayser Y, Beckhoff B, Klapetek P, Marchi D, Cossi M. The thermodynamics of self-assembled monolayer formation: a computational and experimental study of thiols on a flat gold surface. Phys Chem Chem Phys 2024; 26:18799-18807. [PMID: 38938190 DOI: 10.1039/d4cp01322k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024]
Abstract
A methodology based on molecular dynamics simulations is presented to determine the chemical potential of thiol self-assembled monolayers on a gold surface. The thiol de-solvation and then the monolayer formation are described by thermodynamic integration with a gradual decoupling of one molecule from the environment, with the necessary corrections to account for standard state changes. The procedure is applied both to physisorbed undissociated thiol molecules and to chemisorbed dissociated thiyl radicals, considering in the latter case the possible chemical potential of the produced hydrogen. We considered monolayers formed by either 7-mercapto-4-methylcoumarin (MMC) or 3-mercapto-propanoic acid (MPA) on a flat gold surface: the free energy profiles with respect to the monolayer density are consistent with a transition from a very stable lying-down phase at low densities to a standing-up phase at higher densities, as expected. The maximum densities of thermodynamically stable monolayers are compared to experimental measures performed with reference-free grazing-incidence X-ray fluorescence (RF-GIXRF) on the same systems, finding a better agreement in the case of chemisorbed thiyl radicals.
Collapse
Affiliation(s)
- Alberto Zoccante
- Dipartimento di Scienze e Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, I-15121 Alessandria, Italy.
| | - Eleonora Cara
- Istituto Nazionale di Ricerca Metrologica (INRiM), Strada delle Cacce, 91, 10135 Torino, Italy
| | - Federico Ferrarese Lupi
- Istituto Nazionale di Ricerca Metrologica (INRiM), Strada delle Cacce, 91, 10135 Torino, Italy
| | - Philipp Hönicke
- Physikalisch-Technische Bundesanstalt (PTB), Abbestr. 2-12, 10587 Berlin, Germany
- Helmholtz-Zentrum Berlin, Hahn-Meitner Platz 1, 14109 Berlin, Germany
| | - Yves Kayser
- Physikalisch-Technische Bundesanstalt (PTB), Abbestr. 2-12, 10587 Berlin, Germany
- Max-Planck-Institut für chemische Energiekonversion Mulheim an der Ruhr, Nordrhein-Westfalen, DE, Germany
| | - Burkhard Beckhoff
- Physikalisch-Technische Bundesanstalt (PTB), Abbestr. 2-12, 10587 Berlin, Germany
| | - Petr Klapetek
- Department of Nanometrology, Czech Metrology Institute (CMI), Okružní 31, 638 00 Brno, Czech Republic
| | - Davide Marchi
- Dipartimento di Scienze e Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, I-15121 Alessandria, Italy.
| | - Maurizio Cossi
- Dipartimento di Scienze e Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, I-15121 Alessandria, Italy.
| |
Collapse
|
3
|
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
|
4
|
Lagardère L, Maurin L, Adjoua O, El Hage K, Monmarché P, Piquemal JP, Hénin J. Lambda-ABF: Simplified, Portable, Accurate, and Cost-Effective Alchemical Free-Energy Computation. J Chem Theory Comput 2024; 20:4481-4498. [PMID: 38805379 DOI: 10.1021/acs.jctc.3c01249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
We introduce the lambda-Adaptive Biasing Force (lambda-ABF) method for the computation of alchemical free-energy differences. We propose a software implementation and showcase it on biomolecular systems. The method arises from coupling multiple-walker adaptive biasing force with λ-dynamics. The sampling of the alchemical variable is continuous and converges toward a uniform distribution, making manual optimization of the λ schedule unnecessary. Contrary to most other approaches, alchemical free-energy estimates are obtained immediately without any postprocessing. Free diffusion of λ improves orthogonal relaxation compared to fixed-λ thermodynamic integration or free-energy perturbation. Furthermore, multiple walkers provide generic orthogonal space coverage with minimal user input and negligible computational overhead. We show that our high-performance implementations coupling the Colvars library with NAMD and Tinker-HP can address real-world cases including ligand-receptor binding with both fixed-charge and polarizable models, with a demonstrably richer sampling than fixed-λ methods. The implementation is fully open-source, publicly available, and readily usable by practitioners of current alchemical methods. Thanks to the portable Colvars library, lambda-ABF presents a unified user interface regardless of the back-end (NAMD, Tinker-HP, or any software to be interfaced in the future), sparing users the effort of learning multiple interfaces. Finally, the Colvars Dashboard extension of the visual molecular dynamics (VMD) software provides an interactive monitoring and diagnostic tool for lambda-ABF simulations.
Collapse
Affiliation(s)
- Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Institut Parisien de Chimie Physique et Théorique, FR2622 CNRS, 75005 Paris, France
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Lise Maurin
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Laboratoire Jacques-Louis Lions, UMR 7589 CNRS, 75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
| | - Krystel El Hage
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Pierre Monmarché
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Laboratoire Jacques-Louis Lions, UMR 7589 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique, Université Paris Cité, CNRS, UPR 9080, 75005 Paris, France
| |
Collapse
|
5
|
Sandberg JW, Santiago-McRae E, Ennis J, Brannigan G. The density-threshold affinity: Calculating lipid binding affinities from unbiased coarse-grained molecular dynamics simulations. Methods Enzymol 2024; 701:47-82. [PMID: 39025580 DOI: 10.1016/bs.mie.2024.03.008] [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
Many membrane proteins are sensitive to their local lipid environment. As structural methods for membrane proteins have improved, there is growing evidence of direct, specific binding of lipids to protein surfaces. Unfortunately the workhorse of understanding protein-small molecule interactions, the binding affinity for a given site, is experimentally inaccessible for these systems. Coarse-grained molecular dynamics simulations can be used to bridge this gap, and are relatively straightforward to learn. Such simulations allow users to observe spontaneous binding of lipids to membrane proteins and quantify localized densities of individual lipids or lipid fragments. In this chapter we outline a protocol for extracting binding affinities from these localized distributions, known as the "density threshold affinity." The density threshold affinity uses an adaptive and flexible definition of site occupancy that alleviates the need to distinguish between "bound'' lipids and bulk lipids that are simply diffusing through the site. Furthermore, the method allows "bead-level" resolution that is suitable for the case where lipids share binding sites, and circumvents ambiguities about a relevant reference state. This approach provides a convenient and straightforward method for comparing affinities of a single lipid species for multiple sites, multiple lipids for a single site, and/or a single lipid species modeled using multiple forcefields.
Collapse
Affiliation(s)
- Jesse W Sandberg
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ, United States
| | - Ezry Santiago-McRae
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ, United States
| | - Jahmal Ennis
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ, United States
| | - Grace Brannigan
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ, United States; Department of Physics, Rutgers University, Camden, NJ, United States.
| |
Collapse
|
6
|
Fortea E, Lee S, Chadda R, Argyros Y, Sandal P, Mahoney-Kruszka R, Ciftci HD, Falzone ME, Huysmans G, Robertson JL, Boudker O, Accardi A. Structural basis of pH-dependent activation in a CLC transporter. Nat Struct Mol Biol 2024; 31:644-656. [PMID: 38279055 PMCID: PMC11262703 DOI: 10.1038/s41594-023-01210-5] [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] [Received: 08/16/2023] [Accepted: 12/22/2023] [Indexed: 01/28/2024]
Abstract
CLCs are dimeric chloride channels and anion/proton exchangers that regulate processes such as muscle contraction and endo-lysosome acidification. Common gating controls their activity; its closure simultaneously silences both protomers, and its opening allows them to independently transport ions. Mutations affecting common gating in human CLCs cause dominant genetic disorders. The structural rearrangements underlying common gating are unknown. Here, using single-particle cryo-electron microscopy, we show that the prototypical Escherichia coli CLC-ec1 undergoes large-scale rearrangements in activating conditions. The slow, pH-dependent remodeling of the dimer interface leads to the concerted opening of the intracellular H+ pathways and is required for transport. The more frequent formation of short water wires in the open H+ pathway enables Cl- pore openings. Mutations at disease-causing sites favor CLC-ec1 activation and accelerate common gate opening in the human CLC-7 exchanger. We suggest that the pH activation mechanism of CLC-ec1 is related to the common gating of CLC-7.
Collapse
Affiliation(s)
- Eva Fortea
- Department of Physiology and Biophysics, Weill Cornell Medical School, New York, NY, USA
- Department of Anesthesiology, Weill Cornell Medical School, New York, NY, USA
| | - Sangyun Lee
- Department of Anesthesiology, Weill Cornell Medical School, New York, NY, USA
| | - Rahul Chadda
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, USA
| | - Yiorgos Argyros
- Department of Anesthesiology, Weill Cornell Medical School, New York, NY, USA
- Department of Biochemistry, Weill Cornell Medical School, New York, NY, USA
| | - Priyanka Sandal
- Department of Molecular Physiology and Biophysics, The University of Iowa, Iowa City, IA, USA
| | - Robyn Mahoney-Kruszka
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, USA
| | - Hatice Didar Ciftci
- Department of Physiology and Biophysics, Weill Cornell Medical School, New York, NY, USA
- Tri-Institutional Training Program in Chemical Biology, New York, NY, USA
| | - Maria E Falzone
- Department of Anesthesiology, Weill Cornell Medical School, New York, NY, USA
- Department of Biochemistry, Weill Cornell Medical School, New York, NY, USA
| | - Gerard Huysmans
- Department of Physiology and Biophysics, Weill Cornell Medical School, New York, NY, USA
- Erasmus University, Jette, Belgium
| | - Janice L Robertson
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, USA
| | - Olga Boudker
- Department of Physiology and Biophysics, Weill Cornell Medical School, New York, NY, USA.
- Howard Hughes Medical Institute, Chevy Chase, MD, USA.
| | - Alessio Accardi
- Department of Physiology and Biophysics, Weill Cornell Medical School, New York, NY, USA.
- Department of Anesthesiology, Weill Cornell Medical School, New York, NY, USA.
- Department of Biochemistry, Weill Cornell Medical School, New York, NY, USA.
| |
Collapse
|
7
|
Joseph TT, Bu W, Haji-Ghassemi O, Chen YS, Woll K, Allen PD, Brannigan G, van Petegem F, Eckenhoff RG. Propofol directly binds and inhibits skeletal muscle ryanodine receptor 1 (RyR1). BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.10.575040. [PMID: 38260485 PMCID: PMC10802444 DOI: 10.1101/2024.01.10.575040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
As the primary Ca 2+ release channel in skeletal muscle sarcoplasmic reticulum (SR), mutations in the type 1 ryanodine receptor (RyR1) or its binding partners underlie a constellation of muscle disorders, including malignant hyperthermia (MH). In patients with MH mutations, exposure to triggering drugs such as the halogenated volatile anesthetics biases RyR1 to an open state, resulting in uncontrolled Ca 2+ release, sarcomere tension and heat production. Restoration of Ca 2+ into the SR also consumes ATP, generating a further untenable metabolic load. When anesthetizing patients with known MH mutations, the non-triggering intravenous general anesthetic propofol is commonly substituted for triggering anesthetics. Evidence of direct binding of anesthetic agents to RyR1 or its binding partners is scant, and the atomic-level interactions of propofol with RyR1 are entirely unknown. Here, we show that propofol decreases RyR1 opening in heavy SR vesicles and planar lipid bilayers, and that it inhibits activator-induced Ca 2+ release from SR in human skeletal muscle. In addition to confirming direct binding, photoaffinity labeling using m- azipropofol (AziP m ) revealed several putative propofol binding sites on RyR1. Prediction of binding affinity by molecular dynamics simulation suggests that propofol binds at least one of these sites at clinical concentrations. These findings invite the hypothesis that in addition to propofol not triggering MH, it may also be protective against MH by inhibiting induced Ca 2+ flux through RyR1.
Collapse
|
8
|
Naglekar A, Chattopadhyay A, Sengupta D. Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment. J Phys Chem B 2023; 127:11000-11010. [PMID: 38111968 DOI: 10.1021/acs.jpcb.3c05930] [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: 12/20/2023]
Abstract
The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.
Collapse
Affiliation(s)
- Amit Naglekar
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Amitabha Chattopadhyay
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
- CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad 500007, India
| | - Durba Sengupta
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| |
Collapse
|
9
|
Fu H, Chipot C, Shao X, Cai W. Standard Binding Free-Energy Calculations: How Far Are We from Automation? J Phys Chem B 2023; 127:10459-10468. [PMID: 37824848 DOI: 10.1021/acs.jpcb.3c04370] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2023]
Abstract
Recent success stories suggest that in silico protein-ligand binding free-energy calculations are approaching chemical accuracy. However, their widespread application remains limited by the extensive human intervention required, posing challenges for the neophyte. As such, it is critical to develop automated workflows for estimating protein-ligand binding affinities with minimum personal involvement. Key human efforts include setting up and tuning enhanced-sampling or alchemical-transformation algorithms as a preamble to computational binding free-energy estimations. Additionally, preparing input files, bookkeeping, and postprocessing represent nontrivial tasks. In this Perspective, we discuss recent progress in automating standard binding free-energy calculations, featuring the development of adaptive or parameter-free algorithms, standardization of binding free-energy calculation workflows, and the implementation of user-friendly software. We also assess the current state of automated standard binding free-energy calculations and evaluate the limitations of existing methods. Last, we outline the requirements for future algorithms and workflows to facilitate automated free-energy calculations for diverse protein-ligand complexes.
Collapse
Affiliation(s)
- Haohao Fu
- State Key Laboratory of Medicinal Chemical Biology, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| | - Christophe Chipot
- Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign, UMR no. 7019, Université de Lorraine, BP 70239, F-54506 Vandoeuvre-lès-Nancy, France
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
- Department of Chemistry, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
- Department of Chemistry, The University of Hawai'i at Ma̅noa, 2545 McCarthy Mall, Honolulu, Hawaii 96822, United States
| | - Xueguang Shao
- State Key Laboratory of Medicinal Chemical Biology, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| | - Wensheng Cai
- State Key Laboratory of Medicinal Chemical Biology, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| |
Collapse
|
10
|
Jiang W, Wijerathne TD, Zhang H, Lin YC, Jo S, Im W, Lacroix JJ, Luo YL. Structural and thermodynamic framework for PIEZO1 modulation by small molecules. Proc Natl Acad Sci U S A 2023; 120:e2310933120. [PMID: 38060566 PMCID: PMC10723123 DOI: 10.1073/pnas.2310933120] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Accepted: 10/12/2023] [Indexed: 12/17/2023] Open
Abstract
Mechanosensitive PIEZO channels constitute potential pharmacological targets for multiple clinical conditions, spurring the search for potent chemical PIEZO modulators. Among them is Yoda1, a widely used synthetic small molecule PIEZO1 activator discovered through cell-based high-throughput screening. Yoda1 is thought to bind to PIEZO1's mechanosensory arm domain, sandwiched between two transmembrane regions near the channel pore. However, how the binding of Yoda1 to this region promotes channel activation remains elusive. Here, we first demonstrate that cross-linking PIEZO1 repeats A and B with disulfide bridges reduces the effects of Yoda1 in a redox-dependent manner, suggesting that Yoda1 acts by perturbing the contact between these repeats. Using molecular dynamics-based absolute binding free energy simulations, we next show that Yoda1 preferentially occupies a deeper, amphipathic binding site with higher affinity in PIEZO1 open state. Using Yoda1's binding poses in open and closed states, relative binding free energy simulations were conducted in the membrane environment, recapitulating structure-activity relationships of known Yoda1 analogs. Through virtual screening of an 8 million-compound library using computed fragment maps of the Yoda1 binding site, we subsequently identified two chemical scaffolds with agonist activity toward PIEZO1. This study supports a pharmacological model in which Yoda1 activates PIEZO1 by wedging repeats A and B, providing a structural and thermodynamic framework for the rational design of PIEZO1 modulators. Beyond PIEZO channels, the three orthogonal computational approaches employed here represent a promising path toward drug discovery in highly heterogeneous membrane protein systems.
Collapse
Affiliation(s)
- Wenjuan Jiang
- Department of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, CA91766
| | - Tharaka D. Wijerathne
- Department of Basic Medical Sciences, Western University of Health Sciences, Pomona, CA91766
| | - Han Zhang
- Department of Biological Sciences, Lehigh University, Bethlehem, PA18015
- Department of Chemistry, Lehigh University, Bethlehem, PA18015
- Department of Bioengineering, Lehigh University, Bethlehem, PA18015
- Department of Computer Science and Engineering, Lehigh University, Bethlehem, PA18015
| | - Yi-Chun Lin
- Department of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, CA91766
| | - Sunhwan Jo
- Department of Pharmaceutical Sciences, University of Maryland, Baltimore, MD21201
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, PA18015
- Department of Chemistry, Lehigh University, Bethlehem, PA18015
- Department of Bioengineering, Lehigh University, Bethlehem, PA18015
- Department of Computer Science and Engineering, Lehigh University, Bethlehem, PA18015
| | - Jerome J. Lacroix
- Department of Basic Medical Sciences, Western University of Health Sciences, Pomona, CA91766
| | - Yun L. Luo
- Department of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, CA91766
| |
Collapse
|
11
|
Papadourakis M, Sinenka H, Matricon P, Hénin J, Brannigan G, Pérez-Benito L, Pande V, van Vlijmen H, de Graaf C, Deflorian F, Tresadern G, Cecchini M, Cournia Z. Alchemical Free Energy Calculations on Membrane-Associated Proteins. J Chem Theory Comput 2023; 19:7437-7458. [PMID: 37902715 PMCID: PMC11017255 DOI: 10.1021/acs.jctc.3c00365] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Indexed: 10/31/2023]
Abstract
Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.
Collapse
Affiliation(s)
- Michail Papadourakis
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Hryhory Sinenka
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Pierre Matricon
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Jérôme Hénin
- Laboratoire
de Biochimie Théorique UPR 9080, CNRS and Université Paris Cité, 75005 Paris, France
| | - Grace Brannigan
- Center
for Computational and Integrative Biology, Rutgers University−Camden, Camden, New Jersey 08103, United States of America
- Department
of Physics, Rutgers University−Camden, Camden, New Jersey 08102, United States
of America
| | - Laura Pérez-Benito
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Vineet Pande
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Herman van Vlijmen
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Chris de Graaf
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Francesca Deflorian
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Gary Tresadern
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Marco Cecchini
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Zoe Cournia
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| |
Collapse
|
12
|
Ansell TB, Corey RA, Viti LV, Kinnebrew M, Rohatgi R, Siebold C, Sansom MS. The energetics and ion coupling of cholesterol transport through Patched1. SCIENCE ADVANCES 2023; 9:eadh1609. [PMID: 37611095 PMCID: PMC10446486 DOI: 10.1126/sciadv.adh1609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 07/24/2023] [Indexed: 08/25/2023]
Abstract
Patched1 (PTCH1) is a tumor suppressor protein of the mammalian Hedgehog (HH) signaling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating ciliary cholesterol accessibility. Using extensive molecular dynamics simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15 to 20 kilojoule per mole for cholesterol export. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between cation binding sites, transmembrane motions, and PTCH1 activity. Using complementary simulations of Dispatched1, we find that transition between "inward-open" and solvent "occluded" states is accompanied by Na+-induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion coupling stoichiometries of PTCH1 transport mechanisms, whereby one to three Na+ or two to three K+ couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.
Collapse
Affiliation(s)
- T. Bertie Ansell
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | - Robin A. Corey
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
- School of Physiology, Pharmacology and Neuroscience, Bristol University, Bristol BS8 1TD, UK
| | - Lucrezia Vittoria Viti
- Division of Structural Biology, Wellcome Centre for Human Genetics, Roosevelt Drive, Oxford OX3 7BN, UK
| | - Maia Kinnebrew
- Departments of Biochemistry and Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Rajat Rohatgi
- Departments of Biochemistry and Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Christian Siebold
- Division of Structural Biology, Wellcome Centre for Human Genetics, Roosevelt Drive, Oxford OX3 7BN, UK
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| |
Collapse
|
13
|
Languin-Cattoën O, Sterpone F, Stirnemann G. Binding site plasticity regulation of the FimH catch-bond mechanism. Biophys J 2023; 122:2744-2756. [PMID: 37264571 PMCID: PMC10397818 DOI: 10.1016/j.bpj.2023.05.029] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Revised: 04/05/2023] [Accepted: 05/26/2023] [Indexed: 06/03/2023] Open
Abstract
The bacterial fimbrial adhesin FimH is a remarkable and well-studied catch-bond protein found at the tip of E. coli type 1 pili, which allows pathogenic strains involved in urinary tract infections to bind high-mannose glycans exposed on human epithelia. The catch-bond behavior of FimH, where the strength of the interaction increases when a force is applied to separate the two partners, enables the bacteria to resist clearance when they are subjected to shear forces induced by urine flow. Two decades of experimental studies performed at the single-molecule level, as well as x-ray crystallography and modeling studies, have led to a consensus picture whereby force separates the binding domain from an inhibitor domain, effectively triggering an allosteric conformational change in the former. This force-induced allostery is thought to be responsible for an increased binding affinity at the core of the catch-bond mechanism. However, some important questions remain, the most challenging one being that the crystal structures corresponding to these two allosteric states show almost superimposable binding site geometries, which questions the molecular origin for the large difference in affinity. Using molecular dynamics with a combination of enhanced-sampling techniques, we demonstrate that the static picture provided by the crystal structures conceals a variety of binding site conformations that have a key impact on the apparent affinity. Crucially, the respective populations in each of these conformations are very different between the two allosteric states of the binding domain, which can then be related to experimental affinity measurements. We also evidence a previously unappreciated but important effect: in addition to the well-established role of the force as an allosteric regulator via domain separation, application of force tends to directly favor the high-affinity binding site conformations. We hypothesize that this additional "local" catch-bond effect could delay unbinding between the bacteria and the host cell before the "global" allosteric transition occurs, as well as stabilizing the complex even more once in the high-affinity allosteric state.
Collapse
Affiliation(s)
- Olivier Languin-Cattoën
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, Université Paris Cité, PSL University, Paris, France
| | - Fabio Sterpone
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, Université Paris Cité, PSL University, Paris, France.
| | - Guillaume Stirnemann
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, Université Paris Cité, PSL University, Paris, France.
| |
Collapse
|
14
|
Clark F, Robb G, Cole DJ, Michel J. Comparison of Receptor-Ligand Restraint Schemes for Alchemical Absolute Binding Free Energy Calculations. J Chem Theory Comput 2023; 19:3686-3704. [PMID: 37285579 PMCID: PMC10308817 DOI: 10.1021/acs.jctc.3c00139] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2023] [Indexed: 06/09/2023]
Abstract
Alchemical absolute binding free energy calculations are of increasing interest in drug discovery. These calculations require restraints between the receptor and ligand to restrict their relative positions and, optionally, orientations. Boresch restraints are commonly used, but they must be carefully selected in order to sufficiently restrain the ligand and to avoid inherent instabilities. Applying multiple distance restraints between anchor points in the receptor and ligand provides an alternative framework without inherent instabilities which may provide convergence benefits by more strongly restricting the relative movements of the receptor and ligand. However, there is no simple method to calculate the free energy of releasing these restraints due to the coupling of the internal and external degrees of freedom of the receptor and ligand. Here, a method to rigorously calculate free energies of binding with multiple distance restraints by imposing intramolecular restraints on the anchor points is proposed. Absolute binding free energies for the human macrophage migration inhibitory factor/MIF180, system obtained using a variety of Boresch restraints and rigorous and nonrigorous implementations of multiple distance restraints are compared. It is shown that several multiple distance restraint schemes produce estimates in good agreement with Boresch restraints. In contrast, calculations without orientational restraints produce erroneously favorable free energies of binding by up to approximately 4 kcal mol-1. These approaches offer new options for the deployment of alchemical absolute binding free energy calculations.
Collapse
Affiliation(s)
- Finlay Clark
- EaStCHEM
School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Graeme Robb
- Oncology
R&D, AstraZeneca, Cambridge CB4 0WG, United Kingdom
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon Tyne NE1 7RU, United Kingdom
| | - Julien Michel
- EaStCHEM
School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| |
Collapse
|
15
|
Treptow W. Allosteric Modulation of Membrane Proteins by Small Low-Affinity Ligands. J Chem Inf Model 2023; 63:2047-2057. [PMID: 36933226 DOI: 10.1021/acs.jcim.2c01542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023]
Abstract
Membrane proteins may respond to a variety of ligands under an applied external stimulus. These ligands include small low-affinity molecules that account for functional effects in the mM range. Understanding the modulation of protein function by low-affinity ligands requires characterizing their atomic-level interactions under dilution, challenging the current resolution of theoretical and experimental routines. Part of the problem derives from the fact that small low-affinity ligands may interact with multiple sites of a membrane protein in a highly degenerate manner to a degree that it is better conceived as a partition phenomenon, hard to track at the molecular interface of the protein. Looking for new developments in the field, we rely on the classic two-state Boltzmann model to devise a novel theoretical description of the allosteric modulation mechanism of membrane proteins in the presence of small low-affinity ligands and external stimuli. Free energy stability of the partition process and its energetic influence on the protein coupling with the external stimulus are quantified. The outcome is a simple formulation that allows the description of the equilibrium shifts of the protein in terms of the grand-canonical partition function of the ligand at dilute concentrations. The model's predictions of the spatial distribution and response probability shift across a variety of ligand concentrations, and thermodynamic conjugates can be directly compared to macroscopic measurements, making it especially useful to interpret experimental data at the atomic level. Illustration and discussion of the theory is shown in the context of general anesthetics and voltage-gated channels for which structural data are available.
Collapse
Affiliation(s)
- Werner Treptow
- Laboratório de Biologia Teórica e Computacional (LBTC), Universidade de Brasília, Distrito Federal, Brasília CEP 70904-970, Brasil
| |
Collapse
|
16
|
Ansell TB, Corey RA, Viti LV, Kinnebrew M, Rohatgi R, Siebold C, Sansom MSP. The Energetics and Ion Coupling of Cholesterol Transport Through Patched1. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.14.528445. [PMID: 36824746 PMCID: PMC9949057 DOI: 10.1101/2023.02.14.528445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Patched1 (PTCH1) is the principal tumour suppressor protein of the mammalian Hedgehog (HH) signalling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the Class F G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating accessible cholesterol levels within ciliary membranes. Using extensive molecular dynamics (MD) simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15-20 kJ mol -1 for cholesterol export. In simulations we identify cation binding sites within the PTCH1 transmembrane domain (TMD) which may provide the energetic impetus for cholesterol transport. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between transmembrane motions and PTCH1 activity. Using complementary simulations of Dispatched1 (DISP1) we find that transition between 'inward-open' and solvent 'occluded' states is accompanied by Na + induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion-coupling stoichiometries of PTCH1 transport mechanisms, whereby 1-3 Na + or 2-3 K + couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.
Collapse
Affiliation(s)
- T. Bertie Ansell
- Department of Biochemistry, South Parks Road, Oxford, OX1 3QU, UK
| | - Robin A. Corey
- Department of Biochemistry, South Parks Road, Oxford, OX1 3QU, UK
| | - Lucrezia Vittoria Viti
- Division of Structural Biology, Wellcome Centre for Human Genetics, Roosevelt Drive, Oxford, OX3 7BN, UK
| | - Maia Kinnebrew
- Departments of Biochemistry and Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Rajat Rohatgi
- Departments of Biochemistry and Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Christian Siebold
- Division of Structural Biology, Wellcome Centre for Human Genetics, Roosevelt Drive, Oxford, OX3 7BN, UK
| | | |
Collapse
|
17
|
Levental I, Lyman E. Regulation of membrane protein structure and function by their lipid nano-environment. Nat Rev Mol Cell Biol 2023; 24:107-122. [PMID: 36056103 PMCID: PMC9892264 DOI: 10.1038/s41580-022-00524-4] [Citation(s) in RCA: 139] [Impact Index Per Article: 139.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/19/2022] [Indexed: 02/04/2023]
Abstract
Transmembrane proteins comprise ~30% of the mammalian proteome, mediating metabolism, signalling, transport and many other functions required for cellular life. The microenvironment of integral membrane proteins (IMPs) is intrinsically different from that of cytoplasmic proteins, with IMPs solvated by a compositionally and biophysically complex lipid matrix. These solvating lipids affect protein structure and function in a variety of ways, from stereospecific, high-affinity protein-lipid interactions to modulation by bulk membrane properties. Specific examples of functional modulation of IMPs by their solvating membranes have been reported for various transporters, channels and signal receptors; however, generalizable mechanistic principles governing IMP regulation by lipid environments are neither widely appreciated nor completely understood. Here, we review recent insights into the inter-relationships between complex lipidomes of mammalian membranes, the membrane physicochemical properties resulting from such lipid collectives, and the regulation of IMPs by either or both. The recent proliferation of high-resolution methods to study such lipid-protein interactions has led to generalizable insights, which we synthesize into a general framework termed the 'functional paralipidome' to understand the mutual regulation between membrane proteins and their surrounding lipid microenvironments.
Collapse
Affiliation(s)
- Ilya Levental
- Department of Molecular Physiology and Biological Physics, Center for Molecular and Cell Physiology, University of Virginia, Charlottesville, VA, USA.
| | - Ed Lyman
- Department of Physics and Astronomy, Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, USA.
| |
Collapse
|
18
|
Jiang W, Lacroix J, Luo YL. Importance of molecular dynamics equilibrium protocol on protein-lipid interaction near channel pore. BIOPHYSICAL REPORTS 2022; 2:100080. [PMID: 36425669 PMCID: PMC9680783 DOI: 10.1016/j.bpr.2022.100080] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Accepted: 09/27/2022] [Indexed: 06/16/2023]
Abstract
Multiscale molecular dynamics simulations using Martini coarse-grained (CG) and all-atom (AA) force fields are commonly used in membrane protein studies. In particular, reverse mapping an equilibrated CG model to an AA model offers an efficient way for preparing large membrane protein systems with complex protein shapes and lipid compositions. Here, we report that this hybrid CG-equilibrium-AA-production protocol may artificially increase lipid density and decrease hydration in ion channel pores walled with transmembrane gaps. To understand the origin of this conundrum, we conducted replicas of CG, AA, and CG reverse-mapped AA simulations of the pore domain of the mechanosensitive Piezo1 channel in a nonconducting conformation. Lipid/water density analysis and free energy calculations reveal that the lack of initial pore hydration allows excessive lipids to enter the upper pore lumen through gaps between pore helices during CG simulation. Due to the mismatch between CG and AA lipid kinetics, these pore lipids remain trapped in the subsequent AA simulations, despite unfavorable binding free energy. We tested several CG equilibrium protocols and found that a protocol restraining the whole lipid produces pore hydration consistent with AA results, thus eliminating this artifact for further studies of lipid gating and protein-lipid interactions.
Collapse
Affiliation(s)
- Wenjuan Jiang
- College of Pharmacy, Western University of Health Sciences, Pomona, California
| | - Jerome Lacroix
- Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, California
| | - Yun Lyna Luo
- College of Pharmacy, Western University of Health Sciences, Pomona, California
| |
Collapse
|
19
|
Petroff JT, Dietzen NM, Santiago-McRae E, Deng B, Washington MS, Chen LJ, Trent Moreland K, Deng Z, Rau M, Fitzpatrick JAJ, Yuan P, Joseph TT, Hénin J, Brannigan G, Cheng WWL. Open-channel structure of a pentameric ligand-gated ion channel reveals a mechanism of leaflet-specific phospholipid modulation. Nat Commun 2022; 13:7017. [PMID: 36385237 PMCID: PMC9668969 DOI: 10.1038/s41467-022-34813-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Accepted: 11/08/2022] [Indexed: 11/18/2022] Open
Abstract
Pentameric ligand-gated ion channels (pLGICs) mediate synaptic transmission and are sensitive to their lipid environment. The mechanism of phospholipid modulation of any pLGIC is not well understood. We demonstrate that the model pLGIC, ELIC (Erwinia ligand-gated ion channel), is positively modulated by the anionic phospholipid, phosphatidylglycerol, from the outer leaflet of the membrane. To explore the mechanism of phosphatidylglycerol modulation, we determine a structure of ELIC in an open-channel conformation. The structure shows a bound phospholipid in an outer leaflet site, and structural changes in the phospholipid binding site unique to the open-channel. In combination with streamlined alchemical free energy perturbation calculations and functional measurements in asymmetric liposomes, the data support a mechanism by which an anionic phospholipid stabilizes the activated, open-channel state of a pLGIC by specific, state-dependent binding to this site.
Collapse
Affiliation(s)
- John T Petroff
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA
| | - Noah M Dietzen
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA
| | - Ezry Santiago-McRae
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ, USA
| | - Brett Deng
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA
| | - Maya S Washington
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA
| | - Lawrence J Chen
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA
| | - K Trent Moreland
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA
| | - Zengqin Deng
- Department of Cell Biology and Physiology, Washington University School of Medicine, Saint Louis, MO, USA
- Center for the Investigation of Membrane Excitability Diseases, Washington University School of Medicine, Saint Louis, MO, USA
| | - Michael Rau
- Center for Cellular Imaging, Washington University School of Medicine, Saint Louis, MO, USA
| | - James A J Fitzpatrick
- Department of Cell Biology and Physiology, Washington University School of Medicine, Saint Louis, MO, USA
- Center for Cellular Imaging, Washington University School of Medicine, Saint Louis, MO, USA
- Department of Neuroscience, Washington University School of Medicine, Saint Louis, MO, USA
- Department of Biomedical Engineering, Washington University School of Medicine, Saint Louis, MO, USA
| | - Peng Yuan
- Department of Cell Biology and Physiology, Washington University School of Medicine, Saint Louis, MO, USA
- Center for the Investigation of Membrane Excitability Diseases, Washington University School of Medicine, Saint Louis, MO, USA
| | - Thomas T Joseph
- Department of Anesthesiology and Critical Care, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique, Université Paris Cité, CNRS UPR 9080, Paris, France
| | - Grace Brannigan
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ, USA
- Department of Physics, Rutgers University, Camden, NJ, USA
| | - Wayland W L Cheng
- Department of Anesthesiology, Washington University School of Medicine, Saint Louis, MO, USA.
| |
Collapse
|
20
|
Cirqueira L, Stock L, Treptow W. Concentration-Dependent Thermodynamic Analysis of the Partition Process of Small Ligands into Proteins. Comput Struct Biotechnol J 2022; 20:4885-4891. [PMID: 36147679 PMCID: PMC9468351 DOI: 10.1016/j.csbj.2022.08.049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2022] [Revised: 08/19/2022] [Accepted: 08/20/2022] [Indexed: 11/28/2022] Open
Abstract
In the category of functional low-affinity interactions, small ligands may interact with multiple protein sites in a highly degenerate manner. Better conceived as a partition phenomenon at the molecular interface of proteins, such low-affinity interactions appear to be hidden to our current experimental resolution making their structural and functional characterization difficult in the low concentration regime of physiological processes. Characterization of the partition phenomenon under higher chemical forces could be a relevant strategy to tackle the problem provided the results can be scaled back to the low concentration range. Far from being trivial, such scaling demands a concentration-dependent understanding of self-interactions of the ligands, structural perturbations of the protein, among other molecular effects. Accordingly, we elaborate a novel and detailed concentration-dependent thermodynamic analysis of the partition process of small ligands aiming at characterizing the stability and structure of the dilute phenomenon from high concentrations. In analogy to an “aggregate” binding constant of a small molecule over multiple sites of a protein receptor, the model defines the stability of the process as a macroscopic equilibrium constant for the partition number of ligands that can be used to analyze biochemical and functional data of two-component systems driven by low-affinity interactions. Acquisition of such modeling-based structural information is expected to be highly welcome by revealing more traceable protein-binding spots for non-specific ligands.
Collapse
|
21
|
Puyo-Fourtine J, Juillé M, Hénin J, Clavaguéra C, Duboué-Dijon E. Consistent Picture of Phosphate-Divalent Cation Binding from Models with Implicit and Explicit Electronic Polarization. J Phys Chem B 2022; 126:4022-4034. [PMID: 35608554 DOI: 10.1021/acs.jpcb.2c01158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The binding of divalent cations to the ubiquitous phosphate group is essential for a number of key biological processes, such as DNA compaction, RNA folding, or interactions of some proteins with membranes. Yet, probing their binding sites, modes, and associated binding free energy is a challenge for both experiments and simulations. In simulations, standard force fields strongly overestimate the interaction between phosphate groups and divalent cations. Here, we examine how different strategies to include electronic polarization effects in force fields─implicitly, through the use of scaled charges or pair-specific Lennard-Jones parameters, or explicitly, with the polarizable force fields Drude and AMOEBA─capture the interactions of a model phosphate compound, dimethyl phosphate, with calcium and magnesium divalent cations. We show that both implicit and explicit approaches, when carefully parameterized, are successful in capturing the overall binding free energy and that common trends emerge from the comparison of different simulation approaches. Overall, the binding is very moderate, slightly weaker for Ca2+ than Mg2+, and the solvent-shared ion pair is slightly more stable than the contact monodentate ion pair. The bidentate ion pair is higher in energy (or even fully unstable for Mg2+). Our results thus suggest practical ways to capture the divalent cations with biomolecular phosphate groups in complex biochemical systems. In particular, the computational efficiency of implicit models makes them ideally suited for large-scale simulations of biological assemblies, with improved accuracy compared to state-of-the-art fixed-charge force fields.
Collapse
Affiliation(s)
- Julie Puyo-Fourtine
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Marie Juillé
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Jérôme Hénin
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Elise Duboué-Dijon
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| |
Collapse
|
22
|
Ebrahimi M, Hénin J. Symmetry-Adapted Restraints for Binding Free Energy Calculations. J Chem Theory Comput 2022; 18:2494-2502. [PMID: 35230113 DOI: 10.1021/acs.jctc.1c01235] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Binding free energy calculations rely critically on a precise definition of the bound state and well-designed ligand restraints to ensure that binding free energy calculations converge rapidly and yield estimates of well-defined thermodynamic quantities. The distance-to-bound-configuration (DBC) is a single variable that can precisely delineate the bound state of a ligand including translational, rotational and conformational degrees of freedom and has been successfully used to capture binding modes with complex geometries. DBC is defined as the root-mean-square deviation (RMSD) of ligand coordinates in the frame of reference of the binding site. In the special case where the ligand features symmetry-equivalent atoms, a standard RMSD arbitrarily distinguishes equivalent poses, mixing equivalent and nonequivalent degrees of freedom, and preventing the precise delineation of the bound state ensemble, which negates the benefits of defining a flat-bottom binding restraint. To remedy this, we introduce a symmetry-adapted DBC coordinate where the RMSD is minimized over permutations of equivalent ligand atoms. This coordinate is implemented in a portable software library, the Collective Variables Module. We tested the approach by computing the absolute binding free energy of benzene to the engineered site of a mutant lysozyme (L99A/M102H) using alchemical free energy perturbation. We found that the symmetry-adapted restraint leads to well-behaved convergence of both the decoupling free energy in the binding site and the restrained free energy in the gas phase, recovering the affinity computed using a classic center-of-mass restraint. Thus, symmetry-adapted DBC seamlessly generalizes the benefits of DBC restraints to the case of symmetric ligands. The underlying symmetric RMSD coordinate can also be used for analyzing or biasing simulations in other contexts than affinity predictions.
Collapse
Affiliation(s)
- Mina Ebrahimi
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran 1417935840, Iran.,CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris 75005, France.,Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris 75005, France
| | - Jérôme Hénin
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris 75005, France.,Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris 75005, France
| |
Collapse
|
23
|
Hénin J, Lopes LJS, Fiorin G. Human Learning for Molecular Simulations: The Collective Variables Dashboard in VMD. J Chem Theory Comput 2022; 18:1945-1956. [PMID: 35143194 DOI: 10.1021/acs.jctc.1c01081] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Collective Variables Dashboard is a software tool for real-time, seamless exploration of molecular structures and trajectories in a customizable space of collective variables. The Dashboard arises from the integration of the Collective Variables Module (also known as Colvars) with the visualization software VMD, augmented with a fully discoverable graphical interface offering interactive workflows for the design and analysis of collective variables. Typical use cases include a priori design of collective variables for enhanced sampling and free energy simulations as well as analysis of any type of simulation or collection of structures in a collective variable space. A combination of those cases commonly occurs when preliminary simulations, biased or unbiased, reveal that an optimized set of collective variables is necessary to improve sampling in further simulations. Then the Dashboard provides an efficient way to intuitively explore the space of likely collective variables, validate them on existing data, and use the resulting collective variable definitions directly in further biased simulations using the Collective Variables Module. Visualization of biasing energies and forces is proposed to help analyze or plan biased simulations. We illustrate the use of the Dashboard on two applications: discovering coordinates to describe ligand unbinding from a protein binding site and designing volume-based variables to bias the hydration of a transmembrane pore.
Collapse
Affiliation(s)
- Jérôme Hénin
- Laboratoire de Biochimie Théorique UPR 9080, CNRS, Université de Paris, 75005 Paris, France.,Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Laura J S Lopes
- Theoretical Division T-1, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Giacomo Fiorin
- National Institute of Neurological Disorders and Stroke (NINDS) and National Heart, Lung and Blood Institute (NHLBI), Bethesda, Maryland 20892, United States
| |
Collapse
|
24
|
Song W, Corey RA, Ansell TB, Cassidy CK, Horrell MR, Duncan AL, Stansfeld PJ, Sansom MSP. PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:1188-1201. [PMID: 35020380 PMCID: PMC8830038 DOI: 10.1021/acs.jctc.1c00708] [Citation(s) in RCA: 60] [Impact Index Per Article: 30.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Indexed: 12/11/2022]
Abstract
Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.
Collapse
Affiliation(s)
- Wanling Song
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
- Rahko,
Clifton House, 46 Clifton
Terrace, Finsbury Park, London N4 3JP, United Kingdom
| | - Robin A. Corey
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - T. Bertie Ansell
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - C. Keith Cassidy
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Michael R. Horrell
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Anna L. Duncan
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Phillip J. Stansfeld
- School
of Life Sciences & Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Mark S. P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| |
Collapse
|
25
|
Cheng WWL, Arcario MJ, Petroff JT. Druggable Lipid Binding Sites in Pentameric Ligand-Gated Ion Channels and Transient Receptor Potential Channels. Front Physiol 2022; 12:798102. [PMID: 35069257 PMCID: PMC8777383 DOI: 10.3389/fphys.2021.798102] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Accepted: 12/02/2021] [Indexed: 12/17/2022] Open
Abstract
Lipids modulate the function of many ion channels, possibly through direct lipid-protein interactions. The recent outpouring of ion channel structures by cryo-EM has revealed many lipid binding sites. Whether these sites mediate lipid modulation of ion channel function is not firmly established in most cases. However, it is intriguing that many of these lipid binding sites are also known sites for other allosteric modulators or drugs, supporting the notion that lipids act as endogenous allosteric modulators through these sites. Here, we review such lipid-drug binding sites, focusing on pentameric ligand-gated ion channels and transient receptor potential channels. Notable examples include sites for phospholipids and sterols that are shared by anesthetics and vanilloids. We discuss some implications of lipid binding at these sites including the possibility that lipids can alter drug potency or that understanding protein-lipid interactions can guide drug design. Structures are only the first step toward understanding the mechanism of lipid modulation at these sites. Looking forward, we identify knowledge gaps in the field and approaches to address them. These include defining the effects of lipids on channel function in reconstituted systems using asymmetric membranes and measuring lipid binding affinities at specific sites using native mass spectrometry, fluorescence binding assays, and computational approaches.
Collapse
Affiliation(s)
- Wayland W L Cheng
- Department of Anesthesiology, Washington University in St. Louis, St. Louis, MO, United States
| | - Mark J Arcario
- Department of Anesthesiology, Washington University in St. Louis, St. Louis, MO, United States
| | - John T Petroff
- Department of Anesthesiology, Washington University in St. Louis, St. Louis, MO, United States
| |
Collapse
|
26
|
Duboué-Dijon E, Hénin J. Building intuition for binding free energy calculations: Bound state definition, restraints, and symmetry. J Chem Phys 2021; 154:204101. [PMID: 34241173 DOI: 10.1063/5.0046853] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The theory behind computation of absolute binding free energies using explicit-solvent molecular simulations is well-established, yet somewhat complex, with counter-intuitive aspects. This leads to frequent frustration, common misconceptions, and sometimes erroneous numerical treatment. To improve this, we present the main practically relevant segments of the theory with constant reference to physical intuition. We pinpoint the role of the implicit or explicit definition of the bound state (or the binding site) to make a robust link between an experimental measurement and a computational result. We clarify the role of symmetry and discuss cases where symmetry number corrections have been misinterpreted. In particular, we argue that symmetry corrections as classically presented are a source of confusion and could be advantageously replaced by restraint free energy contributions. We establish that contrary to a common intuition, partial or missing sampling of some modes of symmetric bound states does not affect the calculated decoupling free energies. Finally, we review these questions and pitfalls in the context of a few common practical situations: binding to a symmetric receptor (equivalent binding sites), binding of a symmetric ligand (equivalent poses), and formation of a symmetric complex, in the case of homodimerization.
Collapse
Affiliation(s)
- E Duboué-Dijon
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - J Hénin
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
27
|
Sharp L, Brannigan G. Spontaneous lipid binding to the nicotinic acetylcholine receptor in a native membrane. J Chem Phys 2021; 154:185102. [PMID: 34241006 DOI: 10.1063/5.0046333] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
The nicotinic acetylcholine receptor (nAChR) and other pentameric ligand-gated ion channels are native to neuronal membranes with an unusual lipid composition. While it is well-established that these receptors can be significantly modulated by lipids, the underlying mechanisms have been primarily studied in model membranes with few lipid species. Here, we use coarse-grained molecular dynamics simulation to probe specific binding of lipids in a complex quasi-neuronal membrane. We ran a total of 50 μs of simulations of a single nAChR in a membrane composed of 36 species of lipids. Competition between multiple lipid species produces a complex distribution. We find that overall, cholesterol selects for concave inter-subunit sites and polyunsaturated fatty acids select for convex M4 sites, while monounsaturated and saturated lipids are unenriched in the nAChR boundary. We propose the "density-threshold affinity" as a metric calculated from continuous density distributions, which reduces to a standard affinity in two-state binding. We find that the density-threshold affinity for M4 weakens with chain rigidity, which suggests that flexible chains may help relax packing defects caused by the conical protein shape. For any site, PE headgroups have the strongest affinity of all phospholipid headgroups, but anionic lipids still yield moderately high affinities for the M4 sites as expected. We observe cooperative effects between anionic headgroups and saturated chains at the M4 site in the inner leaflet. We also analyze affinities for individual anionic headgroups. When combined, these insights may reconcile several apparently contradictory experiments on the role of anionic phospholipids in modulating nAChR.
Collapse
Affiliation(s)
- Liam Sharp
- Center for Computational and Integrative Biology, Rutgers University-Camden, Camden, New Jersey 08102, USA
| | - Grace Brannigan
- Center for Computational and Integrative Biology, Rutgers University-Camden, Camden, New Jersey 08102, USA
| |
Collapse
|
28
|
Joseph TT, Bu W, Lin W, Zoubak L, Yeliseev A, Liu R, Eckenhoff RG, Brannigan G. Ketamine Metabolite (2 R,6 R)-Hydroxynorketamine Interacts with μ and κ Opioid Receptors. ACS Chem Neurosci 2021; 12:1487-1497. [PMID: 33905229 PMCID: PMC8154314 DOI: 10.1021/acschemneuro.0c00741] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
![]()
Ketamine is an anesthetic,
analgesic, and antidepressant whose
secondary metabolite (2R,6R)-hydroxynorketamine
(HNK) has N-methyl-d-aspartate-receptor-independent
antidepressant activity in a rodent model. In humans, naltrexone attenuates
its antidepressant effect, consistent with opioid pathway involvement.
No detailed biophysical description is available of opioid receptor
binding of ketamine or its metabolites. Using molecular dynamics simulations
with free energy perturbation, we characterize the binding site and
affinities of ketamine and metabolites in μ and κ opioid
receptors, finding a profound effect of the protonation state. G-protein
recruitment assays show that HNK is an inverse agonist, attenuated
by naltrexone, in these receptors with IC50 values congruous
with our simulations. Overall, our findings are consistent with opioid
pathway involvement in ketamine function.
Collapse
Affiliation(s)
- Thomas T. Joseph
- Department of Anesthesiology and Critical Care, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Weiming Bu
- Department of Anesthesiology and Critical Care, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Wenzhen Lin
- Department of Anesthesiology and Critical Care, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Department of Biochemistry and Molecular Biology, Guangxi Medical University, Nanning, Guangxi 530021, China
| | - Lioudmila Zoubak
- National Institute of Alcohol Abuse and Alcoholism, National Institutes of Health, Bethesda, Maryland 20852, United States
| | - Alexei Yeliseev
- National Institute of Alcohol Abuse and Alcoholism, National Institutes of Health, Bethesda, Maryland 20852, United States
| | - Renyu Liu
- Department of Anesthesiology and Critical Care, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Roderic G. Eckenhoff
- Department of Anesthesiology and Critical Care, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Grace Brannigan
- Center for Computational and Integrative Biology and Department of Physics, Rutgers University, Camden, New Jersey 08102, United States
| |
Collapse
|
29
|
Activation of G-protein-coupled receptors is thermodynamically linked to lipid solvation. Biophys J 2021; 120:1777-1787. [PMID: 33640381 DOI: 10.1016/j.bpj.2021.02.029] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 01/11/2021] [Accepted: 02/22/2021] [Indexed: 12/31/2022] Open
Abstract
Preferential lipid solvation of the G-protein-coupled A2A adenosine receptor (A2AR) is evaluated from 35 μs of all-atom molecular dynamics simulation. A coarse-grained transition matrix algorithm is developed to overcome slow equilibration of the first solvation shell, obtaining estimates of the free energy of solvation by different lipids for the receptor in different activation states. Results indicate preference for solvation by unsaturated chains, which favors the active receptor. A model for lipid-dependent G-protein-coupled receptor activity is proposed in which the chemical potential of lipids in the bulk membrane modulates receptor activity. The entropies associated with moving saturated and unsaturated lipids from bulk to A2AR's first solvation shell are evaluated. Overall, the acyl chains are more disordered (i.e., obtain a favorable entropic contribution) when partitioning to the receptor surface, and this effect is augmented for the saturated chains, which are relatively more ordered in bulk.
Collapse
|
30
|
Cholesterol as a modulator of cannabinoid receptor CB 2 signaling. Sci Rep 2021; 11:3706. [PMID: 33580091 PMCID: PMC7881127 DOI: 10.1038/s41598-021-83245-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 01/20/2021] [Indexed: 11/21/2022] Open
Abstract
Signaling through integral membrane G protein-coupled receptors (GPCRs) is influenced by lipid composition of cell membranes. By using novel high affinity ligands of human cannabinoid receptor CB2, we demonstrate that cholesterol increases basal activation levels of the receptor and alters the pharmacological categorization of these ligands. Our results revealed that (2-(6-chloro-2-((2,2,3,3-tetramethylcyclopropane-1-carbonyl)imino)benzo[d]thiazol-3(2H)-yl)ethyl acetate ligand (MRI-2646) acts as a partial agonist of CB2 in membranes devoid of cholesterol and as a neutral antagonist or a partial inverse agonist in cholesterol-containing membranes. The differential effects of a specific ligand on activation of CB2 in different types of membranes may have implications for screening of drug candidates in a search of modulators of GPCR activity. MD simulation suggests that cholesterol exerts an allosteric effect on the intracellular regions of the receptor that interact with the G-protein complex thereby altering the recruitment of G protein.
Collapse
|
31
|
Abstract
Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.
Collapse
|
32
|
Wague A, Joseph TT, Woll KA, Bu W, Vaidya KA, Bhanu NV, Garcia BA, Nimigean CM, Eckenhoff RG, Riegelhaupt PM. Mechanistic insights into volatile anesthetic modulation of K2P channels. eLife 2020; 9:59839. [PMID: 33345771 PMCID: PMC7781597 DOI: 10.7554/elife.59839] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Accepted: 12/19/2020] [Indexed: 01/01/2023] Open
Abstract
K2P potassium channels are known to be modulated by volatile anesthetic (VA) drugs and play important roles in clinically relevant effects that accompany general anesthesia. Here, we utilize a photoaffinity analog of the VA isoflurane to identify a VA-binding site in the TREK1 K2P channel. The functional importance of the identified site was validated by mutagenesis and biochemical modification. Molecular dynamics simulations of TREK1 in the presence of VA found multiple neighboring residues on TREK1 TM2, TM3, and TM4 that contribute to anesthetic binding. The identified VA-binding region contains residues that play roles in the mechanisms by which heat, mechanical stretch, and pharmacological modulators alter TREK1 channel activity and overlaps with positions found to modulate TASK K2P channel VA sensitivity. Our findings define molecular contacts that mediate VA binding to TREK1 channels and suggest a mechanistic basis to explain how K2P channels are modulated by VAs.
Collapse
Affiliation(s)
- Aboubacar Wague
- Department of Anesthesiology, Weill Cornell Medical College, New York City, United States
| | - Thomas T Joseph
- Department of Anesthesiology and Critical Care, University of Pennsylvania, Philadelphia, United States
| | - Kellie A Woll
- Department of Anesthesiology and Critical Care, University of Pennsylvania, Philadelphia, United States
| | - Weiming Bu
- Department of Anesthesiology and Critical Care, University of Pennsylvania, Philadelphia, United States
| | - Kiran A Vaidya
- Department of Anesthesiology, Weill Cornell Medical College, New York City, United States
| | - Natarajan V Bhanu
- Epigenetics Program, Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States
| | - Benjamin A Garcia
- Epigenetics Program, Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States
| | - Crina M Nimigean
- Department of Anesthesiology, Weill Cornell Medical College, New York City, United States.,Department of Physiology and Biophysics, Weill Cornell Medical College, New York City, United States.,Department of Biochemistry, Weill Cornell Medical College, New York City, United States
| | - Roderic G Eckenhoff
- Department of Anesthesiology and Critical Care, University of Pennsylvania, Philadelphia, United States
| | - Paul M Riegelhaupt
- Department of Anesthesiology, Weill Cornell Medical College, New York City, United States
| |
Collapse
|
33
|
Phillips JC, Hardy DJ, Maia JDC, Stone JE, Ribeiro JV, Bernardi RC, Buch R, Fiorin G, Hénin J, Jiang W, McGreevy R, Melo MCR, Radak BK, Skeel RD, Singharoy A, Wang Y, Roux B, Aksimentiev A, Luthey-Schulten Z, Kalé LV, Schulten K, Chipot C, Tajkhorshid E. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J Chem Phys 2020; 153:044130. [PMID: 32752662 PMCID: PMC7395834 DOI: 10.1063/5.0014475] [Citation(s) in RCA: 1368] [Impact Index Per Article: 342.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Accepted: 07/01/2020] [Indexed: 02/06/2023] Open
Abstract
NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.
Collapse
Affiliation(s)
| | - David J. Hardy
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Julio D. C. Maia
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - John E. Stone
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - João V. Ribeiro
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Rafael C. Bernardi
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | | | - Giacomo Fiorin
- National Heart, Lung and Blood Institute, National
Institutes of Health, Bethesda, Maryland 20814,
USA
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique UPR 9080, CNRS
and Université de Paris, Paris, France
| | | | - Ryan McGreevy
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | | | - Brian K. Radak
- NIH Center for Macromolecular Modeling and
Bioinformatics, Theoretical and Computational Biophysics Group, Beckman Institute for
Advanced Science and Technology, University of Illinois at
Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Robert D. Skeel
- School of Mathematical and Statistical Sciences,
Arizona State University, Tempe, Arizona 85281,
USA
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State
University, Tempe, Arizona 85281, USA
| | - Yi Wang
- Department of Physics, The Chinese University of
Hong Kong, Shatin, Hong Kong, China
| | - Benoît Roux
- Department of Biochemistry, University of
Chicago, Chicago, Illinois 60637, USA
| | | | | | | | | | - Christophe Chipot
- Authors to whom correspondence should be addressed:
and . URL: http://www.ks.uiuc.edu
| | - Emad Tajkhorshid
- Authors to whom correspondence should be addressed:
and . URL: http://www.ks.uiuc.edu
| |
Collapse
|
34
|
Distinct binding of cetirizine enantiomers to human serum albumin and the human histamine receptor H 1. J Comput Aided Mol Des 2020; 34:1045-1062. [PMID: 32572668 DOI: 10.1007/s10822-020-00328-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2020] [Accepted: 06/18/2020] [Indexed: 02/02/2023]
Abstract
Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H1 antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation, non-sedating and highly selective antihistamine. Both enantiomers bind tightly to the human histamine H1 receptor (hH1R) and behave as inverse agonists but the affinity and residence time of (R)-cetirizine are greater than those of (S)-cetirizine. In blood plasma, cetirizine exists in the zwitterionic form and more than 90% of the circulating drug is bound to human serum albumin (HSA), which acts as an inactive reservoir. Independent X-ray crystallographic work has solved the structure of the hH1R:doxepin complex and has identified two drug-binding sites for cetirizine on equine serum albumin (ESA). Given this background, we decided to model a membrane-embedded hH1R in complex with either (R)- or (S)-cetirizine and also the complexes of both ESA and HSA with these two enantiomeric drugs to analyze possible differences in binding modes between enantiomers and also among targets. The ensuing molecular dynamics simulations in explicit solvent and additional computational chemistry calculations provided structural and energetic information about all of these complexes that is normally beyond current experimental possibilities. Overall, we found very good agreement between our binding energy estimates and extant biochemical and pharmacological evidence. A much higher degree of solvent exposure in the cetirizine-binding site(s) of HSA and ESA relative to the more occluded orthosteric binding site in hH1R is translated into larger positional fluctuations and considerably lower affinities for these two nonspecific targets. Whereas it is demonstrated that the two known pockets in ESA provide enough stability for cetirizine binding, only one such site does so in HSA due to a number of amino acid replacements. At the histamine-binding site in hH1R, the distinct interactions established between the phenyl and chlorophenyl moieties of the two enantiomers with the amino acids lining up the pocket and between their free carboxylates and Lys179 in the second extracellular loop account for the improved pharmacological profile of (R)-cetirizine.
Collapse
|
35
|
Corey RA, Stansfeld PJ, Sansom MS. The energetics of protein-lipid interactions as viewed by molecular simulations. Biochem Soc Trans 2020; 48:25-37. [PMID: 31872229 PMCID: PMC7054751 DOI: 10.1042/bst20190149] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Revised: 12/06/2019] [Accepted: 12/10/2019] [Indexed: 12/20/2022]
Abstract
Membranes are formed from a bilayer containing diverse lipid species with which membrane proteins interact. Integral, membrane proteins are embedded in this bilayer, where they interact with lipids from their surroundings, whilst peripheral membrane proteins bind to lipids at the surface of membranes. Lipid interactions can influence the function of membrane proteins, either directly or allosterically. Both experimental (structural) and computational approaches can reveal lipid binding sites on membrane proteins. It is, therefore, important to understand the free energies of these interactions. This affords a more complete view of the engagement of a particular protein with the biological membrane surrounding it. Here, we describe many computational approaches currently in use for this purpose, including recent advances using both free energy and unbiased simulation methods. In particular, we focus on interactions of integral membrane proteins with cholesterol, and with anionic lipids such as phosphatidylinositol 4,5-bis-phosphate and cardiolipin. Peripheral membrane proteins are exemplified via interactions of PH domains with phosphoinositide-containing membranes. We summarise the current state of the field and provide an outlook on likely future directions of investigation.
Collapse
Affiliation(s)
- Robin A. Corey
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
| | - Phillip J. Stansfeld
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
- School of Life Sciences and Department of Chemistry, University of Warwick, Coventry CV4 7AL, U.K
| | - Mark S.P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
| |
Collapse
|
36
|
Corey RA, Vickery ON, Sansom MSP, Stansfeld PJ. Insights into Membrane Protein-Lipid Interactions from Free Energy Calculations. J Chem Theory Comput 2019; 15:5727-5736. [PMID: 31476127 PMCID: PMC6785801 DOI: 10.1021/acs.jctc.9b00548] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
![]()
Integral membrane proteins are regulated
by specific interactions
with lipids from the surrounding bilayer. The structures of protein–lipid
complexes can be determined through a combination of experimental
and computational approaches, but the energetic basis of these interactions
is difficult to resolve. Molecular dynamics simulations provide the
primary computational technique to estimate the free energies of these
interactions. We demonstrate that the energetics of protein–lipid
interactions may be reliably and reproducibly calculated using three
simulation-based approaches: potential of mean force calculations,
alchemical free energy perturbation, and well-tempered metadynamics.
We employ these techniques within the framework of a coarse-grained
force field and apply them to both bacterial and mammalian membrane
protein–lipid systems. We demonstrate good agreement between
the different techniques, providing a robust framework for their automated
implementation within a pipeline for annotation of newly determined
membrane protein structures.
Collapse
Affiliation(s)
- Robin A Corey
- Department of Biochemistry , University of Oxford , South Parks Road , Oxford OX1 3QU , U.K
| | - Owen N Vickery
- Department of Biochemistry , University of Oxford , South Parks Road , Oxford OX1 3QU , U.K
| | - Mark S P Sansom
- Department of Biochemistry , University of Oxford , South Parks Road , Oxford OX1 3QU , U.K
| | - Phillip J Stansfeld
- Department of Biochemistry , University of Oxford , South Parks Road , Oxford OX1 3QU , U.K
| |
Collapse
|
37
|
Marrink SJ, Corradi V, Souza PC, Ingólfsson HI, Tieleman DP, Sansom MS. Computational Modeling of Realistic Cell Membranes. Chem Rev 2019; 119:6184-6226. [PMID: 30623647 PMCID: PMC6509646 DOI: 10.1021/acs.chemrev.8b00460] [Citation(s) in RCA: 422] [Impact Index Per Article: 84.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Indexed: 12/15/2022]
Abstract
Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.
Collapse
Affiliation(s)
- Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Paulo C.T. Souza
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Mark S.P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
38
|
Muller MP, Jiang T, Sun C, Lihan M, Pant S, Mahinthichaichan P, Trifan A, Tajkhorshid E. Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation. Chem Rev 2019; 119:6086-6161. [PMID: 30978005 PMCID: PMC6506392 DOI: 10.1021/acs.chemrev.8b00608] [Citation(s) in RCA: 153] [Impact Index Per Article: 30.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.
Collapse
Affiliation(s)
- Melanie P. Muller
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Tao Jiang
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Chang Sun
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Muyun Lihan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Shashank Pant
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Paween Mahinthichaichan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Anda Trifan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Emad Tajkhorshid
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|