1
|
Mahapatra S, Kar P. Computational biophysical characterization of the effect of gatekeeper mutations on the binding of ponatinib to the FGFR kinase. Arch Biochem Biophys 2024; 758:110070. [PMID: 38909834 DOI: 10.1016/j.abb.2024.110070] [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/29/2024] [Revised: 05/15/2024] [Accepted: 06/20/2024] [Indexed: 06/25/2024]
Abstract
Fibroblast Growth Factor Receptor (FGFR) is connected to numerous downstream signalling cascades regulating cellular behavior. Any dysregulation leads to a plethora of illnesses, including cancer. Therapeutics are available, but drug resistance driven by gatekeeper mutation impedes the treatment. Ponatinib is an FDA-approved drug against BCR-ABL kinase and has shown effective results against FGFR-mediated carcinogenesis. Herein, we undertake molecular dynamics simulation-based analysis on ponatinib against all the FGFR isoforms having Val to Met gatekeeper mutations. The results suggest that ponatinib is a potent and selective inhibitor for FGFR1, FGFR2, and FGFR4 gatekeeper mutations. The extensive electrostatic and van der Waals interaction network accounts for its high potency. The FGFR3_VM mutation has shown resistance towards ponatinib, which is supported by their lesser binding affinity than wild-type complexes. The disengaged molecular brake and engaged hydrophobic spine were believed to be the driving factors for weak protein-ligand interaction. Taken together, the inhibitory and structural characteristics exhibited by ponatinib may aid in thwarting resistance based on Val-to-Met gatekeeper mutations at an earlier stage of treatment and advance the design and development of other inhibitors targeted at FGFRs harboring gatekeeper mutations.
Collapse
Affiliation(s)
- Subhasmita Mahapatra
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore, 453552, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore, 453552, Madhya Pradesh, India.
| |
Collapse
|
2
|
Franke L, Peter C. Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding. J Chem Theory Comput 2023; 19:2985-2995. [PMID: 37122117 DOI: 10.1021/acs.jctc.2c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.
Collapse
Affiliation(s)
- Leon Franke
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
- Konstanz Research School Chemical Biology, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| | - Christine Peter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| |
Collapse
|
3
|
Manrique PD, Chakraborty S, Henderson R, Edwards RJ, Mansbach R, Nguyen K, Stalls V, Saunders C, Mansouri K, Acharya P, Korber B, Gnanakaran S. Network analysis uncovers the communication structure of SARS-CoV-2 spike protein identifying sites for immunogen design. iScience 2023; 26:105855. [PMID: 36590900 PMCID: PMC9791713 DOI: 10.1016/j.isci.2022.105855] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Revised: 11/17/2022] [Accepted: 12/19/2022] [Indexed: 12/27/2022] Open
Abstract
The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike's stability. Taken together, the network approach serves as a valuable tool to evaluate communication "hot-spots" in proteins to guide design of stable immunogens.
Collapse
Affiliation(s)
- Pedro D. Manrique
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| | - Srirupa Chakraborty
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| | - Rory Henderson
- Duke Human Vaccine Institute, Duke University, Durham, NC 27710, USA
- Department of Medicine, Duke University, Durham, NC 27710, USA
| | - Robert J. Edwards
- Duke Human Vaccine Institute, Duke University, Durham, NC 27710, USA
- Department of Medicine, Duke University, Durham, NC 27710, USA
| | - Rachael Mansbach
- Physics Department, Concordia University, Montreal, QC H4B IR6, Canada
| | - Kien Nguyen
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| | - Victoria Stalls
- Duke Human Vaccine Institute, Duke University, Durham, NC 27710, USA
| | - Carrie Saunders
- Duke Human Vaccine Institute, Duke University, Durham, NC 27710, USA
| | - Katayoun Mansouri
- Duke Human Vaccine Institute, Duke University, Durham, NC 27710, USA
| | - Priyamvada Acharya
- Duke Human Vaccine Institute, Duke University, Durham, NC 27710, USA
- Department of Surgery, Duke University, Durham, NC 27710, USA
| | - Bette Korber
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| | - S. Gnanakaran
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| |
Collapse
|
4
|
Verkhivker GM, Agajanian S, Oztas D, Gupta G. Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation. J Biomol Struct Dyn 2022; 40:9724-9741. [PMID: 34060425 DOI: 10.1080/07391102.2021.1933594] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA.,Depatment of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, USA
| | - Steve Agajanian
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
| | - Deniz Oztas
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
| | - Grace Gupta
- Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
| |
Collapse
|
5
|
Lata S, Akif M. Probing structural basis for enhanced binding of SARS-CoV-2 P.1 variant spike protein with the human ACE2 receptor. J Cell Biochem 2022; 123:1207-1221. [PMID: 35620980 PMCID: PMC9347910 DOI: 10.1002/jcb.30276] [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: 01/21/2022] [Revised: 04/22/2022] [Accepted: 05/10/2022] [Indexed: 11/25/2022]
Abstract
The initial step of infection by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) involves the binding of receptor binding domain (RBD) of the spike protein to the angiotensin converting enzyme 2 (ACE2) receptor. Each successive wave of SARS‐CoV‐2 reports emergence of many new variants, which is associated with mutations in the RBD as well as other parts of the spike protein. These mutations are reported to have enhanced affinity towards the ACE2 receptor as well as are also crucial for the virus transmission. Many computational and experimental studies have demonstrated the effect of individual mutation on the RBD‐ACE2 binding. However, the cumulative effect of mutations on the RBD and away from the RBD was not investigated in detail. We report here a comparative analysis on the structural communication and dynamics of the RBD and truncated S1 domain of spike protein in complex with the ACE2 receptor from SARS‐CoV‐2 wild type and its P.1 variant. Our integrative network and dynamics approaches highlighted a subtle conformational changes in the RBD as well as truncated S1 domain of spike protein at the protein contact level, responsible for the increased affinity with the ACE2 receptor. Moreover, our study also identified the commonalities and differences in the dynamics of the interactions between spike protein of SARS‐CoV‐2 wild type and its P.1 variant with the ACE2 receptor. Further, our investigation yielded an understanding towards identification of the unique RBD residues crucial for the interaction with the ACE2 host receptor. Overall, the study provides an insight for designing better therapeutics against the circulating P.1 variants as well as other future variants.
Collapse
Affiliation(s)
- Surabhi Lata
- Laboratory of Structural Biology, Department of Biochemistry, School of Life Sciences, University of Hyderabad, Hyderabad, Telangana, India
| | - Mohd Akif
- Laboratory of Structural Biology, Department of Biochemistry, School of Life Sciences, University of Hyderabad, Hyderabad, Telangana, India
| |
Collapse
|
6
|
Behbahani M, Mohabatkar H, Hosseini B. In silico design of quadruplex aptamers against the spike protein of SARS-CoV-2. INFORMATICS IN MEDICINE UNLOCKED 2021; 26:100757. [PMID: 34664030 PMCID: PMC8514331 DOI: 10.1016/j.imu.2021.100757] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Revised: 10/06/2021] [Accepted: 10/11/2021] [Indexed: 12/24/2022] Open
Abstract
Nucleic acid aptamers are short sequences of nucleic acid ligands that bind to a specific target molecule. Aptamers are experimentally nominated using the well-designed SELEX (systematic evolution of ligands by exponential enrichment) method. Here, we designed a new method for diagnosis and blocking SARS-CoV-2 based on G-quadruplex aptamer. This aptamer was developed against the receptor-binding domain (RBD) region of the spike protein. In the current study, ten quadruplex DNA aptamers entitled AP1, AP2, AP3, AP4, AP5, AP6, AP7, AP8, AP9, and AP10 were designed in silico and had high HADDOCK scores. One quadruplex aptamer sequence (AP1) was selected based on the interaction with RBD of SARS-CoV-2. Results showed that AP1 aptamer could be used as an agent in the diagnosis and therapy of SARS-CoV-2, although more works are still needed.
Collapse
Affiliation(s)
- Mandana Behbahani
- Department of Biotechnology, Faculty of Biological Science and Technology, University of Isfahan, Hezar Jareeb St., Isfahan, 81746-73441, Iran
| | - Hassan Mohabatkar
- Department of Biotechnology, Faculty of Biological Science and Technology, University of Isfahan, Hezar Jareeb St., Isfahan, 81746-73441, Iran
| | - Barumand Hosseini
- Department of Biotechnology, Faculty of Biological Science and Technology, University of Isfahan, Hezar Jareeb St., Isfahan, 81746-73441, Iran
| |
Collapse
|
7
|
Verkhivker GM, Agajanian S, Oztas DY, Gupta G. Landscape-Based Mutational Sensitivity Cartography and Network Community Analysis of the SARS-CoV-2 Spike Protein Structures: Quantifying Functional Effects of the Circulating D614G Variant. ACS OMEGA 2021; 6:16216-16233. [PMID: 34179666 PMCID: PMC8223427 DOI: 10.1021/acsomega.1c02336] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Accepted: 06/01/2021] [Indexed: 05/05/2023]
Abstract
We developed and applied a computational approach to simulate functional effects of the global circulating mutation D614G of the SARS-CoV-2 spike protein. All-atom molecular dynamics simulations are combined with deep mutational scanning and analysis of the residue interaction networks to investigate conformational landscapes and energetics of the SARS-CoV-2 spike proteins in different functional states of the D614G mutant. The results of conformational dynamics and analysis of collective motions demonstrated that the D614 site plays a key regulatory role in governing functional transitions between open and closed states. Using mutational scanning and sensitivity analysis of protein residues, we identified the stability hotspots in the SARS-CoV-2 spike structures of the mutant trimers. The results suggest that the D614G mutation can induce the increased stability of the open form acting as a driver of conformational changes, which may result in the increased exposure to the host receptor and promote infectivity of the virus. The network community analysis of the SARS-CoV-2 spike proteins showed that the D614G mutation can enhance long-range couplings between domains and strengthen the interdomain interactions in the open form, supporting the reduced shedding mechanism. This study provides the landscape-based perspective and atomistic view of the allosteric interactions and stability hotspots in the SARS-CoV-2 spike proteins, offering a useful insight into the molecular mechanisms underpinning functional effects of the global circulating mutations.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Steve Agajanian
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange, California 92866, United States
| | - Deniz Yasar Oztas
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange, California 92866, United States
| | - Grace Gupta
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange, California 92866, United States
| |
Collapse
|
8
|
Kumari A, Mittal L, Srivastava M, Asthana S. Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations. J Biomol Struct Dyn 2021; 40:9287-9305. [PMID: 34029506 DOI: 10.1080/07391102.2021.1927844] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The main protease, Mpro/3CLpro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the Mpro mechanism of action. A comparative analysis of both forms of Mpro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of Mpro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of Mpro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent Mpro inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Anita Kumari
- Translational Health Science and Technology Institute (THSTI), Faridabad, Haryana, India
| | - Lovika Mittal
- Translational Health Science and Technology Institute (THSTI), Faridabad, Haryana, India
| | - Mitul Srivastava
- Translational Health Science and Technology Institute (THSTI), Faridabad, Haryana, India
| | - Shailendra Asthana
- Translational Health Science and Technology Institute (THSTI), Faridabad, Haryana, India
| |
Collapse
|