1
|
Demirtaş K, Erman B, Haliloğlu T. Dynamic correlations: exact and approximate methods for mutual information. Bioinformatics 2024; 40:btae076. [PMID: 38341647 PMCID: PMC10898342 DOI: 10.1093/bioinformatics/btae076] [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: 10/13/2023] [Revised: 01/17/2024] [Accepted: 02/08/2024] [Indexed: 02/12/2024] Open
Abstract
MOTIVATION Proteins are dynamic entities that undergo conformational changes critical for their functions. Understanding the communication pathways and information transfer within proteins is crucial for elucidating allosteric interactions in their mechanisms. This study utilizes mutual information (MI) analysis to probe dynamic allostery. Using two cases, Ubiquitin and PLpro, we have evaluated the accuracy and limitations of different approximations including the exact anisotropic and isotropic models, multivariate Gaussian model, isotropic Gaussian model, and the Gaussian Network Model (GNM) in revealing allosteric interactions. RESULTS Our findings emphasize the required trajectory length for capturing accurate mutual information profiles. Long molecular dynamics trajectories, 1 ms for Ubiquitin and 100 µs for PLpro are used as benchmarks, assuming they represent the ground truth. Trajectory lengths of approximately 5 µs for Ubiquitin and 1 µs for PLpro marked the onset of convergence, while the multivariate Gaussian model accurately captured mutual information with trajectories of 5 ns for Ubiquitin and 350 ns for PLpro. However, the isotropic Gaussian model is less successful in representing the anisotropic nature of protein dynamics, particularly in the case of PLpro, highlighting its limitations. The GNM, however, provides reasonable approximations of long-range information exchange as a minimalist network model based on a single crystal structure. Overall, the optimum trajectory lengths for effective Gaussian approximations of long-time dynamic behavior depend on the inherent dynamics within the protein's topology. The GNM, by showcasing dynamics across relatively diverse time scales, can be used either as a standalone method or to gauge the adequacy of MD simulation lengths. AVAILABILITY AND IMPLEMENTATION Mutual information codes are available at https://github.com/kemaldemirtas/prc-MI.git.
Collapse
Affiliation(s)
- Kemal Demirtaş
- Department of Chemical Engineering, Bogazici University, 34342 Istanbul, Turkey
- Polymer Research Center, Bogazici University, 34342 Istanbul, Turkey
| | - Burak Erman
- Department of Chemical and Biological Engineering, Koc University, 34450 Istanbul, Turkey
| | - Türkan Haliloğlu
- Department of Chemical Engineering, Bogazici University, 34342 Istanbul, Turkey
- Polymer Research Center, Bogazici University, 34342 Istanbul, Turkey
| |
Collapse
|
2
|
Soya N, Xu H, Roldan A, Yang Z, Ye H, Jiang F, Premchandar A, Veit G, Cole SPC, Kappes J, Hegedüs T, Lukacs GL. Folding correctors can restore CFTR posttranslational folding landscape by allosteric domain-domain coupling. Nat Commun 2023; 14:6868. [PMID: 37891162 PMCID: PMC10611759 DOI: 10.1038/s41467-023-42586-8] [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: 11/16/2022] [Accepted: 10/16/2023] [Indexed: 10/29/2023] Open
Abstract
The folding/misfolding and pharmacological rescue of multidomain ATP-binding cassette (ABC) C-subfamily transporters, essential for organismal health, remain incompletely understood. The ABCC transporters core consists of two nucleotide binding domains (NBD1,2) and transmembrane domains (TMD1,2). Using molecular dynamic simulations, biochemical and hydrogen deuterium exchange approaches, we show that the mutational uncoupling or stabilization of NBD1-TMD1/2 interfaces can compromise or facilitate the CFTR(ABCC7)-, MRP1(ABCC1)-, and ABCC6-transporters posttranslational coupled domain-folding in the endoplasmic reticulum. Allosteric or orthosteric binding of VX-809 and/or VX-445 folding correctors to TMD1/2 can rescue kinetically trapped CFTR posttranslational folding intermediates of cystic fibrosis (CF) mutants of NBD1 or TMD1 by global rewiring inter-domain allosteric-networks. We propose that dynamic allosteric domain-domain communications not only regulate ABCC-transporters function but are indispensable to tune the folding landscape of their posttranslational intermediates. These allosteric networks can be compromised by CF-mutations, and reinstated by correctors, offering a framework for mechanistic understanding of ABCC-transporters (mis)folding.
Collapse
Affiliation(s)
- Naoto Soya
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Haijin Xu
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Ariel Roldan
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Zhengrong Yang
- Heersink School of Medicine, University of Alabama School of Medicine, Birmingham, AL, USA
| | - Haoxin Ye
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Fan Jiang
- Heersink School of Medicine, University of Alabama School of Medicine, Birmingham, AL, USA
| | - Aiswarya Premchandar
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Guido Veit
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Susan P C Cole
- Division of Cancer Biology and Genetics, Department of Pathology and Molecular Medicine, Queen's University Cancer Research Institute, Kingston, ON, Canada
| | - John Kappes
- Heersink School of Medicine, University of Alabama School of Medicine, Birmingham, AL, USA
| | - Tamás Hegedüs
- Department of Biophysics and Radiation Biology, Semmelweis University, 1085, Budapest, Hungary
- ELKH-SE Biophysical Virology Research Group, Eötvös Loránd Research Network, Budapest, Hungary
| | - Gergely L Lukacs
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada.
| |
Collapse
|
3
|
Soya N, Xu H, Roldan A, Yang Z, Ye H, Jiang F, Premchandar A, Veit G, Cole SPC, Kappes J, Hegedus T, Lukacs GL. Folding correctors can restore CFTR posttranslational folding landscape by allosteric domain-domain coupling. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.19.563107. [PMID: 37905074 PMCID: PMC10614980 DOI: 10.1101/2023.10.19.563107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The folding/misfolding and pharmacological rescue of multidomain ATP-binding cassette (ABC) C-subfamily transporters, essential for organismal health, remain incompletely understood. The ABCC transporters core consists of two nucleotide binding domains (NBD1,2) and transmembrane domains (TMD1,2). Using molecular dynamic simulations, biochemical and hydrogen deuterium exchange approaches, we show that the mutational uncoupling or stabilization of NBD1-TMD1/2 interfaces can compromise or facilitate the CFTR(ABCC7)-, MRP1(ABCC1)-, and ABCC6-transporters posttranslational coupled domain-folding in the endoplasmic reticulum. Allosteric or orthosteric binding of VX-809 and/or VX-445 folding correctors to TMD1/2 can rescue kinetically trapped CFTR post-translational folding intermediates of cystic fibrosis (CF) mutants of NBD1 or TMD1 by global rewiring inter-domain allosteric-networks. We propose that dynamic allosteric domain-domain communications not only regulate ABCC-transporters function but are indispensable to tune the folding landscape of their post-translational intermediates. These allosteric networks can be compromised by CF-mutations, and reinstated by correctors, offering a framework for mechanistic understanding of ABCC-transporters (mis)folding. One-Sentence Summary Allosteric interdomain communication and its modulation are critical determinants of ABCC-transporters post-translational conformational biogenesis, misfolding, and pharmacological rescue.
Collapse
|
4
|
Cookis T, Mattos C. Crystal Structure Reveals the Full Ras-Raf Interface and Advances Mechanistic Understanding of Raf Activation. Biomolecules 2021; 11:996. [PMID: 34356620 PMCID: PMC8301913 DOI: 10.3390/biom11070996] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 07/01/2021] [Accepted: 07/02/2021] [Indexed: 01/17/2023] Open
Abstract
Ras and Raf-kinase interact through the Ras-binding (RBD) and cysteine-rich domains (CRD) of Raf to signal through the mitogen-activated protein kinase pathway, yet the molecular mechanism leading to Raf activation has remained elusive. We present the 2.8 Å crystal structure of the HRas-CRaf-RBD_CRD complex showing the Ras-Raf interface as a continuous surface on Ras, as seen in the KRas-CRaf-RBD_CRD structure. In molecular dynamics simulations of a Ras dimer model formed through the α4-α5 interface, the CRD is dynamic and located between the two Ras protomers, poised for direct or allosteric modulation of functionally relevant regions of Ras and Raf. We propose a molecular model in which Ras binding is involved in the release of Raf autoinhibition while the Ras-Raf complex dimerizes to promote a platform for signal amplification, with Raf-CRD centrally located to impact regulation and function.
Collapse
Affiliation(s)
| | - Carla Mattos
- Department of Chemistry and Chemical Biology, Northeastern University, Boston, MA 02115, USA;
| |
Collapse
|
5
|
Li Y, Guan X, Chaffey PK, Ruan Y, Ma B, Shang S, Himmel ME, Beckham GT, Long H, Tan Z. Carbohydrate-binding module O-mannosylation alters binding selectivity to cellulose and lignin. Chem Sci 2020; 11:9262-9271. [PMID: 34123172 PMCID: PMC8163390 DOI: 10.1039/d0sc01812k] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Improved understanding of the effect of protein glycosylation is expected to provide the foundation for the design of protein glycoengineering strategies. In this study, we examine the impact of O-glycosylation on the binding selectivity of a model Family 1 carbohydrate-binding module (CBM), which has been shown to be one of the primary sub-domains responsible for non-productive lignin binding in multi-modular cellulases. Specifically, we examine the relationship between glycan structure and the binding specificity of the CBM to cellulose and lignin substrates. We find that the glycosylation pattern of the CBM exhibits a strong influence on the binding affinity and the selectivity between both cellulose and lignin. In addition, the large set of binding data collected allows us to examine the relationship between binding affinity and the correlation in motion between pairs of glycosylation sites. Our results suggest that glycoforms displaying highly correlated motion in their glycosylation sites tend to bind cellulose with high affinity and lignin with low affinity. Taken together, this work helps lay the groundwork for future exploitation of glycoengineering as a tool to improve the performance of industrial enzymes. Improved understanding of the effect of protein glycosylation is expected to provide the foundation for the design of protein glycoengineering strategies.![]()
Collapse
Affiliation(s)
- Yaohao Li
- Institute of Materia Medica, Chinese Academy of Medical Sciences, Peking Union Medical College Beijing 100050 China .,Department of Chemistry and Biochemistry and BioFrontiers Institute, University of Colorado Boulder CO 80303 USA
| | - Xiaoyang Guan
- Department of Chemistry and Biochemistry and BioFrontiers Institute, University of Colorado Boulder CO 80303 USA
| | - Patrick K Chaffey
- Department of Chemistry and Biochemistry and BioFrontiers Institute, University of Colorado Boulder CO 80303 USA
| | - Yuan Ruan
- Department of Chemistry and Biochemistry and BioFrontiers Institute, University of Colorado Boulder CO 80303 USA
| | - Bo Ma
- Institute of Materia Medica, Chinese Academy of Medical Sciences, Peking Union Medical College Beijing 100050 China
| | - Shiying Shang
- School of Pharmaceutical Sciences, Tsinghua University Beijing 100084 China
| | - Michael E Himmel
- Biosciences Center, National Renewable Energy Laboratory Golden CO 80401 USA
| | - Gregg T Beckham
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory Golden CO 80401 USA
| | - Hai Long
- Computational Science Center, National Renewable Energy Laboratory Golden CO 80401 USA
| | - Zhongping Tan
- Institute of Materia Medica, Chinese Academy of Medical Sciences, Peking Union Medical College Beijing 100050 China
| |
Collapse
|
6
|
Sajeevan KA, Roy D. Simulation of differential structure and dynamics of disulfide bond isoforms of conopeptide
AuIB
in presence of human
α
3
β
4
nAChR. Pept Sci (Hoboken) 2020. [DOI: 10.1002/pep2.24183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Karuna Anna Sajeevan
- Department of Chemistry Birla Institute of Technology and Science‐Pilani Hyderabad Telangana India
| | - Durba Roy
- Department of Chemistry Birla Institute of Technology and Science‐Pilani Hyderabad Telangana India
| |
Collapse
|
7
|
Lambrughi M, Sanader Maršić Ž, Saez-Jimenez V, Mapelli V, Olsson L, Papaleo E. Conformational gating in ammonia lyases. Biochim Biophys Acta Gen Subj 2020; 1864:129605. [PMID: 32222547 DOI: 10.1016/j.bbagen.2020.129605] [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: 09/09/2019] [Revised: 03/11/2020] [Accepted: 03/23/2020] [Indexed: 11/17/2022]
Abstract
BACKGROUND Ammonia lyases are enzymes of industrial and biomedical interest. Knowledge of structure-dynamics-function relationship in ammonia lyases is instrumental for exploiting the potential of these enzymes in industrial or biomedical applications. METHODS We investigated the conformational changes in the proximity of the catalytic pocket of a 3-methylaspartate ammonia lyase (MAL) as a model system. At this scope, we used microsecond all-atom molecular dynamics simulations, analyzed with dimensionality reduction techniques, as well as in terms of contact networks and correlated motions. RESULTS We identify two regulatory elements in the MAL structure, i.e., the β5-α2 loop and the helix-hairpin-loop subdomain. These regulatory elements undergo conformational changes switching from 'occluded' to 'open' states. The rearrangements are coupled to changes in the accessibility of the active site. The β5-α2 loop and the helix-hairpin-loop subdomain modulate the formation of tunnels from the protein surface to the catalytic site, making the active site more accessible to the substrate when they are in an open state. CONCLUSIONS Our work pinpoints a sequential mechanism, in which the helix-hairpin-loop subdomain of MAL needs to break a subset of intramolecular interactions first to favor the displacement of the β5-α2 loop. The coupled conformational changes of these two elements contribute to modulate the accessibility of the catalytic site. GENERAL SIGNIFICANCE Similar molecular mechanisms can have broad relevance in other ammonia lyases with similar regulatory loops. Our results also imply that it is important to account for protein dynamics in the design of variants of ammonia lyases for industrial and biomedical applications.
Collapse
Affiliation(s)
- Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Željka Sanader Maršić
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Veronica Saez-Jimenez
- Division of Industrial Biotechnology, Department of Biology and Biological Engineering, Chalmers University of Technology, Göteborg, Sweden
| | - Valeria Mapelli
- Division of Industrial Biotechnology, Department of Biology and Biological Engineering, Chalmers University of Technology, Göteborg, Sweden
| | - Lisbeth Olsson
- Division of Industrial Biotechnology, Department of Biology and Biological Engineering, Chalmers University of Technology, Göteborg, Sweden
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark; Translational Disease Systems Biology, Faculty of Health and Medical Sciences, Novo Nordisk Foundation Center for Protein Research University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
8
|
Izmailov SA, Rabdano SO, Hasanbasri Z, Podkorytov IS, Saxena S, Skrynnikov NR. Structural and dynamic origins of ESR lineshapes in spin-labeled GB1 domain: the insights from spin dynamics simulations based on long MD trajectories. Sci Rep 2020; 10:957. [PMID: 31969574 PMCID: PMC6976580 DOI: 10.1038/s41598-019-56750-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2019] [Accepted: 12/06/2019] [Indexed: 12/04/2022] Open
Abstract
Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 μs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1 and χ2 are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.
Collapse
Affiliation(s)
- Sergei A Izmailov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, 199034, Russia
| | - Sevastyan O Rabdano
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, 199034, Russia
| | - Zikri Hasanbasri
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA, 15260, USA
| | - Ivan S Podkorytov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, 199034, Russia
| | - Sunil Saxena
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA, 15260, USA
| | - Nikolai R Skrynnikov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, 199034, Russia.
- Department of Chemistry, Purdue University, West Lafayette, IN, 47907, USA.
| |
Collapse
|
9
|
Ahuja LG, Aoto PC, Kornev AP, Veglia G, Taylor SS. Dynamic allostery-based molecular workings of kinase:peptide complexes. Proc Natl Acad Sci U S A 2019; 116:15052-15061. [PMID: 31285328 PMCID: PMC6660753 DOI: 10.1073/pnas.1900163116] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
A dense interplay between structure and dynamics underlies the working of proteins, especially enzymes. Protein kinases are molecular switches that are optimized for their regulation rather than catalytic turnover rates. Using long-simulations dynamic allostery analysis, this study describes an exploration of the dynamic kinase:peptide complex. We have used protein kinase A (PKA) as a model system as a generic prototype of the protein kinase superfamily of signaling enzymes. Our results explain the role of dynamic coupling of active-site residues that must work in coherence to provide for a successful activation or inhibition response from the kinase. Amino acid networks-based community analysis allows us to ponder the conformational entropy of the kinase:nucleotide:peptide ternary complex. We use a combination of 7 peptides that include 3 types of PKA-binding partners: Substrates, products, and inhibitors. The substrate peptides provide for dynamic insights into the enzyme:substrate complex, while the product phospho-peptide allows for accessing modes of enzyme:product release. Mapping of allosteric communities onto the PKA structure allows us to locate the more unvarying and flexible dynamic regions of the kinase. These distributions, when correlated with the structural elements of the kinase core, allow for a detailed exploration of key dynamics-based signatures that could affect peptide recognition and binding at the kinase active site. These studies provide a unique dynamic allostery-based perspective to kinase:peptide complexes that have previously been explored only in a structural or thermodynamic context.
Collapse
Affiliation(s)
- Lalima G Ahuja
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093;
| | - Phillip C Aoto
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093
| | - Alexandr P Kornev
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093
| | - Gianluigi Veglia
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455
- Department of Biochemistry, Molecular Biology, and Biophysics, University of Minnesota, Minneapolis, MN 55455
| | - Susan S Taylor
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093;
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093
| |
Collapse
|
10
|
Sajeevan KA, Roy D. Principal Component Analysis of a Conotoxin Delineates the Link among Peptide Sequence, Dynamics, and Disulfide Bond Isoforms. J Phys Chem B 2019; 123:5483-5493. [DOI: 10.1021/acs.jpcb.9b04090] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Karuna Anna Sajeevan
- Department of Chemistry, Birla Institute of Technology and Science-Pilani, Hyderabad Campus, Jawahar Nagar,
Kapra Mandal, Hyderabad, Telangana 500078, India
| | - Durba Roy
- Department of Chemistry, Birla Institute of Technology and Science-Pilani, Hyderabad Campus, Jawahar Nagar,
Kapra Mandal, Hyderabad, Telangana 500078, India
| |
Collapse
|
11
|
Lambrughi M, Tiberti M, Allega MF, Sora V, Nygaard M, Toth A, Salamanca Viloria J, Bignon E, Papaleo E. Analyzing Biomolecular Ensembles. Methods Mol Biol 2019; 2022:415-451. [PMID: 31396914 DOI: 10.1007/978-1-4939-9608-7_18] [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: 06/10/2023]
Abstract
Several techniques are available to generate conformational ensembles of proteins and other biomolecules either experimentally or computationally. These methods produce a large amount of data that need to be analyzed to identify structure-dynamics-function relationship. In this chapter, we will cover different tools to unveil the information hidden in conformational ensemble data and to guide toward the rationalization of the data. We included routinely used approaches such as dimensionality reduction, as well as new methods inspired by high-order statistics and graph theory.
Collapse
Affiliation(s)
- Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Matteo Tiberti
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Maria Francesca Allega
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Valentina Sora
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Mads Nygaard
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Agota Toth
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Juan Salamanca Viloria
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Emmanuelle Bignon
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark.
| |
Collapse
|
12
|
Long S, Wang J, Tian P. Significance of triple torsional correlations in proteins. RSC Adv 2019; 9:13949-13958. [PMID: 35519605 PMCID: PMC9064167 DOI: 10.1039/c9ra02191d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Accepted: 04/21/2019] [Indexed: 11/21/2022] Open
Abstract
The free energy landscape (FEL) of a given complex molecular system is fundamentally the joint probability density of its many comprising degrees of freedom (DOFs). Computation of a complete FEL at atomistic scale is unfortunately intractable for a typical biomolecular system. The challenge of entropy calculation comes from various correlations among different DOFs. The common strategy to treat such complexity is expansion of the full correlation into various orders of local correlations. In reality, expansion is usually cut off at the second order (i.e. pairwise interactions) for protein torsional correlations without reliable estimation of the resulting error. Here, we estimated the mutual information of different torsion sets and found that triple correlations were significant for both local/distant residue pairs and consecutive backbone torsional segments. As expected, the third order approximations were found to be consistently better than the second order approximations. These findings were true for all analyzed proteins with different folds, were independent of the two different force fields utilized to generate trajectory sets, and were therefore likely to be of general importance for proteins. Additionally, binning strategies are of universal importance for numerical computation of correlations, we here provided a detailed comparison between equal-width and equal-sample binning for different bin numbers and demonstrated the impact of binning strategies on variances and biases of calculated mutual information. Our observation suggested that caution should be taken when quantitative comparison of correlations were intended between different studies with different binning strategies. Torsional mutual information for 10 typical residue pairs calculated with full joint distributions (MI), second order expansion (MI2), third order expansions (MI3), and their linear recombinations (MILR).![]()
Collapse
Affiliation(s)
| | | | - Pu Tian
- School of Life Science
- School Artificial Intelligence
- Jilin University
- Changchun
- China 130012
| |
Collapse
|
13
|
Tekpinar M, Yildirim A. Only a Subset of Normal Modes is Sufficient to Identify Linear Correlations in Proteins. J Chem Inf Model 2018; 58:1947-1961. [DOI: 10.1021/acs.jcim.8b00486] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
| | - Ahmet Yildirim
- Department of Physics, Siirt University, 56100 Siirt, Turkey
| |
Collapse
|
14
|
Hashem S, Tiberti M, Fornili A. Allosteric modulation of cardiac myosin dynamics by omecamtiv mecarbil. PLoS Comput Biol 2017; 13:e1005826. [PMID: 29108014 PMCID: PMC5690683 DOI: 10.1371/journal.pcbi.1005826] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2017] [Revised: 11/16/2017] [Accepted: 10/16/2017] [Indexed: 01/10/2023] Open
Abstract
New promising avenues for the pharmacological treatment of skeletal and heart muscle diseases rely on direct sarcomeric modulators, which are molecules that can directly bind to sarcomeric proteins and either inhibit or enhance their activity. A recent breakthrough has been the discovery of the myosin activator omecamtiv mecarbil (OM), which has been shown to increase the power output of the cardiac muscle and is currently in clinical trials for the treatment of heart failure. While the overall effect of OM on the mechano-chemical cycle of myosin is to increase the fraction of myosin molecules in the sarcomere that are strongly bound to actin, the molecular basis of its action is still not completely clear. We present here a Molecular Dynamics study of the motor domain of human cardiac myosin bound to OM, where the effects of the drug on the dynamical properties of the protein are investigated for the first time with atomistic resolution. We found that OM has a double effect on myosin dynamics, inducing a) an increased coupling of the motions of the converter and lever arm subdomains to the rest of the protein and b) a rewiring of the network of dynamic correlations, which produces preferential communication pathways between the OM binding site and distant functional regions. The location of the residues responsible for these effects suggests possible strategies for the future development of improved drugs and the targeting of specific cardiomyopathy-related mutations.
Collapse
Affiliation(s)
- Shaima Hashem
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
| | - Matteo Tiberti
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
| | - Arianna Fornili
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
- The Thomas Young Centre for Theory and Simulation of Materials, London, United Kingdom
- * E-mail:
| |
Collapse
|
15
|
Salamanca Viloria J, Allega MF, Lambrughi M, Papaleo E. An optimal distance cutoff for contact-based Protein Structure Networks using side-chain centers of mass. Sci Rep 2017; 7:2838. [PMID: 28588190 PMCID: PMC5460117 DOI: 10.1038/s41598-017-01498-6] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Accepted: 03/28/2017] [Indexed: 02/05/2023] Open
Abstract
Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.
Collapse
Affiliation(s)
- Juan Salamanca Viloria
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Maria Francesca Allega
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100, Copenhagen, Denmark.
| |
Collapse
|
16
|
Nygaard M, Terkelsen T, Vidas Olsen A, Sora V, Salamanca Viloria J, Rizza F, Bergstrand-Poulsen S, Di Marco M, Vistesen M, Tiberti M, Lambrughi M, Jäättelä M, Kallunki T, Papaleo E. The Mutational Landscape of the Oncogenic MZF1 SCAN Domain in Cancer. Front Mol Biosci 2016; 3:78. [PMID: 28018905 PMCID: PMC5156680 DOI: 10.3389/fmolb.2016.00078] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Accepted: 11/17/2016] [Indexed: 11/24/2022] Open
Abstract
SCAN domains in zinc-finger transcription factors are crucial mediators of protein-protein interactions. Up to 240 SCAN-domain encoding genes have been identified throughout the human genome. These include cancer-related genes, such as the myeloid zinc finger 1 (MZF1), an oncogenic transcription factor involved in the progression of many solid cancers. The mechanisms by which SCAN homo- and heterodimers assemble and how they alter the transcriptional activity of zinc-finger transcription factors in cancer and other diseases remain to be investigated. Here, we provide the first description of the conformational ensemble of the MZF1 SCAN domain cross-validated against NMR experimental data, which are probes of structure and dynamics on different timescales. We investigated the protein-protein interaction network of MZF1 and how it is perturbed in different cancer types by the analyses of high-throughput proteomics and RNASeq data. Collectively, we integrated many computational approaches, ranging from simple empirical energy functions to all-atom microsecond molecular dynamics simulations and network analyses to unravel the effects of cancer-related substitutions in relation to MZF1 structure and interactions.
Collapse
Affiliation(s)
- Mads Nygaard
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Thilde Terkelsen
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - André Vidas Olsen
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Valentina Sora
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Juan Salamanca Viloria
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Fabio Rizza
- Department of Biomedical Sciences, University of Padua Padua, Italy
| | - Sanne Bergstrand-Poulsen
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Miriam Di Marco
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Mette Vistesen
- Cell Stress and Survival Unit and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Matteo Tiberti
- Department of Chemistry and Biochemistry, School of Biological and Chemical Sciences, Queen Mary University of London London, UK
| | - Matteo Lambrughi
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Marja Jäättelä
- Unit of Cell Death and Metabolism and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Tuula Kallunki
- Unit of Cell Death and Metabolism and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory and Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center Copenhagen, Denmark
| |
Collapse
|
17
|
Calligari P, Gerolin M, Abergel D, Polimeno A. Decomposition of Proteins into Dynamic Units from Atomic Cross-Correlation Functions. J Chem Theory Comput 2016; 13:309-319. [DOI: 10.1021/acs.jctc.6b00702] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Paolo Calligari
- Dipartimento
di Scienze Chimiche, Università di Padova, via Marzolo, 1, I-35131 Padova, Italy
| | - Marco Gerolin
- Dipartimento
di Scienze Chimiche, Università di Padova, via Marzolo, 1, I-35131 Padova, Italy
- Département
de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules (LBM), 24 rue Lhomond, 75005 Paris, France
| | - Daniel Abergel
- Département
de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules (LBM), 24 rue Lhomond, 75005 Paris, France
| | - Antonino Polimeno
- Dipartimento
di Scienze Chimiche, Università di Padova, via Marzolo, 1, I-35131 Padova, Italy
| |
Collapse
|