SARS-CoV Mpro inhibitory activity of aromatic disulfide compounds: QSAR model

Review and perspective on bioinformatics tools using machine learning and deep learning for predicting antiviral peptides

Viruses constitute a constant threat to global health and have caused millions of human and animal deaths throughout human history. Despite advances in the discovery of antiviral compounds that help fight these pathogens, finding a solution to this problem continues to be a task that consumes time and financial resources. Currently, artificial intelligence (AI) has revolutionized many areas of the biological sciences, making it possible to decipher patterns in amino acid sequences that encode different functions and activities. Within the field of AI, machine learning, and deep learning algorithms have been used to discover antimicrobial peptides. Due to their effectiveness and specificity, antimicrobial peptides (AMPs) hold excellent promise for treating various infections caused by pathogens. Antiviral peptides (AVPs) are a specific type of AMPs that have activity against certain viruses. Unlike the research focused on the development of tools and methods for the prediction of antimicrobial peptides, those related to the prediction of AVPs are still scarce. Given the significance of AVPs as potential pharmaceutical options for human and animal health and the ongoing AI revolution, we have reviewed and summarized the current machine learning and deep learning-based tools and methods available for predicting these types of peptides.

QSAR and molecular docking studies on designing potent inhibitors of SARS-CoVs main protease

Background: Severe acute respiratory syndrome coronavirus (SARS-CoVs) have emerged as a global health threat, which had caused a high rate of mortality. There is an urgent need to find effective drugs against these viruses. Objective: This study aims to predict the activity of unsymmetrical aromatic disulfides by constructing a QSAR model, and to design new compounds according to the structural and physicochemical attributes responsible for higher activity towards SARS-CoVs main protease. Methods: All molecules were constructed in ChemOffice software and molecular descriptors were calculated by CODESSA software. A regression-based linear heuristic method was established by changing descriptors datasets and calculating predicted IC₅₀ values of compounds. Then, some new compounds were designed according to molecular descriptors from the heuristic method model. The compounds with predicted values smaller than a set point were constantly screened out. Finally, the properties analysis and molecular docking were conducted to further understand the structure-activity relationships of these finalized compounds. Results: The heuristic method explored the various descriptors responsible for bioactivity and gained the best linear model with R² 0.87. The success of the model fully passed the testing set validation, proving that the model has both high statistical significance and excellent predictive ability. A total of 5 compounds with ideal predicted IC₅₀ were found from the 96 newly designed derivatives and their properties analyze was carried out. Molecular docking experiments were conducted for the optimal compound 31a, which has the best compound activity with good target protein binding capability. Conclusion: The heuristic method was quite reliable for predicting IC₅₀ values of unsymmetrical aromatic disulfides. The present research provides meaningful guidance for further exploration of the highly active inhibitors for SARS-CoVs.

[The prediction of SARS-CoV-2 main protease inhibition with filtering by position of ligand]

The paper analyzes a set of equations that adequately predict the IC50 value for SARS-CoV-2 main protease inhibitors. The training set was obtained using filtering by criteria independent of prediction of target value. It included 76 compounds, and the test set included nine compounds. We used the values of energy contributions obtained in the calculation of the change of the free energy of complex by MMGBSA method and a number of characteristics of the physical and chemical properties of the inhibitors as independent variables. It is sufficient to use only seven independent variables without loss of prediction quality (Q² = 0.79; R²prediction = 0.89). The maximum error in this case does not exceed 0.92 lg(IC50) units with a full range of observed values from 1.26 to 4.95.

Multi-Condition QSAR Model for the Virtual Design of Chemicals with Dual Pan-Antiviral and Anti-Cytokine Storm Profiles

Respiratory viruses are infectious agents, which can cause pandemics. Although nowadays the danger associated with respiratory viruses continues to be evidenced by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the virus responsible for the current COVID-19 pandemic, other viruses such as SARS-CoV-1, the influenza A and B viruses (IAV and IBV, respectively), and the respiratory syncytial virus (RSV) can lead to globally spread viral diseases. Also, from a biological point of view, most of these viruses can cause an organ-damaging hyperinflammatory response known as the cytokine storm (CS). Computational approaches constitute an essential component of modern drug development campaigns, and therefore, they have the potential to accelerate the discovery of chemicals able to simultaneously inhibit multiple molecular and nonmolecular targets. We report here the first multicondition model based on quantitative structure-activity relationships and an artificial neural network (mtc-QSAR-ANN) for the virtual design and prediction of molecules with dual pan-antiviral and anti-CS profiles. Our mtc-QSAR-ANN model exhibited an accuracy higher than 80%. By interpreting the different descriptors present in the mtc-QSAR-ANN model, we could retrieve several molecular fragments whose assembly led to new molecules with drug-like properties and predicted pan-antiviral and anti-CS activities.

A Review on Herbal Secondary Metabolites Against COVID-19 Focusing on the Genetic Variants of SARS-CoV-2

Context: An outbreak of the new coronavirus disease 2019 (COVID-19) was reported in Wuhan, China, in December 2019, subsequently affecting countries worldwide and causing a pandemic. Although several vaccines, such as mRNA vaccines, inactivated vaccines, and adenovirus vaccines, have been licensed in several countries, the danger of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants persists. To date, Alpha (B.1.1.7), Beta (B.1.351, B.1.351.2, B.1.351.3), Delta (B.1.617.2, AY.1, AY.2, AY. 3), Gamma (P.1, P.1.1, P.1.2), and Iota (B.1 .526) circulating in the United States, Kappa (B.1.617.1) in India, Lambda (C.37) in Peru and Mu (B.1.621) in Colombia are considered the variants of concern and interest. Evidence Acquisition: Data were collected through the end of August 2021 by searching PubMed, Scopus, and Google Scholar databases. There were findings from in silico, in vitro cell-based, and non-cell-based investigations. Results: The potential and safety profile of herbal medicines need clarification to scientifically support future recommendations regarding the benefits and risks of their use. Conclusions: Current research results on natural products against SARS-CoV-2 and variants are discussed, and their specific molecular targets and possible mechanisms of action are summarized.

Synthesis and Biological Activity of Unsymmetrical Monoterpenylhetaryl Disulfides

New unsymmetrical monoterpenylhetaryl disulfides based on heterocyclic disulfides and monoterpene thiols were synthesized for the first time in 48–88% yields. Hydrolysis of disulfides with fragments of methyl esters of 2-mercaptonicotinic acid was carried out in 73–95% yields. The obtained compounds were evaluated for antioxidant, antibacterial, antifungal activity, cytotoxicity and mutagenicity.

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.

SMILES-based 2D-QSAR and similarity search for identification of potential new scaffolds for development of SARS-CoV-2 MPRO inhibitors

COVID-19, whose etiological agent is the SARS-CoV-2 virus, has caused over 537.5 million cases and killed over 6.3 million people since its discovery in 2019. Despite the recent development of the first drugs indicated for treating people already infected, the great need to develop new anti-SARS-CoV-2 drugs still exists, mainly due to the possible emergence of new variants of this virus and resistant strains of the current variants. Thus, this work presents the results of QSAR and similarity search studies based only on 2D structures from a set of 32 bicycloproline derivatives, aiming to quickly, reproducibly, and reliably identify potentially useful compounds as scaffolds of new major protease inhibitors (M^pro) of the virus. The obtained QSAR model is based only on topological molecular descriptors. The model has good internal and external statistics, is robust, and does not present a chance correlation. This model was used as one of the tools to support the virtual screening stage carried out in the SwissADME web tool. Five molecules, from an initial set of 2695 molecules, proved to be the most promising, as they were within the model’s applicability domain and linearity range, with low potential to cause carcinogenic, teratogenic, and reproductive toxicity effects and promising pharmacokinetic properties. These five compounds were then selected as the most competent to generate, in future studies, new anti-SARS-CoV-2 agents with drug-likeness properties suitable for use in therapy.

Suggestion of active 3-chymotrypsin like protease (3CLPro) inhibitors as potential anti-SARS-CoV-2 agents using predictive QSAR model based on the combination of ALASSO with an ANN model

The novel severe acute respiratory syndrome coronavirus (SARS CoV-2) was introduced as an epidemic in 2019 and had millions of deaths worldwide. Given the importance of this disease, the recommendation and design of new active compounds are crucial. 3-chymotrypsin-like protease (3 CL^pro) inhibitors have been identified as potent compounds for treating SARS-CoV-2 disease. So, the design of new 3 CL^pro inhibitors was proposed using a quantitative structure-activity relationship (QSAR) study. In this context, a powerful adaptive least absolute shrinkage and selection operator (ALASSO) penalized variable selection method with inherent advantages coupled with a nonlinear artificial neural network (ANN) modelling method were used to provide a QSAR model with high interpretability and predictability. After evaluating the accuracy and validity of the developed ALASSO-ANN model, new compounds were proposed using effective descriptors, and the biological activity of the new compounds was predicted. Ligand-receptor (LR) interactions were also performed to confirm the interaction strength of the compounds using molecular docking (MD) study. The pharmacokinetics properties and calculated Lipinski's rule of five were applied to all proposed compounds. Due to the ease of synthesis of these suggested new compounds, it is expected that they have acceptable pharmacological properties.

Semi-correlations as a tool to model for skin sensitization

Semi-correlation specifically assesses the correlation between a binary variable and a continuous variable. Semi-correlations were applied to develop binary models for various endpoints. We applied the semi-correlation to develop models of two kinds of skin sensitization one related to animals (local lymph node assay LLNA) and one to human beings (direct peptide reactivity assay DPRA and/or human cell line activation test h-CLAT). The models refer to binary classification for a two-level strategy: the first level (analysis of all compounds) is used in the format "sensitizer or non-sensitizer", and the second level (only sensitizers) is a further classification in the format "strong or weak sensitizer". The ranges of statistical characteristics of the models depend on the endpoint, LLNA or DPRA/h-CLAT: for the first level, sensitivity: 0.69-0.88, specificity: 0.75-0.89, accuracy: 0.77-0.87, Matthew's correlation coefficient (MCC): 0.54-0.57 and for the second level, sensitivity: 0.70-1.0, specificity: 0.78-0.83, accuracy: 0.77-0.87, MCC: 0.54-0.76. Thus, the described approach can be applied to building up models of the skin sensitization potency.

First structure-activity relationship analysis of SARS-CoV-2 virus main protease (Mpro) inhibitors: an endeavor on COVID-19 drug discovery

Main protease (Mpro) of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) intervenes in the replication and transcription processes of the virus. Hence, it is a lucrative target for anti-viral drug development. In this study, molecular modeling analyses were performed on the structure activity data of recently reported diverse SARS-CoV-2 Mpro inhibitors to understand the structural requirements for higher inhibitory activity. The classification-based quantitative structure-activity relationship (QSAR) models were generated between SARS-CoV-2 Mpro inhibitory activities and different descriptors. Identification of structural fingerprints to increase or decrease in the inhibitory activity was mapped for possible inclusion/exclusion of these fingerprints in the lead optimization process. Challenges in ADME properties of protease inhibitors were also discussed to overcome the problems of oral bioavailability. Further, depending on the modeling results, we have proposed novel as well as potent SARS-CoV-2 Mpro inhibitors.

The index of ideality of correlation improves the predictive potential of models of the antioxidant activity of tripeptides from frog skin (Litoria rubella)

It is usually held that good-quality models for the biological activity of peptides must take into account their 3D architecture and descriptors of quantum mechanics. However, the present study shows that it is possible to build up models without these complex calculations. The structure of tripeptides represented by sequences of one-symbol abbreviations of the corresponding amino acids serves to build up quantitative structure-activity relationships for the antioxidant activity of tripeptides from frog skin. The statistical quality of the best model for the validation set is n = 27, r² = 0.93, RMSE = 0.15.

The system of self-consistent semi-correlations as one of the tools of cheminformatics for designing antiviral drugs

The development of antiviral agents against SARS-CoV-2 is necessary.

Application 2D Descriptors and Artificial Neural Networks for Beta-Glucosidase Inhibitors Screening

Beta-glucosidase inhibitors play important medical and biological roles. In this study, simple two-variable artificial neural network (ANN) classification models were developed for beta-glucosidase inhibitors screening. All bioassay data were obtained from the ChEMBL database. The classifiers were generated using 2D molecular descriptors and the data miner tool available in the STATISTICA package (STATISTICA Automated Neural Networks, SANN). In order to evaluate the models' accuracy and select the best classifiers among automatically generated SANNs, the Matthews correlation coefficient (MCC) was used. The application of the combination of maxHBint3 and SpMax8_Bhs descriptors leads to the highest predicting abilities of SANNs, as evidenced by the averaged test set prediction results (MCC = 0.748) calculated for ten different dataset splits. Additionally, the models were analyzed employing receiver operating characteristics (ROC) and cumulative gain charts. The thirteen final classifiers obtained as a result of the model development procedure were applied for a natural compounds collection available in the BIOFACQUIM database. As a result of this beta-glucosidase inhibitors screening, eight compounds were univocally classified as active by all SANNs.