1
|
Gupta G, Verkhivker G. Exploring Binding Pockets in the Conformational States of the SARS-CoV-2 Spike Trimers for the Screening of Allosteric Inhibitors Using Molecular Simulations and Ensemble-Based Ligand Docking. Int J Mol Sci 2024; 25:4955. [PMID: 38732174 PMCID: PMC11084335 DOI: 10.3390/ijms25094955] [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: 03/19/2024] [Revised: 04/24/2024] [Accepted: 04/30/2024] [Indexed: 05/13/2024] Open
Abstract
Understanding mechanisms of allosteric regulation remains elusive for the SARS-CoV-2 spike protein, despite the increasing interest and effort in discovering allosteric inhibitors of the viral activity and interactions with the host receptor ACE2. The challenges of discovering allosteric modulators of the SARS-CoV-2 spike proteins are associated with the diversity of cryptic allosteric sites and complex molecular mechanisms that can be employed by allosteric ligands, including the alteration of the conformational equilibrium of spike protein and preferential stabilization of specific functional states. In the current study, we combine conformational dynamics analysis of distinct forms of the full-length spike protein trimers and machine-learning-based binding pocket detection with the ensemble-based ligand docking and binding free energy analysis to characterize the potential allosteric binding sites and determine structural and energetic determinants of allosteric inhibition for a series of experimentally validated allosteric molecules. The results demonstrate a good agreement between computational and experimental binding affinities, providing support to the predicted binding modes and suggesting key interactions formed by the allosteric ligands to elicit the experimentally observed inhibition. We establish structural and energetic determinants of allosteric binding for the experimentally known allosteric molecules, indicating a potential mechanism of allosteric modulation by targeting the hinges of the inter-protomer movements and blocking conformational changes between the closed and open spike trimer forms. The results of this study demonstrate that combining ensemble-based ligand docking with conformational states of spike protein and rigorous binding energy analysis enables robust characterization of the ligand binding modes, the identification of allosteric binding hotspots, and the prediction of binding affinities for validated allosteric modulators, which is consistent with the experimental data. This study suggested that the conformational adaptability of the protein allosteric sites and the diversity of ligand bound conformations are both in play to enable efficient targeting of allosteric binding sites and interfere with the conformational changes.
Collapse
Affiliation(s)
- Grace Gupta
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, USA;
| | - Gennady Verkhivker
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, USA;
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| |
Collapse
|
2
|
Zeng L, Lu Y, Yan W, Yang Y. A Protein Co-Conservation Network Model Characterizes Mutation Effects on SARS-CoV-2 Spike Protein. Int J Mol Sci 2023; 24:ijms24043255. [PMID: 36834664 PMCID: PMC9960056 DOI: 10.3390/ijms24043255] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Revised: 02/03/2023] [Accepted: 02/04/2023] [Indexed: 02/10/2023] Open
Abstract
The emergence of numerous variants of SARS-CoV-2 has presented challenges to the global efforts to control the COVID-19 pandemic. The major mutation is in the SARS-CoV-2 viral envelope spike protein that is responsible for virus attachment to the host, and is the main target for host antibodies. It is critically important to study the biological effects of the mutations to understand the mechanisms of how mutations alter viral functions. Here, we propose a protein co-conservation weighted network (PCCN) model only based on the protein sequence to characterize the mutation sites by topological features and to investigate the mutation effects on the spike protein from a network view. Frist, we found that the mutation sites on the spike protein had significantly larger centrality than the non-mutation sites. Second, the stability changes and binding free energy changes in the mutation sites were positively significantly correlated with their neighbors' degree and the shortest path length separately. The results indicate that our PCCN model provides new insights into mutations on spike proteins and reflects the mutation effects on protein function alternations.
Collapse
Affiliation(s)
- Lianjie Zeng
- School of Computer Science & Technology, Soochow University, Suzhou 215000, China
- Collaborative Innovation Center of Novel Software Technology and Industrialization, Nanjing 210000, China
| | - Yitan Lu
- Department of Bioinformatics, School of Biology and Basic Medical Sciences, Medical College of Soochow University, Suzhou 215123, China
| | - Wenying Yan
- Department of Bioinformatics, School of Biology and Basic Medical Sciences, Medical College of Soochow University, Suzhou 215123, China
- Jiangsu Province Engineering Research Center of Precision Diagnostics and Therapeutics Development, Suzhou 215123, China
- Correspondence: (W.Y.); (Y.Y.)
| | - Yang Yang
- School of Computer Science & Technology, Soochow University, Suzhou 215000, China
- Collaborative Innovation Center of Novel Software Technology and Industrialization, Nanjing 210000, China
- Correspondence: (W.Y.); (Y.Y.)
| |
Collapse
|
3
|
Erba F, Di Paola L, Di Venere A, Mastrangelo E, Cossu F, Mei G, Minicozzi V. Head or tail? A molecular dynamics approach to the complex structure of TNF-associated factor TRAF2. Biomol Concepts 2023; 14:bmc-2022-0031. [PMID: 37377424 DOI: 10.1515/bmc-2022-0031] [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: 03/31/2023] [Accepted: 06/09/2023] [Indexed: 06/29/2023] Open
Abstract
Tumor necrosis factor receptor-associated factor proteins (TRAFs) are trimeric proteins that play a fundamental role in signaling, acting as intermediaries between the tumor necrosis factor (TNF) receptors and the proteins that transmit the downstream signal. The monomeric subunits of all the TRAF family members share a common tridimensional structure: a C-terminal globular domain and a long coiled-coil tail characterizing the N-terminal section. In this study, the dependence of the TRAF2 dynamics on the length of its tail was analyzed in silico. In particular, we used the available crystallographic structure of a C-terminal fragment of TRAF2 (168 out of 501 a.a.), TRAF2-C, and that of a longer construct, addressed as TRAF2-plus, that we have re-constructed using the AlphaFold2 code. The results indicate that the longer N-terminal tail of TRAF2-plus has a strong influence on the dynamics of the globular regions in the protein C-terminal head. In fact, the quaternary interactions among the TRAF2-C subunits change asymmetrically in time, while the movements of TRAF2-plus monomers are rather limited and more ordered than those of the shorter construct. Such findings shed a new light on the dynamics of TRAF subunits and on the protein mechanism in vivo, since TRAF monomer-trimer equilibrium is crucial for several reasons (receptor recognition, membrane binding, hetero-oligomerization).
Collapse
Affiliation(s)
- Fulvio Erba
- Department of Clinical Science and Translational Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00133 Rome, Italy
| | - Luisa Di Paola
- Unit of Chemical-Physics Fundamentals in Chemical Engineering, Department of Engineering, University Campus Bio-Medico of Rome, Via Álvaro del Portillo 21, 00128 Rome, Italy
| | - Almerinda Di Venere
- Department of Experimental Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00133 Rome, Italy
| | - Eloise Mastrangelo
- National Research Council (IBF-CNR) Milan Unit, Institute of Biophysics, Via Celoria 26, 20133 Milan, Italy
| | - Federica Cossu
- National Research Council (IBF-CNR) Milan Unit, Institute of Biophysics, Via Celoria 26, 20133 Milan, Italy
| | - Giampiero Mei
- Department of Experimental Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00133 Rome, Italy
| | - Velia Minicozzi
- Department of Physics and INFN, Tor Vergata University of Rome, Via Della Ricerca Scientifica 1, 00133 Rome, Italy
| |
Collapse
|
4
|
Barroso da Silva FL, Giron CC, Laaksonen A. Electrostatic Features for the Receptor Binding Domain of SARS-COV-2 Wildtype and Its Variants. Compass to the Severity of the Future Variants with the Charge-Rule. J Phys Chem B 2022; 126:6835-6852. [PMID: 36066414 DOI: 10.1021/acs.jpcb.2c04225] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Electrostatic intermolecular interactions are important in many aspects of biology. We have studied the main electrostatic features involved in the interaction of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with the human receptor Angiotensin-converting enzyme 2 (ACE2). As the principal computational tool, we have used the FORTE approach, capable to model proton fluctuations and computing free energies for a very large number of protein-protein systems under different physical-chemical conditions, here focusing on the RBD-ACE2 interactions. Both the wild-type and all critical variants are included in this study. From our large ensemble of extensive simulations, we obtain, as a function of pH, the binding affinities, charges of the proteins, their charge regulation capacities, and their dipole moments. In addition, we have calculated the pKas for all ionizable residues and mapped the electrostatic coupling between them. We are able to present a simple predictor for the RBD-ACE2 binding based on the data obtained for Alpha, Beta, Gamma, Delta, and Omicron variants, as a linear correlation between the total charge of the RBD and the corresponding binding affinity. This "RBD charge rule" should work as a quick test of the degree of severity of the coming SARS-CoV-2 variants in the future.
Collapse
Affiliation(s)
- Fernando L Barroso da Silva
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. café, s/no-campus da USP, BR-14040-903 Ribeirão Preto, SP, Brazil.,Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Carolina Corrêa Giron
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. café, s/no-campus da USP, BR-14040-903 Ribeirão Preto, SP, Brazil.,Hospital de Clínicas, Universidade Federal do Triângulo Mineiro, Av. Getúlio Guaritá, 38025-440 Uberaba, MG, Brazil
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden.,State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing, 210009, P. R. China.,Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Aleea Grigore Ghica-Voda, 41A, 700487 Iasi, Romania.,Department of Engineering Sciences and Mathematics, Division of Energy Science, Luleå University of Technology, SE-97187 Luleå, Sweden.,Department of Chemical and Geological Sciences, Campus Monserrato, University of Cagliari, SS 554 bivio per Sestu, 09042 Monserrato, Italy
| |
Collapse
|
5
|
In-Silico Characterization of von Willebrand Factor Bound to FVIII. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12157855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Factor VIII belongs to the coagulation cascade and is expressed as a long pre-protein (mature form, 2351 amino acids long). FVIII is deficient or defective in hemophilic A patients, who need to be treated with hemoderivatives or recombinant FVIII substitutes, i.e., biologic drugs. The interaction between FVIII and von Willebrand factor (VWF) influences the pharmacokinetics of FVIII medications. In vivo, full-length FVIII (FL-FVIII) is secreted in a plasma-inactive form, which includes the B domain, which is then proteolyzed by thrombin protease activity, leading to an inactive plasma intermediate. In this work, we analyzed through a computational approach the binding of VWF with two structure models of FVIII (secreted full-length with B domain, and B domain-deleted FVIII). We included in our analysis the atomic model of efanesoctocog alfa, a novel and investigational recombinant FVIII medication, in which the VWF is covalently linked to FVIII. We carried out a structural analysis of VWF/FVIII interfaces by means of protein–protein docking, PISA (Proteins, Interfaces, Structures and Assemblies), and protein contact networks (PCN) analyses. Accordingly, our computational approaches to previously published experimental data demonstrated that the domains A3-C1 of B domain-deleted FVIII (BDD-FVIII) is the preferential binding site for VWF. Overall, our computational approach applied to topological analysis of protein–protein interface can be aimed at the rational design of biologic drugs other than FVIII medications.
Collapse
|
6
|
A Statistical Journey through the Topological Determinants of the β2 Adrenergic Receptor Dynamics. ENTROPY 2022; 24:e24070998. [PMID: 35885221 PMCID: PMC9318934 DOI: 10.3390/e24070998] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/28/2022] [Revised: 07/12/2022] [Accepted: 07/15/2022] [Indexed: 02/06/2023]
Abstract
Activation of G-protein-coupled receptors (GPCRs) is mediated by molecular switches throughout the transmembrane region of the receptor. In this work, we continued along the path of a previous computational study wherein energy transport in the β2 Adrenergic Receptor (β2-AR) was examined and allosteric switches were identified in the molecular structure through the reorganization of energy transport networks during activation. In this work, we further investigated the allosteric properties of β2-AR, using Protein Contact Networks (PCNs). In this paper, we report an extensive statistical analysis of the topological and structural properties of β2-AR along its molecular dynamics trajectory to identify the activation pattern of this molecular system. The results show a distinct character to the activation that both helps to understand the allosteric switching previously identified and confirms the relevance of the network formalism to uncover relevant functional features of protein molecules.
Collapse
|