skip to main content


Title: A multiple-step in silico screening protocol to identify allosteric inhibitors of Spike–hACE2 binding
While the COVID-19 pandemic continues to worsen, effective medicines that target the life cycle of SARS-CoV-2 are still under development. As more highly infective and dangerous variants of the coronavirus emerge, the protective power of vaccines will decrease or vanish. Thus, the development of drugs, which are free of drug resistance is direly needed. The aim of this study is to identify allosteric binding modulators from a large compound library to inhibit the binding between the Spike protein of the SARS-CoV-2 virus and human angiotensin-converting enzyme 2 (hACE2). The binding of the Spike protein to hACE2 is the first step of the infection of host cells by the coronavirus. We first built a compound library containing 77 448 antiviral compounds. Molecular docking was then conducted to preliminarily screen compounds which can potently bind to the Spike protein at two allosteric binding sites. Next, molecular dynamics simulations were performed to accurately calculate the binding affinity between the spike protein and an identified compound from docking screening and to investigate whether the compound can interfere with the binding between the Spike protein and hACE2. We successfully identified two possible drug binding sites on the Spike protein and discovered a series of antiviral compounds which can weaken the interaction between the Spike protein and hACE2 receptor through conformational changes of the key Spike residues at the Spike–hACE2 binding interface induced by the binding of the ligand at the allosteric binding site. We also applied our screening protocol to another compound library which consists of 3407 compounds for which the inhibitory activities of Spike/hACE2 binding were measured. Encouragingly, in vitro data supports that the identified compounds can inhibit the Spike–ACE2 binding. Thus, we developed a promising computational protocol to discover allosteric inhibitors of the binding of the Spike protein of SARS-CoV-2 to the hACE2 receptor, and several promising allosteric modulators were discovered.  more » « less
Award ID(s):
1955260
PAR ID:
10325217
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
24
Issue:
7
ISSN:
1463-9076
Page Range / eLocation ID:
4305 to 4316
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The COVID-19 pandemic was declared due to the spread of the novel coronavirus, SARS-CoV-2. Viral infection is caused by the interaction between the SARS-CoV-2 receptor binding domain (RBD) and the human ACE2 receptor (hACE2). Previous computational studies have identified repurposed small molecules that target the RBD, but very few have screened drugs in the RBD–hACE2 interface. When studies focus solely on the binding affinity between the drug and the RBD, they ignore the effect of hACE2, resulting in an incomplete analysis. We screened ACE inhibitors and previously identified SARS-CoV-2 inhibitors for binding to the RBD—hACE2 interface, and then conducted 500 ns of unrestrained molecular dynamics (MD) simulations of fosinopril, fosinoprilat, lisinopril, emodin, diquafosol, and physcion bound to the interface to assess the binding characteristics of these ligands. Based on MM-GBSA analysis, all six ligands bind favorably in the interface and inhibit the RBD–hACE2 interaction. However, when we repeat our simulation by first binding the drug to the RBD before interacting with hACE2, we find that fosinopril, fosinoprilat, and lisinopril result in a strongly interacting trimeric complex (RBD-drug-hACE2). Hydrogen bonding and pairwise decomposition analyses further suggest that fosinopril is the best RBD inhibitor. However, when lisinopril is bound, it stabilizes the trimeric complex and, therefore, is not an ideal potential drug candidate. Overall, these results reveal important atomistic interactions critical to the binding of the RBD to hACE2 and highlight the significance of including all protein partners in the evaluation of a potential drug candidate.

     
    more » « less
  2. Although COVID-19 transmission has been reduced by the advent of vaccinations and a variety of rapid monitoring techniques, the SARS-CoV-2 virus itself has shown a remarkable ability to mutate and persist. With this long track record of immune escape, researchers are still exploring prophylactic treatments to curtail future SARS-CoV-2 variants. Specifically, much focus has been placed on the antiviral lectin Griffithsin in preventing spike protein-mediated infection via the hACE2 receptor (direct infection). However, an oft-overlooked aspect of SARS-CoV-2 infection is viral capture by attachment receptors such as DC-SIGN, which is thought to facilitate the initial stages of COVID-19 infection in the lung tissue (called trans-infection). In addition, while immune escape is dictated by mutations in the spike protein, coronaviral virions also incorporate M, N, and E structural proteins within the particle. In this paper, we explored how several structural facets of both the SARS-CoV-2 virion and the antiviral lectin Griffithsin can affect and attenuate the infectivity of SARS-CoV-2 pseudovirus. We found that Griffithsin was a better inhibitor of hACE2-mediated direct infection when the coronaviral M protein is present compared to when it is absent (possibly providing an explanation regarding why Griffithsin shows better inhibition against authentic SARS-CoV-2 as opposed to pseudotyped viruses, which generally do not contain M) and that Griffithsin was not an effective inhibitor of DC-SIGN-mediated trans-infection. Furthermore, we found that DC-SIGN appeared to mediate trans-infection exclusively via binding to the SARS-CoV-2 spike protein, with no significant effect observed when other viral proteins (M, N, and/or E) were present. These results provide etiological data that may help to direct the development of novel antiviral treatments, either by leveraging Griffithsin binding to the M protein as a novel strategy to prevent SARS-CoV-2 infection or by narrowing efforts to inhibit trans-infection to focus on DC-SIGN binding to SARS-CoV-2 spike protein.

     
    more » « less
  3. Severe Acute respiratory syndrome coronavirus (SARS-CoV-1) attaches to the host cell surface to initiate the interaction between the receptor-binding domain (RBD) of its spike glycoprotein (S) and the human Angiotensin-converting enzyme (hACE2) receptor. SARS-CoV-1 mutates frequently because of its RNA genome, which challenges the antiviral development. Here, we per-formed computational saturation mutagenesis of the S protein of SARS-CoV-1 to identify the residues crucial for its functions. We used the structure-based energy calculations to analyze the effects of the missense mutations on the SARS-CoV-1 S stability and the binding affinity with hACE2. The sequence and structure alignment showed similarities between the S proteins of SARS-CoV-1 and SARS-CoV-2. Interestingly, we found that target mutations of S protein amino acids generate similar effects on their stabilities between SARS-CoV-1 and SARS-CoV-2. For example, G839W of SARS-CoV-1 corresponds to G857W of SARS-CoV-2, which decrease the stability of their S glycoproteins. The viral mutation analysis of the two different SARS-CoV-1 isolates showed that mutations, T487S and L472P, weakened the S-hACE2 binding of the 2003–2004 SARS-CoV-1 isolate. In addition, the mutations of L472P and F360S destabilized the 2003–2004 viral isolate. We further predicted that many mutations on N-linked glycosylation sites would increase the stability of the S glycoprotein. Our results can be of therapeutic importance in the design of antivirals or vaccines against SARS-CoV-1 and SARS-CoV-2. 
    more » « less
  4. The rapid spread of SARS-CoV-2 required immediate actions to control the transmission of the virus and minimize its impact on humanity. An extensive mutation rate of this viral genome contributes to the virus’ ability to quickly adapt to environmental changes, impacts transmissibility and antigenicity, and may facilitate immune escape. Therefore, it is of great interest for researchers working in vaccine development and drug design to consider the impact of mutations on virus-drug interactions. Here, we propose a multitarget drug discovery pipeline for identifying potential drug candidates which can efficiently inhibit the Receptor Binding Domain (RBD) of spike glycoproteins from different variants of SARS-CoV-2. Eight homology models of RBDs for selected variants were created and validated using reference crystal structures. We then investigated interactions between host receptor ACE2 and RBDs from nine variants of SARS-CoV-2. It led us to conclude that efficient multi-variant targeting drugs should be capable of blocking residues Q(R)493 and N487 in RBDs. Using methods of molecular docking, molecular mechanics, and molecular dynamics, we identified three lead compounds (hesperidin, narirutin, and neohesperidin) suitable for multitarget SARS-CoV-2 inhibition. These compounds are flavanone glycosides found in citrus fruits – an active ingredient of Traditional Chinese Medicines. The developed pipeline can be further used to (1) model mutants for which crystal structures are not yet available and (2) scan a more extensive library of compounds against other mutated viral proteins. 
    more » « less
  5. The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2020 has impacted daily life globally for over a year. While multiple vaccines have been authorized for emergency use and one oral medication has entered clinical trials, we are still seeking antiviral drugs for a long-term treatment for SARS-CoV-2 as well as other coronaviruses. Computational drug screenings of two SARS-CoV-2 protein target candidates are presented in this thesis: the nidoviral RNA uridylate-specific endoribonuclease (Nsp15) and the main protease (Mpro) of SARS-CoV-2. Nonstructural proteins of coronaviruses were selected as targets as they are more conserved across coronavirus strains than structural proteins. High throughput virtual screening of small molecule libraries including DrugBank and ZINC 15 resulted in several promising compounds for each of these targets. Molecular dynamics simulation allowed us to predict the binding energies for these compounds using molecular mechanics with generalized born surface area solvation calculations (MM-GBSA). Four top compounds were discovered for Nsp15 and eight compounds for Mpro. 
    more » « less