1
|
Inan T, Flinko R, Lewis GK, MacKerell AD, Kurkcuoglu O. Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation. J Phys Chem B 2024; 128:5157-5174. [PMID: 38647430 PMCID: PMC11139592 DOI: 10.1021/acs.jpcb.4c00925] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 04/04/2024] [Accepted: 04/08/2024] [Indexed: 04/25/2024]
Abstract
The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.
Collapse
Affiliation(s)
- Tugce Inan
- Department
of Chemical Engineering, Istanbul Technical
University, Istanbul 34469, Turkey
| | - Robin Flinko
- Institute
of Human Virology, University of Maryland
School of Medicine, Baltimore, Maryland 21201, United States
| | - George K. Lewis
- Institute
of Human Virology, University of Maryland
School of Medicine, Baltimore, Maryland 21201, United States
| | - 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
|
2
|
Haliloglu T, Hacisuleyman A, Erman B. Prediction of Allosteric Communication Pathways in Proteins. Bioinformatics 2022; 38:3590-3599. [PMID: 35674396 DOI: 10.1093/bioinformatics/btac380] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 04/12/2022] [Accepted: 06/01/2022] [Indexed: 12/24/2022] Open
Abstract
MOTIVATION Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form. RESULTS Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Turkan Haliloglu
- Polymer Research Center and Chemical Engineering Department, Bogazici University, 34342, Turkey
| | - Aysima Hacisuleyman
- Institute of Bioengineering, Swiss Federal Institute of Technology (EPFL), 1015, Switzerland
| | - Burak Erman
- Chemical and Biological Engineering, Koc University, 34450, Turkey
| |
Collapse
|
3
|
Kutlu Y, Ben-Tal N, Haliloglu T. Global Dynamics Renders Protein Sites with High Functional Response. J Phys Chem B 2021; 125:4734-4745. [PMID: 33914546 DOI: 10.1021/acs.jpcb.1c02511] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Deep mutational scanning enables examination of the effects of many mutations at each amino acid position in a query protein, readily disclosing positions that are particularly sensitive. Mutations in these positions alter protein function the most. Here, on the premise that dynamics underlie function, we explore to what extent the measured sensitivity to mutations could be linked to-perhaps be explained by-the structural dynamics of the protein. We employ a minimalist perturbation-response approach based on the Gaussian Network Model (GNM) on a data set of seven proteins with deep mutational scanning data. The analysis shows that the mutation-sensitive positions are often of capacity to modulate the global dynamics and to intermediate allosteric interactions in the structure. With that, upon strain perturbation, these positions decrease residue fluctuations the most, affecting function via entropy changes. This is particularly relevant for positions that are distant from binding sites or other functional regions of the protein and are sensitive to mutations, nevertheless. Our results indicate that mutations in these positions allosterically manipulate protein function.
Collapse
Affiliation(s)
- Yiǧit Kutlu
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Bebek, Istanbul 34342, Turkey
| | - Nir Ben-Tal
- Department of Biochemistry and Molecular Biology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Turkan Haliloglu
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Bebek, Istanbul 34342, Turkey
| |
Collapse
|
4
|
Sora V, Sanchez D, Papaleo E. Bcl-xL Dynamics under the Lens of Protein Structure Networks. J Phys Chem B 2021; 125:4308-4320. [PMID: 33848145 DOI: 10.1021/acs.jpcb.0c11562] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Understanding the finely orchestrated interactions leading to or preventing programmed cell death (apoptosis) is of utmost importance in cancer research because the failure of these systems could eventually lead to the onset of the disease. In this regard, the maintenance of a delicate balance between the promoters and inhibitors of mitochondrial apoptosis is crucial, as demonstrated by the interplay among the Bcl-2 family members. In particular, B-cell lymphoma extra-large (Bcl-xL) is a target of interest due to the forefront role of its dysfunctions in cancer development. Bcl-xL prevents apoptosis by binding both the pro-apoptotic BH3-only proteins, like PUMA, and the noncanonical partners, such as p53, at different sites. An allosteric communication between the BH3-only protein binding pocket and the p53 binding site, mediating the release of p53 from Bcl-xL upon PUMA binding, has been postulated and supported by nuclear magnetic resonance and other biophysical data. The molecular details of this mechanism, especially at the residue level, remain unclear. In this work, we investigated the distal communication between these two sites in Bcl-xL in its free state and when bound to PUMA. We also evaluated how missense mutations of Bcl-xL found in cancer samples might impair this communication and therefore the allosteric mechanism. We employed all-atom explicit solvent microsecond molecular dynamics simulations, analyzed through a Protein Structure Network approach and integrated with calculations of changes in free energies upon cancer-related mutations identified by genomics studies. We found a subset of candidate residues responsible for both maintaining protein stability and for conveying structural information between the two binding sites and hypothesized possible communication routes between specific residues at both sites.
Collapse
Affiliation(s)
- Valentina Sora
- Computational Biology Laboratory, Danish Cancer Society Research Center, 2100 Copenhagen, Denmark.,Cancer Systems Biology, Section for Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Dionisio Sanchez
- Computational Biology Laboratory, Danish Cancer Society Research Center, 2100 Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, 2100 Copenhagen, Denmark.,Cancer Systems Biology, Section for Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| |
Collapse
|
5
|
Acuner SE, Sumbul F, Torun H, Haliloglu T. Oncogenic mutations on Rac1 affect global intrinsic dynamics underlying GTP and PAK1 binding. Biophys J 2021; 120:866-876. [PMID: 33515600 PMCID: PMC8008323 DOI: 10.1016/j.bpj.2021.01.016] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Revised: 12/31/2020] [Accepted: 01/08/2021] [Indexed: 12/31/2022] Open
Abstract
Rac1 is a small member of the Rho GTPase family. One of the most important downstream effectors of Rac1 is a serine/threonine kinase, p21-activated kinase 1 (PAK1). Mutational activation of PAK1 by Rac1 has oncogenic signaling effects. Here, although we focus on Rac1-PAK1 interaction by atomic-force-microscopy-based single-molecule force spectroscopy experiments, we explore the effect of active mutations on the intrinsic dynamics and binding interactions of Rac1 by Gaussian network model analysis and molecular dynamics simulations. We observe that Rac1 oncogenic mutations are at the hinges of three global modes of motion, suggesting the mechanical changes as potential markers of oncogenicity. Indeed, the dissociation of wild-type Rac1-PAK1 complex shows two distinct unbinding dynamic states that are reduced to one with constitutively active Q61L and oncogenic Y72C mutant Rac1, as revealed by single-molecule force spectroscopy experiments. Q61L and Y72C mutations change the mechanics of the Rac1-PAK1 complex by increasing the elasticity of the protein and slowing down the transition to the unbound state. On the other hand, Rac1's intrinsic dynamics reveal more flexible GTP and PAK1-binding residues on switches I and II with Q61L, Y72C, oncogenic P29S and Q61R, and negative T17N mutations. The cooperativity in the fluctuations of GTP-binding sites around the p-loop and switch I decreases in all mutants, mostly in Q61L, whereas some PAK1-binding residues display enhanced coupling with GTP-binding sites in Q61L and Y72C and within each other in P29S. The predicted binding free energies of the modeled Rac1-PAK1 complexes show that the change in the dynamic behavior likely means a more favorable PAK1 interaction. Overall, these findings suggest that the active mutations affect intrinsic functional dynamic events and alter the mechanics underlying the binding of Rac1 to GTP and upstream and downstream partners including PAK1.
Collapse
Affiliation(s)
- Saliha Ece Acuner
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey
| | - Fidan Sumbul
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey
| | - Hamdi Torun
- Department of Mathematics, Physics and Electrical Engineering, Northumbria University, Newcastle, United Kingdom.
| | - Turkan Haliloglu
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey.
| |
Collapse
|
6
|
Sayılgan JF, Haliloğlu T, Gönen M. Protein dynamics analysis identifies candidate cancer driver genes and mutations in TCGA data. Proteins 2021; 89:721-730. [PMID: 33550612 DOI: 10.1002/prot.26054] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2020] [Revised: 01/04/2021] [Accepted: 01/31/2021] [Indexed: 11/09/2022]
Abstract
Recently, it has been showed that cancer missense mutations selectively target the neighborhood of hinge residues, which are key sites in protein dynamics. Here, we show that this approach can be extended to find previously unknown candidate mutations and genes. To this aim, we developed a computational pipeline to detect significantly enriched three-dimensional (3D) clustering of missense mutations around hinge residues. The hinge residues were detected by applying a Gaussian network model. By systematically analyzing the PanCancer compendium of somatic mutations in nearly 10 000 tumors from the Cancer Genome Atlas, we identified candidate genes and mutations in addition to well known ones. For instance, we found significantly enriched 3D clustering of missense mutations in known cancer genes including CDK4, CDKN2A, TCL1A, and MAPK1. Beside these known genes, we also identified significantly enriched 3D clustering of missense mutations around hinge residues in PLA2G4A, which may lead to excessive phosphorylation of the extracellular signal-regulated kinases. Furthermore, we demonstrated that hinge-based features improves pathogenicity prediction for missense mutations. Our results show that the consideration of clustering around hinge residues can help us explain the functional role of the mutations in known cancer genes and identify candidate genes.
Collapse
Affiliation(s)
- Jan Fehmi Sayılgan
- Graduate School of Sciences and Engineering, Koç University, İstanbul, Turkey
| | - Türkan Haliloğlu
- Department of Chemical Engineering, School of Engineering, Boğaziçi University, İstanbul, Turkey.,Polymer Research Center, Boğaziçi University, İstanbul, Turkey
| | - Mehmet Gönen
- Department of Industrial Engineering, College of Engineering, Koç University, İstanbul, Turkey.,School of Medicine, Koç University, İstanbul, Turkey
| |
Collapse
|
7
|
Abstract
Allostery is a fundamental regulatory mechanism in the majority of biological processes of molecular machines. Allostery is well-known as a dynamic-driven process, and thus, the molecular mechanism of allosteric signal transmission needs to be established. Elastic network models (ENMs) provide efficient methods for investigating the intrinsic dynamics and allosteric communication pathways in proteins. In this chapter, two ENM methods including Gaussian network model (GNM) coupled with Markovian stochastic model, as well as the anisotropic network model (ANM), were introduced to identify allosteric effects in hemoglobins. Techniques on model parameters, scripting and calculation, analysis, and visualization are shown step by step.
Collapse
|
8
|
Yan W, Liu X, Wang Y, Han S, Wang F, Liu X, Xiao F, Hu G. Identifying Drug Targets in Pancreatic Ductal Adenocarcinoma Through Machine Learning, Analyzing Biomolecular Networks, and Structural Modeling. Front Pharmacol 2020; 11:534. [PMID: 32425783 PMCID: PMC7204992 DOI: 10.3389/fphar.2020.00534] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Accepted: 04/06/2020] [Indexed: 12/16/2022] Open
Abstract
Pancreatic ductal adenocarcinoma (PDAC) is one of the leading causes of cancer-related death and has an extremely poor prognosis. Thus, identifying new disease-associated genes and targets for PDAC diagnosis and therapy is urgently needed. This requires investigations into the underlying molecular mechanisms of PDAC at both the systems and molecular levels. Herein, we developed a computational method of predicting cancer genes and anticancer drug targets that combined three independent expression microarray datasets of PDAC patients and protein-protein interaction data. First, Support Vector Machine–Recursive Feature Elimination was applied to the gene expression data to rank the differentially expressed genes (DEGs) between PDAC patients and controls. Then, protein-protein interaction networks were constructed based on the DEGs, and a new score comprising gene expression and network topological information was proposed to identify cancer genes. Finally, these genes were validated by “druggability” prediction, survival and common network analysis, and functional enrichment analysis. Furthermore, two integrins were screened to investigate their structures and dynamics as potential drug targets for PDAC. Collectively, 17 disease genes and some stroma-related pathways including extracellular matrix-receptor interactions were predicted to be potential drug targets and important pathways for treating PDAC. The protein-drug interactions and hinge sites predication of ITGAV and ITGA2 suggest potential drug binding residues in the Thigh domain. These findings provide new possibilities for targeted therapeutic interventions in PDAC, which may have further applications in other cancer types.
Collapse
Affiliation(s)
- Wenying Yan
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Xingyi Liu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Yibo Wang
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Shuqing Han
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Fan Wang
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Xin Liu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Fei Xiao
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China.,State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou, China
| |
Collapse
|
9
|
Tang QY, Kaneko K. Long-range correlation in protein dynamics: Confirmation by structural data and normal mode analysis. PLoS Comput Biol 2020; 16:e1007670. [PMID: 32053592 PMCID: PMC7043781 DOI: 10.1371/journal.pcbi.1007670] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 02/26/2020] [Accepted: 01/21/2020] [Indexed: 11/18/2022] Open
Abstract
Proteins in cellular environments are highly susceptible. Local perturbations to any residue can be sensed by other spatially distal residues in the protein molecule, showing long-range correlations in the native dynamics of proteins. The long-range correlations of proteins contribute to many biological processes such as allostery, catalysis, and transportation. Revealing the structural origin of such long-range correlations is of great significance in understanding the design principle of biologically functional proteins. In this work, based on a large set of globular proteins determined by X-ray crystallography, by conducting normal mode analysis with the elastic network models, we demonstrate that such long-range correlations are encoded in the native topology of the proteins. To understand how native topology defines the structure and the dynamics of the proteins, we conduct scaling analysis on the size dependence of the slowest vibration mode, average path length, and modularity. Our results quantitatively describe how native proteins balance between order and disorder, showing both dense packing and fractal topology. It is suggested that the balance between stability and flexibility acts as an evolutionary constraint for proteins at different sizes. Overall, our result not only gives a new perspective bridging the protein structure and its dynamics but also reveals a universal principle in the evolution of proteins at all different sizes. The long-range correlated fluctuations are closely related to many biological processes of the proteins, such as catalysis, ligand binding, biomolecular recognition, and transportation. In this paper, we elucidate the structural origin of the long-range correlation and describe how native contact topology defines the slow-mode dynamics of the native proteins. Our result suggests an evolutionary constraint for proteins at different sizes, which may shed light on solving many biophysical problems such as structure prediction, multi-scale molecular simulations, and the design of molecular machines. Moreover, in statistical physics, as the long-range correlations are notable signs of the critical point, unveiling the origin of such criticality can extend our understanding of the organizing principle of a large variety of complex systems.
Collapse
Affiliation(s)
- Qian-Yuan Tang
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Tokyo, Japan
- * E-mail:
| | - Kunihiko Kaneko
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Tokyo, Japan
| |
Collapse
|
10
|
Jankauskaite J, Jiménez-García B, Dapkunas J, Fernández-Recio J, Moal IH. SKEMPI 2.0: an updated benchmark of changes in protein-protein binding energy, kinetics and thermodynamics upon mutation. Bioinformatics 2019; 35:462-469. [PMID: 30020414 PMCID: PMC6361233 DOI: 10.1093/bioinformatics/bty635] [Citation(s) in RCA: 141] [Impact Index Per Article: 28.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Accepted: 07/17/2018] [Indexed: 11/18/2022] Open
Abstract
Motivation Understanding the relationship between the sequence, structure, binding energy, binding kinetics and binding thermodynamics of protein–protein interactions is crucial to understanding cellular signaling, the assembly and regulation of molecular complexes, the mechanisms through which mutations lead to disease, and protein engineering. Results We present SKEMPI 2.0, a major update to our database of binding free energy changes upon mutation for structurally resolved protein–protein interactions. This version now contains manually curated binding data for 7085 mutations, an increase of 133%, including changes in kinetics for 1844 mutations, enthalpy and entropy changes for 443 mutations, and 440 mutations, which abolish detectable binding. Availability and implementation The database is available as supplementary data and at https://life.bsc.es/pid/skempi2/. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Justina Jankauskaite
- Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Brian Jiménez-García
- Barcelona Supercomputing Center (BSC), Barcelona, Spain.,Bijvoet Center for Biomolecular Research, Faculty of Science, Utrecht University, Utrecht, the Netherlands
| | - Justas Dapkunas
- Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Juan Fernández-Recio
- Barcelona Supercomputing Center (BSC), Barcelona, Spain.,Institut de Biologia Molecular de Barcelona (IBMB), CSIC, Barcelona, Spain
| | - Iain H Moal
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Hinxton, Cambridge, UK
| |
Collapse
|
11
|
Sayılgan JF, Haliloğlu T, Gönen M. Protein dynamics analysis reveals that missense mutations in cancer‐related genes appear frequently on hinge‐neighboring residues. Proteins 2019; 87:512-519. [DOI: 10.1002/prot.25673] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Revised: 01/09/2019] [Accepted: 02/17/2019] [Indexed: 01/26/2023]
Affiliation(s)
- Jan Fehmi Sayılgan
- Graduate School of Sciences and EngineeringKoç University İstanbul Turkey
| | - Türkan Haliloğlu
- Department of Chemical Engineering, School of EngineeringBoğaziçi University İstanbul Turkey
- Polymer Research CenterBoğaziçi University İstanbul Turkey
| | - Mehmet Gönen
- Department of Industrial Engineering, College of EngineeringKoç University İstanbul Turkey
- School of MedicineKoç University İstanbul Turkey
- Department of Biomedical Engineering, School of MedicineOregon Health and Science University Portland Oregon
| |
Collapse
|
12
|
Liang Z, Hu J, Yan W, Jiang H, Hu G, Luo C. Deciphering the role of dimer interface in intrinsic dynamics and allosteric pathways underlying the functional transformation of DNMT3A. Biochim Biophys Acta Gen Subj 2018; 1862:1667-1679. [DOI: 10.1016/j.bbagen.2018.04.015] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2017] [Revised: 04/03/2018] [Accepted: 04/13/2018] [Indexed: 01/06/2023]
|
13
|
Abstract
An orthosteric site is commonly viewed as the primary, functionally binding pocket on a receptor. Signal molecules, endogenous agonists, and substrates are recognized by and bind to the orthosteric site of a specific target, resulting in a biological effect. A malfunctioning active site on a crucial receptor has been confirmed as the culprit that causes many metabolic disturbances, neurologic disorders, and genetic diseases. A competitive inhibitor that has a stronger binding affinity can outcompete an orthosteric ligand. An allosteric site, which is nonoverlapping and topographically distinct from the active pocket, can emerge as a potential regulatory site on the protein surface. An allosteric modulator interacts with a specific binding site, affecting the atoms of nearby residues, thus eliciting a series of conformational changes in the residues at the active site through propagation pathways. Allosteric regulation can potentiate or inhibit function instead of blocking it, and this is a promising strategy for drug design. In this chapter, we describe the tools and protocols for allosteric site analysis and allosteric ligand design.
Collapse
Affiliation(s)
- Kun Song
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Ministry of Education, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Jian Zhang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Ministry of Education, Shanghai Jiao-Tong University School of Medicine, Shanghai, China.
| |
Collapse
|
14
|
Salamanca Viloria J, Allega MF, Lambrughi M, Papaleo E. An optimal distance cutoff for contact-based Protein Structure Networks using side-chain centers of mass. Sci Rep 2017; 7:2838. [PMID: 28588190 PMCID: PMC5460117 DOI: 10.1038/s41598-017-01498-6] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Accepted: 03/28/2017] [Indexed: 02/05/2023] Open
Abstract
Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.
Collapse
Affiliation(s)
- Juan Salamanca Viloria
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Maria Francesca Allega
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark.
| |
Collapse
|
15
|
Schueler-Furman O, Wodak SJ. Computational approaches to investigating allostery. Curr Opin Struct Biol 2016; 41:159-171. [PMID: 27607077 DOI: 10.1016/j.sbi.2016.06.017] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Accepted: 06/23/2016] [Indexed: 01/01/2023]
Abstract
Allosteric regulation plays a key role in many biological processes, such as signal transduction, transcriptional regulation, and many more. It is rooted in fundamental thermodynamic and dynamic properties of macromolecular systems that are still poorly understood and are moreover modulated by the cellular context. Here we review the computational approaches used in the investigation of allosteric processes in protein systems. We outline how the models of allostery have evolved from their initial formulation in the sixties to the current views, which more fully account for the roles of the thermodynamic and dynamic properties of the system. We then describe the major classes of computational approaches employed to elucidate the mechanisms of allostery, the insights they have provided, as well as their limitations. We complement this analysis by highlighting the role of computational approaches in promising practical applications, such as the engineering of regulatory modules and identifying allosteric binding sites.
Collapse
Affiliation(s)
- Ora Schueler-Furman
- Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada (IMRIC), Hebrew University, Hadassah Medical School, POB 12272, Jerusalem 91120, Israel
| | - Shoshana J Wodak
- VIB Structural Biology Research Center, VUB, Pleinlaan 2, 1050 Brussels, Belgium.
| |
Collapse
|
16
|
Wei G, Xi W, Nussinov R, Ma B. Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell. Chem Rev 2016; 116:6516-51. [PMID: 26807783 PMCID: PMC6407618 DOI: 10.1021/acs.chemrev.5b00562] [Citation(s) in RCA: 253] [Impact Index Per Article: 31.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.
Collapse
Affiliation(s)
- Guanghong Wei
- State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (MOE), and Department of Physics, Fudan University, Shanghai, P. R. China
| | - Wenhui Xi
- State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (MOE), and Department of Physics, Fudan University, Shanghai, P. R. China
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland 21702, USA
- Sackler Inst. of Molecular Medicine Department of Human Genetics and Molecular Medicine Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland 21702, USA
| |
Collapse
|
17
|
Abstract
Allosteric transition, defined as conformational changes induced by ligand binding, is one of the fundamental properties of proteins. Allostery has been observed and characterized in many proteins, and has been recently utilized to control protein function via regulation of protein activity. Here, we review the physical and evolutionary origin of protein allostery, as well as its importance to protein regulation, drug discovery, and biological processes in living systems. We describe recently developed approaches to identify allosteric pathways, connected sets of pairwise interactions that are responsible for propagation of conformational change from the ligand-binding site to a distal functional site. We then present experimental and computational protein engineering approaches for control of protein function by modulation of allosteric sites. As an example of application of these approaches, we describe a synergistic computational and experimental approach to rescue the cystic-fibrosis-associated protein cystic fibrosis transmembrane conductance regulator, which upon deletion of a single residue misfolds and causes disease. This example demonstrates the power of allosteric manipulation in proteins to both elucidate mechanisms of molecular function and to develop therapeutic strategies that rescue those functions. Allosteric control of proteins provides a tool to shine a light on the complex cascades of cellular processes and facilitate unprecedented interrogation of biological systems.
Collapse
Affiliation(s)
- Nikolay V Dokholyan
- Department of Biochemistry and Biophysics, University of North Carolina , Chapel Hill, North Carolina 27599, United States
| |
Collapse
|