Elucidation of Binding Features and Dissociation Pathways of Inhibitors and Modulators in SARS-CoV-2 Main Protease by Multiple Molecular Dynamics Simulations

COVID-19 can cause different neurological symptoms in some people, including smell, inability to taste, dizziness, confusion, delirium, seizures, stroke, etc. Owing to the issue of vaccine effectiveness, update and coverage, we still need one or more diversified strategies as the backstop to manage illness. Characterizing the structural basis of ligand recognition in the main protease (M^pro) of SARS-CoV-2 will facilitate its rational design and development of potential drug candidates with high affinity and selectivity against COVID-19. Up to date, covalent-, non-covalent inhibitors and allosteric modulators have been reported to bind to different active sites of M^pro. In the present work, we applied the molecular dynamics (MD) simulations to systematically characterize the potential binding features of catalytic active site and allosteric binding sites in M^pro using a dataset of 163 3D structures of M^pro-inhibitor complexes, in which our results are consistent with the current studies. In addition, umbrella sampling (US) simulations were used to explore the dissociation processes of substrate pathway and allosteric pathway. All the information provided new insights into the protein features of M^pro and will facilitate its rational drug design for COVID-19.     

Conducting the RBD of SARS-CoV-2 Omicron Variant with Phytoconstituents from Euphorbia dendroides to Repudiate the Binding of Spike Glycoprotein Using Computational Molecular Search and Simulation Approach

(1) Background: Natural constituents are still a preferred route for counteracting the outbreak of COVID-19. Essentially, flavonoids have been found to be among the most promising molecules identified as coronavirus inhibitors. Recently, a new SARS-CoV-2 B.1.1.529 variant has spread in many countries, which has raised awareness of the role of natural constituents in attempts to contribute to therapeutic protocols. (2) Methods: Using various chromatographic techniques, triterpenes (1–7), phenolics (8–11), and flavonoids (12–17) were isolated from Euphorbia dendroides and computationally screened against the receptor-binding domain (RBD) of the SARS-CoV-2 Omicron variant. As a first step, molecular docking calculations were performed for all investigated compounds. Promising compounds were subjected to molecular dynamics simulations (MD) for 200 ns, in addition to molecular mechanics Poisson–Boltzmann surface area calculations (MM/PBSA) to determine binding energy. (3) Results: MM/PBSA binding energy calculations showed that compound 14 (quercetin-3-O-β-D-glucuronopyranoside) and compound 15 (quercetin-3-O-glucuronide 6″-O-methyl ester) exhibited strong inhibition of Omicron, with ΔG_binding of −41.0 and −32.4 kcal/mol, respectively. Finally, drug likeness evaluations based on Lipinski’s rule of five also showed that the discovered compounds exhibited good oral bioavailability. (4) Conclusions: It is foreseeable that these results provide a novel intellectual contribution in light of the decreasing prevalence of SARS-CoV-2 B.1.1.529 and could be a good addition to the therapeutic protocol.     

Insights into the Dynamics and Binding of Two Polyprotein Substrate Cleavage Points in the Context of the SARS-CoV-2 Main and Papain-like Proteases

It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1–PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro–CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.     

The Discovery of Small Allosteric and Active Site Inhibitors of the SARS-CoV-2 Main Protease via Structure-Based Virtual Screening and Biological Evaluation

The main protease enzyme (M^pro) of SARS-CoV-2 is one of the most promising targets for COVID-19 treatment. Accordingly, in this work, a structure-based virtual screening of 3.8 million ligand libraries was carried out. After rigorous filtering, docking, and post screening assessments, 78 compounds were selected for biological evaluation, 3 of which showed promising inhibition of the M^pro enzyme. The obtained hits (CB03, GR04, and GR20) had reasonable potencies with K_i values in the medium to high micromolar range. Interestingly, while our most potent hit, GR20, was suggested to act via a reversible covalent mechanism, GR04 was confirmed as a noncompetitive inhibitor that seems to be one of a kind when compared to the other allosteric inhibitors discovered so far. Moreover, all three compounds have small sizes (~300 Da) with interesting fittings in their relevant binding sites, and they possess lead-like characteristics that can introduce them as very attractive candidates for the future development of COVID-19 treatments.     

The Development of Pharmacophore Models for the Search of New Natural Inhibitors of SARS-CoV-2 Spike RBD–ACE2 Binding Interface

To date, some succeeding variants of SARS-CoV-2 have become more contagious. This virus is known to enter human cells by binding the receptor-binding domain (RBD) of spike protein with the angiotensin-converting enzyme 2 (ACE2), the latter being a membrane protein that regulates the renin–angiotensin system. Since the host cell receptor plays a critical role in viral entry, inhibition of the RBD–ACE2 complex is a promising strategy for preventing COVID-19 infection. In the present communication, we propose and utilize an approach based on the generation of a complex of pharmacophore models and subsequent Induced Fit Docking (IFD) to identify potential inhibitors of the main binding sites of the Omicron SARS-CoV-2 RBD(S1)–ACE2 complex (PDB ID: 7T9L) among a number of natural products of various types and origins. Several natural compounds have been found to provide a high affinity for the receptor of interest. It is expected that the present results will stimulate further research aimed at the development of specialized drugs against this virus.     

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

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

Derivation,Functionalization of (S)-Goniothalamin from Goniothalamus wightii and Their Derivative Targets SARS-CoV-2 MPro,SPro, and RdRp: A Pharmacological Perspective

The tracing of an alternative drug, Phytochemicals is a promising approach to the viral threats that have emerged over the past two years. Across the world, herbal medicine is a better solution against anti-viral diseases during pandemic periods. Goniothalamus wightii is an herbal plant, which has diverse bioactive compounds with anticancer, antioxidant, and anti-viral properties. The aim of the study was to isolate the compound by chromatography studies and functionalization by FT-IR, LC-MS, and NMR (C-NMR, H-NMR). As a result, the current work focuses on whether (S)-Goniathalamin and its analogue could act as natural anti-viral molecules for multiple target proteins viz., M^Pro, RdRp, and S^Pro, which are required for SARS-CoV-2 infection. Overall, 954 compounds were examined and the molecular-docking studies were performed on the maestro platform of Schrodinger software. Molecular-dynamics simulation studies were performed on two complex major compounds to confirm their affinity across 150 simulations. This research suggests that plant-based drugs have high levels of antiviral properties against coronavirus. However, more research is needed to verify its antiviral properties.     

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

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

Molecular and Structural Analysis of Specific Mutations from Saudi Isolates of SARS-CoV-2 RNA-Dependent RNA Polymerase and their Implications on Protein Structure and Drug–Protein Binding

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has stressed the global health system to a significant level, which has not only resulted in high morbidity and mortality but also poses a threat for future pandemics. This situation warrants efforts to develop novel therapeutics to manage SARS-CoV-2 in specific and other emerging viruses in general. This study focuses on SARS-CoV2 RNA-dependent RNA polymerase (RdRp) mutations collected from Saudi Arabia and their impact on protein structure and function. The Saudi SARS-CoV-2 RdRp sequences were compared with the reference Wuhan, China RdRp using a variety of computational and biophysics-based approaches. The results revealed that three mutations—A97V, P323I and Y606C—may affect protein stability, and hence the relationship of protein structure to function. The apo wild RdRp is more dynamically stable with compact secondary structure elements compared to the mutants. Further, the wild type showed stable conformational dynamics and interaction network to remdesivir. The net binding energy of wild-type RdRp with remdesivir is -50.76 kcal/mol, which is more stable than the mutants. The findings of the current study might deliver useful information regarding therapeutic development against the mutant RdRp, which may further furnish our understanding of SARS-CoV-2 biology.     

An Effective MM/GBSA Protocol for Absolute Binding Free Energy Calculations: A Case Study on SARS-CoV-2 Spike Protein and the Human ACE2 Receptor

The binding free energy calculation of protein–ligand complexes is necessary for research into virus–host interactions and the relevant applications in drug discovery. However, many current computational methods of such calculations are either inefficient or inaccurate in practice. Utilizing implicit solvent models in the molecular mechanics generalized Born surface area (MM/GBSA) framework allows for efficient calculations without significant loss of accuracy. Here, GBNSR6, a new flavor of the generalized Born model, is employed in the MM/GBSA framework for measuring the binding affinity between SARS-CoV-2 spike protein and the human ACE2 receptor. A computational protocol is developed based on the widely studied Ras–Raf complex, which has similar binding free energy to SARS-CoV-2/ACE2. Two options for representing the dielectric boundary of the complexes are evaluated: one based on the standard Bondi radii and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein–ligand binding. Predictions based on the two radii sets provide upper and lower bounds on the experimental references: −14.7(ΔGbindBondi)<−10.6(ΔGbindExp.)<−4.1(ΔGbindOPT1) kcal/mol. The consensus estimates of the two bounds show quantitative agreement with the experiment values. This work also presents a novel truncation method and computational strategies for efficient entropy calculations with normal mode analysis. Interestingly, it is observed that a significant decrease in the number of snapshots does not affect the accuracy of entropy calculation, while it does lower computation time appreciably. The proposed MM/GBSA protocol can be used to study the binding mechanism of new variants of SARS-CoV-2, as well as other relevant structures.     

Evidence for Multiple Binding Modes in the Initial Contact Between SARS-CoV-2 Spike S1 Protein and Cell Surface Glycans**

Infection of host cells by SARS-CoV-2 begins with recognition by the virus S (spike) protein of cell surface heparan sulfate (HS), tethering the virus to the extracellular matrix environment, and causing the subunit S1-RBD to undergo a conformational change into the ‘open’ conformation. These two events promote the binding of S1-RBD to the angiotensin converting enzyme 2 (ACE2) receptor, a preliminary step toward viral-cell membrane fusion. Combining ligand-based NMR spectroscopy with molecular dynamics, oligosaccharide analogues were used to explore the interactions between S1-RBD of SARS CoV-2 and HS, revealing several low-specificity binding modes and previously unidentified potential sites for the binding of extended HS polysaccharide chains. The evidence for multiple binding modes also suggest that highly specific inhibitors will not be optimal against protein S but, rather, diverse HS-based structures, characterized by high affinity and including multi-valent compounds, may be required.     

D936Y and Other Mutations in the Fusion Core of the SARS-CoV-2 Spike Protein Heptad Repeat 1: Frequency,Geographical Distribution,and Structural Effect

The crown of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is constituted by its spike (S) glycoprotein. S protein mediates the SARS-CoV-2 entry into the host cells. The “fusion core” of the heptad repeat 1 (HR1) on S plays a crucial role in the virus infectivity, as it is part of a key membrane fusion architecture. While SARS-CoV-2 was becoming a global threat, scientists have been accumulating data on the virus at an impressive pace, both in terms of genomic sequences and of three-dimensional structures. On 15 February 2021, from the SARS-CoV-2 genomic sequences in the GISAID resource, we collected 415,673 complete S protein sequences and identified all the mutations occurring in the HR1 fusion core. This is a 21-residue segment, which, in the post-fusion conformation of the protein, gives many strong interactions with the heptad repeat 2, bringing viral and cellular membranes in proximity for fusion. We investigated the frequency and structural effect of novel mutations accumulated over time in such a crucial region for the virus infectivity. Three mutations were quite frequent, occurring in over 0.1% of the total sequences. These were S929T, D936Y, and S949F, all in the N-terminal half of the HR1 fusion core segment and particularly spread in Europe and USA. The most frequent of them, D936Y, was present in 17% of sequences from Finland and 12% of sequences from Sweden. In the post-fusion conformation of the unmutated S protein, D936 is involved in an inter-monomer salt bridge with R1185. We investigated the effect of the D936Y mutation on the pre-fusion and post-fusion state of the protein by using molecular dynamics, showing how it especially affects the latter one.     

Ultra-Large-Scale Screening of Natural Compounds and Free Energy Calculations Revealed Potential Inhibitors for the Receptor-Binding Domain (RBD) of SARS-CoV-2

The emergence of immune-evading variants of SARS-CoV-2 further aggravated the ongoing pandemic. Despite the deployments of various vaccines, the acquired mutations are capable of escaping both natural and vaccine-induced immune responses. Therefore, further investigation is needed to design a decisive pharmacological treatment that could efficiently block the entry of this virus into cells. Hence, the current study used structure-based methods to target the RBD of the recombinant variant (Deltacron) of SARS-CoV-2, which was used as a model variant. From the virtual drug screenings of various databases, a total of four hits were identified as potential lead molecules. Key residues were blocked by these molecules with favorable structural dynamic features. The binding free energies further validated the potentials of these molecules. The TBE for MNP was calculated to be −32.86 ± 0.10 kcal/mol, for SANC00222 the TBE was −23.41 ± 0.15 kcal/mol, for Liriodenine the TBE was −34.29 ± 0.07 kcal/mol, while for Carviolin the TBE was calculated to be −27.67 ± 0.12 kcal/mol. Moreover, each complex demonstrated distinct internal motion and a free energy profile, indicating a different strategy for the interaction with and inhibition of the RBD. In conclusion, the current study demands further in vivo and in vitro validation for the possible usage of these compounds as potential drugs against SARS-CoV-2 and its variants.     

In Vitro and In Silico Studies for the Identification of Potent Metabolites of Some High-Altitude Medicinal Plants from Nepal Inhibiting SARS-CoV-2 Spike Protein

Despite ongoing vaccination programs against COVID-19 around the world, cases of infection are still rising with new variants. This infers that an effective antiviral drug against COVID-19 is crucial along with vaccinations to decrease cases. A potential target of such antivirals could be the membrane components of the causative pathogen, SARS-CoV-2, for instance spike (S) protein. In our research, we have deployed in vitro screening of crude extracts of seven ethnomedicinal plants against the spike receptor-binding domain (S1-RBD) of SARS-CoV-2 using an enzyme-linked immunosorbent assay (ELISA). Following encouraging in vitro results for Tinospora cordifolia, in silico studies were conducted for the 14 reported antiviral secondary metabolites isolated from T. cordifolia—a species widely cultivated and used as an antiviral drug in the Himalayan country of Nepal—using Genetic Optimization for Ligand Docking (GOLD), Molecular Operating Environment (MOE), and BIOVIA Discovery Studio. The molecular docking and binding energy study revealed that cordifolioside-A had a higher binding affinity and was the most effective in binding to the competitive site of the spike protein. Molecular dynamics (MD) simulation studies using GROMACS 5.4.1 further assayed the interaction between the potent compound and binding sites of the spike protein. It revealed that cordifolioside-A demonstrated better binding affinity and stability, and resulted in a conformational change in S1-RBD, hence hindering the activities of the protein. In addition, ADMET analysis of the secondary metabolites from T. cordifolia revealed promising pharmacokinetic properties. Our study thus recommends that certain secondary metabolites of T. cordifolia are possible medicinal candidates against SARS-CoV-2.     

Plant-Derived Antiviral Compounds as Potential Entry Inhibitors against Spike Protein of SARS-CoV-2 Wild-Type and Delta Variant: An Integrative in SilicoApproach

The wild-type SARS-CoV-2 has continuously evolved into several variants with increased transmissibility and virulence. The Delta variant which was initially identified in India created a devastating impact throughout the country during the second wave. While the efficacy of the existing vaccines against the latest SARS-CoV-2 variants remains unclear, extensive research is being carried out to develop potential antiviral drugs through approaches like in silico screening and drug-repurposing. This study aimed to conduct the docking-based virtual screening of 50 potential phytochemical compounds against a Spike glycoprotein of the wild-type and the Delta SARS-CoV-2 variant. Subsequently, molecular docking was performed for the five best compounds, such as Lupeol, Betulin, Hypericin, Corilagin, and Geraniin, along with synthetic controls. From the results obtained, it was evident that Lupeol exhibited a remarkable binding affinity towards the wild-type Spike protein (−8.54 kcal/mol), while Betulin showed significant binding interactions with the mutated Spike protein (−8.83 kcal/mol), respectively. The binding energy values of the selected plant compounds were slightly higher than that of the controls. Key hydrogen bonding and hydrophobic interactions of the resulting complexes were visualized, which explained their greater binding affinity against the target proteins—the Delta S protein of SARS-CoV-2, in particular. The lower RMSD, the RMSF values of the complexes and the ligands, Rg, H-bonds, and the binding free energies of the complexes together revealed the stability of the complexes and significant binding affinities of the ligands towards the target proteins. Our study suggests that Lupeol and Betulin could be considered as potential ligands for SARS-CoV-2 spike antagonists. Further experimental validations might provide new insights for the possible antiviral therapeutic interventions of the identified lead compounds and their analogs against COVID-19 infection.     

Structural Characterization of N-Linked Glycans in the Receptor Binding Domain of the SARS-CoV-2 Spike Protein and their Interactions with Human Lectins

The glycan structures of the receptor binding domain of the SARS-CoV2 spike glycoprotein expressed in human HEK293F cells have been studied by using NMR. The different possible interacting epitopes have been deeply analysed and characterized, providing evidence of the presence of glycan structures not found in previous MS-based analyses. The interaction of the RBD ¹³C-labelled glycans with different human lectins, which are expressed in different organs and tissues that may be affected during the infection process, has also been evaluated by NMR. In particular, ¹⁵N-labelled galectins (galectins-3, -7 and -8 N-terminal), Siglecs (Siglec-8, Siglec-10), and C-type lectins (DC-SIGN, MGL) have been employed. Complementary experiments from the glycoprotein perspective or from the lectin's point of view have permitted to disentangle the specific interacting epitopes in each case. Based on these findings, 3D models of the interacting complexes have been proposed.     

Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development

Severe acute respiratory syndrome coronavirus (SARS-CoV-2) disease is a global rapidly spreading virus showing very high rates of complications and mortality. Till now, there is no effective specific treatment for the disease. Aloe is a rich source of isolated phytoconstituents that have an enormous range of biological activities. Since there are no available experimental techniques to examine these compounds for antiviral activity against SARS-CoV-2, we employed an in silico approach involving molecular docking, dynamics simulation, and binding free energy calculation using SARS-CoV-2 essential proteins as main protease and spike protein to identify lead compounds from Aloe that may help in novel drug discovery. Results retrieved from docking and molecular dynamics simulation suggested a number of promising inhibitors from Aloe. Root mean square deviation (RMSD) and root mean square fluctuation (RMSF) calculations indicated that compounds 132, 134, and 159 were the best scoring compounds against main protease, while compounds 115, 120, and 131 were the best scoring ones against spike glycoprotein. Compounds 120 and 131 were able to achieve significant stability and binding free energies during molecular dynamics simulation. In addition, the highest scoring compounds were investigated for their pharmacokinetic properties and drug-likeness. The Aloe compounds are promising active phytoconstituents for drug development for SARS-CoV-2.     

Structural Dynamics of the SARS-CoV-2 Spike Protein: A 2-Year Retrospective Analysis of SARS-CoV-2 Variants (from Alpha to Omicron) Reveals an Early Divergence between Conserved and Variable Epitopes

We analyzed the epitope evolution of the spike protein in 1,860,489 SARS-CoV-2 genomes. The structural dynamics of these epitopes was determined by molecular modeling approaches. The D614G mutation, selected in the first months of the pandemic, is still present in currently circulating SARS-CoV-2 strains. This mutation facilitates the conformational change leading to the demasking of the ACE2 binding domain. D614G also abrogated the binding of facilitating antibodies to a linear epitope common to SARS-CoV-1 and SARS-CoV-2. The main neutralizing epitope of the N-terminal domain (NTD) of the spike protein showed extensive structural variability in SARS-CoV-2 variants, especially Delta and Omicron. This epitope is located on the flat surface of the NTD, a large electropositive area which binds to electronegatively charged lipid rafts of host cells. A facilitating epitope located on the lower part of the NTD appeared to be highly conserved among most SARS-CoV-2 variants, which may represent a risk of antibody-dependent enhancement (ADE). Overall, this retrospective analysis revealed an early divergence between conserved (facilitating) and variable (neutralizing) epitopes of the spike protein. These data aid in the designing of new antiviral strategies that could help to control COVID-19 infection by mimicking neutralizing antibodies or by blocking facilitating antibodies.     

Fullerenes against COVID-19: Repurposing C60 and C70 to Clog the Active Site of SARS-CoV-2 Protease

The persistency of COVID-19 in the world and the continuous rise of its variants demand new treatments to complement vaccines. Computational chemistry can assist in the identification of moieties able to lead to new drugs to fight the disease. Fullerenes and carbon nanomaterials can interact with proteins and are considered promising antiviral agents. Here, we propose the possibility to repurpose fullerenes to clog the active site of the SARS-CoV-2 protease, M^pro. Through the use of docking, molecular dynamics, and energy decomposition techniques, it is shown that C₆₀ has a substantial binding energy to the main protease of the SARS-CoV-2 virus, M^pro, higher than masitinib, a known inhibitor of the protein. Furthermore, we suggest the use of C₇₀ as an innovative scaffold for the inhibition of SARS-CoV-2 M^pro. At odds with masitinib, both C₆₀ and C₇₀ interact more strongly with SARS-CoV-2 M^pro when different protonation states of the catalytic dyad are considered. The binding of fullerenes to M^pro is due to shape complementarity, i.e., vdW interactions, and is aspecific. As such, it is not sensitive to mutations that can eliminate or invert the charges of the amino acids composing the binding pocket. Fullerenic cages should therefore be more effective against the SARS-CoV-2 virus than the available inhibitors such as masinitib, where the electrostatic term plays a crucial role in the binding.