1
|
Zheng F, Jiang X, Wen Y, Yang Y, Li M. Systematic investigation of machine learning on limited data: A study on predicting protein-protein binding strength. Comput Struct Biotechnol J 2024; 23:460-472. [PMID: 38235359 PMCID: PMC10792694 DOI: 10.1016/j.csbj.2023.12.018] [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/03/2023] [Revised: 12/14/2023] [Accepted: 12/16/2023] [Indexed: 01/19/2024] Open
Abstract
The application of machine learning techniques in biological research, especially when dealing with limited data availability, poses significant challenges. In this study, we leveraged advancements in method development for predicting protein-protein binding strength to conduct a systematic investigation into the application of machine learning on limited data. The binding strength, quantitatively measured as binding affinity, is vital for understanding the processes of recognition, association, and dysfunction that occur within protein complexes. By incorporating transfer learning, integrating domain knowledge, and employing both deep learning and traditional machine learning algorithms, we mitigated the impact of data limitations and made significant advancements in predicting protein-protein binding affinity. In particular, we developed over 20 models, ultimately selecting three representative best-performing ones that belong to distinct categories. The first model is structure-based, consisting of a random forest regression and thirteen handcrafted features. The second model is sequence-based, employing an architecture that combines transferred embedding features with a multilayer perceptron. Finally, we created an ensemble model by averaging the predictions of the two aforementioned models. The comparison with other predictors on three independent datasets confirms the significant improvements achieved by our models in predicting protein-protein binding affinity. The programs for running these three models are available at https://github.com/minghuilab/BindPPI.
Collapse
Affiliation(s)
- Feifan Zheng
- MOE Key Laboratory of Geriatric Diseases and Immunology, School of Biology and Basic Medical Sciences, Suzhou Medical College of Soochow University, Suzhou, Jiangsu Province 215123, China
| | - Xin Jiang
- MOE Key Laboratory of Geriatric Diseases and Immunology, School of Biology and Basic Medical Sciences, Suzhou Medical College of Soochow University, Suzhou, Jiangsu Province 215123, China
| | - Yuhao Wen
- MOE Key Laboratory of Geriatric Diseases and Immunology, School of Biology and Basic Medical Sciences, Suzhou Medical College of Soochow University, Suzhou, Jiangsu Province 215123, China
| | - Yan Yang
- MOE Key Laboratory of Geriatric Diseases and Immunology, School of Biology and Basic Medical Sciences, Suzhou Medical College of Soochow University, Suzhou, Jiangsu Province 215123, China
| | - Minghui Li
- MOE Key Laboratory of Geriatric Diseases and Immunology, School of Biology and Basic Medical Sciences, Suzhou Medical College of Soochow University, Suzhou, Jiangsu Province 215123, China
| |
Collapse
|
2
|
Cai W, Peng S, Tian Y, Bao Y, Liu Q, Dong Y, Liang Z, Liu Q, Ren Y, Ding P, Liu J, Xu T, Li Y. Hydrophobic core evolution of major histocompatibility complex class I chain-related protein A for dramatic enhancing binding affinity. Int J Biol Macromol 2024; 271:132588. [PMID: 38788878 DOI: 10.1016/j.ijbiomac.2024.132588] [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: 02/20/2024] [Revised: 05/19/2024] [Accepted: 05/21/2024] [Indexed: 05/26/2024]
Abstract
Interface residues at sites of protein-protein interaction (PPI) are the focus for affinity optimisation. However, protein hydrophobic cores (HCs) play critical roles and shape the protein surface. We hypothesise that manipulating protein HCs can enhance PPI interaction affinities. A cell stress molecule, major histocompatibility complex class I chain-related protein A (MICA), binds to the natural killer group 2D (NKG2D) homodimer to form three molecule interactions. MICA was used as a study subject to support our hypothesis. We redesigned MICA HCs by directed mutagenesis and isolated high-affinity variants through a newly designed partial-denature panning (PDP) method. A few mutations in MICA HCs increased the NKG2D-MICA interaction affinity by 325-5613-fold. Crystal structures of the NKG2D-MICA variant complexes indicated that mutagenesis of MICA HCs stabilised helical elements for decreasing intermolecular interactive free energy (ΔG) of the NKG2D-MICA heterotrimer. The repacking of MICA HC mutants maintained overall surface residues and the authentic binding specificity of MICA. In conclusion, this study provides a new method for MICA redesign and affinity optimisation through HC manipulation without mutating PPI interface residues. Our study introduces a novel approach to protein manipulation, potentially expanding the toolkit for protein affinity optimisation.
Collapse
Affiliation(s)
- Wenxuan Cai
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; TIOC Therapeutics Limited, Hangzhou 310018, China
| | - Siqi Peng
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Ye Tian
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; TIOC Therapeutics Limited, Hangzhou 310018, China
| | - Yifeng Bao
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China
| | - Qiang Liu
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China
| | - Yan Dong
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China
| | - Zhaoduan Liang
- Bioland Laboratory, Guangzhou Regenerative Medicine and Health Guangdong Laboratory, Guangzhou 510530, China
| | - Qi Liu
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; TIOC Therapeutics Limited, Hangzhou 310018, China
| | - Yuefei Ren
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; Division of Life Sciences and Medicine, University of Science and Technology of China, Hefei 230026, China
| | - Peng Ding
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China
| | - Jinsong Liu
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; University of Chinese Academy of Sciences, Beijing 100049, China; Division of Life Sciences and Medicine, University of Science and Technology of China, Hefei 230026, China; Guangdong Provincial Key Laboratory of Biocomputing, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China
| | - Tingting Xu
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; Guangdong Provincial Key Laboratory of Biocomputing, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China.
| | - Yi Li
- State Key Laboratory of Respiratory Disease, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou 510530, China; TIOC Therapeutics Limited, Hangzhou 310018, China; University of Chinese Academy of Sciences, Beijing 100049, China; Bioland Laboratory, Guangzhou Regenerative Medicine and Health Guangdong Laboratory, Guangzhou 510530, China; Division of Life Sciences and Medicine, University of Science and Technology of China, Hefei 230026, China.
| |
Collapse
|
3
|
Nikam R, Jemimah S, Gromiha MM. DeepPPAPredMut: deep ensemble method for predicting the binding affinity change in protein-protein complexes upon mutation. Bioinformatics 2024; 40:btae309. [PMID: 38718170 PMCID: PMC11112046 DOI: 10.1093/bioinformatics/btae309] [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: 12/13/2023] [Revised: 04/08/2024] [Accepted: 05/08/2024] [Indexed: 05/24/2024] Open
Abstract
MOTIVATION Protein-protein interactions underpin many cellular processes and their disruption due to mutations can lead to diseases. With the evolution of protein structure prediction methods like AlphaFold2 and the availability of extensive experimental affinity data, there is a pressing need for updated computational tools that can efficiently predict changes in binding affinity caused by mutations in protein-protein complexes. RESULTS We developed a deep ensemble model that leverages protein sequences, predicted structure-based features, and protein functional classes to accurately predict the change in binding affinity due to mutations. The model achieved a correlation of 0.97 and a mean absolute error (MAE) of 0.35 kcal/mol on the training dataset, and maintained robust performance on the test set with a correlation of 0.72 and a MAE of 0.83 kcal/mol. Further validation using Leave-One-Out Complex (LOOC) cross-validation exhibited a correlation of 0.83 and a MAE of 0.51 kcal/mol, indicating consistent performance. AVAILABILITY AND IMPLEMENTATION https://web.iitm.ac.in/bioinfo2/DeepPPAPredMut/index.html.
Collapse
Affiliation(s)
- Rahul Nikam
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sherlyn Jemimah
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
- Department of Biomedical Engineering, Khalifa University, P.O. Box: 127788 , Abu Dhabi, United Arab Emirates
| | - M Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
- Department of Computer Science, Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsutacho, Midori-ku, Yokohama, Kanagawa 226-8501, Japan
| |
Collapse
|
4
|
Tong H, Yang T, Xu S, Li X, Liu L, Zhou G, Yang S, Yin S, Li XJ, Li S. Huntington's Disease: Complex Pathogenesis and Therapeutic Strategies. Int J Mol Sci 2024; 25:3845. [PMID: 38612657 PMCID: PMC11011923 DOI: 10.3390/ijms25073845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Revised: 03/21/2024] [Accepted: 03/26/2024] [Indexed: 04/14/2024] Open
Abstract
Huntington's disease (HD) arises from the abnormal expansion of CAG repeats in the huntingtin gene (HTT), resulting in the production of the mutant huntingtin protein (mHTT) with a polyglutamine stretch in its N-terminus. The pathogenic mechanisms underlying HD are complex and not yet fully elucidated. However, mHTT forms aggregates and accumulates abnormally in neuronal nuclei and processes, leading to disruptions in multiple cellular functions. Although there is currently no effective curative treatment for HD, significant progress has been made in developing various therapeutic strategies to treat HD. In addition to drugs targeting the neuronal toxicity of mHTT, gene therapy approaches that aim to reduce the expression of the mutant HTT gene hold great promise for effective HD therapy. This review provides an overview of current HD treatments, discusses different therapeutic strategies, and aims to facilitate future therapeutic advancements in the field.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | - Xiao-Jiang Li
- Guangdong Key Laboratory of Non-Human Primate Research, Key Laboratory of CNS Regeneration (Ministry of Education), Guangdong-Hongkong-Macau Institute of CNS Regeneration, Jinan University, Guangzhou 510632, China; (H.T.); (T.Y.); (S.X.); (X.L.); (L.L.); (G.Z.); (S.Y.); (S.Y.)
| | - Shihua Li
- Guangdong Key Laboratory of Non-Human Primate Research, Key Laboratory of CNS Regeneration (Ministry of Education), Guangdong-Hongkong-Macau Institute of CNS Regeneration, Jinan University, Guangzhou 510632, China; (H.T.); (T.Y.); (S.X.); (X.L.); (L.L.); (G.Z.); (S.Y.); (S.Y.)
| |
Collapse
|
5
|
Tandiana R, Barletta GP, Soler MA, Fortuna S, Rocchia W. Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔ G Predictions. J Chem Theory Comput 2024; 20:2630-2642. [PMID: 38445482 DOI: 10.1021/acs.jctc.3c01225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
The development of highly potent antibodies and antibody fragments as binding agents holds significant implications in fields such as biosensing and biotherapeutics. Their binding strength is intricately linked to the arrangement and composition of residues at the binding interface. Computational techniques offer a robust means to predict the three-dimensional structure of these complexes and to assess the affinity changes resulting from mutations. Given the interdependence of structure and affinity prediction, our objective here is to disentangle their roles. We aim to evaluate independently six side-chain reconstruction methods and ten binding affinity estimation techniques. This evaluation was pivotal in predicting affinity alterations due to single mutations, a key step in computational affinity maturation protocols. Our analysis focuses on a data set comprising 27 distinct antibody/hen egg white lysozyme complexes, each with crystal structures and experimentally determined binding affinities. Using six different side-chain reconstruction methods, we transformed each structure into its corresponding mutant via in silico single-point mutations. Subsequently, these structures undergo minimization and molecular dynamics simulation. We therefore estimate ΔΔG values based on the original crystal structure, its energy-minimized form, and the ensuing molecular dynamics trajectories. Our research underscores the critical importance of selecting reliable side-chain reconstruction methods and conducting thorough molecular dynamics simulations to accurately predict the impact of mutations. In summary, our study demonstrates that the integration of conformational sampling and scoring is a potent approach to precisely characterizing mutation processes in single-point mutagenesis protocols and crucial for computational antibody design.
Collapse
Affiliation(s)
- Rika Tandiana
- Computational MOdelling of NanosCalE and BioPhysical SysTems─CONCEPT Lab Istituto Italiano di Tecnologia (IIT), Via Melen-83, B Block, 16152 Genoa, Italy
| | - German P Barletta
- Computational MOdelling of NanosCalE and BioPhysical SysTems─CONCEPT Lab Istituto Italiano di Tecnologia (IIT), Via Melen-83, B Block, 16152 Genoa, Italy
- The Abdus Salam International Centre for Theoretical Physics─ICTP, Strada Costiera 11, 34151 Trieste, Italy
| | - Miguel Angel Soler
- Dipartimento di Scienze Matematiche, Informatiche e Fisiche, Universita' di Udine, Via delle Scienze 206, 33100 Udine, Italy
| | - Sara Fortuna
- Computational MOdelling of NanosCalE and BioPhysical SysTems─CONCEPT Lab Istituto Italiano di Tecnologia (IIT), Via Melen-83, B Block, 16152 Genoa, Italy
| | - Walter Rocchia
- Computational MOdelling of NanosCalE and BioPhysical SysTems─CONCEPT Lab Istituto Italiano di Tecnologia (IIT), Via Melen-83, B Block, 16152 Genoa, Italy
| |
Collapse
|
6
|
Nikam R, Yugandhar K, Gromiha MM. Deep learning-based method for predicting and classifying the binding affinity of protein-protein complexes. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2023; 1871:140948. [PMID: 37567456 DOI: 10.1016/j.bbapap.2023.140948] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2023] [Revised: 08/05/2023] [Accepted: 08/08/2023] [Indexed: 08/13/2023]
Abstract
Protein-protein interactions (PPIs) play a critical role in various biological processes. Accurately estimating the binding affinity of PPIs is essential for understanding the underlying molecular recognition mechanisms. In this study, we employed a deep learning approach to predict the binding affinity (ΔG) of protein-protein complexes. To this end, we compiled a dataset of 903 protein-protein complexes, each with its corresponding experimental binding affinity, which belong to six functional classes. We extracted 8 to 20 non-redundant features from the sequence information as well as the predicted three-dimensional structures using feature selection methods for each protein functional class. Our method showed an overall mean absolute error of 1.05 kcal/mol and a correlation of 0.79 between experimental and predicted ΔG values. Additionally, we evaluated our model for discriminating high and low affinity protein-protein complexes and it achieved an accuracy of 87% with an F1 score of 0.86 using 10-fold cross-validation on the selected features. Our approach presents an efficient tool for studying PPIs and provides crucial insights into the underlying mechanisms of the molecular recognition process. The web server can be freely accessed at https://web.iitm.ac.in/bioinfo2/DeepPPAPred/index.html.
Collapse
Affiliation(s)
- Rahul Nikam
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, Tamil Nadu, India
| | - Kumar Yugandhar
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, Tamil Nadu, India; Department of Computational Biology, Cornell University, New York, 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
|
7
|
Appadurai R, Koneru JK, Bonomi M, Robustelli P, Srivastava A. Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding. J Chem Theory Comput 2023; 19:4711-4727. [PMID: 37338049 PMCID: PMC11108026 DOI: 10.1021/acs.jctc.3c00224] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.
Collapse
Affiliation(s)
- Rajeswari Appadurai
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| | | | - Massimiliano Bonomi
- Structural Bioinformatics Unit, Department of Structural Biology and Chemistry. CNRS UMR 3528, C3BI, CNRS USR 3756, Institut Pasteur, Paris, France
| | - Paul Robustelli
- Dartmouth College, Department of Chemistry, Hanover, NH, 03755, USA
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| |
Collapse
|
8
|
Yang YX, Huang JY, Wang P, Zhu BT. AREA-AFFINITY: A Web Server for Machine Learning-Based Prediction of Protein-Protein and Antibody-Protein Antigen Binding Affinities. J Chem Inf Model 2023. [PMID: 37235532 DOI: 10.1021/acs.jcim.2c01499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Protein-Protein binding affinity reflects the binding strength between the binding partners. The prediction of protein-protein binding affinity is important for elucidating protein functions and also for designing protein-based therapeutics. The geometric characteristics such as area (both interface and surface areas) in the structure of a protein-protein complex play an important role in determining protein-protein interactions and their binding affinity. Here, we present a free web server for academic use, AREA-AFFINITY, for prediction of protein-protein or antibody-protein antigen binding affinity based on interface and surface areas in the structure of a protein-protein complex. AREA-AFFINITY implements 60 effective area-based protein-protein affinity predictive models and 37 effective area-based models specific for antibody-protein antigen binding affinity prediction developed in our recent studies. These models take into consideration the roles of interface and surface areas in binding affinity by using areas classified according to different amino acid types with different biophysical nature. The models with the best performances integrate machine learning methods such as neural network or random forest. These newly developed models have superior or comparable performance compared to the commonly used existing methods. AREA-AFFINITY is available for free at: https://affinity.cuhk.edu.cn/.
Collapse
Affiliation(s)
- Yong Xiao Yang
- Shenzhen Key Laboratory of Steroid Drug Discovery and Development, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, China
| | - Jin Yan Huang
- Shenzhen Key Laboratory of Steroid Drug Discovery and Development, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, China
| | - Pan Wang
- Shenzhen Key Laboratory of Steroid Drug Discovery and Development, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, China
| | - Bao Ting Zhu
- Shenzhen Key Laboratory of Steroid Drug Discovery and Development, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, China
- Shenzhen Bay Laboratory, Shenzhen, 518055, China
| |
Collapse
|
9
|
Li Z, Zeng M, Geng K, Lai D, Xu Z, Zhou W. Chemical Constituents and Hypoglycemic Mechanisms of Dendrobium nobile in Treatment of Type 2 Diabetic Rats by UPLC-ESI-Q-Orbitrap, Network Pharmacology and In Vivo Experimental Verification. Molecules 2023; 28:molecules28062683. [PMID: 36985655 PMCID: PMC10057382 DOI: 10.3390/molecules28062683] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/07/2023] [Accepted: 03/10/2023] [Indexed: 03/18/2023] Open
Abstract
This study aimed to systematically explore the chemical constituents of D. nobile and its hypoglycemic effect by UPLC-ESI-Q-Orbitrap, network pharmacology and in vivo experiment. The chemical constituents of D. nobile were qualitatively analyzed, and the hypoglycemic compounds were quickly identified. Network pharmacological analysis and molecular docking technique were applied to assist in the elucidation of the hypoglycemic mechanisms of D. nobile. A type 2 diabetic mellitus (T2DM) rat model was established using the HFD and STZ method for in vivo experimental verification, and these T2DM rats were treated with D. nobile extract and D. nobile polysaccharide for two months by gavage. The results showed that a total of 39 chemical constituents of D. nobile, including alkaloids, bibenzyls, phenanthrenes and other types of compounds, were identified. D. nobile extract and D. nobile polysaccharide could significantly ameliorate the body weight, hyperglycemia, insulin resistance, dyslipidemia and morphological impairment of the liver and pancreas in the T2DM rats. α-Linolenic acid, dihydroconiferyl dihydro-p-coumarate, naringenin, trans-N-feruloyltyramine, gigantol, moscatilin, 4-O-methylpinosylvic acid, venlafaxine, nordendrobin and tristin were regarded as the key hypoglycemic compounds of D. nobile, along with the hypoglycemic effect on the PI3K-AKT signaling pathway, the insulin signaling pathway, the FOXO signaling pathway, the improvement of insulin resistance and the AGE-RAGE signaling pathway. The Western blotting experiment results confirmed that D. nobile activated the PI3K/AKT pathway and insulin signaling pathway, promoted glycogen synthesis via regulating the expression of glycogen synthase kinase 3 beta (GSK-3β) and glucose transporter 4 (GLUT4), and inhibited liver gluconeogenesis by regulating the expression of phosphoenolpyruvate carboxykinase (PEPCK) and glucose 6 phosphatase (G6pase) in the liver. The results suggested that the hypoglycemic mechanism of D. nobile might be associated with liver glycogen synthesis and gluconeogenesis, contributing to improving insulin resistance and abnormal glucose metabolism in the T2DM rats.
Collapse
Affiliation(s)
- Zhaoyang Li
- School of Pharmacy, Guizhou Medical University, Guiyang 550025, China
| | - Meiling Zeng
- School of Pharmacy, Guizhou Medical University, Guiyang 550025, China
| | - Keyong Geng
- School of Pharmacy, Guizhou Medical University, Guiyang 550025, China
| | - Donna Lai
- School of Medicine, Western Sydney University, Penrith, NSW 2751, Australia
- Correspondence: (D.L.); (Z.X.); (W.Z.)
| | - Zhi Xu
- Guizhou Miaoaitang Health Management Co., Ltd., Guiyang 550025, China
- Correspondence: (D.L.); (Z.X.); (W.Z.)
| | - Wei Zhou
- School of Pharmacy, Guizhou Medical University, Guiyang 550025, China
- Correspondence: (D.L.); (Z.X.); (W.Z.)
| |
Collapse
|
10
|
Harnessing Fc/FcRn Affinity Data from Patents with Different Machine Learning Methods. Int J Mol Sci 2023; 24:ijms24065724. [PMID: 36982796 PMCID: PMC10052518 DOI: 10.3390/ijms24065724] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/07/2023] [Accepted: 03/11/2023] [Indexed: 03/19/2023] Open
Abstract
Monoclonal antibodies are biopharmaceuticals with a very long half-life due to the binding of their Fc portion to the neonatal receptor (FcRn), a pharmacokinetic property that can be further improved through engineering of the Fc portion, as demonstrated by the approval of several new drugs. Many Fc variants with increased binding to FcRn have been found using different methods, such as structure-guided design, random mutagenesis, or a combination of both, and are described in the literature as well as in patents. Our hypothesis is that this material could be subjected to a machine learning approach in order to generate new variants with similar properties. We therefore compiled 1323 Fc variants affecting the affinity for FcRn, which were disclosed in twenty patents. These data were used to train several algorithms, with two different models, in order to predict the affinity for FcRn of new randomly generated Fc variants. To determine which algorithm was the most robust, we first assessed the correlation between measured and predicted affinity in a 10-fold cross-validation test. We then generated variants by in silico random mutagenesis and compared the prediction made by the different algorithms. As a final validation, we produced variants, not described in any patents, and compared the predicted affinity with the experimental binding affinities measured by surface plasmon resonance (SPR). The best mean absolute error (MAE) between predicted and experimental values was obtained with a support vector regressor (SVR) using six features and trained on 1251 examples. With this setting, the error on the log(KD) was less than 0.17. The obtained results show that such an approach could be used to find new variants with better half-life properties that are different from those already extensively used in therapeutic antibody development.
Collapse
|
11
|
Ferraz MVF, Neto JCS, Lins RD, Teixeira ES. An artificial neural network model to predict structure-based protein-protein free energy of binding from Rosetta-calculated properties. Phys Chem Chem Phys 2023; 25:7257-7267. [PMID: 36810523 DOI: 10.1039/d2cp05644e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The prediction of the free energy (ΔG) of binding for protein-protein complexes is of general scientific interest as it has a variety of applications in the fields of molecular and chemical biology, materials science, and biotechnology. Despite its centrality in understanding protein association phenomena and protein engineering, the ΔG of binding is a daunting quantity to obtain theoretically. In this work, we devise a novel Artificial Neural Network (ANN) model to predict the ΔG of binding for a given three-dimensional structure of a protein-protein complex with Rosetta-calculated properties. Our model was tested using two data sets, and it presented a root-mean-square error ranging from 1.67 kcal mol-1 to 2.45 kcal mol-1, showing a better performance compared to the available state-of-the-art tools. Validation of the model for a variety of protein-protein complexes is showcased.
Collapse
Affiliation(s)
- Matheus V F Ferraz
- Department of Virology, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, FIOCRUZ, Recife, PE, Brazil.,Department of Fundamental Chemistry, Federal University of Pernambuco, UFPE, Recife, PE, Brazil.,Heidelberg Institute for Theoretical Studies, HITS, Heidelberg, Germany
| | - José C S Neto
- Recife Center for Advanced Studies and Systems, CESAR, Recife, PE, Brazil.
| | - Roberto D Lins
- Department of Virology, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, FIOCRUZ, Recife, PE, Brazil.,Department of Fundamental Chemistry, Federal University of Pernambuco, UFPE, Recife, PE, Brazil
| | - Erico S Teixeira
- Recife Center for Advanced Studies and Systems, CESAR, Recife, PE, Brazil.
| |
Collapse
|
12
|
Mai H, Zimmer MH, Miller TF. Exploring PROTAC Cooperativity with Coarse-Grained Alchemical Methods. J Phys Chem B 2023; 127:446-455. [PMID: 36607139 PMCID: PMC9869335 DOI: 10.1021/acs.jpcb.2c05795] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Revised: 12/18/2022] [Indexed: 01/07/2023]
Abstract
Proteolysis targeting chimera (PROTAC) is a novel drug modality that facilitates the degradation of a target protein by inducing proximity with an E3 ligase. In this work, we present a new computational framework to model the cooperativity between PROTAC-E3 binding and PROTAC-target binding principally through protein-protein interactions (PPIs) induced by the PROTAC. Due to the scarcity and low resolution of experimental measurements, the physical and chemical drivers of these non-native PPIs remain to be elucidated. We develop a coarse-grained (CG) approach to model interactions in the target-PROTAC-E3 complexes, which enables converged thermodynamic estimations using alchemical free energy calculation methods despite an unconventional scale of perturbations. With minimal parametrization, we successfully capture fundamental principles of cooperativity, including the optimality of intermediate PROTAC linker lengths that originates from configurational entropy. We qualitatively characterize the dependency of cooperativity on PROTAC linker lengths and protein charges and shapes. Minimal inclusion of sequence- and conformation-specific features in our current force field, however, limits quantitative modeling to reproduce experimental measurements, but further development of the CG model may allow for efficient computational screening to optimize PROTAC cooperativity.
Collapse
Affiliation(s)
- Huanghao Mai
- Division of Chemistry and Chemical
Engineering, California Institute of Technology, Pasadena, California91125, United States
| | - Matthew H. Zimmer
- Division of Chemistry and Chemical
Engineering, California Institute of Technology, Pasadena, California91125, United States
| | - Thomas F. Miller
- Division of Chemistry and Chemical
Engineering, California Institute of Technology, Pasadena, California91125, United States
| |
Collapse
|
13
|
Jung Y, Geng C, Bonvin AMJJ, Xue LC, Honavar VG. MetaScore: A Novel Machine-Learning-Based Approach to Improve Traditional Scoring Functions for Scoring Protein-Protein Docking Conformations. Biomolecules 2023; 13:121. [PMID: 36671507 PMCID: PMC9855734 DOI: 10.3390/biom13010121] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 12/22/2022] [Accepted: 12/26/2022] [Indexed: 01/11/2023] Open
Abstract
Protein-protein interactions play a ubiquitous role in biological function. Knowledge of the three-dimensional (3D) structures of the complexes they form is essential for understanding the structural basis of those interactions and how they orchestrate key cellular processes. Computational docking has become an indispensable alternative to the expensive and time-consuming experimental approaches for determining the 3D structures of protein complexes. Despite recent progress, identifying near-native models from a large set of conformations sampled by docking-the so-called scoring problem-still has considerable room for improvement. We present MetaScore, a new machine-learning-based approach to improve the scoring of docked conformations. MetaScore utilizes a random forest (RF) classifier trained to distinguish near-native from non-native conformations using their protein-protein interfacial features. The features include physicochemical properties, energy terms, interaction-propensity-based features, geometric properties, interface topology features, evolutionary conservation, and also scores produced by traditional scoring functions (SFs). MetaScore scores docked conformations by simply averaging the score produced by the RF classifier with that produced by any traditional SF. We demonstrate that (i) MetaScore consistently outperforms each of the nine traditional SFs included in this work in terms of success rate and hit rate evaluated over conformations ranked among the top 10; (ii) an ensemble method, MetaScore-Ensemble, that combines 10 variants of MetaScore obtained by combining the RF score with each of the traditional SFs outperforms each of the MetaScore variants. We conclude that the performance of traditional SFs can be improved upon by using machine learning to judiciously leverage protein-protein interfacial features and by using ensemble methods to combine multiple scoring functions.
Collapse
Affiliation(s)
- Yong Jung
- Bioinformatics & Genomics Graduate Program, Pennsylvania State University, University Park, PA 16802, USA
- Artificial Intelligence Research Laboratory, Pennsylvania State University, University Park, PA 16802, USA
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA 16802, USA
| | - Cunliang Geng
- Bijvoet Centre for Biomolecular Research, Faculty of Science—Chemistry, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands
| | - Alexandre M. J. J. Bonvin
- Bijvoet Centre for Biomolecular Research, Faculty of Science—Chemistry, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands
| | - Li C. Xue
- Bijvoet Centre for Biomolecular Research, Faculty of Science—Chemistry, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands
- Center for Molecular and Biomolecular Informatics, Radboudumc, Greet Grooteplein 26-28, 6525 GA Nijmegen, The Netherlands
| | - Vasant G. Honavar
- Bioinformatics & Genomics Graduate Program, Pennsylvania State University, University Park, PA 16802, USA
- Artificial Intelligence Research Laboratory, Pennsylvania State University, University Park, PA 16802, USA
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA 16802, USA
- Clinical and Translational Sciences Institute, Pennsylvania State University, University Park, PA 16802, USA
- College of Information Sciences & Technology, Pennsylvania State University, University Park, PA 16802, USA
- Institute for Computational and Data Sciences, Pennsylvania State University, University Park, PA 16802, USA
- Center for Big Data Analytics and Discovery Informatics, Pennsylvania State University, University Park, PA 16823, USA
| |
Collapse
|
14
|
Lee J, Seok C, Ham S, Chong S. Atomic-level thermodynamics analysis of the binding free energy of SARS-CoV-2 neutralizing antibodies. Proteins 2023; 91:694-704. [PMID: 36564921 PMCID: PMC9880660 DOI: 10.1002/prot.26458] [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] [Revised: 12/20/2022] [Accepted: 12/20/2022] [Indexed: 12/25/2022]
Abstract
Understanding how protein-protein binding affinity is determined from molecular interactions at the interface is essential in developing protein therapeutics such as antibodies, but this has not yet been fully achieved. Among the major difficulties are the facts that it is generally difficult to decompose thermodynamic quantities into contributions from individual molecular interactions and that the solvent effect-dehydration penalty-must also be taken into consideration for every contact formation at the binding interface. Here, we present an atomic-level thermodynamics analysis that overcomes these difficulties and illustrate its utility through application to SARS-CoV-2 neutralizing antibodies. Our analysis is based on the direct interaction energy computed from simulated antibody-protein complex structures and on the decomposition of solvation free energy change upon complex formation. We find that the formation of a single contact such as a hydrogen bond at the interface barely contributes to binding free energy due to the dehydration penalty. On the other hand, the simultaneous formation of multiple contacts between two interface residues favorably contributes to binding affinity. This is because the dehydration penalty is significantly alleviated: the total penalty for multiple contacts is smaller than a sum of what would be expected for individual dehydrations of those contacts. Our results thus provide a new perspective for designing protein therapeutics of improved binding affinity.
Collapse
Affiliation(s)
- Jihyeon Lee
- Department of ChemistrySeoul National UniversitySeoulSouth Korea
| | - Chaok Seok
- Department of ChemistrySeoul National UniversitySeoulSouth Korea
| | - Sihyun Ham
- Department of ChemistrySookmyung Women's UniversitySeoulSouth Korea
| | - Song‐Ho Chong
- Global Center for Natural Resources Sciences, Faculty of Life SciencesKumamoto UniversityKumamotoJapan
| |
Collapse
|
15
|
Harini K, Christoffer C, Gromiha MM, Kihara D. Pairwise and Multi-chain Protein Docking Enhanced Using LZerD Web Server. Methods Mol Biol 2023; 2690:355-373. [PMID: 37450159 PMCID: PMC10561630 DOI: 10.1007/978-1-0716-3327-4_28] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
Interactions of proteins with other macromolecules have important structural and functional roles in the basic processes of living cells. To understand and elucidate the mechanisms of interactions, it is important to know the 3D structures of the complexes. Proteomes contain numerous protein-protein complexes, for which experimentally determined structures often do not exist. Computational techniques can be a practical alternative to obtain useful complex structure models. Here, we present a web server that provides access to the LZerD and Multi-LZerD protein docking tools, which can perform both pairwise and multi-chain docking. The web server is user-friendly, with options to visualize the distribution and structures of binding poses of top-scoring models. The LZerD web server is available at https://lzerd.kiharalab.org . This chapter dictates the algorithm and step-by-step procedure to model the monomeric structures with AttentiveDist, and also provides the detail of pairwise LZerD docking, and multi-LZerD. This also provided case studies for each of the three modules.
Collapse
Affiliation(s)
- Kannan Harini
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, India
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA
| | | | - M Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, India
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA
| | - Daisuke Kihara
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA.
- Department of Computer Science, Purdue University, West Lafayette, IN, USA.
| |
Collapse
|
16
|
Gurusinghe SN, Oppenheimer B, Shifman JM. Cold spots are universal in protein–protein interactions. Protein Sci 2022; 31:e4435. [PMID: 36173158 PMCID: PMC9490803 DOI: 10.1002/pro.4435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2022] [Revised: 07/22/2022] [Accepted: 08/26/2022] [Indexed: 12/02/2022]
Abstract
Proteins interact with each other through binding interfaces that differ greatly in size and physico‐chemical properties. Within the binding interface, a few residues called hot spots contribute the majority of the binding free energy and are hence irreplaceable. In contrast, cold spots are occupied by suboptimal amino acids, providing possibility for affinity enhancement through mutations. In this study, we identify cold spots due to cavities and unfavorable charge interactions in multiple protein–protein interactions (PPIs). For our cold spot analysis, we first use a small affinity database of PPIs with known structures and affinities and then expand our search to nearly 4000 homo‐ and heterodimers in the Protein Data Bank (PDB). We observe that cold spots due to cavities are present in nearly all PPIs unrelated to their binding affinity, while unfavorable charge interactions are relatively rare. We also find that most cold spots are located in the periphery of the binding interface, with high‐affinity complexes showing fewer centrally located colds spots than low‐affinity complexes. A larger number of cold spots is also found in non‐cognate interactions compared to their cognate counterparts. Furthermore, our analysis reveals that cold spots are more frequent in homo‐dimeric complexes compared to hetero‐complexes, likely due to symmetry constraints imposed on sequences of homodimers. Finally, we find that glycines, glutamates, and arginines are the most frequent amino acids appearing at cold spot positions. Our analysis emphasizes the importance of cold spot positions to protein evolution and facilitates protein engineering studies directed at enhancing binding affinity and specificity in a wide range of applications.
Collapse
Affiliation(s)
- Sagara N.S. Gurusinghe
- Department of Biological Chemistry The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem Jerusalem Israel
| | - Ben Oppenheimer
- Department of Biological Chemistry The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem Jerusalem Israel
| | - Julia M. Shifman
- Department of Biological Chemistry The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem Jerusalem Israel
| |
Collapse
|
17
|
In Silico Maturation of a Nanomolar Antibody against the Human CXCR2. Biomolecules 2022; 12:biom12091285. [PMID: 36139124 PMCID: PMC9496334 DOI: 10.3390/biom12091285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2022] [Revised: 08/31/2022] [Accepted: 08/31/2022] [Indexed: 11/17/2022] Open
Abstract
The steady increase in computational power in the last 50 years is opening unprecedented opportunities in biology, as computer simulations of biological systems have become more accessible and can reproduce experimental results more accurately. Here, we wanted to test the ability of computer simulations to replace experiments in the limited but practically useful scope of improving the biochemical characteristics of the abN48 antibody, a nanomolar antagonist of the CXC chemokine receptor 2 (CXCR2) that was initially selected from a combinatorial antibody library. Our results showed a good correlation between the computed binding energies of the antibody to the peptide target and the experimental binding affinities. Moreover, we showed that it is possible to design new antibody sequences in silico with a higher affinity to the desired target using a Monte Carlo Metropolis algorithm. The newly designed sequences had an affinity comparable to the best ones obtained using in vitro affinity maturation and could be obtained within a similar timeframe. The methodology proposed here could represent a valid alternative for improving antibodies in cases in which experiments are too expensive or technically tricky and could open an opportunity for designing antibodies for targets that have been elusive so far.
Collapse
|
18
|
CSatDTA: Prediction of Drug–Target Binding Affinity Using Convolution Model with Self-Attention. Int J Mol Sci 2022; 23:ijms23158453. [PMID: 35955587 PMCID: PMC9369082 DOI: 10.3390/ijms23158453] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Revised: 07/27/2022] [Accepted: 07/27/2022] [Indexed: 12/10/2022] Open
Abstract
Drug discovery, which aids to identify potential novel treatments, entails a broad range of fields of science, including chemistry, pharmacology, and biology. In the early stages of drug development, predicting drug–target affinity is crucial. The proposed model, the prediction of drug–target affinity using a convolution model with self-attention (CSatDTA), applies convolution-based self-attention mechanisms to the molecular drug and target sequences to predict drug–target affinity (DTA) effectively, unlike previous convolution methods, which exhibit significant limitations related to this aspect. The convolutional neural network (CNN) only works on a particular region of information, excluding comprehensive details. Self-attention, on the other hand, is a relatively recent technique for capturing long-range interactions that has been used primarily in sequence modeling tasks. The results of comparative experiments show that CSatDTA surpasses previous sequence-based or other approaches and has outstanding retention abilities.
Collapse
|
19
|
Williams A, Zhan CG. Fast Prediction of Binding Affinities of SARS-CoV-2 Spike Protein and Its Mutants with Antibodies through Intermolecular Interaction Modeling-Based Machine Learning. J Phys Chem B 2022; 126:5194-5206. [PMID: 35817617 PMCID: PMC9301770 DOI: 10.1021/acs.jpcb.2c02123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 06/28/2022] [Indexed: 11/30/2022]
Abstract
Since the introduction of the novel SARS-CoV-2 virus (COVID-19) in late 2019, various new variants have appeared with mutations that confer resistance to the vaccines and monoclonal antibodies that were developed in response to the wild-type virus. As we continue through the pandemic, an accurate and efficient methodology is needed to help predict the effects certain mutations will have on both our currently produced therapeutics and those that are in development. Using published cryo-electron microscopy and X-ray crystallography structures of the spike receptor binding domain region with currently known antibodies, in the present study, we created and cross-validated an intermolecular interaction modeling-based multi-layer perceptron machine learning approach that can accurately predict the mutation-caused shifts in the binding affinity between the spike protein (wild-type or mutant) and various antibodies. This validated artificial intelligence (AI) model was used to predict the binding affinity (Kd) of reported SARS-CoV-2 antibodies with various variants of concern, including the most recently identified "Deltamicron" (or "Deltacron") variant. This AI model may be employed in the future to predict the Kd of developed novel antibody therapeutics to overcome the challenging antibody resistance issue and develop structural bases for the effects of both current and new mutants of the spike protein. In addition, the similar AI strategy and approach based on modeling of the intermolecular interactions may be useful in development of machine learning models predicting binding affinities for other protein-protein binding systems, including other antibodies binding with their antigens.
Collapse
Affiliation(s)
- Alexander
H. Williams
- Molecular
Modeling and Biopharmaceutical Center, University
of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Chang-Guo Zhan
- Molecular
Modeling and Biopharmaceutical Center, University
of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| |
Collapse
|
20
|
Desantis F, Miotto M, Di Rienzo L, Milanetti E, Ruocco G. Spatial organization of hydrophobic and charged residues affects protein thermal stability and binding affinity. Sci Rep 2022; 12:12087. [PMID: 35840609 PMCID: PMC9287411 DOI: 10.1038/s41598-022-16338-5] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2022] [Accepted: 07/08/2022] [Indexed: 11/12/2022] Open
Abstract
What are the molecular determinants of protein–protein binding affinity and whether they are similar to those regulating fold stability are two major questions of molecular biology, whose answers bring important implications both from a theoretical and applicative point of view. Here, we analyze chemical and physical features on a large dataset of protein–protein complexes with reliable experimental binding affinity data and compare them with a set of monomeric proteins for which melting temperature data was available. In particular, we probed the spatial organization of protein (1) intramolecular and intermolecular interaction energies among residues, (2) amino acidic composition, and (3) their hydropathy features. Analyzing the interaction energies, we found that strong Coulombic interactions are preferentially associated with a high protein thermal stability, while strong intermolecular van der Waals energies correlate with stronger protein–protein binding affinity. Statistical analysis of amino acids abundances, exposed to the molecular surface and/or in interaction with the molecular partner, confirmed that hydrophobic residues present on the protein surfaces are preferentially located in the binding regions, while charged residues behave oppositely. Leveraging on the important role of van der Waals interface interactions in binding affinity, we focused on the molecular surfaces in the binding regions and evaluated their shape complementarity, decomposing the molecular patches in the 2D Zernike basis. For the first time, we quantified the correlation between local shape complementarity and binding affinity via the Zernike formalism. In addition, considering the solvent interactions via the residue hydropathy, we found that the hydrophobicity of the binding regions dictates their shape complementary as much as the correlation between van der Waals energy and binding affinity. In turn, these relationships pave the way to the fast and accurate prediction and design of optimal binding regions as the 2D Zernike formalism allows a rapid and superposition-free comparison between possible binding surfaces.
Collapse
Affiliation(s)
- Fausta Desantis
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Viale Regina Elena 291, 00161, Rome, Italy.,The Open University Affiliated Research Centre at Istituto Italiano di Tecnologia, Via Morego, 30, 16163, Genoa, Italy
| | - Mattia Miotto
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Viale Regina Elena 291, 00161, Rome, Italy.
| | - Lorenzo Di Rienzo
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Viale Regina Elena 291, 00161, Rome, Italy
| | - Edoardo Milanetti
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Viale Regina Elena 291, 00161, Rome, Italy.,Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro, 5, 00185, Rome, Italy
| | - Giancarlo Ruocco
- Center for Life Nano and Neuro Science, Istituto Italiano di Tecnologia (IIT), Viale Regina Elena 291, 00161, Rome, Italy.,Department of Physics, Sapienza University of Rome, Piazzale Aldo Moro, 5, 00185, Rome, Italy
| |
Collapse
|
21
|
Fang Z, Peltz G. An automated multi-modal graph-based pipeline for mouse genetic discovery. Bioinformatics 2022; 38:3385-3394. [PMID: 35608290 PMCID: PMC9992076 DOI: 10.1093/bioinformatics/btac356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 04/18/2022] [Accepted: 05/19/2022] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Our ability to identify causative genetic factors for mouse genetic models of human diseases and biomedical traits has been limited by the difficulties associated with identifying true causative factors, which are often obscured by the many false positive genetic associations produced by a GWAS. RESULTS To accelerate the pace of genetic discovery, we developed a graph neural network (GNN)-based automated pipeline (GNNHap) that could rapidly analyze mouse genetic model data and identify high probability causal genetic factors for analyzed traits. After assessing the strength of allelic associations with the strain response pattern; this pipeline analyzes 29M published papers to assess candidate gene-phenotype relationships; and incorporates the information obtained from a protein-protein interaction network and protein sequence features into the analysis. The GNN model produces markedly improved results relative to that of a simple linear neural network. We demonstrate that GNNHap can identify novel causative genetic factors for murine models of diabetes/obesity and for cataract formation, which were validated by the phenotypes appearing in previously analyzed gene knockout mice. The diabetes/obesity results indicate how characterization of the underlying genetic architecture enables new therapies to be discovered and tested by applying 'precision medicine' principles to murine models. AVAILABILITY AND IMPLEMENTATION The GNNHap source code is freely available at https://github.com/zqfang/gnnhap, and the new version of the HBCGM program is available at https://github.com/zqfang/haplomap. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Zhuoqing Fang
- Department of Anesthesia, Pain and Perioperative Medicine, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Gary Peltz
- Department of Anesthesia, Pain and Perioperative Medicine, Stanford University School of Medicine, Stanford, CA 94305, USA
| |
Collapse
|
22
|
The impact of the suppression of highly connected protein interactions on the corona virus infection. Sci Rep 2022; 12:9188. [PMID: 35654986 PMCID: PMC9160517 DOI: 10.1038/s41598-022-13373-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2021] [Accepted: 05/09/2022] [Indexed: 11/28/2022] Open
Abstract
Several highly effective Covid-19 vaccines are in emergency use, although more-infectious coronavirus strains, could delay the end of the pandemic even further. Because of this, it is highly desirable to develop fast antiviral drug treatments to accelerate the lasting immunity against the virus. From a theoretical perspective, computational approaches are useful tools for antiviral drug development based on the data analysis of gene expression, chemical structure, molecular pathway, and protein interaction mapping. This work studies the structural stability of virus–host interactome networks based on the graphical representation of virus–host protein interactions as vertices or nodes connected by commonly shared proteins. These graphical network visualization methods are analogous to those use in the design of artificial neural networks in neuromorphic computing. In standard protein-node-based network representation, virus–host interaction merges with virus–protein and host–protein networks, introducing redundant links associated with the internal virus and host networks. On the contrary, our approach provides a direct geometrical representation of viral infection structure and allows the effective and fast detection of the structural robustness of the virus–host network through proteins removal. This method was validated by applying it to H1N1 and HIV viruses, in which we were able to pinpoint the changes in the Interactome Network produced by known vaccines. The application of this method to the SARS-CoV-2 virus–host protein interactome implies that nonstructural proteins nsp4, nsp12, nsp16, the nuclear pore membrane glycoprotein NUP210, and ubiquitin specific peptidase USP54 play a crucial role in the viral infection, and their removal may provide an efficient therapy. This method may be extended to any new mutations or other viruses for which the Interactome Network is experimentally determined. Since time is of the essence, because of the impact of more-infectious strains on controlling the spread of the virus, this method may be a useful tool for novel antiviral therapies.
Collapse
|
23
|
Kotthoff I, Kundrotas PJ, Vakser IA. Dockground
scoring benchmarks for protein docking. Proteins 2022; 90:1259-1266. [DOI: 10.1002/prot.26306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Revised: 12/06/2021] [Accepted: 01/21/2022] [Indexed: 11/05/2022]
Affiliation(s)
- Ian Kotthoff
- Computational Biology Program The University of Kansas Lawrence Kansas USA
| | | | - Ilya A. Vakser
- Computational Biology Program The University of Kansas Lawrence Kansas USA
- Department of Molecular Biosciences The University of Kansas Lawrence Kansas USA
| |
Collapse
|
24
|
De Lauro A, Di Rienzo L, Miotto M, Olimpieri PP, Milanetti E, Ruocco G. Shape Complementarity Optimization of Antibody–Antigen Interfaces: The Application to SARS-CoV-2 Spike Protein. Front Mol Biosci 2022; 9:874296. [PMID: 35669567 PMCID: PMC9163568 DOI: 10.3389/fmolb.2022.874296] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Accepted: 04/07/2022] [Indexed: 11/17/2022] Open
Abstract
Many factors influence biomolecule binding, and its assessment constitutes an elusive challenge in computational structural biology. In this aspect, the evaluation of shape complementarity at molecular interfaces is one of the main factors to be considered. We focus on the particular case of antibody–antigen complexes to quantify the complementarities occurring at molecular interfaces. We relied on a method we recently developed, which employs the 2D Zernike descriptors, to characterize the investigated regions with an ordered set of numbers summarizing the local shape properties. Collecting a structural dataset of antibody–antigen complexes, we applied this method and we statistically distinguished, in terms of shape complementarity, pairs of the interacting regions from the non-interacting ones. Thus, we set up a novel computational strategy based on in silico mutagenesis of antibody-binding site residues. We developed a Monte Carlo procedure to increase the shape complementarity between the antibody paratope and a given epitope on a target protein surface. We applied our protocol against several molecular targets in SARS-CoV-2 spike protein, known to be indispensable for viral cell invasion. We, therefore, optimized the shape of template antibodies for the interaction with such regions. As the last step of our procedure, we performed an independent molecular docking validation of the results of our Monte Carlo simulations.
Collapse
Affiliation(s)
| | - Lorenzo Di Rienzo
- Center for Life Nano & Neuro-Science, Istituto Italiano di Tecnologia, Rome, Italy
- *Correspondence: Lorenzo Di Rienzo,
| | - Mattia Miotto
- Center for Life Nano & Neuro-Science, Istituto Italiano di Tecnologia, Rome, Italy
| | | | - Edoardo Milanetti
- Center for Life Nano & Neuro-Science, Istituto Italiano di Tecnologia, Rome, Italy
- Department of Physics, Sapienza University, Rome, Italy
| | - Giancarlo Ruocco
- Center for Life Nano & Neuro-Science, Istituto Italiano di Tecnologia, Rome, Italy
- Department of Physics, Sapienza University, Rome, Italy
| |
Collapse
|
25
|
Beshnova D, Fang Y, Du M, Sun Y, Du F, Ye J, Chen ZJ, Li B. Computational approach for binding prediction of SARS-CoV-2 with neutralizing antibodies. Comput Struct Biotechnol J 2022; 20:2212-2222. [PMID: 35530743 PMCID: PMC9059344 DOI: 10.1016/j.csbj.2022.04.038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2022] [Revised: 04/27/2022] [Accepted: 04/28/2022] [Indexed: 12/03/2022] Open
Abstract
Coronavirus disease 2019 (COVID-19) caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide as a severe pandemic and caused enormous global health and economical damage. Since December 2019, more than 197 million cases have been reported, causing 4.2 million deaths. In the settings of pandemic it is an urgent necessity for the development of an effective COVID-19 treatment. While in-vitro screening of hundreds of antibodies isolated from convalescent patients is challenging due to its high cost, use of computational methods may provide an attractive solution in selecting the top candidates. Here, we developed a computational approach (SARS-AB) for binding prediction of spike protein SARS-CoV-2 with monoclonal antibodies. We validated our approach using existing structures in the protein data bank (PDB), and demonstrated its prediction power in antibody-spike protein binding prediction. We further tested its performance using antibody sequences from the literature where crystal structure is not available, and observed a high prediction accuracy (AUC = 99.6%). Finally, we demonstrated that SARS-AB can be used to design effective antibodies against novel SARS-CoV-2 mutants that might escape the current antibody protections. We believe that SARS-AB can significantly accelerate the discovery of neutralizing antibodies against SARS-CoV-2 and its mutants.
Collapse
Affiliation(s)
- Daria Beshnova
- Lyda Hill Department of Bioinformatics, UT Southwestern Medical Center, Dallas, TX 75390, USA
| | - Yan Fang
- Department of Molecular Biology, USA
| | | | - Yehui Sun
- Department of Molecular Biology, USA
| | - Fenghe Du
- Department of Molecular Biology, USA
| | - Jianfeng Ye
- Lyda Hill Department of Bioinformatics, UT Southwestern Medical Center, Dallas, TX 75390, USA
| | | | - Bo Li
- Lyda Hill Department of Bioinformatics, UT Southwestern Medical Center, Dallas, TX 75390, USA
| |
Collapse
|
26
|
Zhang J, Fan F, Liu A, Zhang C, Li Q, Zhang C, He F, Shang M. Icariin: A Potential Molecule for Treatment of Knee Osteoarthritis. Front Pharmacol 2022; 13:811808. [PMID: 35479319 PMCID: PMC9037156 DOI: 10.3389/fphar.2022.811808] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Accepted: 02/21/2022] [Indexed: 01/24/2023] Open
Abstract
Background: Knee osteoarthritis (KOA) is a degenerative disease that develops over time. Icariin (ICA) has a positive effect on KOA, although the mechanism is unknown. To investigate drug-disease connections and processes, network pharmacology is commonly used. The molecular mechanisms of ICA for the treatment of KOA were investigated using network pharmacology, molecular docking and literature research approaches in this study. Methods: We gathered KOA-related genes using the DisGeNET database, the OMIM database, and GEO microarray data. TCMSP database, Pubchem database, TTD database, SwissTargetPrediction database, and Pharmmapper database were used to gather ICA-related data. Following that, a protein-protein interaction (PPI) network was created. Using the Metascape database, we performed GO and KEGG enrichment analyses. After that, we built a targets-pathways network. Furthermore, molecular docking confirms the prediction. Finally, we looked back over the last 5 years of literature on icariin for knee osteoarthritis to see if the findings of this study were accurate. Results: core targets relevant to KOA treatment include TNF, IGF1, MMP9, PTGS2, ESR1, MMP2 and so on. The main biological process involved regulation of inflammatory response, collagen catabolic process, extracellular matrix disassembly and so on. The most likely pathways involved were the IL-17 signaling pathway, TNF signaling pathway, Estrogen signaling pathway. Conclusion: ICA may alleviate KOA by inhibiting inflammation, cartilage breakdown and extracellular matrix degradation. Our study reveals the molecular mechanism of ICA for the treatment of KOA, demonstrating its potential value for further research and as a new drug.
Collapse
Affiliation(s)
- Juntao Zhang
- Academy of Medical Engineering and Traditional Medicine, Tianjin University, Tianjin China.,Orthopedics Department, National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, The First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
| | - Fangyang Fan
- Orthopedics Department, National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, The First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
| | - Aifeng Liu
- Orthopedics Department, National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, The First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
| | - Chao Zhang
- Orthopedics Department, National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, The First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
| | - Qi Li
- Orthopedics Department, National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, The First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
| | - Chenglong Zhang
- Orthopedics Department, National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, The First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
| | - Feng He
- Academy of Medical Engineering and Traditional Medicine, Tianjin University, Tianjin China
| | - Man Shang
- Department of Pharmacology, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
| |
Collapse
|
27
|
Panday S, Alexov E. Protein-Protein Binding Free Energy Predictions with the MM/PBSA Approach Complemented with the Gaussian-Based Method for Entropy Estimation. ACS OMEGA 2022; 7:11057-11067. [PMID: 35415339 PMCID: PMC8991903 DOI: 10.1021/acsomega.1c07037] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 03/10/2022] [Indexed: 06/14/2023]
Abstract
Here, we present a Gaussian-based method for estimation of protein-protein binding entropy to augment the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method for computational prediction of binding free energy (ΔG). The method is termed f5-MM/PBSA/E, where "E" stands for entropy and f5 for five adjustable parameters. The enthalpy components of ΔG (molecular mechanics, polar and non-polar solvation energies) are computed from a single implicit solvent generalized Born (GB) energy minimized structure of a protein-protein complex, while the binding entropy is computed using independently GB energy minimized unbound and bound structures. It should be emphasized that the f5-MM/PBSA/E method does not use snapshots, just energy minimized structures, and is thus very fast and computationally efficient. The method is trained and benchmarked in 5-fold validation test over a data set consisting of 46 protein-protein binding cases with experimentally determined dissociation constant K d values. This data set has been used for benchmarking in recently published protein-protein binding studies that apply conventional MM/PBSA and MM/PBSA with an enhanced sampling method. The f5-MM/PBSA/E tested on the same data set achieves similar or better performance than these computationally demanding approaches, making it an excellent choice for high throughput protein-protein binding affinity prediction studies.
Collapse
|
28
|
Williams AH, Zhan CG. Generalized Methodology for the Quick Prediction of Variant SARS-CoV-2 Spike Protein Binding Affinities with Human Angiotensin-Converting Enzyme II. J Phys Chem B 2022; 126:2353-2360. [PMID: 35315274 PMCID: PMC8982491 DOI: 10.1021/acs.jpcb.1c10718] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 02/24/2022] [Indexed: 12/25/2022]
Abstract
Variants of the SARS-CoV-2 virus continue to remain a threat 2 years from the beginning of the pandemic. As more variants arise, and the B.1.1.529 (Omicron) variant threatens to create another wave of infections, a method is needed to predict the binding affinity of the spike protein quickly and accurately with human angiotensin-converting enzyme II (ACE2). We present an accurate and convenient energy minimization/molecular mechanics Poisson-Boltzmann surface area methodology previously used with engineered ACE2 therapeutics to predict the binding affinity of the Omicron variant. Without any additional data from the variants discovered after the publication of our first model, the methodology can accurately predict the binding of the spike/ACE2 variant complexes. From this methodology, we predicted that the Omicron variant spike has a Kd of ∼22.69 nM (which is very close to the experimental Kd of 20.63 nM published during the review process of the current report) and that spike protein of the new "Stealth" Omicron variant (BA.2) will display a Kd of ∼12.9 nM with the wild-type ACE2 protein. This methodology can be used with as-yet discovered variants, allowing for quick determinations regarding the variant's infectivity versus either the wild-type virus or its variants.
Collapse
Affiliation(s)
- Alexander H. Williams
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
29
|
Sarnowski CP, Bikaki M, Leitner A. Cross-linking and mass spectrometry as a tool for studying the structural biology of ribonucleoproteins. Structure 2022; 30:441-461. [PMID: 35366400 DOI: 10.1016/j.str.2022.03.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2021] [Revised: 02/03/2022] [Accepted: 03/01/2022] [Indexed: 11/17/2022]
Abstract
Cross-linking and mass spectrometry (XL-MS) workflows represent an increasingly popular technique for low-resolution structural studies of macromolecular complexes. Cross-linking reactions take place in the solution state, capturing contact sites between components of a complex that represent the native, functionally relevant structure. Protein-protein XL-MS protocols are widely adopted, providing precise localization of cross-linking sites to single amino acid positions within a pair of cross-linked peptides. In contrast, protein-RNA XL-MS workflows are evolving rapidly and differ in their ability to localize interaction regions within the RNA sequence. Here, we review protein-protein and protein-RNA XL-MS workflows, and discuss their applications in studies of protein-RNA complexes. The examples highlight the complementary value of XL-MS in structural studies of protein-RNA complexes, where more established high-resolution techniques might be unable to produce conclusive data.
Collapse
Affiliation(s)
- Chris P Sarnowski
- Institute of Molecular Systems Biology, Department of Biology, ETH Zürich, 8093 Zurich, Switzerland; Systems Biology PhD Program, University of Zürich and ETH Zürich, Zurich, Switzerland
| | - Maria Bikaki
- Institute of Molecular Systems Biology, Department of Biology, ETH Zürich, 8093 Zurich, Switzerland
| | - Alexander Leitner
- Institute of Molecular Systems Biology, Department of Biology, ETH Zürich, 8093 Zurich, Switzerland.
| |
Collapse
|
30
|
PCA-MutPred: Prediction of binding free energy change upon missense mutation in protein-carbohydrate complexes. J Mol Biol 2022; 434:167526. [DOI: 10.1016/j.jmb.2022.167526] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 02/26/2022] [Accepted: 03/01/2022] [Indexed: 11/22/2022]
|
31
|
Yang YX, Wang P, Zhu BT. Relative importance of interface and surface areas in protein-protein binding affinity prediction: A machine learning analysis based on linear regression and artificial neural network. Biophys Chem 2022; 283:106762. [DOI: 10.1016/j.bpc.2022.106762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Revised: 01/11/2022] [Accepted: 01/14/2022] [Indexed: 11/02/2022]
|
32
|
Zacharias M. Match_Motif: A rapid computational tool to assist in protein-protein interaction design. Protein Sci 2022; 31:147-157. [PMID: 34648221 PMCID: PMC8740833 DOI: 10.1002/pro.4208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Revised: 10/06/2021] [Accepted: 10/12/2021] [Indexed: 11/12/2022]
Abstract
In order to generate protein assemblies with a desired function, the rational design of protein-protein binding interfaces is of significant interest. Approaches based on random mutagenesis or directed evolution may involve complex experimental selection procedures. Also, molecular modeling approaches to design entirely new proteins and interactions with partner molecules can involve large computational efforts and screening steps. In order to simplify at least the initial effort for designing a putative binding interface between two proteins the Match_Motif approach has been developed. It employs the large collection of known protein-protein complex structures to suggest interface modifications that may lead to improved binding for a desired input interaction geometry. The approach extracts interaction motifs based on the backbone structure of short (four residues) segments and the relative arrangement with respect to short segments on the partner protein. The interaction geometry is used to search through a database of such motifs in known stable bound complexes. All matches are rapidly identified (within a few seconds) and collected and can be used to guide changes in the interface that may lead to improved binding. In the output, an alternative interface structure is also proposed based on the frequency of occurrence of side chains at a given interface position in all matches and based on sterical considerations. Applications of the procedure to known complex structures and alternative arrangements are presented and discussed. The program, data files, and example applications can be downloaded from https://www.groups.ph.tum.de/t38/downloads/.
Collapse
Affiliation(s)
- Martin Zacharias
- Center of Functional Protein AssembliesTechnical University of MunichGarchingGermany
| |
Collapse
|
33
|
Renaud N, Geng C, Georgievska S, Ambrosetti F, Ridder L, Marzella DF, Réau MF, Bonvin AMJJ, Xue LC. DeepRank: a deep learning framework for data mining 3D protein-protein interfaces. Nat Commun 2021; 12:7068. [PMID: 34862392 PMCID: PMC8642403 DOI: 10.1038/s41467-021-27396-0] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Accepted: 11/12/2021] [Indexed: 11/08/2022] Open
Abstract
Three-dimensional (3D) structures of protein complexes provide fundamental information to decipher biological processes at the molecular scale. The vast amount of experimentally and computationally resolved protein-protein interfaces (PPIs) offers the possibility of training deep learning models to aid the predictions of their biological relevance. We present here DeepRank, a general, configurable deep learning framework for data mining PPIs using 3D convolutional neural networks (CNNs). DeepRank maps features of PPIs onto 3D grids and trains a user-specified CNN on these 3D grids. DeepRank allows for efficient training of 3D CNNs with data sets containing millions of PPIs and supports both classification and regression. We demonstrate the performance of DeepRank on two distinct challenges: The classification of biological versus crystallographic PPIs, and the ranking of docking models. For both problems DeepRank is competitive with, or outperforms, state-of-the-art methods, demonstrating the versatility of the framework for research in structural biology.
Collapse
Affiliation(s)
- Nicolas Renaud
- Netherlands eScience Center, Science Park 140, 1098 XG, Amsterdam, The Netherlands
| | - Cunliang Geng
- Netherlands eScience Center, Science Park 140, 1098 XG, Amsterdam, The Netherlands
- Bijvoet Centre for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, 3584, Utrecht, CH, The Netherlands
| | - Sonja Georgievska
- Netherlands eScience Center, Science Park 140, 1098 XG, Amsterdam, The Netherlands
| | - Francesco Ambrosetti
- Bijvoet Centre for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, 3584, Utrecht, CH, The Netherlands
| | - Lars Ridder
- Netherlands eScience Center, Science Park 140, 1098 XG, Amsterdam, The Netherlands
| | - Dario F Marzella
- Center for Molecular and Biomolecular Informatics, Radboudumc, Greet Grooteplein 26-28, 6525, Nijmegen, GA, The Netherlands
| | - Manon F Réau
- Bijvoet Centre for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, 3584, Utrecht, CH, The Netherlands
| | - Alexandre M J J Bonvin
- Bijvoet Centre for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, 3584, Utrecht, CH, The Netherlands.
| | - Li C Xue
- Bijvoet Centre for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, 3584, Utrecht, CH, The Netherlands.
- Center for Molecular and Biomolecular Informatics, Radboudumc, Greet Grooteplein 26-28, 6525, Nijmegen, GA, The Netherlands.
| |
Collapse
|
34
|
Wu N, Jiao L, Bütikofer M, Zeng Z, Zenobi R. High-Mass Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry for Absolute Quantitation of Noncovalent Protein-Protein Binding Interactions. Anal Chem 2021; 93:10982-10989. [PMID: 34328720 DOI: 10.1021/acs.analchem.1c02126] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) is a robust and powerful tool for studying biomacromolecules and their interactions. However, quantitative detection of high-mass analytes (kDa to MDa range) remains challenging for MALDI-MS. Herein, we successfully used commercially available purified proteins (β-galactosidase and BSA) as internal standards for high-mass MALDI-MS analysis and achieved absolute quantification of several high-mass analytes. We systematically evaluated four sample deposition methods, and using the sandwich deposition method with saturated sinapinic acid as the top layer, we performed a robust quantitative analysis by high-mass MALDI-MS. Combined with chemical cross-linking, this quantitative strategy was further used to evaluate the affinity of protein-protein interactions (PPIs), specifically of two soluble protein receptors (interleukin 1 receptor and interleukin 2 receptor) and two membrane protein receptors (rhodopsin and angiotensin 2 receptor 1) with their interaction partners. The measured dissociation constants of the protein complexes formed were between 10 nM and 5 μM. We expect this high-throughput, rapid method, which does not require labeling or immobilization of any of the interaction partners, to become a viable alternative to traditional biophysical methods for studying PPIs.
Collapse
Affiliation(s)
- Na Wu
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich CH-8093, Switzerland
| | - Lingyi Jiao
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich CH-8093, Switzerland
| | - Matthias Bütikofer
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich CH-8093, Switzerland
| | - Zhihui Zeng
- School of Materials Science and Engineering, Shandong University, Jinan 250061, P.R. China.,Empa, Swiss Federal Laboratories for Materials Science and Technology, Überlandstrasse 129, Dübendorf CH-8600, Switzerland
| | - Renato Zenobi
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich CH-8093, Switzerland
| |
Collapse
|
35
|
Kurcinski M, Kmiecik S, Zalewski M, Kolinski A. Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations. Int J Mol Sci 2021; 22:ijms22147341. [PMID: 34298961 PMCID: PMC8306105 DOI: 10.3390/ijms22147341] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 07/03/2021] [Accepted: 07/04/2021] [Indexed: 12/21/2022] Open
Abstract
Most of the protein–protein docking methods treat proteins as almost rigid objects. Only the side-chains flexibility is usually taken into account. The few approaches enabling docking with a flexible backbone typically work in two steps, in which the search for protein–protein orientations and structure flexibility are simulated separately. In this work, we propose a new straightforward approach for docking sampling. It consists of a single simulation step during which a protein undergoes large-scale backbone rearrangements, rotations, and translations. Simultaneously, the other protein exhibits small backbone fluctuations. Such extensive sampling was possible using the CABS coarse-grained protein model and Replica Exchange Monte Carlo dynamics at a reasonable computational cost. In our proof-of-concept simulations of 62 protein–protein complexes, we obtained acceptable quality models for a significant number of cases.
Collapse
|
36
|
Mishra SK, Cooper CJ, Parks JM, Mitchell JC. Hotspot Coevolution Is a Key Identifier of Near-Native Protein Complexes. J Phys Chem B 2021; 125:6058-6067. [PMID: 34077660 DOI: 10.1021/acs.jpcb.0c11525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein-protein interactions play a key role in mediating numerous biological functions, with more than half the proteins in living organisms existing as either homo- or hetero-oligomeric assemblies. Protein subunits that form oligomers minimize the free energy of the complex, but exhaustive computational search-based docking methods have not comprehensively addressed the challenge of distinguishing a natively bound complex from non-native forms. Current protein docking approaches address this problem by sampling multiple binding modes in proteins and scoring each mode, with the lowest-energy (or highest scoring) binding mode being regarded as a near-native complex. However, high-scoring modes often match poorly with the true bound form, suggesting a need for improvement of the scoring function. In this study, we propose a scoring function, KFC-E, that accounts for both conservation and coevolution of putative binding hotspot residues at protein-protein interfaces. We tested KFC-E on four benchmark sets of unbound examples and two benchmark sets of bound examples, with the results demonstrating a clear improvement over scores that examine conservation and coevolution across the entire interface.
Collapse
Affiliation(s)
- Sambit K Mishra
- Biosciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37831-6038, United States
| | - Connor J Cooper
- Biosciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37831-6038, United States
| | - Jerry M Parks
- Biosciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37831-6038, United States
| | - Julie C Mitchell
- Biosciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37831-6038, United States
| |
Collapse
|
37
|
Parate S, Rampogu S, Lee G, Hong JC, Lee KW. Exploring the Binding Interaction of Raf Kinase Inhibitory Protein With the N-Terminal of C-Raf Through Molecular Docking and Molecular Dynamics Simulation. Front Mol Biosci 2021; 8:655035. [PMID: 34124147 PMCID: PMC8194344 DOI: 10.3389/fmolb.2021.655035] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Accepted: 05/14/2021] [Indexed: 12/12/2022] Open
Abstract
Protein-protein interactions are indispensable physiological processes regulating several biological functions. Despite the availability of structural information on protein-protein complexes, deciphering their complex topology remains an outstanding challenge. Raf kinase inhibitory protein (RKIP) has gained substantial attention as a favorable molecular target for numerous pathologies including cancer and Alzheimer’s disease. RKIP interferes with the RAF/MEK/ERK signaling cascade by endogenously binding with C-Raf (Raf-1 kinase) and preventing its activation. In the current investigation, the binding of RKIP with C-Raf was explored by knowledge-based protein-protein docking web-servers including HADDOCK and ZDOCK and a consensus binding mode of C-Raf/RKIP structural complex was obtained. Molecular dynamics (MD) simulations were further performed in an explicit solvent to sample the conformations for when RKIP binds to C-Raf. Some of the conserved interface residues were mutated to alanine, phenylalanine and leucine and the impact of mutations was estimated by additional MD simulations and MM/PBSA analysis for the wild-type (WT) and constructed mutant complexes. Substantial decrease in binding free energy was observed for the mutant complexes as compared to the binding free energy of WT C-Raf/RKIP structural complex. Furthermore, a considerable increase in average backbone root mean square deviation and fluctuation was perceived for the mutant complexes. Moreover, per-residue energy contribution analysis of the equilibrated simulation trajectory by HawkDock and ANCHOR web-servers was conducted to characterize the key residues for the complex formation. One residue each from C-Raf (Arg398) and RKIP (Lys80) were identified as the druggable “hot spots” constituting the core of the binding interface and corroborated by additional long-time scale (300 ns) MD simulation of Arg398Ala mutant complex. A notable conformational change in Arg398Ala mutant occurred near the mutation site as compared to the equilibrated C-Raf/RKIP native state conformation and an essential hydrogen bonding interaction was lost. The thirteen binding sites assimilated from the overall analysis were mapped onto the complex as surface and divided into active and allosteric binding sites, depending on their location at the interface. The acquired information on the predicted 3D structural complex and the detected sites aid as promising targets in designing novel inhibitors to block the C-Raf/RKIP interaction.
Collapse
Affiliation(s)
- Shraddha Parate
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), Jinju, Korea
| | - Shailima Rampogu
- Division of Life Sciences, Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Department of Bio and Medical Big Data (BK21 Four Program), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), Jinju, Korea
| | - Gihwan Lee
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), Jinju, Korea
| | - Jong Chan Hong
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), Jinju, Korea
| | - Keun Woo Lee
- Division of Life Sciences, Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Department of Bio and Medical Big Data (BK21 Four Program), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), Jinju, Korea
| |
Collapse
|
38
|
Williams AH, Zhan CG. Fast Prediction of Binding Affinities of the SARS-CoV-2 Spike Protein Mutant N501Y (UK Variant) with ACE2 and Miniprotein Drug Candidates. J Phys Chem B 2021; 125:4330-4336. [PMID: 33881861 PMCID: PMC8084269 DOI: 10.1021/acs.jpcb.1c00869] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Revised: 03/30/2021] [Indexed: 12/28/2022]
Abstract
A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (Kd) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.
Collapse
Affiliation(s)
- Alexander H. Williams
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
39
|
Patel D, Patel JS, Ytreberg FM. Implementing and Assessing an Alchemical Method for Calculating Protein-Protein Binding Free Energy. J Chem Theory Comput 2021; 17:2457-2464. [PMID: 33709712 PMCID: PMC8044032 DOI: 10.1021/acs.jctc.0c01045] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Protein-protein binding is fundamental to most biological processes. It is important to be able to use computation to accurately estimate the change in protein-protein binding free energy due to mutations in order to answer biological questions that would be experimentally challenging, laborious, or time-consuming. Although nonrigorous free-energy methods are faster, rigorous alchemical molecular dynamics-based methods are considerably more accurate and are becoming more feasible with the advancement of computer hardware and molecular simulation software. Even with sufficient computational resources, there are still major challenges to using alchemical free-energy methods for protein-protein complexes, such as generating hybrid structures and topologies, maintaining a neutral net charge of the system when there is a charge-changing mutation, and setting up the simulation. In the current study, we have used the pmx package to generate hybrid structures and topologies, and a double-system/single-box approach to maintain the net charge of the system. To test the approach, we predicted relative binding affinities for two protein-protein complexes using a nonequilibrium alchemical method based on the Crooks fluctuation theorem and compared the results with experimental values. The method correctly identified stabilizing from destabilizing mutations for a small protein-protein complex, and a larger, more challenging antibody complex. Strong correlations were obtained between predicted and experimental relative binding affinities for both protein-protein systems.
Collapse
Affiliation(s)
- Dharmeshkumar Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho 83844, United States
| | - Jagdish Suresh Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho 83844, United States
- Department of Biological Sciences, University of Idaho, Moscow, Idaho 83844, United States
| | - F Marty Ytreberg
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho 83844, United States
- Department of Physics, University of Idaho, Moscow, Idaho 83844, United States
| |
Collapse
|
40
|
Sepahvandi A, Ghaffari M, Bahmanpour AH, Moztarzadeh F, Zarrintaj P, Uludağ H, Mozafari M. COVID-19: insights into virus-receptor interactions. MOLECULAR BIOMEDICINE 2021; 2:10. [PMID: 34766003 PMCID: PMC8035060 DOI: 10.1186/s43556-021-00033-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2020] [Accepted: 03/03/2021] [Indexed: 01/03/2023] Open
Abstract
The recent outbreak of Coronavirus Disease 2019 (COVID-19) calls for rapid mobilization of scientists to probe and explore solutions to this deadly disease. A limited understanding of the high transmissibility of SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) relative to other coronavirus strains guides a deeper investigation into the virus/receptor interactions. The cutting-edge studies in thermodynamic and kinetic properties of interactions such as protein-protein interplays have been reviewed in many modeling and analysis studies. Highlighting the thermodynamic assessments of biological interactions and emphasizing the boosted transmissibility of SARS-CoV-2 despite its high similarity in structure and sequence with other coronavirus strains is an important and highly valuable investigation that can lead scientists to discover analytical and fundamental approaches in studying virus's interactions. Accordingly, we have attempted to describe the crucial factors such as conformational changes and hydrophobicity particularities that influence on thermodynamic potentials in the SARS-COV-2 S-protein adsorption process. Discussing the thermodynamic potentials and the kinetics of the SARS-CoV-2 S-protein in its interaction with the ACE2 receptors of the host cell is a fundamental approach that would be extremely valuable in designing candidate pharmaceutical agents or exploring alternative treatments.
Collapse
Affiliation(s)
- Azadeh Sepahvandi
- Department of Mechanical Engineering College of Engineering and Computing, University of South Carolina, 301 Main St, Columbia, SC 29208 USA
| | - Maryam Ghaffari
- Biomaterial Group, Faculty of Biomedical Engineering (Center of Excellence), Amirkabir University of Technology, Tehran, Iran
| | - Amir Hossein Bahmanpour
- Biomaterial Group, Faculty of Biomedical Engineering (Center of Excellence), Amirkabir University of Technology, Tehran, Iran
| | - Fathollah Moztarzadeh
- Biomaterial Group, Faculty of Biomedical Engineering (Center of Excellence), Amirkabir University of Technology, Tehran, Iran
| | - Payam Zarrintaj
- School of Chemical Engineering, Oklahoma State University, 420 Engineering North, Stillwater, OK 74078 USA
| | - Hasan Uludağ
- Department of Chemical and Material Engineering, Faculty of Engineering, University of Alberta, Edmonton, AB T6G 2V4 Canada
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, AB T6G 2E1 Canada
- Department of Biomedical Engineering, Faculty of Medicine and Dentistry, University of Alberta, Edmonton, AB T6G 2R3 Canada
| | - Masoud Mozafari
- Department of Tissue Engineering & Regenerative Medicine, Faculty of Advanced Technologies in Medicine, Iran University of Medical Sciences (IUMS), Tehran, Iran
- Currently at: Lunenfeld-Tanenbaum Research Institute, Mount Sinai Hospital, University of Toronto, Toronto, ON Canada
| |
Collapse
|
41
|
Lee W, Park JW, Go YJ, Kim WJ, Rhee YM. Considering both small and large scale motions of vascular endothelial growth factor (VEGF) is crucial for reliably predicting its binding affinities to DNA aptamers. RSC Adv 2021; 11:9315-9326. [PMID: 35423456 PMCID: PMC8695334 DOI: 10.1039/d0ra10106k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Accepted: 02/23/2021] [Indexed: 11/21/2022] Open
Abstract
Considering both small and large scale motions of VEGF is crucial to predict its relative binding affinities to DNA aptamer variants with docking.
Collapse
Affiliation(s)
- Wook Lee
- Department of Chemistry
- Pohang University of Science and Technology (POSTECH)
- Pohang 37673
- Korea
- Department of Chemistry
| | - Jae Whee Park
- Department of Chemistry
- Korea Advanced Institute of Science and Technology (KAIST)
- Daejeon 34141
- Korea
| | - Yeon Ju Go
- Department of Chemistry
- Korea Advanced Institute of Science and Technology (KAIST)
- Daejeon 34141
- Korea
| | - Won Jong Kim
- Department of Chemistry
- Pohang University of Science and Technology (POSTECH)
- Pohang 37673
- Korea
| | - Young Min Rhee
- Department of Chemistry
- Korea Advanced Institute of Science and Technology (KAIST)
- Daejeon 34141
- Korea
| |
Collapse
|
42
|
Massoud TF, Paulmurugan R. Molecular Imaging of Protein–Protein Interactions and Protein Folding. Mol Imaging 2021. [DOI: 10.1016/b978-0-12-816386-3.00071-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
|
43
|
Gonzalez TR, Martin KP, Barnes JE, Patel JS, Ytreberg FM. Assessment of software methods for estimating protein-protein relative binding affinities. PLoS One 2020; 15:e0240573. [PMID: 33347442 PMCID: PMC7751979 DOI: 10.1371/journal.pone.0240573] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2020] [Accepted: 12/07/2020] [Indexed: 11/19/2022] Open
Abstract
A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83-98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.
Collapse
Affiliation(s)
- Tawny R. Gonzalez
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
| | - Kyle P. Martin
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| | - Jonathan E. Barnes
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| | - Jagdish Suresh Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Biological Sciences, University of Idaho, Moscow, Idaho, United States of America
| | - F. Marty Ytreberg
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| |
Collapse
|
44
|
Siebenmorgen T, Zacharias M. Efficient Refinement and Free Energy Scoring of Predicted Protein-Protein Complexes Using Replica Exchange with Repulsive Scaling. J Chem Inf Model 2020; 60:5552-5562. [PMID: 33075222 DOI: 10.1021/acs.jcim.0c00853] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate prediction and evaluation of protein-protein complex structures are of major importance to understand the cellular interactome. Typically, putative complexes are predicted based on docking methods, and simple force field or knowledge-based scoring functions are applied to evaluate single complex structures. We have extended a replica-exchange-based scheme employing different levels of a repulsive biasing between partners in each replica simulation (RS-REMD) to simultaneously refine and score protein-protein complexes. The bias acts specifically on the intermolecular interactions based on an increase in effective pairwise van der Waals radii (repulsive scaling (RS)-REMD) without affecting interactions within each protein or with the solvent. The method provides a free energy score that correlates quite well with experimental binding free energies on a set of 36 complexes with correlation coefficients of 0.77 and 0.55 in explicit and implicit solvent simulations, respectively. For a large set of docked decoy complexes, significant improvement of docked complexes was found in many cases with the starting structure in the vicinity (within 20 Å) of the native complex. In the majority of cases (14 out of 20 in explicit solvent), near native docking solutions were identified as the best scoring complexes. The approach is computational demanding but may offer a route for refinement and realistic ranking of predicted protein-protein docking geometries.
Collapse
Affiliation(s)
- Till Siebenmorgen
- Physik-Department T38, Technische Universität München, James-Franck-Str. 1, 85748 Garching, Germany
| | - Martin Zacharias
- Physik-Department T38, Technische Universität München, James-Franck-Str. 1, 85748 Garching, Germany
| |
Collapse
|
45
|
Predicting antibody affinity changes upon mutations by combining multiple predictors. Sci Rep 2020; 10:19533. [PMID: 33177627 PMCID: PMC7658247 DOI: 10.1038/s41598-020-76369-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 10/19/2020] [Indexed: 02/07/2023] Open
Abstract
Antibodies are proteins working in our immune system with high affinity and specificity for target antigens, making them excellent tools for both biotherapeutic and bioengineering applications. The prediction of antibody affinity changes upon mutations (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\Delta \Delta {\mathrm{G}}}}_{\mathrm{binding}}$$\end{document}ΔΔGbinding) is important for antibody engineering. Numerous computational methods have been proposed based on different approaches including molecular mechanics and machine learning. However, the accuracy by each individual predictor is not enough for efficient antibody development. In this study, we develop a new prediction method by combining multiple predictors based on machine learning. Our method was tested on the SiPMAB database, evaluating the Pearson’s correlation coefficient between predicted and experimental \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\Delta \Delta {\mathrm{G}}}}_{\mathrm{binding}}$$\end{document}ΔΔGbinding. Our method achieved higher accuracy (R = 0.69) than previous molecular mechanics or machine-learning based methods (R = 0.59) and the previous method using the average of multiple predictors (R = 0.64). Feature importance analysis indicated that the improved accuracy was obtained by combining predictors with different importance, which have different protocols for calculating energies and for generating mutant and unbound state structures. This study demonstrates that machine learning is a powerful framework for combining different approaches to predict antibody affinity changes.
Collapse
|
46
|
Zhou Y, Zhao R, Schwarz EC, Akbar R, Kaba M, Pattu V, Helms V, Rieger H, Nunes-Hasler P, Qu B. Interorganelle Tethering to Endocytic Organelles Determines Directional Cytokine Transport in CD4 + T Cells. THE JOURNAL OF IMMUNOLOGY 2020; 205:2988-3000. [PMID: 33106338 DOI: 10.4049/jimmunol.2000195] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Accepted: 09/20/2020] [Indexed: 12/24/2022]
Abstract
Delivery of vesicles to their desired destinations plays a central role in maintaining proper cell functionality. In certain scenarios, depending on loaded cargos, the vesicles have spatially distinct destinations. For example, in T cells, some cytokines (e.g., IL-2) are polarized to the T cell-target cell interface, whereas the other cytokines are delivered multidirectionally (e.g., TNF-α). In this study, we show that in primary human CD4+ T cells, both TNF-α+ and IL-2+ vesicles can tether with endocytic organelles (lysosomes/late endosomes) by forming membrane contact sites. Tethered cytokine-containing vesicle (CytV)-endocytic organelle pairs are released sequentially. Only endocytic organelle-tethered CytVs are preferentially transported to their desired destination. Mathematical models suggest that endocytic organelle tethering could regulate the direction of cytokine transport by selectively attaching different microtubule motor proteins (such as kinesin and dynein) to the corresponding CytVs. These findings establish the previously unknown interorganelle tethering to endocytic organelles as a universal solution for directional cytokine transport in CD4+ T cells. Modulating tethering to endocytic organelles can, therefore, coordinately control directionally distinct cytokine transport.
Collapse
Affiliation(s)
- Yan Zhou
- Department of Biophysics, Center for Integrative Physiology and Molecular Medicine, School of Medicine, Saarland University, 66421 Homburg, Germany
| | - Renping Zhao
- Department of Biophysics, Center for Integrative Physiology and Molecular Medicine, School of Medicine, Saarland University, 66421 Homburg, Germany
| | - Eva C Schwarz
- Department of Biophysics, Center for Integrative Physiology and Molecular Medicine, School of Medicine, Saarland University, 66421 Homburg, Germany
| | - Rahmad Akbar
- Center for Bioinformatics, Saarland University, 66123 Saarbrücken, Germany
| | - Mayis Kaba
- Department of Cell Physiology and Metabolism, University Medical Center, University of Geneva, 1211 Geneva, Switzerland
| | - Varsha Pattu
- Department of Physiology, Center for Integrative Physiology and Molecular Medicine, School of Medicine, Saarland University, 66421 Homburg, Germany
| | - Volkhard Helms
- Center for Bioinformatics, Saarland University, 66123 Saarbrücken, Germany
| | - Heiko Rieger
- Department of Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany
| | - Paula Nunes-Hasler
- Department of Cell Physiology and Metabolism, University Medical Center, University of Geneva, 1211 Geneva, Switzerland.,Department of Pathology and Immunology, University Medical Center, University of Geneva, 1211 Geneva, Switzerland; and
| | - Bin Qu
- Department of Biophysics, Center for Integrative Physiology and Molecular Medicine, School of Medicine, Saarland University, 66421 Homburg, Germany; .,Leibniz Institute for New Materials, 66123 Saarbrücken, Germany
| |
Collapse
|
47
|
Di Rienzo L, Milanetti E, Testi C, Montemiglio LC, Baiocco P, Boffi A, Ruocco G. A novel strategy for molecular interfaces optimization: The case of Ferritin-Transferrin receptor interaction. Comput Struct Biotechnol J 2020; 18:2678-2686. [PMID: 33101606 PMCID: PMC7548301 DOI: 10.1016/j.csbj.2020.09.020] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Revised: 09/10/2020] [Accepted: 09/11/2020] [Indexed: 11/24/2022] Open
Abstract
Protein-protein interactions regulate almost all cellular functions and rely on a fine tune of surface amino acids properties involved on both molecular partners. The disruption of a molecular association can be caused even by a single residue mutation, often leading to a pathological modification of a biochemical pathway. Therefore the evaluation of the effects of amino acid substitutions on binding, and the ad hoc design of protein-protein interfaces, is one of the biggest challenges in computational biology. Here, we present a novel strategy for computational mutation and optimization of protein-protein interfaces. Modeling the interaction surface properties using the Zernike polynomials, we describe the shape and electrostatics of binding sites with an ordered set of descriptors, making possible the evaluation of complementarity between interacting surfaces. With a Monte Carlo approach, we obtain protein mutants with controlled molecular complementarities. Applying this strategy to the relevant case of the interaction between Ferritin and Transferrin Receptor, we obtain a set of Ferritin mutants with increased or decreased complementarity. The extensive molecular dynamics validation of the method results confirms its efficacy, showing that this strategy represents a very promising approach in designing correct molecular interfaces.
Collapse
Affiliation(s)
- Lorenzo Di Rienzo
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | - Edoardo Milanetti
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Claudia Testi
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
| | | | - Paola Baiocco
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
- Department of Biochemical Sciences ‘A. Rossi Fanelli’ Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Alberto Boffi
- Department of Biochemical Sciences ‘A. Rossi Fanelli’ Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Giancarlo Ruocco
- Center for Life Nanoscience, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
48
|
Jemimah S, Sekijima M, Gromiha MM. ProAffiMuSeq: sequence-based method to predict the binding free energy change of protein-protein complexes upon mutation using functional classification. Bioinformatics 2020; 36:1725-1730. [PMID: 31713585 DOI: 10.1093/bioinformatics/btz829] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2019] [Revised: 10/23/2019] [Accepted: 11/11/2019] [Indexed: 12/20/2022] Open
Abstract
MOTIVATION Protein-protein interactions are essential for the cell and mediate various functions. However, mutations can disrupt these interactions and may cause diseases. Currently available computational methods require a complex structure as input for predicting the change in binding affinity. Further, they have not included the functional class information for the protein-protein complex. To address this, we have developed a method, ProAffiMuSeq, which predicts the change in binding free energy using sequence-based features and functional class. RESULTS Our method shows an average correlation between predicted and experimentally determined ΔΔG of 0.73 and mean absolute error (MAE) of 0.86 kcal/mol in 10-fold cross-validation and correlation of 0.75 with MAE of 0.94 kcal/mol in the test dataset. ProAffiMuSeq was also tested on an external validation set and showed results comparable to structure-based methods. Our method can be used for large-scale analysis of disease-causing mutations in protein-protein complexes without structural information. AVAILABILITY AND IMPLEMENTATION Users can access the method at https://web.iitm.ac.in/bioinfo2/proaffimuseq/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Sherlyn Jemimah
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Masakazu Sekijima
- Advanced Computational Drug Discovery Unit, Tokyo Institute of Technology, Midori-ku, Kanagawa 226-8503, Yokohama, Japan
| | - M Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.,Advanced Computational Drug Discovery Unit, Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovative Research, Tokyo Institute of Technology, Midori-ku, Kanagawa 226-8503, Yokohama, Japan
| |
Collapse
|
49
|
Insights into changes in binding affinity caused by disease mutations in protein-protein complexes. Comput Biol Med 2020; 123:103829. [DOI: 10.1016/j.compbiomed.2020.103829] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Revised: 05/20/2020] [Accepted: 05/20/2020] [Indexed: 01/11/2023]
|
50
|
Savojardo C, Martelli PL, Casadio R. Protein–Protein Interaction Methods and Protein Phase Separation. Annu Rev Biomed Data Sci 2020. [DOI: 10.1146/annurev-biodatasci-011720-104428] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
In the last decade, newly developed experimental methods have made it possible to highlight that macromolecules in the cell milieu physically interact to support physiology. This has shifted the problem of protein–protein interaction from a microscopic, electron-density scale to a mesoscopic one. Further, nowadays there is increasing evidence that proteins in the nucleus and in the cytoplasm can aggregate in membraneless organelles for different physiological reasons. In this scenario, it is urgent to face the problem of biomolecule functional annotation with efficient computational methods, suited to extract knowledge from reliable data and transfer information across different domains of investigation. Here, we revise the present state of the art of our knowledge of protein–protein interaction and the computational methods that differently implement it. Furthermore, we explore experimental and computational features of a set of proteins involved in phase separation.
Collapse
Affiliation(s)
- Castrense Savojardo
- Biocomputing Group, Department of Pharmacy and Biotechnology and Interdepartmental Center “Luigi Galvani” for Integrated Studies of Bioinformatics, Biophysics, and Biocomplexity, University of Bologna, 40126 Bologna, Italy
| | - Pier Luigi Martelli
- Biocomputing Group, Department of Pharmacy and Biotechnology and Interdepartmental Center “Luigi Galvani” for Integrated Studies of Bioinformatics, Biophysics, and Biocomplexity, University of Bologna, 40126 Bologna, Italy
| | - Rita Casadio
- Biocomputing Group, Department of Pharmacy and Biotechnology and Interdepartmental Center “Luigi Galvani” for Integrated Studies of Bioinformatics, Biophysics, and Biocomplexity, University of Bologna, 40126 Bologna, Italy
- Institute of Biomembranes, Bioenergetics, and Molecular Biotechnologies (IBIOM), Italian National Research Council (CNR), 70126 Bari, Italy
| |
Collapse
|