Computational Study on the Inhibitory Effect of Natural Compounds against the SARS-CoV-2 Proteins.

COVID-19 is more virulent and challenging to human life. In India, the Ministry of AYUSH recommended some strategies through Siddha, homeopathy, and other methods to effectively manage COVID-19 (Guidelines for AYUSH Clinical Studies in COVID-19, 2020). Kabasura Kudineer and homeopathy medicines are in use for the prevention and treatment of COVID-19 infection; however, the mechanism of action is less explored. This study aims to understand the antagonist activity of natural compounds found in Kabasura Kudineer and homeopathy medicines against the SARS-CoV-2 using computational methods. Potential compounds were screened against NSP-12, NSP-13, NSP-14, NSP-15, main protease, and spike proteins. Structure-based virtual screening results shows that, out of 14,682 Kabasura Kudineer compounds, the 250395, 129677029, 44259583, 44259584, and 88583189 compounds and, out of 3,112 homeopathy compounds, the 3802778, 320361, 5315832, 14590080, and 74029795 compounds have good scoring function against the SARS-CoV-2 structural and nonstructural proteins. As a result of docking, homeopathy compounds have a docking score ranging from -5.636 to 13.631 kcal/mol, while Kabasura Kudineer compounds have a docking score varying from -8.290 to -13.759 kcal/mol. It has been found that the selected compounds bind well to the active site of SARS-CoV-2 proteins and form hydrogen bonds. The molecular dynamics simulation study shows that the selected compounds have maintained stable conformation in the simulation period and interact with the target. This study supports the antagonist activity of natural compounds from Kabasura Kudineer and homeopathy against SARS-CoV-2's structural and nonstructural proteins.

Detection of Binding Sites on SARS-CoV-2 Spike Protein Receptor-Binding Domain by Molecular Dynamics Simulations in Mixed Solvents.

The novel SARS-CoV-2 uses ACE2 (Angiotensin-Converting Enzyme 2) receptor as an entry point. Insights on S protein receptor-binding domain (RBD) interaction with ACE2 receptor and drug repurposing has accelerated drug discovery for the novel SARS-CoV-2 infection. Finding small molecule binding sites in S protein and ACE2 interface is crucial in search of effective drugs to prevent viral entry. In this study, we employed molecular dynamics simulations in mixed solvents together with virtual screening to identify small molecules that could be potential inhibitors of S protein -ACE2 interaction. Observation of organic probe molecule localization during the simulations revealed multiple sites at the S protein surface related to small molecule, antibody, and ACE2 binding. In addition, a novel conformation of the S protein was discovered that could be stabilized by small molecules to inhibit attachment to ACE2. The most promising binding site on RBD-ACE2 interface was targeted with virtual screening and top-ranked compounds (DB08248, DB02651, DB03714, and DB14826) are suggested for experimental testing. The protocol described here offers an extremely fast method for characterizing key proteins of a novel pathogen and for the identification of compounds that could inhibit or accelerate spreading of the disease.

Fragment- and negative image-based screening of phosphodiesterase 10A inhibitors.

A novel virtual screening methodology called fragment- and negative image-based (F-NiB) screening is introduced and tested experimentally using phosphodiesterase 10A (PDE10A) as a case study. Potent PDE10A-specific small-molecule inhibitors are actively sought after for their antipsychotic and neuroprotective effects. The F-NiB combines features from both fragment-based drug discovery and negative image-based (NIB) screening methodologies to facilitate rational drug discovery. The selected structural parts of protein-bound ligand(s) are seamlessly combined with the negative image of the target's ligand-binding cavity. This cavity- and fragment-based hybrid model, namely its shape and electrostatics, is used directly in the rigid docking of ab initio generated ligand 3D conformers. In total, 14 compounds were acquired using the F-NiB methodology, 3D quantitative structure-activity relationship modeling, and pharmacophore modeling. Three of the small molecules inhibited PDE10A at ~27 to ~67 µM range in a radiometric assay. In a larger context, the study shows that the F-NiB provides a flexible way to incorporate small-molecule fragments into the drug discovery.

Getting Docking into Shape Using Negative Image-Based Rescoring.

The failure of default scoring functions to ensure virtual screening enrichment is a persistent problem for the molecular docking algorithms used in structure-based drug discovery. To remedy this problem, elaborate rescoring and postprocessing schemes have been developed with a varying degree of success, specificity, and cost. The negative image-based rescoring (R-NiB) has been shown to improve the flexible docking performance markedly with a variety of drug targets. The yield improvement is achieved by comparing the alternative docking poses against the negative image of the target protein's ligand-binding cavity. In other words, the shape and electrostatics of the binding pocket is directly used in the similarity comparison to rank the explicit docking poses. Here, the PANTHER/ShaEP-based R-NiB methodology is tested with six popular docking softwares, including GLIDE, PLANTS, GOLD, DOCK, AUTODOCK, and AUTODOCK VINA, using five validated benchmark sets. Overall, the results indicate that R-NiB outperforms the default docking scoring consistently and inexpensively, demonstrating that the methodology is ready for wide-scale virtual screening usage.

Ultrafast protein structure-based virtual screening with Panther.

Molecular docking is by far the most common method used in protein structure-based virtual screening. This paper presents Panther, a novel ultrafast multipurpose docking tool. In Panther, a simple shape-electrostatic model of the ligand-binding area of the protein is created by utilizing the protein crystal structure. The features of the possible ligands are then compared to the model by using a similarity search algorithm. On average, one ligand can be processed in a few minutes by using classical docking methods, whereas using Panther processing takes <1 s. The presented Panther protocol can be used in several applications, such as speeding up the early phases of drug discovery projects, reducing the number of failures in the clinical phase of the drug development process, and estimating the environmental toxicity of chemicals. Panther-code is available in our web pages (http://www.jyu.fi/panther) free of charge after registration.

Comparison of virtual high-throughput screening methods for the identification of phosphodiesterase-5 inhibitors.

Reliable and effective virtual high-throughput screening (vHTS) methods are desperately needed to minimize the expenses involved in drug discovery projects. Here, we present an improvement to the negative image-based (NIB) screening: the shape, the electrostatics, and the solvation state of the target protein's ligand-binding site are included into the vHTS. Additionally, the initial vHTS results are postprocessed with molecular mechanics/generalized Born surface area (MMGBSA) calculations to estimate the favorability of ligand-protein interactions. The results show that docking produces very good early enrichment for phosphodiesterase-5 (PDE-5); however, in general, the NIB and the ligand-based screening performed better with or without the added electrostatics. Furthermore, the postprocessing of the NIB screening results using MMGBSA calculations improved the early enrichment for the PDE-5 considerably, thus, making hit discovery affordable.

Efficient virtual screening using multiple protein conformations described as negative images of the ligand-binding site.

The protein structure-based virtual screening is typically accomplished using a molecular docking procedure. However, docking is a fairly slow process that is limited by the available scoring functions that cannot reliably distinguish between active and inactive ligands. In contrast, the ligand-based screening methods that are based on shape similarity identify the active ligands with high accuracy. Here, we show that the usage of negative images of the ligand-binding site, together with shape comparison tools, which are typically used in ligand-based virtual screening, improve the discrimination of active molecules from inactives. In contrast to ligand-based shape comparison, the negative image of the binding site allows identification of compounds whose shape complements the shape of the ligand-binding cavity as closely as possible. Furthermore, the use of several target protein conformations allows the identification of active ligands whose shape is not optimal for crystallized protein conformation. Accordingly, the presented virtual screening method improves the identification of novel lead molecules by concentrating on the optimally shaped molecules for the flexible ligand binding site.

alpha 11beta 1 integrin recognizes the GFOGER sequence in interstitial collagens.

The integrins alpha(1)beta(1), alpha(2)beta(1), alpha(10)beta(1), and alpha(11)beta(1) are referred to as a collagen receptor subgroup of the integrin family. Recently, both alpha(1)beta(1) and alpha(2)beta(1) integrins have been shown to recognize triple-helical GFOGER (where single letter amino acid nomenclature is used, O = hydroxyproline) or GFOGER-like motifs found in collagens, despite their distinct binding specificity for various collagen subtypes. In the present study we have investigated the mechanism whereby the latest member in the integrin family, alpha(11)beta(1), recognizes collagens using C2C12 cells transfected with alpha(11) cDNA and the bacterially expressed recombinant alpha(11) I domain. The ligand binding properties of alpha(11)beta(1) were compared with those of alpha(2)beta(1). Mg(2+)-dependent alpha(11)beta(1) binding to type I collagen required micromolar Ca(2+) but was inhibited by 1 mm Ca(2+), whereas alpha(2)beta(1)-mediated binding was refractory to millimolar concentrations of Ca(2+). The bacterially expressed recombinant alpha(11) I domain preference for fibrillar collagens over collagens IV and VI was the same as the alpha(2) I domain. Despite the difference in Ca(2+) sensitivity, alpha(11)beta(1)-expressing cells and the alpha(11) I domain bound to helical GFOGER sequences in a manner similar to alpha(2)beta(1)-expressing cells and the alpha(2) I domain. Modeling of the alpha I domain-collagen peptide complexes could partially explain the observed preference of different I domains for certain GFOGER sequence variations. In summary, our data indicate that the GFOGER sequence in fibrillar collagens is a common recognition motif used by alpha(1)beta(1), alpha(2)beta(1), and also alpha(11)beta(1) integrins. Although alpha(10) and alpha(11) chains show the highest sequence identity, alpha(2) and alpha(11) are more similar with regard to collagen specificity. Future studies will reveal whether alpha(2)beta(1) and alpha(11)beta(1) integrins also show overlapping biological functions.