Screening S protein - ACE2 blockers from natural products: Strategies and advances in the discovery of potential inhibitors of COVID-19.

The Coronavirus disease, 2019 (COVID-19) is caused by severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2), which poses a major threat to human life and health. Given its continued development, limiting the spread of COVID-19 in the population remains a challenging task. Currently, multiple therapies are being tried around the world to deal with SARS-CoV-2 infection, and a variety of studies have shown that natural products have a significant effect on COVID-19 patients. The combination of SARS-CoV-2 S protein with Angiotensin converting enzyme II(ACE2) of host cell to promote membrane fusion is an initial critical step for SARS-CoV-2 infection. Therefore, screening natural products that inhibit the binding of SARS-CoV-2 S protein and ACE2 also provides a feasible strategy for the treatment of COVID-19. Establishment of high throughput screening model is an important basis and key technology for screening S protein-ACE2 blockers. Based on this, the molecular structures of SARS-CoV-2 and ACE2 and their processes in the life cycle of SARS-CoV-2 and host cell infection were firstly reviewed in this paper, with emphasis on the methods and techniques of screening S protein-ACE2 blockers, including Virtual Screening (VS), Surface Plasmon Resonance (SPR), Biochromatography, Biotin-avidin with Enzyme-linked Immunosorbent assay and Gene Chip Technology. Furthermore, the technical principle, advantages and disadvantages and application scope were further elaborated. Combined with the application of the above screening technologies in S protein-ACE2 blockers, a variety of natural products, such as flavonoids, terpenoids, phenols, alkaloids, were summarized, which could be used as S protein-ACE2 blockers, in order to provide ideas for the efficient discovery of S protein-ACE2 blockers from natural sources and contribute to the development of broad-spectrum anti coronavirus drugs.

Introduction

COVID-19 is the most widespread global pandemic to hit humans in the last century, with more than 165 million people infected worldwide as of May 21, 2021. Because of its high infectivity and mortality, it has been declared as an international public health emergency by the World Health Organization, which has brought a major threat to human life safety and health, and led to a huge economic and social burden to the world [1,2]. COVID-19 is caused by SARS-CoV-2, which is structurally related to two other coronaviruses that cause SARS and Middle East respiratory syndrome (MERS) [3]. One of the main ways to prevent SARS-CoV-2 infection is to develop vaccines. At present, major breakthroughs have been made in the development of safe and effective vaccines against SARS-CoV-2 to prevent further spread of COVID-19, with a total of 201 candidate vaccines under development worldwide, 45 of which are in clinical trials and 156 in preclinical studies [4]. However, mutations in the virus have resulted in changes in the infectivity and pathogenicity of SARS-CoV-2. The mutated Alpha strain, which first circulated in the United Kingdom, increased the transmissibility by about 56%. Subsequently, Beta strain circulating in South Africa also increased the transmissibility, and Gamma strain circulating in South America increased the transmissibility by 1.4–2.2 times [5]. Recently, the mutated new coronavirus with at least two key point mutations was found in India and spread to many countries and regions in the world. Virology experts widely believe that the mutation occurred in a key part of SARS-CoV-2, making it significantly more infectious and potentially making existing vaccines less effective [6]. Delta strain that is currently circulating around the world has the characteristics of rapid transmission, certain immune escape ability, and may aggravate the disease, so its threat has attracted the attention of the whole world. Delta strain has multiple mutations in the spike protein that give it both greater ACE2 receptor binding ability and immune escape ability. Delta strain grows more rapidly in the lung and throat of infected people and has a thousand-fold increase in viral load levels, most likely associated with the ACE2 receptor [7]. In view of its continued development, limiting the spread of COVID-19 in population remains a challenging task. Although several highly effective vaccines have recently been widely deployed around the world, widespread adoption and equitable distribution of vaccines face challenges in the context of global inequality in wealth and education. Therefore, we need to pay attention not only to the development of vaccines, but also to the application of therapeutic drugs. (see Table 1 )

Table 1

Natural products inhibiting S protein to bind to ACE2.

Compounds	Chemical class	Source	Effect	Screening method	Reference
Quercetin (1)	Flavonoids	Bupleurum chinense DC.Cyathula officinalis Kuan	Binds to ACE2, impairs the interaction between S protein and ACE2	Net-work pharmacology, Molecular docking, SPR	[123]
Quercetin (1)	Flavonoids	In silico	Targets ACE2 expression and alters the expression of genes encoding protein targets of SARS-CoV-2 in human cells	GSEA, expression profiling experiments	[115]
Quercetin (1)	Flavonoids	In silico	Binds to ACE2 and MPro	Molecular docking, Net-work pharmacology	[62]
Quercetin (1)	Flavonoids	Forsythia suspensa (Thunb.) Vahl Lonicera japonica Thunb. Bupleurum chinense DC.	Binds to the RBD region of S protein and inactivate it and then prevent S protein binding to ACE2 of epithelial cell surface	Net-work pharmacology, SPR, Molecular docking	[86]
Quercetin (1)	Flavonoids	In silico	Decreases ACE2 expression via regulation of transcription factors or miRNAs of ACE2	Net-work pharmacology	[155]
Isorhamnetin (2)	Flavonoids	In silico	Bind to ACE2 and MPro	Molecular docking, Net-work pharmacology	[62]
Isorhamnetin (2)	Flavonoids	Hippophae rhamnoides L.	Binds to S protein and ACE2, inhibites SARS-CoV-2 spike pseudotyped virus entering ACE2hcell	CMC, Molecular docking, SPR	[130]
Puerarin (3)	Flavonoids	Pueraria lobata (Willd.) Ohwi	Binds to ACE2 and S protein, impairs the interaction between S protein and ACE2	Net-work pharmacology, Molecular docking, SPR	[123]
Puerarin (3)	Flavonoids	In silico	Binds to ACE2	Molecular docking	[156]
Rutin (4)	Flavonoids	Carthamus tinctorius L.	Binds to S protein and ACE2	Molecular docking	[157]
Rutin (4)	Flavonoids	Bupleurum chinense DC.Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao	Binds to ACE2	Net-work pharmacology, SPR, Molecular docking	[158]
Oroxylin A (5)	Flavonoids	Scutellaria baicalensis Georgi	Binds to ACE2, suppresses entrance of SARS-CoV-spiked pseudo typed virus into ACE2 cells, inhibits LC3-mediated autophagy of ACE2 cells	CMC, Molecular docking, SPR	[100]
Procyanidin (6)	Flavonoids	In silico	Binds to ACE2 and S protein	Molecular docking, Molecular dynamics simulations	[124]
Irisolidone (7)	Flavonoids	In silico	Binds to ACE2	Molecular docking	[156]
Luteolin (8)	Flavonoids	In silico	Binds to ACE2 and MPro	Molecular docking, Network pharmacology	[62]
Hesperidin (9)	Flavonoids	Valeriana Jatamansi	Binds to S protein and ACE2, destabilize S protein binding to human host ACE2 receptor	Molecular docking, Molecular dynamics simulations	[125]
Chrysin (10)	Flavonoids	Oroxylum indtcum (L.)Vent.	Binds to S protein and ACE2, destabilize S-protein binding to human host ACE2 receptor	Molecular docking, Molecular dynamics simulations	[125]
Anhydrosafflor yellow B (11)	Flavonoids	Carthamus tinctorius L.	Binds to S protein and ACE2	Molecular docking	[157]
Glabridin (12)	Flavonoids	In silico	Decreases ACE2 expression via regulation of transcription factors or miRNAs of ACE2	Net-work pharmacology	[155]
Myricitrin (13)	Flavonoids	Myristica fragrans Houtt.	Strong binding affinity to ACE2 and RNA dependent RNA polymerase	Molecular docking	[158]
Euchrenone (14)	Flavonoids	Glycyrrhiza uralensis Fisch	Binds to ACE2, Mpro and RdRp	Net-work pharmacology, Molecular docking	[66]
Epigallocatechin-3- gallate (15)	Flavonoids	Camellia sinensis (L.) O. Kuntze	Binds to S protein to impair the interaction between S protein and ACE2	Molecular docking, in vitro antiviral experiment	[128,129]
Glycyrrhizic acid (16)	Terpenoids	Glycyrrhiza uralensis Fisch	Binds to ACE2, Mpro and RdRp	Net-work pharmacology, Molecular docking	[66]
Glycyrrhizic acid (16)	Terpenoids	In silico	Binds to RBD domain of the SARS-CoV-2 S protein	Molecular docking, Molecular dynamics simulations	[139]
Glycyrrhizic acid (16)	Terpenoids	In silico	Binds to ACE2	Molecular docking, Molecular dynamics simulations	[60]
Glycyrrhizic acid (16)	Terpenoids	Glycyrrhiza uralensis Fisch	Binds to S protein and ACE2, destroys the interaction between S-RBD and ACE2	SPR, Molecular docking	[140]
Uncaric acid (17)	Terpenoids	Uncaria tomentosa	Binds to RBD/ACE-2 interface and the ACE2 binding site on SARS-CoV-2 RBD viral spike	Molecular docking, Molecular dynamics simulations	[159]
Ursolic acid (18)	Terpenoids	In silico	Binds to ACE2	Molecular docking, Molecular dynamics simulations	[60]
Demethylzeylasteral (19)	Terpenoids	Tripterygium wilfordii Hook.f.	Binds to ACE2 and S protein, impairs the interaction between S protein and ACE2	SPR, Molecular docking	[85]
Maslinic acid (20)	Terpenoids	In silico	Binds to ACE2	Molecular docking, Molecular dynamics simulations	[60]
Obacunone (21)	Terpenoids	In silico	Binds to ACE2	Molecular docking, Molecular dynamics simulations	[60]
Andrographolide (22)	Terpenoids	In silico	Binds to ACE2	Molecular docking, Molecular dynamics simulations	[160]
Atractylenolide III(23)	Terpenoids	Atractylodes lancea (Tunb.) Dc. (Cangzhu)	Binds to ACE2 and has anti-inflammatory effects and antiviral effects in vitro	Molecular docking	[161]
Astragaloside IV(24)	Terpenoids	Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao	Binds to ACE2	Net-work pharmacology, SPR, Molecular docking	[86]
Taraxerol (25)Daturaolone (26)	Terpenoids	Clerodendrum D. innoxia	Bind to ACE2, S protein and Mpro	Molecular docking, Molecular dynamics simulations	[77]
Limonin (27)	Terpenoids	In silico	Binds to ACE2 to impair the interaction between S protein and ACE2	Molecular docking, Molecular dynamics simulations	[162]
Cucurbitacin G 2-glucoside (28)	Terpenoids	Cucurbita pepo L.	Binds to ACE2	Molecular docking	[163]
Citronellol (29)	Terpenoids	Pelargonium graveolens L'Hér.	Downregulates ACE2 expression in epithelial cells	Spectrum-effect relationship analysis	[119]
Limonene (30)	Terpenoids	Lemon essential oils	Downregulates ACE2 expression in epithelial cells	Spectrum-effect relationship analysis	[119]
Salvianolic acid A (31),Salvianolic acid B (32), Salvianolic acid C (33)	Phenols	Salvia miltiorrhiza Bge.	Bind to S protein and ACE2	Molecular docking, SPR	[146]
Neochlorogenic acid (34)	Phenols	Lonicera japonica Thunb.	Binds to ACE2 and inhibits ACE2 activity, binds to the contact surface of ACE2 and spike complex	MBC, SPR, Molecular docking	[95]
Kobophenol A (35)	Phenols	Caragana sinica (Buc'hoz) Rehder	Inhibits SARS-CoV-2 binding to cells through blocking spike RBD to the host ACE2 receptor	Molecular docking, Molecular dynamics simulations, ELISA	[107]
Resveratrol (36)	Phenols	In silico	Binds to ACE2	Molecular docking	[150]
Chlorogenic acid (37)	Phenols	Lonicera japonica Thunb.	Binds to ACE2 to impair the interaction between S-protein and ACE2	Molecular docking, Network pharmacology	[132]
Gallic acid (38)	Phenols	In silico	Decreases ACE2 expression via regulation of transcription factors or miRNAs of ACE2	Net-work pharmacology	[155]
Bis-demethoxycur-cumin (39)	Phenols	Curcuma longa L.	Binds to RBD domain of the SARS-CoV-2 S-protein	Molecular docking, Molecular dynamics simulations	[164]
δ-viniferin (40)	Phenols	Syzygium jambos (L.) Alston	Strong binding affinity to ACE-2 and RNA dependent RNA polymerase	Molecular docking	[158]
Pterostilbene (41)	Phenols	In silico	Binds to ACE2	Molecular docking, Molecular dynamics simulations	[160]
Ephedrine (42), Pseudoephedrine (43),Methylephedrine (44)	Alkaloids	Ephedra sinica Stapf	Bind to ACE2 and spike RBD, inhibit SARS-CoV-2 spike pseudovirus entering ACE2hcells	CMC-HPLC-IT-TOF-MS, SPR,Molecular docking	[99]
Cryptospirolepine (45)Cryptoquindoline (46)	Alkaloids	Cryptolepis sanguinolenta (Lindl.) Schlt (Periplocaceae)	Bind to ACE2 and S-protein, destroy the stability of ACE2-RBD complex	Molecular docking, Molecular dynamics simulations	[150]
Speciophylline (47),Uncarine F (48)	Alkaloids	Uncaria rhynchophylla (Miq.) Miq.	Bind to RBD/ACE-2 interface and the ACE2 binding site on SARS-CoV-2 RBD viral spike	Molecular docking, Molecular dynamics simulations	[159]
Solanine (49)	Alkaloids	Solanum Sps.	Inhibits spike RBD and main protease	Molecular docking combinedMolecular dynamics simulations	[27]
Bicuculline (50)	Alkaloids	In silico	Binds to ACE2	Molecular docking	[156]
Anisodamine(51)	Alkaloids	Anisodus tanguticus (Maxim.) Pascher.	Binds to ACE2	Net-work pharmacology, Molecular docking	[165]
Thymoquinone (52)	Quinones	Nigella sativa L.	Covers active sites of ACE2 and inhibits the activity of ACE2	Molecular docking	[166]
Emodin (53)	Quinones	Rheum palmatum L.	Binds to S-protein and ACE2, destabilizes S-protein binding to human host ACE2 receptor	Molecular docking, Molecular dynamics simulations	[125]
Rhein (54)	Quinones	Rheum palmatum L.	Binds to ACE2 and inhibits ACE2 activity, binds to the contact surface of ACE2 and spike complex	MBC, SPR, Molecular docking	[95]
Forsythoside A (55),Forsythoside I (56)	Glycosides	Forsythia suspensa (Thunb.) Vahl	Bind to ACE2 and inhibits ACE2 activity, bind to the contact surface of ACE2 and spike complex	MBC, SPR, Molecular docking	[95]
Amygdalin (57)	Glycosides	Prunus armeniaca L. var. ansu Maxim.	Binds to ACE2, Mpro and RdRp	Net-work pharmacology, Molecular docking	[66]
Urso-deoxycholic acid (58)	Steroids	Ipomoea obscura (L.)	Binds to ACE2 and MPro	Molecular docking, Molecular dynamics simulations	[167]
Withanolide A (59)	Steroids	Withania somnifera	Binds to ACE2, S-protein and Mpro	Molecular docking, Molecular dynamics simulations	[77]
β-sitosterol (60)	Steroids	Urtica dioica	Binds to ACE2	Molecular docking	[168]
Sesamin (61)Sesaminol (62)Sesamolin (63)	Lignans	Sesamum indicum L.	Bind to ACE2 and S-protein	Molecular docking	[61]
Phillyrin (64)	Lignans	Forsythia suspense (Thunb.) Vahl	Binds to ACE2 to impair the interaction between S-protein and ACE2	Molecular docking, Network pharmacology	[132]
Farnesiferol B (65)	Coumarins	Ferula sinkiangensisK. M. Shen	Binds to ACE2 and S-protein	Molecular docking	[61]

Natural products, widely found in nature, have long been an important source of innovative drug discovery. According to statistics, more than 60% of the drugs on the market are obtained from natural products or natural products as lead compounds after structural optimization and derivation [8]. Natural products are widely used in the treatment of antiviral infection and enhancement of host immune response. In previous coronavirus infections including SARS-CoV or MERS-CoV, natural products have significant therapeutic effects. Natural products also play an important role in the control of COVID-19, as SARS-CoV-2 is similar to SARS-CoV or MERS-CoV in certain genomic sequences [9]. For patients with mild symptoms, early intervention of natural products can effectively prevent the deterioration of the disease. For patients who have developed to advanced stage, natural products can improve symptoms, reduce many complications and mortality [10]. In addition, some natural products with antiviral activity come from a variety of spices, fruits and vegetables, which can reduce the risk and severity of various viral infections by antiviral and enhancing immune response [11]. Therefore, natural products are not only expected to be used to develop effective and less toxic anti SARS-CoV-2 drugs, but also can be used as an alternative supplement for people at risk of COVID-19 infection.

Studies have shown that ACE2 is the key target of SARS-CoV-2 entering host cells and the key point of drug intervention [12]. A trimer fusion protein S protein in SARS-CoV-2 structure is the main target of vaccine design [13], it is responsible for attaching the virus to the receptor of host cell surface, especially ACE2 [14]. ACE2 is widely expressed in the cell membranes of a variety of epithelial cells, including the lung, respiratory tract, heart, kidney, liver, pancreas, and intestine [15]. Receptor-binding domain (RBD) in S protein is a fragment of SARS-CoV-2 that binds to the ACE2 receptor of host cells and promotes the receptor-mediated endocytosis of the virus within host cells [16]. With the disclosure of ACE2 as the target of SARS-CoV-2 invasion into human body, blocking or competitively inhibiting the binding of SARS-CoV-2 and ACE2, so as to disrupt the cell process of virus infection, has become one of the therapeutic directions. Therefore, screening active ingredients from natural products that can inhibit the binding of SARS-CoV-2 and ACE2 also provides a feasible strategy for the treatment of COVID-19. It is worth noting that ACE2 located on the surface of host cells can also bind effectively with S protein of other coronavirus due to the genetic similarity between coronavirus [17]. Therefore, the active components which inhibit the binding of SARS-CoV-2 and ACE2 may also be effective for other types of coronaviruses, and can be used to fight other emerging coronavirus infection. Some of the natural products showed antiviral activity against different viruses in vitro and in vivo [18], so it is more conducive to the development of broad-spectrum anti-coronavirus drugs. At the same time, the establishment of a simple, rapid, sensitive, economic and efficient drug screening model is also an important basis and key technology for screening out active ingredients from natural products that block or competitively inhibit the binding of SARS-CoV-2 and ACE2. To this end, this paper firstly reviews the molecular structure of SARS-CoV-2, ACE2 and their roles in the life cycle of SARS-CoV-2 as well as the infection of host cells, and focuses on the research progress of screening methods and technologies for SARS-CoV-2 and ACE2 binding inhibitors. Finally, the active natural products that inhibit the binding of SARS-CoV-2 and ACE2 were summarized, in order to provide ideas for the efficient discovery of natural S protein-ACE2 blockers, and contribute to the development of broad-spectrum anti-coronavirus drugs.

The characteristics of SARS-CoV-2 and ACE2

SARS-CoV-2

Coronavirus belongs to the genus Coronavirinae, which contains four genera, α, β, δ and γ. Its genome consists of a single-stranded position-sense RNA with a length of 27–32 kb, which is the largest RNA virus identified so far [19,20]. SARS-CoV-2 is the seventh coronavirus of the genus β that can be transmitted from person to person. Sequencing results showed that the genome of SARS-CoV-2 was composed of 29,903 nucleotide bases and 10 open reading frames (ORFs), with a size of about 29.9 kb [21]. ORF1ab encodes large replicase polyprotein PP1ab, which cleaved Papain Like protease (Plpro) and 3-chymotrypsin Like protease (3Clpro) into nonstructural proteins (NSP1-16), ORF2-10 encodes four structural proteins: spike (S), envelope (E), membrane (M), and nucleocapsid (N) [22]. Among these structural proteins, S protein plays the most critical role in viral attachment, fusion and entry into target cells [23].

When imaged by cryo-electron microscopy, SARS-CoV-2 appears as approximately spherical particles with a diameter of about 100 nm. The surface of SARS-CoV-2 was covered with dense viral substance by a prominent S protein trimer lipid bilayer [24]. S protein is a trimeric fusion protein of SARS-CoV-2, which exists in a metastable prefusion conformation after translation [25]. Each monomer contains two functional subunits, S1 and S2. S1 subunit is exposed to the surface of the envelope and is important site for recognizing and binding to the receptor on the surface of host cells, including an N-terminal domain and RBD. SARS-CoV-2 binds to ACE2 of host cells through RBD. S2 subunit is the key structure of virus invasion and fusion between virus envelope and host cell membrane [26,27], which is embedded in the envelope and contains an N-terminal hydrophobic fusion peptide, two heptad repeats (HR1 and HR2)、a transmembrane domain and a cytoplasmic tail [28]. Due to the mutation of the RBD amino acid sequence of SARS-CoV-2 S protein, the Q493 and P499 amino acid residues are more stable to ACE2 [28,29], and the fusion core structure formed by the HR1 and HR2 domains in the S2 subunit is highly stable [30]. Therefore, in the process of virus adsorption and entry, SARS-CoV-2 is more capable than SARS-CoV. At the same time, S protein has multiple N-glycosylation sites, glycosyl is linked to protein through covalent bonds to form glycoprotein, and the existence of a large number of glycosyl can change the spatial structure of protein molecules through glycosylation and block or destroy antigen epitopes, thereby inhibiting the immune response of the body and playing a protective role on the virus [31]. Therefore, SARS-CoV-2 spreads from person to person at a high rate of infection.

The key factor of SARS-CoV-2 entering host cells——ACE2

In 2000, ACE2 attracted wide attention as the first reported ACE homolog. ACE2 is also the anti-regulator of renin-angiotensin system, which plays a key role in the cardiovascular system [32]. Meanwhile, ACE2, as the receptor for coronavirus to enter host cells, plays a certain role in lung diseases. In 2003, it was confirmed that SARS-CoV entered host cells by binding to ACE2 through the S protein on its surface, so ACE2 was confirmed as the SARS-CoV receptor [33]. SARS-CoV has 79.6% sequence identity with SARS-CoV-2 [34], and several studies have confirmed that ACE2 is also a SARS-CoV-2 receptor [[34], [35], [36]].

ACE2 is widely expressed in the cell membranes of a variety of epithelial cells, including the lung, respiratory tract, heart, kidney, liver, pancreas and intestine, with the highest distribution in the small intestine. Some COVID-19 patients show digestive tract symptoms such as vomiting and diarrhea in the early stage, which may also be related to the distribution of ACE2 in intestinal tissues [15,37]. Analysis of RNA sequencing data of human lung cells revealed that 83% of ACE2 was expressed in the type II apical surface of epithelial cells (AT2) [38]. Therefore, intestinal tissue and AT2 may be the main targets of SARS-CoV-2 infection. ACE2 is a transmembrane protein, which is difficult to exist stably in vitro. Yan et al. [37] first found that ACE2 could form a stable complex with an amino acid transporter B0AT1 in the intestine by using cryo-electron microscopy. By analyzing the full-length protein structure of ACE2, it was found that ACE2 existed in the form of dimer and had two conformational changes of open and close. However, both conformations contain interfaces for mutual recognition with SARS-CoV-2. In terms of function and distribution, there are two forms of ACE2 (mACE2 and sACE2). mACE2 is located on the cell membrane and contains a transmembrane anchor and an extracellular domain. The extracellular domain has been shown to be a receptor site for SARS-CoV-2 [21]. sACE2 is a soluble form that lacks a membrane anchor and circulates in the blood in small amounts [39]. Some studies have shown that sACE2 can inhibit the binding between SARS-CoV-2 and mACE2, so high level of sACE2 may have a protective effect on the body [40]. However, a recent study confirmed that sACE2 can interact with S protein of SARS-CoV-2, and the formed sACE2-S protein complex can enter host cells through AT1 receptor mediated endocytosis [41]. Therefore, sACE2 also plays an important role in SARS-CoV-2 infection.

The process of S protein binding to ACE2

SARS-CoV-2 mainly enters cells through endocytosis. Since S protein belongs to type I virus fusion protein, it needs protease cleavage to activate its fusion potential [42,43]. There are two S protein hydrolysis sites: one located at the S1/S2 junction and the other located inside S2 (S2’). S protein hydrolyzed into S1 and S2 subunits and connected through non-covalent bonds to obtain the ability of membrane fusion [44], which is crucial for its mediated SARS-CoV-2 infection. Therefore, the activation of S protein is an important step for SARS-CoV-2 to enter host cells. Transmembrane protease serine 2 (TMPRSS2) is a type II transmembrane serine protease, which is widely expressed in the epithelial cells of the respiratory tract, gastrointestinal tract and urogenital tract. It cleans S protein mainly through the S2’ site and is the key enzyme to activate S protein [45,46]. In addition, the insertion of a “RRAR” sequence recognized by furin enzyme at S1/S2 is a unique structure of SARS-CoV-2 S protein. The experiment verified that the furin enzyme cleaves S protein through S1/S2 sites, inhibiting the furin enzyme, TMPRSS2 or simultaneously inhibiting the two proteases will affect the fusion activity of S protein. Therefore, TMPRSS2 and furin enzyme play an important role in the activation of S protein [47,48].

As mentioned before, S Protein exists in a metastable trimer conformation before activation. RBD of S1 subunit exists in two states, one is “up” conformation with poor stability, the other is stable “down” conformation, which exposes and hides receptor binding sites respectively [28]. The structural study of S protein showed that RBD in the “down” conformation was closer to the center of trimer, with fewer exposed sites, and it was not easy to recognize ACE2 [37]. Thus only “up” RBD can bind ACE2, while “down” RBD is in the receptor inactive binding conformation [28,37]. Studies have shown that ACE2 binds to SARS-CoV-2 S protein with an affinity of about 15 nmol/L, which is about 10–20 times higher than that of SARS-CoV S protein [49]. When RBD binds ACE2 of the host cell with a high affinity, the virus attaches to the cell surface and destroys the stability of the trimer before S protein fusion. S protein cleavage usually occurs sequentially, with furin enzyme cleavage at S1/S2 site first, and followed by TMPRSS2 cleavage at S2’ site [48,50]. Under protease cleavage, the S1 (RBD)-ACE2 complex was separated from the S2 subunit, and the S2 subunit changed dramatically to a stable fusion conformation. The interaction between HR-1 and HR-2 in the S2 subunit formed a stable six-helix bundle structure. These changes expose the FP of the S2 subunit and insert it into the host membrane to initiate the fusion process [51]. Once viral RNA is released into the host cell, ss-positive sense RNA acts as messenger RNA and is translated through the host cells, and the life cycle of SARS-CoV-2 is shown in Fig. 1 .

Fig. 1

The life cycle of SARS-CoV-2. When SARS-CoV-2 S protein binds ACE2 of the host cell with a high affinity, the virus attaches to the cell surface and destroys the stability of the trimer before S protein fusion. S protein cleavage usually occurs sequentially, with furin enzyme cleavage at the S1/S2 site first, followed by TMPRSS2 cleavage at the S2’ site. Under protease cleavage, the S1 (RBD)–ACE2 complex was separated from the S2 subunit, and SARS-CoV-2 enters cells through endocytosis. SARS-CoV-2 releases the single chain positive RNA, then translates into pp1a and pp1ab with the help of ribosomes of host cells, and produces hydrolytic enzymes that can precisely cut polyprotein by means of self-shearing, namely main protease (Mpro) and papain-like protease (Plpro). Under the hydrolysis of Mpro and Plpro, RNA dependent RNA polymerase (RdRp) was formed. RdRp synthesized RNA genome (−) and then synthesized virus genome through genome replication. N protein binds to genomic RNA and replicates, transcribes and synthesizes in the cytoplasm. S, M and E proteins integrate into the membrane of the endoplasmic reticulum (ER). Virus are transported to the host cell membrane and released by exocytosis, and then infect other tissues and cells.

The binding of S protein to ACE2 is the first step of viral infection, in order to prevent SARS-CoV-2 from entering host cells, ACE2 has become one of the ideal targets for the development of novel broad-spectrum anti-coronavirus drugs. Firstly, the natural products that can bind to ACE2 were screened. On the one hand, their affinity for binding to other proteins was reduced biologically, and on the other hand, available sites were saturated to prevent viruses from binding to ACE2 and entering host cells. Secondly, similar to SARS-CoV, targeting SARS-CoV-2 S protein can also interfere with the virus. By binding to S protein through RBD and destroying S protein which directly binds to ACE2, the interaction between S protein and host cells can be prevented and infection can be prevented.

Screening methods for S protein-ACE2 blockers

Virtual screening

Virtual screening (VS) is a common strategy to obtain small molecule active compounds with high affinity for target proteins before biological activity test by preliminary screening multiple target molecules in the database of compounds with the help of computer in the early stage of drug development [52]. It is mainly divided into Structure-Based Virtual Screening (SBVS) and Ligand-Based Virtual Screening (LBVS) [53]. SBVS is primarily used to predict the optimal interaction mode or the formation of stable complexes between two molecules, and to evaluate the force of non-covalent interactions between ligands and molecular targets using the scoring function [54]. LBVS is based on the known structural information and physical and chemical properties of active compounds, and then a new active compound is selected according to the similarity principle [53]. At present, the research foundation of S protein-ACE2 blockers is weak, and there are few drugs with definite clinical efficacy and clinical trials. Therefore, the application of LBVS method is limited to a certain extent. Most studies on natural products as S protein-ACE2 blockers have focused on SBVS.

Molecular docking

SBVS requires to know the 3D structure of the target protein, so that the interaction between the target and each compound can be predicted by computer [55]. Among the SBVS methods, molecular docking has attracted wide attention because of its low cost and good results. Molecular docking is a new drug design method based on the receptor characteristics and the interaction between receptor and drug molecules. As a new research method combining physical and chemical principles with scientific calculation algorithm, it provides a feasible strategy for exploring natural product as S protein-ACE2 blockers [56]. In the process of molecular docking for drug screening, the scoring function can be used to determine the binding site of ligand and the conformation between target and ligand to find out the allosteric site, and can also be used to predict the binding affinity between target and ligand [57]. According to the affinity of the evaluated ligands to the receptor sites, the ligands with some pharmacological activities can be identified more accurately. It is worth noting that even if the compound can specifically bind to the target, it does not mean that it can play a pharmacological role in vivo. Only when the compounds are absorbed by the body and have certain bioavailability to avoid being metabolized and excreted, can they play a good role in vivo. Therefore, in the process of molecular docking, the ADME (absorption, distribution, metabolism, and exception) properties of compounds should be evaluated on the computer with the help of mathematical algorithm [58], At the same time, the most classical Lipinski's rule can also be used to predict the bioavailability of compounds [59]. The process of molecular docking is shown in Fig. 2 .

Fig. 2

Molecular docking process. Cryo-EM and X-ray reveal the structure of SARS-CoV-2 and ACE2. Downloading the structure and manually removing all binding ligands, ions and solvent molecules from the protein database for parameterization, and constructing 3D structure of targets. The ligands are mainly obtained from TCMSP, Drug Bank, ZINC, SANC and other databases. Then, Autodock tools is used to parameterize the structure to add complete hydrogen to the ligands. Finally, the ligands were subjected to molecular docking with SARS-CoV-2 and ACE2. The toxicity and pharmacokinetics of the compounds were analyzed. The ADME (Absorption, Distribution, Metabolism, and Excretion) properties of the compounds were evaluated by means of mathematical algorithm, and the bioavailability of compounds could be predicted by the most classic Lipinski rule to further obtain the potential drug-forming compounds with activity.

Vardhan [60]bases on molecular docking, ADME and drug-likeness prediction to screen potential inhibitors for ACE2 which is the therapeutic protein targets of SARS-CoV-2 from limonoids and triterpenoids. The results show that glycyrrhizic acid (16), ursolic acid (18), maslinic acid (20), and obacunone (21) were found binding at the catalytic site of ACE2 or the RBD site of ACE2. Natesh [61] uses molecular docking, ADME and Lipinski rules to analyze the binding efficiency of culinary spice bioactives to SARS-CoV-2 target protein. Among the 46 culinary spice bioactives, farnesiferol B (65), sesamin (61), sesaminol (62), and sesamolin (63) have good affinity with SARS-CoV-2 S protein and ACE2. Xu [62] screens the natural products acting on ACE2 from Traditional Chinese Medicine Systems Pharmacology Database (TCMSP). According to the binding affinity of natural products with ACE2, ADME and Lipinski rules, and the results showed that quercetin (1), luteolin (8), and isorhamnetin (2) have good binding ability with ACE2. Finally, the research focuses on the food rich in three active ingredients. Through large-scale literature search, it is concluded that the food rich in three active ingredients may have potential anti-COVID-19 bioactive are red wine, Chinese hawthorn, and blackberry.

Molecular docking is committed to identifying compounds with higher affinity for protein active sites from library or database [63]. However, for specific single traditional Chinese medicine (TCM) and compound preparations, it is usually necessary to use network pharmacology to find the material basis and molecular mechanism [64]. Therefore, the network pharmacology integrated the structure, pharmacological action and action pathway of TCM, and then combined with molecular docking technology to further clarify the material basis and action mechanism of TCM and compound preparation. Ephedra sinica Stapf-Prunus armeniaca L. is a common couplet medicine in anti-COVID-19 prescriptions. Gao [65] screed 47 potentially active ingredients from Ephedra sinica Stapf-Prunus armeniaca L. on the basis of the existing literature. They screened potential action targets of Ephedra sinica Stapf-Prunus armeniaca L. against COVID-19 from TCMSP database through network pharmacology, constructed herb-component-target (H–C-T) network and protein-protein interaction (PPI) network. Finally, they used ACE2 as the target to simulate the docking of 19 main compounds in Ephedra sinica Stapf-Prunus armeniaca L. with ACE2. The result showed that luteolin (8), β-sitosterol (60) had high binding activity to ACE2. Li [66] screened the active components and potential targets of Maxing Shigan Decoction through network pharmacology, and revealed its potential material basis and mechanism of action in antiviral and anti-inflammatory aspects. Through molecular docking, they further identified the components with potential binding ability to ACE2. Among them, amygdalin (57) of Prunus armeniaca L. had the strongest binding ability to ACE2.

In addition, using machine learning method in molecular docking can get more reasonable screening model from known data. In recent years, deep learning method based on neural network has shown great potential in the field of VS. SSnet is a deep neural network framework based on end-to-end, which uses protein structure and ligand information to predict protein ligand interaction probability [67]. SSnet is superior to Atomnet, 3D-CNN, Autodock Vina and other methods in identifying protein ligand complexes with high binding affinity [68]. At present, most of the studies focus on drug screening for S protein-ACE2 binding. This method only considers a single ACE2 structure and cannot identify potential allosteric inhibitors formed by ACE2-S protein [37]. Karki [68]first conducted drug screening for ACE2 in the open conformation and ACE2:S1 complex from the approved drug molecular library through SSnet, and verified the screened drugs through smina virtual docking. The results showed that natural products such as avonoids, flavonones, and polyketides have a strong binding ability with ACE2.

As an important method to assist new drug research and development, VS technology based on molecular docking plays a very important role in the early virtual screening of drug design and the search for new targets. Through the docking software, the screening time is effectively reduced and the screening efficiency is greatly improved. However, there are some common problems in molecular docking technology. First of all, the spatial conformation of the receptor structure is complex and changeable, so it may not be possible to obtain its complete and real three-dimensional structure. Secondly, during molecular docking, the molecular model of the receptor protein is almost impossible to be completely rigid and can undergo certain deformation. Therefore, whether the flexibility and solvation effect of receptor protein molecules in biological system can be truly simulated may directly affect the results of molecular docking.

Molecular dynamics simulations (MD)

Proteins are dynamic biomolecules, and their flexibility will affect the ligand recognition process. Molecular docking has limited ability to deal with the inherent flexibility of protein during docking, which greatly limits the expected applicability of the technology [69]. To solve this problem, several strategies have been tried, including flexible butting and local optimization of side chains, but these relatively simple approaches do not allow for extension and rearrangement of protein structures and significant conformational changes at the skeleton level [[70], [71], [72]]. Therefore, the flexibility of receptors remains a major challenge in ligand recognition. MD is a method based on the physical model of interatomic interaction to predict how each atom in a protein or other molecule changes with time. Based on virtual screening, the flexibility of protein is considered to obtain more reasonable ligand binding sites. MD can not only provide abundant dynamic structure information of biological macromolecules, but also provide a lot of energy information of protein ligand interaction, which provides a theoretical basis for understanding the nature of protein ligand interaction and guiding the process of drug design.

The recessive binding sites of protein ligand binding are invisible, which can only be seen in the crystal structure during the binding process. These sites can also be used as effective drug targets, but due to their concealment, it is difficult to find them through experimental screening. Mixing MD is one of the best methods to identify and characterize recessive binding sites [73]. It is a cosolvent simulation technique for identification of binding hots pots such as orthosteric, allosteric and cofactor binding sites as well as PPI sites [74]. Gopinath [75] used MixMD to detect potential inhibitor binding sites on S protein-ACE2 PPI interface. Drug-like organic probe molecules were added to the solvent to observe their localization in the simulation process, so as to detect the possible small molecule binding sites on the surface of S protein-ACE2 interface. Then, the natural products that bind to S protein and interfere with SARS-CoV-2 binding with host cells were screened through the database, and the compounds with strong and stable binding to RBD were further verified by MD.

MD can also be combined with molecular docking, which can be used to evaluate the structural stability of protein ligand complexes, and improve the accuracy of prediction of drug target binding ability, so as to obtain compounds with high activity [76]. Therefore, molecular docking and MD are highly complementary computational methods for drug screening. Mondal [77] selected 60 active components from 22 kinds of medicinal plants for molecular docking with S protein and ACE2 receptor proteins. The results showed that withanolide A (59) had a high affinity towards two targets of SARS-CoV-2 with the binding energies of −7.7 kcal/mol for S protein and −7.0 kcal/mol for ACE2, theoretically, it has stronger binding ability to S protein than ACE2. And the interaction between withanolide A (59) and ACE2 was almost stable after MD verification, indicating that withanolide A (59) could block the ACE2 receptor and make it difficult for SARS-CoV-2 to enter host cells.

Compared with molecular docking, MD takes full account of the flexibility of protein and can obtain a variety of target protein conformations. Combining with molecular docking is also a strategy to improve the enrichment factor of VS. In the study of S protein-ACE2 blockers, molecular docking is often used in combination with MD. MD can also simulate experimental conditions such as temperature, pressure, atomic number, ion concentration and solvent type used for solvating molecules [78]. All factors can be adjusted and controlled by the statistical mechanics integration under the MD, which is more conducive to drug screening and design.

At present, VS has become the main means for researchers to discover the lead compounds of enzyme inhibitors. VS greatly reduces the number of experimental screening compounds, shortens the research cycle, reduces the cost of drug research and development, and does not need the physical existence of compounds, so computational screening can be performed prior to synthesis. In addition, with the continuous improvement of virtual database, the continuous innovation of screening software, and the combination of the latest computer algorithm technology, the function of VS is more perfect and the results are more accurate. On the other hand, the flexibility of protein structure is a major limitation of VS. The binding sites of proteins are flexible, and there are many conformational states, which often involve the structural changes of rings and the secondary structural elements reconstruction caused by ligand binding [79,80], so false positive results will appear. Due to the complexity of receptor binding interaction, it is difficult to accurately predict the correct binding position of the complex, resulting in the best receptor-binding performance with some software, but poor performance with others [54]. Therefore, VS is usually combined with biological activity experiment.

Surface plasmon resonance (SPR)

SPR is a kind of biosensor technology used to study the interaction between various biomolecules. As a new biochemical analysis and detection tool, SPR is widely used in drug screening with the characteristics of label free, high sensitivity and low consumption [81]. In the classical SPR screening method, the target protein is coupled to the chip surface as the receptor, and then the analyte molecular solution flows through the chip surface. The compounds in the analyte solution that can bind with the receptor are used as the ligands. The binding between the ligand and the receptor leads to the increase of the surface mass of the sensor and the change of the refractive index. Through the computer control system, it is transformed into the sensing signal graph of time and resonance signal. By analyzing the sensing signal graph, the binding affinity and kinetic parameters of ligand and protein can be obtained, the affinity of compounds can be evaluated, and the compounds with good affinity can be selected [82]. The process of SPR is shown in Fig. 3 . The detection principle of SPR is based on the change of the refractive index of the chip surface, and it does not need to label the small molecules to be tested. Therefore, it can eliminate the influence of modified labeling on the activity of compounds, and the detection results can more truly reflect the affinity between compounds and proteins [83,84].

Fig. 3

Schematic diagram of SPR. The target protein ACE2 or S protein are coupled to the chip surface as the receptor, and then, the analyte molecular solution flows through the chip surface. The compounds in the analyte solution that can bind with the. The binding between the ligand and the receptor leads to the increase in the surface mass of the sensor. The computer control system is transformed into the sensing signal graph of time and the resonance signal. By analyzing the sensing signal graph, the binding affinity and kinetic parameters of the ligand and protein can be obtained. The affinity of the compounds was evaluated and the compounds with good affinity were screened.

SPR technology can be used to rapidly screen lead compounds by detecting the KD value of the interaction between the target protein and small molecules. Zhu [85]captured the ACE2-His protein as ligand, dissolved S-RBD-mFc in buffer solution as analyte to establish an SPR screening model of S-RBD and ACE2. When the compounds flowed through the chip surface fixed with ACE2 or S-RBD, the KD value was determined to screen the compounds that could bind to ACE2 or S-RBD. Among them, demethylzeylasteral (19) can bind to S-RBD-mFc in mobile buffer solution with the lowest dissociation rate constant and the longest residence time. Demethylzeylaster (18) was placed in the analytical solution containing RBD and continuously flowed over the chip surface of ACE2 to observe its signal response. The results showed that demethylzeylaster (18) could also block the binding of S-RBD to ACE2.

SPR also can be used as in vitro molecular authentication model combined with VS. For complex system screening of active compounds, some non-specific components can be eliminated through VS. Then SPR can monitor the whole process of the selected compounds binding to the target protein in real time to further verify the binding ability of the compounds obtained from VS. Ye [86]used network pharmacology, molecular docking and SPR technology to explore the potential compounds and interaction mechanism of Toujie Quwen granules in the treatment of COVID-19. Firstly, the H–C-T network was constructed by network pharmacology, and the key compounds were screened to dock with ACE2. The selected quercetin (1), astragaloside IV(24), rutin (4) were tested by SPR. S Protein and ACE2 are fixed in CM5 chip by covalent bond, and the screened compound solution flows through the chip respectively to test the binding energy with S Protein and ACE2. The results showed that quercetin (1) and astragaloside IV (24) bound to SARS-CoV-2 S protein and inactivated it, while rutin (4) bound to ACE2.

Compared with the traditional inhibitor screening based on enzyme activity, SPR technology does not depend on the enzyme reaction process and is more efficient. It does not need to label the small molecules to be detected, and it can monitor the whole process of the binding of compounds with protein targets in real time, so as to analyze the binding of compounds with protein and generate sensing images to evaluate the affinity and kinetic information of candidate compounds. In addition, SPR has strong specificity, and it does not require high matrix of samples. It can also detect colored solutions such as traditional Chinese medicine extracts, which further improves the screening efficiency.

Biochromatography

Biochromatography is a new chromatography technology developed by the combination of life science and chromatography in the 1980s. It has been widely used in the screening of active components of drugs and the study of action mechanism [87]. It uses bioactive materials (such as active macromolecules, active cell membranes, living cells, etc.) fixed on chromatographic carriers as stationary phases to simulate the interactions between drugs and biological macromolecules, targets or cells, or to simulate some key steps of biological activity expression of drugs in vivo under physiological and pathological conditions, avoiding the tedious and chemical separation process of the system [88]. The retention of a drug in a biochromatographic column is directly related to its activity, and it can exhibit certain pharmacological effects when it binds to biological macromolecules, targets or cells [89,90]. According to the different stationary phases, biochromatography can be divided into: Molecular-Biology Chromatography (MBC), Cell Membrane Chromatography (CMC), biomembrane chromatography, or artificial biomembrane chromatography [91]. Among them, MBC and CMC are widely used in the study of S protein-ACE2 blockers. The schematic diagram of MBC and CMC is shown in Fig. 4 .

Fig. 4

Schematic diagram of Molecular-Biology Chromatography(MBC) and Cell Membrane Chromatography(CMC). MBC uses enzymes, receptors, DNA, transporters in plasma, and other biomacromolecules with essential physiological functions as stationary phases to separate and purify active compounds. CMC adsorbs the cell membrane containing specific receptors on the surface of SiO2. The self-fusion characteristics of living cell phospholipids and the strong binding characteristics of silicon hydroxyl (Si–OH) on SiO2 form the cell membrane stationary phase, which can identify the specific targets in complex samples that produce the biological effects of ligands.

Molecular-Biology Chromatography (MBC)

MBC is a new technology which can screen, separate and identify the active components in complex system [92]. When drugs are distributed in the body, they can interact specifically with biological macromolecules such as antibodies, enzymes and receptors. Based on the specific interaction of biomacromolecules, MBC uses enzymes, receptors, DNA, transporters in plasma and other biomacromolecules with important physiological functions as stationary phases to separate and purify active compounds and determine their biochemical parameters [93,94]. The ability to identify the chemical structure of bioactive compounds can be enhanced by establishing an integrated system with MBC as the core and combining with NMR, MS and other methods which can provide structural information.

Chen [95] first used UPLC-HRMS to analyze the prototype and metabolite components in plasma and urine after taking Lianhuaqingwen capsule (LHQW), covalently fixed ACE2 protein on silica gel stationary phase, and then used comprehensive 2D ACE2 column/C18 column/TOFMS system to screen potential ACE2 binding components in enriched urine extract. After preliminary screening, 8 active ingredients that may bind ACE2 were obtained. Through SPR, it was further confirmed that the active ingredients screened by biological chromatography comprehensive system had good binding force with ACE2, thus revealing the high efficiency of screening target compounds through this system.

The composition of natural products is complex, and the conventional chromatographic separation has impurity interference. MBC is simple, time-consuming, with high selectivity and good repeatability, which can eliminate most of the impurity interference. In addition, MBC can be directly related to some pharmacological parameters such as plasma protein binding rate or strength, which has certain pharmacological significance [88]. Therefore, MBC is widely used to screen, isolate and identify active components. However, there are still some problems in the research of MBC in drug active ingredient screening. As an in vitro analysis method of active components, MBC cannot fully represent the process of drug active components combined with biological macromolecules in vivo, and the preparation process of biological chromatographic column is relatively complex and its life is short, which hinders its application [91].

Cell membrane chromatography (CMC)

CMC is a biochromatographic method based on the selective interaction between membrane receptors and their ligands. The cell membrane containing specific receptors is adsorbed on the surface of SiO2. The self-fusion characteristics of living cell phospholipids and the strong binding characteristics of silicon hydroxyl (Si–OH) on SiO2 form the cell membrane stationary phase [96]. The binding characteristics between ligands and membrane receptors can effectively identify the specific target components in complex samples that produce ligand cell biological effects [88]. With the rapid development of molecular and cell biology and Bioengineering, CMC can screen active ingredients from complex samples such as natural plants, and use the specific affinity between drugs and membrane receptors to transform the drug interaction process in vivo into chromatographic process, so as to conduct biomimetic research on drug interaction in vitro [97]. Compared with the traditional methods, CMC has the dual functions of chromatography separation and activity characterization, and has a unique advantage in high-throughput screening of target components in complex samples. By combining CMC with HPLC and MS, the specificity, sensitivity and selectivity of target components identified from complex systems are improved, and the screening and identification efficiency of active components are greatly improved [98].

Lv [99] collected ACE2 h cells, separated cell membrane, mixed with pre activated SiO2 to make cell membrane stationary phase, and then loaded it into the standard column to prepare ACE2/CMC column. An ACE2/CMC bioaffinity chromatography model combined with HPLC-IT-TOF-MS system was established to screen three active components——Ephedrine (42), pseudoephedrine (43), and methylephedrine (44) from Ephedra sinica Stapf, which can act on ACE2. Gao [100] used ACE2-HEK293T cell membrane as stationary phase to prepare CMC column to screen potential ACE2 binding inhibitors in Scutellaria baicalensis Georgi. The result showed that oroxylin A (5) had the strongest binding ability to ACE2. At the same time, molecular docking and SPR verified that oroxylin A (5) bonded well to ACE2 from different perspectives. Further studies on autophagy and apoptosis of oroxylin A (5) showed that it could effectively reduce autophagy of ACE2 cells, and had no effect on apoptosis induced by ACE2 cells.

CMC, as a new receptor pharmacology research method, is simple and reliable, with the characteristics of specificity, applicability, saturation, multi-target, etc. It can combine the separation and screening of active ingredients, avoid the complex process of studying drug metabolism reaction in vivo, drug resistance difference between species and individuals, and eliminate the interference of impurity components, which can be used for specific, selective screening of active components in complex systems [101]. CMC has been widely used since its inception, but there are still some problems in the process of experiment and application. First of all, the membrane receptors attached to the silica gel of the cell membrane chromatography column are easy to fall off or lose activity, so the service life is short, and it is difficult to promote and apply in a wide range. There are many kinds of active protein receptors on the membrane, which may be interfered by non-specific binding in the process of CMC research. Therefore, blank carrier column should be used as the control group to eliminate the interference of non-specific binding on the phospholipid bimolecular membrane. In addition, as an in vitro analysis method of active ingredients, CMC cannot completely simulate the influence of complex environment and other systems in vivo on drug efficacy [98,102]. Therefore, it is necessary to carry out sufficient pharmacological evaluation experiments for further verification.

Biotin-avidin with enzyme-linked immunosorbent assay

Enzyme linked immunosorbent assay (ELISA) is a qualitative and quantitative method for immune reaction by adsorbing soluble antigen or antibody onto solid carrier such as polystyrene. It combines the specificity of antigen-antibody binding with the high efficiency of enzyme reaction. It can detect very low content of antigen such as protein, peptide, hormone or antibody in liquid samples, making the detection signal amplified. It has the advantages of high efficiency, high sensitivity and strong specificity, and has been widely used in various fields of life science [103]. According to the principle, it can be divided into many types, and competitive ELISA is the most widely used in the high-throughput screening of small molecule inhibitors. In order to improve the detection sensitivity, the conventional ELISA detection method has been improved, and a new method combining ELISA with other technologies has been developed. Biotin-Avidin with ELISA is one of the most widely used methods.

Biotin-Avidin System is a new type of amplification system of biological reaction, which can be combined with ELISA organically. Through the high affinity and multi-site binding between biotin and avidin, the majority of biological reaction signals can be amplified, and the specificity and sensitivity can be further improved [104]. In addition, the specific binding ability of avidin and biotin is 1 million times of antigen-antibody, and is not affected by the concentration of reaction reagent. The schematic diagram of biotin-avidin with ELISA is shown in Fig. 5 . Once the binding is extremely stable. Biotinylated ELISA has been proved to be able to evaluate the binding ability of SARS-CoV S proteins with ACE2. Ho [105] covalently linked biotin to SARS-CoV S protein for biotinylation and established competitive ELISA. The biotinylated S protein was mixed with different concentrations of drug extracts, and incubated at 37 °C for 2 h. Then the mixture was added to the adsorption plate coated with ACE2 for OD value determination. The results showed that emodin (53), the active component of Rheum palmatum L., blocked the binding of SARS-CoV S protein to ACE2 in a dose-dependent manner. Some scholars [106] selected 5 drugs from 958 FDA approved drug databases that could destroy the interaction between SARS-CoV-2 RBD and ACE2. Biotinylated RBD and 10-fold diluted drugs were added to 96 well plates coated with ACE2 respectively, and the RBD bound to ACE2 was detected by avidin HRP and TMB substrate. The results of ELISA-based high-throughput showed that RBD/ACE2 interaction could be inhibited by the 5 drugs at selected low and high concentrations. At present, the reports about screening drugs with anti-SARS-CoV-2 activity are all based on authentic virus, which limits the unqualified laboratories to carry out drug screening. The high-throughput ELISA based methodology can be better applied to the screening of anti-SARS-CoV-2 activity drugs.

Fig. 5

Schematic diagram of biotin-avidin with ELISA. Biotin was added into S protein, and after the S protein was incubated with biotin, the biotinylated S protein was added to the pores of 96 plate coated with ACE2, and avidin was added to each hole. The binding ability of S protein and ACE2 was evaluated by ELISA. For competitive ELISA, biotinylated S protein was mixed with drug extract and incubated. After incubation, the mixture was added into the pores of 96 plate coated with ACE2. After adding avidin, the absorbance was measured by ELISA.

In addition, ELISA can also be used to quantitatively evaluate the inhibitory activity of compounds in vitro. Gangadevi [107] found that kobophenol A (35) and its three metabolites have strong binding ability in the ACE2/spike interface and the hydrophobic pocket of the ACE2 domain through molecular docking. Therefore, kobophenol A (35) is the primary target of ACE2/spike RBD in vitro. The recombinant 2019-nCoV S1-RBD was added ACE2 receptor protein and a certain concentration of kobophenol A (35) and verified by ELISA in vitro. It was found that kobophenol A (35) inhibited the binding of ACE2 to SARS-CoV-2 S1-RBD, indicating that kobophenol A (35) may inhibit the entry of S SARS-CoV-2 into host cells.

ELISA screening method has high sample handling capacity, high sensitivity and low false positive rate for detecting ligand-receptor interaction complexes on cell surface, and is especially suitable for high-throughput screening [108]. However, due to its high selectivity to reagents, ELISA is difficult to analyze multiple components at the same time, and has a certain degree of cross-reaction to compounds with similar structures. In the screening process, antibodies need to be incubated and washed for many times, which leads to high screening cost. Despite these limitations, ELISA has the potential to rapidly explore ligand-receptor interactions.

Gene Chip Technology

Gene chip Technology, also known as DNA microarray, is a kind of biological Chip [109]. It refers to the technology that high-density DNA fragments are attached to the surface of silicon, glass and other solid phases in a certain order or arrangement by microarray technology, DNA probes are labeled with isotope or fluorescence, and a large number of gene expression and monitoring are carried out by virtue of the principle of base complementary hybridization [110,111]. Among them, the expression profile chip technology is one of the most widely used technology, it can not only be used to analyze gene function and explore the pathogenesis of diseases, but also be used in high-throughput screening of drugs. Through the analysis of the parallel expression profile chip, the normal human cells or abnormal cells compared with pathological changes to detect differentially expressed genes in different specimens, so as to identify drug targets. It greatly improves the traditional experimental method of observing the expression of single or multiple genes, so as to speed up the identification of differentially expressed genes and the construction of differential expression profile [112]. It can also directly analyze the changes of gene expression profiles of different tissues and organs before and after medication, construct gene expression profiles, and screen new drugs or lead compounds efficiently by analyzing gene expression profiles, pathology and physiology [113,114].

Glinsky [115]first identified human genes that regulate the expression and function of ACE2 and FURIN genes through Gene Set Enrichment Analyses (GSEA) to establish a genomic regulatory interaction model that may influence SARS-CoV-2 infection. The expression profiling experiments (EPEs) revealed that ACE2 and FURIN were ubiquitous in human tissues. Then the transcription factor binding sites affecting the expression of ACE2 and FURIN were identified by Gene Expression Omnibus (GEO). The results showed that HMGA2, INSIG1, RUNX1, HNF-4α, JNK1/c-FOS activated the expression of ACE2 and FURIN. The GSEA of the drug perturbations from GEO database focused on downregulated genes identified estradiol and quercetin (1) among the top significantly enriched records, in which quercetin (1) targets the expression of ACE2, which may be an effective inhibitor of SARS-CoV-2 infection.

Compared with traditional biotechnology, gene chip can be used to analyze thousands of genes in the same test process, which can be used to track the changes of thousands of genes in cells after drug action [116]. It has the characteristics of high throughput and high sensitivity, which changes the limitation of traditional methods that can only study one or several genes. However, at present, this technology has high cost and lack of repeatability, which directly affects the evaluation results of differential genes and reduces the credibility of differential gene screening. In addition, gene screening experiments are mostly limited to the cell level, so we should pay more attention to the animal level and clinical level experiments.

Others

“Spectrum-effect relationship analysis” is a potential method to determine the effective components in complex systems. On the basis of modern chromatographic separation and analysis technology, we can obtain the most useful chemical information, evaluate the pharmacodynamic activity of the corresponding components, and finally connect the chemical information of each component with the pharmacodynamic results, so as to determine the compounds related to the pharmacodynamics [117]. At present, the spectrum-effect relationship has been introduced into the screening of enzyme inhibitors [118]. Through the determination method of enzyme catalyzed substrate reaction combined with the separation and identification of components with enzyme inhibitory activity, the target natural products can be quickly determined, so as to screen the lead compounds more effectively. Senthil Kumar [119]confirmed that geranium and lemon oils could significantly inhibit the expression of ACE2 in epithelial cells by Western blotting and qPCR analysis, and identified the main components of the two essential oils by GC-MS analysis. Finally, the inhibitory effect of each component on ACE2 was determined, citronellol (29) and limonene (30) in geranium and lemon oils were screened out to have strong inhibitory effect on ACE2. Amplified luminescent proximity homogeneous assay (AlphaLISA) is a kind of new homogeneous detection technology based on photoluminescence of nano-beads. It mainly uses the principle of intermolecular interaction, such as antigen-antibody specific interaction, and takes silica microspheres as the carrier. When the interaction between molecules makes the microspheres close to each other, the interaction between molecules can be detected by exciting fluorescence [120]. Hanson [121]has developed a method based on AlphaLISA to detect the binding ability of SARS-CoV-2 RBD to ACE2. A simplified Spike RBD-ACE2 interaction model was established by using recombinant SARS-CoV-2 RBD fused with Fc fragment and recombinant human ACE2 fused with biotinylated fragment. Using this assay platform and TruHits counterassay, compounds with inhibitory activity against ACE2-RBD were screened from 3384 small molecule drugs.

Natural products—S protein-ACE2 blockers

Flavonoids

Flavonoids are natural products found in many plants. Modern pharmacological studies have shown that they have antiviral, anti-inflammatory, antioxidant, immunomodulatory and antitumor effects. As a potential antiviral agent, flavonoids can be traced back to the 1950s. In 1990, the antiviral effect of flavonoids on coronavirus was first found [122]. At present, a large number of drug screening studies have shown that flavonoids can block the combination of SARS-CoV-2 S protein and ACE2, thus inhibiting SARS-CoV-2 from entering the host cells. The chemical structure is shown in Fig. 6 . Many VS methods, such as molecular docketing, MD and net-work pharmacology, screened out quercetin (1), isorhamnetin (2), purearin (3), rutin (4), oroxylin A (5), procyanidin (6), irisulidone (7), luteolin (8), and hesperidin (9) as S-protein-ACE2 blockers [62,86,100,[123], [124], [125]]. Anamika [125] selected hesperidin (9), emodin (53), rhein (54) and chrysin (10) to dock with SARS-CoV2 S protein and its human receptor ACE2 molecule. The results showed that when S protein did not bind to ACE2, these natural products could form stable complexes with ACE2. Hesperidin (9) also binds ACE2 noncompetitively with S protein and the binding sites of S protein and hesperidin (9) are located at different sites in ACE2. Due to the existence of hesperidin (9), the binding structure of ACE2 and S protein becomes unstable. It is further confirmed that hesperidin (9) has antiviral activity as a noncompetitive regulator of S protein binding to its host receptor. The dihydroflavone moiety of hesperidin was predicted to be parallel to the β-6 RBD S protein sheet, while the sugar moiety fits into a shallow hole in the direction away from ACE2 [126]. Some of the natural products of flavonoids were tested for their activity in vitro after VS. Quercetin (1) and astragaloside IV(24) have been proved to bind and inactivate SARS-CoV-2 S protein by SPR in vitro, and then prevent S protein from binding to ACE2, while rutin (4) binds to ACE2, thus weakening the ability of SARS-CoV-2 to infect the host to a certain extent [86]. The antiviral studies in vitro have confirmed that isorhamnetin (2) and oxylin A (5) inhibit SARS-CoV-2 spike pseudotyped virus from entering ACE2h cells [100]. Epigallocatechin-3-gallate (EGCG,15) the most abundant component of catechins in tea (Camellia sinensis (L.) O. Kuntze), plays a role against viruses through inhibiting virus invasiveness, restraining gene expression and replication [127]. Molecular docking experiments revealed that EGCG had a higher atomic contact energy value, binding energy, ligand efficiency and surface area than hydroxy-chloroquine (HCQ) during binding with the S protein [128]. In vitro antiviral experiment results showed that EGCG inhibited SARS-CoV-2 virus replication and interfered directly with SARS-CoV-2 receptor binding [129]. Therefore, as a small compound with broad antiviral activity, EGCG could potentially be used as a lead structure to further develop highly effective antiviral drugs. In addition, there are also a variety of screening methods have been used to screen flavonoids of S protein-ACE2 blockers.

Fig. 6

Structures of flavonoids for S protein-ACE2 repressors. (1) Quercetin; (2) Isorhamnetin; (3) Puerarin; (4) Rutin; (5) Oroxylin A; (6) Procyanidin; (7) Irisolidone; (8) Luteolin; (9) Hesperidin; (10) Chrysin; (11) Anhydrosafflor yellow B; (12) Glabridin; (13) Myricitrin; (14) Euchrenone; (15) Epigallocatechin-3- gallate.

Glinsky [115] found that quercetin (1) targeted ACE2 expression through GSEA and EPEs, which could inhibit the expression of c-FOS in human and rat cells, the expression of RUNX1 in rat cells and HNF-4α in human cells (RUNX1, HNF-4α and c-FOS all can activate the expression of ACE2). After quercetin (1) administration during the differentiation of human intestinal cells, the gene expression of ACE2 was significantly decreased. Furthermore, GSEA analysis of 332 human genes encoding SARS-CoV-2 showed that quercetin (1) changed the expression of 98 genes, which may interfere with 23 functions of SARS-CoV-2 protein in human cells. Therefore, quercetin (1) may be an effective inhibitor of SARS-CoV-2 infection. Zhan [130] found that quercetin (1) and isorhamnetin (2) in Hippophae rhamnoides L. had strong binding effect on ACE2-overexpressed HEK293 (ACE2h) cells through CMC analysis. It was confirmed by SPR that isorhamnetin (2) had a good binding effect with ACE2 and further verification in vitro showed that isorhamnetin (2) inhibited the entry of SARS-CoV-2 spike pseudotyped virus into ACE2 h cells. Therefore, isorhamnetin (2) may be a potential therapeutic candidate for COVID-19.

Flavonoids contain a flavan core with a 15-carbon skeleton. There are two benzene rings, connected by a heterocyclic pyran ring. Studies have shown that heterocyclic pyran ring comprises a C–C double bond and carbonyl groups that play an important role in the biological activities. The hydroxyl groups of two benzene rings are known to be responsible for the radical scavenging activity of flavonoids. Quantitative structure–activity relationship modeling was conducted, and the lack of the heterocyclic pyran ring in the flavonoid skeleton was shown to reduce the inhibitory activity of ACE by up to 91%. The absence of carbonyl groups in the heterocyclic pyran ring also reduced the inhibitory activity of ACE by 74% [131]. The resorcinol moiety might play a role in ACE2 inhibition, as this group might disrupt hydrogen bonds between Glu329/Gln325 of ACE2 and Arg426 of the S protein of SARS-CoV-2, which form a salt bridge to stabilize their interaction [132]. The functional groups of flavonoids, such as the pyran moiety in the heterocyclic pyran ring and hydroxyl groups of the two benzene rings, may play an important role in their ACE2 inhibition. Therefore, the application of flavonoid-based scaffolds in the design of new S protein-ACE2 blockers could be a good approach.

Among flavonoids, quercetin (1) is the most widely studied S protein-ACE2 blocker which has potential antiviral effects. In a randomized, double-blind, placebo-controlled trial, 1002 subjects took 500 or 1000 mg quercetin (1) or placebo daily for 12 days, in which 1000 mg/day quercetin (1) significantly reduced the symptoms of upper respiratory tract infection in middle-aged and elderly people [133]. A clinical study conducted in a hospital in Wuhan showed that 54 patients with COVID-19 were treated with TCM including high content of quercetin (1), which improved the symptoms of patients with COVID-19 and shortened the length of hospital stay of patients [134]. Therefore, quercetin (1) is expected to be used as a new anti- SARS-CoV-2 drug in clinical treatment. However, quercetin (1) has low bioavailability owing to biotransformation during digestion, absorption and metabolism, plants store quercetin (1) attached to sugars which substantially interfere with protein interactions but are released in the human digestive tract [135]. Therefore, other routes of delivery could be more effective in disrupting SARS-CoV-2 S protein/ACE2 interactions. For example, nasal spray containing quercetin (1) in a suitable form could deliver the appropriate concentration directly in the active molecular form, free unconjugated quercetin (1) [136].

Terpenoids

Terpenoids is a volatile substance distributed in essential oils from plant. It is mainly composed of monoterpenes, sesquiterpenes, diterpenes and triterpenes. It has been widely concerned about many pharmacological activities such as anti-tumor, antiviral, antibacterial, anti-inflammatory and immunomodulatory [137]. In recent years, many triterpenoids and their derivatives have been proved to have significant broad-spectrum antiviral activities. Cinatl [138] reported that glycyrrhizic acid (16) and its derivatives have anti-SARS-CoV activity, which can not only inhibit the virus replication process, but also inhibit the virus adsorption and penetration at the beginning of the virus life cycle. Several VS methods such as molecular docking, MD and Net-work pharmacology screened glycyrrhizic acid (16) as S protein-ACE2 blocker [60,66,139]. Yu [140] took S protein RBD of SARS-CoV-2 as the target, verified by SPR that the binding of glycyrrhizic acid (16) with S1 subunit of SARS-CoV-2 destroyed the binding between RBD and ACE2, and evaluated the cytotoxicity of glycyrrhizic acid (16) by MTT. The results showed that glycyrrhizic acid (16) was an effective anti-coronavirus active ingredient in vitro with low toxicity, which might be a potential candidate drug for anti-coronavirus. Zhu [85] established the SPR screening model of S-RBD and ACE2, the results showed that demethylzeylasteral (19) could bind to S-RBD and block the binding of S-RBD and ACE2. In addition, there are other triterpenoids such as ursolic acid (18), Maslinic Acid (20), obacunone (21), atractylenolide III (23), and astragaloside IV (24), which were confirmed that they could bind to ACE2 or S protein as S-protein-ACE2 blockers through a variety of VS methods, and the chemical structure is shown in Fig. 7 . Some diterpenoids such as citronellol (29) and limonene (30) were also proved to have ACE2 inhibitory activity. Senthil Kumar [119]found that citronellol (29) and limonene (30) in geranium and lemon oils significantly down regulated the expression of ACE2 in epithelial cells by spectrum-effect relationship. The cytotoxicity of citronellol (29) and limonene (30) was determined by MTT assay. The results showed that the HT-29 cells treated with two diterpenoids showed no cytotoxicity within 48 h at the concentration of 100 μM. Therefore, geranium and lemon oils and their derivatives citronellol (29) and limonene (30) may have potential antiviral effects.

Fig. 7

Structures of Terpenoids for S protein-ACE2 repressors. (16) Glycyrrhizic acid; (17) Uncaric acid; (18) Ursolic acid; (19) Demethylzeylasteral; (20) Maslinic acid; (21) Obacunone; (22) Andrographolide; (23) Atractylenolide III; (24) Astragaloside IV; (25) Taraxerol; (26) Daturaolone; (27) Limonin; (28) Cucurbitacin G 2-glucoside; (29) Citronellol; (30) Limonene.

Terpenoids play an important role in antiviral drugs, but they have a bottleneck as candidate drugs, such as low solubility, poor intestinal permeability and oral bioavailability. New methods should be introduced to improve their biopharmaceutical properties, especially drug delivery properties [141]. Currently, several delivery systems have been developed successfully, such as nanoemulsions, mesoporous silica nanoparticles, solid lipid nanoparticles, and liposomes, niosomal gels, and solid dispersions. Continuous progress will change the pharmacokinetics of terpenoids and increase their solubility, bioavailability, and even oral bioavailability. Glycyrrhizic acid (16) has an antiviral effect against coronavirus in vitro with low toxicity. It is considered as one potential candidate drug for anti-coronavirus to fight off SARS-CoV-2 by testing alone or combining with other drugs [142]. However, its cytotoxicity and poor solubility in water and biological fluids limit its further clinical application, as it may result in low bioavailability [143]. Therefore, to protect against COVID-19, highly biocompatible glycyrrhizic acid (16) nanoparticles (GANPs) based on the active component of glycyrrhizic acid (16) were synthesized. The results showed that GANPs had significant antiviral, anti-inflammatory, and antioxidant effects in vitro and in vivo. GANPs could target areas of severe inflammation through the EPR effect in the surrogate mouse model of COVID-19, which appeared to improve the accumulation of GANPs in the lungs and livers, further increasing the effectiveness of the treatment. Such a novel approach may provide an effective therapeutic solution for the pandemic, as well as treatment of hyperinflammation of other diseases [144].

Phenols

Phenols are the most abundant secondary metabolites in plants, which are widely used in food and medicine fields because of their antioxidant, anti-tumor, liver protection and other functions [145]. Salvianolic acid A (SAA,31), Salvianolic acid B(SAB,32), Salvianolic acid C(SAC,33) are phenolic carboxylic acids obtained from Salvia miltiorrhiza Bge. Hu [146] confirmed that the toxicity of the three compounds in ACE2h cells was low at a given concentration, and evaluated the binding ability of the three compounds with RBD and ACE2 by molecular docking. They further confirmed that SAA, SAB and SAC had good binding ability with RBD and ACE2 by SPR, and the binding ability of SAB was stronger than SAA and SAC. Finally, in vitro antiviral experiments shown that SAA, SAB and SAC can effectively inhibit the entry of 2019-nCoV pseudovirus into cells, and the inhibitory effect of SAB is greater than SAA and SAC, indicating that SAB may inhibit the binding of 2019-nCoV pseudovirus to ACE2. Chen [95] used the comprehensive 2D ACE2 column/C18 column/TOFMS system to screen the active components in LHQW that may bind to ACE2, and further confirmed that the selected active components have good binding force with ACE2 through SPR. In addition, the ACE2 Inhibitor Screening Kit was used to test the inhibitory activity of active components to ACE2 in vitro, and molecular docking technology was used to simulate the binding sites of the active components. The results showed that neochlorogenic acid (34), one of the phenolic acids, was highly exposed to human body and inhibited ACE2 in vivo. Gangadevi [107] found that kobophenol A (35) had a strong binding ability with the ACE2/spike interface and the hydrophobic pocket of the ACE2 domain through molecular docking, and it was verified by ELISA that kobophenol A (35) inhibited the binding of ACE2 to SARS-CoV-2 S1-RBD. The anti-SARS-CoV2 assay of kobophenol A (35) in Veroe6-EGFP cells further suggests that kobophenol A (35) may inhibit virus entry into the host to cause infection without cytotoxicity. In addition, resveratrol (36), chlorogenic acid (37), gallic acid (38), bis-demethoxycurcumin (39), δ-viniferin (40), pterostilbene (41) and other phenols that can be bind to ACE2 and SARS-CoV-2 spike were screened by a variety of VS methods. And the chemical structure is shown in Fig. 8 .

Fig. 8

Structures of Phenols for S protein-ACE2 repressors. (31) Salvianolic acid A; (32) Salvianolic acid B; (33) Salvianolic acid C; (34) Neochlorogenic acid; (35) Kobophenol A; (36) Resveratrol; (37) Chlorogenic acid; (38) Gallic acid; (39) Bis-demethoxycur-cumin; (40) δ-viniferin; (41) Pterostilbene.

Recent studies have shown that the antiviral activity of these compounds is affected by the number of hydroxyl groups, substituents and the relative position of functional groups. Compounds with more hydroxyl groups showed higher antiviral activity. The location of hydroxyl groups in the benzene ring at the meta position (resorcinol structure) in the flavan ring which may disrupt hydrogen bonds between amino acid residues between ACE2 and S protein [132]. In addition, it appears that compounds with a methyl group have a stronger antiviral activity than compounds with a methoxy group, the position of methoxy group on benzene ring also affects antiviral activity [147]. Therefore, optimizing the structure of phenols compounds may be a new strategy for designing antiviral drugs. In addition, metal ions or metal-containing compounds show antiviral activity, and the presence of metal ions significantly affects the affinity of compounds to SARS-CoV-2 [148]. Complexation of phenols compounds with metals favorably changes their lipophilicity, water-solubility, bioavailability, stability and also allows the introduction of important macro-and microelements into the compound [149]. Therefore, metal complexation of phenols compounds may be one of the drug design strategies for SARS-CoV-2 to enhance the antiviral properties of phenols compounds.

Alkaloids

Alkaloids are a class of nitrogen-containing organic compounds with various structures, which have good antiviral activity and are an important source of antiviral drugs. Lv [99] established ACE2/CMC-HPLC-IT-TOF-MS system and screened three active components EP (42), PEP (43) and MEP (44) in Ephedra sinica Stapf to act on ACE2. Molecular docking and SPR were used to evaluate the binding ability of the three active compounds to ACE2 and SARS-CoV-2 RBD, and CCK-8 staining was used to evaluate the cytotoxicity of the selected active compounds to ACE2h cells. The results showed that the three active ingredients had double binding effect on ACE2 and SARS-CoV-2 RBD, and had little cytotoxicity to ACE2h cells at a certain concentration. In vitro simulation of SARS-CoV-2 spike pseudovirus entry into ACE2h cells showed that the selected compounds, especially EP (41), could inhibit the entry of SARS-CoV-2 spike pseudovirus into ACE2h cells. Gyebi [150] screened the alkaloids from African plants by molecular docking and MD simulations to block the receptor function of ACE2 and destroy the fusion of SARS-CoV-2 and host cells. The results showed that cryptospirolepine (45) and cryptoquindoline (46) have high affinity with ACE2, TMPRSS2 and SARS-CoV-2 S protein. It can bind to the complex formed by S protein RBD and ACE2 (RBD-ACE2), even when S1 subunit of S protein binds to ACE2, cryptospirolepine (45) and cryptoquindoline (46) may destroy the cell entry of SARS-CoV-2 and the stability of RBD-ACE2 complex. In addition, there are many alkaloids, such as speciophylline (47), uncarine F (48), solanine (49), bicuculline (50), anisodamine (51), which can bind to ACE2 and S-protein. The chemical structure is shown in Fig. 9 .

Fig. 9

Structures of Alkaloids for S protein-ACE2 repressors. (42) Ephedrine; (43) Pseudoephedrine; (44) Methylephedrine; (45) Cryptospirolepine; (46) Cryptoquindoline; (47) Speciophylline; (48) Uncarine F; (49) Solanine; (50) Bicuculline; (51) Anisodamine.

Bis-benzylisoquinoline alkaloids can effectively inhibit S protein-ACE2 mediated cell membrane fusion by targeting host calcium channels, and the fusion rate is reduced by about 90%.These compounds block host calcium channels, thereby inhibiting Ca2+-mediated fusion and inhibiting virus entry. Although less research has been done on screening S protein-ACE2 blockers for bis-benzylisoquinoline alkaloids, these compounds also have potential as lead agents for the development of SARS-CoV-2 therapies, especially given the effectiveness of calcium channel blockers in controlling hypertension. Bis-benzylisoquinoline alkaloids should be considered for the treatment of SARS-COV-2 infection in hypertensive patients [151]. Due to the structural diversity and the screening methods of the alkaloids, it is not easy to establish a relationship between chemical structure and anti-SARS-CoV-2 activity of the alkaloids. Interestingly, the alkaloids are highly oxygenated due to the presence of several functional groups, including methoxyls, hydroxyls, methylene dioxide, and carbonyls of ketones, esters, and amides. Also, the nitrogens are heterocyclic atoms and are sometimes in ring fusion [152]. These may be the mechanism and influencing factors of alkaloids to anti-SARS-CoV-2.

In addition to the natural products mentioned above, there are also many compounds that can block the binding of S-protein-ACE2, including quinones, glycosides, steroids and lignans, with the structural formula as shown in Fig. 10 . For example, the components of LHQW that may bind to ACE2 and have ACE2 inhibitory activity in vitro were screened by comprehensive 2D ACE2 column/C18 column/TOFMS system, including rhein (54), forsythoside A (55) and forsythoside I (56) [95]. Interestingly, the compounds in cooking spices including sesamin (61), sesaminol (62), sesamolin (63), and farnesifero B (65) have been shown to bind ACE2 and S proteins through molecular docking [61]. Because of its antioxidant, anti-inflammatory, immunomodulatory and antiviral effects, some scholars suggest that health products based on spices can be developed to prevent and treat COVID-19.

Fig. 10

Structures of others S protein-ACE2 repressors. Quinones: (52) Thymoquinone; (53) Emodin; (54) Rhein; Glycosides: (55) Forsythoside A; (56) Forsythoside I; (57) Amygdalin; Steroids: (58) Urso-deoxycholic acid; (59) Withanolide A; (60) β-sitosterol; Lignans: (61) Sesamin; (62) Sesaminol; (63) Sesamolin; (64) Phillyrin; (65) Farnesiferol B.

Conclusion and perspectives

In recent years, with the combination of pharmaceutical analysis, molecular biology, bioinformatics and artificial intelligence, new vitality has been injected into the development of drug screening strategies and technologies, which greatly improves the efficiency of drug research and development from natural products. At present, S protein-ACE2 blocker screening methods mainly include molecular docking, MD and other VS methods. Except VS, other screening methods such as SPR, MBC, CMC and biotin avidin with ELISA are based on biochemical experiments. VS effectively reduces the screening time and improves the screening efficiency. However, due to the complex and changeable conformation of receptor structure, the binding sites of receptor protein are flexible and there are many conformational states, so false positive results will appear. Therefore, in order to improve the accuracy of VS results, it is usually combined with biological activity experiments. In addition to the above screening methods that have been applied to S protein-ACE2 blockers, a variety of screening methods used to detect protein-protein interactions can also be considered for the screening of S protein-ACE2 blockers. Affinity mass spectrometry, for example, is an effective method to efficiently find small molecular ligands with strong affinity to target proteins. It can efficiently find ligands with strong binding to target proteins from complex mixtures [153]. Fluorescence polarization assay is mainly based on the principle of positive correlation between the intensity of polarized fluorescence and the size of fluorescent molecules, which is widely used in high-throughput drug screening and drug analysis [154]. It should be noted that each of the above screening techniques has its own scope of application, and a single technique often has certain limitations in accuracy or efficiency. Therefore, in practical research, it is necessary to carry out multi-technology combination to screen and verify the accuracy of the results. For example, molecular docking combined with MD, on the one hand, can accurately locate the binding sites and recessive binding sites in the dynamic process of drug binding to the target, so as to more accurately explore the potential active compounds. On the other hand, after molecular docking is used to determine the active compounds with good binding ability to the target, MD can be used to further verify the stability of the active compounds and target complexes, and the results of molecular docking can be further verified. For complex systems, SPR can be combined with VS, and some non-specific components can be eliminated by VS. Then the whole process of the selected compounds binding to the target protein can be monitored by SPR in real time to further verify the binding ability of the compounds obtained by VS to the target protein.

As one of the important sources of new drug research and development, natural products have attracted the attention of researchers at home and abroad for a long time. The use of natural products to prevent or treat COVID-19 is largely inspired by the treatment of SARS. Many studies have proved the role of natural products in the treatment or prevention of viral diseases. A variety of natural products, such as flavonoids, terpenoids, phenols, alkaloids, were obtained by modern screening technology, which can be used as S protein-ACE2 blockers. Although a number of studies have proved that natural products have a significant effect on patients with COVID-19, the above natural products are lack of in vitro and clinical verification. Therefore, it is still necessary to increase the number of samples for randomized clinical trials to evaluate the efficacy and safety of natural products on patients with COVID-19.

Based on the history of drug development, the combination of natural products and chemical synthesis can produce more effective drugs. As for the active lead compounds screened for a variety of S Protein-ACE2 blockers, some of them may not have good ADMET parameters, so we can further modify their structures to obtain good ADMET parameters to enhance their biological activity, so as to develop as candidate drugs. Phenols and flavonoids are considered to have broad potential for preventing viral infections and modulating immune responses in nutritional and functional foods. Based on in vitro analysis and molecular docking models, it will be possible to design and chemically synthesize more efficient compounds based on phenols and flavonoid structures. Both natural and synthetic compounds require functional validation and are not limited to in vitro assays based on recombinant viral proteins. In addition, for the key problem of low bioavailability of candidate compounds, they can be administered through different routes from oral (e.g., inhalation, intravenous), which protects the drug molecule from intestinal metabolism and adsorption and improves its pharmacokinetic and pharmacodynamic parameters. Most of the current research work is based on theoretical calculations and in vitro bioassays, and there is still a lot of work to be done in the real sense of drug design and development. Nevertheless, these studies also laid the foundation for follow-up drug discovery and provided guidance for further in-depth in vivo research.

Author contributions

Jun-zhi Lin and Run-chun Xu put forward the idea; Le-le Ma wrote the manuscript and made the figures; Hui-min Liu, Xue-mei Liu, Xiao-yu Yuan, Chao Xu, and Fang Wang collected the literatures; Ding-kun Zhang contributed to the revisions.

Compliance with ethics guidelines

The authors declare that they have no conflict of interest.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.