1
|
Banerjee T, Gosai A, Yousefi N, Garibay OO, Seal S, Balasubramanian G. Examining sialic acid derivatives as potential inhibitors of SARS-CoV-2 spike protein receptor binding domain. J Biomol Struct Dyn 2024; 42:6342-6358. [PMID: 37424217 DOI: 10.1080/07391102.2023.2234044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Accepted: 07/01/2023] [Indexed: 07/11/2023]
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has been the primary reason behind the COVID-19 global pandemic which has affected millions of lives worldwide. The fundamental cause of the infection is the molecular binding of the viral spike protein receptor binding domain (SP-RBD) with the human cell angiotensin-converting enzyme 2 (ACE2) receptor. The infection can be prevented if the binding of RBD-ACE2 is resisted by utilizing certain inhibitors or drugs that demonstrate strong binding affinity towards the SP RBD. Sialic acid based glycans found widely in human cells and tissues have notable propensity of binding to viral proteins of the coronaviridae family. Recent experimental literature have used N-acetyl neuraminic acid (Sialic acid) to create diagnostic sensors for SARS-CoV-2, but a detailed interrogation of the underlying molecular mechanisms is warranted. Here, we perform all atom molecular dynamics (MD) simulations for the complexes of certain Sialic acid-based molecules with that of SP RBD of SARS CoV-2. Our results indicate that Sialic acid not only reproduces a binding affinity comparable to the RBD-ACE2 interactions, it also assumes the longest time to dissociate completely from the protein binding pocket of SP RBD. Our predictions corroborate that a combination of electrostatic and van der Waals energies as well the polar hydrogen bond interactions between the RBD residues and the inhibitors influence free energy of binding.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Tanumoy Banerjee
- Department of Mechanical Engineering and Mechanics, Lehigh University, Bethlehem, PA, USA
| | | | - Niloofar Yousefi
- Industrial Engineering and Management Systems, University of Central Florida, Orlando, FL, USA
| | - Ozlem Ozmen Garibay
- Industrial Engineering and Management Systems, University of Central Florida, Orlando, FL, USA
| | - Sudipta Seal
- College of Medicine, Bionix Cluster, University of Central Florida, Orlando, FL, USA
- Advanced Materials Processing and Analysis Center, Dept. of Materials Science and Engineering, University of Central Florida, Orlando, FL, USA
| | - Ganesh Balasubramanian
- Department of Mechanical Engineering and Mechanics, Lehigh University, Bethlehem, PA, USA
- Institute of Functional Materials & Devices and College of Health, Lehigh University, Bethlehem, PA, USA
| |
Collapse
|
2
|
Aguayo-Acosta A, Oyervides-Muñoz MA, Rodriguez-Aguillón KO, Ovalle-Carcaño A, Romero-Castillo KD, Robles-Zamora A, Johnson M, Parra-Saldívar R, Sosa-Hernández JE. Omicron and Delta variant prevalence detection and identification during the fourth COVID-19 wave in Mexico using wastewater-based epidemiology. IJID REGIONS 2024; 10:44-51. [PMID: 38149263 PMCID: PMC10750064 DOI: 10.1016/j.ijregi.2023.11.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Revised: 11/03/2023] [Accepted: 11/06/2023] [Indexed: 12/28/2023]
Abstract
Objectives To identify the SARS-CoV-2 variants Delta and Omicron during the fourth wave of the COVID-19 pandemic in Mexico using samples taken from 19 locations in 18 out of the 32 states. Methods The genetic material concentration was done with PEG/NaCl precipitation, SARS-CoV-2 presence was confirmed by reverse transcriptase-quantitative polymerase chain reaction assay, the variant detection was carried out using a commercial mutation detection panel kit, and variant/mutation confirmation was done by amplicon sequencing of receptor-binding domain target region. The study used 41 samples. Results The Delta variant was confirmed in two samples during August 2021 (Querétaro and CDMX) and in three samples during November 2021 (Aguascalientes, Ciudad Juárez campuses, and Nuevo Leon). In December 2021, another sample with the Delta variant was confirmed in Nuevo Leon. Between January to March 2022 only the presence of Omicron was confirmed, (variant BA.1). Additionally, in this period six samples were identified with the status "Variant Not Determined". Conclusion To our knowledge, this study is one of the first to identify Omicron and Delta variants with polymerase chain reaction in Mexico and Latin America and its distribution across the country with 56% Mexican states making it a viable alternative for variant detection without conducting a large quantity of sequencing of clinical tests.
Collapse
Affiliation(s)
- Alberto Aguayo-Acosta
- Tecnologico de Monterrey, Institute of Advanced Materials for Sustainable Manufacturing, Monterrey, Mexico
- Tecnologico de Monterrey, School of Engineering and Sciences, Monterrey, Mexico
| | - Mariel Araceli Oyervides-Muñoz
- Tecnologico de Monterrey, Institute of Advanced Materials for Sustainable Manufacturing, Monterrey, Mexico
- Tecnologico de Monterrey, School of Engineering and Sciences, Monterrey, Mexico
| | - Kassandra O. Rodriguez-Aguillón
- Tecnologico de Monterrey, Institute of Advanced Materials for Sustainable Manufacturing, Monterrey, Mexico
- Tecnologico de Monterrey, School of Engineering and Sciences, Monterrey, Mexico
| | - Antonio Ovalle-Carcaño
- Tecnologico de Monterrey, Institute of Advanced Materials for Sustainable Manufacturing, Monterrey, Mexico
- Tecnologico de Monterrey, School of Engineering and Sciences, Monterrey, Mexico
| | | | | | - Marc Johnson
- Department of Molecular Microbiology and Immunology, University of Missouri School of Medicine, Columbia, USA
| | - Roberto Parra-Saldívar
- Tecnologico de Monterrey, Institute of Advanced Materials for Sustainable Manufacturing, Monterrey, Mexico
- Tecnologico de Monterrey, School of Engineering and Sciences, Monterrey, Mexico
| | - Juan Eduardo Sosa-Hernández
- Tecnologico de Monterrey, Institute of Advanced Materials for Sustainable Manufacturing, Monterrey, Mexico
- Tecnologico de Monterrey, School of Engineering and Sciences, Monterrey, Mexico
| |
Collapse
|
3
|
Bhadane R, Salo-Ahen OMH. High-Throughput Molecular Dynamics-Based Alchemical Free Energy Calculations for Predicting the Binding Free Energy Change Associated with the Selected Omicron Mutations in the Spike Receptor-Binding Domain of SARS-CoV-2. Biomedicines 2022; 10:2779. [PMID: 36359299 PMCID: PMC9687918 DOI: 10.3390/biomedicines10112779] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2022] [Revised: 10/17/2022] [Accepted: 10/25/2022] [Indexed: 11/10/2023] Open
Abstract
The ongoing pandemic caused by SARS-CoV-2 has gone through various phases. Since the initial outbreak, the virus has mutated several times, with some lineages showing even stronger infectivity and faster spread than the original virus. Among all the variants, omicron is currently classified as a variant of concern (VOC) by the World Health Organization, as the previously circulating variants have been replaced by it. In this work, we have focused on the mutations observed in omicron sub lineages BA.1, BA.2, BA.4 and BA.5, particularly at the receptor-binding domain (RBD) of the spike protein that is responsible for the interactions with the host ACE2 receptor and binding of antibodies. Studying such mutations is particularly important for understanding the viral infectivity, spread of the disease and for tracking the escape routes of this virus from antibodies. Molecular dynamics (MD) based alchemical free energy calculations have been shown to be very accurate in predicting the free energy change, due to a mutation that could have a deleterious or a stabilizing effect on either the protein itself or its binding affinity to another protein. Here, we investigated the significance of five spike RBD mutations on the stability of the spike protein binding to ACE2 by free energy calculations using high throughput MD simulations. For comparison, we also used conventional MD simulations combined with a Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) based approach, and compared our results with the available experimental data. Overall, the alchemical free energy calculations performed far better than the MM-GBSA approach in predicting the individual impact of the mutations. When considering the experimental variation, the alchemical free energy method was able to produce a relatively accurate prediction for N501Y, the mutant that has previously been reported to increase the binding affinity to hACE2. On the other hand, the other individual mutations seem not to have a significant effect on the spike RBD binding affinity towards hACE2.
Collapse
Affiliation(s)
- Rajendra Bhadane
- Structural Bioinformatics Laboratory, Faculty of Science and Engineering, Biochemistry, Åbo Akademi University, FI-20520 Turku, Finland
- Pharmaceutical Sciences Laboratory, Faculty of Science and Engineering, Pharmacy, Åbo Akademi University, FI-20520 Turku, Finland
| | - Outi M. H. Salo-Ahen
- Structural Bioinformatics Laboratory, Faculty of Science and Engineering, Biochemistry, Åbo Akademi University, FI-20520 Turku, Finland
- Pharmaceutical Sciences Laboratory, Faculty of Science and Engineering, Pharmacy, Åbo Akademi University, FI-20520 Turku, Finland
| |
Collapse
|
4
|
Coderc de Lacam EG, Blazhynska M, Chen H, Gumbart JC, Chipot C. When the Dust Has Settled: Calculation of Binding Affinities from First Principles for SARS-CoV-2 Variants with Quantitative Accuracy. J Chem Theory Comput 2022; 18:5890-5900. [PMID: 36108303 PMCID: PMC9518821 DOI: 10.1021/acs.jctc.2c00604] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Indexed: 11/30/2022]
Abstract
Accurate determination of binding free energy is pivotal for the study of many biological processes and has been applied in a number of theoretical investigations to compare the affinity of severe acute respiratory syndrome coronavirus 2 variants toward the host cell. Diversity of these variants challenges the development of effective general therapies, their transmissibility relying either on an increased affinity toward their dedicated human receptor, the angiotensin-converting enzyme 2 (ACE2), or on escaping the immune response. Now that robust structural data are available, we have determined with utmost accuracy the standard binding free energy of the receptor-binding domain to the most widespread variants, namely, Alpha, Beta, Delta, and Omicron BA.2, as well as the wild type (WT) in complex either with ACE2 or with antibodies, namely, S2E12 and H11-D4, using a rigorous theoretical framework that combines molecular dynamics and potential-of-mean-force calculations. Our results show that an appropriate starting structure is crucial to ensure appropriate reproduction of the binding affinity, allowing the variants to be compared. They also emphasize the necessity to apply the relevant methodology, bereft of any shortcut, to account for all the contributions to the standard binding free energy. Our estimates of the binding affinities support the view that while the Alpha and Beta variants lean on an increased affinity toward the host cell, the Delta and Omicron BA.2 variants choose immune escape. Moreover, the S2E12 antibody, already known to be active against the WT (Starr et al., 2021; Mlcochova et al., 2021), proved to be equally effective against the Delta variant. In stark contrast, H11-D4 retains a low affinity toward the WT compared to that of ACE2 for the latter. Assuming robust structural information, the methodology employed herein successfully addresses the challenging protein-protein binding problem in the context of coronavirus disease 2019 while offering promising perspectives for predictive studies of ever-emerging variants.
Collapse
Affiliation(s)
- Emma Goulard Coderc de Lacam
- Laboratoire International Associé Centre
National de la Recherche Scientifique et University of Illinois at Urbana-Champaign,
Unité Mixte de Recherche No 7019, Université de
Lorraine, B.P. 70239, Vandœuvre-lès-Nancy Cedex54506,
France
| | - Marharyta Blazhynska
- Laboratoire International Associé Centre
National de la Recherche Scientifique et University of Illinois at Urbana-Champaign,
Unité Mixte de Recherche No 7019, Université de
Lorraine, B.P. 70239, Vandœuvre-lès-Nancy Cedex54506,
France
| | - Haochuan Chen
- Laboratoire International Associé Centre
National de la Recherche Scientifique et University of Illinois at Urbana-Champaign,
Unité Mixte de Recherche No 7019, Université de
Lorraine, B.P. 70239, Vandœuvre-lès-Nancy Cedex54506,
France
| | - James C. Gumbart
- School of Physics, Georgia Institute of
Technology, Atlanta, Georgia30332, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre
National de la Recherche Scientifique et University of Illinois at Urbana-Champaign,
Unité Mixte de Recherche No 7019, Université de
Lorraine, B.P. 70239, Vandœuvre-lès-Nancy Cedex54506,
France
- Theoretical and Computational Biophysics Group, Beckman
Institute, and Department of Physics, University of Illinois at
Urbana-Champaign, UrbanaIllinois61802, United
States
- Department of Biochemistry and Molecular Biology,
The University of Chicago, 929 E. 57th Street W225, Chicago,
Illinois60637, United States
| |
Collapse
|
5
|
Rodríguez Y, Cardoze SM, Obineche OW, Melo C, Persaud A, Fernández Romero JA. Small Molecules Targeting SARS-CoV-2 Spike Glycoprotein Receptor-Binding Domain. ACS OMEGA 2022; 7:28779-28789. [PMID: 35991504 PMCID: PMC9380817 DOI: 10.1021/acsomega.2c00844] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/10/2022] [Accepted: 07/06/2022] [Indexed: 06/15/2023]
Abstract
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the coronavirus disease 2019 (COVID-19) pandemic. Several variants of SARS-CoV-2 have emerged worldwide. These variants show different transmissibility infectivity due to mutations in the viral spike (S) glycoprotein that interacts with the human angiotensin-converting enzyme 2 (hACE2) receptor and facilitates viral entry into target cells. Despite the effective SARS-CoV-2 vaccines, we still need to identify selective antivirals, and the S glycoprotein is a key target to neutralize the virus. We hypothesize that small molecules could disrupt the interaction of S glycoprotein with hACE2 and inhibit viral entry. We analyzed the S glycoprotein-hACE2 complex structure (PDB: 7DF4) and created models for different viral variants using visual molecular dynamics (VMD) and molecular operating environment (MOE) programs. Moreover, we started the hits search by performing structure-based molecular docking virtual screening of commercially available small molecules against S glycoprotein models using OEDocking FRED-4.0.0.0 software. The FRED-4.0.0.0 Chemguass4 scoring function was used to rank the small molecules based on their affinities. The best candidate compounds were purchased and tested using a standard SARS-CoV-2 pseudotyped cell-based bioassay to investigate their antiviral activity. Three of these compounds, alone or in combination, showed antiviral selectivity. These small molecules may lead to an effective antiviral treatment or serve as probes to better understand the biology of SARS-CoV-2.
Collapse
Affiliation(s)
- Yoel Rodríguez
- Department
of Natural Sciences, Hostos Community College
of The City University of New York, 500 Grand Concourse, Bronx, New York, New York 10451, United States
- Department
of Pharmacological Sciences, Icahn School
of Medicine at Mount Sinai, 1425 Madison Avenue, New York, New York 10029, United
States
| | - Scarlet Martínez Cardoze
- Department
of Natural Sciences, Hostos Community College
of The City University of New York, 500 Grand Concourse, Bronx, New York, New York 10451, United States
| | - Onyinyechi W. Obineche
- Department
of Natural Sciences, Hostos Community College
of The City University of New York, 500 Grand Concourse, Bronx, New York, New York 10451, United States
| | - Claudia Melo
- Department
of Science, Borough of Manhattan Community
College of The City University of New York, 199 Chambers St, New York, New York 10007, United States
- Brooklyn
College of The City University of New York, 2900 Bedford Avenue, New
York, New York 11210, United States
| | - Ashanna Persaud
- Department
of Science, Borough of Manhattan Community
College of The City University of New York, 199 Chambers St, New York, New York 10007, United States
| | - José A. Fernández Romero
- Department
of Science, Borough of Manhattan Community
College of The City University of New York, 199 Chambers St, New York, New York 10007, United States
- Center
for Biomedical Research, The Population
Council, 1230 York Avenue, New York, New York 10065, United States
| |
Collapse
|
6
|
Ghasemitarei M, Privat-Maldonado A, Yusupov M, Rahnama S, Bogaerts A, Ejtehadi MR. Effect of Cysteine Oxidation in SARS-CoV-2 Receptor-Binding Domain on Its Interaction with Two Cell Receptors: Insights from Atomistic Simulations. J Chem Inf Model 2022; 62:129-141. [PMID: 34965734 PMCID: PMC8751020 DOI: 10.1021/acs.jcim.1c00853] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Indexed: 12/15/2022]
Abstract
Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.
Collapse
Affiliation(s)
- Maryam Ghasemitarei
- Department
of Physics, Sharif University of Technology, Tehran 14588-89694, Iran
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium
| | - Angela Privat-Maldonado
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium
| | - Maksudbek Yusupov
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium
- Laboratory
of Thermal Physics of Multiphase Systems, Arifov Institute of Ion-Plasma
and Laser Technologies, Academy of Sciences
of Uzbekistan, Durmon
yuli str. 33, 100125 Tashkent, Uzbekistan
| | - Shadi Rahnama
- Institute
for Nanoscience & Nanotechnology (INST), Sharif University of Technology, Azadi Avenue, Tehran 14588-89694, Iran
| | - Annemie Bogaerts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium
| | | |
Collapse
|
7
|
Ngo VA, Jha RK. Correction: Identifying key determinants and dynamics of SARS-CoV-2/ACE2 tight interaction. PLoS One 2021; 16:e0259705. [PMID: 34731214 PMCID: PMC8565735 DOI: 10.1371/journal.pone.0259705] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|