Reversal of the unique Q493R mutation increases the affinity of Omicron S1-RBD for ACE2

The SARS-CoV-2 Omicron variant containing 15 mutations, including the unique Q493R, in the spike protein receptor binding domain (S1-RBD) is highly infectious. While comparison with previously reported mutations provide some insights, the mechanism underlying the increased infections and the impact of the reversal of the unique Q493R mutation seen in BA.4, BA.5, BA.2.75, BQ.1 and XBB lineages is not yet completely understood. Here, using structural modelling and molecular dynamics (MD) simulations, we show that the Omicron mutations increases the affinity of S1-RBD for ACE2, and a reversal of the unique Q493R mutation further increases the ACE2-S1-RBD affinity. Specifically, we performed all atom, explicit solvent MD simulations using a modelled structure of the Omicron S1-RBD-ACE2 and compared the trajectories with the WT complex revealing a substantial reduction in the Cα-atom fluctuation in the Omicron S1-RBD and increased hydrogen bond and other interactions. Residue level analysis revealed an alteration in the interaction between several residues including a switch in the interaction of ACE2 D38 from S1-RBD Y449 in the WT complex to the mutated R residue (Q493R) in Omicron complex. Importantly, simulations with Revertant (Omicron without the Q493R mutation) complex revealed further enhancement of the interaction between S1-RBD and ACE2. Thus, results presented here not only provide insights into the increased infectious potential of the Omicron variant but also a mechanistic basis for the reversal of the Q493R mutation seen in some Omicron lineages and will aid in understanding the impact of mutations in SARS-CoV-2 evolution.

Phytochemical rich Himalayan Rhododendron arboreum petals inhibit SARS-CoV-2 infection in vitro.

Phytochemicals with potential to competitively bind to the host receptors or inhibit SARS-CoV-2 replication, may prove to be useful as adjunct therapeutics for COVID-19. We profiled and investigated the phytochemicals of Rhododendron arboreum petals sourced from Himalayan flora, undertook in vitro studies and found it as a promising candidate against SARS-CoV-2. The phytochemicals were reported in various scientific investigations to act against a range of virus in vitro and in vivo, which prompted us to test against SARS-CoV-2. In vitro assays of R. arboreum petals hot aqueous extract confirmed dose dependent reduction in SARS-CoV-2 viral load in infected Vero E6 cells (80% inhibition at 1 mg/ml; IC50 = 173 µg/ml) and phytochemicals profiled were subjected to molecular docking studies against SARS CoV-2 target proteins. The molecules 5-O-Feruloyl-quinic acid, 3-Caffeoyl-quinic acid, 5-O-Coumaroyl-D-quinic acid, Epicatechin and Catechin showed promising binding affinity with SARS-CoV-2 Main protease (MPro; PDB ID: 6LU7; responsible for viral replication) and Human Angiotensin Converting Enzyme-2 (ACE2; PDB ID: 1R4L; mediate viral entry in the host). Molecular dynamics (MD) simulation of 5-O-Feruloyl-quinic acid, an abundant molecule in the extract complexed with the target proteins showed stable interactions. Taken together, the phytochemical profiling, in silico analysis and in vitro anti-viral assay revealed that the petals extract act upon MPro and may be inhibiting SARS-CoV-2 replication. This is the first report highlighting R. arboreum petals as a reservoir of antiviral phytochemicals with potential anti-SARS-CoV-2 activity using an in vitro system.

Synthesis, X-ray crystal structure, Hirshfeld surface analysis and computational investigation into the potential inhibitory action of novel 6-(p-tolyl)-2-((p-tolyl)thio)methyl-7H-[1.2.4]triazolo[5,1-b][1,3,4]thiadiazine inhibits the main protease of COVID-19.

In recent times, the novel coronavirus disease (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has now become a worldwide pandemic. With over 71 million confirmed cases, even though the effectiveness and side effects of the specific drugs and vaccines approved for this disease are still limited. Scientists and researchers from all across the world are working to find a vaccine and a cure for COVID-19 by using large-scale drug discovery and analysis. Heterocyclic compounds are regarded to be valuable sources for the discovery of new antiviral medications against SARS-CoV-2 because virus occurrences are still on the rise, and infectivity and mortality may also rise shortly. In this regard, we have synthesized a new triazolothiadiazine derivative. The structure was characterized by NMR spectra and confirmed by X-ray diffraction analysis. The structural geometry coordinates of the title compound are well reproduced by DFT calculations. NBO and NPA analyses have been performed to determine the interaction energies between bonding and antibonding orbital, and natural atomic charges of heavy atoms. Molecular docking suggests that the compounds may have good affinity for SAR-CoV-2 main protease, RNA-dependent RNA polymerase and nucleocapsid enzymes, particularly the main protease enzyme (binding energy of -11.9 kcal mol-1). The predicted docked pose of the compound is dynamically stable and reports a major van der Waals contribution (-62.00 kcal mol-1) to overall net energy.Communicated by Ramaswamy H. Sarma.

Modelling the Transitioning of SARS-CoV-2 nsp3 and nsp4 Lumenal Regions towards a More Stable State on Complex Formation.

During coronavirus infection, three non-structural proteins, nsp3, nsp4, and nsp6, are of great importance as they induce the formation of double-membrane vesicles where the replication and transcription of viral gRNA takes place, and the interaction of nsp3 and nsp4 lumenal regions triggers membrane pairing. However, their structural states are not well-understood. We investigated the interactions between nsp3 and nsp4 by predicting the structures of their lumenal regions individually and in complex using AlphaFold2 as implemented in ColabFold. The ColabFold prediction accuracy of the nsp3-nsp4 complex was increased compared to nsp3 alone and nsp4 alone. All cysteine residues in both lumenal regions were modelled to be involved in intramolecular disulphide bonds. A linker region in the nsp4 lumenal region emerged as crucial for the interaction, transitioning to a structured state when predicted in complex. The key interactions modelled between nsp3 and nsp4 appeared stable when the transmembrane regions of nsp3 and nsp4 were added to the modelling either alone or together. While molecular dynamics simulations (MD) demonstrated that the proposed model of the nsp3 lumenal region on its own is not stable, key interactions between nsp and nsp4 in the proposed complex model appeared stable after MD. Together, these observations suggest that the interaction is robust to different modelling conditions. Understanding the functional importance of the nsp4 linker region may have implications for the targeting of double membrane vesicle formation in controlling coronavirus infection.

Identifying non-nucleoside inhibitors of RNA-dependent RNA-polymerase of SARS-CoV-2 through per-residue energy decomposition-based pharmacophore modeling, molecular docking, and molecular dynamics simulation.

BACKGROUND AND OBJECTIVE: The current coronavirus disease-2019 (COVID-19) pandemic has triggered a worldwide health and economic crisis. The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes the disease and completes its life cycle using the RNA-dependent RNA-polymerase (RdRp) enzyme, a prominent target for antivirals. In this study, we have computationally screened ∼690 million compounds from the ZINC20 database and 11,698 small molecule inhibitors from DrugBank to find existing and novel non-nucleoside inhibitors for SARS-CoV-2 RdRp. METHODS: Herein, a combination of the structure-based pharmacophore modeling and hybrid virtual screening methods, including per-residue energy decomposition-based pharmacophore screening, molecular docking, pharmacokinetics, and toxicity evaluation were employed to retrieve novel as well as existing RdRp non-nucleoside inhibitors from large chemical databases. Besides, molecular dynamics simulation and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) method were used to investigate the binding stability and calculate the binding free energy of RdRp-inhibitor complexes. RESULTS: Based on docking scores and significant binding interactions with crucial residues (Lys553, Arg557, Lys623, Cys815, and Ser816) in the RNA binding site of RdRp, three existing drugs, ZINC285540154, ZINC98208626, ZINC28467879, and five compounds from ZINC20 (ZINC739681614, ZINC1166211307, ZINC611516532, ZINC1602963057, and ZINC1398350200) were selected, and the conformational stability of RdRp due to their binding was confirmed through molecular dynamics simulation. The free energy calculations revealed these compounds possess strong binding affinities for RdRp. In addition, these novel inhibitors exhibited drug-like features, good absorption, distribution, metabolism, and excretion profile and were found to be non-toxic. CONCLUSION: The compounds identified in the study by multifold computational strategy can be validated in vitro as potential non-nucleoside inhibitors of SARS-CoV-2 RdRp and holds promise for the discovery of novel drugs against COVID-19 in future.

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations.

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.

A new approach for extracting information from protein dynamics.

Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.

An insight into SARS-CoV-2 membrane protein interaction with spike, envelope, and nucleocapsid proteins.

Intraviral protein-protein interactions are crucial for replication, pathogenicity, and viral assembly. Among these, virus assembly is a critical step as it regulates the arrangements of viral structural proteins and helps in the encapsulation of genomic material. SARS-CoV-2 structural proteins play an essential role in the self-rearrangement, RNA encapsulation, and mature virus particle formation. In SARS-CoV, the membrane protein interacts with the envelope and spike protein in Endoplasmic Reticulum Golgi Intermediate Complex (ERGIC) to form an assembly in the lipid bilayer, followed by membrane-ribonucleoprotein (nucleocapsid) interaction. In this study, we tried to understand the interaction of membrane protein's interaction with envelope, spike, and nucleocapsid proteins using protein-protein docking. Further, simulation studies were performed up to 100 ns to examine the stability of protein-protein complexes of Membrane-Envelope, Membrane-Spike, and Membrane-Nucleocapsid proteins. Prime MM-GBSA showed high binding energy calculations for the simulated structures than the docked complex. The interactions identified in our study will be of great importance, as it provides valuable insight into the protein-protein complex, which could be the potential drug targets for future studies.Communicated by Ramaswamy H. Sarma.

Aloin A inhibits SARS CoV-2 replication by targeting its binding with ACE2 - Evidence from modeling-supported molecular dynamics simulation.

The current study aimed to expand on the recently published results and assess the inhibitory efficacy of aloin A against SARS CoV-2. In vitro testing of aloin A against SARS CoV-2 proteases (i.e., MPro and PLPro) showed weak to moderate activity (IC50 = 68.56 ± 1.13 µM and 24.77 ± 1.57 µM, respectively). However, aloin A was able to inhibit the replication of SARS CoV-2 in Vero E6 cells efficiently with an IC50 of 0.095 ± 0.022 µM. Depending on the reported poor permeability of aloin A alongside its insignificant protease inhibitory activities presented in this study, we ran a number of extensive virtual screenings and physics-based simulations to determine the compound's potential mode of action. As a result, RBD-ACE2 was identified as a key target for aloin A. Results from 600 ns-long molecular dynamics (MD) simulation experiments pointed to aloin A's role as an RBD-ACE2 destabilizer. Therefore, the results of this work may pave the way for further development of this scaffold and the eventual production of innovative anti-SARS CoV-2 medicines with several mechanisms of action.Communicated by Ramaswamy H. Sarma.

Identification of Pyrazole Derivatives of Usnic Acid as Novel Inhibitor of SARS-CoV-2 Main Protease Through Virtual Screening Approaches.

The infection produced by the SARS-CoV-2 virus remains a significant health crisis worldwide. The lack of specific medications for COVID-19 necessitates a concerted effort to find the much-desired therapies for this condition. The main protease (Mpro) of SARS-CoV-2 is a promising target, vital for virus replication and transcription. In this study, fifty pyrazole derivatives were tested for their pharmacokinetics and drugability, resulting in eight hit compounds. Subsequent molecular docking simulations on SARS-CoV-2 main protease afforded two lead compounds with strong affinity at the active site. Additionally, the molecular dynamics (MD) simulations of lead compounds (17 and 39), along with binding free energy calculations, were accomplished to validate the stability of the docked complexes and the binding poses achieved in docking experiments. Based on these findings, compound 17 and 39, with their favorable projected pharmacokinetics and pharmacological characteristics, are the proposed potential antiviral candidates which require further investigation to be used as anti-SARS-CoV-2 medication.

Exploring novel targets of sitagliptin for type 2 diabetes mellitus: Network pharmacology, molecular docking, molecular dynamics simulation, and SPR approaches.

Background: Diabetes has become a serious global public health problem. With the increasing prevalence of type 2 diabetes mellitus (T2DM), the incidence of complications of T2DM is also on the rise. Sitagliptin, as a targeted drug of DPP4, has good therapeutic effect for T2DM. It is well known that sitagliptin can specifically inhibit the activity of DPP4 to promote insulin secretion, inhibit islet ß cell apoptosis and reduce blood glucose levels, while other pharmacological mechanisms are still unclear, such as improving insulin resistance, anti-inflammatory, anti-oxidative stress, and anti-fibrosis. The aim of this study was to explore novel targets and potential signaling pathways of sitagliptin for T2DM. Methods: Firstly, network pharmacology was applied to find the novel target most closely related to DPP4. Semi-flexible molecular docking was performed to confirm the binding ability between sitagliptin and the novel target, and molecular dynamics simulation (MD) was carried to verify the stability of the complex formed by sitagliptin and the novel target. Furthermore, surface-plasmon resonance (SPR) was used to explored the affinity and kinetic characteristics of sitagliptin with the novel target. Finally, the molecular mechanism of sitagliptin for T2DM was predicted by the enrichment analysis of GO function and KEGG pathway. Results: In this study, we found the cell surface receptor-angiotensin-converting enzyme 2 (ACE2) most closely related to DPP4. Then, we confirmed that sitagliptin had strong binding ability with ACE2 from a static perspective, and the stability of sitagliptin-ACE2 complex had better stability and longer binding time than BAR708-ACE2 in simulated aqueous solution within 50 ns. Significantly, we have demonstrated a strong affinity between sitagliptin and ACE2 on SPR biosensor, and their kinetic characteristics were "fast binding/fast dissociation". The guiding significance of clinical administration: low dose can reach saturation, but repeated administration was needed. Finally, there was certain relationship between COVID-19 and T2DM, and ACE2/Ang-(1-7)/Mas receptor (MasR) axis may be the important pathway of sitagliptin targeting ACE2 for T2DM. Conclusion: This study used different methods to prove that ACE2 may be another novel target of sitagliptin for T2DM, which extended the application of ACE2 in improving diabetes mellitus.

Core fucose-specific Pholiota squarrosa lectin (PhoSL) as a potent broad-spectrum inhibitor of SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S protein) is highly N-glycosylated, and a "glycan shield" is formed to limit the access of other molecules; however, a small open area coincides with the interface to the host's receptor and also neutralising antibodies. Most of the variants of concern have mutations in this area, which could reduce the efficacy of existing antibodies. In contrast, N-glycosylation sites are relatively invariant, and some are essential for infection. Here, we observed that the S proteins of the ancestral (Wuhan) and Omicron strains bind with Pholiota squarrosa lectin (PhoSL), a 40-amino-acid chemically synthesised peptide specific to core-fucosylated N-glycans. The affinities were at a low nanomolar level, which were ~ 1000-fold stronger than those between PhoSL and the core-fucosylated N-glycans at the micromolar level. We demonstrated that PhoSL inhibited infection by both strains at similar submicromolar levels, suggesting its broad-spectrum effect on SARS-CoV-2 variants. Cryogenic electron microscopy revealed that PhoSL caused an aggregation of the S protein, which was likely due to the multivalence of both the trimeric PhoSL and S protein. This characteristic is likely relevant to the inhibitory mechanism. Structural modelling of the PhoSL-S protein complex indicated that PhoSL was in contact with the amino acids of the S protein, which explains the enhanced affinity with S protein and also indicates the significant potential for developing specific binders by the engineering of PhoSL.

Predictive modeling and therapeutic repurposing of natural compounds against the receptor-binding domain of SARS-CoV-2.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a member of the Coronaviridae family, causing major destructions to human life directly and indirectly to the economic crisis around the world. Although there is significant reporting on the whole genome sequences and updated data for the different receptors are widely analyzed and screened to find a proper medication. Only a few bioassay experiments were completed against SARS-CoV-2 spike protein. We collected the compounds dataset from the PubChem Bioassay database having 1786 compounds and split it into the ratio of 80-20% for model training and testing purposes, respectively. Initially, we have created 11 models and validated them using a fivefold validation strategy. The hybrid consensus model shows a predictive accuracy of 95.5% for training and 94% for the test dataset. The model was applied to screen a virtual chemical library of Natural products of 2598 compounds. Our consensus model has successfully identified 75 compounds with an accuracy range of 70-100% as active compounds against SARS-CoV-2 RBD protein. The output of ML data (75 compounds) was taken for the molecular docking and dynamics simulation studies. In the complete analysis, the Epirubicin and Daunorubicin have shown the docking score of -9.937 and -9.812, respectively, and performed well in the molecular dynamics simulation studies. Also, Pirarubicin, an analogue of anthracycline, has widely been used due to its lower cardiotoxicity. It shows the docking score of -9.658, which also performed well during the complete analysis. Hence, after the following comprehensive pipeline-based study, these drugs can be further tested in vivo for further human utilization.Communicated by Ramaswamy H. Sarma.

DFT, ADMET, molecular docking and molecular dynamics studies of pyridoxal

We report in silico studies of pyridoxal, which is of interest both as a precursor for further functionalization due to the presence of the aldehyde functionality, as well as a bioactive compound. So far, the crystal structure of pyridoxal has not been reported and, thus, we have optimized its structure both under water solvation and in gas phase using the DFT calculations. Under water solvation conditions the optimized structure of pyridoxal is 7.62 kcal/mol more favorable in comparison to that in gas phase. The DFT calculations were also applied to verify optical and electronic properties of the optimized structure of pyridoxal in water. The HOMO and LUMO were revealed to subtract a set of descriptors of the so-called global chemical reactivity as well as to probe pyridoxal as a potential corrosion inhibitor for some important metals used in implants. The title compound exhibits the best electron charge transfer from the molecule to the surface of Ni and Co. Some biological properties of pyridoxal were evaluated using the respective on-line tools. Molecular docking was additionally applied to study interaction of pyridoxal with some SARS-CoV-2 proteins as well as one of the monkeypox proteins. It was established that the title compound is active against all the applied proteins with the most efficient interaction with nonstructural protein 15 (endoribonuclease) and Omicron Spike protein of SARS-CoV-2. Pyridoxal was found to be also active against the studied monkeypox protein. Interaction of pyridoxal with nonstructural protein 15 (endoribonuclease) was further studied using molecular dynamics simulation.

Unveiling the Potentiality of Shikonin Derivatives Inhibiting SARS-CoV-2 Main Protease by Molecular Dynamic Simulation Studies.

Shikonin, a phytochemical present in the roots of Lithospermum erythrorhizon, is well-known for its broad-spectrum activity against cancer, oxidative stress, inflammation, viruses, and anti-COVID-19 agents. A recent report based on a crystallographic study revealed a distinct conformation of shikonin binding to the SARS-CoV-2 main protease (Mpro), suggesting the possibility of designing potential inhibitors based on shikonin derivatives. The present study aimed to identify potential shikonin derivatives targeting the Mpro of COVID-19 by using molecular docking and molecular dynamics simulations. A total of 20 shikonin derivatives were screened, of which few derivatives showed higher binding affinity than shikonin. Following the MM-GBSA binding energy calculations using the docked structures, four derivatives were retained with the highest binding energy and subjected to molecular dynamics simulation. Molecular dynamics simulation studies suggested that alpha-methyl-n-butyl shikonin, beta-hydroxyisovaleryl shikonin, and lithospermidin-B interacted with two conserved residues, His41 and Cys145, through multiple bonding in the catalytic sites. This suggests that these residues may effectively suppress SARS-CoV-2 progression by inhibiting Mpro. Taken together, the present in silico study concluded that shikonin derivatives may play an influential role in Mpro inhibition.

Computational design of candidate multi-epitope vaccine against SARS-CoV-2 targeting structural (S and N) and non-structural (NSP3 and NSP12) proteins.

The COVID-19 pandemic caused by SARS-CoV-2 virus has created a global damage and has exposed the vulnerable side of scientific research towards novel diseases. The intensity of the pandemic is huge, with mortality rates of more than 6 million people worldwide in a span of 2 years. Considering the gravity of the situation, scientists all across the world are continuously attempting to create successful therapeutic solutions to combat the virus. Various vaccination strategies are being devised to ensure effective immunization against SARS-CoV-2 infection. SARS-CoV-2 spreads very rapidly, and the infection rate is remarkably high than other respiratory tract viruses. The viral entry and recognition of the host cell is facilitated by S protein of the virus. N protein along with NSP3 is majorly responsible for viral genome assembly and NSP12 performs polymerase activity for RNA synthesis. In this study, we have designed a multi-epitope, chimeric vaccine considering the two structural (S and N protein) and two non-structural proteins (NSP3 and NSP12) of SARS-CoV-2 virus. The aim is to induce immune response by generating antibodies against these proteins to target the viral entry and viral replication in the host cell. In this study, computational tools were used, and the reliability of the vaccine was verified using molecular docking, molecular dynamics simulation and immune simulation studies in silico. These studies demonstrate that the vaccine designed shows steady interaction with Toll like receptors with good stability and will be effective in inducing a strong and specific immune response in the body.Communicated by Ramaswamy H. Sarma.

Computational approaches to define poncirin from Magnolia champaka leaves as a novel multi-target inhibitor of SARS-CoV-2.

Phytochemical-based drug discovery against the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has been the focus of the current scenario. In this context, we aimed to perform the phytochemical profiling of Magnolia champaka, an evergreen tree from the Magnoliaceae family, in order to perform a virtual screening of its phytoconstituents against different biological targets of SARS-CoV-2. The phytochemicals identified from the ethanol extract of M. champaka leaves using liquid chromatography-mass spectroscopy (LC-MS) technique were screened against SARS-CoV-2 spike glycoprotein (PDB ID: 6M0J), main protease/Mpro (PDB ID: 6LU7), and papain-like protease/PLpro (PDB ID: 7CMD) through computational tools. The experimentation design included molecular docking simulation, molecular dynamics simulation, and binding free energy calculations. Through molecular docking simulation, we identified poncirin as a common potential inhibitor of all the above-mentioned target proteins. In addition, molecular dynamics simulations, binding free energy calculations, and PCA analysis also supported the outcomes of the virtual screening. By the virtue of all the in silico results obtained, poncirin could be taken for in vitro and in vivo studies in near future.Communicated by Ramaswamy H. Sarma.

Identification of potent COVID-19 main protease inhibitors by loading of favipiravir on Mg12O12 and Zn12O12 nanoclusters: an in silico strategy for COVID-19 treatment.

Pandemic new severe acute respiratory syndrome coronavirus (SARS-CoV-2) virus has increased throughout the world. There is no effective treatment against this virus until now. Since its appearance in Wuhan, China in December 2019, SARS-CoV-2 becomes the largest challenge the world is opposite today, including the discovery of an antiviral drug for this virus. Several viral proteins have been prioritized as SARS-CoV-2 antiviral drug targets, among them the papain-like protease (PLpro) and the main protease (Mpro). Inhibition of these proteases would target viral replication, viral maturation and suppression of host innate immune responses. Potential candidates have been identified to show inhibitory effects against Mpro, both in biochemical assays and viral replication in cells. There are different molecules such as lopinavir and favipiravir considerably inhibit the activity of Mpro in vitro. Different studies have shown that structurally improved favipiravir and other similar compounds can inhibit SARS-CoV-2 main protease. In this work, we study the interactions between favipiravir with Mg12O12 and Zn12O12 nanoclusters by density functional theory (DFT) and quantum mechanics atoms in molecules (QMAIM) methods to summarize the ability to load favipiravir onto Mg12O12 and Zn12O12 nanoclusters. Favipiravir-Mg12O12 and favipiravir-Zn12O12 lowest structures complexes were chosen to dock inside the SARS-CoV-2 main protease by molecular docking study. The molecular docking analysis revealed that the binding affinity of Mg12O12 and Zn12O12 nanoclusters inside the Mpro receptor is larger than that of favipiravir. Also, the loading of favipiravir on the surface of Mg12O12 and Zn12O12 nanoclusters increased the binding affinity against the Mpro receptor. Subsequently, 100 ns molecular dynamics simulation of the favipiravir-Mg12O12, and favipiravir-Zn12O12 docked inside the Mpro complexes established that favipiravir-Mg12O12, forms the most stable complex with the Mpro. Further molecular mechanics Poisson Boltzmann surface area (MMPBSA) analyses using the MD trajectories also demonstrated the higher binding affinity of favipiravir-Mg12O12 inside the Mpro. In summary, this study demonstrates a new way to characterize leads for novel anti-viral drugs against SARS-CoV-2, by improving the drug ability of favipiravir via loading it on Mg12O12 and Zn12O12 nanoclusters.Communicated by Ramaswamy H. Sarma.

Modeling Substrate Coordination to Zn-Bound Angiotensin Converting Enzyme 2

The Angiotensin Converting Enzyme 2 (ACE2) is a crucial regulator for the renin-angiotensin system. ACE2 converts the Angiotensin (Ang) II peptide into Ang 1-7 and thus promotes various anti-proliferative, anti-inflammatory and cardioprotective effects. In this study, we computationally designed several Ang II mutants to find a strong binding sequence to ACE2 receptor and examined the role of ligand substitution in the docking of native as well as mutant Ang II to the receptor. The peptide in the ACE2-peptide complex was coordinated to zinc (Zn) in the ACE2 cleft. The MD-generated root-mean-square deviation values were mostly similar between the native and mutant peptides considered in this work. The initial peptide-ACE2 poses were generated by molecular docking. The MD simulations used were post-processed by MM-PBSA to generate the binding free energies. All of the peptides studied here demonstrated negative binding free energies, which suggest that all the tested peptides form stable complexes with ACE2. Additionally, by examining the trends in the binding free energies calculated with different internal dielectric constants, it is evident that native Ang II and two of its variants have strongest binding to ACE2 receptor. Even though free energy measurements through classical MD simulation have certain limitations, in the absence of the availability of crystal structures of ACE2-peptide complexes, our work provides some structural insights for various Ang II analogs and how they may interact with a zinc atom within the active site of the enzyme. Copyright © 2022, The Author(s), under exclusive licence to Springer Nature B.V.

Virtual screening for finding inhibitor against the main protease of SARS-Cov-2 from the FDA-approved drugs database

At the end of 2019,a new coronavirus suddenly broke out all over the world.To date, there is still no targeted medicine available for the treatment of this disease. Vaccineis essential for controlling the epidemicofSARS-CoV-2. But the effective ofvaccine was reduced because of the SARS-CoV-2constant mutation. It is gratifying that scientistuncover theinfection mechanisms of the SARS-CoV-2. The main protease of SARS-CoV-2 is highly conserved and plays an important role of the life cycle of virus. Therefore, we executed virtual screening on the FDA-approved database and hoped to find a potential candidate against the main protease. As a result, we obtained eight available active compounds derived from the database through molecular dynamics simulations. As antiviral treatment candidates, the drugs can also be used to clinical emergencies. © 2022 SPIE. All rights reserved.