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
|
Chikhale RV, Abdelghani HTM, Deka H, Pawar AD, Patil PC, Bhowmick S. Machine learning assisted methods for the identification of low toxicity inhibitors of Enoyl-Acyl Carrier Protein Reductase (InhA). Comput Biol Chem 2024; 110:108034. [PMID: 38430612 DOI: 10.1016/j.compbiolchem.2024.108034] [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/30/2023] [Revised: 01/20/2024] [Accepted: 02/18/2024] [Indexed: 03/05/2024]
Abstract
Tuberculosis (TB) is one of the life-threatening infectious diseases with prehistoric origins and occurs in almost all habitable parts of the world. TB mainly affects the lungs, and its etiological agent is Mycobacterium tuberculosis (Mtb). In 2022, more than 10 million people were infected worldwide, and 1.3 million were children. The current study considered the in-silico and machine learning (ML) approaches to explore the potential anti-TB molecules from the SelleckChem database against Enoyl-Acyl Carrier Protein Reductase (InhA). Initially, the entire database of ∼ 119000 molecules was sorted out through drug-likeness. Further, the molecular docking study was conducted to reduce the chemical space. The standard TB drug molecule's binding energy was considered a threshold, and molecules found with lower affinity were removed for further analyses. Finally, the molecules were checked for the pharmacokinetic and toxicity studies, and compounds found to have acceptable pharmacokinetic parameters and were non-toxic were considered as final promising molecules for InhA. The above approach further evaluated five molecules for ML-based toxicity and synthetic accessibility assessment. Not a single molecule was found toxic and each of them was revealed as easy to synthesise. The complex between InhA and proposed and standard molecules was considered for molecular dynamics simulation. Several statistical parameters showed the stability between InhA and the proposed molecule. The high binding affinity was also found for each of the molecules towards InhA using the MM-GBSA approach. Hence, the above approaches and findings exposed the potentiality of the proposed molecules against InhA.
Collapse
Affiliation(s)
- Rupesh V Chikhale
- Department of Pharmaceutical and Biological Chemistry, School of Pharmacy, University College London, London, UK
| | - Heba Taha M Abdelghani
- Department of Exercise Physiology, College of Sport Sciences and Physical Activity, King Saud University, Riyadh 11451, Saudi Arabia
| | - Hemchandra Deka
- SilicoScientia Private Limited, Nagananda Commercial Complex, No. 07/3, 15/1, 18th Main Road, Jayanagar 9th Block, Bengaluru 560041, India
| | - Atul Darasing Pawar
- SilicoScientia Private Limited, Nagananda Commercial Complex, No. 07/3, 15/1, 18th Main Road, Jayanagar 9th Block, Bengaluru 560041, India
| | - Pritee Chunarkar Patil
- Department of Bioinformatics, Rajiv Gandhi Institute of IT and Biotechnology, Bharati Vidyapeeth (Deemed to be University), Pune-Satara Road, Pune, India
| | - Shovonlal Bhowmick
- SilicoScientia Private Limited, Nagananda Commercial Complex, No. 07/3, 15/1, 18th Main Road, Jayanagar 9th Block, Bengaluru 560041, India.
| |
Collapse
|
3
|
Trueba-Gómez R, Rosenfeld-Mann F, Estrada-Juárez H. Prediction of the antigenic regions in eight RhD variants identified by computational biology. Vox Sang 2024; 119:590-597. [PMID: 38523363 DOI: 10.1111/vox.13620] [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: 12/14/2023] [Revised: 02/23/2024] [Accepted: 03/08/2024] [Indexed: 03/26/2024]
Abstract
BACKGROUND AND OBJECTIVES Changes in RHD generate variations in protein structure that lead to antigenic variants. The classical model divides them into quantitative (weak and Del) and qualitative (partial D). There are two types of protein antigens: linear and conformational. Computational biology analyses the theoretical assembly of tertiary protein structures and allows us to identify the 'topological' differences between isoforms. Our aim was to determine the theoretical antigenic differences between weak RhD variants compared with normal RhD based on structural analysis using bioinformatic techniques. MATERIALS AND METHODS We analysed the variations in secondary structures and hydrophobicity of RHD*01, RHD*01W.1, W2, W3, RHD*09.03.01, RHD*09.04, RHD*11, RHD*15 and RHD*21. We then modelled the tertiary structure and calculated their probable antigenic regions, intra-protein interactions, displacement and membrane width and compared them with Rhce. RESULTS The 10 proteins are similar in their secondary structure and hydrophobicity, with the main differences observed in the exofacial coils. We identified six potential antigenic regions: one that is unique to RhD (R3), one that is common to all D (R6), three that are highly variable among RhD isoforms (R1, R2 and R4), one that they share with Rhce (R5) and two that are unique to Rhce (Ra and Rbc). CONCLUSION The alloimmunization capacity of these subjects could be explained by the variability of the antigen pattern, which is not necessarily recognized or recognized with lower intensity by the commercially available antibodies, and not because they have a lower protein concentration in the membrane.
Collapse
Affiliation(s)
- Rocio Trueba-Gómez
- Instituto Nacional de Perinatología "Isidro Espinosa de los Reyes," Coordinación de Hematología Perinatal, Mexico City, Mexico
- Comité de Trombosis y Hemostasia AMEH-CLAHT, A.C., Mexico City, Mexico
| | | | - Higinio Estrada-Juárez
- Instituto Nacional de Perinatología "Isidro Espinosa de los Reyes," Coordinación de Hematología Perinatal, Mexico City, Mexico
- Comité de Trombosis y Hemostasia AMEH-CLAHT, A.C., Mexico City, Mexico
| |
Collapse
|
4
|
Chitluri KK, Emerson IA. The importance of protein domain mutations in cancer therapy. Heliyon 2024; 10:e27655. [PMID: 38509890 PMCID: PMC10950675 DOI: 10.1016/j.heliyon.2024.e27655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Revised: 02/28/2024] [Accepted: 03/05/2024] [Indexed: 03/22/2024] Open
Abstract
Cancer is a complex disease that is caused by multiple genetic factors. Researchers have been studying protein domain mutations to understand how they affect the progression and treatment of cancer. These mutations can significantly impact the development and spread of cancer by changing the protein structure, function, and signalling pathways. As a result, there is a growing interest in how these mutations can be used as prognostic indicators for cancer prognosis. Recent studies have shown that protein domain mutations can provide valuable information about the severity of the disease and the patient's response to treatment. They may also be used to predict the response and resistance to targeted therapy in cancer treatment. The clinical implications of protein domain mutations in cancer are significant, and they are regarded as essential biomarkers in oncology. However, additional techniques and approaches are required to characterize changes in protein domains and predict their functional effects. Machine learning and other computational tools offer promising solutions to this challenge, enabling the prediction of the impact of mutations on protein structure and function. Such predictions can aid in the clinical interpretation of genetic information. Furthermore, the development of genome editing tools like CRISPR/Cas9 has made it possible to validate the functional significance of mutants more efficiently and accurately. In conclusion, protein domain mutations hold great promise as prognostic and predictive biomarkers in cancer. Overall, considerable research is still needed to better define genetic and molecular heterogeneity and to resolve the challenges that remain, so that their full potential can be realized.
Collapse
Affiliation(s)
- Kiran Kumar Chitluri
- Bioinformatics Programming Lab, Department of Bio-Sciences, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, TN, 632014, India
| | - Isaac Arnold Emerson
- Bioinformatics Programming Lab, Department of Bio-Sciences, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, TN, 632014, India
| |
Collapse
|
5
|
Hait S, Kundu S. Revisiting structural organization of proteins at high temperature from a network perspective. Comput Biol Chem 2024; 108:107978. [PMID: 37956471 DOI: 10.1016/j.compbiolchem.2023.107978] [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: 08/01/2023] [Revised: 10/08/2023] [Accepted: 10/29/2023] [Indexed: 11/15/2023]
Abstract
Interactions between distantly placed amino acids in the primary chain (long-range) play a very crucial role in the formation and stabilization of the tertiary structure of a protein, while interactions between closely placed amino acids in the primary chain (short-range) mostly stabilize the secondary structures. Every protein needs to maintain marginal stability in order to perform its physiological functions in its native environment. The requirements for this stability in mesophilic and thermophilic proteins are different. Thermophilic proteins need to form more interactions as well as more stable interactions to survive in the extreme environment, they live in. Here, we aim to find out how the interacting amino acids in three-dimensional space are positioned in the primary chains in thermophilic and mesophilic. How does this arrangement help thermophiles to maintain their structural integrity at high temperatures? Working on a dataset of 1560 orthologous pairs we perceive that thermophiles are not only enriched with long-range interactions, they feature bigger connected clusters and higher network densities compared to their mesophilic orthologs, at higher interaction strengths between the amino acids. Moreover, we have observed the enrichment of different types of interactions at different secondary structural regions.
Collapse
Affiliation(s)
- Suman Hait
- Department of Biophysics, Molecular Biology and Bioinformatics, 92, Acharya Prafulla Chandra Road, Kolkata 700009, India.
| | - Sudip Kundu
- Department of Biophysics, Molecular Biology and Bioinformatics, 92, Acharya Prafulla Chandra Road, Kolkata 700009, India.
| |
Collapse
|
6
|
Gong J, Jiang L, Chen Y, Zhang Y, Li X, Ma Z, Fu Z, He F, Sun P, Ren Z, Tian M. THPLM: a sequence-based deep learning framework for protein stability changes prediction upon point variations using pretrained protein language model. Bioinformatics 2023; 39:btad646. [PMID: 37874953 PMCID: PMC10627365 DOI: 10.1093/bioinformatics/btad646] [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/29/2023] [Revised: 09/25/2023] [Accepted: 10/22/2023] [Indexed: 10/26/2023] Open
Abstract
MOTIVATION Quantitative determination of protein thermodynamic stability is a critical step in protein and drug design. Reliable prediction of protein stability changes caused by point variations contributes to developing-related fields. Over the past decades, dozens of structure-based and sequence-based methods have been proposed, showing good prediction performance. Despite the impressive progress, it is necessary to explore wild-type and variant protein representations to address the problem of how to represent the protein stability change in view of global sequence. With the development of structure prediction using learning-based methods, protein language models (PLMs) have shown accurate and high-quality predictions of protein structure. Because PLM captures the atomic-level structural information, it can help to understand how single-point variations cause functional changes. RESULTS Here, we proposed THPLM, a sequence-based deep learning model for stability change prediction using Meta's ESM-2. With ESM-2 and a simple convolutional neural network, THPLM achieved comparable or even better performance than most methods, including sequence-based and structure-based methods. Furthermore, the experimental results indicate that the PLM's ability to generate representations of sequence can effectively improve the ability of protein function prediction. AVAILABILITY AND IMPLEMENTATION The source code of THPLM and the testing data can be accessible through the following links: https://github.com/FPPGroup/THPLM.
Collapse
Affiliation(s)
- Jianting Gong
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Lili Jiang
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Yongbing Chen
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Yixiang Zhang
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Xue Li
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Zhiqiang Ma
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Department of Computer Science, College of Humanities and Sciences of Northeast Normal University, Changchun 130117, China
| | - Zhiguo Fu
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
| | - Fei He
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
| | - Pingping Sun
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
| | - Zilin Ren
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Mingyao Tian
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| |
Collapse
|
7
|
Ramakrishna Reddy P, Kulandaisamy A, Michael Gromiha M. TMH Stab-pred: Predicting the stability of α-helical membrane proteins using sequence and structural features. Methods 2023; 218:118-124. [PMID: 37572768 DOI: 10.1016/j.ymeth.2023.08.005] [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: 04/08/2023] [Revised: 08/02/2023] [Accepted: 08/04/2023] [Indexed: 08/14/2023] Open
Abstract
The folding and stability of transmembrane proteins (TMPs) are governed by the insertion of secondary structural elements into the cell membrane followed by their assembly. Understanding the important features that dictate the stability of TMPs is important for elucidating their functions. In this work, we related sequence and structure-based parameters with free energy (ΔG0) of α-helical membrane proteins. Our results showed that the free energy transfer of hydrophobic peptides, relative contact order, total interaction energy, number of hydrogen bonds and lipid accessibility of transmembrane regions are important for stability. Further, we have developed multiple-regression models to predict the stability of α-helical membrane proteins using these features and our method can predict the stability with a correlation and mean absolute error (MAE) of 0.89 and 1.21 kcal/mol, respectively, on jack-knife test. The method was validated with a blind test set of three recently reported experimental ΔG0, which could predict the stability within an average MAE of 0.51 kcal/mol. Further, we developed a webserver for predicting the stability and it is freely available at (https://web.iitm.ac.in/bioinfo2/TMHS/). The importance of selected parameters and limitations are discussed.
Collapse
Affiliation(s)
- P Ramakrishna Reddy
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, Tamil Nadu, India
| | - A Kulandaisamy
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, Tamil Nadu, India; Basic and Translational Research Division, Department of Cardiology, Boston Children's Hospital, Boston, MA 02115, USA
| | - M Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, Tamil Nadu, India; Department of Computer Science, Tokyo Institute of Technology, Yokohama, Japan; Department of Computer Science, National University of Singapore, Singapore.
| |
Collapse
|
8
|
Han J, Jiang S, Zhou Z, Lin M, Wang J. Artificial Proteins Designed from G3LEA Contribute to Enhancement of Oxidation Tolerance in E. coli in a Chaperone-like Manner. Antioxidants (Basel) 2023; 12:1147. [PMID: 37371877 DOI: 10.3390/antiox12061147] [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: 03/07/2023] [Revised: 04/12/2023] [Accepted: 05/22/2023] [Indexed: 06/29/2023] Open
Abstract
G3LEA is a family of proteins that exhibit chaperone-like activity when under distinct stress. In previous research, DosH was identified as a G3LEA protein from model extremophile-Deinococcus radiodurans R1 with a crucial core HD domain consisting of eight 11-mer motifs. However, the roles of motifs participating in the process of resistance to stress and their underlying mechanisms remain unclear. Here, eight different proteins with tandem repeats of the same motif were synthesized, named Motif1-8, respectively, whose function and structure were discussed. In this way, the role of each motif in the HD domain can be comprehensively analyzed, which can help in finding possibly crucial amino acid sites. Circular dichroism results showed that all proteins were intrinsically ordered in phosphate buffer, and changed into more α-helical ordered structures with the addition of trifluoroethanol and glycerol. Transformants expressing artificial proteins had significantly higher stress resistance to oxidation, desiccation, salinity and freezing compared with the control group; E. coli with Motif1 and Motif8 had more outstanding performance in particular. Moreover, enzymes and membrane protein protection viability suggested that Motif1 and Motif8 had more positive influences on various molecules, demonstrating a protective role in a chaperone-like manner. Based on these results, the artificial proteins synthesized according to the rule of 11-mer motifs have a similar function to wildtype protein. Regarding the sequence in all motifs, there are more amino acids to produce H bonds and α-helices, and more amino acids to promote interaction between proteins in Motif1 and Motif8; in addition, considering linkers, there are possibly more amino acids forming α-helix and binding substrates in these two proteins, which potentially provides some ideas for us to design potential ideal stress-response elements for synthetic biology. Therefore, the amino acid composition of the 11-mer motif and linker is likely responsible for its biological function.
Collapse
Affiliation(s)
- Jiahui Han
- Key Laboratory of Agricultural Microbiome (MARA), Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Shijie Jiang
- School of Life Science and Engineering, Southwest University of Science and Technology, Mianyang 621010, China
| | - Zhengfu Zhou
- Key Laboratory of Agricultural Microbiome (MARA), Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Min Lin
- Key Laboratory of Agricultural Microbiome (MARA), Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Jin Wang
- Key Laboratory of Agricultural Microbiome (MARA), Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| |
Collapse
|
9
|
PSnpBind-ML: predicting the effect of binding site mutations on protein-ligand binding affinity. J Cheminform 2023; 15:31. [PMID: 36864534 PMCID: PMC9983232 DOI: 10.1186/s13321-023-00701-3] [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: 10/21/2022] [Accepted: 02/17/2023] [Indexed: 03/04/2023] Open
Abstract
Protein mutations, especially those which occur in the binding site, play an important role in inter-individual drug response and may alter binding affinity and thus impact the drug's efficacy and side effects. Unfortunately, large-scale experimental screening of ligand-binding against protein variants is still time-consuming and expensive. Alternatively, in silico approaches can play a role in guiding those experiments. Methods ranging from computationally cheaper machine learning (ML) to the more expensive molecular dynamics have been applied to accurately predict the mutation effects. However, these effects have been mostly studied on limited and small datasets, while ideally a large dataset of binding affinity changes due to binding site mutations is needed. In this work, we used the PSnpBind database with six hundred thousand docking experiments to train a machine learning model predicting protein-ligand binding affinity for both wild-type proteins and their variants with a single-point mutation in the binding site. A numerical representation of the protein, binding site, mutation, and ligand information was encoded using 256 features, half of them were manually selected based on domain knowledge. A machine learning approach composed of two regression models is proposed, the first predicting wild-type protein-ligand binding affinity while the second predicting the mutated protein-ligand binding affinity. The best performing models reported an RMSE value within 0.5 [Formula: see text] 0.6 kcal/mol-1 on an independent test set with an R2 value of 0.87 [Formula: see text] 0.90. We report an improvement in the prediction performance compared to several reported models developed for protein-ligand binding affinity prediction. The obtained models can be used as a complementary method in early-stage drug discovery. They can be applied to rapidly obtain a better overview of the ligand binding affinity changes across protein variants carried by people in the population and narrow down the search space where more time-demanding methods can be used to identify potential leads that achieve a better affinity for all protein variants.
Collapse
|
10
|
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
|
11
|
Evangelista-Falcón W, Denhez C, Baena-Moncada A, Ponce-Vargas M. Revisiting the Sweet Taste Receptor T1R2-T1R3 through Molecular Dynamics Simulations Coupled with a Noncovalent Interactions Analysis. J Phys Chem B 2023; 127:1110-1119. [PMID: 36705604 DOI: 10.1021/acs.jpcb.2c07180] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
It is nowadays widely accepted that sweet taste perception is elicited by the activation of the heterodimeric complex T1R2-T1R3, customarily known as sweet taste receptor (STR). However, the interplay between STR and sweeteners has not yet been fully clarified. Here through a methodology coupling molecular dynamics and the independent gradient model (igm) approach we determine the main interacting signatures of the closed (active) conformation of the T1R2 Venus flytrap domain (VFD) toward aspartame. The igm methodology provides a rapid and reliable quantification of noncovalent interactions through a score (Δginter) based on the attenuation of the electronic density gradient when two molecular fragments approach each other. Herein, this approach is coupled to a 100 ns molecular dynamics simulation (MD-igm) to explore the ligand-cavity contacts on a per-residue basis as well as a series of key inter-residue interactions that stabilize the closed form of VFD. We also apply an atomic decomposition scheme of noncovalent interactions to quantify the contribution of the ligand segments to the noncovalent interplay. Finally, a series of structural modification on aspartame are conducted in order to obtain guidelines for the rational design of novel sweeteners. Given that innovative methodologies to reliably quantify the extent of ligand-protein coupling are strongly demanded, this approach combining a noncovalent analysis and MD simulations represents a valuable contribution, that can be easily applied to other relevant biomolecular systems.
Collapse
Affiliation(s)
- Wilfredo Evangelista-Falcón
- Laboratory of Biomolecules, Faculty of Health Sciences, Universidad Peruana de Ciencias Aplicadas, Lima15023, Perú
| | - Clément Denhez
- Institut de Chimie Moléculaire de Reims UMR CNRS 7312, Université de Reims Champagne-Ardenne, Moulin de la Housse 51687, ReimsCedex 02 BP39, France
| | - Angélica Baena-Moncada
- Laboratorio de Investigación de Electroquímica Aplicada, Facultad de Ciencias de la Universidad Nacional de Ingeniería, Av. Túpac Amaru 210, Rímac, Lima31-139, Perú
| | - Miguel Ponce-Vargas
- Institut de Chimie Moléculaire de Reims UMR CNRS 7312, Université de Reims Champagne-Ardenne, Moulin de la Housse 51687, ReimsCedex 02 BP39, France
| |
Collapse
|
12
|
Gomari MM, Tarighi P, Choupani E, Abkhiz S, Mohamadzadeh M, Rostami N, Sadroddiny E, Baammi S, Uversky VN, Dokholyan NV. Structural evolution of Delta lineage of SARS-CoV-2. Int J Biol Macromol 2023; 226:1116-1140. [PMID: 36435470 PMCID: PMC9683856 DOI: 10.1016/j.ijbiomac.2022.11.227] [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: 08/17/2022] [Revised: 11/19/2022] [Accepted: 11/21/2022] [Indexed: 11/27/2022]
Abstract
One of the main obstacles in prevention and treatment of COVID-19 is the rapid evolution of the SARS-CoV-2 Spike protein. Given that Spike is the main target of common treatments of COVID-19, mutations occurring at this virulent factor can affect the effectiveness of treatments. The B.1.617.2 lineage of SARS-CoV-2, being characterized by many Spike mutations inside and outside of its receptor-binding domain (RBD), shows high infectivity and relative resistance to existing cures. Here, utilizing a wide range of computational biology approaches, such as immunoinformatics, molecular dynamics (MD), analysis of intrinsically disordered regions (IDRs), protein-protein interaction analyses, residue scanning, and free energy calculations, we examine the structural and biological attributes of the B.1.617.2 Spike protein. Furthermore, the antibody design protocol of Rosetta was implemented for evaluation the stability and affinity improvement of the Bamlanivimab (LY-CoV55) antibody, which is not capable of interactions with the B.1.617.2 Spike. We observed that the detected mutations in the Spike of the B1.617.2 variant of concern can cause extensive structural changes compatible with the described variation in immunogenicity, secondary and tertiary structure, oligomerization potency, Furin cleavability, and drug targetability. Compared to the Spike of Wuhan lineage, the B.1.617.2 Spike is more stable and binds to the Angiotensin-converting enzyme 2 (ACE2) with higher affinity.
Collapse
Affiliation(s)
- Mohammad Mahmoudi Gomari
- Student Research Committee, Iran University of Medical Sciences, Tehran 1449614535, Iran,Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Parastoo Tarighi
- Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Edris Choupani
- Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Shadi Abkhiz
- Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Masoud Mohamadzadeh
- Department of Chemistry, Faculty of Sciences, University of Hormozgan, Bandar Abbas 7916193145, Iran
| | - Neda Rostami
- Department of Chemical Engineering, Faculty of Engineering, Arak University, Arak 3848177584, Iran
| | - Esmaeil Sadroddiny
- Medical Biotechnology Department, School of Advanced Technologies in Medicine, Tehran University of Medical Sciences, Tehran 1417613151, Iran
| | - Soukayna Baammi
- African Genome Centre (AGC), Mohammed VI Polytechnic University, Benguerir 43150, Morocco
| | - Vladimir N. Uversky
- Department of Molecular Medicine and USF Health Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, Tampa, FL 33620, USA,Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia,Correspondence to: V.N. Uversky, Department of Molecular Medicine and USF Health Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, Tampa, FL 33620, USA
| | - Nikolay V. Dokholyan
- Department of Pharmacology, Department of Biochemistry & Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 16802, USA,Corresponding author
| |
Collapse
|
13
|
Ma J, Ma Y, Li Y, Sun Z, Sun X, Padmakumar V, Cheng Y, Zhu W. Characterization of feruloyl esterases from Pecoramyces sp. F1 and the synergistic effect in biomass degradation. World J Microbiol Biotechnol 2022; 39:17. [PMID: 36409385 DOI: 10.1007/s11274-022-03466-3] [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: 06/01/2022] [Accepted: 11/10/2022] [Indexed: 11/22/2022]
Abstract
Feruloyl esterase (FAE; EC 3.1.1.73) cleaves the ester bond between ferulic acid (FA) and sugar, to assist the release of FAs and degradation of plant cell walls. In this study, two FAEs (Fae13961 and Fae16537) from the anaerobic fungus Pecoramyces sp. F1 were heterologously expressed in Pichia pastoris (P. pastoris). Compared with Fae16537, Fae13961 had higher catalytic efficiency. The optimum temperature and pH of both the FAEs were 45 ℃ and 7.0, respectively. They showed good stability-Fae16537 retained up to 80% activity after incubation at 37 ℃ for 24 h. The FAEs activity was enhanced by Ca2+ and reduced by Zn2+, Mn2+, Fe2+ and Fe3+. Additionally, the effect of FAEs on the hydrolytic efficiency of xylanase and cellulase was also determined. The FAE Fae13961 had synergistic effect with xylanase and it promoted the degradation of xylan substrates by xylanase, but it did not affect the degradation of cellulose substrates by cellulase. When Fae13961 was added in a mixture of xylanase and cellulase to degrade complex agricultural biomass, it significantly enhanced the mixture's ability to disintegrate complex substrates. These FAEs could serve as superior auxiliary enzymes for other lignocellulosic enzymes in the process of degradation of agricultural residues for industrial applications.
Collapse
Affiliation(s)
- Jing Ma
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China
| | - Yuping Ma
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China
| | - Yuqi Li
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China
| | - Zhanying Sun
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China
| | - Xiaoni Sun
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China
| | | | - Yanfen Cheng
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China.
| | - Weiyun Zhu
- Laboratory of Gastrointestinal Microbiology, National Center for International Research on Animal Gut Nutrition, Nanjing Agricultural University, Nanjing, 210095, China
| |
Collapse
|
14
|
Yakiyama Y. Molecular-Shape-Organized Stimuli-Responsive Functional Crystalline Systems. J SYN ORG CHEM JPN 2022. [DOI: 10.5059/yukigoseikyokaishi.80.1036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Yumi Yakiyama
- Division of Applied Chemistry, Graduate School of Engineering, Osaka University
| |
Collapse
|
15
|
Ramasamy P, Vandermarliere E, Vranken WF, Martens L. Panoramic Perspective on Human Phosphosites. J Proteome Res 2022; 21:1894-1915. [PMID: 35793420 DOI: 10.1021/acs.jproteome.2c00164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein phosphorylation is the most common reversible post-translational modification of proteins and is key in the regulation of many cellular processes. Due to this importance, phosphorylation is extensively studied, resulting in the availability of a large amount of mass spectrometry-based phospho-proteomics data. Here, we leverage the information in these large-scale phospho-proteomics data sets, as contained in Scop3P, to analyze and characterize proteome-wide protein phosphorylation sites (P-sites). First, we set out to differentiate correctly observed P-sites from false-positive sites using five complementary site properties. We then describe the context of these P-sites in terms of the protein structure, solvent accessibility, structural transitions and disorder, and biophysical properties. We also investigate the relative prevalence of disease-linked mutations on and around P-sites. Moreover, we assess the structural dynamics of P-sites in their phosphorylated and unphosphorylated states. As a result, we show how large-scale reprocessing of available proteomics experiments can enable a more reliable view on proteome-wide P-sites. Furthermore, adding the structural context of proteins around P-sites helps uncover possible conformational switches upon phosphorylation. Moreover, by placing sites in different biophysical contexts, we show the differential preference in protein dynamics at phosphorylated sites when compared to the nonphosphorylated counterparts.
Collapse
Affiliation(s)
- Pathmanaban Ramasamy
- VIB-UGent Center for Medical Biotechnology, VIB, 9000 Ghent, Belgium.,Department of Biomolecular Medicine, Faculty of Health Sciences and Medicine, Ghent University, 9000 Ghent, Belgium.,Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, 1050 Brussels, Belgium.,Structural Biology Brussels, Vrije Universiteit Brussel, 1050 Brussels, Belgium.,Centre for Structural Biology, VIB, 1050 Brussels, Belgium
| | | | - Wim F Vranken
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, 1050 Brussels, Belgium.,Structural Biology Brussels, Vrije Universiteit Brussel, 1050 Brussels, Belgium.,Centre for Structural Biology, VIB, 1050 Brussels, Belgium
| | - Lennart Martens
- VIB-UGent Center for Medical Biotechnology, VIB, 9000 Ghent, Belgium.,Department of Biomolecular Medicine, Faculty of Health Sciences and Medicine, Ghent University, 9000 Ghent, Belgium
| |
Collapse
|
16
|
Padányi R, Farkas B, Tordai H, Kiss B, Grubmüller H, Soya N, Lukács GL, Kellermayer M, Hegedűs T. Nanomechanics combined with HDX reveals allosteric drug binding sites of CFTR NBD1. Comput Struct Biotechnol J 2022; 20:2587-2599. [PMID: 35685375 PMCID: PMC9160490 DOI: 10.1016/j.csbj.2022.05.036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Revised: 05/17/2022] [Accepted: 05/19/2022] [Indexed: 11/29/2022] Open
Abstract
Cystic fibrosis (CF) is a frequent genetic disease in Caucasians that is caused by the deletion of F508 (ΔF508) in the nucleotide binding domain 1 (NBD1) of the CF transmembrane conductance regulator (CFTR). The ΔF508 compromises the folding energetics of the NBD1, as well as the folding of three other CFTR domains. Combination of FDA approved corrector molecules can efficiently but incompletely rescue the ΔF508-CFTR folding and stability defect. Thus, new pharmacophores that would reinstate the wild-type-like conformational stability of the ΔF508-NBD1 would be highly beneficial. The most prominent molecule, 5-bromoindole-3-acetic acid (BIA) that can thermally stabilize the NBD1 has low potency and efficacy. To gain insights into the NBD1 (un)folding dynamics and BIA binding site localization, we combined molecular dynamics (MD) simulations, atomic force spectroscopy (AFM) and hydrogen-deuterium exchange (HDX) experiments. We found that the NBD1 α-subdomain with three adjacent strands from the β-subdomain plays an important role in early folding steps, when crucial non-native interactions are formed via residue F508. Our AFM and HDX experiments showed that BIA associates with this α-core region and increases the resistance of the ΔF508-NBD1 against mechanical unfolding, a phenomenon that could be exploited in future developments of folding correctors.
Collapse
Affiliation(s)
- Rita Padányi
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Bianka Farkas
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
- Faculty of Information Technology and Bionics, Pázmány Péter Catholic University, Budapest, Hungary
| | - Hedvig Tordai
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Bálint Kiss
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Helmut Grubmüller
- Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Naoto Soya
- Department of Physiology and Biochemistry, McGill University, Montréal, Quebec, Canada
| | - Gergely L. Lukács
- Department of Physiology and Biochemistry, McGill University, Montréal, Quebec, Canada
| | - Miklós Kellermayer
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Tamás Hegedűs
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
- ELKH-SE Molecular Biophysics Research Group, ELKH, Budapest, Hungary
- Corresponding author at: Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary.
| |
Collapse
|
17
|
Prabantu VM, Gadiyaram V, Vishveshwara S, Srinivasan N. Understanding structural variability in proteins using protein structural networks. Curr Res Struct Biol 2022; 4:134-145. [PMID: 35586857 PMCID: PMC9108755 DOI: 10.1016/j.crstbi.2022.04.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Revised: 03/01/2022] [Accepted: 04/09/2022] [Indexed: 11/13/2022] Open
Abstract
Proteins perform their function by accessing a suitable conformer from the ensemble of available conformations. The conformational diversity of a chosen protein structure can be obtained by experimental methods under different conditions. A key issue is the accurate comparison of different conformations. A gold standard used for such a comparison is the root mean square deviation (RMSD) between the two structures. While extensive refinements of RMSD evaluation at the backbone level are available, a comprehensive framework including the side chain interaction is not well understood. Here we employ protein structure network (PSN) formalism, with the non-covalent interactions of side chain, explicitly treated. The PSNs thus constructed are compared through graph spectral method, which provides a comparison at the local and at the global structural level. In this work, PSNs of multiple crystal conformers of single-chain, single-domain proteins, are subject to pair-wise analysis to examine the dissimilarity in their network topologies and in order to determine the conformational diversity of their native structures. This information is utilized to classify the structural domains of proteins into different categories. It is observed that proteins typically tend to retain structure and interactions at the backbone level. However, some of them also depict variability in either their overall structure or only in their inter-residue connectivity at the sidechain level, or both. Variability of sub-networks based on solvent accessibility and secondary structure is studied. The types of specific interactions are found to contribute differently to structure variability. An ensemble analysis by computing the mathematical variance of edge-weights across multiple conformers provided information on the contribution to overall variability from each edge of the PSN. Interactions that are highly variable are identified and their impact on structure variability has been discussed with the help of a case study. The classification based on the present side-chain network-based studies provides a framework to correlate the structure-function relationships in protein structures. Monomeric, single domain protein structures can exhibit non-rigid behaviour and be highly variable. The comparison of protein structural networks can better discriminate conformations with similar backbones. Specific interactions between solvent accessible and inaccessible residues are poorly preserved. Network edge-variation offers insights on which interacting residues are likely to influence their dynamics and function. These side-chain network-based studies provide a framework to correlate protein structure-function relationships.
Collapse
|
18
|
Cao M, Zhao P, Liu C, Xia J, Xu H. When Dynamic Diselenide Bonds meet Dynamic Imine Bonds in Polymeric Materials. Macromol Rapid Commun 2022; 43:e2200083. [PMID: 35257443 DOI: 10.1002/marc.202200083] [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: 01/30/2022] [Revised: 02/27/2022] [Indexed: 11/09/2022]
Abstract
In both natural and artificial functional systems, the cooperation between different dynamic interactions is of vital importance for realizing complicated functions. Dynamic covalent bonds are one kind of relatively stable dynamic interactions, and have shown synergistic effect in natural systems such as functional proteins. However, synergistic interactions between different dynamic covalent bonds in polymeric materials are still unclear. Herein, polymeric materials containing diselenide and imine bonds are prepared, and then the synergistic effect between the two dynamic covalent bonds is quantitatively evaluated in typical processes of dynamic materials. The results reveal that dynamic covalent bonds show weak synergistic effect in the degradation process, and have strong synergistic effect in stress relaxation process. Therefore, introducing multiple dynamic covalent bonds in polymeric materials could extensively enhance their dynamic properties. This article is protected by copyright. All rights reserved.
Collapse
Affiliation(s)
- Muqing Cao
- Key Lab of Organic Optoelectronics & Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, 100084, People's Republic of China
| | - Peng Zhao
- Key Lab of Organic Optoelectronics & Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, 100084, People's Republic of China
| | - Cheng Liu
- Key Lab of Organic Optoelectronics & Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, 100084, People's Republic of China
| | - Jiahao Xia
- Key Lab of Organic Optoelectronics & Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, 100084, People's Republic of China
| | - Huaping Xu
- Key Lab of Organic Optoelectronics & Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, 100084, People's Republic of China
| |
Collapse
|
19
|
Yadav S, Bharti S, Srivastava P, Mathur P. PepEngine: A Manually Curated Structural Database of Peptides Containing α, β- Dehydrophenylalanine (ΔPhe) and α-Amino Isobutyric Acid (Aib). Int J Pept Res Ther 2022. [DOI: 10.1007/s10989-022-10362-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
20
|
rsRNASP: A residue-separation-based statistical potential for RNA 3D structure evaluation. Biophys J 2022; 121:142-156. [PMID: 34798137 PMCID: PMC8758408 DOI: 10.1016/j.bpj.2021.11.016] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/23/2021] [Accepted: 11/10/2021] [Indexed: 01/07/2023] Open
Abstract
Knowledge-based statistical potentials have been shown to be rather effective in protein 3-dimensional (3D) structure evaluation and prediction. Recently, several statistical potentials have been developed for RNA 3D structure evaluation, while their performances are either still at a low level for the test datasets from structure prediction models or dependent on the "black-box" process through neural networks. In this work, we have developed an all-atom distance-dependent statistical potential based on residue separation for RNA 3D structure evaluation, namely rsRNASP, which is composed of short- and long-ranged potentials distinguished by residue separation. The extensive examinations against available RNA test datasets show that rsRNASP has apparently higher performance than the existing statistical potentials for the realistic test datasets with large RNAs from structure prediction models, including the newly released RNA-Puzzles dataset, and is comparable to the existing top statistical potentials for the test datasets with small RNAs or near-native decoys. In addition, rsRNASP is superior to RNA3DCNN, a recently developed scoring function through 3D convolutional neural networks. rsRNASP and the relevant databases are available to the public.
Collapse
|
21
|
Casier R, Duhamel J. Effects of Glycine on the Local Conformation and Internal Backbone Dynamics of Polypeptides. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c01479] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Remi Casier
- Institute for Polymer Research, Waterloo Institute for Nanotechnology, Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L3G1, Canada
| | - Jean Duhamel
- Institute for Polymer Research, Waterloo Institute for Nanotechnology, Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L3G1, Canada
| |
Collapse
|
22
|
The SARS-CoV-2 spike protein is vulnerable to moderate electric fields. Nat Commun 2021; 12:5407. [PMID: 34518528 PMCID: PMC8437970 DOI: 10.1038/s41467-021-25478-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 08/12/2021] [Indexed: 12/23/2022] Open
Abstract
Most of the ongoing projects aimed at the development of specific therapies and vaccines against COVID-19 use the SARS-CoV-2 spike (S) protein as the main target. The binding of the spike protein with the ACE2 receptor (ACE2) of the host cell constitutes the first and key step for virus entry. During this process, the receptor binding domain (RBD) of the S protein plays an essential role, since it contains the receptor binding motif (RBM), responsible for the docking to the receptor. So far, mostly biochemical methods are being tested in order to prevent binding of the virus to ACE2. Here we show, with the help of atomistic simulations, that external electric fields of easily achievable and moderate strengths can dramatically destabilise the S protein, inducing long-lasting structural damage. One striking field-induced conformational change occurs at the level of the recognition loop L3 of the RBD where two parallel beta sheets, believed to be responsible for a high affinity to ACE2, undergo a change into an unstructured coil, which exhibits almost no binding possibilities to the ACE2 receptor. We also show that these severe structural changes upon electric-field application also occur in the mutant RBDs corresponding to the variants of concern (VOC) B.1.1.7 (UK), B.1.351 (South Africa) and P.1 (Brazil). Remarkably, while the structural flexibility of S allows the virus to improve its probability of entering the cell, it is also the origin of the surprising vulnerability of S upon application of electric fields of strengths at least two orders of magnitude smaller than those required for damaging most proteins. Our findings suggest the existence of a clean physical method to weaken the SARS-CoV-2 virus without further biochemical processing. Moreover, the effect could be used for infection prevention purposes and also to develop technologies for in-vitro structural manipulation of S. Since the method is largely unspecific, it can be suitable for application to other mutations in S, to other proteins of SARS-CoV-2 and in general to membrane proteins of other virus types. The SARS-CoV-2 Spike protein is essential for viral infectivity and binds to the host receptor ACE2. Here, the authors present MD simulations of the Spike protein and its variants of concern and observe that the Spike protein is destabilised by moderate static electric fields, and undergoes field-induced conformational changes that hinder binding to ACE2.
Collapse
|
23
|
Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances. Life (Basel) 2021; 11:life11060570. [PMID: 34204410 PMCID: PMC8234428 DOI: 10.3390/life11060570] [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: 04/14/2021] [Revised: 05/30/2021] [Accepted: 06/09/2021] [Indexed: 11/17/2022] Open
Abstract
Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.
Collapse
|
24
|
Floch A, Pirenne F, Barrault A, Chami B, Toly-Ndour C, Tournamille C, de Brevern AG. Insights into anti-D formation in carriers of RhD variants through studies of 3D intraprotein interactions. Transfusion 2021; 61:1286-1301. [PMID: 33586199 DOI: 10.1111/trf.16301] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 11/04/2020] [Accepted: 01/13/2021] [Indexed: 01/11/2023]
Abstract
BACKGROUND Many RhD variants associated with anti-D formation (partial D) in carriers exposed to the conventional D antigen carry mutations affecting extracellular loop residues. Surprisingly, some carry mutations affecting transmembrane or intracellular domains, positions not thought likely to have a major impact on D epitopes. STUDY DESIGN AND METHODS A wild-type Rh trimer (RhD1 RhAG2 ) was modeled by comparative modeling with the human RhCG structure. Taking trimer conformation, residue accessibility, and position relative to the lipid bilayer into account, we redefine the domains of the RhD protein. We generated models for RhD variants carrying one or two amino acid substitutions associated with anti-D formation in published articles (25 variants) or abstracts (12 variants) and for RHD*weak D type 38. We determined the extracellular substitutions and compared the interactions of the variants with those of the standard RhD. RESULTS The findings of the three-dimensional (3D) analysis were correlated with anti-D formation for 76% of RhD variants: 15 substitutions associated with anti-D formation concerned extracellular residues, and structural differences in intraprotein interactions relative to standard RhD were observed in the others. We discuss the mechanisms by which D epitopes may be modified in variants in which the extracellular residues are identical to those of standard RhD and provide arguments for the benignity of p.T379M (RHD*DAU0) and p.G278D (RHD*weak D type 38) in transfusion medicine. CONCLUSION The study of RhD intraprotein interactions and the precise redefinition of residue accessibility provide insight into the mechanisms through which RhD point mutations may lead to anti-D formation in carriers.
Collapse
Affiliation(s)
- Aline Floch
- Univ Paris Est Creteil, INSERM U955, Transfusion et Maladies du Globule Rouge, IMRB, Creteil, France.,Etablissement francais du sang Ile-de-France, Creteil, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| | - France Pirenne
- Univ Paris Est Creteil, INSERM U955, Transfusion et Maladies du Globule Rouge, IMRB, Creteil, France.,Etablissement francais du sang Ile-de-France, Creteil, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| | - Aurélie Barrault
- Univ Paris Est Creteil, INSERM U955, Transfusion et Maladies du Globule Rouge, IMRB, Creteil, France.,Etablissement francais du sang Ile-de-France, Creteil, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| | - Btissam Chami
- Etablissement francais du sang Ile-de-France, Creteil, France
| | - Cécile Toly-Ndour
- Unité Fonctionnelle d'expertise en Immuno-Hémobiologie Périnatale, Centre National de Référence en Hémobiologie Périnatale (CNRHP), Service de Médecine Fœtale, Pôle Périnatalité, Hôpital Trousseau, GH HUEP, APHP, Paris, France
| | - Christophe Tournamille
- Univ Paris Est Creteil, INSERM U955, Transfusion et Maladies du Globule Rouge, IMRB, Creteil, France.,Etablissement francais du sang Ile-de-France, Creteil, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| | - Alexandre G de Brevern
- Laboratoire d'Excellence GR-Ex, Paris, France.,Université de Paris, Biologie Intégrée du Globule Rouge UMR_S1134, Inserm, Université de la Réunion, Université des Antilles, Paris, France.,Institut National de la Transfusion Sanguine (INTS), Paris, France
| |
Collapse
|
25
|
Chowdhury A, Chatterjee S, Pongen A, Sarania D, Tripathi NM, Bandyopadhyay A. nSite-Selective, Chemical Modification of Protein at Aromatic Side Chain and Their Emergent Applications. Protein Pept Lett 2021; 28:788-808. [PMID: 33511938 DOI: 10.2174/0929866528666210129152535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 12/11/2020] [Accepted: 12/20/2020] [Indexed: 11/22/2022]
Abstract
Site-selective chemical modification of protein side chain has probed enormous opportunities in the fundamental understanding of cellular biology and therapeutic applications. Primarily, in the field of biopharmaceutical where formulation of bioconjugates is found to be potential medicine than an individual constituent. In this regard, Lysine and Cysteine are the most widely used endogenous amino acid for these purposes. Recently, the aromatic side chain residues (Trp, Tyr, and His) that are low abundant in protein have gained more attention in therapeutic applications due to their advantages of chemical reactivity and specificity. This review discusses the site-selective bioconjugation methods for aromatic side chains (Trp, Tyr and His) and highlights the developed strategies in the last three years, along with their applications. Also, the review highlights the prevalent methods published earlier. We have examined that metal-catalyzed and photocatalytic reactions are gaining more attention for bioconjugation, though their practical operation is under development. The review has been summarized with the future perspective of protein and peptide conjugations contemplating therapeutic applications and challenges.
Collapse
Affiliation(s)
- Arnab Chowdhury
- Biomimetic Peptide Engineering Laboratory, Department of Chemistry, Indian Institute of Technology, Ropar, Birla Farms, Punjab-781039. India
| | - Saurav Chatterjee
- Biomimetic Peptide Engineering Laboratory, Department of Chemistry, Indian Institute of Technology, Ropar, Birla Farms, Punjab-781039. India
| | - Akumlong Pongen
- Biomimetic Peptide Engineering Laboratory, Department of Chemistry, Indian Institute of Technology, Ropar, Birla Farms, Punjab-781039. India
| | - Dhanjit Sarania
- Biomimetic Peptide Engineering Laboratory, Department of Chemistry, Indian Institute of Technology, Ropar, Birla Farms, Punjab-781039. India
| | - Nitesh Mani Tripathi
- Biomimetic Peptide Engineering Laboratory, Department of Chemistry, Indian Institute of Technology, Ropar, Birla Farms, Punjab-781039. India
| | - Anupam Bandyopadhyay
- Biomimetic Peptide Engineering Laboratory, Department of Chemistry, Indian Institute of Technology, Ropar, Birla Farms, Punjab-781039. India
| |
Collapse
|
26
|
Casier R, Duhamel J. Blob-Based Predictions of Protein Folding Times from the Amino Acid-Dependent Conformation of Polypeptides in Solution. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c02617] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Remi Casier
- Institute for Polymer Research, Waterloo Institute for Nanotechnology, Department of Chemistry, University of Waterloo, Waterloo, ON N2L3G1, Canada
| | - Jean Duhamel
- Institute for Polymer Research, Waterloo Institute for Nanotechnology, Department of Chemistry, University of Waterloo, Waterloo, ON N2L3G1, Canada
| |
Collapse
|
27
|
Casier R, Duhamel J. Blob-Based Approach to Estimate the Folding Time of Proteins Supported by Pyrene Excimer Fluorescence Experiments. Macromolecules 2020. [DOI: 10.1021/acs.macromol.0c02201] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Remi Casier
- Institute for Polymer Research, Waterloo Institute for Nanotechnology, Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Jean Duhamel
- Institute for Polymer Research, Waterloo Institute for Nanotechnology, Department of Chemistry, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
28
|
Kulandaisamy A, Zaucha J, Frishman D, Gromiha MM. MPTherm-pred: Analysis and Prediction of Thermal Stability Changes upon Mutations in Transmembrane Proteins. J Mol Biol 2020; 433:166646. [PMID: 32920050 DOI: 10.1016/j.jmb.2020.09.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2020] [Revised: 09/04/2020] [Accepted: 09/04/2020] [Indexed: 01/06/2023]
Abstract
The stability of membrane proteins differs from globular proteins due to the presence of nonpolar membrane-spanning regions. Using a dataset of 929 membrane protein mutations whose effects on thermal stability (ΔTm) were experimentally determined, we found that the average ΔTm due to 190 stabilizing and 232 destabilizing mutations occurring in membrane-spanning regions are 2.43(3.1) °C and -5.48(5.5) °C, respectively. The ΔTm values for mutations occurring in solvent-exposed regions are 2.56(2.82) and - 6.8(7.2) °C. We have systematically analyzed the factors influencing the stability of mutants and observed that changes in hydrophobicity, number of contacts between Cα atoms and frequency of aliphatic residues are important determinants of the stability change induced by mutations occurring in membrane-spanning regions. We have developed structure- and sequence-based machine learning predictors of ΔTm due to mutations specifically for membrane proteins. They showed a correlation and mean absolute error (MAE) of 0.72 and 2.85 °C, respectively, between experimental and predicted ΔTm for mutations in membrane-spanning regions on 10-fold group-wise cross-validation. The average correlation and MAE for mutations in aqueous regions are 0.73 and 3.7 °C, respectively. These MAE values are about 50% lower than standard deviations from the mean ΔTm values. The reliability of the method was affirmed on a test set of mutations occurring in evolutionary independent protein sequences. The developed MPTherm-pred server for predicting thermal stability changes upon mutations in membrane proteins is available at https://web.iitm.ac.in/bioinfo2/mpthermpred/. Our results provide insights into factors influencing the stability of membrane proteins and can aid in designing mutants that are more resistant to thermal stress.
Collapse
Affiliation(s)
- A Kulandaisamy
- Department of Biotechnology, Bhupat and Jyoti Mehta School of BioSciences, Indian Institute of Technology Madras, Chennai 600 036, Tamilnadu, India
| | - Jan Zaucha
- Department of Bioinformatics, Technische Universität München, Wissenschaftszentrum Weihenstephan, Freising, Germany
| | - Dmitrij Frishman
- Department of Bioinformatics, Technische Universität München, Wissenschaftszentrum Weihenstephan, Freising, Germany; Department of Bioinformatics, Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russian Federation
| | - M Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of BioSciences, Indian Institute of Technology Madras, Chennai 600 036, Tamilnadu, India.
| |
Collapse
|
29
|
Shao J, Liu B. ProtFold-DFG: protein fold recognition by combining Directed Fusion Graph and PageRank algorithm. Brief Bioinform 2020; 22:5901980. [PMID: 32892224 DOI: 10.1093/bib/bbaa192] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Revised: 07/16/2020] [Accepted: 07/28/2020] [Indexed: 12/27/2022] Open
Abstract
As one of the most important tasks in protein structure prediction, protein fold recognition has attracted more and more attention. In this regard, some computational predictors have been proposed with the development of machine learning and artificial intelligence techniques. However, these existing computational methods are still suffering from some disadvantages. In this regard, we propose a new network-based predictor called ProtFold-DFG for protein fold recognition. We propose the Directed Fusion Graph (DFG) to fuse the ranking lists generated by different methods, which employs the transitive closure to incorporate more relationships among proteins and uses the KL divergence to calculate the relationship between two proteins so as to improve its generalization ability. Finally, the PageRank algorithm is performed on the DFG to accurately recognize the protein folds by considering the global interactions among proteins in the DFG. Tested on a widely used and rigorous benchmark data set, LINDAHL dataset, experimental results show that the ProtFold-DFG outperforms the other 35 competing methods, indicating that ProtFold-DFG will be a useful method for protein fold recognition. The source code and data of ProtFold-DFG can be downloaded from http://bliulab.net/ProtFold-DFG/download.
Collapse
Affiliation(s)
- Jiangyi Shao
- School of Computer Science and Technology, Beijing Institute of Technology, China
| | - Bin Liu
- School of Computer Science and Technology, Beijing Institute of Technology, Beijing, China
| |
Collapse
|
30
|
Darmawan KK, Karagiannis TC, Hughes JG, Small DM, Hung A. High temperature induced structural changes of apo-lactoferrin and interactions with β-lactoglobulin and α-lactalbumin for potential encapsulation strategies. Food Hydrocoll 2020. [DOI: 10.1016/j.foodhyd.2020.105817] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
31
|
Liu XR, Zhang MM, Gross ML. Mass Spectrometry-Based Protein Footprinting for Higher-Order Structure Analysis: Fundamentals and Applications. Chem Rev 2020; 120:4355-4454. [PMID: 32319757 PMCID: PMC7531764 DOI: 10.1021/acs.chemrev.9b00815] [Citation(s) in RCA: 130] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Proteins adopt different higher-order structures (HOS) to enable their unique biological functions. Understanding the complexities of protein higher-order structures and dynamics requires integrated approaches, where mass spectrometry (MS) is now positioned to play a key role. One of those approaches is protein footprinting. Although the initial demonstration of footprinting was for the HOS determination of protein/nucleic acid binding, the concept was later adapted to MS-based protein HOS analysis, through which different covalent labeling approaches "mark" the solvent accessible surface area (SASA) of proteins to reflect protein HOS. Hydrogen-deuterium exchange (HDX), where deuterium in D2O replaces hydrogen of the backbone amides, is the most common example of footprinting. Its advantage is that the footprint reflects SASA and hydrogen bonding, whereas one drawback is the labeling is reversible. Another example of footprinting is slow irreversible labeling of functional groups on amino acid side chains by targeted reagents with high specificity, probing structural changes at selected sites. A third footprinting approach is by reactions with fast, irreversible labeling species that are highly reactive and footprint broadly several amino acid residue side chains on the time scale of submilliseconds. All of these covalent labeling approaches combine to constitute a problem-solving toolbox that enables mass spectrometry as a valuable tool for HOS elucidation. As there has been a growing need for MS-based protein footprinting in both academia and industry owing to its high throughput capability, prompt availability, and high spatial resolution, we present a summary of the history, descriptions, principles, mechanisms, and applications of these covalent labeling approaches. Moreover, their applications are highlighted according to the biological questions they can answer. This review is intended as a tutorial for MS-based protein HOS elucidation and as a reference for investigators seeking a MS-based tool to address structural questions in protein science.
Collapse
Affiliation(s)
| | | | - Michael L. Gross
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO, USA, 63130
| |
Collapse
|
32
|
Sheik Amamuddy O, Musyoka TM, Boateng RA, Zabo S, Tastan Bishop Ö. Determining the unbinding events and conserved motions associated with the pyrazinamide release due to resistance mutations of Mycobacterium tuberculosis pyrazinamidase. Comput Struct Biotechnol J 2020; 18:1103-1120. [PMID: 32489525 PMCID: PMC7251373 DOI: 10.1016/j.csbj.2020.05.009] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Revised: 04/23/2020] [Accepted: 05/06/2020] [Indexed: 01/04/2023] Open
Abstract
Pyrazinamide (PZA) is the only first-line antitubercular drug active against latent Mycobacterium tuberculosis (Mtb). It is activated to pyrazinoic acid by the pncA-encoded pyrazinamidase enzyme (PZase). Despite the emergence of PZA drug resistance, the underlying mechanisms of resistance remain unclear. This study investigated part of these mechanisms by modelling a PZA-bound wild type and 82 mutant PZase structures before applying molecular dynamics (MD) with an accurate Fe2+ cofactor coordination geometry. After observing nanosecond-scale PZA unbinding from several PZase mutants, an algorithm was developed to systematically detect ligand release via centre of mass distances (COM) and ligand average speed calculations, before applying the statistically guided network analysis (SGNA) method to investigate conserved protein motions associated with ligand unbinding. Ligand and cofactor perspectives were also investigated. A conserved pair of lid-destabilising motions was found. These consisted of (1) antiparallel lid and side flap motions; (2) the contractions of a flanking region within the same flap and residue 74 towards the core. Mutations affecting the hinge residues (H51 and H71), nearby residues or L19 were found to destabilise the lid. Additionally, other metal binding site (MBS) mutations delocalised the Fe2+ cofactor, also facilitating lid opening. In the early stages of unbinding, a wider variety of PZA poses were observed, suggesting multiple exit pathways. These findings provide insights into the late events preceding PZA unbinding, which we found to occur in some resistant PZase mutants. Further, the algorithm developed here to identify unbinding events coupled with SGNA can be applicable to other similar problems.
Collapse
Key Words
- 3D, Three-dimensional
- ACPYPE, AnteChamber Python Parser interface
- Amber force field parameters
- CHPC, Center for High Performance Computing
- COM, Center of mass
- Drug resistance
- Drug unbinding
- FDA, Food and Drug Administration
- HTMD, High throughput molecular dynamics
- INH, Isoniazid
- MBS, Metal binding site
- MCBP, Metal Center Parameter Builder
- MD, Molecular dynamics
- MDR-TB, Multidrug-resistant tuberculosis
- Missense mutations
- Molecular dynamics simulations
- PBC, Periodic boundary conditions
- PDB, Protein Data bank
- POA, Pyrazinoic acid
- PZA, Pyrazinamide
- PZase, Pyrazinamidase
- QM, Quantum Mechanics
- RIF, Rifampicin
- SGNA, Statistically guided network analysis
- Statistically guided network analysis
- TB, Tuberculosis
- VAPOR, Variant Analysis Portal
- WHO, World Health Organization
- WT, Wild type
Collapse
Affiliation(s)
| | | | - Rita Afriyie Boateng
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| | - Sophakama Zabo
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| |
Collapse
|
33
|
Jing X, Zeng H, Wang S, Xu J. A Web-Based Protocol for Interprotein Contact Prediction by Deep Learning. Methods Mol Biol 2020; 2074:67-80. [PMID: 31583631 DOI: 10.1007/978-1-4939-9873-9_6] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Identifying residue-residue contacts in protein-protein interactions or complex is crucial for understanding protein and cell functions. DCA (direct-coupling analysis) methods shed some light on this, but they need many sequence homologs to yield accurate prediction. Inspired by the success of our deep-learning method for intraprotein contact prediction, we have developed RaptorX-ComplexContact, a web server for interprotein residue-residue contact prediction. Given a pair of interacting protein sequences, RaptorX-ComplexContact first searches for their sequence homologs and builds two paired multiple sequence alignments (MSA) based on genomic distance and phylogeny information, respectively. Then, RaptorX-ComplexContact uses two deep convolutional residual neural networks (ResNet) to predict interprotein contacts from sequential features and coevolution information of paired MSAs. RaptorX-ComplexContact shall be useful for protein docking, protein-protein interaction prediction, and protein interaction network construction.
Collapse
Affiliation(s)
- Xiaoyang Jing
- Toyota Technological Institute at Chicago, Chicago, IL, USA
- School of Computer Science, Fudan University, Shanghai, China
| | - Hong Zeng
- School of Computer Science and Technology, Hangzhou Dianzi University, Hangzhou, China
| | - Sheng Wang
- Toyota Technological Institute at Chicago, Chicago, IL, USA
| | - Jinbo Xu
- Toyota Technological Institute at Chicago, Chicago, IL, USA.
| |
Collapse
|
34
|
Summers TJ, Daniel BP, Cheng Q, DeYonker NJ. Quantifying Inter-Residue Contacts through Interaction Energies. J Chem Inf Model 2019; 59:5034-5044. [DOI: 10.1021/acs.jcim.9b00804] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Thomas J. Summers
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Baty P. Daniel
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Qianyi Cheng
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Nathan J. DeYonker
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| |
Collapse
|
35
|
Zhang H, Zhang Q, Ju F, Zhu J, Gao Y, Xie Z, Deng M, Sun S, Zheng WM, Bu D. Predicting protein inter-residue contacts using composite likelihood maximization and deep learning. BMC Bioinformatics 2019; 20:537. [PMID: 31664895 PMCID: PMC6821021 DOI: 10.1186/s12859-019-3051-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Accepted: 08/22/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Accurate prediction of inter-residue contacts of a protein is important to calculating its tertiary structure. Analysis of co-evolutionary events among residues has been proved effective in inferring inter-residue contacts. The Markov random field (MRF) technique, although being widely used for contact prediction, suffers from the following dilemma: the actual likelihood function of MRF is accurate but time-consuming to calculate; in contrast, approximations to the actual likelihood, say pseudo-likelihood, are efficient to calculate but inaccurate. Thus, how to achieve both accuracy and efficiency simultaneously remains a challenge. RESULTS In this study, we present such an approach (called clmDCA) for contact prediction. Unlike plmDCA using pseudo-likelihood, i.e., the product of conditional probability of individual residues, our approach uses composite-likelihood, i.e., the product of conditional probability of all residue pairs. Composite likelihood has been theoretically proved as a better approximation to the actual likelihood function than pseudo-likelihood. Meanwhile, composite likelihood is still efficient to maximize, thus ensuring the efficiency of clmDCA. We present comprehensive experiments on popular benchmark datasets, including PSICOV dataset and CASP-11 dataset, to show that: i) clmDCA alone outperforms the existing MRF-based approaches in prediction accuracy. ii) When equipped with deep learning technique for refinement, the prediction accuracy of clmDCA was further significantly improved, suggesting the suitability of clmDCA for subsequent refinement procedure. We further present a successful application of the predicted contacts to accurately build tertiary structures for proteins in the PSICOV dataset. CONCLUSIONS Composite likelihood maximization algorithm can efficiently estimate the parameters of Markov Random Fields and can improve the prediction accuracy of protein inter-residue contacts.
Collapse
Affiliation(s)
- Haicang Zhang
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Qi Zhang
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Fusong Ju
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Jianwei Zhu
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Yujuan Gao
- Center for Quantitative Biology, School of Mathematical Sciences, Center for Statistical Sciences, Peking University, Beijing, China
| | - Ziwei Xie
- College of Life Science and Technology, Huazhong University of Science and Technology, Wuhan, China
| | - Minghua Deng
- Center for Quantitative Biology, School of Mathematical Sciences, Center for Statistical Sciences, Peking University, Beijing, China
| | - Shiwei Sun
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China.
| | - Wei-Mou Zheng
- Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China.
| | - Dongbo Bu
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China. .,University of Chinese Academy of Sciences, Beijing, China.
| |
Collapse
|
36
|
Chandrasekaran P, Santosh Kumar C, Rangachari K, Sekar K. Disassociation of β1-α1-β2 from the α2-α3 domain of prion protein (PrP) is a prerequisite for the conformational conversion of PrPC into PrPSc: Driven by the free energy landscape. Int J Biol Macromol 2019; 136:368-376. [DOI: 10.1016/j.ijbiomac.2019.06.099] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2019] [Revised: 05/29/2019] [Accepted: 06/13/2019] [Indexed: 12/16/2022]
|
37
|
Mayol E, Campillo M, Cordomí A, Olivella M. Inter-residue interactions in alpha-helical transmembrane proteins. Bioinformatics 2019; 35:2578-2584. [PMID: 30566615 DOI: 10.1093/bioinformatics/bty978] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2018] [Revised: 10/19/2018] [Accepted: 12/17/2018] [Indexed: 01/23/2023] Open
Abstract
MOTIVATION The number of available membrane protein structures has markedly increased in the last years and, in parallel, the reliability of the methods to detect transmembrane (TM) segments. In the present report, we characterized inter-residue interactions in α-helical membrane proteins using a dataset of 3462 TM helices from 430 proteins. This is by far the largest analysis published to date. RESULTS Our analysis of residue-residue interactions in TM segments of membrane proteins shows that almost all interactions involve aliphatic residues and Phe. There is lack of polar-polar, polar-charged and charged-charged interactions except for those between Thr or Ser sidechains and the backbone carbonyl of aliphatic and Phe residues. The results are discussed in the context of the preferences of amino acids to be in the protein core or exposed to the lipid bilayer and to occupy specific positions along the TM segment. Comparison to datasets of β-barrel membrane proteins and of α-helical globular proteins unveils the specific patterns of interactions and residue composition characteristic of α-helical membrane proteins that are the clue to understanding their structure. AVAILABILITY AND IMPLEMENTATION Results data and datasets used are available at http://lmc.uab.cat/TMalphaDB/interactions.php. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Eduardo Mayol
- Laboratori de Medicina Computacional, Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Mercedes Campillo
- Laboratori de Medicina Computacional, Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Arnau Cordomí
- Laboratori de Medicina Computacional, Unitat de Bioestadística, Facultat de Medicina, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Mireia Olivella
- Bioinformatics Area, School of International Studies, ESCI-UPF, Barcelona, Spain.,Bioinformatics and Medical Statistics Group, U Science Tech, Central University of Catalonia, Vic, Barcelona, Spain
| |
Collapse
|
38
|
The Hsp70 Chaperone System Stabilizes a Thermo-sensitive Subproteome in E. coli. Cell Rep 2019; 28:1335-1345.e6. [DOI: 10.1016/j.celrep.2019.06.081] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2019] [Revised: 06/01/2019] [Accepted: 06/21/2019] [Indexed: 01/05/2023] Open
|
39
|
Luttrell J, Liu T, Zhang C, Wang Z. Predicting protein residue-residue contacts using random forests and deep networks. BMC Bioinformatics 2019; 20:100. [PMID: 30871477 PMCID: PMC6419322 DOI: 10.1186/s12859-019-2627-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The ability to predict which pairs of amino acid residues in a protein are in contact with each other offers many advantages for various areas of research that focus on proteins. For example, contact prediction can be used to reduce the computational complexity of predicting the structure of proteins and even to help identify functionally important regions of proteins. These predictions are becoming especially important given the relatively low number of experimentally determined protein structures compared to the amount of available protein sequence data. RESULTS Here we have developed and benchmarked a set of machine learning methods for performing residue-residue contact prediction, including random forests, direct-coupling analysis, support vector machines, and deep networks (stacked denoising autoencoders). These methods are able to predict contacting residue pairs given only the amino acid sequence of a protein. According to our own evaluations performed at a resolution of +/- two residues, the predictors we trained with the random forest algorithm were our top performing methods with average top 10 prediction accuracy scores of 85.13% (short range), 74.49% (medium range), and 54.49% (long range). Our ensemble models (stacked denoising autoencoders combined with support vector machines) were our best performing deep network predictors and achieved top 10 prediction accuracy scores of 75.51% (short range), 60.26% (medium range), and 43.85% (long range) using the same evaluation. These tests were blindly performed on targets from the CASP11 dataset; and the results suggested that our models achieved comparable performance to contact predictors developed by groups that participated in CASP11. CONCLUSIONS Due to the challenging nature of contact prediction, it is beneficial to develop and benchmark a variety of different prediction methods. Our work has produced useful tools with a simple interface that can provide contact predictions to users without requiring a lengthy installation process. In addition to this, we have released our C++ implementation of the direct-coupling analysis method as a standalone software package. Both this tool and our RFcon web server are freely available to the public at http://dna.cs.miami.edu/RFcon /.
Collapse
Affiliation(s)
- Joseph Luttrell
- School of Computing Sciences and Computer Engineering, University of Southern Mississippi, 118 College Drive, Hattiesburg, MS, 39406, USA
| | - Tong Liu
- Department of Computer Science, University of Miami, 1365 Memorial Drive, Coral Gables, FL, 33124, USA
| | - Chaoyang Zhang
- School of Computing Sciences and Computer Engineering, University of Southern Mississippi, 118 College Drive, Hattiesburg, MS, 39406, USA
| | - Zheng Wang
- Department of Computer Science, University of Miami, 1365 Memorial Drive, Coral Gables, FL, 33124, USA.
| |
Collapse
|
40
|
Wuyun Q, Zheng W, Peng Z, Yang J. A large-scale comparative assessment of methods for residue-residue contact prediction. Brief Bioinform 2019; 19:219-230. [PMID: 27802931 DOI: 10.1093/bib/bbw106] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2016] [Indexed: 11/14/2022] Open
Abstract
Sequence-based prediction of residue-residue contact in proteins becomes increasingly more important for improving protein structure prediction in the big data era. In this study, we performed a large-scale comparative assessment of 15 locally installed contact predictors. To assess these methods, we collected a big data set consisting of 680 nonredundant proteins covering different structural classes and target difficulties. We investigated a wide range of factors that may influence the precision of contact prediction, including target difficulty, structural class, the alignment depth and distribution of contact pairs in a protein structure. We found that: (1) the machine learning-based methods outperform the direct-coupling-based methods for short-range contact prediction, while the latter are significantly better for long-range contact prediction. The consensus-based methods, which combine machine learning and direct-coupling methods, perform the best. (2) The target difficulty does not have clear influence on the machine learning-based methods, while it does affect the direct-coupling and consensus-based methods significantly. (3) The alignment depth has relatively weak effect on the machine learning-based methods. However, for the direct-coupling-based methods and consensus-based methods, the predicted contacts for targets with deeper alignment tend to be more accurate. (4) All methods perform relatively better on β and α + β proteins than on α proteins. (5) Residues buried in the core of protein structure are more prone to be in contact than residues on the surface (22 versus 6%). We believe these are useful results for guiding future development of new approach to contact prediction.
Collapse
Affiliation(s)
- Qiqige Wuyun
- School of Mathematical Sciences, Nankai University, Tianjin, China
| | - Wei Zheng
- School of Mathematical Sciences, Nankai University, Tianjin, China
| | - Zhenling Peng
- Center for Applied Mathematics, Tianjin University, Tianjin, China
| | - Jianyi Yang
- School of Mathematical Sciences, Nankai University, Tianjin, China
| |
Collapse
|
41
|
Joshi PN, Rai V. Single-site labeling of histidine in proteins, on-demand reversibility, and traceless metal-free protein purification. Chem Commun (Camb) 2019; 55:1100-1103. [PMID: 30620346 DOI: 10.1039/c8cc08733d] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
A precision methodology distinguishes one His from all the nucleophilic residues and its multiple copies. An easy-to-operate C-N bond formation labels diverse proteins without adversely affecting their structure and function. The late-stage transformation allows installation of distinct probes. The chemically triggered reversibility enables traceless metal-free purification of proteins with a His-tag.
Collapse
Affiliation(s)
- Pralhad Namdev Joshi
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhauri, Bhopal, MP 462 066, India.
| | - Vishal Rai
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhauri, Bhopal, MP 462 066, India.
| |
Collapse
|
42
|
Amala A, Emerson IA. Understanding contact patterns of protein structures from protein contact map and investigation of unique patterns in the globin-like folded domains. J Cell Biochem 2018; 120:9877-9886. [PMID: 30525229 DOI: 10.1002/jcb.28270] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2018] [Accepted: 10/24/2018] [Indexed: 11/06/2022]
Abstract
Proteins are biochemical compounds made up of one or more polypeptides in a specific order, typically folded into a functionally active form. Proteins are categorized into four different structural classes according to the topology of α-helices and β-strands. In this study, we modeled these four structural classes as an undirected network depicting amino acids as nodes and interaction between them as edges. Results infer that basic protein classes can be easily recognized as well as distinguished by utilizing protein contact maps (PCM). Toward studying the globin-like fold, the helix-loop-helix region contacts were seen to be of a unique pattern, and these remained in all the folds. Further, the averaged diagonal contacts were analyzed and identified those contacts in α/β proteins were higher in comparison with the other class. Interesting, we noticed that anti-parallel beta sheets were dominant in all-β and α + β classes that lead to similar diagonal patterns. Network properties of all four basic classes were analyzed and found to possess small-world property. Findings infer that PCM may assist classify protein structure classes and it also helps in evaluating the predicted protein structures.
Collapse
Affiliation(s)
- Arumugam Amala
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Tamil Nadu, India
| | - Isaac Arnold Emerson
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Tamil Nadu, India
| |
Collapse
|
43
|
Kalhor H, Poorebrahim M, Rahimi H, Shabani AA, Karimipoor M, Akbari Eidgahi MR, Teimoori-Toolabi L. Structural and dynamic characterization of human Wnt2-Fzd7 complex using computational approaches. J Mol Model 2018; 24:274. [PMID: 30191337 DOI: 10.1007/s00894-018-3788-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Accepted: 08/09/2018] [Indexed: 12/20/2022]
Abstract
Wnt and Frizzled (Fzd) family members play crucial roles in the self-renewal of tumor-initiating cells. Until now, only a few studies have addressed the distinct mechanism of Wnt-Fzd interactions. In this study, we suggest a possible interaction mode of Wnt2 with the Fzd7 cysteine-rich domain (CRD)-both of which are up-regulated in some types of cancer. A combination of homology modeling, molecular docking and molecular dynamics (MD) simulations was carried out to study this ligand-receptor complex in great detail. The results demonstrated the unique dynamic behavior of Wnt2 upon binding to Fzd7. Interestingly, the β-strand content of the C-terminal binding site of Wnt2 was obviously reduced when bound to Fzd7 CRD. Moreover, the N-terminal and C-terminal binding sites of Wnt2 appeared to interact with the C-terminal and N-terminal binding sites of Fzd7, respectively. Calculation of the binding energies uncovered the pivotal role of electrostatic and hydrophobic interactions in the binding of Wnt2 to Fzd7 CRD. In conclusion, this study provides valuable insights into the mechanism of the Wnt2-Fzd7 CRD interaction for application in colorectal cancer prevention programs. Graphical abstract Flowchart representation of different steps used in this study.
Collapse
Affiliation(s)
- Hourieh Kalhor
- Department and Biotechnology Research Center, Semnan University of Medical Sciences, Semnan, Iran.,Molecular Medicine Department, Pasteur Institute of Iran, Tehran, Iran
| | | | - Hamzeh Rahimi
- Molecular Medicine Department, Pasteur Institute of Iran, Tehran, Iran
| | - Ali Akbar Shabani
- Department and Biotechnology Research Center, Semnan University of Medical Sciences, Semnan, Iran
| | | | | | | |
Collapse
|
44
|
Li B, Fooksa M, Heinze S, Meiler J. Finding the needle in the haystack: towards solving the protein-folding problem computationally. Crit Rev Biochem Mol Biol 2018; 53:1-28. [PMID: 28976219 PMCID: PMC6790072 DOI: 10.1080/10409238.2017.1380596] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2017] [Revised: 08/22/2017] [Accepted: 09/13/2017] [Indexed: 12/22/2022]
Abstract
Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as "the protein folding problem," has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.
Collapse
Affiliation(s)
- Bian Li
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Michaela Fooksa
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
- Chemical and Physical Biology Graduate Program, Vanderbilt University, Nashville, TN, USA
| | - Sten Heinze
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
45
|
Bittner HJ, Guixà-González R, Hildebrand PW. Structural basis for the interaction of the beta-secretase with copper. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2018; 1860:1105-1113. [PMID: 29391167 DOI: 10.1016/j.bbamem.2018.01.019] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Revised: 01/21/2018] [Accepted: 01/23/2018] [Indexed: 12/18/2022]
Abstract
The β-secretase (BACE1) features a unique sulfur rich motif (M462xxxC466xxxM470xxxC474xxxC478) in its transmembrane helix (BACE1-TM) which is characteristic for proteins involved in copper ion storage and transport. While this motif has been shown to promote BACE1-TM trimerization and binding of copper ions in vitro, the structural basis for the interaction of copper ions with the BACE1-TM is still not well understood. Using molecular dynamics (MD) simulations, we show that membrane embedded BACE1-TMs adopt a flexible trimeric structure that binds and conducts copper ions through variable coordination. In coarse-grained (CG) MD simulations, the spontaneous assembly of BACE1-TMs trimers results in a right-handed helix packing arrangement. In subsequent atomistic MD simulations the sulfur rich motif defines characteristic copper ion coordination sites along a constricted partially solvated axial pore. Sliding and tilting of BACE1-TMs along smooth A459xxxA463/464xxA467 surfaces, facilitated by a central P472 induced kink, enables copper ions to alternate between different coordination sites, including the prominent C466 and M470. We shed light into the structural arrangement of BACE1-TM trimers and propose a mechanism for copper ion conduction that might also apply to other proteins involved in metal ion transport.
Collapse
Affiliation(s)
- Heiko J Bittner
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institute of Medical Physics and Biophysics {Institut für Medizinische Physik und Biophysik}, AG ProteInFormatics, Charitéplatz 1, 10117 Berlin, Germany; Faculty of Medicine, Leipzig University, Institute for Medical Physics and Biophysics {Institut für Medizinische Physik und Biophysik}, Härtelstr. 16-18, 04107 Leipzig, Germany.
| | - Ramon Guixà-González
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institute of Medical Physics and Biophysics {Institut für Medizinische Physik und Biophysik}, AG ProteInFormatics, Charitéplatz 1, 10117 Berlin, Germany.
| | - Peter W Hildebrand
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institute of Medical Physics and Biophysics {Institut für Medizinische Physik und Biophysik}, AG ProteInFormatics, Charitéplatz 1, 10117 Berlin, Germany; Faculty of Medicine, Leipzig University, Institute for Medical Physics and Biophysics {Institut für Medizinische Physik und Biophysik}, Härtelstr. 16-18, 04107 Leipzig, Germany.
| |
Collapse
|
46
|
Molecular Chaperones: Structure-Function Relationship and their Role in Protein Folding. REGULATION OF HEAT SHOCK PROTEIN RESPONSES 2018. [DOI: 10.1007/978-3-319-74715-6_8] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
47
|
GRAFENE: Graphlet-based alignment-free network approach integrates 3D structural and sequence (residue order) data to improve protein structural comparison. Sci Rep 2017; 7:14890. [PMID: 29097661 PMCID: PMC5668259 DOI: 10.1038/s41598-017-14411-y] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2016] [Accepted: 10/11/2017] [Indexed: 12/26/2022] Open
Abstract
Initial protein structural comparisons were sequence-based. Since amino acids that are distant in the sequence can be close in the 3-dimensional (3D) structure, 3D contact approaches can complement sequence approaches. Traditional 3D contact approaches study 3D structures directly and are alignment-based. Instead, 3D structures can be modeled as protein structure networks (PSNs). Then, network approaches can compare proteins by comparing their PSNs. These can be alignment-based or alignment-free. We focus on the latter. Existing network alignment-free approaches have drawbacks: 1) They rely on naive measures of network topology. 2) They are not robust to PSN size. They cannot integrate 3) multiple PSN measures or 4) PSN data with sequence data, although this could improve comparison because the different data types capture complementary aspects of the protein structure. We address this by: 1) exploiting well-established graphlet measures via a new network alignment-free approach, 2) introducing normalized graphlet measures to remove the bias of PSN size, 3) allowing for integrating multiple PSN measures, and 4) using ordered graphlets to combine the complementary PSN data and sequence (specifically, residue order) data. We compare synthetic networks and real-world PSNs more accurately and faster than existing network (alignment-free and alignment-based), 3D contact, or sequence approaches.
Collapse
|
48
|
Borrel A, Camproux AC, Xhaard H. Characterization of Ionizable Groups' Environments in Proteins and Protein-Ligand Complexes through a Statistical Analysis of the Protein Data Bank. ACS OMEGA 2017; 2:7359-7374. [PMID: 31457307 PMCID: PMC6645025 DOI: 10.1021/acsomega.7b00739] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Accepted: 10/11/2017] [Indexed: 06/10/2023]
Abstract
We conduct a statistical analysis of the molecular environment of common ionizable functional groups in both protein-ligand complexes and inside proteins from the Protein Data Bank (PDB). In particular, we characterize the frequency, type, and density of the interacting atoms as well as the presence of a potential counterion. We found that for ligands, most guanidinium groups, half of primary and secondary amines, and one-fourth of imidazole neighbor a carboxylate group. Tertiary amines bind more rarely near carboxylate groups, which may be explained by a crowded neighborhood and hydrophobic character. In comparison to the environment seen by the ligands, inside proteins, an environment enriched in main-chain atoms is found, and the prevalence of direct charge neutralization by carboxylate groups is different. When the ionizable character of water molecules and phenolic or hydroxyl groups is accounted, considering a high-resolution dataset (less than 1.5 Å), charge neutralization could occur for well above 80% of the ligand functional groups considered, but for tertiary amines.
Collapse
Affiliation(s)
- Alexandre Borrel
- Molécules
Thérapeutiques in silico (MTi), INSERM UMRS-973, University Paris Diderot, Sorbonne Paris Cité, 75205 Paris Cedex 13, France
- Faculty
of Pharmacy, Division of Pharmaceutical Chemistry and Technology, University of Helsinki, Viikinkaari 5E, P.O. Box 56, FI-00014 Helsinki, Finland
| | - Anne-Claude Camproux
- Molécules
Thérapeutiques in silico (MTi), INSERM UMRS-973, University Paris Diderot, Sorbonne Paris Cité, 75205 Paris Cedex 13, France
| | - Henri Xhaard
- Faculty
of Pharmacy, Division of Pharmaceutical Chemistry and Technology, University of Helsinki, Viikinkaari 5E, P.O. Box 56, FI-00014 Helsinki, Finland
| |
Collapse
|
49
|
Khan I, Ansari IA, Singh P, Dass J FP. Prediction of functionally significant single nucleotide polymorphisms in PTEN tumor suppressor gene: An in silico approach. Biotechnol Appl Biochem 2017; 64:657-666. [PMID: 26800850 DOI: 10.1002/bab.1483] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2015] [Accepted: 01/16/2016] [Indexed: 11/06/2022]
Abstract
The phosphatase and tensin homolog (PTEN) gene plays a crucial role in signal transduction by negatively regulating the PI3K signaling pathway. It is the most frequent mutated gene in many human-related cancers. Considering its critical role, a functional analysis of missense mutations of PTEN gene was undertaken in this study. Thirty five nonsynonymous single nucleotide polymorphisms (nsSNPs) within the coding region of the PTEN gene were selected for our in silico investigation, and five nsSNPs (G129E, C124R, D252G, H61D, and R130G) were found to be deleterious based on combinatorial predictions of different computational tools. Moreover, molecular dynamics (MD) simulation was performed to investigate the conformational variation between native and all the five mutant PTEN proteins having predicted deleterious nsSNPs. The results of MD simulation of all mutant models illustrated variation in structural attributes such as root-mean-square deviation, root-mean-square fluctuation, radius of gyration, and total energy; which depicts the structural stability of PTEN protein. Furthermore, mutant PTEN protein structures also showed a significant variation in the solvent accessible surface area and hydrogen bond frequencies from the native PTEN structure. In conclusion, results of this study have established the deleterious effect of the all the five predicted nsSNPs on the PTEN protein structure. Thus, results of the current study can pave a new platform to sort out nsSNPs that can be undertaken for the confirmation of their phenotype and their correlation with diseased status in case of control studies.
Collapse
Affiliation(s)
- Imran Khan
- Department of Biosciences, Integral University, Lucknow, India
| | - Irfan A Ansari
- Department of Biosciences, Integral University, Lucknow, India
| | - Pratichi Singh
- School of Biosciences and Technology, Vellore Institute of Technology, Vellore, Tamilnadu, India
| | - Febin Prabhu Dass J
- School of Biosciences and Technology, Vellore Institute of Technology, Vellore, Tamilnadu, India
| |
Collapse
|
50
|
Rahamim G, Amir D, Haas E. Simultaneous Determination of Two Subdomain Folding Rates Using the "Transfer-Quench" Method. Biophys J 2017; 112:1786-1796. [PMID: 28494950 DOI: 10.1016/j.bpj.2017.01.037] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Revised: 12/21/2016] [Accepted: 01/06/2017] [Indexed: 11/29/2022] Open
Abstract
The investigation of the mechanism of protein folding is complicated by the context dependence of the rates of intramolecular contact formation. Methods based on site-specific labeling and ultrafast spectroscopic detection of fluorescence signals were developed for monitoring the rates of individual subdomain folding transitions in situ, in the context of the whole molecule. However, each site-specific labeling modification might affect rates of folding of near-neighbor structural elements, and thus limit the ability to resolve fine differences in rates of folding of these elements. Therefore, it is highly desirable to be able to study the rates of folding of two or more neighboring subdomain structures using a single mutant to facilitate resolution of the order and interdependence of such steps. Here, we report the development of the "Transfer-Quench" method for measuring the rate of formation of two structural elements using a single triple-labeled mutant. This method is based on Förster resonance energy transfer combined with fluorescence quenching. We placed the donor and acceptor at the loop ends, and a quencher at an α-helical element involved in the node forming the loop. The folding of the triple-labeled mutant is monitored by the acceptor emission. The formation of nonlocal contact (loop closure) increases the time-dependent acceptor emission, while the closure of the labeled helix turn reduces this emission. The method was applied in a study of the folding mechanism of the common model protein, the B domain of staphylococcal protein A. Only natural amino acids were used as probes, and thus possible structural perturbations were minimized. Tyr and Trp residues served as donor and acceptor at the ends of a long loop between helices I and II, and a Cys residue as a quencher for the acceptor. We found that the closure of the loop (segment 14-33) occurs with the same rate constant as the nucleation of helix HII (segment 33-29), in line with the nucleation-condensation model.
Collapse
Affiliation(s)
- Gil Rahamim
- The Goodman Faculty of Life Sciences Bar Ilan University, Ramat Gan, Israel
| | - Dan Amir
- The Goodman Faculty of Life Sciences Bar Ilan University, Ramat Gan, Israel
| | - Elisha Haas
- The Goodman Faculty of Life Sciences Bar Ilan University, Ramat Gan, Israel.
| |
Collapse
|