Structural analysis, virtual screening and molecular simulation to identify potential inhibitors targeting 2'-O-ribose methyltransferase of SARS-CoV-2 coronavirus

Hesperidin, a Potential Antiviral Agent against SARS-CoV-2: The Influence of Citrus Consumption on COVID-19 Incidence and Severity in China

This review examines hesperidin, a citrus bioflavonoid, as a potential antiviral agent against SARS-CoV-2. The COVID-19 pandemic has demanded an urgent need to search for effective antiviral compounds, including those of natural origin, such as hesperidin. The review provides a comprehensive analysis of the chemical properties, bioavailability and antiviral mechanisms of hesperidin, particularly its potential efficacy against SARS-CoV-2. A review of databases, including PubMedPico, Scopus and Web of Science, was conducted using specific keywords and search criteria in accordance with PRISMA (Re-porting Items for Systematic Reviews and Meta-Analysis) guidelines between 2020 and 2024. Of the 207 articles, 37 were selected for the review. A key aspect is the correlation of in vitro, in silico and clinical studies on the antiviral effects of hesperidin with epidemiological data on citrus consumption in China during 2020-2024. The importance of integrating laboratory findings with actual consumption patterns to better understand the role of hesperidin in mitigating COVID-19 was highlighted, and an attempt was made to analyze epidemiological studies to examine the association between citrus juice consumption as a source of hesperidin and the incidence and severity of COVID-19 using China as an example. The review identifies consistencies and discrepancies between experimental and epidemiological data, highlighting the need to correlate the two fields to better understand the potential of hesperidin as an agent against SARS-CoV-2. Challenges and limitations in interpreting the results and future research perspectives in this area are discussed. The aim of this comprehensive review is to bridge the gap between experimental studies and epidemiological evidence and to contribute to the understanding of their correlation.

A review on in silico virtual screening methods in COVID-19 using anticancer drugs and other natural/chemical inhibitors

The present coronavirus disease 2019 (COVID-19) pandemic scenario has posed a difficulty for cancer treatment. Even under ideal conditions, malignancies like small cell lung cancer (SCLC) are challenging to treat because of their fast development and early metastases. The treatment of these patients must not be jeopardized, and they must be protected as much as possible from the continuous spread of the COVID-19 infection. Initially identified in December 2019 in Wuhan, China, the contagious coronavirus illness 2019 (COVID-19) is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Finding inhibitors against the druggable targets of SARS-CoV-2 has been a significant focus of research efforts across the globe. The primary motivation for using molecular modeling tools against SARS-CoV-2 was to identify candidates for use as therapeutic targets from a pharmacological database. In the published study, scientists used a combination of medication repurposing and virtual drug screening methodologies to target many structures of SARS-CoV-2. This virus plays an essential part in the maturation and replication of other viruses. In addition, the total binding free energy and molecular dynamics (MD) modeling findings showed that the dynamics of various medications and substances were stable; some of them have been tested experimentally against SARS-CoV-2. Different virtual screening (VS) methods have been discussed as potential means by which the evaluated medications that show strong binding to the active site might be repurposed for use against SARS-CoV-2.

Coronavirus 2'-O-methyltransferase: A promising therapeutic target

Coronaviruses (CoVs) have been the source of multiple epidemics and a global pandemic since the start of century, and there is an urgent need to understand CoV biology and develop better therapeutics. Here, we review the role of NSP16 in CoV replication, specifically its importance to 2'-O-methylation and CoV RNA capping. We describe the attenuation phenotypes of NSP16-mutant CoVs, the roles of MDA5 and IFITs in sensing and antagonizing viral RNA lacking 2'O methylation, and the dependence on 2'-O-methylation in other virus families. We also detail the growing body of research into targeting 2'-O-methylation for therapeutics or as a platform for live attenuated vaccines. Beyond its role in RNA capping, NSP16 may have yet uncharacterized importance to CoV replication, highlighting the need for continued studies into NSP16 functions. Understanding the full contribution of NSP16 to the replicative fitness of CoVs will better inform the development of treatments against future CoV outbreaks.

Identifying inhibitors of NSP16-NSP10 of SARS-CoV-2 from large databases

The COVID-19 pandemic, which has already claimed millions of lives, continues to pose a serious threat to human health, requiring the development of new effective drugs. Non-structural proteins of SARS-CoV-2 play an important role in viral replication and infection. Among them, NSP16 (non-structured protein 16) and its cofactor NSP10 (non-structured protein 10) perform C2'-O methylation at the 5' end of the viral RNA, which promotes efficient virus replication. Therefore, the NSP16-NSP10 complex becomes an attractive target for drug development. Using a multi-step virtual screening protocol which includes Lipinski's rule, docking, steered molecular dynamics and umbrella sampling, we searched for potential inhibitors from the PubChem and anti-HIV databases. It has been shown that CID 135566620 compound from PubChem is the best candidate with an inhibition constant in the sub-μM range. The Van der Waals interaction was found to be more important than the electrostatic interaction in the binding affinity of this compound to NSP16-NSP10. Further in vitro and in vivo studies are needed to test the activity of the identified compound against COVID-19.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 Uses Nonstructural Protein 16 To Evade Restriction by IFIT1 and IFIT3

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'-O-methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'-O-MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo, using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive than wild-type SARS-CoV-2 to type I interferon (IFN-I) in vitro. Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'-O-methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, an MTase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment and attenuates viral replication. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a target for future antiviral therapies. IMPORTANCE Similar to other coronaviruses, disruption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo, our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1 but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2'-O-methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

High throughput bioluminescent assay to characterize and monitor the activity of SARS-CoV-2 methyltransferases

The fast rate of viral mutations of SARS CoV-2 result in decrease in the efficacy of the vaccines that have been developed before the emergence of these mutations. Thus, it is believed that using additional measures to combat the virus is not only advisable but also beneficial. Two antiviral drugs were authorized for emergency use by the FDA, namely Pfizer's two-drug regimen sold under the brand name Paxlovid, and Merck's drug Lagevrio. Pfizer's two-drug combination consists of nirmatrelvir, a protease inhibitor that blocks coronavirus ability to multiply and another antiviral, ritonavir, that lowers the rate of drug clearance to boost the longevity and activity of the protease inhibitor. Merck's drug Lagevrio (molnupiravir) is a nucleoside analogue with a mechanism of action that aims to introduce errors into the genetic code of the virus. We believe the armament against the virus can be augmented by the addition of another class of enzyme inhibitors that are required for viral survival and its ability to replicate. Enzymes like nsp14 and nsp10/16 methyltransferases (MTases) represent another class of drug targets since they are required for viral RNA translation and evading the host immune system. In this communication, we have successfully verified that the MTase-Glo, which is universal and homogeneous MTase assay can be used to screen for inhibitors of the two pivotal enzymes nsp14 and nsp16 of SARS CoV-2. Furthermore, we have carried out extensive studies on those enzymes using different RNA substrates and tested their activity using various inhibitors and verified the utility of this assay for use in drug screening programs. We anticipate our work will be pursued further to screen for large libraries to discover new and selective inhibitors for the viral enzymes particularly that these enzymes are structurally different from their mammalian counterparts.

SARS-CoV-2 Uses Nonstructural Protein 16 to Evade Restriction by IFIT1 and IFIT3

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'- O methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'- O MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo , using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive to type I interferon (IFN-I) in vitro . Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'- O methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, a methyltransferase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a possible target for future antiviral therapies.

Importance

Similar to other coronaviruses, disruption of SARS-CoV-2 NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo , our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1, but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2'- O methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

African derived phytocompounds may interfere with SARS-CoV-2 RNA capping machinery via inhibition of 2'-O-ribose methyltransferase: An in silico perspective

Despite the ongoing vaccination against the life-threatening COVID-19, there is need for viable therapeutic interventions. The S-adenosyl-l-Methionine (SAM) dependent 2-O'-ribose methyltransferase (2'-O-MTase) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a therapeutic target against COVID-19 infection. In a bid to profile bioactive principles from natural sources, a custom-made library of 226 phytochemicals from African medicinal plants with especially anti-malarial activity was screened for direct interactions with SARS-CoV-2 2'-O-MTase (S2RMT) using molecular docking and molecular dynamics (MD) simulations as well as binding free energies methods. Based on minimal binding energy lower than sinefungin (a reference methyl-transferase inhibitor) and binding mode analysis at the catalytic site of S2RMT, a list of 26 hit phytocompounds was defined. The interaction of these phytocompounds was compared with the 2'-O-MTase of SARS-CoV and MERS-CoV. Among these compounds, the lead phytocompounds (LPs) viz: mulberrofuran F, 24-methylene cycloartenol, ferulate, 3-benzoylhosloppone and 10-hydroxyusambarensine interacted strongly with the conserved KDKE tetrad within the substrate binding pocket of the 2'-O-MTase of the coronavirus strains which is critical for substrate binding. The thermodynamic parameters analyzed from the MD simulation trajectories of the LPs-S2RMT complexes presented an eminent structural stability and compactness. These LPs demonstrated favorable druggability and in silico ADMET properties over a diverse array of molecular computing descriptors. The LPs show promising prospects in the disruption of S2RMT capping machinery in silico. However, these LPs should be validated via in vitro and in vivo experimental models.

Chemical biology and medicinal chemistry of RNA methyltransferases

RNA methyltransferases (MTases) are ubiquitous enzymes whose hitherto low profile in medicinal chemistry, contrasts with the surging interest in RNA methylation, the arguably most important aspect of the new field of epitranscriptomics. As MTases become validated as drug targets in all major fields of biomedicine, the development of small molecule compounds as tools and inhibitors is picking up considerable momentum, in academia as well as in biotech. Here we discuss the development of small molecules for two related aspects of chemical biology. Firstly, derivates of the ubiquitous cofactor S-adenosyl-l-methionine (SAM) are being developed as bioconjugation tools for targeted transfer of functional groups and labels to increasingly visible targets. Secondly, SAM-derived compounds are being investigated for their ability to act as inhibitors of RNA MTases. Drug development is moving from derivatives of cosubstrates towards higher generation compounds that may address allosteric sites in addition to the catalytic centre. Progress in assay development and screening techniques from medicinal chemistry have led to recent breakthroughs, e.g. in addressing human enzymes targeted for their role in cancer. Spurred by the current pandemic, new inhibitors against coronaviral MTases have emerged at a spectacular rate, including a repurposed drug which is now in clinical trial.

Energy and reactivity profile and proton affinity analysis of rimegepant with special reference to its potential activity against SARS-CoV-2 virus proteins using molecular dynamics

Rimegepant is a new medicine developed for the management of chronic headache due to migraine. This manuscript is an attempt to study the various structural, physical, and chemical properties of the molecules. The molecule was optimized using B3LYP functional with 6-311G + (2d,p) basis set. Excited state properties of the compound were studied using CAM-B3LYP functional with same basis sets using IEFPCM model in methanol for the implicit solvent atmosphere. The various electronic descriptors helped to identify the reactivity behavior and stability. The compound is found to possess good nonlinear optical properties in the gas phase. The various intramolecular electronic delocalizations and non-covalent interactions were analyzed and explained. As the compound contain several heterocyclic nitrogen atoms, they have potential proton abstraction features, which was analyzed energetically. The most important result from this study is from the molecular docking analysis which indicates that rimegepant binds irreversibly with three established SARS-CoV-2 proteins with ID 6LU7, 6M03, and 6W63 with docking scores − 9.2988, − 8.3629, and − 9.5421 kcal/mol respectively. Further assessment of docked complexes with molecular dynamics simulations revealed that hydrophobic interactions, water bridges, and π–π interactions play a significant role in stabilizing the ligand within the binding region of respective proteins. MMGBSA-free energies further demonstrated that rimegepant is more stable when complexed with 6LU7 among the selected PDB models. As the pharmacology and pharmacokinetics of this molecule are already established, rimegepant can be considered as an ideal candidate with potential for use in the treatment of COVID patients after clinical studies.

Computational Chemistry to Repurposing Drugs for the Control of COVID-19

Thus far, in 2021, 219 countries with over 175 million people have been infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). SARS-CoV-2 is a positive sense, single-stranded RNA virus, and is the causal agent for coronavirus disease (COVID-19). Due to the urgency of the situation, virtual screening as a computational modeling method offers a fast and effective modality of identifying drugs that may be effective against SARS-CoV-2. There has been an overwhelming abundance of molecular docking against SARS-CoV-2 in the last year. Due to the massive volume of computational studies, this systematic review has been created to evaluate and summarize the findings of existing studies. Herein, we report on computational articles of drugs which target, (1) viral protease, (2) Spike protein-ACE 2 interaction, (3) RNA-dependent RNA polymerase, and (4) other proteins and nonstructural proteins of SARS-CoV-2. Based on the studies presented, there are 55 identified natural or drug compounds with potential anti-viral activity. The next step is to show anti-viral activity in vitro and translation to determine effectiveness into human clinical trials.

Evolutionary selectivity of amino acid is inspired from the enhanced structural stability and flexibility of the folded protein

AIM

The present study attempts to decipher the site-specific amino acid alterations at certain positions experiencing preferential selectivity and their effect on proteins' stability and flexibility. The study examines the selection preferences by considering pair-wise non-bonded interaction energies of adjacent and interacting amino acids present at the interacting site, along with their evolutionary history.

MATERIALS AND METHODS

For the study, variations in the interacting residues of spike protein (S-Protein) receptor-binding domain (RBD) of different coronaviruses were examined. The MD simulation trajectory analysis revealed that, though all the variants studied were structurally stable at their native and bound confirmations, the RBD of 2019-nCoV/SARS-CoV-2 was found to be more flexible and more dynamic. Furthermore, a noticeable change observed in the non-bonded interaction energies of the amino acids interacting with the receptor corroborated their selection at respective positions.

KEY FINDINGS

The conformational changes exerted by the altered amino acids could be the reason for a broader range of interacting receptors among the selected proteins.

SIGNIFICANCE

The results envisage a strong indication that the residue selection at certain positions is governed by a well-orchestrated feedback mechanism, which follows increased stability and flexibility in the folded structure compared to its evolutionary predecessor.

The Role of Coronavirus RNA-Processing Enzymes in Innate Immune Evasion

Viral RNA sensing triggers innate antiviral responses in humans by stimulating signaling pathways that include crucial antiviral genes such as interferon. RNA viruses have evolved strategies to inhibit or escape these mechanisms. Coronaviruses use multiple enzymes to synthesize, modify, and process their genomic RNA and sub-genomic RNAs. These include Nsp15 and Nsp16, whose respective roles in RNA capping and dsRNA degradation play a crucial role in coronavirus escape from immune surveillance. Evolutionary studies on coronaviruses demonstrate that genome expansion in Nidoviruses was promoted by the emergence of Nsp14-ExoN activity and led to the acquisition of Nsp15- and Nsp16-RNA-processing activities. In this review, we discuss the main RNA-sensing mechanisms in humans as well as recent structural, functional, and evolutionary insights into coronavirus Nsp15 and Nsp16 with a view to potential antiviral strategies.

In-silico evaluation of bioactive compounds from tea as potential SARS-CoV-2 nonstructural protein 16 inhibitors

Background and aim

A novel coronavirus, called the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been found to cause COVID-19 in humans and some other mammals. The nonstructural protein 16 (NSP16) of SARS-CoV-2 plays a significant part in the replication of viruses and suppresses the ability of innate immune system to detect the virus. Therefore, inhibiting NSP16 can be a secure path towards identifying a potent medication against SARS-CoV-2. Tea (Camellia sinensis) polyphenols have been reported to exhibit potential treatment options against various viral diseases.

Methods

We conducted molecular docking and structural dynamics studies with a set of 65 Tea bioactive compounds to illustrate their ability to inhibit NSP16 of SARS-CoV-2. Moreover, post-simulations end state thermodynamic free energy calculations were estimated to strengthen our results.

Results and conclusion

Six bioactive tea molecules showed better docking scores than the standard molecule sinefungin. These results were further validated by MD simulations, where Theaflavin compound demonstrated lower binding free energy in comparison to the standard molecule sinefungin. The compound theaflavin could be considered as a novel lead compound for further evaluation by in-vitro and in-vivo studies.

Molecular Simulation-Based Investigation of Highly Potent Natural Products to Abrogate Formation of the nsp10-nsp16 Complex of SARS-CoV-2

The SARS-CoV-2 non-structural protein (nsp) nsp10–nsp16 complex is essential for the 2′-O-methylation of viral mRNA, a crucial step for evading the innate immune system, and it is an essential process in SARS-CoV-2 life cycle. Therefore, detecting molecules that can disrupt the nsp10–nsp16 interaction are prospective antiviral drugs. In this study, we screened the North African Natural Products database (NANPDB) for molecules that can interact with the nsp10 interface and disturb the nsp10–nsp16 complex formation. Following rigorous screening and validation steps, in addition to toxic side effects, drug interactions and a risk /benefit assessment, we identified four compounds (genkwanin-6-C-beta-glucopyranoside, paraliane diterpene, 4,5-di-p-trans-coumaroylquinic acid and citrinamide A) that showed the best binding affinity and most favourable interaction with nsp10 interface residues. To understand the conformational stability and dynamic features of nsp10 bound to the four selected compounds, we subjected each complex to 200 ns molecular dynamics simulations. We then calculated the free binding energies of compounds interacting with nsp10 structure using the molecular mechanics-generalised Born surface area (MMGBSA). Of the four compounds, genkwanin-6-C-beta-glucopyranoside demonstrated the most stable complex with nsp10, in addition to a tighter binding affinity of −37.4 ± 1.3 Kcal/mol. This potential to disrupt the nsp10–nsp16 interface interaction and inhibit it now sets the path for functional studies.

A review on the interaction of nucleoside analogues with SARS-CoV-2 RNA dependent RNA polymerase

The outbreaks of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) in 2019, have highlighted the concerns about the lack of potential vaccines or antivirals approved for inhibition of CoVs infection. SARS-CoV-2 RNA dependent RNA polymerase (RdRp) which is almost preserved across different viral species can be a potential target for development of antiviral drugs, including nucleoside analogues (NA). However, ExoN proofreading activity of CoVs leads to their protection from several NAs. Therefore, potential platforms based on the development of efficient NAs with broad-spectrum efficacy against human CoVs should be explored. This study was then aimed to present an overview on the development of NAs-based drug repurposing for targeting SARS-CoV-2 RdRp by computational analysis. Afterwards, the clinical development of some NAs including Favipiravir, Sofosbuvir, Ribavirin, Tenofovir, and Remdesivir as potential inhibitors of RdRp, were surveyed. Overall, exploring broad-spectrum NAs as promising inhibitors of RdRp may provide useful information about the identification of potential antiviral repurposed drugs against SARS-CoV-2.

Aptamer BC 007 - Efficient binder of spreading-crucial SARS-CoV-2 proteins

Corona virus disease 2019 (COVID-19) is a respiratory disease caused by a new coronavirus (SARS-CoV-2) which causes significant morbidity and mortality. The emergence of this novel and highly pathogenic SARS-CoV-2 and its rapid international spread poses a serious global public health emergency. To date 32,174,627 cases, of which 962,613 (2.99%) have died, have been reported (https://www.who.int/westernpacific/health-topics/coronavirus, accessed 23 Sep 2020). The outbreak was declared a Public Health Emergency of International Concern on 30 January 2020.

There are still not many SARS-CoV-2-specific and effective treatments or vaccines available. A second round of infection is obviously unavoidable.

Aptamers had already been at the centre of interest in the fight against viruses before now. The selection and development of a new aptamer is, however, a time-consuming process. We therefore checked whether a clinically developed aptamer, BC 007, which is currently in phase 2 of clinical testing for a different indication, would also be able to efficiently bind DNA-susceptible peptide structures from SARS-CoV-2-spreading crucial proteins, such as the receptor binding domain (RBD) of the spike protein and the RNA dependent RNA polymerase of SARS-CoV-2 (re-purposing). Indeed, several such sequence-sections have been identified. In particular for two of these sequences, BC 007 showed specific binding in a therapy-relevant concentration range, as shown in Nuclear magnetic resonance (NMR)- and Circular dicroism (CD)-spectroscopy and isothermal titration calorimetry (ITC). The excellent clinical toxicity and tolerability profile of this substance opens up an opportunity for rapid clinical testing of its COVID-19 effectiveness.

COVID-19 and SARS-CoV-2. Modeling the present, looking at the future

Since December 2019 the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has produced an outbreak of pulmonary disease which has soon become a global pandemic, known as COronaVIrus Disease-19 (COVID-19). The new coronavirus shares about 82% of its genome with the one which produced the 2003 outbreak (SARS CoV-1). Both coronaviruses also share the same cellular receptor, which is the angiotensin-converting enzyme 2 (ACE2) one. In spite of these similarities, the new coronavirus has expanded more widely, more faster and more lethally than the previous one. Many researchers across the disciplines have used diverse modeling tools to analyze the impact of this pandemic at global and local scales. This includes a wide range of approaches - deterministic, data-driven, stochastic, agent-based, and their combinations - to forecast the progression of the epidemic as well as the effects of non-pharmaceutical interventions to stop or mitigate its impact on the world population. The physical complexities of modern society need to be captured by these models. This includes the many ways of social contacts - (multiplex) social contact networks, (multilayers) transport systems, metapopulations, etc. - that may act as a framework for the virus propagation. But modeling not only plays a fundamental role in analyzing and forecasting epidemiological variables, but it also plays an important role in helping to find cures for the disease and in preventing contagion by means of new vaccines. The necessity for answering swiftly and effectively the questions: could existing drugs work against SARS CoV-2? and can new vaccines be developed in time? demands the use of physical modeling of proteins, protein-inhibitors interactions, virtual screening of drugs against virus targets, predicting immunogenicity of small peptides, modeling vaccinomics and vaccine design, to mention just a few. Here, we review these three main areas of modeling research against SARS CoV-2 and COVID-19: (1) epidemiology; (2) drug repurposing; and (3) vaccine design. Therefore, we compile the most relevant existing literature about modeling strategies against the virus to help modelers to navigate this fast-growing literature. We also keep an eye on future outbreaks, where the modelers can find the most relevant strategies used in an emergency situation as the current one to help in fighting future pandemics.