1
|
Xu Y, Das P, McCord RP, Shen T. Node features of chromosome structure networks and their connections to genome annotation. Comput Struct Biotechnol J 2024; 23:2240-2250. [PMID: 38827231 PMCID: PMC11140560 DOI: 10.1016/j.csbj.2024.05.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Revised: 05/14/2024] [Accepted: 05/15/2024] [Indexed: 06/04/2024] Open
Abstract
The 3D conformations of chromosomes can encode biological significance, and the implications of such structures have been increasingly appreciated recently. Certain chromosome structural features, such as A/B compartmentalization, are frequently extracted from Hi-C pairwise genome contact information (physical association between different regions of the genome) and compared with linear annotations of the genome, such as histone modifications and lamina association. We investigate how additional properties of chromosome structure can be deduced using an abstract graph representation of the contact heatmap, and describe specific network properties that can have a strong connection with some of these biological annotations. We constructed chromosome structure networks (CSNs) from bulk Hi-C data and calculated a set of site-resolved (node-based) network properties. These properties are useful for characterizing certain aspects of chromosomal structure. We examined the ability of network properties to differentiate several scenarios, such as haploid vs diploid cells, partially inverted nuclei vs conventional architecture, depletion of chromosome architectural proteins, and structural changes during cell development. We also examined the connection between network properties and a series of other linear annotations, such as histone modifications and chromatin states including poised promoter and enhancer labels. We found that semi-local network properties exhibit greater capability in characterizing genome annotations compared to diffusive or ultra-local node features. For example, the local square clustering coefficient can be a strong classifier of lamina-associated domains. We demonstrated that network properties can be useful for highlighting large-scale chromosome structure differences that emerge in different biological situations.
Collapse
Affiliation(s)
- Yingjie Xu
- Graduate School of Genome Science & Technology, University of Tennessee, Knoxville, TN 37996, USA
| | - Priyojit Das
- Department of Genetics, Harvard Medical School, Boston, MA 02115, USA
- Department of Molecular Biology, Massachusetts General Hospital, Boston, MA 02114, USA
| | - Rachel Patton McCord
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Tongye Shen
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| |
Collapse
|
2
|
Inan T, Yuce M, MacKerell AD, Kurkcuoglu O. Exploring Druggable Binding Sites on the Class A GPCRs Using the Residue Interaction Network and Site Identification by Ligand Competitive Saturation. ACS OMEGA 2024; 9:40154-40171. [PMID: 39346853 PMCID: PMC11425613 DOI: 10.1021/acsomega.4c06172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/03/2024] [Revised: 09/04/2024] [Accepted: 09/06/2024] [Indexed: 10/01/2024]
Abstract
G protein-coupled receptors (GPCRs) play a central role in cellular signaling and are linked to many diseases. Accordingly, computational methods to explore potential allosteric sites for this class of proteins to facilitate the identification of potential modulators are needed. Importantly, the availability of rich structural data providing the locations of the orthosteric ligands and allosteric modulators targeting different GPCRs allows for the validation of approaches to identify new allosteric binding sites. Here, we validate the combination of two computational techniques, the residue interaction network (RIN) model and the site identification by ligand competitive saturation (SILCS) method, to predict putative allosteric binding sites of class A GPCRs. RIN analysis identifies hub residues that mediate allosteric signaling within a receptor and have a high capacity to alter receptor dynamics upon ligand binding. The known orthosteric (and allosteric) binding sites of 18 distinct class A GPCRs were successfully predicted by RIN through a dataset of 105 crystal structures (91 ligand-bound, 14 unbound) with up to 77.8% (76.9%) sensitivity, 92.5% (95.3%) specificity, 51.9% (50%) precision, and 86.2% (92.4%) accuracy based on the experimental and theoretical binding site data. Moreover, graph spectral analysis of the residue networks revealed that the proposed sites were located at the interfaces of highly interconnected residue clusters with a high ability to coordinate the functional dynamics. Then, we employed the SILCS-Hotspots method to assess the druggability of the novel sites predicted for 7 distinct class A GPCRs that are critical for a variety of diseases. While the known orthosteric and allosteric binding sites are successfully explored by our approach, numerous putative allosteric sites with the potential to bind drug-like molecules are proposed. The computational approach presented here promises to be a highly effective tool to predict putative allosteric sites of GPCRs to facilitate the design of effective modulators.
Collapse
Affiliation(s)
- Tugce Inan
- Department of Chemical Engineering, Istanbul Technical University, Istanbul 34469, Turkey
- Chemical Engineering Department, Faculty of Engineering & Architecture, Istanbul Beykent University, Istanbul 34396, Turkey
| | - Merve Yuce
- Department of Chemical Engineering, Istanbul Technical University, Istanbul 34469, Turkey
| | - Alexander D MacKerell
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, United States
| | - Ozge Kurkcuoglu
- Department of Chemical Engineering, Istanbul Technical University, Istanbul 34469, Turkey
| |
Collapse
|
3
|
Konovalov KA, Wu CG, Qiu Y, Balakrishnan VK, Parihar PS, O’Connor MS, Xing Y, Huang X. Disease mutations and phosphorylation alter the allosteric pathways involved in autoinhibition of protein phosphatase 2A. J Chem Phys 2023; 158:215101. [PMID: 37260014 PMCID: PMC10238128 DOI: 10.1063/5.0150272] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Accepted: 05/16/2023] [Indexed: 06/02/2023] Open
Abstract
Mutations in protein phosphatase 2A (PP2A) are connected to intellectual disability and cancer. It has been hypothesized that these mutations might disrupt the autoinhibition and phosphorylation-induced activation of PP2A. Since they are located far from both the active and substrate binding sites, it is unclear how they exert their effect. We performed allosteric pathway analysis based on molecular dynamics simulations and combined it with biochemical experiments to investigate the autoinhibition of PP2A. In the wild type (WT), the C-arm of the regulatory subunit B56δ obstructs the active and substrate binding sites exerting a dual autoinhibition effect. We find that the disease mutant, E198K, severely weakens the allosteric pathways that stabilize the C-arm in the WT. Instead, the strongest allosteric pathways in E198K take a different route that promotes exposure of the substrate binding site. To facilitate the allosteric pathway analysis, we introduce a path clustering algorithm for lumping pathways into channels. We reveal remarkable similarities between the allosteric channels of E198K and those in phosphorylation-activated WT, suggesting that the autoinhibition can be alleviated through a conserved mechanism. In contrast, we find that another disease mutant, E200K, which is in spatial proximity of E198, does not repartition the allosteric pathways leading to the substrate binding site; however, it may still induce exposure of the active site. This finding agrees with our biochemical data, allowing us to predict the activity of PP2A with the phosphorylated B56δ and provide insight into how disease mutations in spatial proximity alter the enzymatic activity in surprisingly different mechanisms.
Collapse
Affiliation(s)
- Kirill A. Konovalov
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | | | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Vijaya Kumar Balakrishnan
- McArdle Laboratory for Cancer Research, Department of Oncology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin 53705, USA
| | - Pankaj Singh Parihar
- McArdle Laboratory for Cancer Research, Department of Oncology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin 53705, USA
| | - Michael S. O’Connor
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yongna Xing
- Authors to whom correspondence should be addressed: and
| | - Xuhui Huang
- Authors to whom correspondence should be addressed: and
| |
Collapse
|
4
|
Summers TJ, Hemmati R, Miller JE, Agbaglo DA, Cheng Q, DeYonker NJ. Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase. J Chem Phys 2023; 158:065101. [PMID: 36792523 DOI: 10.1063/5.0127106] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.
Collapse
Affiliation(s)
- Thomas J Summers
- Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, USA
| | - Reza Hemmati
- Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, USA
| | - Justin E Miller
- Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, USA
| | - Donatus A Agbaglo
- Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, USA
| | - Qianyi Cheng
- Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, USA
| | - Nathan J DeYonker
- Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, USA
| |
Collapse
|
5
|
Yao XQ, Hamelberg D. From Distinct to Differential Conformational Dynamics to Map Allosteric Communication Pathways in Proteins. J Phys Chem B 2022; 126:2612-2620. [PMID: 35319195 DOI: 10.1021/acs.jpcb.2c00199] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Initiation of biological processes involving protein-ligand binding, transient protein-protein interactions, or amino acid modifications alters the conformational dynamics of proteins. Accompanying these biological processes are ensuing coupled atomic level conformational changes within the proteins. These conformational changes collectively connect multiple amino acid residues at distal allosteric, binding, and/or active sites. Local changes due to, for example, binding of a regulatory ligand at an allosteric site initiate the allosteric regulation. The allosteric signal propagates throughout the protein structure, causing changes at distal sites, activating, deactivating, or modifying the function of the protein. Hence, dynamical responses within protein structures to stimuli contain critical information on protein function. In this Perspective, we examine the description of allosteric regulation from protein dynamical responses and associated alternative and emerging computational approaches to map allosteric communication pathways between distal sites in proteins at the atomic level.
Collapse
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
6
|
Foutch D, Pham B, Shen T. Protein conformational switch discerned via network centrality properties. Comput Struct Biotechnol J 2021; 19:3599-3608. [PMID: 34257839 PMCID: PMC8246261 DOI: 10.1016/j.csbj.2021.06.004] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Revised: 06/01/2021] [Accepted: 06/02/2021] [Indexed: 11/17/2022] Open
Abstract
Network analysis has emerged as a powerful tool for examining structural biology systems. The spatial organization of the components of a biomolecular structure has been rendered as a graph representation and analyses have been performed to deduce the biophysical and mechanistic properties of these components. For proteins, the analysis of protein structure networks (PSNs), especially via network centrality measurements and cluster coefficients, has led to identifying amino acid residues that play key functional roles and classifying amino acid residues in general. Whether these network properties examined in various studies are sensitive to subtle (yet biologically significant) conformational changes remained to be addressed. Here, we focused on four types of network centrality properties (betweenness, closeness, degree, and eigenvector centralities) for conformational changes upon ligand binding of a sensor protein (constitutive androstane receptor) and an allosteric enzyme (ribonucleotide reductase). We found that eigenvector centrality is sensitive and can distinguish salient structural features between protein conformational states while other centrality measures, especially closeness centrality, are less sensitive and rather generic with respect to the structural specificity. We also demonstrated that an ensemble-informed, modified PSN with static edges removed (which we term PSN*) has enhanced sensitivity at discerning structural changes.
Collapse
Affiliation(s)
- David Foutch
- Department of Biomedical Informatics, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Bill Pham
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Tongye Shen
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA.,UT-ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| |
Collapse
|
7
|
Lindsay RJ, Mansbach RA, Gnanakaran S, Shen T. Effects of pH on an IDP conformational ensemble explored by molecular dynamics simulation. Biophys Chem 2021; 271:106552. [PMID: 33581430 PMCID: PMC8024028 DOI: 10.1016/j.bpc.2021.106552] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 01/15/2021] [Accepted: 01/20/2021] [Indexed: 01/03/2023]
Abstract
The conformational ensemble of intrinsically disordered proteins, such as α-synuclein, are responsible for their function and malfunction. Misfolding of α-synuclein can lead to neurodegenerative diseases, and the ability to study their conformations and those of other intrinsically disordered proteins under varying physiological conditions can be crucial to understanding and preventing pathologies. In contrast to well-folded peptides, a consensus feature of IDPs is their low hydropathy and high charge, which makes their conformations sensitive to pH perturbation. We examine a prominent member of this subset of IDPs, α-synuclein, using a divide-and-conquer scheme that provides enhanced sampling of IDP structural ensembles. We constructed conformational ensembles of α-synuclein under neutral (pH ~ 7) and low (pH ~ 3) pH conditions and compared our results with available information obtained from smFRET, SAXS, and NMR studies. Specifically, α-synuclein has been found to in a more compact state at low pH conditions and the structural changes observed are consistent with those from experiments. We also characterize the conformational and dynamic differences between these ensembles and discussed the implication on promoting pathogenic fibril formation. We find that under low pH conditions, neutralization of negatively charged residues leads to compaction of the C-terminal portion of α-synuclein while internal reorganization allows α-synuclein to maintain its overall end-to-end distance. We also observe different levels of intra-protein interaction between three regions of α-synuclein at varying pH and a shift towards more hydrophilic interactions with decreasing pH.
Collapse
Affiliation(s)
- Richard J Lindsay
- UT- ORNL Graduate School of Genome Science and Technology, Knoxville, TN, 37996, USA.
| | - Rachael A Mansbach
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM, 87544, USA; Department of Physics, Concordia University, Montreal, Quebec, Canada.
| | - S Gnanakaran
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM, 87544, USA.
| | - Tongye Shen
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN, 37996, USA.
| |
Collapse
|
8
|
Summers TJ, Daniel BP, Cheng Q, DeYonker NJ. Quantifying Inter-Residue Contacts through Interaction Energies. J Chem Inf Model 2019; 59:5034-5044. [DOI: 10.1021/acs.jcim.9b00804] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Thomas J. Summers
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Baty P. Daniel
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Qianyi Cheng
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Nathan J. DeYonker
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| |
Collapse
|
9
|
Nussinov R, Tsai CJ, Jang H. Does Ras Activate Raf and PI3K Allosterically? Front Oncol 2019; 9:1231. [PMID: 31799192 PMCID: PMC6874141 DOI: 10.3389/fonc.2019.01231] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Accepted: 10/28/2019] [Indexed: 12/11/2022] Open
Abstract
The mechanism through which oncogenic Ras activates its effectors is vastly important to resolve. If allostery is at play, then targeting allosteric pathways could help in quelling activation of MAPK (Raf/MEK/ERK) and PI3K (PI3K/Akt/mTOR) cell proliferation pathways. On the face of it, allosteric activation is reasonable: Ras binding perturbs the conformational ensembles of its effectors. Here, however, we suggest that at least for Raf, PI3K, and NORE1A (RASSF5), that is unlikely. Raf's long disordered linker dampens effective allosteric activation. Instead, we suggest that the high-affinity Ras–Raf binding relieves Raf's autoinhibition, shifting Raf's ensemble from the inactive to the nanocluster-mediated dimerized active state, as Ras also does for NORE1A. PI3K is recruited and allosterically activated by RTK (e.g., EGFR) at the membrane. Ras restrains PI3K's distribution and active site orientation. It stabilizes and facilitates PIP2 binding at the active site and increases the PI3K residence time at the membrane. Thus, RTKs allosterically activate PI3Kα; however, merging their action with Ras accomplishes full activation. Here we review their activation mechanisms in this light and draw attention to implications for their pharmacology.
Collapse
Affiliation(s)
- Ruth Nussinov
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, MD, United States.,Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Chung-Jung Tsai
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, MD, United States
| | - Hyunbum Jang
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, MD, United States
| |
Collapse
|
10
|
A framework for measuring association of random vectors via collapsed random variables. J MULTIVARIATE ANAL 2019. [DOI: 10.1016/j.jmva.2019.02.012] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
11
|
Smith IN, Thacker S, Seyfi M, Cheng F, Eng C. Conformational Dynamics and Allosteric Regulation Landscapes of Germline PTEN Mutations Associated with Autism Compared to Those Associated with Cancer. Am J Hum Genet 2019; 104:861-878. [PMID: 31006514 PMCID: PMC6506791 DOI: 10.1016/j.ajhg.2019.03.009] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2018] [Accepted: 03/08/2019] [Indexed: 01/07/2023] Open
Abstract
Individuals with germline PTEN tumor-suppressor variants have PTEN hamartoma tumor syndrome (PHTS). Clinically, PHTS has variable presentations; there are distinct subsets of PHTS-affected individuals, such as those diagnosed with autism spectrum disorder (ASD) or cancer. It remains unclear why mutations in one gene can lead to such seemingly disparate phenotypes. Therefore, we sought to determine whether it is possible to predict a given PHTS-affected individual's a priori risk of ASD, cancer, or the co-occurrence of both phenotypes. By integrating network proximity analysis performed on the human interactome, molecular simulations, and residue-interaction networks, we demonstrate the role of conformational dynamics in the structural communication and long-range allosteric regulation of germline PTEN variants associated with ASD or cancer. We show that the PTEN interactome shares significant overlap with the ASD and cancer interactomes, providing network-based evidence that PTEN is a crucial player in the biology of both disorders. Importantly, this finding suggests that a germline PTEN variant might perturb the ASD or cancer networks differently, thus favoring one disease outcome at any one time. Furthermore, protein-dynamic structural-network analysis reveals small-world structural communication mediated by highly conserved functional residues and potential allosteric regulation of PTEN. We identified a salient structural-communication pathway that extends across the inter-domain interface for cancer-only mutations. In contrast, the structural-communication pathway is predominantly restricted to the phosphatase domain for ASD-only mutations. Our integrative approach supports the prediction and potential modulation of the relevant conformational states that influence structural communication and long-range perturbations associated with mutational effects that lead to PTEN-ASD or PTEN-cancer phenotypes.
Collapse
Affiliation(s)
- Iris Nira Smith
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Stetson Thacker
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA
| | - Marilyn Seyfi
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Feixiong Cheng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA; Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Charis Eng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA; Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA; Taussig Cancer Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Department of Genetics and Genome Sciences, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA.
| |
Collapse
|
12
|
Hodges M, Barahona M, Yaliraki SN. Allostery and cooperativity in multimeric proteins: bond-to-bond propensities in ATCase. Sci Rep 2018; 8:11079. [PMID: 30038211 PMCID: PMC6056424 DOI: 10.1038/s41598-018-27992-z] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2018] [Accepted: 06/13/2018] [Indexed: 11/08/2022] Open
Abstract
Aspartate carbamoyltransferase (ATCase) is a large dodecameric enzyme with six active sites that exhibits allostery: its catalytic rate is modulated by the binding of various substrates at distal points from the active sites. A recently developed method, bond-to-bond propensity analysis, has proven capable of predicting allosteric sites in a wide range of proteins using an energy-weighted atomistic graph obtained from the protein structure and given knowledge only of the location of the active site. Bond-to-bond propensity establishes if energy fluctuations at given bonds have significant effects on any other bond in the protein, by considering their propagation through the protein graph. In this work, we use bond-to-bond propensity analysis to study different aspects of ATCase activity using three different protein structures and sources of fluctuations. First, we predict key residues and bonds involved in the transition between inactive (T) and active (R) states of ATCase by analysing allosteric substrate binding as a source of energy perturbations in the protein graph. Our computational results also indicate that the effect of multiple allosteric binding is non linear: a switching effect is observed after a particular number and arrangement of substrates is bound suggesting a form of long range communication between the distantly arranged allosteric sites. Second, cooperativity is explored by considering a bisubstrate analogue as the source of energy fluctuations at the active site, also leading to the identification of highly significant residues to the T ↔ R transition that enhance cooperativity across active sites. Finally, the inactive (T) structure is shown to exhibit a strong, non linear communication between the allosteric sites and the interface between catalytic subunits, rather than the active site. Bond-to-bond propensity thus offers an alternative route to explain allosteric and cooperative effects in terms of detailed atomistic changes to individual bonds within the protein, rather than through phenomenological, global thermodynamic arguments.
Collapse
Affiliation(s)
- Maxwell Hodges
- Department of Chemistry, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
- Institute of Chemical Biology, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
| | - Mauricio Barahona
- Department of Mathematics, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
- Institute of Chemical Biology, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom
| | - Sophia N Yaliraki
- Department of Chemistry, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom.
- Institute of Chemical Biology, Imperial College London, South Kensington Campus, London, SW7 2AZ, United Kingdom.
| |
Collapse
|
13
|
Yao XQ, Momin M, Hamelberg D. Elucidating Allosteric Communications in Proteins with Difference Contact Network Analysis. J Chem Inf Model 2018; 58:1325-1330. [DOI: 10.1021/acs.jcim.8b00250] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Mohamed Momin
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
14
|
Dutta S, Eckmann JP, Libchaber A, Tlusty T. Green function of correlated genes in a minimal mechanical model of protein evolution. Proc Natl Acad Sci U S A 2018; 115:E4559-E4568. [PMID: 29712824 PMCID: PMC5960285 DOI: 10.1073/pnas.1716215115] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The function of proteins arises from cooperative interactions and rearrangements of their amino acids, which exhibit large-scale dynamical modes. Long-range correlations have also been revealed in protein sequences, and this has motivated the search for physical links between the observed genetic and dynamic cooperativity. We outline here a simplified theory of protein, which relates sequence correlations to physical interactions and to the emergence of mechanical function. Our protein is modeled as a strongly coupled amino acid network with interactions and motions that are captured by the mechanical propagator, the Green function. The propagator describes how the gene determines the connectivity of the amino acids and thereby, the transmission of forces. Mutations introduce localized perturbations to the propagator that scatter the force field. The emergence of function is manifested by a topological transition when a band of such perturbations divides the protein into subdomains. We find that epistasis-the interaction among mutations in the gene-is related to the nonlinearity of the Green function, which can be interpreted as a sum over multiple scattering paths. We apply this mechanical framework to simulations of protein evolution and observe long-range epistasis, which facilitates collective functional modes.
Collapse
Affiliation(s)
- Sandipan Dutta
- Center for Soft and Living Matter, Institute for Basic Science, Ulsan 44919, Korea
| | - Jean-Pierre Eckmann
- Département de Physique Théorique and Section de Mathématiques, Université de Genève, CH-1211 Geneva 4, Switzerland
| | - Albert Libchaber
- Center for Studies in Physics and Biology, The Rockefeller University, New York, NY 10021;
| | - Tsvi Tlusty
- Center for Soft and Living Matter, Institute for Basic Science, Ulsan 44919, Korea;
- Department of Physics, Ulsan National Institute of Science and Technology, Ulsan 44919, Korea
| |
Collapse
|
15
|
Self-organization, entropy and allostery. Biochem Soc Trans 2018; 46:587-597. [PMID: 29678954 DOI: 10.1042/bst20160144] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Revised: 03/14/2018] [Accepted: 03/19/2018] [Indexed: 12/12/2022]
Abstract
Allostery is a fundamental regulatory mechanism in biology. Although generally accepted that it is a dynamics-driven process, the exact molecular mechanism of allosteric signal transmission is hotly debated. We argue that allostery is as a part of a bigger picture that also includes fractal-like properties of protein interior, hierarchical protein folding and entropy-driven molecular recognition. Although so far all these phenomena were studied separately, they stem from the same common root: self-organization of polypeptide chains and, thus, has to be studied collectively. This merge will allow the cross-referencing of a broad spectrum of multi-disciplinary data facilitating progress in all these fields.
Collapse
|
16
|
Design of Elastic Networks with Evolutionary Optimized Long-Range Communication as Mechanical Models of Allosteric Proteins. Biophys J 2017; 113:558-571. [PMID: 28793211 PMCID: PMC5550307 DOI: 10.1016/j.bpj.2017.06.043] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Revised: 05/25/2017] [Accepted: 06/14/2017] [Indexed: 02/08/2023] Open
Abstract
Allosteric effects often underlie the activity of proteins, and elucidating generic design aspects and functional principles unique to allosteric phenomena represent a major challenge. Here an approach consisting of the in silico design of synthetic structures, which, as the principal element of allostery, encode dynamical long-range coupling among two sites, is presented. The structures are represented by elastic networks, similar to coarse-grained models of real proteins. A strategy of evolutionary optimization was implemented to iteratively improve allosteric coupling. In the designed structures, allosteric interactions were analyzed in terms of strain propagation, and simple pathways that emerged during evolution were identified as signatures through which long-range communication was established. Moreover, robustness of allosteric performance with respect to mutations was demonstrated. As it turned out, the designed prototype structures reveal dynamical properties resembling those found in real allosteric proteins. Hence, they may serve as toy models of complex allosteric systems, such as proteins. Application of the developed modeling scheme to the allosteric transition in the myosin V molecular motor was also demonstrated.
Collapse
|
17
|
Kapetis D, Sassone J, Yang Y, Galbardi B, Xenakis MN, Westra RL, Szklarczyk R, Lindsey P, Faber CG, Gerrits M, Merkies ISJ, Dib-Hajj SD, Mantegazza M, Waxman SG, Lauria G. Network topology of NaV1.7 mutations in sodium channel-related painful disorders. BMC SYSTEMS BIOLOGY 2017; 11:28. [PMID: 28235406 PMCID: PMC5324268 DOI: 10.1186/s12918-016-0382-0] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2016] [Accepted: 12/20/2016] [Indexed: 11/30/2022]
Abstract
BACKGROUND Gain-of-function mutations in SCN9A gene that encodes the voltage-gated sodium channel NaV1.7 have been associated with a wide spectrum of painful syndromes in humans including inherited erythromelalgia, paroxysmal extreme pain disorder and small fibre neuropathy. These mutations change the biophysical properties of NaV1.7 channels leading to hyperexcitability of dorsal root ganglion nociceptors and pain symptoms. There is a need for better understanding of how gain-of-function mutations alter the atomic structure of Nav1.7. RESULTS We used homology modeling to build an atomic model of NaV1.7 and a network-based theoretical approach, which can predict interatomic interactions and connectivity arrangements, to investigate how pain-related NaV1.7 mutations may alter specific interatomic bonds and cause connectivity rearrangement, compared to benign variants and polymorphisms. For each amino acid substitution, we calculated the topological parameters betweenness centrality (B ct ), degree (D), clustering coefficient (CC ct ), closeness (C ct ), and eccentricity (E ct ), and calculated their variation (Δ value = mutant value -WT value ). Pathogenic NaV1.7 mutations showed significantly higher variation of |ΔB ct | compared to benign variants and polymorphisms. Using the cut-off value ±0.26 calculated by receiver operating curve analysis, we found that ΔB ct correctly differentiated pathogenic NaV1.7 mutations from variants not causing biophysical abnormalities (nABN) and homologous SNPs (hSNPs) with 76% sensitivity and 83% specificity. CONCLUSIONS Our in-silico analyses predict that pain-related pathogenic NaV1.7 mutations may affect the network topological properties of the protein and suggest |ΔB ct | value as a potential in-silico marker.
Collapse
Affiliation(s)
- Dimos Kapetis
- Bioinformatics Unit, IRCCS Foundation “Carlo Besta” Neurological Institute, Milan, Italy
- Neuroalgology Unit, IRCCS Foundation “Carlo Besta” Neurological Institute, Milan, Italy
| | - Jenny Sassone
- Neuroalgology Unit, IRCCS Foundation “Carlo Besta” Neurological Institute, Milan, Italy
- Present address: San Raffaele Scientific Institute and Vita-Salute University, Milan, Italy
| | - Yang Yang
- Department of Neurology, Yale University School of Medicine, New Haven, USA
- Center for Neuroscience and Regeneration Research, Yale University School of Medicine, New Haven, USA
| | - Barbara Galbardi
- Bioinformatics Unit, IRCCS Foundation “Carlo Besta” Neurological Institute, Milan, Italy
| | - Markos N. Xenakis
- Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands
- Department of Knowledge Engineering, Maastricht University, Maastricht, The Netherlands
| | - Ronald L. Westra
- Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands
- Department of Knowledge Engineering, Maastricht University, Maastricht, The Netherlands
| | - Radek Szklarczyk
- Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Patrick Lindsey
- Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Catharina G. Faber
- Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands
- Department of Neurology, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Monique Gerrits
- Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Ingemar S. J. Merkies
- Department of Neurology, Maastricht University Medical Center, Maastricht, The Netherlands
- Department of Neurology, Spaarne Hospital, Hoofddorp, The Netherlands
| | - Sulayman D. Dib-Hajj
- Department of Neurology, Yale University School of Medicine, New Haven, USA
- Center for Neuroscience and Regeneration Research, Yale University School of Medicine, New Haven, USA
| | - Massimo Mantegazza
- Laboratory of Excellence Ion Channel Science and Therapeutics, Institute of Molecular and Cellular Pharmacology, CNRS UMR7275 & University of Nice-Sophia Antipolis, Valbonne, France
| | - Stephen G. Waxman
- Department of Neurology, Yale University School of Medicine, New Haven, USA
- Center for Neuroscience and Regeneration Research, Yale University School of Medicine, New Haven, USA
| | - Giuseppe Lauria
- Neuroalgology Unit, IRCCS Foundation “Carlo Besta” Neurological Institute, Milan, Italy
| |
Collapse
|
18
|
Tamaki FK, Souza DP, Souza VP, Ikegami CM, Farah CS, Marana SR. Using the Amino Acid Network to Modulate the Hydrolytic Activity of β-Glycosidases. PLoS One 2016; 11:e0167978. [PMID: 27936116 PMCID: PMC5148593 DOI: 10.1371/journal.pone.0167978] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2016] [Accepted: 11/24/2016] [Indexed: 11/18/2022] Open
Abstract
The active site residues in GH1 β-glycosidases are compartmentalized into 3 functional regions, involved in catalysis or binding of glycone and aglycone motifs from substrate. However, it still remains unclear how residues outside the active site modulate the enzymatic activity. To tackle this question, we solved the crystal structure of the GH1 β-glycosidase from Spodoptera frugiperda (Sfβgly) to systematically map its residue contact network and correlate effects of mutations within and outside the active site. External mutations neighbouring the functional residues involved in catalysis and glycone-binding are deleterious, whereas mutations neighbouring the aglycone-binding site are less detrimental or even beneficial. The large dataset of new and previously characterized Sfβgly mutants supports that external perturbations are coherently transmitted to active site residues possibly through contacts and specifically disturb functional regions they interact to, reproducing the effects observed for direct mutations of functional residues. This allowed us to suggest that positions related to the aglycone-binding site are preferential targets for introduction of mutations aiming to further improve the hydrolytic activity of β-glycosidases.
Collapse
Affiliation(s)
- Fábio K. Tamaki
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
- * E-mail:
| | - Diorge P. Souza
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Valquiria P. Souza
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Cecilia M. Ikegami
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Chuck S. Farah
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Sandro R. Marana
- Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| |
Collapse
|
19
|
Barman A, Hamelberg D. Fe(II)/Fe(III) Redox Process Can Significantly Modulate the Conformational Dynamics and Electrostatics of Pirin in NF-κB Regulation. ACS OMEGA 2016; 1:837-842. [PMID: 31457166 PMCID: PMC6640773 DOI: 10.1021/acsomega.6b00231] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2016] [Accepted: 10/21/2016] [Indexed: 06/10/2023]
Abstract
Pirin is an iron (Fe)-dependent regulatory protein of nuclear factor κB (NF-κB) transcription factors. Binding studies have suggested that the oxidative state of iron plays a crucial role in modulating the binding of Pirin to NF-κB p65, in turn enhancing the binding of p65 to DNA. The Fe(III) form of Pirin is the active form and binds to NF-κB, whereas the Fe(II) form does not bind to NF-κB. However, the surprising consequence of a single charge perturbation in the functional modulation of NF-κB is not well understood. Here, we use quantum mechanical calculations and microsecond-long molecular dynamics simulations to explore the free-energy landscapes of the Fe(II) and Fe(III) forms of Pirin. We show that the restricted conformational space and electrostatic complementarity of the Fe(III) form of Pirin are crucial for binding and regulation of NF-κB. Our results suggest that a subtle single-electron redox trigger could significantly modulate the conformational dynamics and electrostatics of proteins in subcellular allosteric regulatory processes.
Collapse
|
20
|
Lambrughi M, De Gioia L, Gervasio FL, Lindorff-Larsen K, Nussinov R, Urani C, Bruschi M, Papaleo E. DNA-binding protects p53 from interactions with cofactors involved in transcription-independent functions. Nucleic Acids Res 2016; 44:9096-9109. [PMID: 27604871 PMCID: PMC5100575 DOI: 10.1093/nar/gkw770] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2015] [Revised: 08/19/2016] [Accepted: 08/23/2016] [Indexed: 12/15/2022] Open
Abstract
Binding-induced conformational changes of a protein at regions distant from the binding site may play crucial roles in protein function and regulation. The p53 tumour suppressor is an example of such an allosterically regulated protein. Little is known, however, about how DNA binding can affect distal sites for transcription factors. Furthermore, the molecular details of how a local perturbation is transmitted through a protein structure are generally elusive and occur on timescales hard to explore by simulations. Thus, we employed state-of-the-art enhanced sampling atomistic simulations to unveil DNA-induced effects on p53 structure and dynamics that modulate the recruitment of cofactors and the impact of phosphorylation at Ser215. We show that DNA interaction promotes a conformational change in a region 3 nm away from the DNA binding site. Specifically, binding to DNA increases the population of an occluded minor state at this distal site by more than 4-fold, whereas phosphorylation traps the protein in its major state. In the minor conformation, the interface of p53 that binds biological partners related to p53 transcription-independent functions is not accessible. Significantly, our study reveals a mechanism of DNA-mediated protection of p53 from interactions with partners involved in the p53 transcription-independent signalling. This also suggests that conformational dynamics is tightly related to p53 signalling.
Collapse
Affiliation(s)
- Matteo Lambrughi
- Computational Biology Laboratory, Unit of Statistics, Bioinformatics and Registry, Strandboulevarden 49, 2100, Copenhagen, Denmark
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Luca De Gioia
- Department of Biotechnology and Biosciences, University of Milano-Bicocca, Piazza della Scienza 2, 20126, Milan, Italy
| | - Francesco Luigi Gervasio
- Department of Chemistry and Institute of Structural and Molecular Biology, University College London, London WC1H 0AJ, UK
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Ruth Nussinov
- Cancer and Inflammation Program, Leidos Biomedical Research Inc., Frederick National laboratory, National Cancer Institute, Frederick, MD 21702, USA
- Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Chiara Urani
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 1, 20126, Milan, Italy
| | - Maurizio Bruschi
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 1, 20126, Milan, Italy
| | - Elena Papaleo
- Computational Biology Laboratory, Unit of Statistics, Bioinformatics and Registry, Strandboulevarden 49, 2100, Copenhagen, Denmark
- Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
21
|
Huculeci R, Cilia E, Lyczek A, Buts L, Houben K, Seeliger MA, van Nuland N, Lenaerts T. Dynamically Coupled Residues within the SH2 Domain of FYN Are Key to Unlocking Its Activity. Structure 2016; 24:1947-1959. [PMID: 27692963 DOI: 10.1016/j.str.2016.08.016] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Revised: 07/13/2016] [Accepted: 08/10/2016] [Indexed: 12/30/2022]
Abstract
Src kinase activity is controlled by various mechanisms involving a coordinated movement of kinase and regulatory domains. Notwithstanding the extensive knowledge related to the backbone dynamics, little is known about the more subtle side-chain dynamics within the regulatory domains and their role in the activation process. Here, we show through experimental methyl dynamic results and predicted changes in side-chain conformational couplings that the SH2 structure of Fyn contains a dynamic network capable of propagating binding information. We reveal that binding the phosphorylated tail of Fyn perturbs a residue cluster near the linker connecting the SH2 and SH3 domains of Fyn, which is known to be relevant in the regulation of the activity of Fyn. Biochemical perturbation experiments validate that those residues are essential for inhibition of Fyn, leading to a gain of function upon mutation. These findings reveal how side-chain dynamics may facilitate the allosteric regulation of the different members of the Src kinase family.
Collapse
Affiliation(s)
- Radu Huculeci
- Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussel, Belgium; Structural Biology Research Center, VIB, Pleinlaan 2, 1050 Brussel, Belgium
| | - Elisa Cilia
- MLG, Départment d'Informatique, Université Libre de Bruxelles, Boulevard du Triomphe CP212, 1050 Brussels, Belgium; Interuniversity Institute of Bioinformatics Brussels (IB(2)), ULB-VUB, La Plaine Campus, Boulevard du Triomphe CP 263, 1050 Brussels, Belgium
| | - Agatha Lyczek
- Department of Pharmacological Sciences, Stony Brook University School of Medicine, BST 8-140, Stony Brook, NY 11794-8651, USA
| | - Lieven Buts
- Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussel, Belgium; Structural Biology Research Center, VIB, Pleinlaan 2, 1050 Brussel, Belgium
| | - Klaartje Houben
- NMR spectroscopy, Bijvoet Center for Biomolecular Research, Utrecht University, Padualaan 8, 3584 CH Utrecht, the Netherlands
| | - Markus A Seeliger
- Department of Pharmacological Sciences, Stony Brook University School of Medicine, BST 8-140, Stony Brook, NY 11794-8651, USA
| | - Nico van Nuland
- Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussel, Belgium; Structural Biology Research Center, VIB, Pleinlaan 2, 1050 Brussel, Belgium.
| | - Tom Lenaerts
- MLG, Départment d'Informatique, Université Libre de Bruxelles, Boulevard du Triomphe CP212, 1050 Brussels, Belgium; Interuniversity Institute of Bioinformatics Brussels (IB(2)), ULB-VUB, La Plaine Campus, Boulevard du Triomphe CP 263, 1050 Brussels, Belgium; AI-lab, Vakgroep Computerwetenschappen, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium.
| |
Collapse
|
22
|
Strain analysis of protein structures and low dimensionality of mechanical allosteric couplings. Proc Natl Acad Sci U S A 2016; 113:E5847-E5855. [PMID: 27655887 DOI: 10.1073/pnas.1609462113] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
In many proteins, especially allosteric proteins that communicate regulatory states from allosteric to active sites, structural deformations are functionally important. To understand these deformations, dynamical experiments are ideal but challenging. Using static structural information, although more limited than dynamical analysis, is much more accessible. Underused for protein analysis, strain is the natural quantity for studying local deformations. We calculate strain tensor fields for proteins deformed by ligands or thermal fluctuations using crystal and NMR structure ensembles. Strains-primarily shears-show deformations around binding sites. These deformations can be induced solely by ligand binding at distant allosteric sites. Shears reveal quasi-2D paths of mechanical coupling between allosteric and active sites that may constitute a widespread mechanism of allostery. We argue that strain-particularly shear-is the most appropriate quantity for analysis of local protein deformations. This analysis can reveal mechanical and biological properties of many proteins.
Collapse
|
23
|
Vaidehi N, Bhattacharya S. Allosteric communication pipelines in G-protein-coupled receptors. Curr Opin Pharmacol 2016; 30:76-83. [PMID: 27497048 DOI: 10.1016/j.coph.2016.07.010] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2016] [Revised: 07/12/2016] [Accepted: 07/14/2016] [Indexed: 11/29/2022]
Abstract
The binding of ligands to G-protein-coupled receptors (GPCRs) in the extracellular region transmits the signal to the intracellular region to initiate coupling to effector proteins. The mechanism of this allosteric communication remains largely unexplored. Knowledge of the residues involved in the pipeline of the allosteric communication from the extracellular to the intracellular region will provide means to (a) design ligands with bias in potency towards one signaling pathway over others, and (b) design allosteric modulators that show subtype selectivity in GPCRs. In this review we describe the current state of the computational methods that provide insights into the allosteric communication in GPCRs and elucidate how this information can be used to design allosteric modulators.
Collapse
Affiliation(s)
- Nagarajan Vaidehi
- Department of Molecular Immunology, Beckman Research Institute of the City of Hope, 1500 East Duarte Road, Duarte, CA 91010, USA.
| | - Supriyo Bhattacharya
- Department of Molecular Immunology, Beckman Research Institute of the City of Hope, 1500 East Duarte Road, Duarte, CA 91010, USA
| |
Collapse
|
24
|
Westerhold LE, Adams SL, Bergman HL, Zeczycki TN. Pyruvate Occupancy in the Carboxyl Transferase Domain of Pyruvate Carboxylase Facilitates Product Release from the Biotin Carboxylase Domain through an Intermolecular Mechanism. Biochemistry 2016; 55:3447-60. [PMID: 27254467 DOI: 10.1021/acs.biochem.6b00372] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Protein structure, ligand binding, and catalytic turnover contributes to the governance of catalytic events occurring at spatially distinct domains in multifunctional enzymes. Coordination of these catalytic events partially rests on the ability of spatially discrete active sites to communicate with other allosteric and active sites on the same polypeptide chain (intramolecular) or on different polypeptide chains (intermolecular) within the holoenzyme. Often, communication results in long-range effects on substrate binding or product release. For example, pyruvate binding to the carboxyl transferase (CT) domain of pyruvate carboxylase (PC) increases the rate of product release in the biotin carboxylase (BC) domain. In order to address how CT domain ligand occupancy is "sensed" by other domains, we generated functional, mixed hybrid tetramers using the E218A (inactive BC domain) and T882S (low pyruvate binding, low activity) mutant forms of PC. The apparent Ka pyruvate for the pyruvate-stimulated release of Pi catalyzed by the T882S:E218A[1:1] hybrid tetramer was comparable to the wild-type enzyme and nearly 10-fold lower than that for the T882S homotetramer. In addition, the ratio of the rates of oxaloacetate formation to Pi release for the WT:T882S[1:1] and E218A:T882S[1:1] hybrid tetramer-catalyzed reactions was 0.5 and 0.6, respectively, while the T882S homotetramer exhibited a near 1:1 coupling of the two domains, suggesting that the mechanisms coordinating catalytic events is more complicated that we initially assumed. The results presented here are consistent with an intermolecular communication mechanism, where pyruvate binding to the CT domain is "sensed" by domains on a different polypeptide chain within the tetramer.
Collapse
Affiliation(s)
- Lauren E Westerhold
- Department of Biochemistry and Molecular Biology and the ‡East Carolina Diabetes and Obesity Institute, Brody School of Medicine at East Carolina University , Greenville, North Carolina 27834, United States
| | - Stephanie L Adams
- Department of Biochemistry and Molecular Biology and the ‡East Carolina Diabetes and Obesity Institute, Brody School of Medicine at East Carolina University , Greenville, North Carolina 27834, United States
| | - Hanna L Bergman
- Department of Biochemistry and Molecular Biology and the ‡East Carolina Diabetes and Obesity Institute, Brody School of Medicine at East Carolina University , Greenville, North Carolina 27834, United States
| | - Tonya N Zeczycki
- Department of Biochemistry and Molecular Biology and the ‡East Carolina Diabetes and Obesity Institute, Brody School of Medicine at East Carolina University , Greenville, North Carolina 27834, United States
| |
Collapse
|
25
|
Barman A, Hamelberg D. Coupled Dynamics and Entropic Contribution to the Allosteric Mechanism of Pin1. J Phys Chem B 2016; 120:8405-15. [DOI: 10.1021/acs.jpcb.6b02123] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Arghya Barman
- Department
of Chemistry and
the Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department
of Chemistry and
the Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
26
|
Maximova T, Moffatt R, Ma B, Nussinov R, Shehu A. Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics. PLoS Comput Biol 2016; 12:e1004619. [PMID: 27124275 PMCID: PMC4849799 DOI: 10.1371/journal.pcbi.1004619] [Citation(s) in RCA: 132] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.
Collapse
Affiliation(s)
- Tatiana Maximova
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Ryan Moffatt
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
- Department of Biongineering, George Mason University, Fairfax, Virginia, United States of America
- School of Systems Biology, George Mason University, Manassas, Virginia, United States of America
| |
Collapse
|
27
|
Affiliation(s)
- Andre A. S. T. Ribeiro
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Vanessa Ortiz
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| |
Collapse
|
28
|
Bhattacharya S, Vaidehi N. Differences in allosteric communication pipelines in the inactive and active states of a GPCR. Biophys J 2015; 107:422-434. [PMID: 25028884 DOI: 10.1016/j.bpj.2014.06.015] [Citation(s) in RCA: 86] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2014] [Revised: 06/02/2014] [Accepted: 06/09/2014] [Indexed: 01/01/2023] Open
Abstract
G-protein-coupled receptors (GPCRs) are membrane proteins that allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to proteins in the intracellular (IC) domain. However, the complete pathway of allosteric communication from the EC to the IC domain, including the role of individual amino acids in the pathway is not known. Using the correlation in torsion angle movements calculated from microseconds-long molecular-dynamics simulations, we elucidated the allosteric pathways in three different conformational states of β2-adrenergic receptor (β2AR): 1), the inverse-agonist-bound inactive state; 2), the agonist-bound intermediate state; and (3), the agonist- and G-protein-bound fully active state. The inactive state is less dynamic compared with the intermediate and active states, showing dense clusters of allosteric pathways (allosteric pipelines) connecting the EC with the IC domain. The allosteric pipelines from the EC domain to the IC domain are weakened in the intermediate state, thus decoupling the EC domain from the IC domain and making the receptor more dynamic compared with the other states. Also, the orthosteric ligand-binding site becomes the initiator region for allosteric communication in the intermediate state. This finding agrees with the paradigm that the nature of the agonist governs the specific signaling state of the receptor. These results provide an understanding of the mechanism of allosteric communication in class A GPCRs. In addition, our analysis shows that mutations that affect the ligand efficacy, but not the binding affinity, are located in the allosteric pipelines. This clarifies the role of such mutations, which has hitherto been unexplained.
Collapse
Affiliation(s)
- Supriyo Bhattacharya
- Division of Immunology, Beckman Research Institute of the City of Hope, Duarte, California
| | - Nagarajan Vaidehi
- Division of Immunology, Beckman Research Institute of the City of Hope, Duarte, California.
| |
Collapse
|
29
|
Soltan Ghoraie L, Burkowski F, Zhu M. Using kernelized partial canonical correlation analysis to study directly coupled side chains and allostery in small G proteins. Bioinformatics 2015; 31:i124-32. [PMID: 26072474 PMCID: PMC4765857 DOI: 10.1093/bioinformatics/btv241] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Motivation: Inferring structural dependencies among a protein’s side chains helps us understand their coupled motions. It is known that coupled fluctuations can reveal pathways of communication used for information propagation in a molecule. Side-chain conformations are commonly represented by multivariate angular variables, but existing partial correlation methods that can be applied to this inference task are not capable of handling multivariate angular data. We propose a novel method to infer direct couplings from this type of data, and show that this method is useful for identifying functional regions and their interactions in allosteric proteins. Results: We developed a novel extension of canonical correlation analysis (CCA), which we call ‘kernelized partial CCA’ (or simply KPCCA), and used it to infer direct couplings between side chains, while disentangling these couplings from indirect ones. Using the conformational information and fluctuations of the inactive structure alone for allosteric proteins in the Ras and other Ras-like families, our method identified allosterically important residues not only as strongly coupled ones but also in densely connected regions of the interaction graph formed by the inferred couplings. Our results were in good agreement with other empirical findings. By studying distinct members of the Ras, Rho and Rab sub-families, we show further that KPCCA was capable of inferring common allosteric characteristics in the small G protein super-family. Availability and implementation:https://github.com/lsgh/ismb15 Contact:lsoltang@uwaterloo.ca
Collapse
Affiliation(s)
- Laleh Soltan Ghoraie
- Department of Computer Science and Department of Statistics and Actuarial Science, University of Waterloo, Waterloo, ON, Canada
| | - Forbes Burkowski
- Department of Computer Science and Department of Statistics and Actuarial Science, University of Waterloo, Waterloo, ON, Canada
| | - Mu Zhu
- Department of Computer Science and Department of Statistics and Actuarial Science, University of Waterloo, Waterloo, ON, Canada
| |
Collapse
|
30
|
Soltan Ghoraie L, Burkowski F, Zhu M. Sparse networks of directly coupled, polymorphic, and functional side chains in allosteric proteins. Proteins 2015; 83:497-516. [DOI: 10.1002/prot.24752] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2014] [Revised: 12/05/2014] [Accepted: 12/13/2014] [Indexed: 02/05/2023]
Affiliation(s)
| | - Forbes Burkowski
- School of Computer Science, University of Waterloo; Waterloo Ontario Canada
| | - Mu Zhu
- Department of Statistics and Actuarial Science; University of Waterloo; Waterloo Ontario Canada
| |
Collapse
|
31
|
Farabella I, Pham T, Henderson NS, Geibel S, Phan G, Thanassi DG, Delcour AH, Waksman G, Topf M. Allosteric signalling in the outer membrane translocation domain of PapC usher. eLife 2014; 3:e03532. [PMID: 25271373 PMCID: PMC4356140 DOI: 10.7554/elife.03532] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2014] [Accepted: 09/29/2014] [Indexed: 11/13/2022] Open
Abstract
PapC ushers are outer-membrane proteins enabling assembly and secretion of P pili in uropathogenic E. coli. Their translocation domain is a large β-barrel occluded by a plug domain, which is displaced to allow the translocation of pilus subunits across the membrane. Previous studies suggested that this gating mechanism is controlled by a β-hairpin and an α-helix. To investigate the role of these elements in allosteric signal communication, we developed a method combining evolutionary and molecular dynamics studies of the native translocation domain and mutants lacking the β-hairpin and/or the α-helix. Analysis of a hybrid residue interaction network suggests distinct regions (residue 'communities') within the translocation domain (especially around β12-β14) linking these elements, thereby modulating PapC gating. Antibiotic sensitivity and electrophysiology experiments on a set of alanine-substitution mutants confirmed functional roles for four of these communities. This study illuminates the gating mechanism of PapC ushers and its importance in maintaining outer-membrane permeability.
Collapse
Affiliation(s)
- Irene Farabella
- Institute of Structural and Molecular Biology, Birkbeck College and University College London, London, United Kingdom
| | - Thieng Pham
- Department of Biology and Biochemistry, University of Houston, Houston, United States
| | - Nadine S Henderson
- Department of Molecular Genetics and Microbiology, Center for Infectious Diseases, Stony Brook University, Stony Brook, United States
| | - Sebastian Geibel
- Institute of Structural and Molecular Biology, Birkbeck College and University College London, London, United Kingdom
| | - Gilles Phan
- Institute of Structural and Molecular Biology, Birkbeck College and University College London, London, United Kingdom
| | - David G Thanassi
- Department of Molecular Genetics and Microbiology, Center for Infectious Diseases, Stony Brook University, Stony Brook, United States
| | - Anne H Delcour
- Department of Biology and Biochemistry, University of Houston, Houston, United States
| | - Gabriel Waksman
- Institute of Structural and Molecular Biology, Birkbeck College and University College London, London, United Kingdom
| | - Maya Topf
- Institute of Structural and Molecular Biology, Birkbeck College and University College London, London, United Kingdom
| |
Collapse
|
32
|
Engineering allostery. Trends Genet 2014; 30:521-8. [PMID: 25306102 DOI: 10.1016/j.tig.2014.09.004] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2014] [Revised: 09/04/2014] [Accepted: 09/05/2014] [Indexed: 02/04/2023]
Abstract
Allosteric proteins have great potential in synthetic biology, but our limited understanding of the molecular underpinnings of allostery has hindered the development of designer molecules, including transcription factors with new DNA-binding or ligand-binding specificities that respond appropriately to inducers. Such allosteric proteins could function as novel switches in complex circuits, metabolite sensors, or as orthogonal regulators for independent, inducible control of multiple genes. Advances in DNA synthesis and next-generation sequencing technologies have enabled the assessment of millions of mutants in a single experiment, providing new opportunities to study allostery. Using the classic LacI protein as an example, we describe a genetic selection system using a bidirectional reporter to capture mutants in both allosteric states, allowing the positions most crucial for allostery to be identified. This approach is not limited to bacterial transcription factors, and could reveal new mechanistic insights and facilitate engineering of other major classes of allosteric proteins such as nuclear receptors, two-component systems, G protein-coupled receptors, and protein kinases.
Collapse
|
33
|
Tiberti M, Invernizzi G, Lambrughi M, Inbar Y, Schreiber G, Papaleo E. PyInteraph: a framework for the analysis of interaction networks in structural ensembles of proteins. J Chem Inf Model 2014; 54:1537-51. [PMID: 24702124 DOI: 10.1021/ci400639r] [Citation(s) in RCA: 75] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In the last years, a growing interest has been gathering around the ability of Molecular Dynamics (MD) to provide insight into the paths of long-range structural communication in biomolecules. The knowledge of the mechanisms related to structural communication helps in the rationalization in atomistic details of the effects induced by mutations, ligand binding, and the intrinsic dynamics of proteins. We here present PyInteraph, a tool for the analysis of structural ensembles inspired by graph theory. PyInteraph is a software suite designed to analyze MD and structural ensembles with attention to binary interactions between residues, such as hydrogen bonds, salt bridges, and hydrophobic interactions. PyInteraph also allows the different classes of intra- and intermolecular interactions to be represented, combined or alone, in the form of interaction graphs, along with performing network analysis on the resulting interaction graphs. The program also integrates the network description with a knowledge-based force field to estimate the interaction energies between side chains in the protein. It can be used alone or together with the recently developed xPyder PyMOL plugin through an xPyder-compatible format. The software capabilities and associated protocols are here illustrated by biologically relevant cases of study. The program is available free of charge as Open Source software via the GPL v3 license at http://linux.btbs.unimib.it/pyinteraph/.
Collapse
Affiliation(s)
- Matteo Tiberti
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | | | | | | | | | | |
Collapse
|
34
|
Unraveling structural mechanisms of allosteric drug action. Trends Pharmacol Sci 2014; 35:256-64. [PMID: 24742712 DOI: 10.1016/j.tips.2014.03.006] [Citation(s) in RCA: 101] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2014] [Revised: 03/03/2014] [Accepted: 03/07/2014] [Indexed: 11/20/2022]
Abstract
Orthosteric drugs block the active site to obstruct function; allosteric drugs modify the population of the active state, to modulate function. Available data lead us to propose that allosteric drugs can constitute anchors and drivers. The anchor docks into an allosteric pocket. The conformation with which it interacts is unchanged during the transition between the inactive and active states. The anchor provides the foundation that allows the driver to exert a 'pull' and/or 'push' action that shifts the receptor population from the inactive to the active state. The presence or absence of driver atom in an allosteric drug can exert opposite agonism. We map a strategy for driver identification and expect the allosteric trigger concept to transform agonist/antagonist drug discovery.
Collapse
|
35
|
Ito F, Miyake M, Fushinobu S, Nakamura S, Shimizu K, Wakagi T. Engineering the allosteric properties of archaeal non-phosphorylating glyceraldehyde-3-phosphate dehydrogenases. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2014; 1844:759-66. [DOI: 10.1016/j.bbapap.2014.01.017] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/22/2013] [Revised: 01/22/2014] [Accepted: 01/25/2014] [Indexed: 11/25/2022]
|
36
|
Feher VA, Durrant JD, Van Wart AT, Amaro RE. Computational approaches to mapping allosteric pathways. Curr Opin Struct Biol 2014; 25:98-103. [PMID: 24667124 DOI: 10.1016/j.sbi.2014.02.004] [Citation(s) in RCA: 101] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2014] [Revised: 02/20/2014] [Accepted: 02/24/2014] [Indexed: 01/17/2023]
Abstract
Allosteric signaling occurs when chemical and/or physical changes at an allosteric site alter the activity of a primary orthosteric site often many Ångströms distant. A number of recently developed computational techniques, including dynamical network analysis, novel topological and molecular dynamics methods, and hybrids of these methods, are useful for elucidating allosteric signaling pathways at the atomistic level. No single method prevails as best to identify allosteric signal propagation path(s), rather each has particular strengths in characterizing signals that occur over specific timescale ranges and magnitudes of conformational fluctuation. With continued improvement in accuracy and predictive power, these computational techniques aim to become useful drug discovery tools that will allow researchers to identify allostery critical residues for subsequent pharmacological targeting.
Collapse
Affiliation(s)
- Victoria A Feher
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, USA
| | - Jacob D Durrant
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, USA
| | - Adam T Van Wart
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, USA
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, USA.
| |
Collapse
|
37
|
Verkhivker GM. Computational Studies of Allosteric Regulation in the Hsp90 Molecular Chaperone: From Functional Dynamics and Protein Structure Networks to Allosteric Communications and Targeted Anti-Cancer Modulators. Isr J Chem 2014. [DOI: 10.1002/ijch.201300143] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
38
|
Blacklock K, Verkhivker GM. Allosteric regulation of the Hsp90 dynamics and stability by client recruiter cochaperones: protein structure network modeling. PLoS One 2014; 9:e86547. [PMID: 24466147 PMCID: PMC3896489 DOI: 10.1371/journal.pone.0086547] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2013] [Accepted: 12/06/2013] [Indexed: 12/29/2022] Open
Abstract
The fundamental role of the Hsp90 chaperone in supporting functional activity of diverse protein clients is anchored by specific cochaperones. A family of immune sensing client proteins is delivered to the Hsp90 system with the aid of cochaperones Sgt1 and Rar1 that act cooperatively with Hsp90 to form allosterically regulated dynamic complexes. In this work, functional dynamics and protein structure network modeling are combined to dissect molecular mechanisms of Hsp90 regulation by the client recruiter cochaperones. Dynamic signatures of the Hsp90-cochaperone complexes are manifested in differential modulation of the conformational mobility in the Hsp90 lid motif. Consistent with the experiments, we have determined that targeted reorganization of the lid dynamics is a unifying characteristic of the client recruiter cochaperones. Protein network analysis of the essential conformational space of the Hsp90-cochaperone motions has identified structurally stable interaction communities, interfacial hubs and key mediating residues of allosteric communication pathways that act concertedly with the shifts in conformational equilibrium. The results have shown that client recruiter cochaperones can orchestrate global changes in the dynamics and stability of the interaction networks that could enhance the ATPase activity and assist in the client recruitment. The network analysis has recapitulated a broad range of structural and mutagenesis experiments, particularly clarifying the elusive role of Rar1 as a regulator of the Hsp90 interactions and a stability enhancer of the Hsp90-cochaperone complexes. Small-world organization of the interaction networks in the Hsp90 regulatory complexes gives rise to a strong correspondence between highly connected local interfacial hubs, global mediator residues of allosteric interactions and key functional hot spots of the Hsp90 activity. We have found that cochaperone-induced conformational changes in Hsp90 may be determined by specific interaction networks that can inhibit or promote progression of the ATPase cycle and thus control the recruitment of client proteins.
Collapse
Affiliation(s)
- Kristin Blacklock
- School of Computational Sciences and Crean School of Health and Life Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
| | - Gennady M. Verkhivker
- School of Computational Sciences and Crean School of Health and Life Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
- Department of Pharmacology, University of California San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
39
|
Decipher the mechanisms of protein conformational changes induced by nucleotide binding through free-energy landscape analysis: ATP binding to Hsp70. PLoS Comput Biol 2013; 9:e1003379. [PMID: 24348227 PMCID: PMC3861046 DOI: 10.1371/journal.pcbi.1003379] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2013] [Accepted: 10/15/2013] [Indexed: 11/30/2022] Open
Abstract
ATP regulates the function of many proteins in the cell by transducing its binding and hydrolysis energies into protein conformational changes by mechanisms which are challenging to identify at the atomic scale. Based on molecular dynamics (MD) simulations, a method is proposed to analyze the structural changes induced by ATP binding to a protein by computing the effective free-energy landscape (FEL) of a subset of its coordinates along its amino-acid sequence. The method is applied to characterize the mechanism by which the binding of ATP to the nucleotide-binding domain (NBD) of Hsp70 propagates a signal to its substrate-binding domain (SBD). Unbiased MD simulations were performed for Hsp70-DnaK chaperone in nucleotide-free, ADP-bound and ATP-bound states. The simulations revealed that the SBD does not interact with the NBD for DnaK in its nucleotide-free and ADP-bound states whereas the docking of the SBD was found in the ATP-bound state. The docked state induced by ATP binding found in MD is an intermediate state between the initial nucleotide-free and final ATP-bound states of Hsp70. The analysis of the FEL projected along the amino-acid sequence permitted to identify a subset of 27 protein internal coordinates corresponding to a network of 91 key residues involved in the conformational change induced by ATP binding. Among the 91 residues, 26 are identified for the first time, whereas the others were shown relevant for the allosteric communication of Hsp70 s in several experiments and bioinformatics analysis. The FEL analysis revealed also the origin of the ATP-induced structural modifications of the SBD recently measured by Electron Paramagnetic Resonance. The pathway between the nucleotide-free and the intermediate state of DnaK was extracted by applying principal component analysis to the subset of internal coordinates describing the transition. The methodology proposed is general and could be applied to analyze allosteric communication in other proteins. The precise biophysical characterization of the mechanisms of the protein conformational changes controlled by a nucleotide remains a challenge in biology. Molecular dynamics simulations of proteins in different nucleotide-binding states contain information on the nucleotide-dependent conformational dynamics. However, it is difficult to extract relevant information about the conformation-induced mechanism from the raw molecular dynamics data. Herein, we addressed this issue for the major ATP-dependent molecular chaperones Hsp70 s, which contribute to crucial cellular processes and are involved in several neurodegenerative diseases and in cancer. To function, Hsp70 undergoes several conformational changes controlled by the state of its nucleotide-binding domain. We demonstrated that the analysis of the effective free-energy landscape of the protein projected along the amino-acid sequence and computed from the molecular dynamics simulations of Hsp70 in different nucleotide-binding states, holds the key to identify the key residues of the conformational induced pathway. Identification of the key residues involved in the propagation of the structural changes induced by ATP binding offer alternative druggable specific sites other than the ligand binding clefts. The methodology developed for Hsp70 is general and can be adapted to any ligand induced conformational change in proteins.
Collapse
|
40
|
Csermely P, Korcsmáros T, Kiss HJM, London G, Nussinov R. Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review. Pharmacol Ther 2013; 138:333-408. [PMID: 23384594 PMCID: PMC3647006 DOI: 10.1016/j.pharmthera.2013.01.016] [Citation(s) in RCA: 512] [Impact Index Per Article: 46.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2013] [Accepted: 01/22/2013] [Indexed: 02/02/2023]
Abstract
Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.
Collapse
Affiliation(s)
- Peter Csermely
- Department of Medical Chemistry, Semmelweis University, P.O. Box 260, H-1444 Budapest 8, Hungary.
| | | | | | | | | |
Collapse
|
41
|
A mechanistic understanding of allosteric immune escape pathways in the HIV-1 envelope glycoprotein. PLoS Comput Biol 2013; 9:e1003046. [PMID: 23696718 PMCID: PMC3656115 DOI: 10.1371/journal.pcbi.1003046] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2012] [Accepted: 03/15/2013] [Indexed: 11/19/2022] Open
Abstract
The HIV-1 envelope (Env) spike, which consists of a compact, heterodimeric trimer of the glycoproteins gp120 and gp41, is the target of neutralizing antibodies. However, the high mutation rate of HIV-1 and plasticity of Env facilitates viral evasion from neutralizing antibodies through various mechanisms. Mutations that are distant from the antibody binding site can lead to escape, probably by changing the conformation or dynamics of Env; however, these changes are difficult to identify and define mechanistically. Here we describe a network analysis-based approach to identify potential allosteric immune evasion mechanisms using three known HIV-1 Env gp120 protein structures from two different clades, B and C. First, correlation and principal component analyses of molecular dynamics (MD) simulations identified a high degree of long-distance coupled motions that exist between functionally distant regions within the intrinsic dynamics of the gp120 core, supporting the presence of long-distance communication in the protein. Then, by integrating MD simulations with network theory, we identified the optimal and suboptimal communication pathways and modules within the gp120 core. The results unveil both strain-dependent and -independent characteristics of the communication pathways in gp120. We show that within the context of three structurally homologous gp120 cores, the optimal pathway for communication is sequence sensitive, i.e. a suboptimal pathway in one strain becomes the optimal pathway in another strain. Yet the identification of conserved elements within these communication pathways, termed inter-modular hotspots, could present a new opportunity for immunogen design, as this could be an additional mechanism that HIV-1 uses to shield vulnerable antibody targets in Env that induce neutralizing antibody breadth.
Collapse
|
42
|
Erman B. A fast approximate method of identifying paths of allosteric communication in proteins. Proteins 2013; 81:1097-101. [PMID: 23508936 DOI: 10.1002/prot.24284] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2013] [Revised: 02/22/2013] [Accepted: 03/05/2013] [Indexed: 11/07/2022]
Abstract
Fluctuations of the distance between a pair of residues i and j may be correlated with the fluctuations of the distance between another pair k and l. In this case, information may be transmitted among these four residues. Allosteric activity is postulated to proceed through such correlated paths. In this short communication a fast method for calculating correlations among all possible pairs ij and kl leading to a pathway of correlated residues of a protein is proposed. The method is based on the alpha carbon centered Gaussian Network Model. The model is applied to Glutamine Amidotransferase and pathways of allosteric activity are identified and compared with literature.
Collapse
Affiliation(s)
- Burak Erman
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey.
| |
Collapse
|
43
|
Abstract
Allostery is a biological phenomenon of fundamental importance in regulation and signaling, and efforts to understand this process have led to the development of numerous models. In spite of individual successes in understanding the structural determinants of allostery in well-documented systems, much less success has been achieved in identifying a set of quantitative and transferable ground rules that provide an understanding of how allostery works. Are there organizing principles that allow us to relate structurally different proteins, or are the determinants of allostery unique to each system? Using an ensemble-based model, we show that allosteric phenomena can be formulated in terms of conformational free energies of the cooperative elements in a protein and the coupling interactions between them. Interestingly, the resulting allosteric ground rules provide a framework to reconcile observations that challenge purely structural models of site-to-site coupling, including (a) allostery in the absence of pathways of structural distortions, (b) allostery in the absence of any structural change, and (c) the ability of allosteric ligands to act as agonists under some circumstances and antagonists under others. The ensemble view of allostery that emerges provides insights into the energetic prerequisites of site-to-site coupling and thus into how allostery works.
Collapse
Affiliation(s)
- Vincent J Hilser
- Department of Biology, The Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | | | | |
Collapse
|
44
|
Abstract
Allostery is a biological phenomenon of critical importance in metabolic regulation and cell signalling. The fundamental premise of classical models that describe allostery is that structure mediates 'action at a distance'. Recently, this paradigm has been challenged by the enrichment of IDPs (intrinsically disordered proteins) or ID (intrinsically disordered) segments in transcription factors and signalling pathways of higher organisms, where an allosteric response from external signals is requisite for regulated function. This observation strongly suggests that IDPs elicit the capacity for finely tunable allosteric regulation. Is there a set of transferable ground rules that reconcile these disparate allosteric phenomena? We focus on findings from the human GR (glucocorticoid receptor) which is a nuclear transcription factor in the SHR (steroid hormone receptor) family. GR contains an intrinsically disordered NTD (N-terminal domain) that is obligatory for transcription activity. Different GR translational isoforms have various lengths of NTD and by studying these isoforms we found that the full-length ID NTD consists of two thermodynamically distinct coupled regions. The data are interpreted in the context of an EAM (ensemble allosteric model) that considers only the intrinsic and measurable energetics of allosteric systems. Expansion of the EAM is able to reconcile the paradox that ligands for SHRs can be agonists and antagonists in a cell-context-dependent manner. These findings suggest a mechanism by which SHRs in particular, and IDPs in general, may have evolved to couple thermodynamically distinct ID segments. The ensemble view of allostery that is illuminated provides organizing principles to unify the description of all allosteric systems and insight into 'how' allostery works.
Collapse
|
45
|
Toward an understanding of the sequence and structural basis of allosteric proteins. J Mol Graph Model 2013; 40:30-9. [PMID: 23337573 DOI: 10.1016/j.jmgm.2012.12.011] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2012] [Revised: 12/18/2012] [Accepted: 12/19/2012] [Indexed: 11/19/2022]
Abstract
Allostery is the most efficient means of regulating protein functions, ranging from the control of metabolic mechanisms to signal transduction pathways. Although allosteric regulation has been recognized for half a century, our knowledge is limited to the characteristics of allosteric proteins and the structural coupling of allosteric sites and modulators. In this paper, we present a comprehensive analysis of allosteric proteins that provides insight into the foundation of allosteric interactions by revealing a series of common features in the allosteric proteins. Allosteric proteins mainly appear in transferases, and phosphorylation is the most common type of modification found in allosteric proteins. Disorders related to allosteric proteins primarily comprise metabolic diseases and cancers. In general, allosteric proteins prefer to exist as monomers or even-numbered multimers. Greater stability and hydrophobicity are observed in allosteric proteins than in general proteins. Further analysis of the allosteric sites reveals a series of buried and compact pockets composed of significantly greater hydrophobic surface area than the corresponding orthosteric sites. The hydrophobicity of the allosteric sites plays a dominant role in the binding of allosteric modulators as observed in the analysis of 106 diverse allosteric protein-modulator pairs. These results may be of great significance in predicting which proteins are allosteric and in designing novel triggers to inhibit or activate proteins of interest.
Collapse
|
46
|
Eargle J, Luthey-Schulten Z. NetworkView: 3D display and analysis of protein·RNA interaction networks. Bioinformatics 2012; 28:3000-1. [PMID: 22982572 DOI: 10.1093/bioinformatics/bts546] [Citation(s) in RCA: 156] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
SUMMARY NetworkView is an application for the display and analysis of protein·RNA interaction networks derived from structure and/or dynamics. These networks typically model individual protein residues and nucleic acid monomers as nodes and their pairwise contacts as edges with associated weights. NetworkView projects the network onto the underlying 3D molecular structure so that visualization and analysis of the network can be coupled to physical and biological properties. NetworkView is implemented as a plugin to the molecular visualization software VMD. AVAILABILITY AND IMPLEMENTATION NetworkView is included with VMD, which is available at http://www.ks.uiuc.edu/Research/vmd/. Documentation, tutorials and supporting programs are available at http://www.scs.illinois.edu/schulten/software/. CONTACT networkview@scs.illinois.edu
Collapse
Affiliation(s)
- John Eargle
- Department of Chemistry, University of Illinois at Urbana-Champaign, 505 South Mathews Avenue Urbana, IL 61801, USA
| | | |
Collapse
|
47
|
Qi Y, Wang Q, Tang B, Lai L. Identifying Allosteric Binding Sites in Proteins with a Two-State Go̅ Model for Novel Allosteric Effector Discovery. J Chem Theory Comput 2012; 8:2962-71. [DOI: 10.1021/ct300395h] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Yifei Qi
- BNLMS,
State Key Laboratory for Structural Chemistry of Unstable and Stable
Species, and Peking−Tsinghua Center for Life Sciences at College
of Chemistry and Molecular Engineering and ‡Center for Quantitative Biology, Peking University, Beijing 100871, China
| | - Qian Wang
- BNLMS,
State Key Laboratory for Structural Chemistry of Unstable and Stable
Species, and Peking−Tsinghua Center for Life Sciences at College
of Chemistry and Molecular Engineering and ‡Center for Quantitative Biology, Peking University, Beijing 100871, China
| | - Bo Tang
- BNLMS,
State Key Laboratory for Structural Chemistry of Unstable and Stable
Species, and Peking−Tsinghua Center for Life Sciences at College
of Chemistry and Molecular Engineering and ‡Center for Quantitative Biology, Peking University, Beijing 100871, China
| | - Luhua Lai
- BNLMS,
State Key Laboratory for Structural Chemistry of Unstable and Stable
Species, and Peking−Tsinghua Center for Life Sciences at College
of Chemistry and Molecular Engineering and ‡Center for Quantitative Biology, Peking University, Beijing 100871, China
| |
Collapse
|
48
|
Rader AJ, Yennamalli RM, Harter AK, Sen TZ. A rigid network of long-range contacts increases thermostability in a mutant endoglucanase. J Biomol Struct Dyn 2012; 30:628-37. [PMID: 22731517 DOI: 10.1080/07391102.2012.689696] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35 V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35 V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.
Collapse
Affiliation(s)
- A J Rader
- Department of Physics, Indiana University-Purdue University Indianapolis, Indianapolis, IN 46202, USA.
| | | | | | | |
Collapse
|
49
|
Dixit A, Verkhivker GM. Probing molecular mechanisms of the Hsp90 chaperone: biophysical modeling identifies key regulators of functional dynamics. PLoS One 2012; 7:e37605. [PMID: 22624053 PMCID: PMC3356286 DOI: 10.1371/journal.pone.0037605] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2011] [Accepted: 04/24/2012] [Indexed: 12/26/2022] Open
Abstract
Deciphering functional mechanisms of the Hsp90 chaperone machinery is an important objective in cancer biology aiming to facilitate discovery of targeted anti-cancer therapies. Despite significant advances in understanding structure and function of molecular chaperones, organizing molecular principles that control the relationship between conformational diversity and functional mechanisms of the Hsp90 activity lack a sufficient quantitative characterization. We combined molecular dynamics simulations, principal component analysis, the energy landscape model and structure-functional analysis of Hsp90 regulatory interactions to systematically investigate functional dynamics of the molecular chaperone. This approach has identified a network of conserved regions common to the Hsp90 chaperones that could play a universal role in coordinating functional dynamics, principal collective motions and allosteric signaling of Hsp90. We have found that these functional motifs may be utilized by the molecular chaperone machinery to act collectively as central regulators of Hsp90 dynamics and activity, including the inter-domain communications, control of ATP hydrolysis, and protein client binding. These findings have provided support to a long-standing assertion that allosteric regulation and catalysis may have emerged via common evolutionary routes. The interaction networks regulating functional motions of Hsp90 may be determined by the inherent structural architecture of the molecular chaperone. At the same time, the thermodynamics-based "conformational selection" of functional states is likely to be activated based on the nature of the binding partner. This mechanistic model of Hsp90 dynamics and function is consistent with the notion that allosteric networks orchestrating cooperative protein motions can be formed by evolutionary conserved and sparsely connected residue clusters. Hence, allosteric signaling through a small network of distantly connected residue clusters may be a rather general functional requirement encoded across molecular chaperones. The obtained insights may be useful in guiding discovery of allosteric Hsp90 inhibitors targeting protein interfaces with co-chaperones and protein binding clients.
Collapse
Affiliation(s)
- Anshuman Dixit
- Department of Pharmaceutical Chemistry, School of Pharmacy, The University of Kansas, Lawrence, Kansas, United States of America
| | - Gennady M. Verkhivker
- School of Computational Sciences and Crean School of Health and Life Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
- Department of Pharmacology, University of California San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
50
|
Pandini A, Fornili A, Fraternali F, Kleinjung J. Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics. FASEB J 2012; 26:868-81. [PMID: 22071506 PMCID: PMC3290435 DOI: 10.1096/fj.11-190868] [Citation(s) in RCA: 82] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.
Collapse
Affiliation(s)
- Alessandro Pandini
- Division of Mathematical Biology, Medical Research Council National Institute for Medical Research, London, UK; ,Randall Division of Cell and Molecular Biophysics, King's College London, London, UK; and , Correspondence: Division of Mathematical Biology, MRC National Institute for Medical Research, The Ridgeway, Mill Hill, NW7 1AA London, UK. E-mail: A.P., ; J.K.,
| | - Arianna Fornili
- Randall Division of Cell and Molecular Biophysics, King's College London, London, UK; and
| | - Franca Fraternali
- Randall Division of Cell and Molecular Biophysics, King's College London, London, UK; and ,The Thomas Young Centre for Theory and Simulation of Materials, London, UK
| | - Jens Kleinjung
- Division of Mathematical Biology, Medical Research Council National Institute for Medical Research, London, UK; , Correspondence: Division of Mathematical Biology, MRC National Institute for Medical Research, The Ridgeway, Mill Hill, NW7 1AA London, UK. E-mail: A.P., ; J.K.,
| |
Collapse
|