Computational analysis of dynamic allostery and control in the SARS-CoV-2 main protease

PyPCN: protein contact networks in PyMOL

MOTIVATION

Protein contact networks (PCNs) represent the 3D structure of a protein using network formalism. Inter-residue contacts are described as binary adjacency matrices, which are derived from the graph representation of residues (as α-carbons, β-carbons or centroids) and Euclidean distances according to defined thresholds. Functional characterization algorithms are computed on binary adjacency matrices to unveil allosteric, dynamic, and interaction mechanisms in proteins. Such strategies are usually applied in a combinatorial manner, although rarely in seamless and user-friendly implementations.

RESULTS

PyPCN is a plugin for PyMOL wrapping more than twenty PCN algorithms and metrics in an easy-to-use graphical user interface, to support PCN analysis. The plugin accepts 3D structures from the Protein Data Bank, user-provided PDBs, or precomputed adjacency matrices. The results are directly mapped to 3D protein structures and organized into interactive diagrams for their visualization. A dedicated graphical user interface combined with PyMOL visual support makes analysis more intuitive and easier, extending the applicability of PCNs.

AVAILABILITY AND IMPLEMENTATION

https://github.com/pcnproject/PyPCN.

Structural insights into SARS-CoV-2 main protease conformational plasticity

The spread of different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants underscores the need for insights into the structural properties of its structural and non-structural proteins. The highly conserved homo-dimeric chymotrypsin-like protease (3CL M^PRO ), belonging to the class of cysteine hydrolases, plays an indispensable role in processing viral polyproteins that are involved in viral replication and transcription. Studies have successfully demonstrated the role of M^PRO as an attractive drug target for designing antiviral treatments because of its importance in the viral life cycle. Herein, we report the structural dynamics of six experimentally solved structures of M^PRO (i.e., 6LU7, 6M03, 6WQF, 6Y2E, 6Y84, and 7BUY including both ligand-free and ligand-bound states) at different resolutions. We have employed a structure-based balanced forcefield, CHARMM36m through state-of-the-art all-atoms molecular dynamics simulations at µ-seconds scale at room temperature (303K) and pH 7.0 to explore their structure-function relationship. The helical domain-III responsible for dimerization mostly contributes to the altered conformational states and destabilization of M^PRO . A keen observation of the high degree of flexibility in the P5 binding pocket adjoining domain II-III highlights the reason for observation of conformational heterogeneity among the structural ensembles of M^PRO . We also observe a differential dynamics of the catalytic pocket residues His41, Cys145, and Asp187, which may lead to catalytic impairment of the monomeric proteases. Among the highly populated conformational states of the six systems, 6LU7 and 7M03 forms the most stable and compact M^PRO conformation with intact catalytic site and structural integrity. Altogether, our findings from this extensive study provides a benchmark to identify physiologically relevant structures of such promising drug targets for structure-based drug design and discovery of potent drug-like compounds having clinical potential.

Computational analysis of the sequence-structure relation in SARS-CoV-2 spike protein using protein contact networks

The structure of proteins impacts directly on the function they perform. Mutations in the primary sequence can provoke structural changes with consequent modification of functional properties. SARS-CoV-2 proteins have been extensively studied during the pandemic. This wide dataset, related to sequence and structure, has enabled joint sequence-structure analysis. In this work, we focus on the SARS-CoV-2 S (Spike) protein and the relations between sequence mutations and structure variations, in order to shed light on the structural changes stemming from the position of mutated amino acid residues in three different SARS-CoV-2 strains. We propose the use of protein contact network (PCN) formalism to: (i) obtain a global metric space and compare various molecular entities, (ii) give a structural explanation of the observed phenotype, and (iii) provide context dependent descriptors of single mutations. PCNs have been used to compare sequence and structure of the Alpha, Delta, and Omicron SARS-CoV-2 variants, and we found that omicron has a unique mutational pattern leading to different structural consequences from mutations of other strains. The non-random distribution of changes in network centrality along the chain has allowed to shed light on the structural (and functional) consequences of mutations.

Potential Self-Peptide Inhibitors of the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in viral replication, cleaving viral polyproteins into functional proteins. This makes Mpro an important drug target. Mpro consists of an N-terminal catalytic domain and a C-terminal α-helical domain (MproC). Previous studies have shown that peptides derived from a given protein sequence (self-peptides) can affect the folding and, in turn, the function of that protein. Since the SARS-CoV-1 MproC is known to stabilize its Mpro and regulate its function, we hypothesized that SARS-CoV-2 MproC-derived self-peptides may modulate the folding and the function of SARS-CoV-2 Mpro. To test this, we studied the folding of MproC in the presence of various self-peptides using coarse-grained structure-based models and molecular dynamics simulations. In these simulations of MproC and one self-peptide, we found that two self-peptides, the α1-helix and the loop between α4 and α5 (loop4), could replace the equivalent native sequences in the MproC structure. Replacement of either sequence in full-length Mpro should, in principle, be able to perturb Mpro function albeit through different mechanisms. Some general principles for the rational design of self-peptide inhibitors emerge: The simulations show that prefolded self-peptides are more likely to replace native sequences than those which do not possess structure. Additionally, the α1-helix self-peptide is kinetically stable and once inserted rarely exchanges with the native α1-helix, while the loop4 self-peptide is easily replaced by the native loop4, making it less useful for modulating function. In summary, a prefolded α1-derived peptide should be able to inhibit SARS-CoV-2 Mpro function.

SARS-CoV-2 proteases Mpro and PLpro: Design of inhibitors with predicted high potency and low mammalian toxicity using artificial neural networks, ligand-protein docking, molecular dynamics simulations, and ADMET calculations

The main (Mpro) and papain-like (PLpro) proteases are highly conserved viral proteins essential for replication of the COVID-19 virus, SARS-COV-2. Therefore, a logical plan for producing new drugs against this pathogen is to discover inhibitors of these enzymes. Accordingly, the goal of the present work was to devise a computational approach to design, characterize, and select compounds predicted to be potent dual inhibitors - effective against both Mpro and PLpro. The first step employed LigDream, an artificial neural network, to create a virtual ligand library. Ligands with computed ADMET profiles indicating drug-like properties and low mammalian toxicity were selected for further study. Initial docking of these ligands into the active sites of Mpro and PLpro was done with GOLD, and the highest-scoring ligands were redocked with AutoDock Vina to determine binding free energies (ΔG). Compounds 89-00, 89-07, 89-32, and 89-38 exhibited favorable ΔG values for Mpro (-7.6 to -8.7 kcal/mol) and PLpro (-9.1 to -9.7 kcal/mol). Global docking of selected compounds with the Mpro dimer identified prospective allosteric inhibitors 89-00, 89-27, and 89-40 (ΔG -8.2 to -8.9 kcal/mol). Molecular dynamics simulations performed on Mpro and PLpro active site complexes with the four top-scoring ligands from Vina demonstrated that the most stable complexes were formed with compounds 89-32 and 89-38. Overall, the present computational strategy generated new compounds with predicted drug-like characteristics, low mammalian toxicity, and high inhibitory potencies against both target proteases to form stable complexes. Further preclinical studies will be required to validate the in silico findings before the lead compounds could be considered for clinical trials.

Dual-Reporter System for Real-Time Monitoring of SARS-CoV-2 Main Protease Activity in Live Cells Enables Identification of an Allosteric Inhibition Path

The SARS-CoV-2 pandemic is an ongoing threat to global health, and the continuing emergence of contagious variants highlights the urgent need for additional antiviral therapy to attenuate COVID-19 disease. The SARS-CoV-2 main protease (3CL^pro) presents an attractive target for such therapy due to its high sequence conservation and key role in the viral life cycle. In this study, we designed a fluorescent-luminescent cell-based reporter for the detection and quantification of 3CL^pro intracellular activity. Employing this platform, we examined the efficiency of known protease inhibitors against 3CL^pro and further identified potent inhibitors through high-throughput chemical screening. Computational analysis confirmed a direct interaction of the lead compounds with the protease catalytic site and identified a prototype for efficient allosteric inhibition. These developments address a pressing need for a convenient sensor and specific targets for both virus detection and rapid discovery of potential inhibitors.

Progress of cGAS-STING signaling in response to SARS-CoV-2 infection

Coronavirus disease 2019 (COVID-19) is an epidemic respiratory disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that can cause infections in millions of individuals, who can develop lung injury, organ failure, and subsequent death. As the first line of host defense, the innate immune system is involved in initiating the immune response to SARS-CoV-2 infection and the hyperinflammatory phenotype of COVID-19. However, the interplay between SARS-CoV-2 and host innate immunity is not yet well understood. It had become known that the cGAS-STING pathway is involved in the detection of cytosolic DNA, which elicits an innate immune response involving a robust type I interferon response against viral and bacterial infections. Nevertheless, several lines of evidence indicate that SARS-CoV-2, a single-stranded positive-sense RNA virus, triggered the cGAS-STING signaling pathway. Therefore, understanding the molecular and cellular details of cGAS-STING signaling upon SARS-CoV-2 infection is of considerable biomedical importance. In this review, we discuss the role of cGAS-STING signaling in SARS-CoV-2 infection and summarize the potential therapeutics of STING agonists as virus vaccine adjuvants.

Promising SARS-CoV-2 main protease inhibitor ligand-binding modes evaluated using LB-PaCS-MD/FMO

Parallel cascade selection molecular dynamics-based ligand binding-path sampling (LB-PaCS-MD) was combined with fragment molecular orbital (FMO) calculations to reveal the ligand path from an aqueous solution to the SARS-CoV-2 main protease (M^pro) active site and to customise a ligand-binding pocket suitable for delivering a potent inhibitor. Rubraxanthone exhibited mixed-inhibition antiviral activity against SARS-CoV-2 M^pro, relatively low cytotoxicity, and high cellular inhibition. However, the atomic inhibition mechanism remains ambiguous. LB-PaCS-MD/FMO is a hybrid ligand-binding evaluation method elucidating how rubraxanthone interacts with SARS-CoV-2 M^pro. In the first step, LB-PaCS-MD, which is regarded as a flexible docking, efficiently samples a set of ligand-binding pathways. After that, a reasonable docking pose of LB-PaCS-MD is evaluated by the FMO calculation to elucidate a set of protein-ligand interactions, enabling one to know the binding affinity of a specified ligand with respect to a target protein. A possible conformation was proposed for rubraxanthone binding to the SARS-CoV-2 M^pro active site, and allosteric inhibition was elucidated by combining blind docking with k-means clustering. The interaction profile, key binding residues, and considerable interaction were elucidated for rubraxanthone binding to both M^pro sites. Integrated LB-PaCS-MD/FMO provided a more reasonable complex structure for ligand binding at the SARS-CoV-2 M^pro active site, which is vital for discovering and designing antiviral drugs.

Allosteric Hotspots in the Main Protease of SARS-CoV-2

Inhibiting the main protease of SARS-CoV-2 is of great interest in tackling the COVID-19 pandemic caused by the virus. Most efforts have been centred on inhibiting the binding site of the enzyme. However, considering allosteric sites, distant from the active or orthosteric site, broadens the search space for drug candidates and confers the advantages of allosteric drug targeting. Here, we report the allosteric communication pathways in the main protease dimer by using two novel fully atomistic graph-theoretical methods: Bond-to-bond propensity, which has been previously successful in identifying allosteric sites in extensive benchmark data sets without a priori knowledge, and Markov transient analysis, which has previously aided in finding novel drug targets in catalytic protein families. Using statistical bootstrapping, we score the highest ranking sites against random sites at similar distances, and we identify four statistically significant putative allosteric sites as good candidates for alternative drug targeting.

PCN-Miner: an open-source extensible tool for the analysis of Protein Contact Networks

MOTIVATION

Protein Contact Network (PCN) is a powerful method for analysing the structure and function of proteins, with a specific focus on disclosing the molecular features of allosteric regulation through the discovery of modular substructures. The importance of PCN analysis has been shown in many contexts, such as the analysis of SARS-CoV-2 Spike protein and its complexes with the Angiotensin Converting Enzyme 2 (ACE2) human receptors. Even if there exist many software tools implementing such methods, there is a growing need for the introduction of tools integrating existing approaches.

RESULTS

We present PCN-Miner, a software tool implemented in the Python programming language, able to (i) import protein structures from the Protein Data Bank; (ii) generate the corresponding PCN; (iii) model, analyse and visualize PCNs and related protein structures by using a set of known algorithms and metrics. The PCN-Miner can cover a large set of applications: from clustering to embedding and subsequent analysis.

AVAILABILITY AND IMPLEMENTATION

The PCN-Miner tool is freely available at the following GitHub repository: https://github.com/hguzzi/ProteinContactNetworks. It is also available in the Python Package Index (PyPI) repository.

Allosteric Binding Sites of the SARS-CoV-2 Main Protease: Potential Targets for Broad-Spectrum Anti-Coronavirus Agents

The current pandemic caused by the COVID-19 disease has reached everywhere in the world and has affected every aspect of our lives. As of the current data, the World Health Organization (WHO) has reported more than 300 million confirmed COVID-19 cases worldwide and more than 5 million deaths. M^pro is an enzyme that plays a key role in the life cycle of the SARS-CoV-2 virus, and it is vital for the disease progression. The M^pro enzyme seems to have several allosteric sites that can hinder the enzyme catalytic activity. Furthermore, some of these allosteric sites are located at or nearby the dimerization interface which is essential for the overall M^pro activity. In this review paper, we investigate the potential of the M^pro allosteric site to act as a drug target, especially since they interestingly appear to be resistant to mutation. The work is illustrated through three subsequent sections: First, the two main categories of M^pro allosteric sites have been explained and discussed. Second, a total of six pockets have been studied and evaluated for their druggability and cavity characteristics. Third, the experimental and computational attempts for the discovery of new allosteric inhibitors have been illustrated and discussed. To sum up, this review paper gives a detailed insight into the feasibility of developing new M^pro inhibitors to act as a potential treatment for the COVID-19 disease.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Optimising Elastic Network Models for Protein Dynamics and Allostery: Spatial and Modal Cut-offs and Backbone Stiffness

The family of coarse-grained models for protein dynamics known as Elastic Network Models (ENMs) require careful choice of parameters to represent well experimental measurements or fully-atomistic simulations. The most basic ENM that represents each protein residue by a node at the position of its C-alpha atom, all connected by springs of equal stiffness, up to a cut-off in distance. Even at this level a choice is required of the optimum cut-off distance and the upper limit of elastic normal modes taken in any sum for physical properties, such as dynamic correlation or allosteric effects on binding. Additionally, backbone-enhanced ENM (BENM) may improve the model by allocating a higher stiffness to springs that connect along the protein backbone. This work reports on the effect of varying these three parameters (distance and mode cutoffs, backbone stiffness) on the dynamical structure of three proteins, Catabolite Activator Protein (CAP), Glutathione S-transferase (GST), and the SARS-CoV-2 Main Protease (M pro ). Our main results are: (1) balancing B-factor and dispersion-relation predictions, a near-universal optimal value of 8.5 Å is advisable for ENMs; (2) inhomogeneity in elasticity brings the first mode containing spatial structure not well-resolved by the ENM typically within the first 20; (3) the BENM only affects modes in the upper third of the distribution, and, additionally to the ENM, is only able to model the dispersion curve better in this vicinity; (4) BENM does not typically affect fluctuation-allostery, which also requires careful treatment of the effector binding to the host protein to capture.

Allosteric Modulation of the Main Protease (MPro) of SARS-CoV-2 by Casticin—Insights from Molecular Dynamics Simulations

Inhibition of the main protease (M^pro) of SARS-CoV-2 has been suggested to be vital in shutting down viral replication in a host. Most efforts aimed at inhibiting M^Pro activity have been channeled into competitive inhibition at the active site, but this strategy will require a high inhibitor concentration and impressive inhibitor-M^Pro binding affinity. Allosteric inhibition can potentially serve as an effective strategy for alleviating these limitations. In this study, the ability of antiviral natural products to inhibit M^Pro in an allosteric fashion was explored with in silico techniques. Molecular docking revealed a strong interaction between casticin, an antiviral flavonoid, and M^pro at a site distant from the active site. This site, characterized as a distal site, has been shown to have an interdependent dynamic effect with the active site region. M^pro only, M^pro-peptide (binary) and M^pro-peptide-casticin (ternary) complexes were subjected to molecular dynamics simulations for 50 ns to investigate the modulatory activity of casticin binding on M^pro. Molecular dynamic simulations revealed that binding of casticin at the distal site interferes with the proper orientation of the peptide substrate in the oxyanion hole of the active site, and this could lead to a halt or decrease in catalytic activity. This study therefore highlights casticin as a potential allosteric modulator of the SARS-CoV-2 main protease, which could be optimized and developed into a potential lead compound for anti-SARS-CoV-2 chemotherapy.

Protocetraric and Salazinic Acids as Potential Inhibitors of SARS-CoV-2 3CL Protease: Biochemical, Cytotoxic, and Computational Characterization of Depsidones as Slow-Binding Inactivators

The study investigated the inhibitory activity of protocetraric and salazinic acids against SARS-CoV-2 3CL^pro. The kinetic parameters were determined by microtiter plate-reading fluorimeter using a fluorogenic substrate. The cytotoxic activity was tested on murine Sertoli TM4 cells. In silico analysis was performed to ascertain the nature of the binding with the 3CL^pro. The compounds are slow-binding inactivators of 3CL^pro with a K_i of 3.95 μM and 3.77 μM for protocetraric and salazinic acid, respectively, and inhibitory efficiency k_inact/K_i at about 3 × 10⁻⁵ s⁻¹µM⁻¹. The mechanism of inhibition shows that both compounds act as competitive inhibitors with the formation of a stable covalent adduct. The viability assay on epithelial cells revealed that none of them shows cytotoxicity up to 80 μM, which is well below the K_i values. By molecular modelling, we predicted that the catalytic Cys145 makes a nucleophilic attack on the carbonyl carbon of the cyclic ester common to both inhibitors, forming a stably acyl-enzyme complex. The computational and kinetic analyses confirm the formation of a stable acyl-enzyme complex with 3CL^pro. The results obtained enrich the knowledge of the already numerous biological activities exhibited by lichen secondary metabolites, paving the way for developing promising scaffolds for the design of cysteine enzyme inhibitors.

Kinematic Vibrational Entropy Assessment and Analysis of SARS CoV-2 Main Protease

The three-dimensional conformations of a protein influence its function and select for the ligands it can interact with. The total free energy change during protein-ligand complex formation includes enthalphic and entropic components, which together report on the binding affinity and conformational states of the complex. However, determining the entropic contribution is computationally burdensome. Here, we apply kinematic flexibility analysis (KFA) to efficiently estimate vibrational frequencies from static protein and protein-ligand structures. The vibrational frequencies, in turn, determine the vibrational entropies of the structures and their complexes. Our estimates of the vibrational entropy change caused by ligand binding compare favorably to values obtained from a dynamic Normal Mode Analysis (NMA). Higher correlation factors can be achieved by increasing the distance cutoff in the potential energy model. Furthermore, we apply our new method to analyze the entropy changes of the SARS CoV-2 main protease when binding with different ligand inhibitors, which is relevant for the design of potential drugs.

Dynamic Allostery Highlights the Evolutionary Differences between the CoV-1 and CoV-2 Main Proteases

The SARS-CoV-2 coronavirus has become one of the most immediate and widely-studied systems since its identification and subsequent global outbreak from 2019-2021. In an effort to understand the biophysical changes as a result of mutations, the mechanistic details of multiple different proteins within the SARS-CoV-2 virus have been studied and compared with SARS-CoV-1. Focusing on the main protease (mPro), we first explored the long-range dynamics using the Dynamic Coupling Index (DCI) to investigate the dynamic coupling between the catalytic site residues and the rest of the protein, both inter and intra chain, for the CoV-1 and CoV-2 mPro. We found that there is significant cross-chain coupling between these active sites and specific distal residues in the CoV-2 mPro not present in CoV-1. The enhanced long distance interactions, particularly between the two chains, suggest subsequently enhanced cooperativity for CoV-2. A further comparative analysis of the dynamic flexibility using the Dynamic Flexibility Index (DFI) between the CoV-1 and CoV-2 mPros shows that the inhibitor binding near active sites induces change in flexibility to a distal region of the protein, opposite in behavior between the two systems; this region becomes more flexible upon inhibitor binding in CoV-1 while it becomes less flexible in the CoV-2 mPro. Upon inspection, we show that, on average, the dynamic flexibility of the sites substituted from CoV-1 to CoV-2 changes significantly less than the average calculated across all residues within the structure, indicating that the differences in behaviors between the two systems is likely the result of allosteric influence, where the new substitutions in CoV-2 induce flexibility and dynamical changes elsewhere in the structure.

Gaussian network model revisited: Effects of mutation and ligand binding on protein behavior

The coarse-grained Gaussian Network model, GNM, considers only the alpha carbons of the folded protein. Therefore it is not directly applicable to the study of mutation or ligand binding problems where atomic detail is required. This shortcoming is improved by including all atom pairs within the coordination shell of each other into the Kirchoff Adjacency Matrix. Counting all contacts rather than only alpha carbon contacts diminishes the magnitude of fluctuations in the system. But more importantly, it changes the graph-like connectivity structure, i.e., the Kirchoff Adjacency Matrix of the protein. This change depends on amino acid type which introduces amino acid specific and position specific information into the classical coarse-grained GNM which was originally modelled in analogy with the phantom network model of rubber elasticity. With this modification, it is now possible to explain the consequences of mutation and ligand binding on residue fluctuations, their pair-correlations and mutual information (MI) shared by each pair. We refer to the new model as 'all-atom GNM'. Using examples from published data we show that the all-atom GNM gives B-factors that are in better agreement with experiment, can explain effects of mutation on long range communication in PDZ domains and can predict effects of GDP and GTP binding on the dimerization of KRAS.

Structural deformability induced in proteins of potential interest associated with COVID-19 by binding of homologues present in ivermectin: Comparative study based in elastic networks models

The COVID-19 pandemic has accelerated the study of the potential of multi-target drugs (MTDs). The mixture of homologues called ivermectin (avermectin-B1a + avermectin-B1b) has been shown to be a MTD with potential antiviral activity against SARS-CoV-2 in vitro. However, there are few reports on the effect of each homologue on the flexibility and stiffness of proteins associated with COVID-19, described as ivermectin targets. We observed that each homologue was stably bound to the proteins studied and was able to induce detectable changes with Elastic Network Models (ENM). The perturbations induced by each homologue were characteristic of each compound and, in turn, were represented by a disruption of native intramolecular networks (interactions between residues). The homologues were able to slightly modify the conformation and stability of the connection points between the Cα atoms of the residues that make up the structural network of proteins (nodes), compared to free proteins. Each homologue was able to modified differently the distribution of quasi-rigid regions of the proteins, which could theoretically alter their biological activities. These results could provide a biophysical-computational view of the potential MTD mechanism that has been reported for ivermectin.

Natural Polyphenols Inhibit the Dimerization of the SARS-CoV-2 Main Protease: The Case of Fortunellin and Its Structural Analogs

3CL-Pro is the SARS-CoV-2 main protease (MPro). It acts as a homodimer to cleave the large polyprotein 1ab transcript into proteins that are necessary for viral growth and replication. 3CL-Pro has been one of the most studied SARS-CoV-2 proteins and a main target of therapeutics. A number of drug candidates have been reported, including natural products. Here, we employ elaborate computational methods to explore the dimerization of the 3CL-Pro protein, and we formulate a computational context to identify potential inhibitors of this process. We report that fortunellin (acacetin 7-O-neohesperidoside), a natural flavonoid O-glycoside, and its structural analogs are potent inhibitors of 3CL-Pro dimerization, inhibiting viral plaque formation in vitro. We thus propose a novel basis for the search of pharmaceuticals as well as dietary supplements in the fight against SARS-CoV-2 and COVID-19.

Effectiveness of Natural Antioxidants against SARS-CoV-2? Insights from the In-Silico World

The SARS CoV-2 pandemic has affected millions of people around the globe. Despite many efforts to find some effective medicines against SARS CoV-2, no established therapeutics are available yet. The use of phytochemicals as antiviral agents provides hope against the proliferation of SARS-CoV-2. Several natural compounds were analyzed by virtual screening against six SARS CoV-2 protein targets using molecular docking simulations in the present study. More than a hundred plant-derived secondary metabolites have been docked, including alkaloids, flavonoids, coumarins, and steroids. SARS CoV-2 protein targets include Main protease (MPro), Papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp), Spike glycoprotein (S), Helicase (Nsp13), and E-Channel protein. Phytochemicals were evaluated by molecular docking, and MD simulations were performed using the YASARA structure using a modified genetic algorithm and AMBER03 force field. Binding energies and dissociation constants allowed the identification of potentially active compounds. Ligand-protein interactions provide an insight into the mechanism and potential of identified compounds. Glycyrrhizin and its metabolite 18-β-glycyrrhetinic acid have shown a strong binding affinity for MPro, helicase, RdRp, spike, and E-channel proteins, while a flavonoid Baicalin also strongly binds against PLpro and RdRp. The use of identified phytochemicals may help to speed up the drug development and provide natural protection against SARS-CoV-2.

Accelerating COVID-19 Research Using Molecular Dynamics Simulation

The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.

Elucidation of Cryptic and Allosteric Pockets within the SARS-CoV-2 Main Protease

The SARS-CoV-2 pandemic has rapidly spread across the globe, posing an urgent health concern. Many quests to computationally identify treatments against the virus rely on in silico small molecule docking to experimentally determined structures of viral proteins. One limit to these approaches is that protein dynamics are often unaccounted for, leading to overlooking transient, druggable conformational states. Using Gaussian accelerated molecular dynamics to enhance sampling of conformational space, we identified cryptic pockets within the SARS-CoV-2 main protease, including some within regions far from the active site. These simulations sampled comparable dynamics and pocket volumes to conventional brute force simulations carried out on two orders of magnitude greater timescales.

Dynamic Profiling of β-Coronavirus 3CL Mpro Protease Ligand-Binding Sites

β-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensable role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand-protein interactions. Herein, we provide a comprehensive comparison of all nonredundant ligand-binding sites available for SARS-CoV2, SARS-CoV, and MERS-CoV Mpro. Extensive adaptive sampling has been used to investigate structural conservation of ligand-binding sites using Markov state models (MSMs) and compare conformational dynamics employing convolutional variational auto-encoder-based deep learning. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across β-CoV homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor.

Molecular dynamics and in silico mutagenesis on the reversible inhibitor-bound SARS-CoV-2 main protease complexes reveal the role of lateral pocket in enhancing the ligand affinity

The 2019 novel coronavirus pandemic caused by SARS-CoV-2 remains a serious health threat to humans and there is an urgent need to develop therapeutics against this deadly virus. Recent scientific evidences have suggested that the main protease (Mpro) enzyme in SARS-CoV-2 can be an ideal drug target due to its crucial role in the viral replication and transcription processes. Therefore, there are ongoing research efforts to identify drug candidates against SARS-CoV-2 Mpro that resulted in hundreds of X-ray crystal structures of ligand-bound Mpro complexes in the Protein Data Bank (PDB) describing the interactions of different fragment chemotypes within different sites of the Mpro. In this work, we performed rigorous molecular dynamics (MD) simulation of 62 reversible ligand-Mpro complexes in the PDB to gain mechanistic insights about their interactions at the atomic level. Using a total of over 3 µs long MD trajectories, we characterized different pockets in the apo Mpro structure, and analyzed the dynamic interactions and binding affinity of ligands within those pockets. Our results identified the key residues that stabilize the ligands in the catalytic sites and other pockets of Mpro. Our analyses unraveled the role of a lateral pocket in the catalytic site in Mpro that is critical for enhancing the ligand binding to the enzyme. We also highlighted the important contribution from HIS163 in the lateral pocket towards ligand binding and affinity against Mpro through computational mutation analyses. Further, we revealed the effects of explicit water molecules and Mpro dimerization in the ligand association with the target. Thus, comprehensive molecular-level insights gained from this work can be useful to identify or design potent small molecule inhibitors against SARS-CoV-2 Mpro.

The impact of structural bioinformatics tools and resources on SARS-CoV-2 research and therapeutic strategies

SARS-CoV-2 is the causative agent of COVID-19, the ongoing global pandemic. It has posed a worldwide challenge to human health as no effective treatment is currently available to combat the disease. Its severity has led to unprecedented collaborative initiatives for therapeutic solutions against COVID-19. Studies resorting to structure-based drug design for COVID-19 are plethoric and show good promise. Structural biology provides key insights into 3D structures, critical residues/mutations in SARS-CoV-2 proteins, implicated in infectivity, molecular recognition and susceptibility to a broad range of host species. The detailed understanding of viral proteins and their complexes with host receptors and candidate epitope/lead compounds is the key to developing a structure-guided therapeutic design. Since the discovery of SARS-CoV-2, several structures of its proteins have been determined experimentally at an unprecedented speed and deposited in the Protein Data Bank. Further, specialized structural bioinformatics tools and resources have been developed for theoretical models, data on protein dynamics from computer simulations, impact of variants/mutations and molecular therapeutics. Here, we provide an overview of ongoing efforts on developing structural bioinformatics tools and resources for COVID-19 research. We also discuss the impact of these resources and structure-based studies, to understand various aspects of SARS-CoV-2 infection and therapeutic development. These include (i) understanding differences between SARS-CoV-2 and SARS-CoV, leading to increased infectivity of SARS-CoV-2, (ii) deciphering key residues in the SARS-CoV-2 involved in receptor-antibody recognition, (iii) analysis of variants in host proteins that affect host susceptibility to infection and (iv) analyses facilitating structure-based drug and vaccine design against SARS-CoV-2.

Impact of dimerization and N3 binding on molecular dynamics of SARS-CoV and SARS-CoV-2 main proteases

SARS-CoV-2 main protease is one of the major targets in drug development efforts against Covid-19. Even though several structures were reported to date, its dynamics is not understood well. In particular, impact of dimerization and ligand binding on the dynamics is an important issue to investigate. In this study, we performed molecular dynamics simulations of SARS-CoV and SARS-CoV-2 main proteases to investigate influence of dimerization on the dynamics by modeling monomeric and dimeric apo and holo forms. The dimerization causes an organization of the interdomain dynamics as well as some local structural changes. Moreover, we investigated impact of a peptide mimetic (N3) on the dynamics of SARS-CoV and SARS-CoV-2 Mpro. The ligand binding to the dimeric forms causes some key local changes at the dimer interface and it causes an allosteric interaction between the active sites of two protomers. Our results support the idea that only one protomer is active on SARS-CoV-2 due to this allosteric interaction. Additionally, we analyzed the molecular dynamics trajectories from graph theoretical perspective and found that the most influential residues – as measured by eigenvector centrality – are a group of residues in active site and dimeric interface of the protease. This study may form a bridge between what we know about the dynamics of SARS-CoV and SARS-CoV-2 Mpro. We think that enlightening allosteric communication of the active sites and the role of dimerization in SARS-CoV-2 Mpro can contribute to development of novel drugs against this global health problem as well as other similar proteases.

Communicated by Ramaswamy H. Sarma

A multiscale model of virus pandemic: Heterogeneous interactive entities in a globally connected world

This paper is devoted to the multidisciplinary modelling of a pandemic initiated by an aggressive virus, specifically the so-called SARS-CoV-2 Severe Acute Respiratory Syndrome, corona virus n.2. The study is developed within a multiscale framework accounting for the interaction of different spatial scales, from the small scale of the virus itself and cells, to the large scale of individuals and further up to the collective behaviour of populations. An interdisciplinary vision is developed thanks to the contributions of epidemiologists, immunologists and economists as well as those of mathematical modellers. The first part of the contents is devoted to understanding the complex features of the system and to the design of a modelling rationale. The modelling approach is treated in the second part of the paper by showing both how the virus propagates into infected individuals, successfully and not successfully recovered, and also the spatial patterns, which are subsequently studied by kinetic and lattice models. The third part reports the contribution of research in the fields of virology, epidemiology, immune competition, and economy focussed also on social behaviours. Finally, a critical analysis is proposed looking ahead to research perspectives.