Computational characterization and design of SARS coronavirus receptor recognition and antibody neutralization

The sequential determination of crystal structures of the SARS coronavirus spike receptor-binding domain (RBD) in complex with its cellular receptor or neutralizing antibody opened a door for the design and development of antiviral competitive inhibitors. Based on those complex structures, we conduct computational characterization and design of RBD-mediated receptor recognition and antibody neutralization. The comparisons between computational predictions and experimental evidences validate our structural bioinformatics protocols. And the calculations predict a number of single substitutions on RBD, receptor or antibody that could remarkably elevate the binding affinities of those complexes. It is reasonable to anticipate our structure-based computation-derived hypotheses could be informative to the future biochemical and immunological tests.     

Binding of human ACE2 and RBD of Omicron enhanced by unique interaction patterns among SARS-CoV-2 variants of concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.     

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.     

Possible Transmission Flow of SARS-CoV-2 Based on ACE2 Features

Angiotensin-converting enzyme 2 (ACE2) is the cellular receptor for the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) that is engendering the severe coronavirus disease 2019 (COVID-19) pandemic. The spike (S) protein receptor-binding domain (RBD) of SARS-CoV-2 binds to the three sub-domains viz. amino acids (aa) 22–42, aa 79–84, and aa 330–393 of ACE2 on human cells to initiate entry. It was reported earlier that the receptor utilization capacity of ACE2 proteins from different species, such as cats, chimpanzees, dogs, and cattle, are different. A comprehensive analysis of ACE2 receptors of nineteen species was carried out in this study, and the findings propose a possible SARS-CoV-2 transmission flow across these nineteen species.     

Fruit Bromelain-Derived Peptide Potentially Restrains the Attachment of SARS-CoV-2 Variants to hACE2: A Pharmacoinformatics Approach

Before entering the cell, the SARS-CoV-2 spike glycoprotein receptor-binding domain (RBD) binds to the human angiotensin-converting enzyme 2 (hACE2) receptor. Hence, this RBD is a critical target for the development of antiviral agents. Recent studies have discovered that SARS-CoV-2 variants with mutations in the RBD have spread globally. The purpose of this in silico study was to determine the potential of a fruit bromelain-derived peptide. DYGAVNEVK. to inhibit the entry of various SARS-CoV-2 variants into human cells by targeting the hACE binding site within the RBD. Molecular docking analysis revealed that DYGAVNEVK interacts with several critical RBD binding residues responsible for the adhesion of the RBD to hACE2. Moreover, 100 ns MD simulations revealed stable interactions between DYGAVNEVK and RBD variants derived from the trajectory of root-mean-square deviation (RMSD), radius of gyration (Rg), and root-mean-square fluctuation (RMSF) analysis, as well as free binding energy calculations. Overall, our computational results indicate that DYGAVNEVK warrants further investigation as a candidate for preventing SARS-CoV-2 due to its interaction with the RBD of SARS-CoV-2 variants.     

Gold Metallodrugs to Target Coronavirus Proteins: Inhibitory Effects on the Spike-ACE2 Interaction and on PLpro Protease Activity by Auranofin and Gold Organometallics**

Gold complexes have a long tradition in medicine and for many examples antirheumatic, anticancer or anti-infective effects have been confirmed. Herein, we evaluated the lead compound Auranofin and five selected gold organometallics as inhibitors of two relevant drug targets of severe acute respiratory syndrome coronaviruses (SARS-CoV). The gold metallodrugs were effective inhibitors of the interaction of the SARS-CoV-2 spike protein with the angiotensin converting enzyme 2 (ACE2) host receptor and might thus interfere with the viral entry process. The gold metallodrugs were also efficient inhibitors of the papain-like protease (PLpro) of SARS-CoV-1 and SARS-CoV-2, which is a key enzyme in the viral replication. Regarding PLpro from SARS-CoV-2, the here reported inhibitors are among the very first experimentally confirmed examples with activity against this target enzyme. Importantly, the activity of the complexes against both PLpro enzymes correlated with the ability of the inhibitors to remove zinc ions from the labile zinc center of the enzyme. Taken together, the results of this pilot study suggest further evaluation of gold complexes as SARS-CoV antiviral drugs.     

A molecular docking model of SARS-CoV S1 protein in complex with its receptor, human ACE2

The exact residues within severe acute respiratory syndrome coronavirus (SARS-CoV) S1 protein and its receptor, human ACE2, involved in their interaction still remain largely undetermined. Identification of exact amino acid residues that are crucial for the interaction of S1 with ACE2 could provide working hypotheses for experimental studies and might be helpful for the development of antiviral inhibitor. In this paper, a molecular docking model of SARS-CoV S1 protein in complex with human ACE2 was constructed. The interacting residue pairs within this complex model and their contact types were also identified. Our model, supported by significant biochemical evidence, suggested receptor-binding residues were concentrated in two segments of S1 protein. In contrast, the interfacial residues in ACE2, though close to each other in tertiary structure, were found to be widely scattered in the primary sequence. In particular, the S1 residue ARG453 and ACE2 residue LYS341 might be the key residues in the complex formation.     

Prediction and mitigation of mutation threats to COVID-19 vaccines and antibody therapies

Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic. They, however, are prone to over 5000 mutations on the spike (S) protein uncovered by a Mutation Tracker based on over 200 000 genome isolates. It is imperative to understand how mutations will impact vaccines and antibodies in development. In this work, we first study the mechanism, frequency, and ratio of mutations on the S protein which is the common target of most COVID-19 vaccines and antibody therapies. Additionally, we build a library of 56 antibody structures and analyze their 2D and 3D characteristics. Moreover, we predict the mutation-induced binding free energy (BFE) changes for the complexes of S protein and antibodies or ACE2. By integrating genetics, biophysics, deep learning, and algebraic topology, we reveal that most of the 462 mutations on the receptor-binding domain (RBD) will weaken the binding of S protein and antibodies and disrupt the efficacy and reliability of antibody therapies and vaccines. A list of 31 antibody disrupting mutants is identified, while many other disruptive mutations are detailed as well. We also unveil that about 65% of the existing RBD mutations, including those variants recently found in the United Kingdom (UK) and South Africa, will strengthen the binding between the S protein and human angiotensin-converting enzyme 2 (ACE2), resulting in more infectious COVID-19 variants. We discover the disparity between the extreme values of RBD mutation-induced BFE strengthening and weakening of the bindings with antibodies and angiotensin-converting enzyme 2 (ACE2), suggesting that SARS-CoV-2 is at an advanced stage of evolution for human infection, while the human immune system is able to produce optimized antibodies. This discovery, unfortunately, implies the vulnerability of current vaccines and antibody drugs to new mutations. Our predictions were validated by comparison with more than 1400 deep mutations on the S protein RBD. Our results show the urgent need to develop new mutation-resistant vaccines and antibodies and to prepare for seasonal vaccinations.

Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic.     

Structural behavior of monomer of SARS-CoV-2 spike protein during initial stage of adsorption on graphene

Spike glycoprotein of the SARS-CoV-2 virus and its structure play a crucial role in the infections of cells containing angiotensin-converting enzyme 2 (ACE2) as well as in the interactions of this virus with surfaces. Protection against viruses and often even their deactivation is one of the great varieties of graphene applications. The structural changes of the non-glycosylated monomer of the spike glycoprotein trimer (denoted as S-protein in this work) triggered by its adsorption onto graphene at the initial stage are investigated by means of atomistic molecular dynamics simulations. The adsorption of the S-protein happens readily during the first 10 ns. The shape of the S-protein becomes more prolate during the adsorption, but this trend, albeit less pronounced, is observed also for the freely relaxing S-protein in water. The receptor-binding domain (RBD) of the free and adsorbed S-protein manifests itself as the most rigid fragment of the whole S-protein. The adsorption even enhances the rigidity of the whole S-protein as well as its subunits. Only one residue of the RBD involved in the specific interactions with ACE2 during the cell infection is involved in the direct contact of the adsorbed S-protein with the graphene. The new intramolecular hydrogen bonds formed during the S-protein adsorption replace the S-protein-water hydrogen bonds; this trend, although less apparent, is observed also during the relaxation of the free S-protein in water. In the initial phase, the secondary structure of the RBD fragment specifically interacting with ACE2 receptor is not affected during the S-protein adsorption onto the graphene.     

Identification of novel small-molecule inhibitors of severe acute respiratory syndrome-associated coronavirus by chemical genetics

The severe acute respiratory syndrome-associated coronavirus (SARS-CoV) infected more than 8,000 people across 29 countries and caused more than 900 fatalities. Based on the concept of chemical genetics, we screened 50,240 structurally diverse small molecules from which we identified 104 compounds with anti-SARS-CoV activity. Of these 104 compounds, 2 target the SARS-CoV main protease (M(pro)), 7 target helicase (Hel), and 18 target spike (S) protein-angiotensin-converting enzyme 2 (ACE2)-mediated viral entry. The EC(50) of the majority of the 104 compounds determined by SARS-CoV plaque reduction assay were found to be at low micromolar range. Three selected compounds, MP576, HE602, and VE607, validated to be inhibitors of SARS-CoV M(pro), Hel, and viral entry, respectively, exhibited potent antiviral activity (EC(50) < 10 microM) and comparable inhibitory activities in target-specific in vitro assays.     

Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics

Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.     

Native Structure-Based Peptides as Potential Protein–Protein Interaction Inhibitors of SARS-CoV-2 Spike Protein and Human ACE2 Receptor

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is a positive-strand RNA virus that causes severe respiratory syndrome in humans, which is now referred to as coronavirus disease 2019 (COVID-19). Since December 2019, the new pathogen has rapidly spread globally, with over 65 million cases reported to the beginning of December 2020, including over 1.5 million deaths. Unfortunately, currently, there is no specific and effective treatment for COVID-19. As SARS-CoV-2 relies on its spike proteins (S) to bind to a host cell-surface receptor angiotensin-converting enzyme-2(ACE2), and this interaction is proved to be responsible for entering a virus into host cells, it makes an ideal target for antiviral drug development. In this work, we design three very short peptides based on the ACE2 sequence/structure fragments, which may effectively bind to the receptor-binding domain (RBD) of S protein and may, in turn, disrupt the important virus-host protein–protein interactions, blocking early steps of SARS-CoV-2 infection. Two of our peptides bind to virus protein with affinity in nanomolar range, and as very short peptides have great potential for drug development.     

In Silico Discovery of Antimicrobial Peptides as an Alternative to Control SARS-CoV-2

A serious pandemic has been caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The interaction between spike surface viral protein (Sgp) and the angiotensin-converting enzyme 2 (ACE2) cellular receptor is essential to understand the SARS-CoV-2 infectivity and pathogenicity. Currently, no drugs are available to treat the infection caused by this coronavirus and the use of antimicrobial peptides (AMPs) may be a promising alternative therapeutic strategy to control SARS-CoV-2. In this study, we investigated the in silico interaction of AMPs with viral structural proteins and host cell receptors. We screened the antimicrobial peptide database (APD3) and selected 15 peptides based on their physicochemical and antiviral properties. The interactions of AMPs with Sgp and ACE2 were performed by docking analysis. The results revealed that two amphibian AMPs, caerin 1.6 and caerin 1.10, had the highest affinity for Sgp proteins while interaction with the ACE2 receptor was reduced. The effective AMPs interacted particularly with Arg995 located in the S2 subunits of Sgp, which is key subunit that plays an essential role in viral fusion and entry into the host cell through ACE2. Given these computational findings, new potentially effective AMPs with antiviral properties for SARS-CoV-2 were identified, but they need experimental validation for their therapeutic effectiveness.     

Identification of SARS-CoV-2 Receptor Binding Inhibitors by In Vitro Screening of Drug Libraries

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the coronavirus disease 2019 (COVID-19) global pandemic. The first step of viral infection is cell attachment, which is mediated by the binding of the SARS-CoV-2 receptor binding domain (RBD), part of the virus spike protein, to human angiotensin-converting enzyme 2 (ACE2). Therefore, drug repurposing to discover RBD-ACE2 binding inhibitors may provide a rapid and safe approach for COVID-19 therapy. Here, we describe the development of an in vitro RBD-ACE2 binding assay and its application to identify inhibitors of the interaction of the SARS-CoV-2 RBD to ACE2 by the high-throughput screening of two compound libraries (LOPAC^®1280 and DiscoveryProbeTM). Three compounds, heparin sodium, aurintricarboxylic acid (ATA), and ellagic acid, were found to exert an effective binding inhibition, with IC50 values ranging from 0.6 to 5.5 µg/mL. A plaque reduction assay in Vero E6 cells infected with a SARS-CoV-2 surrogate virus confirmed the inhibition efficacy of heparin sodium and ATA. Molecular docking analysis located potential binding sites of these compounds in the RBD. In light of these findings, the screening system described herein can be applied to other drug libraries to discover potent SARS-CoV-2 inhibitors.     

Atomistic simulation reveals structural mechanisms underlying D614G spike glycoprotein-enhanced fitness in SARS-COV-2

D614G spike glycoprotein (sgp) mutation in rapidly spreading severe acute respiratory syndrome coronavirus-2 (SARS-COV-2) is associated with enhanced fitness and higher transmissibility in new cases of COVID-19 but the underlying mechanism is unknown. Here, using atomistic simulation, a plausible mechanism has been delineated. In G614 sgp but not wild type, increased D(G)614-T859 Cα-distance within 65 ns is interpreted as S1/S2 protomer dissociation. Overall, ACE2-binding, post-fusion core, open-state and sub-optimal antibody-binding conformations were preferentially sampled by the G614 mutant, but not wild type. Furthermore, in the wild type, only one of the three sgp chains has optimal communication route between residue 614 and the receptor-binding domain (RBD); whereas, two of the three chains communicated directly in G614 mutant. These data provide evidence that D614G sgp mutant is more available for receptor binding, cellular invasion and reduced antibody interaction; thus, providing framework for enhanced fitness and higher transmissibility in D614G SARS-COV-2 mutant.     

Interaction of a peptide from the pre-transmembrane domain of the severe acute respiratory syndrome coronavirus spike protein with phospholipid membranes

The severe acute respiratory syndrome coronavirus (SARS-CoV) envelope spike (S) glycoprotein, a Class I viral fusion protein, is responsible for the fusion between the membranes of the virus and the target cell. In order to gain new insight into the protein membrane alteration leading to the viral fusion mechanism, a peptide pertaining to the putative pre-transmembrane domain (PTM) of the S glycoprotein has been studied by infrared and fluorescence spectroscopies regarding its structure, its ability to induce membrane leakage, aggregation, and fusion, as well as its affinity toward specific phospholipids. We demonstrate that the SARS-CoV PTM peptide binds to and interacts with phospholipid model membranes, and, at the same time, it adopts different conformations when bound to membranes of different compositions. As it has been already suggested for other viral fusion proteins such as HIV gp41, the region of the SARS-CoV protein where the PTM peptide resides could be involved in the merging of the viral and target cell membranes working synergistically with other membrane-active regions of the SARS-CoV S glycoprotein to heighten the fusion process and therefore might be essential for the assistance and enhancement of the viral and cell fusion process.     

Novel pyrazole-biphenyl-carboxamides for SARS-CoV2 entry-level restriction and microbial infections

Microbial diseases including viral infection are big issues globally. Effective medicinal discovery for them is the need for the day. In this study, we report pyrazole-biphenyl-carboxamides ( 4a-l ) validated for their SARS-CoV2 entry-level restriction effect over studying the protein–protein interaction of SARS-CoV2 with human ACE protein. Their extended antimicrobial properties were also evaluated. Online and offline software tools predicted MD simulation and ADMET druggability in silico. The antimicrobial efficacy of all compounds was also evaluated against Gram+^ve Streptococcus pneumoniae (MTCC 1936), Staphylococcus aureus (MTCC 737) and Gram-^ve Escherichia coli (MTCC 443), Pseudomonas aeruginosa (MTCC 424) (bacteria). In the results, compounds 4g and 4i were evenly active against both bacteria at a low concentration range (MIC: 1.00 to 9.5 μg/mL) and displayed lesser toxicity to tested mammalian cells (EC₁₀₀ = 75 μg/mL). Furthermore, it was able to kill metabolically inactive bacterial cells and eradicate established biofilms of Methicillin-resistant Staphylococcus aureus (MRSA). Both the compounds inhibited DNA gyrase well with an IC₅₀ 0.25 μM (96% relative activity) and 0.52 μM (97% relative activity) respectively. Compounds ( 4a-l ) showed restrictive efficiency of SARS-CoV2 spike protein (SC2SP) and human angiotensin-converting enzyme 2 (hACE2) entry-level association in COVID-19 in silico. To assess this ability, firstly, we identified the crucial amino acid residues involved in the interface of SARS-CoV-2 and hACE2 virtually. We recognized the ability of 4a-l binding to the binding interface to SARS-CoV2; thus, the interaction of SC2SP-hACE2 was effectively inhibited.     

Screening and identification of linear B-cell epitopes and entry-blocking peptide of severe acute respiratory syndrome (SARS)-associated coronavirus using synthetic overlapping peptide library

A 10-mer overlapping peptide library has been synthesized for screening and identification of linear B-cell epitopes of severe acute respiratory syndrome associated coronavirus (SARS-CoV), which spanned the major structural proteins of SARS-CoV. One hundred and eleven candidate peptides were positive according to the result of PEPscan, which were assembled into 22 longer peptides. Five of these peptides showed high cross-immunoreactivities (approximately 66.7 to 90.5%) to SARS convalescent patients' sera from the severest epidemic regions of the China mainland. Most interestingly, S(471-503), a peptide located at the receptor binding domain (RBD) of SARS-CoV, could specifically block the binding between the RBD and angiotensin-converting enzyme 2, resulting in the inhibition of SARS-CoV entrance into host cells in vitro. The study demonstrated that S(471-503) peptide was a potential immunoantigen for the development of peptide-based vaccine or a candidate for further drug evaluation against the SARS-CoV virus-cell fusion.     

Comparative study of synonymous codon usage variations between the nucleocapsid and spike genes of coronavirus, and C-type lectin domain genes of human and mouse

Coronaviruses (CoVs) are single-stranded RNA viruses which contain the largest RNA genomes, and severe acute respiratory syndrome coronavirus (SARS-CoV), a newly found group 2 CoV, emerged as infectious disease with high mortality rate. In this study, we compared the synonymous codon usage patterns between the nucleocapsid and spike genes of CoVs, and C-type lectin domain (CTLD) genes of human and mouse on the codon basis. Findings indicate that the nucleocapsid genes of CoVs were affected from the synonymous codon usage bias than spike genes, and the CTLDs of human and mouse partially overlapped with the nucleocapsid genes of CoVs. In addition, we observed that CTLDs which showed the similar relative synonymous codon usage (RSCU) patterns with CoVs were commonly derived from the human chromosome 12, and mouse chromosome 6 and 12, suggesting that there might be a specific genomic region or chromosomes which show a more similar synonymous codon usage pattern with viral genes. Our findings contribute to developing the codon-optimization method in DNA vaccines, and further study is needed to determine a specific correlation between the codon usage patterns and the chromosomal locations in higher organisms.     

Computational Insights into the Allosteric Effect and Dynamic Structural Features of the SARS-COV-2 Spike Protein

COVID-19 caused by SARS-COV-2 is continuing to surge globally. The spike (S) protein is the key protein of SARS-COV-2 that recognizes and binds to the host target ACE2. In this study, molecular dynamics simulation was used to elucidate the allosteric effect of the S protein. Binding of ACE2 caused a centripetal movement of the receptor-binding domain of the S protein. The dihedral changes in Phe329 and Phe515 played a key role in this process. Two potential cleavage sites S1/S2 and S2′ were exposed on the surface after the binding of ACE2. The binding affinity of SARS-COV-2 S protein and ACE2 was higher than that of SARS-COV. This was mainly due to the mutation of Asp480 in SARS-COV to Ser494 in SARS-COV-2, which greatly weakened the electrostatic repulsion. The result provides a theoretical basis for the SARS-COV-2 infection and aids the development of biosensors and detection reagents.