1
|
Diessner EM, Freites JA, Tobias DJ, Butts CT. Network Hamiltonian Models for Unstructured Protein Aggregates, with Application to γD-Crystallin. J Phys Chem B 2023; 127:685-697. [PMID: 36637342 PMCID: PMC10437096 DOI: 10.1021/acs.jpcb.2c07672] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Network Hamiltonian models (NHMs) are a framework for topological coarse-graining of protein-protein interactions, in which each node corresponds to a protein, and edges are drawn between nodes representing proteins that are noncovalently bound. Here, this framework is applied to aggregates of γD-crystallin, a structural protein of the eye lens implicated in cataract disease. The NHMs in this study are generated from atomistic simulations of equilibrium distributions of wild-type and the cataract-causing variant W42R in solution, performed by Wong, E. K.; Prytkova, V.; Freites, J. A.; Butts, C. T.; Tobias, D. J. Molecular Mechanism of Aggregation of the Cataract-Related γD-Crystallin W42R Variant from Multiscale Atomistic Simulations. Biochemistry2019, 58 (35), 3691-3699. Network models are shown to successfully reproduce the aggregate size and structure observed in the atomistic simulation, and provide information about the transient protein-protein interactions therein. The system size is scaled from the original 375 monomers to a system of 10000 monomers, revealing a lowering of the upper tail of the aggregate size distribution of the W42R variant. Extrapolation to higher and lower concentrations is also performed. These results provide an example of the utility of NHMs for coarse-grained simulation of protein systems, as well as their ability to scale to large system sizes and high concentrations, reducing computational costs while retaining topological information about the system.
Collapse
Affiliation(s)
- Elizabeth M Diessner
- Department of Chemistry, University of California, Irvine, California92697, United States
| | - J Alfredo Freites
- Department of Chemistry, University of California, Irvine, California92697, United States
| | - Douglas J Tobias
- Department of Chemistry, University of California, Irvine, California92697, United States
| | - Carter T Butts
- Departments of Sociology, Statistics, Computer Science, and EECS, University of California, Irvine, California92697, United States
| |
Collapse
|
2
|
Chaudhuri S, Srivastava A. Network approach to understand biological systems: From single to multilayer networks. J Biosci 2022. [DOI: 10.1007/s12038-022-00285-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
3
|
Fanelli F, Felline A. Uncovering GPCR and G Protein Function by Protein Structure Network Analysis. COMPUTATIONAL TOOLS FOR CHEMICAL BIOLOGY 2017. [DOI: 10.1039/9781788010139-00198] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used for investigating structural communication in biomolecular systems. Information on the system's dynamics can be provided by atomistic molecular dynamics (MD) simulations or coarse grained elastic network models paired with normal mode analysis (ENM-NMA). This chapter reports on selected applications of PSN analysis to uncover the structural communication in G protein coupled receptors (GPCRs) and G proteins. Strategies to highlight changes in structural communication caused by mutations, ligand and protein binding are described. Conserved amino acids, sites of misfolding mutations, or ligands acting as functional switches tend to behave as hubs in the native structure networks. Densely linked regions in the protein structure graphs could be identified as playing central roles in protein stability and function. Changes in the communication pathway fingerprints depending on the bound ligand or following amino acid mutation could be highlighted as well. A bridge between misfolding and misrouting could be established in rhodopsin mutants linked to inherited blindness. The analysis of native network perturbations by misfolding mutations served to infer key structural elements of protein responsiveness to small chaperones with implications for drug discovery.
Collapse
Affiliation(s)
- Francesca Fanelli
- Department of Life Sciences University of Modena and Reggio Emilia Italy
- Center for Neuroscience and Neurotechnology University of Modena and Reggio Emilia Italy
| | - Angelo Felline
- Department of Life Sciences University of Modena and Reggio Emilia Italy
| |
Collapse
|
4
|
Srivastava A, Sinha S. Uncoupling of an ammonia channel as a mechanism of allosteric inhibition in anthranilate synthase of Serratia marcescens: dynamic and graph theoretical analysis. MOLECULAR BIOSYSTEMS 2017; 13:142-155. [DOI: 10.1039/c6mb00646a] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Network modeling and molecular dynamic studies reveal the perturbation in communication pathways as a mechanism of allosteric inhibition in anthranilate synthase.
Collapse
Affiliation(s)
- Ashutosh Srivastava
- Centre for Protein Science
- Design
- Engineering (CPSDE)
- Department of Biological Sciences
- Indian Institute of Science Education Research Mohali
| | - Somdatta Sinha
- Centre for Protein Science
- Design
- Engineering (CPSDE)
- Department of Biological Sciences
- Indian Institute of Science Education Research Mohali
| |
Collapse
|
5
|
Tasdighian S, Di Paola L, De Ruvo M, Paci P, Santoni D, Palumbo P, Mei G, Di Venere A, Giuliani A. Modules Identification in Protein Structures: The Topological and Geometrical Solutions. J Chem Inf Model 2013; 54:159-68. [DOI: 10.1021/ci400218v] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Setareh Tasdighian
- Department
of Plant Systems Biology, Flanders Interuniversity Institute for Biotechnology, Ghent University, K.L. Ledeganckstraat 35, B-9000 Ghent, Belgium
| | - Luisa Di Paola
- Faculty
of Engineering, Università CAMPUS BioMedico, Via A. del
Portillo, 21, 00128 Roma, Italy
| | - Micol De Ruvo
- CNR-Institute of Systems Analysis and Computer Science (IASI), viale Manzoni 30, 00185 Roma, Italy
| | - Paola Paci
- CNR-Institute of Systems Analysis and Computer Science (IASI), viale Manzoni 30, 00185 Roma, Italy
| | - Daniele Santoni
- Department
of Experimental Medicine and Surgery, University of Rome “Tor Vergata”, via Montpellier 1, 00133 Rome, Italy
| | - Pasquale Palumbo
- CNR-Institute of Systems Analysis and Computer Science (IASI), viale Manzoni 30, 00185 Roma, Italy
- Department
of Experimental Medicine and Surgery, University of Rome “Tor Vergata”, via Montpellier 1, 00133 Rome, Italy
| | - Giampiero Mei
- Department
of Experimental Medicine and Surgery, University of Rome “Tor Vergata”, via Montpellier 1, 00133 Rome, Italy
| | - Almerinda Di Venere
- Environment
and Health Department, Istituto Superiore di Sanità, Viale
Regina Elena 299, 00161, Roma, Italy
| | - Alessandro Giuliani
- Environment
and Health Department, Istituto Superiore di Sanità, Viale
Regina Elena 299, 00161, Roma, Italy
| |
Collapse
|
6
|
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
|
7
|
Di Paola L, De Ruvo M, Paci P, Santoni D, Giuliani A. Protein Contact Networks: An Emerging Paradigm in Chemistry. Chem Rev 2012. [DOI: 10.1021/cr3002356] [Citation(s) in RCA: 173] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- L. Di Paola
- Faculty of Engineering, Università CAMPUS BioMedico, Via A. del Portillo,
21, 00128 Roma, Italy
| | | | | | - D. Santoni
- BioMathLab, CNR-Institute of Systems Analysis and Computer Science (IASI), viale Manzoni 30, 00185
Roma, Italy
| | - A. Giuliani
- Environment
and Health Department, Istituto Superiore di Sanità, Viale Regina Elena
299, 00161, Roma, Italy
| |
Collapse
|
8
|
Turgut D, Atilgan AR, Atilgan C. Assortative mixing in close-packed spatial networks. PLoS One 2010; 5:e15551. [PMID: 21179578 PMCID: PMC3002975 DOI: 10.1371/journal.pone.0015551] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2010] [Accepted: 11/11/2010] [Indexed: 12/02/2022] Open
Abstract
BACKGROUND In recent years, there is aroused interest in expressing complex systems as networks of interacting nodes. Using descriptors from graph theory, it has been possible to classify many diverse systems derived from social and physical sciences alike. In particular, folded proteins as examples of self-assembled complex molecules have also been investigated intensely using these tools. However, we need to develop additional measures to classify different systems, in order to dissect the underlying hierarchy. METHODOLOGY AND PRINCIPAL FINDINGS In this study, a general analytical relation for the dependence of nearest neighbor degree correlations on degree is derived. Dependence of local clustering on degree is shown to be the sole determining factor of assortative versus disassortative mixing in networks. The characteristics of networks constructed from spatial atomic/molecular systems exemplified by self-organized residue networks built from folded protein structures and block copolymers, atomic clusters and well-compressed polymeric melts are studied. Distributions of statistical properties of the networks are presented. For these densely-packed systems, assortative mixing in the network construction is found to apply, and conditions are derived for a simple linear dependence. CONCLUSIONS Our analyses (i) reveal patterns that are common to close-packed clusters of atoms/molecules, (ii) identify the type of surface effects prominent in different close-packed systems, and (iii) associate fingerprints that may be used to classify networks with varying types of correlations.
Collapse
Affiliation(s)
- Deniz Turgut
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Ali Rana Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| |
Collapse
|
9
|
Ghosh A, Vishveshwara S. Variations in clique and community patterns in protein structures during allosteric communication: investigation of dynamically equilibrated structures of methionyl tRNA synthetase complexes. Biochemistry 2008; 47:11398-407. [PMID: 18842003 DOI: 10.1021/bi8007559] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The allosteric concept has played a key role in understanding the biological functions of proteins. The rigidity or plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. Although molecular insights have been gained from a large number of structures, a precise assessment of the ligand-induced conformational changes in proteins at different levels, ranging from gross topology to intricate details, remains a challenge. In this study, we have explored the conformational changes in the complexes of methionyl tRNA synthetase (MetRS) through novel network parameters such as cliques and communities, which identify the rigid regions in the protein structure networks (PSNs) constructed from the noncovalent interactions of amino acid side chains. MetRS belongs to the aminoacyl tRNA synthetase (aaRS) family that plays a crucial role in the translation of genetic code. These enzymes are modular with distinct domains from which extensive genetic, kinetic, and structural data are available, highlighting the role of interdomain communication. The network parameters evaluated here on the conformational ensembles of MetRS complexes, generated from molecular dynamics simulations, have enabled us to understand the interdomain communication in detail. Additionally, the characterization of conformational changes in terms of cliques and communities has also become possible, which had eluded conventional analyses. Furthermore, we find that most of the residues participating in cliques and communities are strikingly different from those that take part in long-range communication. The cliques and communities evaluated here for the first time on PSNs have beautifully captured the local geometries in detail within the framework of global topology. Here the allosteric effect is revealed at the residue level via identification of the important residues specific for structural rigidity and functional flexibility in MetRS. This ought to enhance our understanding of the functioning of aaRS in general.
Collapse
Affiliation(s)
- Amit Ghosh
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India 560012
| | | |
Collapse
|
10
|
Goodey NM, Benkovic SJ. Allosteric regulation and catalysis emerge via a common route. Nat Chem Biol 2008; 4:474-82. [PMID: 18641628 DOI: 10.1038/nchembio.98] [Citation(s) in RCA: 531] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Allosteric regulation of protein function is a mechanism by which an event in one place of a protein structure causes an effect at another site, much like the behavior of a telecommunications network in which a collection of transmitters, receivers and transceivers communicate with each other across long distances. For example, ligand binding or an amino acid mutation at an allosteric site can alter enzymatic activity or binding affinity in a distal region such as the active site or a second binding site. The mechanism of this site-to-site communication is of great interest, especially since allosteric effects must be considered in drug design and protein engineering. In this review, conformational mobility as the common route between allosteric regulation and catalysis is discussed. We summarize recent experimental data and the resulting insights into allostery within proteins, and we discuss the nature of future studies and the new applications that may result from increased understanding of this regulatory mechanism.
Collapse
Affiliation(s)
- Nina M Goodey
- Montclair State University, Department of Chemistry and Biochemistry, 1 Normal Avenue, Montclair, New Jersey 07043, USA
| | | |
Collapse
|
11
|
Tsai CJ, Sol AD, Nussinov R. Allostery: absence of a change in shape does not imply that allostery is not at play. J Mol Biol 2008; 378:1-11. [PMID: 18353365 PMCID: PMC2684958 DOI: 10.1016/j.jmb.2008.02.034] [Citation(s) in RCA: 361] [Impact Index Per Article: 22.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2007] [Revised: 02/15/2008] [Accepted: 02/15/2008] [Indexed: 11/17/2022]
Abstract
Allostery is essential for controlled catalysis, signal transmission, receptor trafficking, turning genes on and off, and apoptosis. It governs the organism's response to environmental and metabolic cues, dictating transient partner interactions in the cellular network. Textbooks taught us that allostery is a change of shape at one site on the protein surface brought about by ligand binding to another. For several years, it has been broadly accepted that the change of shape is not induced; rather, it is observed simply because a larger protein population presents it. Current data indicate that while side chains can reorient and rewire, allostery may not even involve a change of (backbone) shape. Assuming that the enthalpy change does not reverse the free-energy change due to the change in entropy, entropy is mainly responsible for binding.
Collapse
Affiliation(s)
- Chung-Jung Tsai
- Basic Research Program, SAIC-Frederick, Inc., Center for Cancer Research Nanobiology Program, NCI-Frederick, Frederick, MD 21702
| | - Antonio del Sol
- Bioinformatics Research Unit, Research and Development Division, Fujirebio Inc., 51 Komiya-cho, Hachioji-shi, Tokyo 192-0031, Japan
| | - Ruth Nussinov
- Basic Research Program, SAIC-Frederick, Inc., Center for Cancer Research Nanobiology Program, NCI-Frederick, Frederick, MD 21702
- Sackler Inst. of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
12
|
González-Díaz H, González-Díaz Y, Santana L, Ubeira FM, Uriarte E. Proteomics, networks and connectivity indices. Proteomics 2008; 8:750-78. [DOI: 10.1002/pmic.200700638] [Citation(s) in RCA: 170] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
13
|
Ghosh A, Vishveshwara S. A study of communication pathways in methionyl- tRNA synthetase by molecular dynamics simulations and structure network analysis. Proc Natl Acad Sci U S A 2007; 104:15711-6. [PMID: 17898174 PMCID: PMC2000407 DOI: 10.1073/pnas.0704459104] [Citation(s) in RCA: 186] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by approximately 70 A in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.
Collapse
Affiliation(s)
- Amit Ghosh
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India
| | - Saraswathi Vishveshwara
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India
- *To whom correspondence may be addressed. E-mail:
| |
Collapse
|