1
|
Kutlu Y, Axel G, Kolodny R, Ben-Tal N, Haliloglu T. Reused Protein Segments Linked to Functional Dynamics. Mol Biol Evol 2024; 41:msae184. [PMID: 39226145 PMCID: PMC11412252 DOI: 10.1093/molbev/msae184] [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: 02/22/2024] [Revised: 08/10/2024] [Accepted: 08/26/2024] [Indexed: 09/05/2024] Open
Abstract
Protein space is characterized by extensive recurrence, or "reuse," of parts, suggesting that new proteins and domains can evolve by mixing-and-matching of existing segments. From an evolutionary perspective, for a given combination to persist, the protein segments should presumably not only match geometrically but also dynamically communicate with each other to allow concerted motions that are key to function. Evidence from protein space supports the premise that domains indeed combine in this manner; we explore whether a similar phenomenon can be observed at the sub-domain level. To this end, we use Gaussian Network Models (GNMs) to calculate the so-called soft modes, or low-frequency modes of motion for a dataset of 150 protein domains. Modes of motion can be used to decompose a domain into segments of consecutive amino acids that we call "dynamic elements", each of which belongs to one of two parts that move in opposite senses. We find that, in many cases, the dynamic elements, detected based on GNM analysis, correspond to established "themes": Sub-domain-level segments that have been shown to recur in protein space, and which were detected in previous research using sequence similarity alone (i.e. completely independently of the GNM analysis). This statistically significant correlation hints at the importance of dynamics in evolution. Overall, the results are consistent with an evolutionary scenario where proteins have emerged from themes that need to match each other both geometrically and dynamically, e.g. to facilitate allosteric regulation.
Collapse
Affiliation(s)
- Yiğit Kutlu
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey
| | - Gabriel Axel
- School of Neurobiology, Biochemistry & Biophysics, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv, Israel
| | - Rachel Kolodny
- Department of Computer Science, University of Haifa, Haifa, Israel
| | - Nir Ben-Tal
- School of Neurobiology, Biochemistry & Biophysics, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv, Israel
| | - Turkan Haliloglu
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey
| |
Collapse
|
2
|
McBride JM, Tlusty T. AI-Predicted Protein Deformation Encodes Energy Landscape Perturbation. PHYSICAL REVIEW LETTERS 2024; 133:098401. [PMID: 39270162 DOI: 10.1103/physrevlett.133.098401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 02/27/2024] [Accepted: 07/24/2024] [Indexed: 09/15/2024]
Abstract
AI algorithms have proven to be excellent predictors of protein structure, but whether and how much these algorithms can capture the underlying physics remains an open question. Here, we aim to test this question using the Alphafold2 (AF) algorithm: We use AF to predict the subtle structural deformation induced by single mutations, quantified by strain, and compare with experimental datasets of corresponding perturbations in folding free energy ΔΔG. Unexpectedly, we find that physical strain alone-without any additional data or computation-correlates almost as well with ΔΔG as state-of-the-art energy-based and machine-learning predictors. This indicates that the AF-predicted structures alone encode fine details about the energy landscape. In particular, the structures encode significant information on stability, enough to estimate (de-)stabilizing effects of mutations, thus paving the way for the development of novel, structure-based stability predictors for protein design and evolution.
Collapse
Affiliation(s)
- John M McBride
- Center for Algorithmic and Robotized Synthesis, Institute for Basic Science, Ulsan 44919, South Korea
| | | |
Collapse
|
3
|
Kohsokabe T, Kuratanai S, Kaneko K. Developmental hourglass: Verification by numerical evolution and elucidation by dynamical-systems theory. PLoS Comput Biol 2024; 20:e1011867. [PMID: 38422161 PMCID: PMC10903806 DOI: 10.1371/journal.pcbi.1011867] [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/17/2023] [Accepted: 01/29/2024] [Indexed: 03/02/2024] Open
Abstract
Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.
Collapse
Affiliation(s)
| | | | - Kunihiko Kaneko
- Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
4
|
Martin NS, Ahnert SE. The Boltzmann distributions of molecular structures predict likely changes through random mutations. Biophys J 2023; 122:4467-4475. [PMID: 37897043 PMCID: PMC10698324 DOI: 10.1016/j.bpj.2023.10.024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 08/19/2023] [Accepted: 10/20/2023] [Indexed: 10/29/2023] Open
Abstract
New folded molecular structures can only evolve after arising through mutations. This aspect is modeled using genotype-phenotype maps, which connect sequence changes through mutations to changes in molecular structures. Previous work has shown that the likelihood of appearing through mutations can differ by orders of magnitude from structure to structure and that this can affect the outcomes of evolutionary processes. Thus, we focus on the phenotypic mutation probabilities φqp, i.e., the likelihood that a random mutation changes structure p into structure q. For both RNA secondary structures and the HP protein model, we show that a simple biophysical principle can explain and predict how this likelihood depends on the new structure q: φqp is high if sequences that fold into p as the minimum-free-energy structure are likely to have q as an alternative structure with high Boltzmann frequency. This generalizes the existing concept of plastogenetic congruence from individual sequences to the entire neutral spaces of structures. Our result helps us understand why some structural changes are more likely than others, may be useful for estimating these likelihoods via sampling and makes a connection to alternative structures with high Boltzmann frequency, which could be relevant in evolutionary processes.
Collapse
Affiliation(s)
- Nora S Martin
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, United Kingdom; Theory of Condensed Matter Group, Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom; Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom.
| | - Sebastian E Ahnert
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Cambridge, United Kingdom; The Alan Turing Institute, London, United Kingdom
| |
Collapse
|
5
|
Jordan DJ, Miska EA. Canalisation and plasticity on the developmental manifold of Caenorhabditis elegans. Mol Syst Biol 2023; 19:e11835. [PMID: 37850520 PMCID: PMC10632735 DOI: 10.15252/msb.202311835] [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: 06/21/2023] [Revised: 09/26/2023] [Accepted: 10/05/2023] [Indexed: 10/19/2023] Open
Abstract
How do the same mechanisms that faithfully regenerate complex developmental programmes in spite of environmental and genetic perturbations also allow responsiveness to environmental signals, adaptation and genetic evolution? Using the nematode Caenorhabditis elegans as a model, we explore the phenotypic space of growth and development in various genetic and environmental contexts. Our data are growth curves and developmental parameters obtained by automated microscopy. Using these, we show that among the traits that make up the developmental space, correlations within a particular context are predictive of correlations among different contexts. Furthermore, we find that the developmental variability of this animal can be captured on a relatively low-dimensional phenotypic manifold and that on this manifold, genetic and environmental contributions to plasticity can be deconvolved independently. Our perspective offers a new way of understanding the relationship between robustness and flexibility in complex systems, suggesting that projection and concentration of dimension can naturally align these forces as complementary rather than competing.
Collapse
Affiliation(s)
- David J Jordan
- Department of BiochemistryUniversity of CambridgeCambridgeUK
| | - Eric A Miska
- Department of BiochemistryUniversity of CambridgeCambridgeUK
| |
Collapse
|
6
|
Rahimi M, Taghdir M, Abasi Joozdani F. Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease. Sci Rep 2023; 13:14179. [PMID: 37648682 PMCID: PMC10469195 DOI: 10.1038/s41598-023-40818-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Accepted: 08/17/2023] [Indexed: 09/01/2023] Open
Abstract
Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.
Collapse
Affiliation(s)
- Mohammad Rahimi
- Department of Biophysics, Faculty of Biological Science, Tarbiat Modares University, Tehran, 14115_111, Iran
| | - Majid Taghdir
- Department of Biophysics, Faculty of Biological Science, Tarbiat Modares University, Tehran, 14115_111, Iran.
| | - Farzane Abasi Joozdani
- Department of Biophysics, Faculty of Biological Science, Tarbiat Modares University, Tehran, 14115_111, Iran
| |
Collapse
|
7
|
Kelly MS, Macke AC, Kahawatte S, Stump JE, Miller AR, Dima RI. The quaternary question: Determining allostery in spastin through dynamics classification learning and bioinformatics. J Chem Phys 2023; 158:125102. [PMID: 37003743 DOI: 10.1063/5.0139273] [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: 03/09/2023] Open
Abstract
The nanomachine from the ATPases associated with various cellular activities superfamily, called spastin, severs microtubules during cellular processes. To characterize the functionally important allostery in spastin, we employed methods from evolutionary information, to graph-based networks, to machine learning applied to atomistic molecular dynamics simulations of spastin in its monomeric and the functional hexameric forms, in the presence or absence of ligands. Feature selection, using machine learning approaches, for transitions between spastin states recognizes all the regions that have been proposed as allosteric or functional in the literature. The analysis of the composition of the Markov State Model macrostates in the spastin monomer, and the analysis of the direction of change in the top machine learning features for the transitions, indicate that the monomer favors the binding of ATP, which primes the regions involved in the formation of the inter-protomer interfaces for binding to other protomer(s). Allosteric path analysis of graph networks, built based on the cross-correlations between residues in simulations, shows that perturbations to a hub specific for the pre-hydrolysis hexamer propagate throughout the structure by passing through two obligatory regions: the ATP binding pocket, and pore loop 3, which connects the substrate binding site to the ATP binding site. Our findings support a model where the changes in the terminal protomers due to the binding of ligands play an active role in the force generation in spastin. The secondary structures in spastin, which are found to be highly degenerative within the network paths, are also critical for feature transitions of the classification models, which can guide the design of allosteric effectors to enhance or block allosteric signaling.
Collapse
Affiliation(s)
- Maria S Kelly
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | - Amanda C Macke
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | - Shehani Kahawatte
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | - Jacob E Stump
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | - Abigail R Miller
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| | - Ruxandra I Dima
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, USA
| |
Collapse
|
8
|
Jia K, Kilinc M, Jernigan RL. Functional Protein Dynamics Directly from Sequences. J Phys Chem B 2023; 127:1914-1921. [PMID: 36848294 PMCID: PMC10009744 DOI: 10.1021/acs.jpcb.2c05766] [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/11/2022] [Revised: 02/15/2023] [Indexed: 03/01/2023]
Abstract
The sequence correlations within a protein multiple sequence alignment are routinely being used to predict contacts within its structure, but here we point out that these data can also be used to predict a protein's dynamics directly. The elastic network protein dynamics models rely directly upon the contacts, and the normal modes of motion are obtained from the decomposition of the inverse of the contact map. To make the direct connection between sequence and dynamics, it is necessary to apply coarse-graining to the structure at the level of one point per amino acid, which has often been done, and protein coarse-grained dynamics from elastic network models has been highly successful, particularly in representing the large-scale motions of proteins that usually relate closely to their functions. The interesting implication of this is that it is not necessary to know the structure itself to obtain its dynamics and instead to use the sequence information directly to obtain the dynamics.
Collapse
Affiliation(s)
- Kejue Jia
- Bioinformatics and Computational
Biology Program and Roy J. Carver Department of Biochemistry, Biophysics
and Molecular Biology Iowa State University, Ames, Iowa 50011, United States
| | - Mesih Kilinc
- Bioinformatics and Computational
Biology Program and Roy J. Carver Department of Biochemistry, Biophysics
and Molecular Biology Iowa State University, Ames, Iowa 50011, United States
| | - Robert L. Jernigan
- Bioinformatics and Computational
Biology Program and Roy J. Carver Department of Biochemistry, Biophysics
and Molecular Biology Iowa State University, Ames, Iowa 50011, United States
| |
Collapse
|
9
|
Yildirim A, Tekpinar M. Building Quantitative Bridges between Dynamics and Sequences of SARS-CoV-2 Main Protease and a Diverse Set of Thirty-Two Proteins. J Chem Inf Model 2023; 63:9-19. [PMID: 36513349 DOI: 10.1021/acs.jcim.2c01206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Proteases are major drug targets for many viral diseases. However, mutations can render several antiprotease drugs inefficient rapidly even though these mutations may not alter protein structures significantly. Understanding relations between quickly mutating residues, protease structures, and the dynamics of the proteases is crucial for designing potent drugs. Due to this reason, we studied relations between the evolutionary information on residues in the amino acid sequences and protein dynamics for SARS-CoV-2 main protease. More precisely, we analyzed three dynamical quantities (Schlitter entropy, root-mean-square fluctuations, and dynamical flexibility index) and their relation to the amino acid conservation extracted from multiple sequence alignments of the main protease. We showed that a quantifiable similarity can be built between a sequence-based quantity called Jensen-Shannon conservation and those three dynamical quantities. We validated this similarity for a diverse set of 32 different proteins, other than the SARS-CoV-2 main protease. We believe that establishing these kinds of quantitative bridges will have larger implications for all viral proteases as well as all proteins.
Collapse
Affiliation(s)
- Ahmet Yildirim
- Department of Biology, Siirt University, 56100Siirt, Turkey
| | - Mustafa Tekpinar
- CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative - UMR 7238, Sorbonne University, 75005Paris, France
| |
Collapse
|
10
|
Maeda Y, Mizuuchi R, Shigenobu S, Shibai A, Kotani H, Furusawa C, Ichihashi N. Experimental evidence for the correlation between RNA structural fluctuations and the frequency of beneficial mutations. RNA (NEW YORK, N.Y.) 2022; 28:1659-1667. [PMID: 36195345 PMCID: PMC9670806 DOI: 10.1261/rna.079291.122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/05/2022] [Accepted: 09/28/2022] [Indexed: 06/16/2023]
Abstract
RNA has been used as a model molecule to understand the adaptive evolution process owing to the simple relationship between the structure (i.e., phenotype) and sequence (i.e., genotype). RNA usually forms multiple substructures with similar thermodynamic stabilities, called structural fluctuations. Ancel and Fontana theoretically proposed that structural fluctuation is directly related to the ease of change in structures by mutations and thus works as a source of adaptive evolution; however, experimental verification is limited. Here, we analyzed 76 RNA genotypes that appeared in our previous in vitro evolution to examine whether (i) RNA fluctuation decreases as adaptive evolution proceeds and (ii) RNAs that have larger fluctuations tend to have higher frequencies of beneficial mutations. We first computationally estimated the structural fluctuations of all RNAs and observed that they tended to decrease as their fitness increased. We next measured the frequency of beneficial mutations for 10 RNA genotypes and observed that the total number of beneficial mutations was correlated with the size of the structural fluctuations. These results consistently support the idea that the structural fluctuation of RNA, at least those evaluated in our study, works as a source of adaptive evolution.
Collapse
Affiliation(s)
- Yutaro Maeda
- Department of Life Science, Graduate School of Arts and Science, The University of Tokyo, Meguro-ku, Tokyo, 153-8902, Japan
| | - Ryo Mizuuchi
- Komaba Institute for Science, The University of Tokyo, Meguro-ku, Tokyo 153-8902, Japan
- JST, PRESTO, Kawaguchi, Saitama 332-0012, Japan
| | - Shuji Shigenobu
- National Institute for Basic Biology, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Atsushi Shibai
- Center for Biosystems Dynamics Research, RIKEN, Suita, Osaka 565-0871, Japan
| | - Hazuki Kotani
- Center for Biosystems Dynamics Research, RIKEN, Suita, Osaka 565-0871, Japan
| | - Chikara Furusawa
- Center for Biosystems Dynamics Research, RIKEN, Suita, Osaka 565-0871, Japan
- Universal Biology Institute, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Norikazu Ichihashi
- Department of Life Science, Graduate School of Arts and Science, The University of Tokyo, Meguro-ku, Tokyo, 153-8902, Japan
- Komaba Institute for Science, The University of Tokyo, Meguro-ku, Tokyo 153-8902, Japan
- Universal Biology Institute, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
11
|
McBride JM, Eckmann JP, Tlusty T. General Theory of Specific Binding: Insights from a Genetic-Mechano-Chemical Protein Model. Mol Biol Evol 2022; 39:msac217. [PMID: 36208205 PMCID: PMC9641994 DOI: 10.1093/molbev/msac217] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Proteins need to selectively interact with specific targets among a multitude of similar molecules in the cell. However, despite a firm physical understanding of binding interactions, we lack a general theory of how proteins evolve high specificity. Here, we present such a model that combines chemistry, mechanics, and genetics and explains how their interplay governs the evolution of specific protein-ligand interactions. The model shows that there are many routes to achieving molecular discrimination-by varying degrees of flexibility and shape/chemistry complementarity-but the key ingredient is precision. Harder discrimination tasks require more collective and precise coaction of structure, forces, and movements. Proteins can achieve this through correlated mutations extending far from a binding site, which fine-tune the localized interaction with the ligand. Thus, the solution of more complicated tasks is enabled by increasing the protein size, and proteins become more evolvable and robust when they are larger than the bare minimum required for discrimination. The model makes testable, specific predictions about the role of flexibility and shape mismatch in discrimination, and how evolution can independently tune affinity and specificity. Thus, the proposed theory of specific binding addresses the natural question of "why are proteins so big?". A possible answer is that molecular discrimination is often a hard task best performed by adding more layers to the protein.
Collapse
Affiliation(s)
- John M McBride
- Center for Soft and Living Matter, Institute for Basic Science, Ulsan 44919, South Korea
| | - Jean-Pierre Eckmann
- Département de Physique Théorique and Section de Mathématiques, University of Geneva, Geneva, Switzerland
| | - Tsvi Tlusty
- Center for Soft and Living Matter, Institute for Basic Science, Ulsan 44919, South Korea
- Departments of Physics and Chemistry, Ulsan National Institute of Science and Technology, Ulsan 44919, South Korea
| |
Collapse
|
12
|
Tang QY, Ren W, Wang J, Kaneko K. The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database. Mol Biol Evol 2022; 39:msac197. [PMID: 36108094 PMCID: PMC9550990 DOI: 10.1093/molbev/msac197] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.
Collapse
Affiliation(s)
- Qian-Yuan Tang
- Laboratory for Neural Computation and Adaptation, RIKEN Center for Brain Science, 2-1 Hirosawa, Wako, Saitama 351-0106, Japan
| | - Weitong Ren
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Jun Wang
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Kunihiko Kaneko
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Komaba, Meguro, Tokyo 153-8902, Japan
- The Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, Copenhagen 2100-DK, Denmark
| |
Collapse
|
13
|
Dubanevics I, McLeish TCB. Optimising Elastic Network Models for Protein Dynamics and Allostery: Spatial and Modal Cut-offs and Backbone Stiffness. J Mol Biol 2022; 434:167696. [PMID: 35810792 DOI: 10.1016/j.jmb.2022.167696] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 06/16/2022] [Accepted: 06/19/2022] [Indexed: 01/15/2023]
Abstract
The family of coarse-grained models for protein dynamics known as Elastic Network Models (ENMs) require careful choice of parameters to represent well experimental measurements or fully-atomistic simulations. The most basic ENM that represents each protein residue by a node at the position of its C-alpha atom, all connected by springs of equal stiffness, up to a cut-off in distance. Even at this level a choice is required of the optimum cut-off distance and the upper limit of elastic normal modes taken in any sum for physical properties, such as dynamic correlation or allosteric effects on binding. Additionally, backbone-enhanced ENM (BENM) may improve the model by allocating a higher stiffness to springs that connect along the protein backbone. This work reports on the effect of varying these three parameters (distance and mode cutoffs, backbone stiffness) on the dynamical structure of three proteins, Catabolite Activator Protein (CAP), Glutathione S-transferase (GST), and the SARS-CoV-2 Main Protease (M pro ). Our main results are: (1) balancing B-factor and dispersion-relation predictions, a near-universal optimal value of 8.5 Å is advisable for ENMs; (2) inhomogeneity in elasticity brings the first mode containing spatial structure not well-resolved by the ENM typically within the first 20; (3) the BENM only affects modes in the upper third of the distribution, and, additionally to the ENM, is only able to model the dispersion curve better in this vicinity; (4) BENM does not typically affect fluctuation-allostery, which also requires careful treatment of the effector binding to the host protein to capture.
Collapse
|
14
|
Muhammad S, Saba A, Khera RA, Al-Sehemi AG, Algarni H, Iqbal J, Alshahrani MY, Chaudhry AR. Virtual screening of potential inhibitor against breast cancer-causing estrogen receptor alpha (ERα): molecular docking and dynamic simulations. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2072840] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Shabbir Muhammad
- Department of Chemistry, College of Science, King Khalid University, Abha 61413, P.O. Box 9004, Saudi Arabia
| | - Afsheen Saba
- Department of Chemistry, College of Science, University of Agriculture, Faisalabad, Pakistan
| | - Rasheed Ahmad Khera
- Department of Chemistry, College of Science, University of Agriculture, Faisalabad, Pakistan
| | - Abdullah. G. Al-Sehemi
- Department of Chemistry, College of Science, King Khalid University, Abha 61413, P.O. Box 9004, Saudi Arabia
| | - H. Algarni
- Department of Physics, College of Science, King Khalid University, Abha, Saudi Arabia
| | - Javed Iqbal
- Department of Chemistry, College of Science, University of Agriculture, Faisalabad, Pakistan
| | - Mohammad Y. Alshahrani
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, King Khalid University, Abha, Saudi Arabia
| | | |
Collapse
|
15
|
Dong F, Zhang M, Ma R, Lu C, Xu F. Insights of conformational dynamics on catalytic activity in the computational stability design of Bacillus subtilis LipA. Arch Biochem Biophys 2022; 722:109196. [PMID: 35339426 DOI: 10.1016/j.abb.2022.109196] [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: 11/23/2021] [Revised: 03/19/2022] [Accepted: 03/21/2022] [Indexed: 11/18/2022]
Abstract
In protein engineering, the contributions of individual mutations to designed combinatorial mutants are unpredictable. Screening designed mutations that affect enzyme catalytic activity enables evolutions towards efficient activities. Here, Bacillus subtilis LipA (BSLA) was selected as a model protein for thermostabilization designs, and the circular dichroism measurements showed six combinatorial designs with improved stability (from 5.81 °C to 13.61 °C). Based on molecular dynamic simulations, the conformational dynamics of the mutants revealed that mutations alter the populations of conformational states and the increased ensembles of inactive conformations might lead to a reduction in activity. We further demonstrated that the mutations responsible for the reduced enzyme catalytic activity involved a short dynamic correlation path to disturbing the equilibrium conformation of active sites. By removing N82V, which had a close dynamic correlation to the active sites in mutant D3, the redesigned mutant RD3 had an increased activity of 57.6%. By combining computational simulation with experimental verification, this work established that essential sites to counteract the activity-stability trade-off in multipoint combinatorial mutants could be computationally predicted and thus provide a possible strategy by which to indirectly or directly guide protein design.
Collapse
Affiliation(s)
- Fangying Dong
- Ministry of Education Key Laboratory of Carbohydrate Chemistry and Biotechnology, School of Biotechnology, Jiangnan University, 214122, Wuxi, China
| | - Meng Zhang
- Ministry of Education Key Laboratory of Carbohydrate Chemistry and Biotechnology, School of Biotechnology, Jiangnan University, 214122, Wuxi, China
| | - Rui Ma
- Ministry of Education Key Laboratory of Carbohydrate Chemistry and Biotechnology, School of Biotechnology, Jiangnan University, 214122, Wuxi, China
| | - Cheng Lu
- Ministry of Education Key Laboratory of Carbohydrate Chemistry and Biotechnology, School of Biotechnology, Jiangnan University, 214122, Wuxi, China.
| | - Fei Xu
- Ministry of Education Key Laboratory of Carbohydrate Chemistry and Biotechnology, School of Biotechnology, Jiangnan University, 214122, Wuxi, China.
| |
Collapse
|