1
|
Moldovean-Cioroianu NS. Reviewing the Structure-Function Paradigm in Polyglutamine Disorders: A Synergistic Perspective on Theoretical and Experimental Approaches. Int J Mol Sci 2024; 25:6789. [PMID: 38928495 PMCID: PMC11204371 DOI: 10.3390/ijms25126789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2024] [Revised: 06/13/2024] [Accepted: 06/17/2024] [Indexed: 06/28/2024] Open
Abstract
Polyglutamine (polyQ) disorders are a group of neurodegenerative diseases characterized by the excessive expansion of CAG (cytosine, adenine, guanine) repeats within host proteins. The quest to unravel the complex diseases mechanism has led researchers to adopt both theoretical and experimental methods, each offering unique insights into the underlying pathogenesis. This review emphasizes the significance of combining multiple approaches in the study of polyQ disorders, focusing on the structure-function correlations and the relevance of polyQ-related protein dynamics in neurodegeneration. By integrating computational/theoretical predictions with experimental observations, one can establish robust structure-function correlations, aiding in the identification of key molecular targets for therapeutic interventions. PolyQ proteins' dynamics, influenced by their length and interactions with other molecular partners, play a pivotal role in the polyQ-related pathogenic cascade. Moreover, conformational dynamics of polyQ proteins can trigger aggregation, leading to toxic assembles that hinder proper cellular homeostasis. Understanding these intricacies offers new avenues for therapeutic strategies by fine-tuning polyQ kinetics, in order to prevent and control disease progression. Last but not least, this review highlights the importance of integrating multidisciplinary efforts to advancing research in this field, bringing us closer to the ultimate goal of finding effective treatments against polyQ disorders.
Collapse
Affiliation(s)
- Nastasia Sanda Moldovean-Cioroianu
- Institute of Materials Science, Bioinspired Materials and Biosensor Technologies, Kiel University, Kaiserstraße 2, 24143 Kiel, Germany;
- Faculty of Physics, Babeș-Bolyai University, Kogălniceanu 1, RO-400084 Cluj-Napoca, Romania
| |
Collapse
|
2
|
Gao CY, Yang GY, Ding XW, Xu JH, Cheng X, Zheng GW, Chen Q. Engineering of Halide Methyltransferase BxHMT through Dynamic Cross-Correlation Network Analysis. Angew Chem Int Ed Engl 2024; 63:e202401235. [PMID: 38623716 DOI: 10.1002/anie.202401235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2024] [Revised: 03/18/2024] [Accepted: 04/15/2024] [Indexed: 04/17/2024]
Abstract
Halide methyltransferases (HMTs) provide an effective way to regenerate S-adenosyl methionine (SAM) from S-adenosyl homocysteine and reactive electrophiles, such as methyl iodide (MeI) and methyl toluene sulfonate (MeOTs). As compared with MeI, the cost-effective unnatural substrate MeOTs can be accessed directly from cheap and abundant alcohols, but shows only limited reactivity in SAM production. In this study, we developed a dynamic cross-correlation network analysis (DCCNA) strategy for quickly identifying hot spots influencing the catalytic efficiency of the enzyme, and applied it to the evolution of HMT from Paraburkholderia xenovorans. Finally, the optimal mutant, M4 (V55T/C125S/L127T/L129P), exhibited remarkable improvement, with a specific activity of 4.08 U/mg towards MeOTs, representing an 82-fold increase as compared to the wild-type (WT) enzyme. Notably, M4 also demonstrated a positive impact on the catalytic ability with other methyl donors. The structural mechanism behind the enhanced enzyme activity was uncovered by molecular dynamics simulations. Our work not only contributes a promising biocatalyst for the regeneration of SAM, but also offers a strategy for efficient enzyme engineering.
Collapse
Affiliation(s)
- Chun-Yu Gao
- State Key Laboratory of Bioreactor Engineering and Shanghai Collaborative Innovation Center for Biomanufacturing, East China University of Science and Technology, Shanghai, 200237, China
| | - Gui-Ying Yang
- State Key Laboratory of Bioreactor Engineering and Shanghai Collaborative Innovation Center for Biomanufacturing, East China University of Science and Technology, Shanghai, 200237, China
| | - Xu-Wei Ding
- State Key Laboratory of Bioreactor Engineering and Shanghai Collaborative Innovation Center for Biomanufacturing, East China University of Science and Technology, Shanghai, 200237, China
| | - Jian-He Xu
- State Key Laboratory of Bioreactor Engineering and Shanghai Collaborative Innovation Center for Biomanufacturing, East China University of Science and Technology, Shanghai, 200237, China
| | - Xiaolin Cheng
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy, The Ohio State University, Columbus, OH 43210, United States
| | - Gao-Wei Zheng
- State Key Laboratory of Bioreactor Engineering and Shanghai Collaborative Innovation Center for Biomanufacturing, East China University of Science and Technology, Shanghai, 200237, China
| | - Qi Chen
- State Key Laboratory of Bioreactor Engineering and Shanghai Collaborative Innovation Center for Biomanufacturing, East China University of Science and Technology, Shanghai, 200237, China
| |
Collapse
|
3
|
Gonzales J, Kim I, Hwang W, Cho JH. Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.24.595776. [PMID: 38826371 PMCID: PMC11142230 DOI: 10.1101/2024.05.24.595776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85β. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85β-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.
Collapse
|
4
|
Agbaglo DA, Summers TJ, Cheng Q, DeYonker NJ. The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study. Phys Chem Chem Phys 2024; 26:12467-12482. [PMID: 38618904 PMCID: PMC11090134 DOI: 10.1039/d3cp06100k] [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] [Indexed: 04/16/2024]
Abstract
Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.
Collapse
Affiliation(s)
- Donatus A Agbaglo
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| | - Thomas J Summers
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| | - Qianyi Cheng
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| | - Nathan J DeYonker
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| |
Collapse
|
5
|
Gorantla R, Kubincová A, Weiße AY, Mey ASJS. From Proteins to Ligands: Decoding Deep Learning Methods for Binding Affinity Prediction. J Chem Inf Model 2024; 64:2496-2507. [PMID: 37983381 PMCID: PMC11005465 DOI: 10.1021/acs.jcim.3c01208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Revised: 10/26/2023] [Accepted: 10/27/2023] [Indexed: 11/22/2023]
Abstract
Accurate in silico prediction of protein-ligand binding affinity is important in the early stages of drug discovery. Deep learning-based methods exist but have yet to overtake more conventional methods such as giga-docking largely due to their lack of generalizability. To improve generalizability, we need to understand what these models learn from input protein and ligand data. We systematically investigated a sequence-based deep learning framework to assess the impact of protein and ligand encodings on predicting binding affinities for commonly used kinase data sets. The role of proteins is studied using convolutional neural network-based encodings obtained from sequences and graph neural network-based encodings enriched with structural information from contact maps. Ligand-based encodings are generated from graph-neural networks. We test different ligand perturbations by randomizing node and edge properties. For proteins, we make use of 3 different protein contact generation methods (AlphaFold2, Pconsc4, and ESM-1b) and compare these with a random control. Our investigation shows that protein encodings do not substantially impact the binding predictions, with no statistically significant difference in binding affinity for KIBA in the investigated metrics (concordance index, Pearson's R Spearman's Rank, and RMSE). Significant differences are seen for ligand encodings with random ligands and random ligand node properties, suggesting a much bigger reliance on ligand data for the learning tasks. Using different ways to combine protein and ligand encodings did not show a significant change in performance.
Collapse
Affiliation(s)
- Rohan Gorantla
- School
of Informatics, University of Edinburgh, Edinburgh, EH8 9AB, U.K.
- EaStCHEM
School of Chemistry, University of Edinburgh, Edinburgh, EH9 3FJ, U.K.
| | | | - Andrea Y. Weiße
- School
of Informatics, University of Edinburgh, Edinburgh, EH8 9AB, U.K.
- School
of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FF, U.K.
| | | |
Collapse
|
6
|
Mahapatra S, Jonniya NA, Koirala S, Kar P. Molecular dynamics simulations reveal phosphorylation-induced conformational dynamics of the fibroblast growth factor receptor 1 kinase. J Biomol Struct Dyn 2024; 42:2929-2941. [PMID: 37160693 DOI: 10.1080/07391102.2023.2209189] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Accepted: 04/26/2023] [Indexed: 05/11/2023]
Abstract
The Fibroblast Growth Factor Receptor1 (FGFR1) kinase wields exquisite control on cell fate, proliferation, differentiation, and homeostasis. An imbalance of FGFR1 signaling leads to several pathogeneses of diseases ranging from multiple cancers to allergic and neurodegenerative disorders. In this study, we investigated the phosphorylation-induced conformational dynamics of FGFR1 in apo and ATP-bound states via all-atom molecular dynamics simulations. All simulations were performed for 2 × 2 µs. We have also investigated the energetics of the binding of ATP to FGFR1 using the molecular mechanics Poisson-Boltzmann scheme. Our study reveals that the FGFR1 kinase can reach a fully active configuration through phosphorylation and ATP binding. A 3-10 helix formation in the activation loop signifies its rearrangement leading to stability upon ATP binding. The interaction of phosphorylated tyrosine (pTyr654) with positively charged residues forms strong salt-bridge interactions, driving the compactness of the structure. The dynamic cross-correlation map reveals phosphorylation enhances correlated motions and reduces anti-correlated motions between different domains. We believe that the mechanistic understanding of large-conformational changes upon the activation of the FGFR1 kinase will aid the development of novel targeted therapeutics.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Subhasmita Mahapatra
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Suman Koirala
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| |
Collapse
|
7
|
Nicolle A, Deng S, Ihme M, Kuzhagaliyeva N, Ibrahim EA, Farooq A. Mixtures Recomposition by Neural Nets: A Multidisciplinary Overview. J Chem Inf Model 2024; 64:597-620. [PMID: 38284618 DOI: 10.1021/acs.jcim.3c01633] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2024]
Abstract
Artificial Neural Networks (ANNs) are transforming how we understand chemical mixtures, providing an expressive view of the chemical space and multiscale processes. Their hybridization with physical knowledge can bridge the gap between predictivity and understanding of the underlying processes. This overview explores recent progress in ANNs, particularly their potential in the 'recomposition' of chemical mixtures. Graph-based representations reveal patterns among mixture components, and deep learning models excel in capturing complexity and symmetries when compared to traditional Quantitative Structure-Property Relationship models. Key components, such as Hamiltonian networks and convolution operations, play a central role in representing multiscale mixtures. The integration of ANNs with Chemical Reaction Networks and Physics-Informed Neural Networks for inverse chemical kinetic problems is also examined. The combination of sensors with ANNs shows promise in optical and biomimetic applications. A common ground is identified in the context of statistical physics, where ANN-based methods iteratively adapt their models by blending their initial states with training data. The concept of mixture recomposition unveils a reciprocal inspiration between ANNs and reactive mixtures, highlighting learning behaviors influenced by the training environment.
Collapse
Affiliation(s)
- Andre Nicolle
- Aramco Fuel Research Center, Rueil-Malmaison 92852, France
| | - Sili Deng
- Massachusetts Institute of Technology, Cambridge 02139, Massachusetts, United States
| | - Matthias Ihme
- Stanford University, Stanford 94305, California, United States
| | | | - Emad Al Ibrahim
- King Abdullah University of Science and Technology, Thuwal 23955, Saudi Arabia
| | - Aamir Farooq
- King Abdullah University of Science and Technology, Thuwal 23955, Saudi Arabia
| |
Collapse
|
8
|
Kornev AP, Weng JH, Maillard RA, Taylor SS. Gauging Dynamics-driven Allostery Using a New Computational Tool: A CAP Case Study. J Mol Biol 2024; 436:168395. [PMID: 38097109 PMCID: PMC10851786 DOI: 10.1016/j.jmb.2023.168395] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 11/22/2023] [Accepted: 12/08/2023] [Indexed: 12/24/2023]
Abstract
In this study, we utilize Protein Residue Networks (PRNs), constructed using Local Spatial Pattern (LSP) alignment, to explore the dynamic behavior of Catabolite Activator Protein (CAP) upon the sequential binding of cAMP. We employed the Degree Centrality of these PRNs to investigate protein dynamics on a sub-nanosecond time scale, hypothesizing that it would reflect changes in CAP's entropy related to its thermal motions. We show that the binding of the first cAMP led to an increase in stability in the Cyclic-Nucleotide Binding Domain A (CNBD-A) and destabilization in CNBD-B, agreeing with previous reports explaining the negative cooperativity of cAMP binding in terms of an entropy-driven allostery. LSP-based PRNs also allow for the study of Betweenness Centrality, another graph-theoretical characteristic of PRNs, providing insights into global residue connectivity within CAP. Using this approach, we were able to correctly identify amino acids that were shown to be critical in mediating allosteric interactions in CAP. The agreement between our studies and previous experimental reports validates our method, particularly with respect to the reliability of Degree Centrality as a proxy for entropy related to protein thermal dynamics. Because LSP-based PRNs can be easily extended to include dynamics of small organic molecules, polynucleotides, or other allosteric proteins, the methods presented here mark a significant advancement in the field, positioning them as vital tools for a fast, cost-effective, and accurate analysis of entropy-driven allostery and identification of allosteric hotspots.
Collapse
Affiliation(s)
- Alexandr P Kornev
- Departmen of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA.
| | - Jui-Hung Weng
- Departmen of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA
| | - Rodrigo A Maillard
- Department of Chemistry, Georgetown University, Washington, DC 20007, USA
| | - Susan S Taylor
- Departmen of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA; Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| |
Collapse
|
9
|
Draškovič-Bračun A, Potisk T, Svenšek D. Modeling ultrasonic metafluids: The significance of discrete oscillators. Phys Rev E 2024; 109:014604. [PMID: 38366442 DOI: 10.1103/physreve.109.014604] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 12/11/2023] [Indexed: 02/18/2024]
Abstract
Significant changes in the acoustic response of a fluid can be induced by the suspension of tiny, subwavelength-size discrete micro-oscillators in the fluid. We investigate how the topological properties of these oscillators, such as the mass distribution and connectivity of the oscillator parts, influence the effective dynamic density and compressibility of the fluid in which they are embedded. We demonstrate a superior, metamaterial-like response of the suspension when using micro-oscillators with a high density of low-frequency modes. Such low-frequency modes occur in loosely connected microstructures and make the system much more experimentally feasible due to the larger ultrasonic attenuation length at these frequencies. In addition, the absence of need for an intricately designed structure brings experimental implementation within reach.
Collapse
Affiliation(s)
- Aljaž Draškovič-Bračun
- Laboratory of Molecular Modeling, National Institute of Chemistry, SI-1001 Ljubljana, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| | - Tilen Potisk
- Laboratory of Molecular Modeling, National Institute of Chemistry, SI-1001 Ljubljana, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| | - Daniel Svenšek
- Laboratory of Molecular Modeling, National Institute of Chemistry, SI-1001 Ljubljana, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
10
|
Huang A, Lu F, Liu F. Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network. Protein Sci 2023; 32:e4837. [PMID: 37984374 PMCID: PMC10682693 DOI: 10.1002/pro.4837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2023] [Revised: 10/18/2023] [Accepted: 11/09/2023] [Indexed: 11/22/2023]
Abstract
Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.
Collapse
Affiliation(s)
- Ailan Huang
- College of BiotechnologyTianjin University of Science & TechnologyTianjinChina
| | - Fuping Lu
- College of BiotechnologyTianjin University of Science & TechnologyTianjinChina
- Key Laboratory of Industrial Fermentation Microbiology, Ministry of EducationTianjin Key Laboratory of Industrial MicrobiologyTianjinChina
| | - Fufeng Liu
- College of BiotechnologyTianjin University of Science & TechnologyTianjinChina
- Key Laboratory of Industrial Fermentation Microbiology, Ministry of EducationTianjin Key Laboratory of Industrial MicrobiologyTianjinChina
| |
Collapse
|
11
|
Rosignoli S, di Paola L, Paiardini A. PyPCN: protein contact networks in PyMOL. Bioinformatics 2023; 39:btad675. [PMID: 37941462 PMCID: PMC10641099 DOI: 10.1093/bioinformatics/btad675] [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: 04/27/2023] [Revised: 09/25/2023] [Accepted: 11/03/2023] [Indexed: 11/10/2023] Open
Abstract
MOTIVATION Protein contact networks (PCNs) represent the 3D structure of a protein using network formalism. Inter-residue contacts are described as binary adjacency matrices, which are derived from the graph representation of residues (as α-carbons, β-carbons or centroids) and Euclidean distances according to defined thresholds. Functional characterization algorithms are computed on binary adjacency matrices to unveil allosteric, dynamic, and interaction mechanisms in proteins. Such strategies are usually applied in a combinatorial manner, although rarely in seamless and user-friendly implementations. RESULTS PyPCN is a plugin for PyMOL wrapping more than twenty PCN algorithms and metrics in an easy-to-use graphical user interface, to support PCN analysis. The plugin accepts 3D structures from the Protein Data Bank, user-provided PDBs, or precomputed adjacency matrices. The results are directly mapped to 3D protein structures and organized into interactive diagrams for their visualization. A dedicated graphical user interface combined with PyMOL visual support makes analysis more intuitive and easier, extending the applicability of PCNs. AVAILABILITY AND IMPLEMENTATION https://github.com/pcnproject/PyPCN.
Collapse
Affiliation(s)
- Serena Rosignoli
- Department of Biochemical Sciences “A. Rossi Fanelli”, Sapienza University of Rome, 00185 Rome, Italy
| | - Luisa di Paola
- Unit of Chemical-Physics Fundamentals in Chemical Engineering, Department of Engineering, Università Campus Bio-Medico di Roma, 00128 Rome, Italy
| | - Alessandro Paiardini
- Department of Biochemical Sciences “A. Rossi Fanelli”, Sapienza University of Rome, 00185 Rome, Italy
| |
Collapse
|
12
|
Burns D, Venditti V, Potoyan DA. Temperature sensitive contact modes allosterically gate TRPV3. PLoS Comput Biol 2023; 19:e1011545. [PMID: 37831724 PMCID: PMC10599574 DOI: 10.1371/journal.pcbi.1011545] [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: 06/30/2023] [Revised: 10/25/2023] [Accepted: 09/27/2023] [Indexed: 10/15/2023] Open
Abstract
TRPV Ion channels are sophisticated molecular sensors designed to respond to distinct temperature thresholds. The recent surge in cryo-EM structures has provided numerous insights into the structural rearrangements accompanying their opening and closing; however, the molecular mechanisms by which TRPV channels establish precise and robust temperature sensing remain elusive. In this work we employ molecular simulations, multi-ensemble contact analysis, graph theory, and machine learning techniques to reveal the temperature-sensitive residue-residue interactions driving allostery in TRPV3. We find that groups of residues exhibiting similar temperature-dependent contact frequency profiles cluster at specific regions of the channel. The dominant mode clusters on the ankyrin repeat domain and displays a linear melting trend while others display non-linear trends. These modes describe the residue-level temperature response patterns that underlie the channel's functional dynamics. With network analysis, we find that the community structure of the channel changes with temperature. And that a network of high centrality contacts connects distant regions of the protomer to the gate, serving as a means for the temperature-sensitive contact modes to allosterically regulate channel gating. Using a random forest model, we show that the contact states of specific temperature-sensitive modes are indeed predictive of the channel gate's state. Supporting the physical validity of these modes and networks are several residues identified with our analyses that are reported in literature to be functionally critical. Our results offer high resolution insight into thermo-TRP channel function and demonstrate the utility of temperature-sensitive contact analysis.
Collapse
Affiliation(s)
- Daniel Burns
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, Iowa, United States of America
| | - Vincenzo Venditti
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, Iowa, United States of America
| | - Davit A. Potoyan
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, Iowa, United States of America
| |
Collapse
|
13
|
Gutland C. The shift from quantitative to qualitative thinking-problems and prospects as viewed from Husserl's and Hegel's philosophy. Front Psychol 2023; 14:1232420. [PMID: 37790237 PMCID: PMC10544961 DOI: 10.3389/fpsyg.2023.1232420] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 08/23/2023] [Indexed: 10/05/2023] Open
Abstract
This article contrasts the views of the philosophers Husserl and Hegel on quantification in science and compares their proposals for conducting rigorous qualitative research. Both deem quantification integral to science, but furthermore proposed methodologies to investigate qualitative necessities achieved by a shift in conscious activity and awareness. However, their methodologies differ significantly. While Husserl rejects idealization and instead proposes intuitive means to ideate qualitative essential relations, Hegel suggests idealizing less one-sidedly, namely, qualitatively over and above quantitatively. The article first examines how quantification is achieved and how it contrasts with measuring. This contrast reveals that measuring implies knowledge of qualities. These qualities, however, thus far remain oddly external to the mathematical relations linking the various established equations. The article then follows Husserl's reconstruction of the development of science to illustrate the dismissal of many experiential qualities and how philosophy further amplified skepticism about science on qualities. Husserl's notion of the life-world and the method of eidetic variation are then introduced as means to counterbalance mathematical proceedings in science. However, this method reveals both eidetic qualitative structures and psychical structures without being able to distinguish between them. It is thus susceptible to idiosyncratic, traditional, and cultural biases. Subsequently, Hegel's description of the shift in conscious experience that sets qualitative from quantitative thinking apart is introduced. This shift may overcome the biases, but it faces skepticism that calls for further investigation of the experience of different kinds of thinking.
Collapse
|
14
|
Burns D, Venditti V, Potoyan DA. Temperature-Sensitive Contact Modes Allosterically Gate TRPV3. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.02.522497. [PMID: 36711981 PMCID: PMC9881879 DOI: 10.1101/2023.01.02.522497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
TRPV Ion channels are sophisticated molecular sensors designed to respond to distinct temperature thresholds. The recent surge in cryo-EM structures has provided numerous insights into the structural rearrangements accompanying their opening and closing; however, the molecular mechanisms by which TRPV channels establish precise and robust temperature sensing remain elusive. In this work we employ molecular simulations, multi-ensemble contact analysis, graph theory, and machine learning techniques to reveal the temperature-sensitive residue-residue interactions driving allostery in TRPV3. We find that groups of residues exhibiting similar temperature-dependent contact frequency profiles cluster at specific regions of the channel. The dominant mode clusters on the ankyrin repeat domain and displays a linear melting trend while others display non-linear trends. These modes describe the residue-level temperature response patterns that underlie the channel's functional dynamics. With network analysis, we find that the community structure of the channel changes with temperature. And that a network of high centrality contacts connects distant regions of the protomer to the gate, serving as a means for the temperature-sensitive contact modes to allosterically regulate channel gating. Using a random forest model, we show that the contact states of specific temperature-sensitive modes are indeed predictive of the channel gate's state. Supporting the physical validity of these modes and networks are several residues identified with our analyses that are reported in literature to be functionally critical. Our results offer high resolution insight into thermo-TRP channel function and demonstrate the utility of temperature-sensitive contact analysis.
Collapse
|
15
|
Gheeraert A, Lesieur C, Batista VS, Vuillon L, Rivalta I. Connected Component Analysis of Dynamical Perturbation Contact Networks. J Phys Chem B 2023; 127:7571-7580. [PMID: 37641933 PMCID: PMC10493978 DOI: 10.1021/acs.jpcb.3c04592] [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: 07/07/2023] [Revised: 08/02/2023] [Indexed: 08/31/2023]
Abstract
Describing protein dynamical networks through amino acid contacts is a powerful way to analyze complex biomolecular systems. However, due to the size of the systems, identifying the relevant features of protein-weighted graphs can be a difficult task. To address this issue, we present the connected component analysis (CCA) approach that allows for fast, robust, and unbiased analysis of dynamical perturbation contact networks (DPCNs). We first illustrate the CCA method as applied to a prototypical allosteric enzyme, the imidazoleglycerol phosphate synthase (IGPS) enzyme from Thermotoga maritima bacteria. This approach was shown to outperform the clustering methods applied to DPCNs, which could not capture the propagation of the allosteric signal within the protein graph. On the other hand, CCA reduced the DPCN size, providing connected components that nicely describe the allosteric propagation of the signal from the effector to the active sites of the protein. By applying the CCA to the IGPS enzyme in different conditions, i.e., at high temperature and from another organism (yeast IGPS), and to a different enzyme, i.e., a protein kinase, we demonstrated how CCA of DPCNs is an effective and transferable tool that facilitates the analysis of protein-weighted networks.
Collapse
Affiliation(s)
- Aria Gheeraert
- Laboratoire
de Mathématiques (LAMA), Université
Savoie Mont Blanc, CNRS, 73376 Le Bourget du Lac, France
- Dipartimento
di Chimica Industriale “Toso Montanari”, Alma Mater
Studiorum, Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Claire Lesieur
- Univ.
Lyon, CNRS, INSA Lyon, Université Claude Bernard Lyon 1, Ecole
Centrale de Lyon, Ampère UMR5005, Villeurbanne 69622, France
- Institut
Rhônalpin des Systèmes Complexes, IXXI-ENS-Lyon, Lyon 69007, France
| | - Victor S. Batista
- Department
of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Laurent Vuillon
- Laboratoire
de Mathématiques (LAMA), Université
Savoie Mont Blanc, CNRS, 73376 Le Bourget du Lac, France
- Institut
Rhônalpin des Systèmes Complexes, IXXI-ENS-Lyon, Lyon 69007, France
| | - Ivan Rivalta
- Dipartimento
di Chimica Industriale “Toso Montanari”, Alma Mater
Studiorum, Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
- ENS
de Lyon,
CNRS, Laboratoire de Chimie UMR 5182, 69364 Lyon, France
| |
Collapse
|
16
|
Jones AC, Kornev AP, Weng JH, Manning G, Taylor SS, Newton AC. Single-residue mutation in protein kinase C toggles between cancer and neurodegeneration. Biochem J 2023; 480:1299-1316. [PMID: 37551632 PMCID: PMC10586763 DOI: 10.1042/bcj20220397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Revised: 08/04/2023] [Accepted: 08/07/2023] [Indexed: 08/09/2023]
Abstract
Conventional protein kinase C (cPKC) isozymes tune the signaling output of cells, with loss-of-function somatic mutations associated with cancer and gain-of-function germline mutations identified in neurodegeneration. PKC with impaired autoinhibition is removed from the cell by quality-control mechanisms to prevent the accumulation of aberrantly active enzyme. Here, we examine how a highly conserved residue in the C1A domain of cPKC isozymes permits quality-control degradation when mutated to histidine in cancer (PKCβ-R42H) and blocks down-regulation when mutated to proline in the neurodegenerative disease spinocerebellar ataxia (PKCγ-R41P). Using FRET-based biosensors, we determined that mutation of R42 to any residue, including lysine, resulted in reduced autoinhibition as indicated by higher basal activity and faster agonist-induced plasma membrane translocation. R42 is predicted to form a stabilizing salt bridge with E655 in the C-tail and mutation of E655, but not neighboring E657, also reduced autoinhibition. Western blot analysis revealed that whereas R42H had reduced stability, the R42P mutant was stable and insensitive to activator-induced ubiquitination and down-regulation, an effect previously observed by deletion of the entire C1A domain. Molecular dynamics (MD) simulations and analysis of stable regions of the domain using local spatial pattern (LSP) alignment suggested that P42 interacts with Q66 to impair mobility and conformation of one of the ligand-binding loops. Additional mutation of Q66 to the smaller asparagine (R42P/Q66N), to remove conformational constraints, restored degradation sensitivity. Our results unveil how disease-associated mutations of the same residue in the C1A domain can toggle between gain- or loss-of-function of PKC.
Collapse
Affiliation(s)
- Alexander C. Jones
- Department of Pharmacology, University of California, La Jolla, CA 92093, U.S.A
- Biomedical Sciences Graduate Program, University of California, La Jolla, CA 92093, U.S.A
| | - Alexandr P. Kornev
- Department of Pharmacology, University of California, La Jolla, CA 92093, U.S.A
| | - Jui-Hung Weng
- Department of Pharmacology, University of California, La Jolla, CA 92093, U.S.A
| | | | - Susan S. Taylor
- Department of Pharmacology, University of California, La Jolla, CA 92093, U.S.A
| | - Alexandra C. Newton
- Department of Pharmacology, University of California, La Jolla, CA 92093, U.S.A
| |
Collapse
|
17
|
Sora V, Tiberti M, Beltrame L, Dogan D, Robbani SM, Rubin J, Papaleo E. PyInteraph2 and PyInKnife2 to Analyze Networks in Protein Structural Ensembles. J Chem Inf Model 2023; 63:4237-4245. [PMID: 37437128 DOI: 10.1021/acs.jcim.3c00574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/14/2023]
Abstract
Due to the complex nature of noncovalent interactions and their long-range effects, analyzing protein conformations using network theory can be enlightening. Protein Structure Networks (PSNs) provide a convenient formalism to study protein structures in relation to essential properties such as key residues for structural stability, allosteric communication, and the effects of modifications of the protein. PSNs can be defined according to very different principles, and the available tools have limitations in input formats, supported models, and version control. Other outstanding problems are related to the definition of network cutoffs and the assessment of the stability of the network properties. The protein science community could benefit from a common framework to carry out these analyses and make them easier to reproduce, reuse, and evaluate. We here provide two open-source software packages, PyInteraph2 and PyInKnife2, to implement and analyze PSNs in a reproducible and documented manner. PyInteraph2 interfaces with multiple formats for protein ensembles and incorporates different network models with the possibility of integrating them into a macronetwork and performing various downstream analyses, including hubs, connected components, and several other centrality measures, and visualizes the networks or further analyzes them thanks to compatibility with Cytoscape.PyInKnife2 that supports the network models implemented in PyInteraph2. It employs a jackknife resampling approach to estimate the convergence of network properties and streamline the selection of distance cutoffs. We foresee that the modular structure of the code and the supported version control system will promote the transition to a community-driven effort, boost reproducibility, and establish common protocols in the PSN field. As developers, we will guarantee the introduction of new functionalities and maintenance, assistance, and training of new contributors.
Collapse
Affiliation(s)
- Valentina Sora
- Cancer Structural Biology, Danish Cancer Institute, Strandboulevarden 49, 2100 Copenhagen, Denmark
- Cancer Systems Biology, Section of Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Matteo Tiberti
- Cancer Structural Biology, Danish Cancer Institute, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Ludovica Beltrame
- Cancer Systems Biology, Section of Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Deniz Dogan
- Cancer Structural Biology, Danish Cancer Institute, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Shahriyar Mahdi Robbani
- Cancer Structural Biology, Danish Cancer Institute, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Joshua Rubin
- Cancer Structural Biology, Danish Cancer Institute, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Elena Papaleo
- Cancer Structural Biology, Danish Cancer Institute, Strandboulevarden 49, 2100 Copenhagen, Denmark
- Cancer Systems Biology, Section of Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| |
Collapse
|
18
|
Lomoio U, Puccio B, Tradigo G, Guzzi PH, Veltri P. SARS-CoV-2 protein structure and sequence mutations: Evolutionary analysis and effects on virus variants. PLoS One 2023; 18:e0283400. [PMID: 37471335 PMCID: PMC10358949 DOI: 10.1371/journal.pone.0283400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 07/04/2023] [Indexed: 07/22/2023] Open
Abstract
The structure and sequence of proteins strongly influence their biological functions. New models and algorithms can help researchers in understanding how the evolution of sequences and structures is related to changes in functions. Recently, studies of SARS-CoV-2 Spike (S) protein structures have been performed to predict binding receptors and infection activity in COVID-19, hence the scientific interest in the effects of virus mutations due to sequence, structure and vaccination arises. However, there is the need for models and tools to study the links between the evolution of S protein sequence, structure and functions, and virus transmissibility and the effects of vaccination. As studies on S protein have been generated a large amount of relevant information, we propose in this work to use Protein Contact Networks (PCNs) to relate protein structures with biological properties by means of network topology properties. Topological properties are used to compare the structural changes with sequence changes. We find that both node centrality and community extraction analysis can be used to relate protein stability and functionality with sequence mutations. Starting from this we compare structural evolution to sequence changes and study mutations from a temporal perspective focusing on virus variants. Finally by applying our model to the Omicron variant we report a timeline correlation between Omicron and the vaccination campaign.
Collapse
Affiliation(s)
- Ugo Lomoio
- Department of Surgical and Medical Sciences, University of Catanzaro, Catanzaro, Italy
| | - Barbara Puccio
- Department of Surgical and Medical Sciences, University of Catanzaro, Catanzaro, Italy
| | | | - Pietro Hiram Guzzi
- Department of Surgical and Medical Sciences, University of Catanzaro, Catanzaro, Italy
| | | |
Collapse
|
19
|
Maschietto F, Allen B, Kyro GW, Batista VS. MDiGest: A Python package for describing allostery from molecular dynamics simulations. J Chem Phys 2023; 158:215103. [PMID: 37272574 PMCID: PMC10769569 DOI: 10.1063/5.0140453] [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: 12/28/2022] [Accepted: 04/04/2023] [Indexed: 06/06/2023] Open
Abstract
Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.
Collapse
Affiliation(s)
- Federica Maschietto
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Brandon Allen
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Gregory W. Kyro
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Victor S. Batista
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| |
Collapse
|
20
|
Franke L, Peter C. Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding. J Chem Theory Comput 2023; 19:2985-2995. [PMID: 37122117 DOI: 10.1021/acs.jctc.2c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.
Collapse
Affiliation(s)
- Leon Franke
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
- Konstanz Research School Chemical Biology, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| | - Christine Peter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| |
Collapse
|
21
|
Sawang N, Phongphanphanee S, Wong-ekkabut J, Sutthibutpong T. Biophysical Interpretation of Evolutionary Consequences on the SARS-CoV2 Main Protease through Molecular Dynamics Simulations and Network Topology Analysis. J Phys Chem B 2023; 127:2331-2343. [PMID: 36913683 PMCID: PMC10022058 DOI: 10.1021/acs.jpcb.2c08312] [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: 11/27/2022] [Revised: 02/25/2023] [Indexed: 03/14/2023]
Abstract
In this study, we present a combined analysis procedure between atomistic molecular dynamics (MD) simulations and network topology to obtain more understanding on the evolutionary consequences on protein stability and substrate binding of the main protease enzyme of SARS-CoV2. Communicability matrices of the protein residue networks (PRNs) were extracted from MD trajectories of both Mpro enzymes in complex with the nsp8/9 peptide substrate to compare the local communicability within both proteases that would affect the enzyme function, along with biophysical details on global protein conformation, flexibility, and contribution of amino acid side chains to both intramolecular and intermolecular interactions. The analysis displayed the significance of the mutated residue 46 with the highest communicability gain to the binding pocket closure. Interestingly, the mutated residue 134 with the highest communicability loss corresponded to a local structural disruption of the adjacent peptide loop. The enhanced flexibility of the disrupted loop connecting to the catalytic residue Cys145 introduced an extra binding mode that brought the substrate in proximity and could facilitate the reaction. This understanding might provide further help in the drug development strategy against SARS-CoV2 and prove the capability of the combined techniques of MD simulations and network topology analysis as a "reverse" protein engineering tool.
Collapse
Affiliation(s)
- Nuttawat Sawang
- Theoretical
and Computational Physics Group, Department of Physics, King Mongkut’s University of Technology Thonburi
(KMUTT), 126 Pracha-Uthit Road, Bang Mod, Thrung Khru, Bangkok 10140, Thailand
- Center
of Excellence in Theoretical and Computational Science (TaCS-CoE),
Faculty of Science, King Mongkut’s
University of Technology Thonburi (KMUTT), 126 Pracha Uthit Rd., Bang Mod, Thung Khru, Bangkok 10140, Thailand
| | - Saree Phongphanphanee
- Computational
Biomodelling Laboratory for Agricultural Science and Technology (CBLAST),
Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
- Thailand
Center of Excellence in Physics (ThEP Center), Ministry of Higher Education, Science, Research and Innovation, Bangkok 10400, Thailand
- Department
of Materials Science, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Jirasak Wong-ekkabut
- Computational
Biomodelling Laboratory for Agricultural Science and Technology (CBLAST),
Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
- Thailand
Center of Excellence in Physics (ThEP Center), Ministry of Higher Education, Science, Research and Innovation, Bangkok 10400, Thailand
- Department
of Physics, Faculty of Science, Kasetsart
University, Bangkok 10900, Thailand
| | - Thana Sutthibutpong
- Theoretical
and Computational Physics Group, Department of Physics, King Mongkut’s University of Technology Thonburi
(KMUTT), 126 Pracha-Uthit Road, Bang Mod, Thrung Khru, Bangkok 10140, Thailand
- Center
of Excellence in Theoretical and Computational Science (TaCS-CoE),
Faculty of Science, King Mongkut’s
University of Technology Thonburi (KMUTT), 126 Pracha Uthit Rd., Bang Mod, Thung Khru, Bangkok 10140, Thailand
- Computational
Biomodelling Laboratory for Agricultural Science and Technology (CBLAST),
Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| |
Collapse
|
22
|
Jones AC, Kornev AP, Weng JH, Manning G, Taylor SS, Newton AC. Single-residue mutation in protein kinase C toggles between cancer and neurodegeneration. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.16.532226. [PMID: 36993163 PMCID: PMC10055082 DOI: 10.1101/2023.03.16.532226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Conventional protein kinase C (PKC) isozymes tune the signaling output of cells, with loss-of-function somatic mutations associated with cancer and gain-of-function germline mutations identified in neurodegeneration. PKC with impaired autoinhibition is removed from the cell by quality-control mechanisms to prevent accumulation of aberrantly active enzyme. Here, we examine how a single residue in the C1A domain of PKCβ, arginine 42 (R42), permits quality-control degradation when mutated to histidine in cancer (R42H) and blocks downregulation when mutated to proline in the neurodegenerative disease spinocerebellar ataxia (R42P). Using FRET-based biosensors, we determined that mutation of R42 to any residue, including lysine, resulted in reduced autoinhibition as indicated by higher basal activity and faster agonist-induced plasma membrane translocation. R42 is predicted to form a stabilizing salt bridge with E655 in the C-tail and mutation of E655, but not neighboring E657, also reduced autoinhibition. Western blot analysis revealed that whereas R42H had reduced stability, the R42P mutant was stable and insensitive to activator-induced ubiquitination and downregulation, an effect previously observed by deletion of the entire C1A domain. Molecular dynamics (MD) simulations and analysis of stable regions of the domain using local spatial pattern (LSP) alignment suggested that P42 interacts with Q66 to impair mobility and conformation of one of the ligand-binding loops. Additional mutation of Q66 to the smaller asparagine (R42P/Q66N), to remove conformational constraints, restored degradation sensitivity to that of WT. Our results unveil how disease-associated mutations of the same residue in the C1A domain can toggle between gain- or loss-of-function of PKC.
Collapse
|
23
|
Boral A, Mitra D. Heterogeneity in winged helix-turn-helix and substrate DNA interactions: Insights from theory and experiments. J Cell Biochem 2023; 124:337-358. [PMID: 36715571 DOI: 10.1002/jcb.30369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2022] [Revised: 12/29/2022] [Accepted: 01/02/2023] [Indexed: 01/31/2023]
Abstract
Specific interactions between transcription factors (TFs) and substrate DNA constitute the fundamental basis of gene expression. Unlike in TFs like basic helix-loop-helix or basic leucine zippers, prediction of substrate DNA is extremely challenging for helix-turn-helix (HTH). Experimental techniques like chromatin immunoprecipitation combined with massively parallel DNA sequencing remains a viable option. We characterize the molecular basis of heterogeneity in HTH-DNA interaction using in silico tools and thence validate them experimentally. Given the profound functional diversity in HTH, we focus primarily on winged-HTH (wHTH). We consider 180 wHTH TFs, whose experimental three-dimensional structures are available in DNA bound/unbound conformations. Starting with PDB-wide scanning and curation of data, we construct a phylogenetic tree, which distributes 180 wHTH sequences under multiple sub-groups. Structure-sequence alignment followed by detailed intra/intergroup analysis, covariation studies and extensive network theory analysis help us to gain deep insight into heterogeneous wHTH-substrate DNA interactions. A central aim of this study is to find a consensus to predict the substrate DNA sequence for wHTH, amidst heterogeneity. The strength of our exhaustive theoretical investigations including molecular docking are successfully tested through experimental characterization of wHTH TF from Sulfurimonas denitrificans.
Collapse
Affiliation(s)
- Aparna Boral
- Department of Life Sciences, Presidency University, Kolkata, West Bengal, India
| | - Devrani Mitra
- Department of Life Sciences, Presidency University, Kolkata, West Bengal, India
| |
Collapse
|
24
|
Guzzi PH, di Paola L, Puccio B, Lomoio U, Giuliani A, Veltri P. Computational analysis of the sequence-structure relation in SARS-CoV-2 spike protein using protein contact networks. Sci Rep 2023; 13:2837. [PMID: 36808182 PMCID: PMC9936485 DOI: 10.1038/s41598-023-30052-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Accepted: 02/15/2023] [Indexed: 02/19/2023] Open
Abstract
The structure of proteins impacts directly on the function they perform. Mutations in the primary sequence can provoke structural changes with consequent modification of functional properties. SARS-CoV-2 proteins have been extensively studied during the pandemic. This wide dataset, related to sequence and structure, has enabled joint sequence-structure analysis. In this work, we focus on the SARS-CoV-2 S (Spike) protein and the relations between sequence mutations and structure variations, in order to shed light on the structural changes stemming from the position of mutated amino acid residues in three different SARS-CoV-2 strains. We propose the use of protein contact network (PCN) formalism to: (i) obtain a global metric space and compare various molecular entities, (ii) give a structural explanation of the observed phenotype, and (iii) provide context dependent descriptors of single mutations. PCNs have been used to compare sequence and structure of the Alpha, Delta, and Omicron SARS-CoV-2 variants, and we found that omicron has a unique mutational pattern leading to different structural consequences from mutations of other strains. The non-random distribution of changes in network centrality along the chain has allowed to shed light on the structural (and functional) consequences of mutations.
Collapse
Affiliation(s)
- Pietro Hiram Guzzi
- Department of Surgical and Medical Sciences, Magna Graecia University of Catanzaro, Catanzaro, Italy.
| | - Luisa di Paola
- grid.9657.d0000 0004 1757 5329Unit of Chemical-Physics Fundamentals in Chemical Engineering, Department of Engineering, Universita Campus Bio-Medico di Roma, via Alvaro del Portillo 21, 00128 Rome, Italy
| | - Barbara Puccio
- grid.411489.10000 0001 2168 2547Department of Surgical and Medical Sciences, Magna Graecia University of Catanzaro, Catanzaro, Italy
| | - Ugo Lomoio
- grid.411489.10000 0001 2168 2547Department of Surgical and Medical Sciences, Magna Graecia University of Catanzaro, Catanzaro, Italy
| | - Alessandro Giuliani
- grid.416651.10000 0000 9120 6856Environment and Health Department, Istituto Superiore di Sanita, Rome, Italy
| | - Pierangelo Veltri
- grid.411489.10000 0001 2168 2547Department of Surgical and Medical Sciences, Magna Graecia University of Catanzaro, Catanzaro, Italy ,grid.7778.f0000 0004 1937 0319Department of Computer, Modeling, Electronics and System Engineering, University of Calabria, Rende, Italy
| |
Collapse
|
25
|
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
|
26
|
Mathew AT, Baidya ATK, Das B, Devi B, Kumar R. N-glycosylation induced changes in tau protein dynamics reveal its role in tau misfolding and aggregation: A microsecond long molecular dynamics study. Proteins 2023; 91:147-160. [PMID: 36029032 DOI: 10.1002/prot.26417] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 08/18/2022] [Accepted: 08/23/2022] [Indexed: 01/07/2023]
Abstract
Various posttranslational modifications like hyperphosphorylation, O-GlcNAcylation, and acetylation have been attributed to induce the abnormal folding in tau protein. Recent in vitro studies revealed the possible involvement of N-glycosylation of tau protein in the abnormal folding and tau aggregation. Hence, in this study, we performed a microsecond long all atom molecular dynamics simulation to gain insights into the effects of N-glycosylation on Asn-359 residue which forms part of the microtubule binding region. Trajectory analysis of the stimulations coupled with essential dynamics and free energy landscape analysis suggested that tau, in its N-glycosylated form tends to exist in a largely folded conformation having high beta sheet propensity as compared to unmodified tau which exists in a large extended form with very less beta sheet propensity. Residue interaction network analysis of the lowest energy conformations further revealed that Phe378 and Lys353 are the functionally important residues in the peptide which helped in initiating the folding process and Phe378, Lys347, and Lys370 helped to maintain the stability of the protein in the folded state.
Collapse
Affiliation(s)
- Alen T Mathew
- Department of Pharmaceutical Engineering and Technology, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh, India
| | - Anurag T K Baidya
- Department of Pharmaceutical Engineering and Technology, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh, India
| | - Bhanuranjan Das
- Department of Pharmaceutical Engineering and Technology, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh, India
| | - Bharti Devi
- Department of Pharmaceutical Engineering and Technology, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh, India
| | - Rajnish Kumar
- Department of Pharmaceutical Engineering and Technology, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh, India
| |
Collapse
|
27
|
Liu J, Amaral LAN, Keten S. A new approach for extracting information from protein dynamics. Proteins 2023; 91:183-195. [PMID: 36094321 PMCID: PMC9844508 DOI: 10.1002/prot.26421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 08/25/2022] [Accepted: 09/06/2022] [Indexed: 01/19/2023]
Abstract
Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.
Collapse
Affiliation(s)
- Jenny Liu
- Department of Mechanical Engineering, Northwestern University
| | - Luís A. N. Amaral
- Department of Chemical and Biological Engineering, Northwestern University
| | - Sinan Keten
- Department of Mechanical Engineering, Northwestern University
| |
Collapse
|
28
|
Structural analysis of SARS-CoV-2 Spike protein variants through graph embedding. NETWORK MODELING AND ANALYSIS IN HEALTH INFORMATICS AND BIOINFORMATICS 2023; 12:3. [PMID: 36506261 PMCID: PMC9718452 DOI: 10.1007/s13721-022-00397-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 10/21/2022] [Accepted: 11/16/2022] [Indexed: 12/03/2022]
Abstract
Since December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected almost all countries. The unprecedented spreading of this virus has led to the insurgence of many variants that impact protein sequence and structure that need continuous monitoring and analysis of the sequences to understand the genetic evolution and to prevent possible dangerous outcomes. Some variants causing the modification of the structure of the proteins, such as the Spike protein S, need to be monitored. Protein contact networks (PCNs) have been recently proposed as a modelling framework for protein structures. In such a framework, the protein structure is represented as an unweighted graph whose nodes are the central atoms of the backbones (C- α ), and edges connect two atoms falling in the spatial distance between 4 and 7 Å. PCN may also be a data-rich representation since we may add to each node/atom biological and topological information. Such formalism enables the possibility of using algorithms from graph theory to analyze the graph. In particular, we refer to graph embedding methods enabling the analysis of such graphs with deep learning methods. In this work, we explore the possibility of embedding PCN using Graph Neural Networks and then analyze in the embedded space each residue to distinguish mutated residues from non-mutated ones. In particular, we analyzed the structure of the Spike protein of the coronavirus. First, we obtained the PCNs of the Spike protein for the wild-type, α , β , and δ variants. Then we used the GraphSage embedding algorithm to obtain an unsupervised embedding. Then we analyzed the point of mutation in the embedded space. Results show the characteristics of the mutation point in the embedding space.
Collapse
|
29
|
Erba F, Di Paola L, Di Venere A, Mastrangelo E, Cossu F, Mei G, Minicozzi V. Head or tail? A molecular dynamics approach to the complex structure of TNF-associated factor TRAF2. Biomol Concepts 2023; 14:bmc-2022-0031. [PMID: 37377424 DOI: 10.1515/bmc-2022-0031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 06/09/2023] [Indexed: 06/29/2023] Open
Abstract
Tumor necrosis factor receptor-associated factor proteins (TRAFs) are trimeric proteins that play a fundamental role in signaling, acting as intermediaries between the tumor necrosis factor (TNF) receptors and the proteins that transmit the downstream signal. The monomeric subunits of all the TRAF family members share a common tridimensional structure: a C-terminal globular domain and a long coiled-coil tail characterizing the N-terminal section. In this study, the dependence of the TRAF2 dynamics on the length of its tail was analyzed in silico. In particular, we used the available crystallographic structure of a C-terminal fragment of TRAF2 (168 out of 501 a.a.), TRAF2-C, and that of a longer construct, addressed as TRAF2-plus, that we have re-constructed using the AlphaFold2 code. The results indicate that the longer N-terminal tail of TRAF2-plus has a strong influence on the dynamics of the globular regions in the protein C-terminal head. In fact, the quaternary interactions among the TRAF2-C subunits change asymmetrically in time, while the movements of TRAF2-plus monomers are rather limited and more ordered than those of the shorter construct. Such findings shed a new light on the dynamics of TRAF subunits and on the protein mechanism in vivo, since TRAF monomer-trimer equilibrium is crucial for several reasons (receptor recognition, membrane binding, hetero-oligomerization).
Collapse
Affiliation(s)
- Fulvio Erba
- Department of Clinical Science and Translational Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00133 Rome, Italy
| | - Luisa Di Paola
- Unit of Chemical-Physics Fundamentals in Chemical Engineering, Department of Engineering, University Campus Bio-Medico of Rome, Via Álvaro del Portillo 21, 00128 Rome, Italy
| | - Almerinda Di Venere
- Department of Experimental Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00133 Rome, Italy
| | - Eloise Mastrangelo
- National Research Council (IBF-CNR) Milan Unit, Institute of Biophysics, Via Celoria 26, 20133 Milan, Italy
| | - Federica Cossu
- National Research Council (IBF-CNR) Milan Unit, Institute of Biophysics, Via Celoria 26, 20133 Milan, Italy
| | - Giampiero Mei
- Department of Experimental Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00133 Rome, Italy
| | - Velia Minicozzi
- Department of Physics and INFN, Tor Vergata University of Rome, Via Della Ricerca Scientifica 1, 00133 Rome, Italy
| |
Collapse
|
30
|
Chen L, Gong W, Han Z, Zhou W, Yang S, Li C. Key Residues in δ Opioid Receptor Allostery Explored by the Elastic Network Model and the Complex Network Model Combined with the Perturbation Method. J Chem Inf Model 2022; 62:6727-6738. [PMID: 36073904 DOI: 10.1021/acs.jcim.2c00513] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Opioid receptors, a kind of G protein-coupled receptors (GPCRs), mainly mediate an analgesic response via allosterically transducing the signal of endogenous ligand binding in the extracellular domain to couple to effector proteins in the intracellular domain. The δ opioid receptor (DOP) is associated with emotional control besides pain control, which makes it an attractive therapeutic target. However, its allosteric mechanism and key residues responsible for the structural stability and signal communication are not completely clear. Here we utilize the Gaussian network model (GNM) and amino acid network (AAN) combined with perturbation methods to explore the issues. The constructed fcfGNMMD, where the force constants are optimized with the inverse covariance estimation based on the correlated fluctuations from the available DOP molecular dynamics (MD) ensemble, shows a better performance than traditional GNM in reproducing residue fluctuations and cross-correlations and in capturing functionally low-frequency modes. Additionally, fcfGNMMD can consider implicitly the environmental effects to some extent. The lowest mode can well divide DOP segments and identify the two sodium ion (important allosteric regulator) binding coordination shells, and from the fastest modes, the key residues important for structure stabilization are identified. Using fcfGNMMD combined with a dynamic perturbation-response method, we explore the key residues related to the sodium ion binding. Interestingly, we identify not only the key residues in sodium ion binding shells but also the ones far away from the perturbation sites, which are involved in binding with DOP ligands, suggesting the possible long-range allosteric modulation of sodium binding for the ligand binding to DOP. Furthermore, utilizing the weighted AAN combined with attack perturbations, we identify the key residues for allosteric communication. This work helps strengthen the understanding of the allosteric communication mechanism in δ opioid receptor and can provide valuable information for drug design.
Collapse
Affiliation(s)
- Lei Chen
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Wenxue Zhou
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Shuang Yang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
31
|
Tajana M, Trovato A, Tiana G. Key interaction patterns in proteins revealed by cluster expansion of the partition function. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2022; 45:95. [PMID: 36447074 DOI: 10.1140/epje/s10189-022-00250-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Accepted: 11/19/2022] [Indexed: 06/16/2023]
Abstract
The native conformation of structured proteins is stabilized by a complex network of interactions. We analyzed the elementary patterns that constitute such network and ranked them according to their importance in shaping protein sequence design. To achieve this goal, we employed a cluster expansion of the partition function in the space of sequences and evaluated numerically the statistical importance of each cluster. An important feature of this procedure is that it is applied to a dense finite system. We found that patterns that contribute most to the partition function are cycles with even numbers of nodes, while cliques are typically detrimental. Each cluster also gives a contribute to the sequence entropy, which is a measure of the evolutionary designability of a fold. We compared the entropies associated with different interaction patterns to their abundances in the native structures of real proteins.
Collapse
Affiliation(s)
- Matteo Tajana
- Department of Physics, Università degli Studi di Milano, Via Celoria 16, 20133, Milan, Italy
| | - Antonio Trovato
- Department of Physics and Astronomy "G. Galilei", Università degli Studi di Padova and INFN, Via Marzolo 8, 35121, Padova, Italy
| | - Guido Tiana
- Department of Physics and Center for Complexity and Biosystems, Università degli Studi di Milano and INFN, Via Celoria 16, 20133, Milan, Italy.
| |
Collapse
|
32
|
Rosignoli S, Paiardini A. Boosting the Full Potential of PyMOL with Structural Biology Plugins. Biomolecules 2022; 12:biom12121764. [PMID: 36551192 PMCID: PMC9775141 DOI: 10.3390/biom12121764] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Revised: 11/23/2022] [Accepted: 11/24/2022] [Indexed: 11/29/2022] Open
Abstract
Over the past few decades, the number of available structural bioinformatics pipelines, libraries, plugins, web resources and software has increased exponentially and become accessible to the broad realm of life scientists. This expansion has shaped the field as a tangled network of methods, algorithms and user interfaces. In recent years PyMOL, widely used software for biomolecules visualization and analysis, has started to play a key role in providing an open platform for the successful implementation of expert knowledge into an easy-to-use molecular graphics tool. This review outlines the plugins and features that make PyMOL an eligible environment for supporting structural bioinformatics analyses.
Collapse
|
33
|
Contreras-Torres E, Marrero-Ponce Y, Terán JE, Agüero-Chapin G, Antunes A, García-Jacas CR. Fuzzy spherical truncation-based multi-linear protein descriptors: From their definition to application in structural-related predictions. Front Chem 2022; 10:959143. [PMID: 36277354 PMCID: PMC9585278 DOI: 10.3389/fchem.2022.959143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2022] [Accepted: 08/15/2022] [Indexed: 11/13/2022] Open
Abstract
This study introduces a set of fuzzy spherically truncated three-dimensional (3D) multi-linear descriptors for proteins. These indices codify geometric structural information from kth spherically truncated spatial-(dis)similarity two-tuple and three-tuple tensors. The coefficients of these truncated tensors are calculated by applying a smoothing value to the 3D structural encoding based on the relationships between two and three amino acids of a protein embedded into a sphere. At considering, the geometrical center of the protein matches with center of the sphere, the distance between each amino acid involved in any specific interaction and the geometrical center of the protein can be computed. Then, the fuzzy membership degree of each amino acid from an spherical region of interest is computed by fuzzy membership functions (FMFs). The truncation value is finally a combination of the membership degrees from interacting amino acids, by applying the arithmetic mean as fusion rule. Several fuzzy membership functions with diverse biases on the calculation of amino acids memberships (e.g., Z-shaped (close to the center), PI-shaped (middle region), and A-Gaussian (far from the center)) were considered as well as traditional truncation functions (e.g., Switching). Such truncation functions were comparatively evaluated by exploring: 1) the frequency of membership degrees, 2) the variability and orthogonality analyses among them based on the Shannon Entropy’s and Principal Component’s methods, respectively, and 3) the prediction performance of alignment-free prediction of protein folding rates and structural classes. These analyses unraveled the singularity of the proposed fuzzy spherically truncated MDs with respect to the classical (non-truncated) ones and respect to the MDs truncated with traditional functions. They also showed an improved prediction power by attaining an external correlation coefficient of 95.82% in the folding rate modelling and an accuracy of 100% in distinguishing structural protein classes. These outcomes are better than the ones attained by existing approaches, justifying the theoretical contribution of this report. Thus, the fuzzy spherically truncated-based protein descriptors from MuLiMs-MCoMPAs (http://tomocomd.com/mulims-mcompas) are promising alignment-free predictors for modeling protein functions and properties.
Collapse
Affiliation(s)
- Ernesto Contreras-Torres
- Grupo de Medicina Molecular y Traslacional (MeM&T), Colegio de Ciencias de la Salud (COCSA), Escuela de Medicina, Universidad San Francisco de Quito (USFQ), Quito, Pichincha, Ecuador
- Instituto de Simulación Computacional (ISC-USFQ), Quito, Pichincha, Ecuador
- BCAM—Basque Center for Applied Mathematics, Bilbao, Spain
| | - Yovani Marrero-Ponce
- Grupo de Medicina Molecular y Traslacional (MeM&T), Colegio de Ciencias de la Salud (COCSA), Escuela de Medicina, Universidad San Francisco de Quito (USFQ), Quito, Pichincha, Ecuador
- Instituto de Simulación Computacional (ISC-USFQ), Quito, Pichincha, Ecuador
- Computer-Aided Molecular “Biosilico” Discovery and Bioinformatics Research International Network (CAMD-BIR IN), Quito, Ecuador
- *Correspondence: Yovani Marrero-Ponce, , , César R. García-Jacas, , ,
| | - Julio E. Terán
- Grupo de Medicina Molecular y Traslacional (MeM&T), Colegio de Ciencias de la Salud (COCSA), Escuela de Medicina, Universidad San Francisco de Quito (USFQ), Quito, Pichincha, Ecuador
- Instituto de Simulación Computacional (ISC-USFQ), Quito, Pichincha, Ecuador
- Department of Textile Engineering, Chemistry and Science, College of Textiles, North Carolina State University, Raleigh, NC, United States
| | - Guillermin Agüero-Chapin
- CIIMAR—Centro Interdisciplinar de Investigação Marinha e Ambiental, Universidade do Porto, Porto, Portugal
- Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, Porto, Portugal
| | - Agostinho Antunes
- CIIMAR—Centro Interdisciplinar de Investigação Marinha e Ambiental, Universidade do Porto, Porto, Portugal
- Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, Porto, Portugal
| | - César R. García-Jacas
- Cátedras Conacyt—Departamento de Ciencias de la Computación, Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE), Ensenada, Baja California, Mexico
- *Correspondence: Yovani Marrero-Ponce, , , César R. García-Jacas, , ,
| |
Collapse
|
34
|
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]
|
35
|
Guzzi PH, Di Paola L, Giuliani A, Veltri P. PCN-Miner: an open-source extensible tool for the analysis of Protein Contact Networks. Bioinformatics 2022; 38:4235-4237. [PMID: 35799364 DOI: 10.1093/bioinformatics/btac450] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 06/14/2022] [Accepted: 07/04/2022] [Indexed: 12/24/2022] Open
Abstract
MOTIVATION Protein Contact Network (PCN) is a powerful method for analysing the structure and function of proteins, with a specific focus on disclosing the molecular features of allosteric regulation through the discovery of modular substructures. The importance of PCN analysis has been shown in many contexts, such as the analysis of SARS-CoV-2 Spike protein and its complexes with the Angiotensin Converting Enzyme 2 (ACE2) human receptors. Even if there exist many software tools implementing such methods, there is a growing need for the introduction of tools integrating existing approaches. RESULTS We present PCN-Miner, a software tool implemented in the Python programming language, able to (i) import protein structures from the Protein Data Bank; (ii) generate the corresponding PCN; (iii) model, analyse and visualize PCNs and related protein structures by using a set of known algorithms and metrics. The PCN-Miner can cover a large set of applications: from clustering to embedding and subsequent analysis. AVAILABILITY AND IMPLEMENTATION The PCN-Miner tool is freely available at the following GitHub repository: https://github.com/hguzzi/ProteinContactNetworks. It is also available in the Python Package Index (PyPI) repository.
Collapse
Affiliation(s)
- Pietro Hiram Guzzi
- Department of Surgical and Medical Sciences, Magna Graecia University of Catanzaro, 88100 Catanzaro, Italy
| | - Luisa Di Paola
- Unit of Chemical-Physics Fundamentals in Chemical Engineering, Department of Engineering, Università Campus Bio-Medico di Roma, 00128 Rome, Italy
| | - Alessandro Giuliani
- Environment and Health Department, Istituto Superiore di Sanità, 00161Rome, Italy
| | - Pierangelo Veltri
- Department of Surgical and Medical Sciences, Magna Graecia University of Catanzaro, 88100 Catanzaro, Italy
| |
Collapse
|
36
|
Han Z, Wu Z, Gong W, Zhou W, Chen L, Li C. Allosteric mechanism for SL RNA recognition by polypyrimidine tract binding protein RRM1: An atomistic MD simulation and network-based study. Int J Biol Macromol 2022; 221:763-772. [PMID: 36058398 DOI: 10.1016/j.ijbiomac.2022.08.181] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 08/20/2022] [Accepted: 08/27/2022] [Indexed: 12/01/2022]
Abstract
Polypyrimidine tract-binding protein (PTB), an RNA-binding protein, is involved in the regulation of diverse processes in mRNA metabolism. However, the allosteric modulation of its binding with RNA remains unclear. We explore the dynamic characteristics of PTB RNA recognition motif 1 (RRM1) in its RNA-free and wild-type/mutant RNA-bound states to understand the issues using molecular dynamics (MD) simulation, perturbation response scanning (PRS) and protein structure network (PSN) models. It is found that RNA binding strengthens RRM1 stability, while L151G mutation in α3 helix far away from the interface makes the complex unstable. The latter is caused by long-distance dynamic couplings, which makes intermolecular electrostatic and entropy energies unfavorable. The weakened couplings between interface β sheets and C-terminal parts upon mutation reveal RNA recognition is co-regulated by these regions. Interestingly, PRS analysis reveals the allostery caused by the perturbation on α3 helix has already been pre-encoded in the equilibrium dynamics of the protein structure. PSN analysis shows the details of the allosteric signal transmission, revealing the necessity of strong couplings between α3 helix and interface for maintaining the high binding affinity. This study sheds light on the mechanisms of PTB allostery and RNA recognition and can provide important information for drug design.
Collapse
Affiliation(s)
- Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Zhixiang Wu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Wenxue Zhou
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Lei Chen
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.
| |
Collapse
|
37
|
A Statistical Journey through the Topological Determinants of the β2 Adrenergic Receptor Dynamics. ENTROPY 2022; 24:e24070998. [PMID: 35885221 PMCID: PMC9318934 DOI: 10.3390/e24070998] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/28/2022] [Revised: 07/12/2022] [Accepted: 07/15/2022] [Indexed: 02/06/2023]
Abstract
Activation of G-protein-coupled receptors (GPCRs) is mediated by molecular switches throughout the transmembrane region of the receptor. In this work, we continued along the path of a previous computational study wherein energy transport in the β2 Adrenergic Receptor (β2-AR) was examined and allosteric switches were identified in the molecular structure through the reorganization of energy transport networks during activation. In this work, we further investigated the allosteric properties of β2-AR, using Protein Contact Networks (PCNs). In this paper, we report an extensive statistical analysis of the topological and structural properties of β2-AR along its molecular dynamics trajectory to identify the activation pattern of this molecular system. The results show a distinct character to the activation that both helps to understand the allosteric switching previously identified and confirms the relevance of the network formalism to uncover relevant functional features of protein molecules.
Collapse
|
38
|
Hadi-Alijanvand H, Di Paola L, Hu G, Leitner DM, Verkhivker GM, Sun P, Poudel H, Giuliani A. Biophysical Insight into the SARS-CoV2 Spike-ACE2 Interaction and Its Modulation by Hepcidin through a Multifaceted Computational Approach. ACS OMEGA 2022; 7:17024-17042. [PMID: 35600142 PMCID: PMC9113007 DOI: 10.1021/acsomega.2c00154] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2022] [Accepted: 04/15/2022] [Indexed: 05/08/2023]
Abstract
At the center of the SARS-CoV2 infection, the spike protein and its interaction with the human receptor ACE2 play a central role in the molecular machinery of SARS-CoV2 infection of human cells. Vaccine therapies are a valuable barrier to the worst effects of the virus and to its diffusion, but the need of purposed drugs is emerging as a core target of the fight against COVID19. In this respect, the repurposing of drugs has already led to discovery of drugs thought to reduce the effects of the cytokine storm, but still a drug targeting the spike protein, in the infection stage, is missing. In this work, we present a multifaceted computational approach strongly grounded on a biophysical modeling of biological systems, so to disclose the interaction of the SARS-CoV2 spike protein with ACE2 with a special focus to an allosteric regulation of the spike-ACE2 interaction. Our approach includes the following methodologies: Protein Contact Networks and Network Clustering, Targeted Molecular Dynamics, Elastic Network Modeling, Perturbation Response Scanning, and a computational analysis of energy flow and SEPAS as a protein-softness and monomer-based affinity predictor. We applied this approach to free (closed and open) states of spike protein and spike-ACE2 complexes. Eventually, we analyzed the interactions of free and bound forms of spike with hepcidin (HPC), the major hormone in iron regulation, recently addressed as a central player in the COVID19 pathogenesis, with a special emphasis to the most severe outcomes. Our results demonstrate that, compared with closed and open states, the spike protein in the ACE2-bound state shows higher allosteric potential. The correspondence between hinge sites and the Allosteric Modulation Region (AMR) in the S-ACE complex suggests a molecular basis for hepcidin involvement in COVID19 pathogenesis. We verify the importance of AMR in different states of spike and then study its interactions with HPC and the consequence of the HPC-AMR interaction on spike dynamics and its affinity for ACE2. We propose two complementary mechanisms for HPC effects on spike of SARS-CoV-2; (a) HPC acts as a competitive inhibitor when spike is in a preinfection state (open and with no ACE2), (b) the HPC-AMR interaction pushes the spike structure into the safer closed state. These findings need clear molecular in vivo verification beside clinical observations.
Collapse
Affiliation(s)
- Hamid Hadi-Alijanvand
- Department
of Biological Sciences, Institute for Advanced
Studies in Basic Sciences, Zanjan 45137-66731, Iran
| | - Luisa Di Paola
- Unit
of Chemical-Physics Fundamentals in Chemical Engineering, Department
of Engineering, Università Campus
Bio-Medico di Roma, via
Álvaro del Portillo 21, Rome 00128, Italy
| | - Guang Hu
- Center
for Systems Biology, Department of Bioinformatics, School of Biology
and Basic Medical Sciences, Soochow University, Suzhou 215123, China
- . Phone: +39 (06) 225419634
| | - David M. Leitner
- Department
of Chemistry, University of Nevada, Reno 89557, Nevada, United States
| | - Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange 92866, California, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine 92618, California, United States
| | - Peixin Sun
- Center
for Systems Biology, Department of Bioinformatics, School of Biology
and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Humanath Poudel
- Department
of Chemistry, University of Nevada, Reno 89557, Nevada, United States
| | - Alessandro Giuliani
- Environmental
and Health Department, Istituto Superiore
di Sanità, Rome 00161, Italy
| |
Collapse
|
39
|
Sk MF, Jonniya NA, Roy R, Kar P. Phosphorylation-Induced Conformational Dynamics and Inhibition of Janus Kinase 1 by Suppressors of Cytokine Signaling 1. J Phys Chem B 2022; 126:3224-3239. [PMID: 35443129 DOI: 10.1021/acs.jpcb.1c10733] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The dysfunction of the JAK/STAT (Janus kinase/signal transducers and activators of transcription) pathway results in several pathophysiological conditions, including autoimmune disorders. The negative feedback regulators of the JAK/STAT signaling pathway, suppressors of cytokine signaling (SOCS), act as a natural inhibitor of JAK and inhibit aberrant activity. SOCS1 is the most potent member of the SOCS family, whose kinase inhibitory region targets the substrate-binding groove of JAK with high affinity and blocks the phosphorylation of JAK kinases. Overall, we performed an aggregate of 13 μs molecular dynamics simulations on the activation loop's three different phosphorylation (double and single) states. Results from our simulations show that the single Tyr1034 phosphorylation could stabilize the JAK1/SOCS1 complex as well as the flexible activation segment. The phosphate-binding loop (P-loop) shows conformational variability at dual and single phosphorylated states. Principal component analysis and protein structure network (PSN) analysis reveal that the different phosphorylation states and SOCS1 binding induce intermediate inactive conformations of JAK1, which could be a better target for future JAK1 selective drug design. PSN analysis suggests that the com-pY1034 system is stabilized due to higher values of network hubs than the other two complex systems. Moreover, the binding free energy calculations suggest that pTyr1034 states show a higher affinity toward SOCS1 than the dual and pTyr1035 states. We believe that the mechanistic understanding of JAK1/SOCS1 complexation will aid future studies related to peptide inhibitors based on SOCS1.
Collapse
Affiliation(s)
- Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Rajarshi Roy
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| |
Collapse
|
40
|
Cheng Q, DeYonker NJ. A Case Study of the Glycoside Hydrolase Enzyme Mechanism Using an Automated QM-Cluster Model Building Toolkit. Front Chem 2022; 10:854318. [PMID: 35402371 PMCID: PMC8987026 DOI: 10.3389/fchem.2022.854318] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2022] [Accepted: 03/08/2022] [Indexed: 12/03/2022] Open
Abstract
Glycoside hydrolase enzymes are important for hydrolyzing the β-1,4 glycosidic bond in polysaccharides for deconstruction of carbohydrates. The two-step retaining reaction mechanism of Glycoside Hydrolase Family 7 (GH7) was explored with different sized QM-cluster models built by the Residue Interaction Network ResidUe Selector (RINRUS) software using both the wild-type protein and its E217Q mutant. The first step is the glycosylation, in which the acidic residue 217 donates a proton to the glycosidic oxygen leading to bond cleavage. In the subsequent deglycosylation step, one water molecule migrates into the active site and attacks the anomeric carbon. Residue interaction-based QM-cluster models lead to reliable structural and energetic results for proposed glycoside hydrolase mechanisms. The free energies of activation for glycosylation in the largest QM-cluster models were predicted to be 19.5 and 31.4 kcal mol−1 for the wild-type protein and its E217Q mutant, which agree with experimental trends that mutation of the acidic residue Glu217 to Gln will slow down the reaction; and are higher in free energy than the deglycosylation transition states (13.8 and 25.5 kcal mol−1 for the wild-type protein and its mutant, respectively). For the mutated protein, glycosylation led to a low-energy product. This thermodynamic sink may correspond to the intermediate state which was isolated in the X-ray crystal structure. Hence, the glycosylation is validated to be the rate-limiting step in both the wild-type and mutated enzyme.
Collapse
Affiliation(s)
- Qianyi Cheng
- *Correspondence: Qianyi Cheng, ; Nathan John DeYonker,
| | | |
Collapse
|
41
|
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
|
42
|
Symmetric versus Asymmetric Features of Homologous Homodimeric Amine Oxidases: When Water and Cavities Make the Difference. Symmetry (Basel) 2022. [DOI: 10.3390/sym14030522] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/10/2022] Open
Abstract
Symmetry is an intrinsic property of homo-oligomers. Amine oxidases are multidomain homodimeric enzymes that contain one catalytic site per subunit, and that share a high homology degree. In this paper, we investigated, by fluorescence spectroscopy measurements, the conformational dynamics and resiliency in solutions of two amine oxidases, one from lentil seedlings, and one from Euphorbia characias latex, of which the crystallographic structure is still unknown. The data demonstrate that slight but significant differences exist at the level of the local tridimensional structure, which arise from the presence of large internal cavities, which are characterized by different hydration extents. Molecular dynamics and a contact network methodology were also used to further explore, in silico, the structural features of the two proteins. The analysis demonstrates that the two proteins show similar long-range symmetrical connectivities, but that they differ in their local (intra-subunit) contact networks, which appear mostly asymmetric. These features have been interpreted to suggest a new rationale for the functioning of amino oxidases as obligate homodimers.
Collapse
|
43
|
Linking protein structural and functional change to mutation using amino acid networks. PLoS One 2022; 17:e0261829. [PMID: 35061689 PMCID: PMC8782487 DOI: 10.1371/journal.pone.0261829] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Accepted: 12/11/2021] [Indexed: 11/30/2022] Open
Abstract
The function of a protein is strongly dependent on its structure. During evolution, proteins acquire new functions through mutations in the amino-acid sequence. Given the advance in deep mutational scanning, recent findings have found functional change to be position dependent, notwithstanding the chemical properties of mutant and mutated amino acids. This could indicate that structural properties of a given position are potentially responsible for the functional relevance of a mutation. Here, we looked at the relation between structure and function of positions using five proteins with experimental data of functional change available. In order to measure structural change, we modeled mutated proteins via amino-acid networks and quantified the perturbation of each mutation. We found that structural change is position dependent, and strongly related to functional change. Strong changes in protein structure correlate with functional loss, and positions with functional gain due to mutations tend to be structurally robust. Finally, we constructed a computational method to predict functionally sensitive positions to mutations using structural change that performs well on all five proteins with a mean precision of 74.7% and recall of 69.3% of all functional positions.
Collapse
|
44
|
Yao XQ, Hamelberg D. Residue–Residue Contact Changes during Functional Processes Define Allosteric Communication Pathways. J Chem Theory Comput 2022; 18:1173-1187. [DOI: 10.1021/acs.jctc.1c00669] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
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
|
45
|
Sequence and Structure-Based Analyses of Human Ankyrin Repeats. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27020423. [PMID: 35056738 PMCID: PMC8781854 DOI: 10.3390/molecules27020423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Revised: 12/27/2021] [Accepted: 01/04/2022] [Indexed: 11/24/2022]
Abstract
Ankyrin is one of the most abundant protein repeat families found across all forms of life. It is found in a variety of multi-domain and single domain proteins in humans with diverse number of repeating units. They are observed to occur in several functionally diverse proteins, such as transcriptional initiators, cell cycle regulators, cytoskeletal organizers, ion transporters, signal transducers, developmental regulators, and toxins, and, consequently, defects in ankyrin repeat proteins have been associated with a number of human diseases. In this study, we have classified the human ankyrin proteins into clusters based on the sequence similarity in their ankyrin repeat domains. We analyzed the amino acid compositional bias and consensus ankyrin motif sequence of the clusters to understand the diversity of the human ankyrin proteins. We carried out network-based structural analysis of human ankyrin proteins across different clusters and showed the association of conserved residues with topologically important residues identified by network centrality measures. The analysis of conserved and structurally important residues helps in understanding their role in structural stability and function of these proteins. In this paper, we also discuss the significance of these conserved residues in disease association across the human ankyrin protein clusters.
Collapse
|
46
|
Baldini L, Martino A, Rizzi A. A class-specific metric learning approach for graph embedding by information granulation. Appl Soft Comput 2022. [DOI: 10.1016/j.asoc.2021.108199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
|
47
|
Abstract
The multilevel organization of nature is self-evident: proteins do interact among them to give rise to an organized metabolism and the same hierarchical organization is in action for gene expression, tissue and organ architectures, and ecological systems.The still more common approach to such state of affairs is to think that causally relevant events originate from the lower level in the form of perturbations, that climb up the hierarchy reaching the ultimate layer of macroscopic behavior (e.g., causing a specific disease). Such rigid bottom-up causative model is unable to offer realistic models of many biological phenomena.Complex network approach allows to uncover the nature of multilevel organization, but in order to operationally define the organization principles of biological systems, we need to go further and complement network approach with sensible measures of order and organization. These measures, while keeping their original physical meaning, must not impose theoretical premises not verifiable in biological frameworks. We will show here how relatively simple and largely hypothesis-free multidimensional statistics tools can satisfactorily meet these criteria.
Collapse
Affiliation(s)
- Mariano Bizzarri
- Istituto Superiore di Sanità AND Sapienza University, Environment and Health Department AND Department of Experimental Medicine, Rome, Italy
| | - Alessandro Giuliani
- Istituto Superiore di Sanità AND Sapienza University, Environment and Health Department AND Department of Experimental Medicine, Rome, Italy.
| |
Collapse
|
48
|
Martinez-Seidel F, Hsieh YC, Walther D, Kopka J, Pereira Firmino AA. [Formula: see text]: ComplexOme-Structural Network Interpreter used to study spatial enrichment in metazoan ribosomes. BMC Bioinformatics 2021; 22:605. [PMID: 34930116 PMCID: PMC8686616 DOI: 10.1186/s12859-021-04510-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2021] [Accepted: 12/01/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Upon environmental stimuli, ribosomes are surmised to undergo compositional rearrangements due to abundance changes among proteins assembled into the complex, leading to modulated structural and functional characteristics. Here, we present the ComplexOme-Structural Network Interpreter ([Formula: see text]), a computational method to allow testing whether ribosomal proteins (rProteins) that exhibit abundance changes under specific conditions are spatially confined to particular regions within the large ribosomal complex. RESULTS [Formula: see text] translates experimentally determined structures into graphs, with nodes representing proteins and edges the spatial proximity between them. In its first implementation, [Formula: see text] considers rProteins and ignores rRNA and other objects. Spatial regions are defined using a random walk with restart methodology, followed by a procedure to obtain a minimum set of regions that cover all proteins in the complex. Structural coherence is achieved by applying weights to the edges reflecting the physical proximity between purportedly contacting proteins. The weighting probabilistically guides the random-walk path trajectory. Parameter tuning during region selection provides the option to tailor the method to specific biological questions by yielding regions of different sizes with minimum overlaps. In addition, other graph community detection algorithms may be used for the [Formula: see text] workflow, considering that they yield different sized, non-overlapping regions. All tested algorithms result in the same node kernels under equivalent regions. Based on the defined regions, available abundance change information of proteins is mapped onto the graph and subsequently tested for enrichment in any of the defined spatial regions. We applied [Formula: see text] to the cytosolic ribosome structures of Saccharomyces cerevisiae, Oryctolagus cuniculus, and Triticum aestivum using datasets with available quantitative protein abundance change information. We found that in yeast, substoichiometric rProteins depleted from translating polysomes are significantly constrained to a ribosomal region close to the tRNA entry and exit sites. CONCLUSIONS [Formula: see text] offers a computational method to partition multi-protein complexes into structural regions and a statistical approach to test for spatial enrichments of any given subsets of proteins. [Formula: see text] is applicable to any multi-protein complex given appropriate structural and abundance-change data. [Formula: see text] is publicly available as a GitHub repository https://github.com/MSeidelFed/COSNet_i and can be installed using the python installer pip.
Collapse
Affiliation(s)
- Federico Martinez-Seidel
- Willmitzer Department, Max-Planck-Institute of Molecular Plant Physiology, 14476 Potsdam-Golm, Germany
- School of BioSciences, University of Melbourne, Parkville, VC 3010 Australia
| | - Yin-Chen Hsieh
- Willmitzer Department, Max-Planck-Institute of Molecular Plant Physiology, 14476 Potsdam-Golm, Germany
- Institute for Arctic and Marine Biology, UiT Arctic University of Norway, 9037 Tromsø, Norway
| | - Dirk Walther
- Willmitzer Department, Max-Planck-Institute of Molecular Plant Physiology, 14476 Potsdam-Golm, Germany
| | - Joachim Kopka
- Willmitzer Department, Max-Planck-Institute of Molecular Plant Physiology, 14476 Potsdam-Golm, Germany
| | | |
Collapse
|
49
|
Mechanical couplings of protein backbone and side chains exhibit scale-free network properties and specific hotspots for function. Comput Struct Biotechnol J 2021; 19:5309-5320. [PMID: 34765086 PMCID: PMC8554173 DOI: 10.1016/j.csbj.2021.09.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Revised: 09/02/2021] [Accepted: 09/05/2021] [Indexed: 11/23/2022] Open
Abstract
Statistical learning from protein dynamics unravels rigidities in interaction network. Backbone and side-chain mechanical couplings exhibit scale-free network properties. Graphical depiction of network rigidities captures sequence co-evolution patterns. Functional sites at secondary structure peripheries are mechanical hotspots. Our rigidity scores are compelling metrics for residue biological significance.
A backbone-side-chain elastic network model (bsENM) is devised in this contribution to decipher the network of molecular interactions during protein dynamics. The chemical details in 5 μs all-atom molecular dynamics (MD) simulation are mapped onto the bsENM spring constants by self-consistent iterations. The elastic parameters obtained by this structure-mechanics statistical learning are then used to construct inter-residue rigidity graphs for the chemical components in protein amino acids. A key discovery is that the mechanical coupling strengths of both backbone and side chains exhibit heavy-tailed distributions and scale-free network properties. In both rat trypsin and PDZ3 proteins, the statistically prominent modes of rigidity graphs uncover the sequence-specific coupling patterns and mechanical hotspots. Based on the contributions to graphical modes, our residue rigidity scores in backbone and side chains are found to be very useful metrics for the biological significance. Most functional sites have high residue rigidity scores in side chains while the biologically important glycines are generally next to mechanical hotspots. Furthermore, prominent modes in the rigidity graphs involving side chains oftentimes coincide with the co-evolution patterns due to evolutionary restraints. The bsENM specifically devised to resolve the protein chemical character thus provides useful means for extracting functional information from all-atom MD.
Collapse
|
50
|
Modulating Glycoside Hydrolase Activity between Hydrolysis and Transfer Reactions Using an Evolutionary Approach. Molecules 2021; 26:molecules26216586. [PMID: 34770995 PMCID: PMC8587830 DOI: 10.3390/molecules26216586] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2021] [Revised: 10/27/2021] [Accepted: 10/28/2021] [Indexed: 01/02/2023] Open
Abstract
The proteins within the CAZy glycoside hydrolase family GH13 catalyze the hydrolysis of polysaccharides such as glycogen and starch. Many of these enzymes also perform transglycosylation in various degrees, ranging from secondary to predominant reactions. Identifying structural determinants associated with GH13 family reaction specificity is key to modifying and designing enzymes with increased specificity towards individual reactions for further applications in industrial, chemical, or biomedical fields. This work proposes a computational approach for decoding the determinant structural composition defining the reaction specificity. This method is based on the conservation of coevolving residues in spatial contacts associated with reaction specificity. To evaluate the algorithm, mutants of α-amylase (TmAmyA) and glucanotransferase (TmGTase) from Thermotoga maritima were constructed to modify the reaction specificity. The K98P/D99A/H222Q variant from TmAmyA doubled the transglycosydation/hydrolysis (T/H) ratio while the M279N variant from TmGTase increased the hydrolysis/transglycosidation ratio five-fold. Molecular dynamic simulations of the variants indicated changes in flexibility that can account for the modified T/H ratio. An essential contribution of the presented computational approach is its capacity to identify residues outside of the active center that affect the reaction specificity.
Collapse
|