Discovery of anti-MERS-CoV small covalent inhibitors through pharmacophore modeling, covalent docking and molecular dynamics simulation

Revealing innovative JAK1 and JAK3 inhibitors: a comprehensive study utilizing QSAR, 3D-Pharmacophore screening, molecular docking, molecular dynamics, and MM/GBSA analyses

The heterocycle compounds, with their diverse functionalities, are particularly effective in inhibiting Janus kinases (JAKs). Therefore, it is crucial to identify the correlation between their complex structures and biological activities for the development of new drugs for the treatment of rheumatoid arthritis (RA) and cancer. In this study, a diverse set of 28 heterocyclic compounds selective for JAK1 and JAK3 was employed to construct quantitative structure-activity relationship (QSAR) models using multiple linear regression (MLR). Artificial neural network (ANN) models were employed in the development of QSAR models. The robustness and stability of the models were assessed through internal and external methodologies, including the domain of applicability (DoA). The molecular descriptors incorporated into the model exhibited a satisfactory correlation with the receptor-ligand complex structures of JAKs observed in X-ray crystallography, making the model interpretable and predictive. Furthermore, pharmacophore models ADRRR and ADHRR were designed for each JAK1 and JAK3, proving effective in discriminating between active compounds and decoys. Both models demonstrated good performance in identifying new compounds, with an ROC of 0.83 for the ADRRR model and an ROC of 0.75 for the ADHRR model. Using a pharmacophore model, the most promising compounds were selected based on their strong affinity compared to the most active compounds in the studied series each JAK1 and JAK3. Notably, the pharmacokinetic, physicochemical properties, and biological activities of the selected compounds (As compounds ZINC79189223 and ZINC66252348) were found to be consistent with their therapeutic effects in RA, owing to their non-toxic, cholinergic nature, absence of P-glycoprotein, high gastrointestinal absorption, and ability to penetrate the blood-brain barrier. Furthermore, ADMET properties were assessed, and molecular dynamics and MM/GBSA analysis revealed stability in these molecules.

Discovery of pyrimidoindol and benzylpyrrolyl inhibitors targeting SARS-CoV-2 main protease (Mpro) through pharmacophore modelling, covalent docking, and biological evaluation

The main protease (Mpro) enzyme has an imperative function in disease progression and the life cycle of the SARS-CoV-2 virus. Although the orally active drug nirmatrelvir (co-administered with ritonavir as paxlovid) has been approved for emergency use as the frontline antiviral agent, there are a number of limitations that necessitate the discovery of new drug scaffolds, such as poor pharmacokinetics and susceptibility to proteolytic degradation due to its peptidomimetic nature. This study utilized a novel virtual screening workflow that combines pharmacophore modelling, multiple-receptor covalent docking, and biological evaluation in order to find new Mpro inhibitors. After filtering and analysing ∼66,000 ligands from three different electrophilic libraries, 29 compounds were shortlisted for experimental testing, and two of them exhibited ≥20% inhibition at 100 μM. Our top candidate, GF04, is a benzylpyrrolyl compound that exhibited the highest inhibition activity of 38.3%, with a relatively small size (<350 Da) and leadlike character. Interestingly, our approach also identified another hit, DR07, a pyrimidoindol with a non-peptide character, and a molecular weight of 438.9 Da, reporting an inhibition of 26.3%. The established approach detailed in this study, in conjunction with the discovered inhibitors, has the capacity to yield novel perspectives for devising covalent inhibitors targeting the COVID-19 Mpro enzyme and other comparable targets.

An in silico approach to develop potential therapies against Middle East Respiratory Syndrome Coronavirus (MERS-CoV)

A deadly respiratory disease Middle East Respiratory Syndrome (MERS) is caused by a perilous virus known as MERS-CoV, which has a severe impact on human health. Currently, there is no approved vaccine, prophylaxis, or antiviral therapeutics for preventing MERS-CoV infection. Due to its inexorable and integral role in the maturation and replication of the MERS-CoV virus, the 3C-like protease is unavoidly a viable therapeutic target. In this study, 2369 phytoconstituents were enlisted from Japanese medicinal plants, and these compounds were screened against 3C-like protease to identify feasible inhibitors. The best three compounds were identified as Kihadanin B, Robustaflavone, and 3-beta-O- (trans-p-Coumaroyl) maslinic acid, with binding energies of -9.8, -9.4, and -9.2 kcal/mol, respectively. The top three potential candidates interacted with several active site residues in the targeted protein, including Cys145, Met168, Glu169, Ala171, and Gln192. The best three compounds were assessed by in silico technique to determine their drug-likeness properties, and they exhibited the least harmful features and the greatest drug-like qualities. Various descriptors, such as solvent-accessible surface area, root-mean-square fluctuation, root-mean-square deviation, hydrogen bond, and radius of gyration, validated the stability and firmness of the protein-ligand complexes throughout the 100ns molecular dynamics simulation. Moreover, the top three compounds exhibited better binding energy along with better stability and firmness than the inhibitor (Nafamostat), which was further confirmed by the binding free energy calculation. Therefore, this computational investigation could aid in the development of efficient therapeutics for life-threatening MERS-CoV infections.

Integrated virtual screening and molecular dynamics simulation approaches revealed potential natural inhibitors for DNMT1 as therapeutic solution for triple negative breast cancer

Triple negative breast cancers (TNBC) are clinically heterogeneous but mostly aggressive malignancies devoid of expression of the estrogen, progesterone, and HER2 (ERBB2 or NEU) receptors. It accounts for 15-20% of all cases. Altered epigenetic regulation including DNA hypermethylation by DNA methyltransferase 1 (DNMT1) has been implicated as one of the causes of TNBC tumorigenesis. The antitumor effect of DNMT1 has also been explored in TNBC that currently lacks targeted therapies. However, the actual treatment for TNBC is yet to be discovered. This study is attributed to the identification of novel drug targets against TNBC. A comprehensive docking and simulation analysis was performed to optimize promising new compounds by estimating their binding affinity to the target protein. Molecular dynamics simulation of 500 ns well complemented the binding affinity of the compound and revealed strong stability of predicted compounds at the docked site. Calculation of binding free energies using MMPBSA and MMGBSA validated the strong binding affinity between compound and binding pockets of DNMT1. In a nutshell, our study uncovered that Beta-Mangostin, Gancaonin Z, 5-hydroxysophoranone, Sophoraflavanone L, and Dorsmanin H showed maximum binding affinity with the active sites of DNMT1. Furthermore, all of these compounds depict maximum drug-like properties. Therefore, the proposed compounds can be a potential candidate for patients with TNBC, but, experimental validation is needed to ensure their safety.Communicated by Ramaswamy H. Sarma.

Identification of novel inhibitors against VP40 protein of Marburg virus by integrating molecular modeling and dynamics approaches

Marburg virus (MV) is a highly etiological agent of haemorrhagic fever in humans and has spread across the world. Its outbreaks caused a 23-90% human death rate. However, there are currently no authorized preventive or curative measures yet. VP40 is the MV matrix protein, which builds protein shell underneath the viral envelope and confers hallmark filamentous. VP40 alone is able to induce assembly and budding of filamentous virus-like particles (VLPs), which resemble authentic virions. As a result, this research is credited with clarifying the function of VP40 and leading to the discovery of new therapeutic targets effective in combating MV disease (MVD). Virtual screening, molecular docking and molecular dynamics (MD) simulation were used to find the putative active chemicals based on a 3D pharmacophore model of the protein's active site cavity. Initially, andrographidine-C, a potent inhibitor was selected for the development of the pharmacophore model. Later, a library of 30,000 compounds along with the andrographidine-C was docked against VP40 protein. Three best hits including avanafil, diuvaretin and macrourone were subjected to further MD simulation analysis, as these compounds had better binding affinities as compared to andrographidine-C. Furthermore, throughout the 100 ns simulations, the back bone of VP40 protein in presence of avanafil, diuvaretin and macrourone remained stable which was further validated by MM-PBSA analysis. Additionally, all of these compounds depict maximum drug-like properties. The predicted drugs based on the ligand, avanafil, diuvaretin and macrourone could be exploited and developed as an alternative or complementary therapy for the treatment of MVD.Communicated by Ramaswamy H. Sarma.

Structural Dynamics of P-Rex1 Complexed with Natural Leads Establishes the Protein as an Attractive Target for Therapeutics to Suppress Cancer Metastasis

Phosphatidylinositol 3,4,5-trisphosphate- (PIP3-) dependent Rac exchanger 1 (P-Rex1) functions as Rho guanine nucleotide exchange factor and is activated by synergistic activity of Gβγ and PIP3 of the heterotrimeric G protein. P-Rex1 activates Rac GTPases for regulating cell invasion and migration and promotes metastasis in several human cancers including breast, prostate, and skin cancer. The protein is a promising therapeutic target because of its multifunction roles in human cancers. Herein, the present study attempts to identify selective P-Rex1 natural inhibitors by targeting PIP3-binding pocket using large-size multiple natural molecule libraries. Each library was filtered subsequently in FAF-Drugs4 based on Lipinski's rule of five (RO5), toxicity, and filter pan assay interference compounds (PAINS). The output hits were virtually screened at the PIP3-binding pocket through PyRx AutoDock Vina and cross-checked by GOLD. The best binders at the PIP3-binding pocket were prioritized using a comparative analysis of the docking scores. Top-ranked two compounds with high GOLD fitness score (>80) and lowest AutoDock binding energy (< -12.7 kcal/mol) were complexed and deciphered for molecular dynamics along with control-P-Rex1 complex to validate compound binding conformation and disclosed binding interaction pattern. Both the systems were seen in good equilibrium, and along the simulation time, the compounds are in strong contact with the P-Rex1 PIP3-binding site. Hydrogen bonding analysis towards simulation end identified the formation of 16 and 22 short- and long-distance hydrogen bonds with different percent of occupancy to the PIP3 residues for compound I and compound 2, respectively. Radial distribution function (RDF) analysis of the key hydrogen bonds between the compound and the PIP3 residues demonstrated a strong affinity of the compounds to the mentioned PIP3 pocket. Additionally, MMGB/PBSA energies were performed that confirmed the dominance of Van der Waals energy in complex formation along with favorable contribution from hydrogen bonding. These findings were also cross-validated by a more robust WaterSwap binding energy predictor, and the results are in good agreement with a strong binding affinity of the compounds for the protein. Lastly, the key contribution of residues in interaction with the compounds was understood by binding free energy decomposition and alanine scanning methods. In short, the results of this study suggest that P-Rex1 is a good druggable target to suppress cancer metastasis; therefore, the screened druglike molecules of this study need in vitro and in vivo anti-P-Rex1 validation and may serve as potent leads to fight cancer.

A Multifaceted Computational Approach to Understanding the MERS-CoV Main Protease and Brown Algae Compounds' Interaction

Middle East Respiratory Syndrome (MERS) is a viral respiratory disease caused b a special type of coronavirus called MERS-CoV. In the search for effective substances against the MERS-CoV main protease, we looked into compounds from brown algae, known for their medicinal benefits. From a set of 1212 such compounds, our computer-based screening highlighted four-CMNPD27819, CMNPD1843, CMNPD4184, and CMNPD3156. These showed good potential in how they might attach to the MERS-CoV protease, comparable to a known inhibitor. We confirmed these results with multiple computer tests. Studies on the dynamics and steadiness of these compounds with the MERS-CoV protease were performed using molecular dynamics (MD) simulations. Metrics like RMSD and RMSF showed their stability. We also studied how these compounds and the protease interact in detail. An analysis technique, PCA, showed changes in atomic positions over time. Overall, our computer studies suggest brown algae compounds could be valuable in fighting MERS. However, experimental validation is needed to prove their real-world effectiveness.

Docking covalent targets for drug discovery: stimulating the computer-aided drug design community of possible pitfalls and erroneous practices

The continuous approval of covalent drugs in recent years for the treatment of diseases has led to an increased search for covalent agents by medicinal chemists and computational scientists worldwide. In the computational parlance, molecular docking which is a popular tool to investigate the interaction of a ligand and a protein target, does not account for the formation of covalent bond, and the increasing application of these conventional programs to covalent targets in early drug discovery practice is a matter of utmost concern. Thus, in this comprehensive review, we sought to educate the docking community about the realization of covalent docking and the existence of suitable programs to make their future virtual-screening events on covalent targets worthwhile and scientifically rational. More interestingly, we went beyond the classical description of the functionality of covalent-docking programs down to selecting the 'best' program to consult with during a virtual-screening campaign based on receptor class and covalent warhead chemistry. In addition, we made a highlight on how covalent docking could be achieved using random conventional docking software. And lastly, we raised an alert on the growing erroneous molecular docking practices with covalent targets. Our aim is to guide scientists in the rational docking pursuit when dealing with covalent targets, as this will reduce false-positive results and also increase the reliability of their work for translational research.

Integrated molecular modeling and dynamics approaches revealed potential natural inhibitors of NF-κB transcription factor as breast cancer therapeutics

Breast cancer is a silent killer malady among women and a serious economic burden in health care management. A case of breast cancer is diagnosed among women every 19 s, and every 74 s, a woman dies of breast cancer somewhere in the world. Despite the pop-up of progressive research, advanced treatment approaches, and preventive measures, breast cancer remains amplifying ailment. The nuclear factor kappa B (NF-κB) is a key transcription factor that links inflammation with cancer and is demonstrated as being involved in the tumorigenesis of breast cancer. The NF-κB transcription factor family in mammals consists of five proteins; c-Rel, RelA(p65), RelB, NF-κB1(p50), and NF-κB2(p52). The antitumor effect of NF-κB has also been explored in breast cancer, however, the actual treatment for breast cancer is yet to be discovered. This study is attributed to the identification of novel drug targets against breast cancer by targeting c-Rel, RelA(p65), RelB, NF-κB1(p50), and NF-κB2(p52) proteins. To identify the putative active compounds, a structure-based 3D pharmacophore model to the protein active site cavity was generated followed by virtual screening, molecular docking, and molecular dynamics (MD) simulation. Initially, a library of 45000 compounds were docked against the target protein and five compounds namely Z56811101, Z653426226, Z1097341967, Z92743432, and Z464101066 were selected for further analysis. The relative binding affinity of Z56811101, Z653426226, Z1097341967, Z92743432, and Z464101066 with NF-κB1 (p50), NF-κB2 (p52), RelA (p65), RelB, and c-Rel proteins were -6.8, -8, -7.0, -6.9, and -7.2 kcal/mol, respectively which remained stable throughout the simulations of 200 ns. Furthermore, all of these compounds depict maximum drug-like properties. Therefore, the proposed compounds can be a potential candidate for patients with breast cancer, but, experimental validation is needed to ensure their safety.Communicated by Ramaswamy H. Sarma.

Identification of molecular mechanisms underlying the therapeutic effects of Celosia Cristata on immunoglobulin nephropathy

Immunoglobulin A (IgA) nephropathy also known as Berger's disease, is a silent monster and perhaps the most prevalent glomerulonephritis that often accounts for end-stage kidney failure, thereby signifying a growing public health problem worldwide. The limited amount of available data and a broad spectrum of dysregulated physiological processes of IgAN make it a challenging task and a disproportionate economic load on the community health sector. Celosia cristata is an Amaranthaceous plant with attractive colorful inflorescences that are used in various regions of earth for the treatment of numerous ailments. A list of studies evidences the therapeutic efficacy of C. cristata against complicated disorders, but the precise molecular mechanism is yet to be discovered. This study is attributed to the identification of bioactive compounds, pathways, and target genes for the better treatment of IgAN. In the current analysis, compound-target genes-pathway networks were explored which uncovered that isorhamnetin, stigmasterol, luteolin, amaranthin, and β-sitosterol may serve as a magic bullet against IgAN by influencing the targets genes involved in the disease pathogenesis. Later, the expression of hub genes was then further analyzed using a microarray dataset (GSE93798). Through expression analysis, it is worth noting that FOS, JUN, and EGFR were considerably upregulated, and at the same moment, AKT1 was considerably downregulated in IgAN patients. Lastly, docking analysis further strengthened the current findings by validating the effective activity of the active ingredients against putative target genes. In summary, we propose that five key compounds including, isorhamnetin, stigmasterol, luteolin, amaranthin, and β-sitosterol, aid in the regulation of JUN, FOS, AKT1, and EGFR, which may serve as a promising and enthralling therapeutic option for IgAN. The overall integration of network pharmacology with molecular docking unveiled the multi-target pharmacological mechanisms of C. cristata against IgAN. This study provides convincing evidence that C. cristata might partially alleviate the IgAN and ultimately lays a foundation for further experimental research on the anti-IgAN activity of C. cristata.

The Discovery of Small Allosteric and Active Site Inhibitors of the SARS-CoV-2 Main Protease via Structure-Based Virtual Screening and Biological Evaluation

The main protease enzyme (M^pro) of SARS-CoV-2 is one of the most promising targets for COVID-19 treatment. Accordingly, in this work, a structure-based virtual screening of 3.8 million ligand libraries was carried out. After rigorous filtering, docking, and post screening assessments, 78 compounds were selected for biological evaluation, 3 of which showed promising inhibition of the M^pro enzyme. The obtained hits (CB03, GR04, and GR20) had reasonable potencies with K_i values in the medium to high micromolar range. Interestingly, while our most potent hit, GR20, was suggested to act via a reversible covalent mechanism, GR04 was confirmed as a noncompetitive inhibitor that seems to be one of a kind when compared to the other allosteric inhibitors discovered so far. Moreover, all three compounds have small sizes (~300 Da) with interesting fittings in their relevant binding sites, and they possess lead-like characteristics that can introduce them as very attractive candidates for the future development of COVID-19 treatments.

Computational exploration of maternal embryonic leucine zipper kinase (MELK) as a cancer drug target

Maternal embryonic leucine zipper kinase (MELK) is of vital importance due to its significant role in cancer development and its association with poor prognosis in different cancers. Here, we employed several computer aided drug design approaches to shortlist potential binding molecules of MELK. For virtual screening, asinex oncology library (containing 6334 drugs) and comprehensive marine natural products database (containing approximately 32,000 drugs) were used. The study identified two drug molecules: Top-2 and Top-3 as high affinity binding MELK molecules compared to the control co-crystalized Top-1 inhibitor. Both the shortlisted compounds and the control showed high stable binding free energy and high GOLD score. The compounds and control also reported stable dynamics with root mean square deviations (RMSD) value ∼ 2 Å in 500 ns. Similarly, the MELK active site residues were observed in good stability with the compounds. Further, it was noticed the compounds/control formed multiple hydrogen bonds with the MELK active pocket residues which is the main reason of high intermolecular stability. Atomic level binding free energies determined van der Waals and electrostatic energies to play vital role in stable complex formation. From drug likeness and pharmacokinetics perspective, the compounds are ideal molecules for further investigation. Overall, the results are promising and might be tested in in vivo and in vitro studies against MELK.

Potential Therapeutic Candidates against Chlamydia pneumonia Discovered and Developed In Silico Using Core Proteomics and Molecular Docking and Simulation-Based Approaches

Chlamydia pneumonia, a species of the family Chlamydiacea, is a leading cause of pneumonia. Failure to eradicate C. pneumoniae can lead to chronic infection, which is why it is also considered responsible for chronic inflammatory disorders such as asthma, arthritis, etc. There is an urgent need to tackle the major concerns arising due to persistent infections caused by C. pneumoniae as no FDA-approved drug is available against this chronic infection. In the present study, an approach named subtractive proteomics was employed to the core proteomes of five strains of C. pneumonia using various bioinformatic tools, servers, and software. However, 958 non-redundant proteins were predicted from the 4754 core proteins of the core proteome. BLASTp was used to analyze the non-redundant genes against the proteome of humans, and the number of potential genes was reduced to 681. Furthermore, based on subcellular localization prediction, 313 proteins with cytoplasmic localization were selected for metabolic pathway analysis. Upon subsequent analysis, only three cytoplasmic proteins, namely 30S ribosomal protein S4, 4-hydroxybenzoate decarboxylase subunit C, and oligopeptide binding protein, were identified, which have the potential to be novel drug target candidates. The Swiss Model server was used to predict the target proteins’ three-dimensional (3D) structure. The molecular docking technique was employed using MOE software for the virtual screening of a library of 15,000 phytochemicals against the interacting residues of the target proteins. Molecular docking experiments were also evaluated using molecular dynamics simulations and the widely used MM-GBSA and MM-PBSA binding free energy techniques. The findings revealed a promising candidate as a novel target against C. pneumonia infections.

Structural Elucidation of Rift Valley Fever Virus L Protein towards the Discovery of Its Potential Inhibitors

Rift valley fever virus (RVFV) is the causative agent of a viral zoonosis that causes a significant clinical burden in domestic and wild ruminants. Major outbreaks of the virus occur in livestock, and contaminated animal products or arthropod vectors can transmit the virus to humans. The viral RNA-dependent RNA polymerase (RdRp; L protein) of the RVFV is responsible for viral replication and is thus an appealing drug target because no effective and specific vaccine against this virus is available. The current study reported the structural elucidation of the RVFV-L protein by in-depth homology modeling since no crystal structure is available yet. The inhibitory binding modes of known potent L protein inhibitors were analyzed. Based on the results, further molecular docking-based virtual screening of Selleckchem Nucleoside Analogue Library (156 compounds) was performed to find potential new inhibitors against the RVFV L protein. ADME (Absorption, Distribution, Metabolism, and Excretion) and toxicity analysis of these compounds was also performed. Besides, the binding mechanism and stability of identified compounds were confirmed by a 50 ns molecular dynamic (MD) simulation followed by MM/PBSA binding free energy calculations. Homology modeling determined a stable multi-domain structure of L protein. An analysis of known L protein inhibitors, including Monensin, Mycophenolic acid, and Ribavirin, provide insights into the binding mechanism and reveals key residues of the L protein binding pocket. The screening results revealed that the top three compounds, A-317491, Khasianine, and VER155008, exhibited a high affinity at the L protein binding pocket. ADME analysis revealed good pharmacodynamics and pharmacokinetic profiles of these compounds. Furthermore, MD simulation and binding free energy analysis endorsed the binding stability of potential compounds with L protein. In a nutshell, the present study determined potential compounds that may aid in the rational design of novel inhibitors of the RVFV L protein as anti-RVFV drugs.

Integrated computational approaches to screen gene expression data to determine key genes and therapeutic targets for type-2 diabetes mellitus

There is a rapid rise in cases of Type-2-diabetes mellitus (T2DM) globally, irrespective of the geography, ethnicity or any other variable factors. The molecular mechanisms that could cause the condition of T2DM need to be more thoroughly analysed to understand the clinical manifestations and to derive better therapeutic regimes. Tools in bioinformatics are used to trace out key gene elements and to identify the key causative gene elements and their possible therapeutic agents. Microarray datasets were retrieved from the Gene expression omnibus database and studied using R to derive different expressed gene (DEG) elements. With the comparison of the expressed genes with disease specific genes in DisGeNET, the final annotated genes were taken for analysis. Gene Ontology studies, Protein-protein interaction (PPI), Co-expression analysis, Gene-drug interactions were performed to scale down the hub genes and to identify the novelty across the genes analysed so far. In vivo and invitro analysis of key genes and the trace of interaction pathway is crucial to better understand the unique outcomes from the novel genes, forming the basis to understand the pathway that ends up causing T2DM. Afterwards, docking was executed enabling recognition of interacting residues involved in inhibition. The complex CCL5-265 and CD8A-40585 thus docked showed best results as is evident from its PCA analysis and MMGBSA calculation. There is now scope for deriving candidate drugs that could possibly detect personalized therapies for T2DM.

Integrated computer-aided drug design and biophysical simulation approaches to determine natural anti-bacterial compounds for Acinetobacter baumannii

Acinetobacter baumannii is a nosocomial bacterial pathogen and is responsible for a wide range of diseases including pneumonia, necrotizing fasciitis, meningitis, and sepsis. The enzyme 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase (encoded by aroA gene) in ESKAPE pathogens catalyzes the sixth step of shikimate pathway. The shikimate pathway is an attractive drug targets pathway as it is present in bacteria but absent in humans. As EPSP is essential for the A. baumannii growth and needed during the infection process, therefore it was used as a drug target herein for high-throughput screening of a comprehensive marine natural products database (CMNPD). The objective was to identify natural molecules that fit best at the substrate binding pocket of the enzyme and interact with functionally critical residues. Comparative assessment of the docking scores allowed selection of three compounds namely CMNPD31561, CMNPD28986, and CMNPD28985 as best binding molecules. The molecules established a balanced network of hydrophobic and hydrophilic interactions, and the binding pose remained in equilibrium throughout the length of molecular simulation time. Radial distribution function (RDF) analysis projected key residues from enzyme active pocket which actively engaged the inhibitors. Further validation is performed through binding free energies estimation that affirms very low delta energy of <-22 kcal/mol in MM-GBSA method and <-12 kcal/mol in MM-PBSA method. Lastly, the most important active site residues were mutated and their ligand binding potential was re-investigated. The molecules also possess good druglike properties and better pharmacokinetics. Together, these findings suggest the potential biological potency of the leads and thus can be used by experimentalists in vivo and in vitro studies.

Peptide candidates for the development of therapeutics and vaccines against β-coronavirus infection

Betacoronaviruses (β-CoVs) have caused major viral outbreaks in the last two decades in the world. The mutation and recombination abilities in β-CoVs resulted in zoonotic diseases in humans. Proteins responsible for viral attachment and replication are highly conserved in β-CoVs. These conserved proteins have been extensively studied as targets for preventing infection and the spread of β-CoVs. Peptides are among the most promising candidates for developing vaccines and therapeutics against viral pathogens. The immunostimulatory and viral inhibitory potential of natural and synthetic peptides has been extensively studied since the SARS-CoV outbreak. Food-derived peptides demonstrating high antiviral activity can be used to develop effective therapeutics against β-CoVs. Specificity, tolerability, and customizability of peptides can be explored to develop potent drugs against β-CoVs. However, the proteolytic susceptibility and low bioavailability of peptides pose challenges for the development of therapeutics. This review illustrates the potential role of peptides in eliciting an adaptive immune response and inhibiting different stages of the β-CoV life cycle. Further, the challenges and future directions associated with developing peptide-based therapeutics and vaccines against existing and future β-CoV pathogens have been discussed.

In Silico Drug Repositioning to Target the SARS-CoV-2 Main Protease as Covalent Inhibitors Employing a Combined Structure-Based Virtual Screening Strategy of Pharmacophore Models and Covalent Docking

The epidemic caused by the SARS-CoV-2 coronavirus, which has spread rapidly throughout the world, requires urgent and effective treatments considering that the appearance of viral variants limits the efficacy of vaccines. The main protease of SARS-CoV-2 (Mpro) is a highly conserved cysteine proteinase, fundamental for the replication of the coronavirus and with a specific cleavage mechanism that positions it as an attractive therapeutic target for the proposal of irreversible inhibitors. A structure-based strategy combining 3D pharmacophoric modeling, virtual screening, and covalent docking was employed to identify the interactions required for molecular recognition, as well as the spatial orientation of the electrophilic warhead, of various drugs, to achieve a covalent interaction with Cys145 of Mpro. The virtual screening on the structure-based pharmacophoric map of the SARS-CoV-2 Mpro in complex with an inhibitor N3 (reference compound) provided high efficiency by identifying 53 drugs (FDA and DrugBank databases) with probabilities of covalent binding, including N3 (Michael acceptor) and others with a variety of electrophilic warheads. Adding the energy contributions of affinity for non-covalent and covalent docking, 16 promising drugs were obtained. Our findings suggest that the FDA-approved drugs Vaborbactam, Cimetidine, Ixazomib, Scopolamine, and Bicalutamide, as well as the other investigational peptide-like drugs (DB04234, DB03456, DB07224, DB7252, and CMX-2043) are potential covalent inhibitors of SARS-CoV-2 Mpro.

Structural and Dynamic Insights into the W68L, L85P, and T87A Mutations of Mycobacterium tuberculosis Inducing Resistance to Pyrazinamide

Tuberculosis (TB), the most frequent bacterium-mediated infectious disease caused by Mycobacterium tuberculosis, has been known to infect humans since ancient times. Although TB is common worldwide, the most recent report by the WHO (World Health Organization) listed the three countries of India, China, and Russia with 27%, 14%, and 8% of the global burden of TB, respectively. It has been reported that resistance to TB drugs, particularly by the pncA gene to the pyrazinamide drug due to mutations, significantly affects the effective treatment of TB. Understanding the mechanism of drug resistance using computational methods is of great interest to design effective TB treatment, exploring the structural features with these tools. Thus, keeping in view the importance of these methods, we employed state-of-the-art computational methods to study the mechanism of resistance caused by the W68L, L85P, and T87A mutations recently reported in 2021. We employed a molecular docking approach to predict the binding conformation and studied the dynamic properties of each complex using molecular dynamics simulation approaches. Our analysis revealed that compared to the wildtype, these three mutations altered the binding pattern and reduced the binding affinity. Moreover, the structural dynamic features also showed that these mutations significantly reduced the structural stability and packing, particularly by the W68L and L85P mutations. Moreover, principal component analysis, free energy landscape, and the binding free energy results revealed variation in the protein's motion and the binding energy. The total binding free energy was for the wildtype -9.61 kcal/mol, W68L -7.57 kcal/mol, L85P -6.99 kcal/mol, and T87A -7.77 kcal/mol. Our findings can help to design a structure-based drug against the MDR (multiple drug-resistant) TB.

Identification of lactoferrin-derived peptides as potential inhibitors against the main protease of SARS-CoV-2

COVID-19 is a global health emergency that causes serious concerns. A global effort is underway to identify drugs for the treatment of COVID-19. One possible solution to the present problem is to develop drugs that can inhibit SARS-CoV-2 main protease (M^pro), a coronavirus protein that been considered as one among many drug targets. In this work, lactoferrin from Bos taurus L. was in silico hydrolyzed. The bioactivity, water solubility, and ADMET properties of the generated peptides were predicted using various online tools. The molecular interactions between M^pro and the peptides were studied using molecular docking and molecular dynamic simulation. The results demonstrated that peptide GSRY was predicted to have better physicochemical properties, and the value of '-C DOCKER interaction energy' between peptide GSRY and M^pro was 80.8505 kcal/mol. The interaction between the peptide GSRY and the native ligand N3 co-crystallized with M^pro had overlapped amino acids, i.e., HIS163, GlY143, GLU166, GLN189 and MET165. Molecular dynamic simulation revealed that M^pro/GSRY complexes were stable. Collectively, the peptide GSRY may be a potential candidate drug against M^pro of SARS-CoV-2.

Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches

The perennial emergence of SARS-CoV-2 and its new variants causing upper respiratory complexities since December 2019 has aggravated the pandemic situation around the world. SARS-CoV-2 encodes several proteins among which ORF8 is a novel factor that is unique to SARS-CoV-2 only and is reported to help the virus in disease severity and immune evasion. ORF8-IRF3 complex induces endoplasmic reticulum stress, thus helps in the evasion of immune response. Consequently, targeting the ORF8-IRF3 complex is considered as a prime target for the discovery of novel drugs against SARS-CoV-2. In this regard, computational methods are of great interest to fast track the identification and development of novel drugs. Virtual screening of South African Natural Compounds Database (SANCDB), followed by docking and molecular dynamics (MD) simulation analysis, were performed to determine novel natural compounds. Computational molecular search and rescoring of the SANCDB database followed by induced-fit docking (IFD) protocol identified Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside (SANC00850), Tribuloside (SANC01050), and Rutin (SANC00867) are the best scoring compounds. Structural-dynamic properties assessment revealed that these three compounds have stable dynamics, compactness, and a higher number of hydrogen bonds. For validation, we used MM/GBSA, in silico bioactivity estimation and dissociation constant (KD) approaches, which revealed that these compounds are the more potent inhibitors of the ORF8-IRF3 complex and would rescue the host immune system potentially. These compounds need further in vitro and in vivo validations to be used as therapeutics against SARS-CoV-2 to rescue the host immune system during COVID-19 infection.

Determination of Novel Anti-Cancer Agents by Targeting OGG1 Enzyme Using Integrated Bioinformatics Methods

The 8-oxoguanine DNA glycosylase (OGG1) enzyme is a key DNA glycosylase mediating the excision of 7,8-dihydro-8-oxoguanine (8-oxoG) from DNA molecule to the start base excision repair pathway. The OGG1 glycosylase function depletion has been seen to obstruct pathological conditions such as inflammation, A3 T-cell lymphoblastic acute leukemia growth, and neurodegenerative diseases, thus warranting OGG1 as an attractive anti-cancer enzyme. Herein, we employed several drug libraries intending to screen non-toxic inhibitory molecules against the active pocket of the enzyme that achieved stable binding mode in dynamics. Two anti-cancer compounds ([O-]C1=C(CC2=CC=CC=C2)SC(=[N+]1CC(=O)NC3=NC=C(CC4=CC=CC=C4)S3)S and CCCN(CCC)[S]-(=O)(=O)C1=CC=C(C=C1)C(=O)NNC2=NC3=CC=C(Br)C=C3C(=N2)C4=CC=CC=C4) from Selleckchem.com were identified to occupy the active pocket of OGG1 and bind with greater affinity than Control TH5487. The binding affinity of Top-1 is -11.6 kcal/mol while that of Top-2 is -10.7 kcal/mol in contrast to TH5487 Control (-9 kcal/mol). During molecular dynamic simulations versus time, the said compounds are tightly held by the enzyme with no minor structural deviations reported except flexible loops in particular those present at the N and C-terminal. Both the compounds produced extensive hydrophobic interactions with the enzyme along with stable hydrogen bonding. The docking and molecular dynamics simulations predictions were further validated by molecular mechanics with generalized Born and surface area solvation (MM/GBSA) and Poisson Boltzmann surface area (MM/PBSA), and WaterSwap binding energies that validated strong binding of the compounds to the enzyme. The MM/GBSA binding free energy for Top-1 complex is -28.10 kcal/mol, Top-2 complex is -50.14 kcal/mol) and Control is -46.91 kcal/mol while MM/PBSA value for Top-1, Top-2 and Control is -23.38 kcal/mol, -35.29 kcal/mol and -38.20 kcal/mol, respectively. Computational pharmacokinetics support good druglike candidacy of the compounds with acceptable profile of pharmacokinetics and very little toxicity. All these findings support the notion that the compounds can be used in experiments to test their anti-cancer activities.

Insights into the Binding of Receptor-Binding Domain (RBD) of SARS-CoV-2 Wild Type and B.1.620 Variant with hACE2 Using Molecular Docking and Simulation Approaches

Recently, a new variant, B.1620, with mutations (S477N-E484K) in the spike protein's receptor-binding domain (RBD) has been reported in Europe. In order to design therapeutic strategies suitable for B.1.620, further studies are required. A detailed investigation of the structural features and variations caused by these substitutions, that is, a molecular level investigation, is essential to uncover the role of these changes. To determine whether and how the binding affinity of ACE2-RBD is affected, we used protein-protein docking and all-atom simulation approaches. Our analysis revealed that B.1.620 binds more strongly than the wild type and alters the hydrogen bonding network. The docking score for the wild type was reported to be -122.6 +/- 0.7 kcal/mol, while for B.1.620, the docking score was -124.9 +/- 3.8 kcal/mol. A comparative binding investigation showed that the wild-type complex has 11 hydrogen bonds and one salt bridge, while the B.1.620 complex has 14 hydrogen bonds and one salt bridge, among which most of the interactions are preserved between the wild type and B.1.620. A dynamic analysis of the two complexes revealed stable dynamics, which corroborated the global stability trend, compactness, and flexibility of the three essential loops, providing a better conformational optimization opportunity and binding. Furthermore, binding free energy revealed that the wild type had a total binding energy of -51.14 kcal/mol, while for B.1.628, the total binding energy was -68.25 kcal/mol. The current findings based on protein complex modeling and bio-simulation methods revealed the atomic features of the B.1.620 variant harboring S477N and E484K mutations in the RBD and the basis for infectivity. In conclusion, the current study presents distinguishing features of B.1.620, which can be used to design structure-based drugs against the B.1.620 variant.

Allergic Diseases: A Comprehensive Review on Risk Factors, Immunological Mechanisms, Link with COVID-19, Potential Treatments, and Role of Allergen Bioinformatics

The prevalence of allergic diseases is regarded as one of the key challenges in health worldwide. Although the precise mechanisms underlying this rapid increase in prevalence are unknown, emerging evidence suggests that genetic and environmental factors play a significant role. The immune system, microbiota, viruses, and bacteria have all been linked to the onset of allergy disorders in recent years. Avoiding allergen exposure is the best treatment option; however, steroids, antihistamines, and other symptom-relieving drugs are also used. Allergen bioinformatics encompasses both computational tools/methods and allergen-related data resources for managing, archiving, and analyzing allergological data. This study highlights allergy-promoting mechanisms, algorithms, and concepts in allergen bioinformatics, as well as major areas for future research in the field of allergology.

Discovery of Novel Inhibitors From Medicinal Plants for V-Domain Ig Suppressor of T-Cell Activation

V-domain Ig suppressor of T cell activation (VISTA) is an immune checkpoint and is a type I transmembrane protein. VISTA is linked to immunotherapy resistance, and it is a potential immune therapeutic target, especially for triple-negative breast cancer. It expresses at a high concentration in regulatory T cells and myeloid-derived suppressor cells, and its functional blockade is found to delay tumor growth. A useful medicinal plant database for drug designing (MPD3), which is a collection of phytochemicals from diverse plant families, was employed in virtual screening against VISTA to prioritize natural inhibitors against VISTA. Three compounds, Paratocarpin K (PubChem ID: 14187087), 3-(1H-Indol-3-yl)-2-(trimethylazaniumyl)propanoate (PubChem ID: 3861164), and 2-[(5-Benzyl-4-ethyl-1,2,4-triazol-3-yl)sulfanylmethyl]-5-methyl-1,3,4-oxadiazole (PubChem ID: 6494266), having binding energies stronger than -6 kcal/mol were found to have two common hydrogen bond interactions with VISTA active site residues: Arg54 and Arg127. The dynamics of the compound-VISTA complexes were further explored to infer binding stability of the systems. Results revealed that the compound 14187087 and 6494266 systems are highly stable with an average RMSD of 1.31 Å. Further affirmation on the results was achieved by running MM-GBSA on the MD simulation trajectories, which re-ranked 14187087 as the top-binder with a net binding energy value of -33.33 kcal/mol. In conclusion, the present study successfully predicted natural compounds that have the potential to block the function of VISTA and therefore can be utilized further in experimental studies to validate their real anti-VISTA activity.

Molecular Insights into Binding Mode and Interactions of Structure-Based Virtually Screened Inhibitors for Pseudomonas aeruginosa Multiple Virulence Factor Regulator (MvfR)

Multi-drug resistance (MDR) bacterial pathogens pose a threat to global health and warrant the discovery of new therapeutic molecules, particularly those that can neutralize their virulence and stop the evolution of new resistant mechanisms. The superbug nosocomial pathogen, Pseudomonas aeruginosa, uses a multiple virulence factor regulator (MvfR) to regulate the expression of multiple virulence proteins during acute and persistent infections. The present study targeted MvfR with the intention of designing novel anti-virulent compounds, which will function in two ways: first, they will block the virulence and pathogenesis P. aeruginosa by disrupting the quorum-sensing network of the bacteria, and second, they will stop the evolution of new resistant mechanisms. A structure-based virtual screening (SBVS) method was used to screen druglike compounds from the Asinex antibacterial library (~5968 molecules) and the comprehensive marine natural products database (CMNPD) (~32 thousand compounds), against the ligand-binding domain (LBD) of MvfR, to identify molecules that show high binding potential for the relevant pocket. In this way, two compounds were identified: Top-1 (4-((carbamoyloxy)methyl)-10,10-dihydroxy-2,6-diiminiodecahydropyrrolo[1,2-c]purin-9-yl sulfate) and Top-2 (10,10-dihydroxy-2,6-diiminio-4-(((sulfonatocarbamoyl)oxy)methyl)decahydropyrrolo[1,2-c]purin-9-yl sulfate), in contrast to the co-crystallized M64 control. Both of the screened leads were found to show deep pocket binding and interactions with several key residues through a network of hydrophobic and hydrophilic interactions. The docking results were validated by a long run of 200 ns of molecular dynamics simulation and MM-PB/GBSA binding free energies. All of these analyses confirmed the presence of strong complex formation and rigorous intermolecular interactions. An additional analysis of normal mode entropy and a WaterSwap assay were also performed to complement the aforementioned studies. Lastly, the compounds were found to show an acceptable range of pharmacokinetic properties, making both compounds potential candidates for further experimental studies to decipher their real biological potency.

In Silico Core Proteomics and Molecular Docking Approaches for the Identification of Novel Inhibitors against Streptococcus pyogenes

Streptococcus pyogenes is a significant pathogen that causes skin and upper respiratory tract infections and non-suppurative complications, such as acute rheumatic fever and post-strep glomerulonephritis. Multidrug resistance has emerged in S. pyogenes strains, making them more dangerous and pathogenic. Hence, it is necessary to identify and develop therapeutic methods that would present novel approaches to S. pyogenes infections. In the current study, a subtractive proteomics approach was employed to core proteomes of four strains of S. pyogenes using several bioinformatic software tools and servers. The core proteome consists of 1324 proteins, and 302 essential proteins were predicted from them. These essential proteins were analyzed using BLASTp against human proteome, and the number of potential targets was reduced to 145. Based on subcellular localization prediction, 46 proteins with cytoplasmic localization were chosen for metabolic pathway analysis. Only two cytoplasmic proteins, i.e., chromosomal replication initiator protein DnaA and two-component response regulator (TCR), were discovered to have the potential to be novel drug target candidates. Three-dimensional (3D) structure prediction of target proteins was carried out via the Swiss Model server. Molecular docking approach was employed to screen the library of 1000 phytochemicals against the interacting residues of the target proteins through the MOE software. Further, the docking studies were validated by running molecular dynamics simulation and highly popular binding free energy approaches of MM-GBSA and MM-PBSA. The findings revealed a promising candidate as a novel target against S. pyogenes infections.

Molecular Docking and Simulation Studies of Antidiabetic Agents Devised from Hypoglycemic Polypeptide-P of Momordica charantia

Diabetes mellitus termed as metabolic disorder is a collection of interlinked diseases and mainly body's inability to manage glucose level which leads to cardiovascular diseases, renal failure, neurological disorders, and many others. The drugs contemporarily used for diabetes have many inevitable side effects, and many of them have become less responsive to this multifactorial disorder. Momordica charantia commonly known as bitter gourd has many bioactive compounds with antidiabetic properties. The current study was designed to use computational methods to discover the best antidiabetic peptides devised from hypoglycemic polypeptide-P of M. charantia. The binding affinity and interaction patterns of peptides were evaluated against four receptor proteins (i.e., as agonists of insulin receptor and inhibitors of sodium-glucose cotransporter 1, dipeptidyl peptidase-IV, and glucose transporter 2) using molecular docking approach. A total of thirty-seven peptides were docked against these receptors. Out of which, top five peptides against each receptor were shortlisted based on their S-scores and binding affinities. Finally, the eight best ligands (i.e., LIVA, TSEP, EKAI, LKHA, EALF, VAEK, DFGAS, and EPGGGG) were selected as these ligands strictly followed Lipinski's rule of five and exhibited good ADMET profiling. One peptide EPGGGG showed activity towards insulin and SGLT1 receptor proteins. The top complex for both these targets was subjected to 50 ns of molecular dynamics simulations and MM-GBSA binding energy test that concluded both complexes as highly stable, and the intermolecular interactions were dominated by van der Waals and electrostatic energies. Overall, the selected ligands strongly fulfilled the drug-like evaluation criterion and proved to have good antidiabetic properties.

Structural insight into the binding pattern and interaction mechanism of chemotherapeutic agents with Sorcin by docking and molecular dynamic simulation

Sorcin (SOluble Resistance-related Calcium bInding proteiN) is a calcium binding protein that plays a key role in multidrug resistance (MDR) in human cancers. This study aimed at understanding the binding mechanism and structural basis for the interaction of structurally and functionally unrelated chemotherapeutic agent, namely doxorubicin, etoposide, omacetaxine mepesuccinate and paclitaxel with Sorcin by utilizing docking and molecular dynamic simulation approaches. The docking evaluation of etoposide, omacetaxine mepesuccinate and paclitaxel have shown a high affinity binding with Sorcin at the Ca²⁺-binding C-terminal domain (SCBD) in a comparable mode and affinity of binding to doxorubicin. Moreover, all of the docked compounds were shown to interact both hydrophilically and hydrophobically with the same residues within the active pocket which is located at interface of the Sorcin and collectively formed by EF5 loop, G helix and EF4 loop. However, the MD simulations revealed that the dynamics of Sorcin structure is different in the presence of the compounds when compared and contrasted to the Apo Sorcin, particularly in the first 25 ns, after which each system gained considerable structure stability. The difference in dynamics might be the outcome of high N and C-terminal flexibility that seem not to disturb compounds binding conformation but more likely is affecting chemical interaction network by breaking and establishing old and new hydrogen bonds, respectively. This detailed mechanistic understanding of different chemotherapeutic agents binding to Sorcin might be useful to open windows for designing and developing new inhibitors that are potentially capable of reversing the MDR in human cancers.

Structural basis of UDP-N-acetylglucosamine pyrophosphorylase and identification of promising terpenes to control Aedes aegypti

According to the world health organization (WHO) 2020 report, vector borne diseases account for 17 % of all infections with reported 700 thousand death each year. They are of considerable importance to health professionals as they are posing a serious health threat and include dengue fever, Zika fever, chikungunya, yellow fever, and other disease agents. Aedes aegypti serve as a vector for transmitting several of these tropical fevers. In the present study, UDP-N-acetylglucosamine pyrophosphorylase enzyme (Aa-UAP) of A. aegypti which plays a significant contribution in chitin metabolism is targeted with natural terpenes to propose an eco-friendly and novel candidates for the development of new insecticides. The three dimensional Aa-UAP structure was constructed via a comparative homology approach and validated, followed by structure-based virtual screening against 1000 terpenes collected from natural MDP3 and NPACT databases. Top hits were subjected to molecular dynamics (MD) simulations and binding free energies analysis to elucidate complex intermolecular stability and affinity over simulated time. The results demonstrated that Aa-UAP possesses a homodimer state and its active site residues are well conserved. Three compounds (NPACT00138, NPACT00452, and NPACT00839) were prioritized as they are establishing conserved and stable interactions with the active binding-site residues of Aa-UAP. Conclusively, the reported Aa-UAP specific terpenes could serve as promising leads in order to develop potential insecticides. Importantly, the FDA approved drug NPACT00839 (Paclitaxel) could be used further in the fast-track experimental testing pipeline for biological optimization.