1
|
Kunzmann P, Krumbach JH, Saponaro A, Moroni A, Thiel G, Hamacher K. Anisotropic Network Analysis of Open/Closed HCN4 Channel Advocates Asymmetric Subunit Cooperativity in cAMP Modulation of Gating. J Chem Inf Model 2024; 64:4727-4738. [PMID: 38830626 PMCID: PMC11203669 DOI: 10.1021/acs.jcim.4c00360] [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: 03/01/2024] [Revised: 05/07/2024] [Accepted: 05/10/2024] [Indexed: 06/05/2024]
Abstract
Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are opened in an allosteric manner by membrane hyperpolarization and cyclic nucleotides such as cAMP. Because of conflicting reports from experimental studies on whether cAMP binding to the four available binding sites in the channel tetramer operates cooperatively in gating, we employ here a computational approach as a promising route to examine ligand-induced conformational changes after binding to individual sites. By combining an elastic network model (ENM) with linear response theory (LRT) for modeling the apo-holo transition of the cyclic nucleotide-binding domain (CNBD) in HCN channels, we observe a distinct pattern of cooperativity matching the "positive-negative-positive" cooperativity reported from functional studies. This cooperativity pattern is highly conserved among HCN subtypes (HCN4, HCN1), but only to a lesser extent visible in structurally related channels, which are only gated by voltage (KAT1) or cyclic nucleotides (TAX4). This suggests an inherent cooperativity between subunits in HCN channels as part of a ligand-triggered gating mechanism in these channels.
Collapse
Affiliation(s)
- Patrick Kunzmann
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Jan H. Krumbach
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Andrea Saponaro
- Department
of Pharmacology and Biomolecular Sciences, University of Milano, via Balzaretti 9, 20133 Milano, Italy
| | - Anna Moroni
- Department
of Biosciences, Ion Channel Biophysics, University of Milan, via Celoria 26, 20133 Milan, Italy
| | - Gerhard Thiel
- Department
of Biology, Membrane Biophysics, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
- Centre
for Synthetic Biology, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Kay Hamacher
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
- Centre
for Synthetic Biology, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| |
Collapse
|
2
|
Fontana ACK, Poli AN, Gour J, Srikanth YV, Anastasi N, Ashok D, Khatiwada A, Reeb KL, Cheng MH, Bahar I, Rawls SM, Salvino JM. Synthesis and Structure-Activity Relationships for Glutamate Transporter Allosteric Modulators. J Med Chem 2024; 67:6119-6143. [PMID: 38626917 PMCID: PMC11056993 DOI: 10.1021/acs.jmedchem.3c01909] [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: 10/12/2023] [Revised: 03/23/2024] [Accepted: 03/27/2024] [Indexed: 04/19/2024]
Abstract
Excitatory amino acid transporters (EAATs) are essential CNS proteins that regulate glutamate levels. Excess glutamate release and alteration in EAAT expression are associated with several CNS disorders. Previously, we identified positive allosteric modulators (PAM) of EAAT2, the main CNS transporter, and have demonstrated their neuroprotective properties in vitro. Herein, we report on the structure-activity relationships (SAR) for the analogs identified from virtual screening and from our medicinal chemistry campaign. This work identified several selective EAAT2 positive allosteric modulators (PAMs) such as compounds 4 (DA-023) and 40 (NA-014) from a library of analogs inspired by GT949, an early generation compound. This series also provides nonselective EAAT PAMs, EAAT inhibitors, and inactive compounds that may be useful for elucidating the mechanism of EAAT allosteric modulation.
Collapse
Affiliation(s)
- Andréia C. K. Fontana
- Department
of Pharmacology and Physiology, Drexel University
College of Medicine, Philadelphia, Pennsylvania 19102, United States
| | - Adi N.R. Poli
- Medicinal
Chemistry, Molecular and Cellular Oncogenesis (MCO) Program, The Wistar Institute, Philadelphia, Pennsylvania 19104, United States
| | - Jitendra Gour
- Medicinal
Chemistry, Molecular and Cellular Oncogenesis (MCO) Program, The Wistar Institute, Philadelphia, Pennsylvania 19104, United States
| | - Yellamelli V.V. Srikanth
- Medicinal
Chemistry, Molecular and Cellular Oncogenesis (MCO) Program, The Wistar Institute, Philadelphia, Pennsylvania 19104, United States
| | - Nicholas Anastasi
- Medicinal
Chemistry, Molecular and Cellular Oncogenesis (MCO) Program, The Wistar Institute, Philadelphia, Pennsylvania 19104, United States
| | - Devipriya Ashok
- Medicinal
Chemistry, Molecular and Cellular Oncogenesis (MCO) Program, The Wistar Institute, Philadelphia, Pennsylvania 19104, United States
| | - Apeksha Khatiwada
- Department
of Pharmacology and Physiology, Drexel University
College of Medicine, Philadelphia, Pennsylvania 19102, United States
| | - Katelyn L. Reeb
- Department
of Pharmacology and Physiology, Drexel University
College of Medicine, Philadelphia, Pennsylvania 19102, United States
| | - Mary Hongying Cheng
- Laufer
Center for Physical & Quantitative Biology, Stony Brook University, Stony
Brook, New York 11794, United States
| | - Ivet Bahar
- Department
of Biochemistry and Cell Biology, College of Arts & Sciences and
School of Medicine, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer
Center for Physical & Quantitative Biology, Stony Brook University, Stony
Brook, New York 11794, United States
| | - Scott M. Rawls
- Center
for Substance Abuse Research, Lewis Katz School of Medicine, Temple University, Philadelphia, Pennsylvania 19140United States
| | - Joseph M. Salvino
- Medicinal
Chemistry, Molecular and Cellular Oncogenesis (MCO) Program, The Wistar Institute, Philadelphia, Pennsylvania 19104, United States
- The
Wistar
Cancer Center Molecular Screening, The Wistar
Institute, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
3
|
Lam JH, Nakano A, Katritch V. Scalable computation of anisotropic vibrations for large macromolecular assemblies. Nat Commun 2024; 15:3479. [PMID: 38658556 PMCID: PMC11043083 DOI: 10.1038/s41467-024-47685-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 04/02/2024] [Indexed: 04/26/2024] Open
Abstract
The Normal Mode Analysis (NMA) is a standard approach to elucidate the anisotropic vibrations of macromolecules at their folded states, where low-frequency collective motions can reveal rearrangements of domains and changes in the exposed surface of macromolecules. Recent advances in structural biology have enabled the resolution of megascale macromolecules with millions of atoms. However, the calculation of their vibrational modes remains elusive due to the prohibitive cost associated with constructing and diagonalizing the underlying eigenproblem and the current approaches to NMA are not readily adaptable for efficient parallel computing on graphic processing unit (GPU). Here, we present eigenproblem construction and diagonalization approach that implements level-structure bandwidth-reducing algorithms to transform the sparse computation in NMA to a globally-sparse-yet-locally-dense computation, allowing batched tensor products to be most efficiently executed on GPU. We map, optimize, and compare several low-complexity Krylov-subspace eigensolvers, supplemented by techniques such as Chebyshev filtering, sum decomposition, external explicit deflation and shift-and-inverse, to allow fast GPU-resident calculations. The method allows accurate calculation of the first 1000 vibrational modes of some largest structures in PDB ( > 2.4 million atoms) at least 250 times faster than existing methods.
Collapse
Affiliation(s)
- Jordy Homing Lam
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA
- Bridge Institute and Michelson Center for Convergent Biosciences, University of Southern California, Los Angeles, CA, USA
- Center for New Technologies in Drug Discovery and Development, University of Southern California, Los Angeles, CA, USA
| | - Aiichiro Nakano
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA.
- Department of Physics and Astronomy, University of Southern California, Los Angeles, CA, USA.
- Department of Computer Science, University of Southern California, Los Angeles, CA, USA.
| | - Vsevolod Katritch
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA.
- Bridge Institute and Michelson Center for Convergent Biosciences, University of Southern California, Los Angeles, CA, USA.
- Center for New Technologies in Drug Discovery and Development, University of Southern California, Los Angeles, CA, USA.
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA.
| |
Collapse
|
4
|
Zhao Y, Zhang X, Liu L, Hu F, Chang F, Han Z, Li C. Insights into Activation Dynamics and Functional Sites of Inwardly Rectifying Potassium Channel Kir3.2 by an Elastic Network Model Combined with Perturbation Methods. J Phys Chem B 2024; 128:1360-1370. [PMID: 38308647 DOI: 10.1021/acs.jpcb.3c06739] [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: 02/05/2024]
Abstract
The inwardly rectifying potassium channel Kir3.2, a member of the inward rectifier potassium (Kir) channel family, exerts important biological functions through transporting potassium ions outside of the cell, during which a large-scale synergistic movement occurs among its different domains. Currently, it is not fully understood how the binding of the ligand to the Kir3.2 channel leads to the structural changes and which key residues are responsible for the channel gating and allosteric dynamics. Here, we construct the Gaussian network model (GNM) of the Kir3.2 channel with the secondary structure and covalent interaction information considered (sscGNM), which shows a better performance in reproducing the channel's flexibility compared with the traditional GNM. In addition, the sscANM-based perturbation method is used to simulate the channel's conformational transition caused by the activator PIP2's binding. By applying certain forces to the PIP2 binding pocket, the coarse-grained calculations generate the similar conformational changes to the experimental observation, suggesting that the topology structure as well as PIP2 binding are crucial to the allosteric activation of the Kir3.2 channel. We also utilize the sscGNM-based thermodynamic cycle method developed by us to identify the key residues whose mutations significantly alter the channel's binding free energy with PIP2. We identify not only the residues important for the specific binding but also the ones critical for the allosteric transition coupled with PIP2 binding. This study is helpful for understanding the working mechanism of Kir3.2 channels and can provide important information for related drug design.
Collapse
Affiliation(s)
- Yingchun Zhao
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Xinyu Zhang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Lamei Liu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Fangrui Hu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Fubin Chang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
5
|
Punia R, Goel G. Free Energy Surface and Molecular Mechanism of Slow Structural Transitions in Lipid Bilayers. J Chem Theory Comput 2023; 19:8245-8257. [PMID: 37947833 DOI: 10.1021/acs.jctc.3c00856] [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: 11/12/2023]
Abstract
Lipid membrane remodeling, crucial for many cellular processes, is governed by the coupling of membrane structure and shape fluctuations. Given the importance of the ∼ nm length scale, details of the transition intermediates for conformational change are not fully captured by a continuum-mechanical description. Slow dynamics and the lack of knowledge of reaction coordinates (RCs) for biasing methods pose a challenge for all-atom (AA) simulations. Here, we map system dynamics on Langevin dynamics in a normal mode space determined from an elastic network model representation for the lipid-water Hamiltonian. AA molecular dynamics (MD) simulations are used to determine model parameters, and Langevin dynamics predictions for bilayer structural, mechanical, and dynamic properties are validated against MD simulations and experiments. Transferability to describe the dynamics of a larger lipid bilayer and a heterogeneous membrane-protein system is assessed. A set of generic RCs for pore formation in two tensionless bilayers is obtained by coupling Langevin dynamics to the underlying energy landscape for membrane deformations. Structure evolution is carried out by AA MD, wherein the generic RCs are used in a path metadynamics or an umbrella sampling simulation to determine the thermodynamics of pore formation and its molecular determinants, such as the role of distinct bilayer motions, lipid solvation, and lipid packing.
Collapse
Affiliation(s)
- Rajat Punia
- Department of Chemical Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India
| | - Gaurav Goel
- Department of Chemical Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India
| |
Collapse
|
6
|
Krieger JM, Sorzano COS, Carazo JM. Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing. Int J Mol Sci 2023; 24:14245. [PMID: 37762547 PMCID: PMC10532346 DOI: 10.3390/ijms241814245] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 09/10/2023] [Accepted: 09/15/2023] [Indexed: 09/29/2023] Open
Abstract
Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.
Collapse
Affiliation(s)
- James M. Krieger
- Biocomputing Unit, National Centre for Biotechnology (CNB CSIC), Campus Universidad Autónoma de Madrid, Darwin 3, Cantoblanco, 28049 Madrid, Spain
| | | | - Jose Maria Carazo
- Biocomputing Unit, National Centre for Biotechnology (CNB CSIC), Campus Universidad Autónoma de Madrid, Darwin 3, Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
7
|
Vilardaga JP, Clark LJ, White AD, Sutkeviciute I, Lee JY, Bahar I. Molecular Mechanisms of PTH/PTHrP Class B GPCR Signaling and Pharmacological Implications. Endocr Rev 2023; 44:474-491. [PMID: 36503956 PMCID: PMC10461325 DOI: 10.1210/endrev/bnac032] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 11/14/2022] [Accepted: 12/05/2022] [Indexed: 12/14/2022]
Abstract
The classical paradigm of G protein-coupled receptor (GPCR) signaling via G proteins is grounded in a view that downstream responses are relatively transient and confined to the cell surface, but this notion has been revised in recent years following the identification of several receptors that engage in sustained signaling responses from subcellular compartments following internalization of the ligand-receptor complex. This phenomenon was initially discovered for the parathyroid hormone (PTH) type 1 receptor (PTH1R), a vital GPCR for maintaining normal calcium and phosphate levels in the body with the paradoxical ability to build or break down bone in response to PTH binding. The diverse biological processes regulated by this receptor are thought to depend on its capacity to mediate diverse modes of cyclic adenosine monophosphate (cAMP) signaling. These include transient signaling at the plasma membrane and sustained signaling from internalized PTH1R within early endosomes mediated by PTH. Here we discuss recent structural, cell signaling, and in vivo studies that unveil potential pharmacological outputs of the spatial versus temporal dimension of PTH1R signaling via cAMP. Notably, the combination of molecular dynamics simulations and elastic network model-based methods revealed how precise modulation of PTH signaling responses is achieved through structure-encoded allosteric coupling within the receptor and between the peptide hormone binding site and the G protein coupling interface. The implications of recent findings are now being explored for addressing key questions on how location bias in receptor signaling contributes to pharmacological functions, and how to drug a difficult target such as the PTH1R toward discovering nonpeptidic small molecule candidates for the treatment of metabolic bone and mineral diseases.
Collapse
Affiliation(s)
- Jean-Pierre Vilardaga
- Laboratory for GPCR Biology, Department of Pharmacology and Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Lisa J Clark
- Laboratory for GPCR Biology, Department of Pharmacology and Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Alex D White
- Laboratory for GPCR Biology, Department of Pharmacology and Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Ieva Sutkeviciute
- Laboratory for GPCR Biology, Department of Pharmacology and Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Ji Young Lee
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| |
Collapse
|
8
|
A structure and evolutionary-based classification of solute carriers. iScience 2022; 25:105096. [PMID: 36164651 PMCID: PMC9508557 DOI: 10.1016/j.isci.2022.105096] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Revised: 07/22/2022] [Accepted: 09/04/2022] [Indexed: 11/22/2022] Open
Abstract
Solute carriers are an operationally defined diverse family of membrane proteins involved in the transport of nutrients, metabolites, xenobiotics, and drugs. Here, we provide an integrative classification of solute carriers by combining evolutionary information with proteome-wide structure models recently made available through the AlphaFold resource. Analyses of orthologous relations among 455 protein-coding genes currently classified as human solute carriers, over the fully sequenced genomes of 2,100 species, suggest no more than approximately 180 independent evolutionary origins. Structural comparative analyses provided further insight revealing a total of 24 structurally distinct transmembrane folds, increasing by approximately 40% the number of previously described SLC structural folds. In addition, a structural comparative analysis identified a new human solute carrier member and revealed details of noncanonical ones. Our analyses uncover new ancestral relations between solute carrier genes, provide insights into the evolution of remote homologs and a platform to test hypotheses of functional deorphanization.
Collapse
|
9
|
Sun MY, Zhang X, Yu PC, Liu D, Yang Y, Cui WW, Yang XN, Lei YT, Li XH, Wang WH, Cao P, Wang HS, Zhu MX, Li CZ, Wang R, Fan YZ, Yu Y. Vanilloid agonist-mediated activation of TRPV1 channels requires coordinated movement of the S1-S4 bundle rather than a quiescent state. Sci Bull (Beijing) 2022; 67:1062-1076. [PMID: 36546250 DOI: 10.1016/j.scib.2022.02.016] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Revised: 01/12/2022] [Accepted: 02/16/2022] [Indexed: 01/07/2023]
Abstract
Transient receptor potential vanilloid1 (TRPV1) channel plays an important role in a wide range of physiological and pathological processes, and a comprehensive understanding of TRPV1 gating will create opportunities for therapeutic intervention. Recent incredible advances in cryo-electron microscopy (cryo-EM) have yielded high-resolution structures of all TRPV subtypes (TRPV1-6) and all of them share highly conserved six transmembrane (TM) domains (S1-S6). As revealed by the open structures of TRPV1 in the presence of a bound vanilloid agonist (capsaicin or resiniferatoxin), TM helicesS1 to S4 form a bundle that remains quiescent during channel activation, highlighting differences in the gating mechanism of TRPV1 and voltage-gated ion channels. Here, however, we argue that the structural dynamics rather than quiescence of S1-S4 domains is necessary for capsaicin-mediated activation of TRPV1. Using fluorescent unnatural amino acid (flUAA) incorporation and voltage-clamp fluorometry (VCF) analysis, we directly observed allostery of the S1-S4 bundle upon capsaicin binding. Covalent occupation of VCF-identified sites, single-channel recording, cell apoptosis analysis, and exploration of the role of PSFL828, a novel non-vanilloid agonist we identified, have collectively confirmed the essential role of this coordinated S1-S4 motility in capsaicin-mediated activation of TRPV1. This study concludes that, in contrast to cryo-EM structural studies, vanilloid agonists are also required for S1-S4 movement during TRPV1 activation. Redefining the gating process of vanilloid agonists and the discovery of new non-vanilloid agonists will allow the evaluation of new strategies aimed at the development of TRPV1 modulators.
Collapse
Affiliation(s)
- Meng-Yang Sun
- School of Life Sciences and Key Laboratory of Preclinical Study for New Drugs of Gansu Province, School of Basic Medical Sciences, Lanzhou University, Lanzhou 730000, China; Department of Pharmacology and Chemical Biology, Institute of Medical Sciences, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China; Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Xue Zhang
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Peng-Cheng Yu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Di Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Yang Yang
- Department of Pharmacology and Chemical Biology, Institute of Medical Sciences, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China; Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Wen-Wen Cui
- Department of Pharmacology and Chemical Biology, Institute of Medical Sciences, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Xiao-Na Yang
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China; College of Bioscience and Biotechnology, Hunan Agricultural University, Changsha 410128, China
| | - Yun-Tao Lei
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Xing-Hua Li
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Wen-Hui Wang
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
| | - Peng Cao
- Hospital of Integrated Traditional Chinese and Western Medicine, Nanjing University of Chinese Medicine, Nanjing 210023, China
| | - Heng-Shan Wang
- State Key Laboratory for Chemistry and Molecular Engineering of Medicinal Resources, Collaborative Innovation Center for Guangxi Ethnic Medicine, School of Chemistry and Pharmaceutical Sciences, Guangxi Normal University, Guilin 541004, China
| | - Michael X Zhu
- Department of Integrative Biology and Pharmacology, McGovern Medical School, The University of Texas Health Science Center at Houston, Houston, TX 77030, USA
| | - Chang-Zhu Li
- State Key Laboratory of Utilization of Woody Oil Resource, Hunan Academy of Forestry, Changsha 410004, China
| | - Rui Wang
- School of Life Sciences and Key Laboratory of Preclinical Study for New Drugs of Gansu Province, School of Basic Medical Sciences, Lanzhou University, Lanzhou 730000, China.
| | - Ying-Zhe Fan
- Putuo Hospital, Shanghai University of Chinese Traditional Medicine, Shanghai 200062, China.
| | - Ye Yu
- Department of Pharmacology and Chemical Biology, Institute of Medical Sciences, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China; Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China.
| |
Collapse
|
10
|
Zimmermann MT. Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research. Compr Physiol 2022; 12:3141-3166. [PMID: 35578963 DOI: 10.1002/cphy.c190047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.
Collapse
Affiliation(s)
- Michael T Zimmermann
- Bioinformatics Research and Development Laboratory, Genomic Sciences and Precision Medicine Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Clinical and Translational Sciences Institute, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Biochemistry, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| |
Collapse
|
11
|
Krieger JM, Sorzano COS, Carazo JM, Bahar I. Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0. ACTA CRYSTALLOGRAPHICA SECTION D STRUCTURAL BIOLOGY 2022; 78:399-409. [PMID: 35362464 PMCID: PMC8972803 DOI: 10.1107/s2059798322001966] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Accepted: 02/18/2022] [Indexed: 11/24/2022]
Abstract
New computational biophysics pipelines for analysing the global dynamics of structural ensembles and large, dynamic complexes resolved by cryoEM are reviewed. Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.
Collapse
|
12
|
Tuerkova A, Ungvári O, Laczkó-Rigó R, Mernyák E, Szakács G, Özvegy-Laczka C, Zdrazil B. Data-Driven Ensemble Docking to Map Molecular Interactions of Steroid Analogs with Hepatic Organic Anion Transporting Polypeptides. J Chem Inf Model 2021; 61:3109-3127. [PMID: 34105971 PMCID: PMC8243326 DOI: 10.1021/acs.jcim.1c00362] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Hepatic organic anion transporting polypeptides—OATP1B1,
OATP1B3, and OATP2B1—are expressed at the basolateral membrane
of hepatocytes, being responsible for the uptake of a wide range of
natural substrates and structurally unrelated pharmaceuticals. Impaired
function of hepatic OATPs has been linked to clinically relevant drug–drug
interactions leading to altered pharmacokinetics of administered drugs.
Therefore, understanding the commonalities and differences across
the three transporters represents useful knowledge to guide the drug
discovery process at an early stage. Unfortunately, such efforts remain
challenging because of the lack of experimentally resolved protein
structures for any member of the OATP family. In this study, we established
a rigorous computational protocol to generate and validate structural
models for hepatic OATPs. The multistep procedure is based on the
systematic exploration of available protein structures with shared
protein folding using normal-mode analysis, the calculation of multiple
template backbones from elastic network models, the utilization of
multiple template conformations to generate OATP structural models
with various degrees of conformational flexibility, and the prioritization
of models on the basis of enrichment docking. We employed the resulting
OATP models of OATP1B1, OATP1B3, and OATP2B1 to elucidate binding
modes of steroid analogs in the three transporters. Steroid conjugates
have been recognized as endogenous substrates of these transporters.
Thus, investigating this data set delivers insights into mechanisms
of substrate recognition. In silico predictions were complemented
with in vitro studies measuring the bioactivity of a compound set
on OATP expressing cell lines. Important structural determinants conferring
shared and distinct binding patterns of steroid analogs in the three
transporters have been identified. Overall, this comparative study
provides novel insights into hepatic OATP-ligand interactions and
selectivity. Furthermore, the integrative computational workflow for
structure-based modeling can be leveraged for other pharmaceutical
targets of interest.
Collapse
Affiliation(s)
- Alzbeta Tuerkova
- University of Vienna, Department of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, Althanstraße 14, A-1090 Vienna, Austria
| | - Orsolya Ungvári
- Drug Resistance Research Group, Institute of Enzymology, RCNS, Eötvös Loránd Research Network, H-1117, Budapest, Magyar tudósok krt. 2, Hungary
| | - Réka Laczkó-Rigó
- Drug Resistance Research Group, Institute of Enzymology, RCNS, Eötvös Loránd Research Network, H-1117, Budapest, Magyar tudósok krt. 2, Hungary
| | - Erzsébet Mernyák
- Department of Organic Chemistry, University of Szeged, Dóm tér 8, H-6720 Szeged, Hungary
| | - Gergely Szakács
- Drug Resistance Research Group, Institute of Enzymology, RCNS, Eötvös Loránd Research Network, H-1117, Budapest, Magyar tudósok krt. 2, Hungary.,Department of Medicine I, Institute of Cancer Research, Comprehensive Cancer Center, Medical University of Vienna, A-1090 Vienna, Austria
| | - Csilla Özvegy-Laczka
- Drug Resistance Research Group, Institute of Enzymology, RCNS, Eötvös Loránd Research Network, H-1117, Budapest, Magyar tudósok krt. 2, Hungary
| | - Barbara Zdrazil
- University of Vienna, Department of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, Althanstraße 14, A-1090 Vienna, Austria
| |
Collapse
|
13
|
Zhang Y, Zhang S, Xing J, Bahar I. Normal mode analysis of membrane protein dynamics using the vibrational subsystem analysis. J Chem Phys 2021; 154:195102. [PMID: 34240914 DOI: 10.1063/5.0046710] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The vibrational subsystem analysis is a useful approach that allows for evaluating the spectrum of modes of a given system by integrating out the degrees of freedom accessible to the environment. The approach could be utilized for exploring the collective dynamics of a membrane protein (system) coupled to the lipid bilayer (environment). However, the application to membrane proteins is limited due to high computational costs of modeling a sufficiently large membrane environment unbiased by end effects, which drastically increases the size of the investigated system. We derived a recursive formula for calculating the reduced Hessian of a membrane protein embedded in a lipid bilayer by decomposing the membrane into concentric cylindrical domains with the protein located at the center. The approach allows for the design of a time- and memory-efficient algorithm and a mathematical understanding of the convergence of the reduced Hessian with respect to increasing membrane sizes. The application to the archaeal aspartate transporter GltPh illustrates its utility and efficiency in capturing the transporter's elevator-like movement during its transition between outward-facing and inward-facing states.
Collapse
Affiliation(s)
- Yan Zhang
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Bldg., 3420 Forbes Avenue, Pittsburgh, Pennsylvania 15260, USA
| | - She Zhang
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Bldg., 3420 Forbes Avenue, Pittsburgh, Pennsylvania 15260, USA
| | - Jianhua Xing
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Bldg., 3420 Forbes Avenue, Pittsburgh, Pennsylvania 15260, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Bldg., 3420 Forbes Avenue, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
14
|
Zhang S, Krieger JM, Zhang Y, Kaya C, Kaynak B, Mikulska-Ruminska K, Doruker P, Li H, Bahar I. ProDy 2.0: Increased Scale and Scope after 10 Years of Protein Dynamics Modelling with Python. Bioinformatics 2021; 37:3657-3659. [PMID: 33822884 PMCID: PMC8545336 DOI: 10.1093/bioinformatics/btab187] [Citation(s) in RCA: 68] [Impact Index Per Article: 22.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 03/16/2021] [Indexed: 12/31/2022] Open
Abstract
Summary ProDy, an integrated application programming interface developed for modelling and analysing protein dynamics, has significantly evolved in recent years in response to the growing data and needs of the computational biology community. We present major developments that led to ProDy 2.0: (i) improved interfacing with databases and parsing new file formats, (ii) SignDy for signature dynamics of protein families, (iii) CryoDy for collective dynamics of supramolecular systems using cryo-EM density maps and (iv) essential site scanning analysis for identifying sites essential to modulating global dynamics. Availability and implementation ProDy is open-source and freely available under MIT License from https://github.com/prody/ProDy. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- She Zhang
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - James M Krieger
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Yan Zhang
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Cihan Kaya
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Burak Kaynak
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Karolina Mikulska-Ruminska
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Pemra Doruker
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Hongchun Li
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
15
|
Zhang Y, Krieger J, Mikulska-Ruminska K, Kaynak B, Sorzano COS, Carazo JM, Xing J, Bahar I. State-dependent sequential allostery exhibited by chaperonin TRiC/CCT revealed by network analysis of Cryo-EM maps. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2021; 160:104-120. [PMID: 32866476 PMCID: PMC7914283 DOI: 10.1016/j.pbiomolbio.2020.08.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/17/2019] [Revised: 06/25/2020] [Accepted: 08/16/2020] [Indexed: 12/17/2022]
Abstract
The eukaryotic chaperonin TRiC/CCT plays a major role in assisting the folding of many proteins through an ATP-driven allosteric cycle. Recent structures elucidated by cryo-electron microscopy provide a broad view of the conformations visited at various stages of the chaperonin cycle, including a sequential activation of its subunits in response to nucleotide binding. But we lack a thorough mechanistic understanding of the structure-based dynamics and communication properties that underlie the TRiC/CCT machinery. In this study, we present a computational methodology based on elastic network models adapted to cryo-EM density maps to gain a deeper understanding of the structure-encoded allosteric dynamics of this hexadecameric machine. We have analysed several structures of the chaperonin resolved in different states toward mapping its conformational landscape. Our study indicates that the overall architecture intrinsically favours cooperative movements that comply with the structural variabilities observed in experiments. Furthermore, the individual subunits CCT1-CCT8 exhibit state-dependent sequential events at different states of the allosteric cycle. For example, in the ATP-bound state, subunits CCT5 and CCT4 selectively initiate the lid closure motions favoured by the overall architecture; whereas in the apo form of the heteromer, the subunit CCT7 exhibits the highest predisposition to structural change. The changes then propagate through parallel fluxes of allosteric signals to neighbours on both rings. The predicted state-dependent mechanisms of sequential activation provide new insights into TRiC/CCT intra- and inter-ring signal transduction events.
Collapse
Affiliation(s)
- Yan Zhang
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA, 15261, USA
| | - James Krieger
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA, 15261, USA
| | - Karolina Mikulska-Ruminska
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA, 15261, USA
| | - Burak Kaynak
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA, 15261, USA
| | | | - José-María Carazo
- Centro Nacional de Biotecnología (CSIC), Darwin, 3, 28049, Madrid, Spain
| | - Jianhua Xing
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA, 15261, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA, 15261, USA.
| |
Collapse
|
16
|
Allosteric interactions in the parathyroid hormone GPCR-arrestin complex formation. Nat Chem Biol 2020; 16:1096-1104. [PMID: 32632293 PMCID: PMC7502484 DOI: 10.1038/s41589-020-0567-0] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Accepted: 05/13/2020] [Indexed: 11/09/2022]
Abstract
Peptide ligands of class B G-protein-coupled receptors act via a two-step binding process, but the essential mechanisms that link their extracellular binding to intracellular receptor-arrestin interactions are not fully understood. Using NMR, crosslinking coupled to mass spectrometry, signaling experiments and computational approaches on the parathyroid hormone (PTH) type 1 receptor (PTHR), we show that initial binding of the PTH C-terminal part constrains the conformation of the flexible PTH N-terminal signaling epitope before a second binding event occurs. A 'hot-spot' PTH residue, His9, that inserts into the PTHR transmembrane domain at this second step allosterically engages receptor-arrestin coupling. A conformational change in PTHR intracellular loop 3 permits favorable interactions with β-arrestin's finger loop. These results unveil structural determinants for PTHR-arrestin complex formation and reveal that the two-step binding mechanism proceeds via cooperative fluctuations between ligand and receptor, which extend to other class B G-protein-coupled receptors.
Collapse
|
17
|
Li H, Chang YY, Lee JY, Bahar I, Yang LW. DynOmics: dynamics of structural proteome and beyond. Nucleic Acids Res 2019; 45:W374-W380. [PMID: 28472330 PMCID: PMC5793847 DOI: 10.1093/nar/gkx385] [Citation(s) in RCA: 99] [Impact Index Per Article: 19.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2017] [Accepted: 04/25/2017] [Indexed: 01/07/2023] Open
Abstract
DynOmics (dynomics.pitt.edu) is a portal developed to leverage rapidly growing structural proteomics data by efficiently and accurately evaluating the dynamics of structurally resolved systems, from individual molecules to large complexes and assemblies, in the context of their physiological environment. At the core of the portal is a newly developed server, ENM 1.0, which permits users to efficiently generate information on the collective dynamics of any structure in PDB format, user-uploaded or database-retrieved. ENM 1.0 integrates two widely used elastic network models (ENMs)—the Gaussian Network Model (GNM) and the Anisotropic Network Model (ANM), extended to take account of molecular environment. It enables users to assess potentially functional sites, signal transduction or allosteric communication mechanisms, and protein–protein and protein–DNA interaction poses, in addition to delivering ensembles of accessible conformers reconstructed at atomic details based on the global modes of motions predicted by the ANM. The ‘environment’ is defined in a flexible manner, from lipid bilayer and crystal contacts, to substrate or ligands bound to a protein, or surrounding subunits in a multimeric structure or assembly. User-friendly interactive features permit users to easily visualize how the environment alter the intrinsic dynamics of the query systems. ENM 1.0 can be accessed at http://enm.pitt.edu/ or http://dyn.life.nthu.edu.tw/oENM/.
Collapse
Affiliation(s)
- Hongchun Li
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh PA, 15213, USA
| | - Yuan-Yu Chang
- Institute of Bioinformatics and Structural Biology, National Tsing-Hua University, Taiwan
| | - Ji Young Lee
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh PA, 15213, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh PA, 15213, USA
| | - Lee-Wei Yang
- Institute of Bioinformatics and Structural Biology, National Tsing-Hua University, Taiwan
| |
Collapse
|
18
|
Ponzoni L, Zhang S, Cheng MH, Bahar I. Shared dynamics of LeuT superfamily members and allosteric differentiation by structural irregularities and multimerization. Philos Trans R Soc Lond B Biol Sci 2019; 373:rstb.2017.0177. [PMID: 29735731 PMCID: PMC5941172 DOI: 10.1098/rstb.2017.0177] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/19/2017] [Indexed: 12/14/2022] Open
Abstract
The LeuT-fold superfamily includes secondary active transporters from different functional families, which share a common tertiary structure, despite having a remarkably low sequence similarity. By identifying the common structural and dynamical features upon principal component analysis of a comprehensive ensemble of 90 experimentally resolved structures and anisotropic network model evaluation of collective motions, we provide a unified point of view for understanding the reasons why this particular fold has been selected by evolution to accomplish such a broad spectrum of functions. The parallel identification of conserved sequence features, localized at specific sites of transmembrane helices, sheds light on the role of broken helices (TM1 and TM6 in LeuT) in promoting ion/substrate binding and allosteric interconversion between the outward- and inward-facing conformations of transporters. Finally, the determination of the dynamics landscape for the structural ensemble provides a promising framework for the classification of transporters based on their dynamics, and the characterization of the collective movements that favour multimerization.This article is part of a discussion meeting issue 'Allostery and molecular machines'.
Collapse
Affiliation(s)
- Luca Ponzoni
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - She Zhang
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Mary Hongying Cheng
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| |
Collapse
|
19
|
Enkavi G, Javanainen M, Kulig W, Róg T, Vattulainen I. Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance. Chem Rev 2019; 119:5607-5774. [PMID: 30859819 PMCID: PMC6727218 DOI: 10.1021/acs.chemrev.8b00538] [Citation(s) in RCA: 173] [Impact Index Per Article: 34.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
![]()
Biological
membranes are tricky to investigate. They are complex
in terms of molecular composition and structure, functional
over a wide range of time scales, and characterized
by nonequilibrium conditions. Because of all of these
features, simulations are a great technique to study biomembrane
behavior. A significant part of the functional processes
in biological membranes takes place at the molecular
level; thus computer simulations are the method of
choice to explore how their properties emerge from specific
molecular features and how the interplay among the numerous
molecules gives rise to function over spatial and
time scales larger than the molecular ones. In this
review, we focus on this broad theme. We discuss the current
state-of-the-art of biomembrane simulations that, until
now, have largely focused on a rather narrow picture
of the complexity of the membranes. Given this, we
also discuss the challenges that we should unravel in the
foreseeable future. Numerous features such as the actin-cytoskeleton
network, the glycocalyx network, and nonequilibrium
transport under ATP-driven conditions have so far
received very little attention; however, the potential
of simulations to solve them would be exceptionally high. A
major milestone for this research would be that one day
we could say that computer simulations genuinely research
biological membranes, not just lipid bilayers.
Collapse
Affiliation(s)
- Giray Enkavi
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland
| | - Matti Javanainen
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland.,Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo naḿesti 542/2 , 16610 Prague , Czech Republic.,Computational Physics Laboratory , Tampere University , P.O. Box 692, FI-33014 Tampere , Finland
| | - Waldemar Kulig
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland
| | - Tomasz Róg
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland.,Computational Physics Laboratory , Tampere University , P.O. Box 692, FI-33014 Tampere , Finland
| | - Ilpo Vattulainen
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland.,Computational Physics Laboratory , Tampere University , P.O. Box 692, FI-33014 Tampere , Finland.,MEMPHYS-Center for Biomembrane Physics
| |
Collapse
|
20
|
van Dijk L, Giladi M, Refaeli B, Hiller R, Cheng MH, Bahar I, Khananshvili D. Key residues controlling bidirectional ion movements in Na +/Ca 2+ exchanger. Cell Calcium 2018; 76:10-22. [PMID: 30248574 PMCID: PMC6688843 DOI: 10.1016/j.ceca.2018.09.004] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Revised: 09/13/2018] [Accepted: 09/13/2018] [Indexed: 01/06/2023]
Abstract
Prokaryotic and eukaryotic Na+/Ca2+ exchangers (NCX) control Ca2+ homeostasis. NCX orthologs exhibit up to 104-fold differences in their turnover rates (kcat), whereas the ratios between the cytosolic (cyt) and extracellular (ext) Km values (Kint = KmCyt/KmExt) are highly asymmetric and alike (Kint ≤ 0.1) among NCXs. The structural determinants controlling a huge divergence in kcat at comparable Kint remain unclear, although 11 (out of 12) ion-coordinating residues are highly conserved among NCXs. The crystal structure of the archaeal NCX (NCX_Mj) was explored for testing the mutational effects of pore-allied and loop residues on kcat and Kint. Among 55 tested residues, 26 mutations affect either kcat or Kint, where two major groups can be distinguished. The first group of mutations (14 residues) affect kcat rather than Kint. The majority of these residues (10 out of 14) are located within the extracellular vestibule near the pore center. The second group of mutations (12 residues) affect Kint rather than kcat, whereas the majority of residues (9 out 12) are randomly dispersed within the extracellular vestibule. In conjunction with computational modeling-simulations and hydrogen-deuterium exchange mass-spectrometry (HDX-MS), the present mutational analysis highlights structural elements that differentially govern the intrinsic asymmetry and transport rates. The key residues, located at specific segments, can affect the characteristic features of local backbone dynamics and thus, the conformational flexibility of ion-transporting helices contributing to critical conformational transitions. The underlying mechanisms might have a physiological relevance for matching the response modes of NCX variants to cell-specific Ca2+ and Na+ signaling.
Collapse
Affiliation(s)
- Liat van Dijk
- Department of Physiology and Pharmacology, Sackler School of Medicine, Tel-Aviv University, Ramat-Aviv, Tel-Aviv, 69978, Israel
| | - Moshe Giladi
- Department of Physiology and Pharmacology, Sackler School of Medicine, Tel-Aviv University, Ramat-Aviv, Tel-Aviv, 69978, Israel
| | - Bosmat Refaeli
- Department of Physiology and Pharmacology, Sackler School of Medicine, Tel-Aviv University, Ramat-Aviv, Tel-Aviv, 69978, Israel
| | - Reuben Hiller
- Department of Physiology and Pharmacology, Sackler School of Medicine, Tel-Aviv University, Ramat-Aviv, Tel-Aviv, 69978, Israel
| | - Mary Hongying Cheng
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, 15213, USA
| | - Ivet Bahar
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, 15213, USA.
| | - Daniel Khananshvili
- Department of Physiology and Pharmacology, Sackler School of Medicine, Tel-Aviv University, Ramat-Aviv, Tel-Aviv, 69978, Israel.
| |
Collapse
|
21
|
Allosteric Modulation of Intact γ-Secretase Structural Dynamics. Biophys J 2018; 113:2634-2649. [PMID: 29262358 DOI: 10.1016/j.bpj.2017.10.012] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 09/26/2017] [Accepted: 10/10/2017] [Indexed: 12/20/2022] Open
Abstract
As a protease complex involved in the cleavage of amyloid precursor proteins that lead to the formation of amyloid β fibrils implicated in Alzheimer's disease, γ-secretase is an important target for developing therapeutics against Alzheimer's disease. γ-secretase is composed of four subunits: nicastrin (NCT) in the extracellular (EC) domain, presenilin-1 (PS1), anterior pharynx defective 1, and presenilin enhancer 2 in the transmembrane (TM) domain. NCT and PS1 play important roles in binding amyloid β precursor proteins and modulating PS1 catalytic activity. Yet, the molecular mechanisms of coupling between substrate/modulator binding and catalytic activity remain to be elucidated. Recent determination of intact human γ-secretase cryo-electron microscopy structure has opened the way for a detailed investigation of the structural dynamics of this complex. Our analysis, based on a membrane-coupled anisotropic network model, reveals two types of NCT motions, bending and twisting, with respect to PS1. These underlie the fluctuations between the "open" and "closed" states of the lid-like NCT with respect to a hydrophilic loop 1 (HL1) on PS1, thus allowing or blocking access of the substrate peptide (EC portion) to HL1 and to the neighboring helix TM2. In addition to this alternating access mechanism, fluctuations in the volume of the PS1 central cavity facilitate the exposure of the catalytic site for substrate cleavage. Druggability simulations show that γ-secretase presents several hot spots for either orthosteric or allosteric inhibition of catalytic activity, consistent with experimental data. In particular, a hinge region at the interface between the EC and TM domains, near the interlobe groove of NCT, emerges as an allo-targeting site that would impact the coupling between HL1/TM2 and the catalytic pocket, opening, to our knowledge, new avenues for structure-based design of novel allosteric modulators of γ-secretase protease activity.
Collapse
|
22
|
Wang J, Albers T, Grewer C. Energy Landscape of the Substrate Translocation Equilibrium of Plasma-Membrane Glutamate Transporters. J Phys Chem B 2017; 122:28-39. [PMID: 29218993 DOI: 10.1021/acs.jpcb.7b09059] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Glutamate transporters maintain a large glutamate concentration gradient across synaptic membranes and are, thus, critical for functioning of the excitatory synapse. Mammalian glutamate transporters concentrate glutamate inside cells through energetic coupling of glutamate flux to the transmembrane concentration gradient of Na+. Structural models based on an archeal homologue, GltPh, suggest an elevator-like carrier mechanism. However, the energetic determinants of this carrier-based movement are not well understood. Although electrostatics play an important role in governing these energetics, their implication on transport dynamics has not been studied. Here, we combine a pre-steady-state kinetic analysis of the translocation equilibrium with electrostatic computations to gain insight into the energetics of the translocation process. Our results show the biphasic nature of translocation, consistent with the existence of an intermediate on the translocation pathway. In the absence of voltage, the equilibrium is shifted to the outward-facing configuration. Electrostatic computations confirm the intermediate state and show that the elevator-like movement is energetically feasible in the presence of bound Na+ ions, whereas a substrate-hopping model is energetically prohibitive. Our results highlight the critical contribution of charge compensation to transport and add to results from previous molecular dynamics simulations for improved understanding of the glutamate translocation process.
Collapse
Affiliation(s)
- Jiali Wang
- Department of Chemistry, Binghamton University , Binghamton, New York 13902, United States
| | - Thomas Albers
- Department of Chemistry, Binghamton University , Binghamton, New York 13902, United States
| | - Christof Grewer
- Department of Chemistry, Binghamton University , Binghamton, New York 13902, United States
| |
Collapse
|
23
|
Kaynak BT, Findik D, Doruker P. RESPEC Incorporates Residue Specificity and the Ligand Effect into the Elastic Network Model. J Phys Chem B 2017; 122:5347-5355. [DOI: 10.1021/acs.jpcb.7b10325] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Burak T. Kaynak
- Department of Physics, Bogazici University, 34342, Bebek, Istanbul, Turkey
| | - Doga Findik
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, 34342, Bebek, Istanbul, Turkey
| | - Pemra Doruker
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, 34342, Bebek, Istanbul, Turkey
| |
Collapse
|
24
|
Li H, Sharma N, General IJ, Schreiber G, Bahar I. Dynamic Modulation of Binding Affinity as a Mechanism for Regulating Interferon Signaling. J Mol Biol 2017; 429:2571-2589. [PMID: 28648616 PMCID: PMC5545807 DOI: 10.1016/j.jmb.2017.06.011] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2017] [Revised: 06/15/2017] [Accepted: 06/16/2017] [Indexed: 12/22/2022]
Abstract
How structural dynamics affects cytokine signaling is under debate. Here, we investigated the dynamics of the type I interferon (IFN) receptor, IFNAR1, and its effect on signaling upon binding IFN and IFNAR2 using a combination of structure-based mechanistic studies, in situ binding, and gene induction assays. Our study reveals that IFNAR1 flexibility modulates ligand-binding affinity, which, in turn, regulates biological signaling. We identified the hinge sites and key interactions implicated in IFNAR1 inter-subdomain (SD1-SD4) movements. We showed that the predicted cooperative movements are essential to accommodate intermolecular interactions. Engineered disulfide bridges, computationally predicted to interfere with IFNAR1 dynamics, were experimentally confirmed. Notably, introducing disulfide bonds between subdomains SD2 and SD3 modulated IFN binding and activity in accordance with the relative attenuation of cooperative movements with varying distance from the hinge center, whereas locking the SD3-SD4 interface flexibility in favor of an extended conformer increased activity.
Collapse
Affiliation(s)
- Hongchun Li
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Nanaocha Sharma
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Ignacio J General
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA; School of Science and Technology, and CONICET, Universidad Nacional de San Martin, San Martin, Buenos Aires 1650, Argentina
| | - Gideon Schreiber
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot 76100, Israel.
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA.
| |
Collapse
|
25
|
Cheng MH, Torres-Salazar D, Gonzalez-Suarez AD, Amara SG, Bahar I. Substrate transport and anion permeation proceed through distinct pathways in glutamate transporters. eLife 2017; 6. [PMID: 28569666 PMCID: PMC5472439 DOI: 10.7554/elife.25850] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2017] [Accepted: 05/10/2017] [Indexed: 11/13/2022] Open
Abstract
Advances in structure-function analyses and computational biology have enabled a deeper understanding of how excitatory amino acid transporters (EAATs) mediate chloride permeation and substrate transport. However, the mechanism of structural coupling between these functions remains to be established. Using a combination of molecular modeling, substituted cysteine accessibility, electrophysiology and glutamate uptake assays, we identified a chloride-channeling conformer, iChS, transiently accessible as EAAT1 reconfigures from substrate/ion-loaded into a substrate-releasing conformer. Opening of the anion permeation path in this iChS is controlled by the elevator-like movement of the substrate-binding core, along with its wall that simultaneously lines the anion permeation path (global); and repacking of a cluster of hydrophobic residues near the extracellular vestibule (local). Moreover, our results demonstrate that stabilization of iChS by chemical modifications favors anion channeling at the expense of substrate transport, suggesting a mutually exclusive regulation mediated by the movement of the flexible wall lining the two regions.
Collapse
Affiliation(s)
- Mary Hongying Cheng
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, United States
| | - Delany Torres-Salazar
- Laboratory of Molecular and Cellular Neurobiology, National Institute of Mental Health, National Institutes of Health, Bethesda, United States
| | - Aneysis D Gonzalez-Suarez
- Laboratory of Molecular and Cellular Neurobiology, National Institute of Mental Health, National Institutes of Health, Bethesda, United States
| | - Susan G Amara
- Laboratory of Molecular and Cellular Neurobiology, National Institute of Mental Health, National Institutes of Health, Bethesda, United States
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, United States
| |
Collapse
|
26
|
Bergman S, Lezon TR. Modeling global changes induced by local perturbations to the HIV-1 capsid. J Mol Graph Model 2016; 71:218-226. [PMID: 27951510 DOI: 10.1016/j.jmgm.2016.12.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2016] [Revised: 11/04/2016] [Accepted: 12/01/2016] [Indexed: 02/02/2023]
Abstract
The HIV-1 capsid is a conical protein shell made up of hexamers and pentamers of the capsid protein. The capsid houses the viral genome and replication machinery, and its opening, or uncoating, within the host cell marks a critical step in the HIV-1 lifecycle. Binding of host factors such as TRIM5α and cyclophilin A (CypA) can alter the capsid's stability, accelerating or delaying the onset of uncoating and disrupting infectivity. We employ coarse-grained computational modeling to investigate the effects of point mutations and host factor binding on HIV-1 capsid stability. We find that the largest fluctuations occur in the low-curvature regions of the capsid, and that its structural dynamics are affected by perturbations at the inter-hexamer interfaces and near the CypA binding loop, suggesting roles for these features in capsid stability. Our models show that linking capsid proteins across hexamers attenuates vibration in the low-curvature regions of the capsid, but that linking within hexamers does not. These results indicate a possible mechanism through which CypA binding alters capsid stability and highlight the utility of coarse-grained network modeling for understanding capsid mechanics.
Collapse
Affiliation(s)
- Shana Bergman
- Department of Computational and Systems Biology, University of Pittsburgh, Suite 3064 Biomedical Science Tower 3, 3501 Fifth Avenue, Pittsburgh, PA 15260, USA.
| | - Timothy R Lezon
- Department of Computational and Systems Biology, University of Pittsburgh, Suite 3064 Biomedical Science Tower 3, 3501 Fifth Avenue, Pittsburgh, PA 15260, USA; University of Pittsburgh Drug Discovery Institute, University of Pittsburgh, W965 Biomedical Science Tower, 200 Lothrop Street, Pittsburgh, PA 15261, USA.
| |
Collapse
|
27
|
LeVine MV, Cuendet MA, Khelashvili G, Weinstein H. Allosteric Mechanisms of Molecular Machines at the Membrane: Transport by Sodium-Coupled Symporters. Chem Rev 2016; 116:6552-87. [PMID: 26892914 DOI: 10.1021/acs.chemrev.5b00627] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Solute transport across cell membranes is ubiquitous in biology as an essential physiological process. Secondary active transporters couple the unfavorable process of solute transport against its concentration gradient to the energetically favorable transport of one or several ions. The study of such transporters over several decades indicates that their function involves complex allosteric mechanisms that are progressively being revealed in atomistic detail. We focus on two well-characterized sodium-coupled symporters: the bacterial amino acid transporter LeuT, which is the prototype for the "gated pore" mechanism in the mammalian synaptic monoamine transporters, and the archaeal GltPh, which is the prototype for the "elevator" mechanism in the mammalian excitatory amino acid transporters. We present the evidence for the role of allostery in the context of a quantitative formalism that can reconcile biochemical and biophysical data and thereby connects directly to recent insights into the molecular structure and dynamics of these proteins. We demonstrate that, while the structures and mechanisms of these transporters are very different, the available data suggest a common role of specific models of allostery in their functions. We argue that such allosteric mechanisms appear essential not only for sodium-coupled symport in general but also for the function of other types of molecular machines in the membrane.
Collapse
Affiliation(s)
- Michael V LeVine
- Department of Physiology and Biophysics, ‡HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Cornell Medical College of Cornell University , New York, New York 10065, United States
| | - Michel A Cuendet
- Department of Physiology and Biophysics, ‡HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Cornell Medical College of Cornell University , New York, New York 10065, United States
| | - George Khelashvili
- Department of Physiology and Biophysics, ‡HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Cornell Medical College of Cornell University , New York, New York 10065, United States
| | - Harel Weinstein
- Department of Physiology and Biophysics, ‡HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Cornell Medical College of Cornell University , New York, New York 10065, United States
| |
Collapse
|
28
|
Computational approaches to detect allosteric pathways in transmembrane molecular machines. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:1652-62. [PMID: 26806157 DOI: 10.1016/j.bbamem.2016.01.010] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2015] [Revised: 01/13/2016] [Accepted: 01/14/2016] [Indexed: 01/05/2023]
Abstract
Many of the functions of transmembrane proteins involved in signal processing and transduction across the cell membrane are determined by allosteric couplings that propagate the functional effects well beyond the original site of activation. Data gathered from breakthroughs in biochemistry, crystallography, and single molecule fluorescence have established a rich basis of information for the study of molecular mechanisms in the allosteric couplings of such transmembrane proteins. The mechanistic details of these couplings, many of which have therapeutic implications, however, have only become accessible in synergy with molecular modeling and simulations. Here, we review some recent computational approaches that analyze allosteric coupling networks (ACNs) in transmembrane proteins, and in particular the recently developed Protein Interaction Analyzer (PIA) designed to study ACNs in the structural ensembles sampled by molecular dynamics simulations. The power of these computational approaches in interrogating the functional mechanisms of transmembrane proteins is illustrated with selected examples of recent experimental and computational studies pursued synergistically in the investigation of secondary active transporters and GPCRs. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
Collapse
|
29
|
Haliloglu T, Bahar I. Adaptability of protein structures to enable functional interactions and evolutionary implications. Curr Opin Struct Biol 2015; 35:17-23. [PMID: 26254902 DOI: 10.1016/j.sbi.2015.07.007] [Citation(s) in RCA: 75] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2015] [Revised: 07/15/2015] [Accepted: 07/20/2015] [Indexed: 12/21/2022]
Abstract
Several studies in recent years have drawn attention to the ability of proteins to adapt to intermolecular interactions by conformational changes along structure-encoded collective modes of motions. These so-called soft modes, primarily driven by entropic effects, facilitate, if not enable, functional interactions. They represent excursions on the conformational space along principal low-ascent directions/paths away from the original free energy minimum, and they are accessible to the protein even before protein-protein/ligand interactions. An emerging concept from these studies is the evolution of structures or modular domains to favor such modes of motion that will be recruited or integrated for enabling functional interactions. Structural dynamics, including the allosteric switches in conformation that are often stabilized upon formation of complexes and multimeric assemblies, emerge as key properties that are evolutionarily maintained to accomplish biological activities, consistent with the paradigm sequence→structure→dynamics→function where 'dynamics' bridges structure and function.
Collapse
Affiliation(s)
- Turkan Haliloglu
- Department of Chemical Engineering and Polymer Research Center, and Center for Life Sciences and Technologies, Bogazici University, 34342 Istanbul, Turkey; Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA.
| |
Collapse
|
30
|
Cooperative Dynamics of Intact AMPA and NMDA Glutamate Receptors: Similarities and Subfamily-Specific Differences. Structure 2015; 23:1692-1704. [PMID: 26256538 DOI: 10.1016/j.str.2015.07.002] [Citation(s) in RCA: 65] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2015] [Revised: 06/19/2015] [Accepted: 07/01/2015] [Indexed: 01/03/2023]
Abstract
Ionotropic glutamate receptors (iGluRs) are tetrameric ion channels that mediate excitatory neurotransmission. Recent structures of α-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) and N-methyl-D-aspartate (NMDA) receptors permit a comparative analysis of whole-receptor dynamics for the first time. Despite substantial differences in the packing of their two-domain extracellular region, the two iGluRs share similar dynamics, elucidated by elastic network models. Motions accessible to either structure enable conformational interconversion, such as compression of the AMPA receptor toward the more tightly packed NMDA receptor conformation, which has been linked to allosteric regulation. Pivoting motions coupled to concerted rotations of the transmembrane ion channel are prominent between dimers of distal N-terminal domains in the loosely packed AMPA receptor. The occurrence and functional relevance of these motions is verified by cross-linking experiments designed to probe the computationally predicted distance changes. Together with the identification of hotspot residues acting as mediators of allosteric communication, our data provide a glimpse into the dynamic spectrum of iGluRs.
Collapse
|
31
|
Bahar I, Cheng MH, Lee JY, Kaya C, Zhang S. Structure-Encoded Global Motions and Their Role in Mediating Protein-Substrate Interactions. Biophys J 2015; 109:1101-9. [PMID: 26143655 DOI: 10.1016/j.bpj.2015.06.004] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2015] [Revised: 06/02/2015] [Accepted: 06/03/2015] [Indexed: 12/22/2022] Open
Abstract
Recent structure-based computational studies suggest that, in contrast to the classical description of equilibrium fluctuations as wigglings and jigglings, proteins have access to well-defined spectra of collective motions, called intrinsic dynamics, encoded by their structure under native state conditions. In particular, the global modes of motions (at the low frequency end of the spectrum) are shown by multiple studies to be highly robust to minor differences in the structure or to detailed interactions at the atomic level. These modes, encoded by the overall fold, usually define the mechanisms of interactions with substrates. They can be estimated by low-resolution models such as the elastic network models (ENMs) exclusively based on interresidue contact topology. The ability of ENMs to efficiently assess the global motions intrinsically favored by the overall fold as well as the relevance of these predictions to the dominant changes in structure experimentally observed for a given protein in the presence of different substrates suggest that the intrinsic dynamics plays a role in mediating protein-substrate interactions. These observations underscore the functional significance of structure-encoded dynamics, or the importance of the predisposition to favor functional global modes in the evolutionary selection of structures.
Collapse
Affiliation(s)
- Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania.
| | - Mary Hongying Cheng
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania
| | - Ji Young Lee
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania
| | - Cihan Kaya
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania
| | - She Zhang
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania
| |
Collapse
|
32
|
Functional mechanisms of neurotransmitter transporters regulated by lipid-protein interactions of their terminal loops. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2015; 1848:1765-74. [PMID: 25847498 DOI: 10.1016/j.bbamem.2015.03.025] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2015] [Revised: 03/16/2015] [Accepted: 03/23/2015] [Indexed: 01/28/2023]
Abstract
The physiological functions of neurotransmitter:sodium symporters (NSS) in reuptake of neurotransmitters from the synapse into the presynaptic nerve have been shown to be complemented by their involvement, together with non-plasma membrane neurotransmitter transporters, in the reverse transport of substrate (efflux) in response to psychostimulants. Recent experimental evidence implicates highly anionic phosphatidylinositol 4,5-biphosphate (PIP(2)) lipids in such functions of the serotonin (SERT) and dopamine (DAT) transporters. Thus, for both SERT and DAT, neurotransmitter efflux has been shown to be strongly regulated by the presence of PIP(2) lipids in the plasma membrane, and the electrostatic interaction of the N-terminal region of DAT with the negatively charged PIP(2) lipids. We examine the experimentally established phenotypes in a structural context obtained from computational modeling based on recent crystallographic data. The results are shown to set the stage for a mechanistic understanding of physiological actions of neurotransmitter transporters in the NSS family of membrane proteins. This article is part of a Special Issue entitled: Lipid-protein interactions.
Collapse
|
33
|
McIlwain BC, Vandenberg RJ, Ryan RM. Transport rates of a glutamate transporter homologue are influenced by the lipid bilayer. J Biol Chem 2015; 290:9780-8. [PMID: 25713135 DOI: 10.1074/jbc.m114.630590] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2014] [Indexed: 11/06/2022] Open
Abstract
The aspartate transporter from Pyrococcus horikoshii (GltPh) is a model for the structure of the SLC1 family of amino acid transporters. Crystal structures of GltPh provide insight into mechanisms of ion coupling and substrate transport; however, structures have been solved in the absence of a lipid bilayer so they provide limited information regarding interactions that occur between the protein and lipids of the membrane. Here, we investigated the effect of the lipid environment on aspartate transport by reconstituting GltPh into liposomes of defined lipid composition where the primary lipid is phosphatidylethanolamine (PE) or its methyl derivatives. We showed that the rate of aspartate transport and the transmembrane orientation of GltPh were influenced by the primary lipid in the liposomes. In PE liposomes, we observed the highest transport rate and showed that 85% of the transporters were orientated right-side out, whereas in trimethyl PE liposomes, 50% of transporters were right-side out, and we observed a 4-fold reduction in transport rate. Differences in orientation can only partially explain the lipid composition effect on transport rate. Crystal structures of GltPh revealed a tyrosine residue (Tyr-33) that we propose interacts with lipid headgroups during the transport cycle. Based on site-directed mutagenesis, we propose that a cation-π interaction between Tyr-33 and the lipid headgroups can influence conformational flexibility of the trimerization domain and thus the rate of transport. These results provide a specific example of how interactions between membrane lipids and membrane-bound proteins can influence function and highlight the importance of the role of the membrane in transporter function.
Collapse
Affiliation(s)
- Benjamin C McIlwain
- From the Transporter Biology Group, Discipline of Pharmacology and Bosch Institute, University of Sydney, Sydney, New South Wales 2006, Australia
| | - Robert J Vandenberg
- From the Transporter Biology Group, Discipline of Pharmacology and Bosch Institute, University of Sydney, Sydney, New South Wales 2006, Australia
| | - Renae M Ryan
- From the Transporter Biology Group, Discipline of Pharmacology and Bosch Institute, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
34
|
Coupling between neurotransmitter translocation and protonation state of a titratable residue during Na ⁺-coupled transport. Biophys J 2015; 106:2547-8. [PMID: 24940770 DOI: 10.1016/j.bpj.2014.05.011] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2014] [Revised: 05/05/2014] [Accepted: 05/05/2014] [Indexed: 12/31/2022] Open
|
35
|
Mondal S, Khelashvili G, Weinstein H. Not just an oil slick: how the energetics of protein-membrane interactions impacts the function and organization of transmembrane proteins. Biophys J 2014; 106:2305-16. [PMID: 24896109 PMCID: PMC4052241 DOI: 10.1016/j.bpj.2014.04.032] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2014] [Revised: 03/31/2014] [Accepted: 04/23/2014] [Indexed: 01/12/2023] Open
Abstract
The membrane environment, its composition, dynamics, and remodeling, have been shown to participate in the function and organization of a wide variety of transmembrane (TM) proteins, making it necessary to study the molecular mechanisms of such proteins in the context of their membrane settings. We review some recent conceptual advances enabling such studies, and corresponding computational models and tools designed to facilitate the concerted experimental and computational investigation of protein-membrane interactions. To connect productively with the high resolution achieved by cognate experimental approaches, the computational methods must offer quantitative data at an atomistically detailed level. We show how such a quantitative method illuminated the mechanistic importance of a structural characteristic of multihelical TM proteins, that is, the likely presence of adjacent polar and hydrophobic residues at the protein-membrane interface. Such adjacency can preclude the complete alleviation of the well-known hydrophobic mismatch between TM proteins and the surrounding membrane, giving rise to an energy cost of residual hydrophobic mismatch. The energy cost and biophysical formulation of hydrophobic mismatch and residual hydrophobic mismatch are reviewed in the context of their mechanistic role in the function of prototypical members of multihelical TM protein families: 1), LeuT, a bacterial homolog of mammalian neurotransmitter sodium symporters; and 2), rhodopsin and the β1- and β2-adrenergic receptors from the G-protein coupled receptor family. The type of computational analysis provided by these examples is poised to translate the rapidly growing structural data for the many TM protein families that are of great importance to cell function into ever more incisive insights into mechanisms driven by protein-ligand and protein-protein interactions in the membrane environment.
Collapse
Affiliation(s)
- Sayan Mondal
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York
| | - George Khelashvili
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York
| | - Harel Weinstein
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York; The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Cornell Medical College of Cornell University, New York, New York.
| |
Collapse
|
36
|
Bakan A, Dutta A, Mao W, Liu Y, Chennubhotla C, Lezon TR, Bahar I. Evol and ProDy for bridging protein sequence evolution and structural dynamics. Bioinformatics 2014; 30:2681-3. [PMID: 24849577 DOI: 10.1093/bioinformatics/btu336] [Citation(s) in RCA: 153] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
UNLABELLED Correlations between sequence evolution and structural dynamics are of utmost importance in understanding the molecular mechanisms of function and their evolution. We have integrated Evol, a new package for fast and efficient comparative analysis of evolutionary patterns and conformational dynamics, into ProDy, a computational toolbox designed for inferring protein dynamics from experimental and theoretical data. Using information-theoretic approaches, Evol coanalyzes conservation and coevolution profiles extracted from multiple sequence alignments of protein families with their inferred dynamics. AVAILABILITY AND IMPLEMENTATION ProDy and Evol are open-source and freely available under MIT License from http://prody.csb.pitt.edu/.
Collapse
Affiliation(s)
- Ahmet Bakan
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Anindita Dutta
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Wenzhi Mao
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Ying Liu
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Chakra Chennubhotla
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Timothy R Lezon
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| | - Ivet Bahar
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
| |
Collapse
|
37
|
Gur M, Zomot E, Bahar I. Global motions exhibited by proteins in micro- to milliseconds simulations concur with anisotropic network model predictions. J Chem Phys 2014; 139:121912. [PMID: 24089724 DOI: 10.1063/1.4816375] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Anton supercomputing technology recently developed for efficient molecular dynamics simulations permits us to examine micro- to milli-second events at full atomic resolution for proteins in explicit water and lipid bilayer. It also permits us to investigate to what extent the collective motions predicted by network models (that have found broad use in molecular biophysics) agree with those exhibited by full-atomic long simulations. The present study focuses on Anton trajectories generated for two systems: the bovine pancreatic trypsin inhibitor, and an archaeal aspartate transporter, GltPh. The former, a thoroughly studied system, helps benchmark the method of comparative analysis, and the latter provides new insights into the mechanism of function of glutamate transporters. The principal modes of motion derived from both simulations closely overlap with those predicted for each system by the anisotropic network model (ANM). Notably, the ANM modes define the collective mechanisms, or the pathways on conformational energy landscape, that underlie the passage between the crystal structure and substates visited in simulations. In particular, the lowest frequency ANM modes facilitate the conversion between the most probable substates, lending support to the view that easy access to functional substates is a robust determinant of evolutionarily selected native contact topology.
Collapse
Affiliation(s)
- M Gur
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, 3501 Fifth Ave, Suite 3064 BST3, Pittsburgh, Pennsylvania 15260, USA
| | | | | |
Collapse
|
38
|
Bastolla U. Computing protein dynamics from protein structure with elastic network models. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1186] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Ugo Bastolla
- Centro de Biologa Molecular Severo Ochoa (CSIC‐UAM)Universidad Autónoma de MadridMadridSpain
| |
Collapse
|
39
|
Das A, Gur M, Cheng MH, Jo S, Bahar I, Roux B. Exploring the conformational transitions of biomolecular systems using a simple two-state anisotropic network model. PLoS Comput Biol 2014; 10:e1003521. [PMID: 24699246 PMCID: PMC3974643 DOI: 10.1371/journal.pcbi.1003521] [Citation(s) in RCA: 96] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2013] [Accepted: 02/01/2014] [Indexed: 11/19/2022] Open
Abstract
Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results. Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.
Collapse
Affiliation(s)
- Avisek Das
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America
| | - Mert Gur
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Mary Hongying Cheng
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Sunhwan Jo
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America
| | - Ivet Bahar
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America
- * E-mail:
| |
Collapse
|
40
|
Neurotransmitter transporters: structure meets function. Structure 2014; 21:694-705. [PMID: 23664361 DOI: 10.1016/j.str.2013.03.002] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2012] [Revised: 02/15/2013] [Accepted: 03/06/2013] [Indexed: 12/22/2022]
Abstract
At synapses, sodium-coupled transporters remove released neurotransmitters, thereby recycling them and maintaining a low extracellular concentration of the neurotransmitter. The molecular mechanism underlying sodium-coupled neurotransmitter uptake is not completely understood. Several structures of homologs of human neurotransmitter transporters have been solved with X-ray crystallography. These crystal structures have spurred a plethora of computational and experimental work to elucidate the molecular mechanism underlying sodium-coupled transport. Here, we compare the structures of GltPh, a glutamate transporter homolog, and LeuT, a homolog of neurotransmitter transporters for the biogenic amines and inhibitory molecules GABA and glycine. We relate these structures to data obtained from experiments and computational simulations, to draw conclusions about the mechanism of uptake by sodium-coupled neurotransmitter transporters. Here, we propose how sodium and substrate binding is coupled and how binding of sodium and substrate opens and closes the gates in these transporters, thereby leading to an efficient coupled transport.
Collapse
|
41
|
Rauen T, Tanui R, Grewer C. Structural and functional dynamics of Excitatory Amino Acid Transporters (EAAT). AIMS MOLECULAR SCIENCE 2014. [DOI: 10.3934/molsci.2014.3.99] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
42
|
Ozcan O, Uyar A, Doruker P, Akten ED. Effect of intracellular loop 3 on intrinsic dynamics of human β2-adrenergic receptor. BMC STRUCTURAL BIOLOGY 2013; 13:29. [PMID: 24206668 PMCID: PMC3834532 DOI: 10.1186/1472-6807-13-29] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/03/2013] [Accepted: 11/06/2013] [Indexed: 11/10/2022]
Abstract
Background To understand the effect of the long intracellular loop 3 (ICL3) on the intrinsic dynamics of human β2-adrenergic receptor, molecular dynamics (MD) simulations were performed on two different models, both of which were based on the inactive crystal structure in complex with carazolol (after removal of carazolol and T4-lysozyme). In the so-called loop model, the ICL3 region that is missing in available crystal structures was modeled as an unstructured loop of 32-residues length, whereas in the clipped model, the two open ends were covalently bonded to each other. The latter model without ICL3 was taken as a reference, which has also been commonly used in recent computational studies. Each model was embedded into POPC bilayer membrane with explicit water and subjected to a 1 μs molecular dynamics (MD) simulation at 310 K. Results After around 600 ns, the loop model started a transition to a “very inactive” conformation, which is characterized by a further movement of the intracellular half of transmembrane helix 6 (TM6) towards the receptor core, and a close packing of ICL3 underneath the membrane completely blocking the G-protein’s binding site. Concurrently, the binding site at the extracellular part of the receptor expanded slightly with the Ser207-Asp113 distance increasing to 18 Å from 11 Å, which was further elaborated by docking studies. Conclusions The essential dynamics analysis indicated a strong coupling between the extracellular and intracellular parts of the intact receptor, implicating a functional relevance for allosteric regulation. In contrast, no such transition to the “very inactive” state, nor any structural correlation, was observed in the clipped model without ICL3. Furthermore, elastic network analysis using different conformers for the loop model indicated a consistent picture on the specific ICL3 conformational change being driven by global modes.
Collapse
Affiliation(s)
| | | | | | - Ebru Demet Akten
- Department of Bioinformatics and Genetics, Faculty of Natural Sciences and Engineering, Kadir Has University, Cibali 34083, Istanbul, Turkey.
| |
Collapse
|
43
|
Lee CH, Chen YW, Huang YT, Pan YJ, Lee CH, Lin SM, Huang LK, Lo YY, Huang YF, Hsu YD, Yen SC, Hwang JK, Pan RL. Functional Investigation of Transmembrane Helix 3 in H+-Translocating Pyrophosphatase. J Membr Biol 2013; 246:959-66. [DOI: 10.1007/s00232-013-9599-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Accepted: 09/24/2013] [Indexed: 10/26/2022]
|
44
|
Protein mechanics: how force regulates molecular function. Biochim Biophys Acta Gen Subj 2013; 1830:4762-8. [PMID: 23791949 DOI: 10.1016/j.bbagen.2013.06.005] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2012] [Revised: 05/26/2013] [Accepted: 06/04/2013] [Indexed: 12/13/2022]
Abstract
BACKGROUND Regulation of proteins is ubiquitous and vital for any organism. Protein activity can be altered chemically, by covalent modifications or non-covalent binding of co-factors. Mechanical forces are emerging as an additional way of regulating proteins, by inducing a conformational change or by partial unfolding. SCOPE We review some advances in experimental and theoretical techniques to study protein allostery driven by mechanical forces, as opposed to the more conventional ligand driven allostery. In this respect, we discuss recent single molecule pulling experiments as they have substantially augmented our view on the protein allostery by mechanical signals in recent years. Finally, we present a computational analysis technique, Force Distribution Analysis, that we developed to reveal allosteric pathways in proteins. MAJOR CONCLUSIONS Any kind of external perturbation, being it ligand binding or mechanical stretching, can be viewed as an external force acting on the macromolecule, rendering force-based experimental or computational techniques, a very general approach to the mechanics involved in protein allostery. GENERAL SIGNIFICANCE This unifying view might aid to decipher how complex allosteric protein machineries are regulated on the single molecular level.
Collapse
|
45
|
Structure of a bacterial energy-coupling factor transporter. Nature 2013; 497:272-6. [DOI: 10.1038/nature12045] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2012] [Accepted: 02/26/2013] [Indexed: 01/11/2023]
|
46
|
Structure and mechanism of a nitrate transporter. Cell Rep 2013; 3:716-23. [PMID: 23523348 DOI: 10.1016/j.celrep.2013.03.007] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2012] [Revised: 01/20/2013] [Accepted: 03/06/2013] [Indexed: 11/23/2022] Open
Abstract
The nitrate/nitrite transporters NarK and NarU play an important role in nitrogen homeostasis in bacteria and belong to the nitrate/nitrite porter family (NNP) of the major facilitator superfamily (MFS) fold. The structure and functional mechanism of NarK and NarU remain unknown. Here, we report the crystal structure of NarU at a resolution of 3.1 Å and systematic biochemical characterization. The two molecules of NarU in an asymmetric unit exhibit two distinct conformational states: occluded and partially inward-open. The substrate molecule nitrate appears to be coordinated by four highly conserved, charged, or polar amino acids. Structural and biochemical analyses allowed the identification of key amino acids that are involved in substrate gating and transport. The observed conformational differences of NarU, together with unique sequence features of the NNP family transporters, suggest a transport mechanism that might deviate from the canonical rocker-switch model.
Collapse
|