A Review on Applications of Computational Methods in Drug Screening and Design

Enhanced Sampling in Molecular Dynamics Simulations: How Many MD Snapshots can be Needed to Reproduce the Biological Behavior?

Since its early days in the 19th century, medicinal chemistry has concentrated its efforts on the treatment of diseases, using tools from areas such as chemistry, pharmacology, and molecular biology. The understanding of biological mechanisms and signaling pathways is crucial information for the development of potential agents for the treatment of diseases mainly because they are such complex processes. Given the limitations that the experimental approach presents, computational chemistry is a valuable alternative for the study of these systems and their behavior. Thus, classical molecular dynamics, based on Newton's laws, is considered a technique of great accuracy, when appropriated force fields are used, and provides satisfactory contributions to the scientific community. However, as many configurations are generated in a large MD simulation, methods such as Statistical Inefficiency and Optimal Wavelet Signal Compression Algorithm are great tools that can reduce the number of subsequent QM calculations. Accordingly, this review aims to briefly discuss the importance and relevance of medicinal chemistry allied to computational chemistry as well as to present a case study where, through a molecular dynamics simulation of AMPK protein (50 ns) and explicit solvent (TIP3P model), a minimum number of snapshots necessary to describe the oscillation profile of the protein behavior was proposed. For this purpose, the RMSD calculation, together with the sophisticated OWSCA method was used to propose the minimum number of snapshots.

The Discovery of Novel Agents against Staphylococcus aureus by Targeting Sortase A: A Combination of Virtual Screening and Experimental Validation

Staphylococcus aureus (S. aureus), commonly known as "superbugs", is a highly pathogenic bacterium that poses a serious threat to human health. There is an urgent need to replace traditional antibiotics with novel drugs to combat S. aureus. Sortase A (SrtA) is a crucial transpeptidase involved in the adhesion process of S. aureus. The reduction in virulence and prevention of S. aureus infections have made it a significant target for antimicrobial drugs. In this study, we combined virtual screening with experimental validation to identify potential drug candidates from a drug library. Three hits, referred to as Naldemedine, Telmisartan, and Azilsartan, were identified based on docking binding energy and the ratio of occupied functional sites of SrtA. The stability analysis manifests that Naldemedine and Telmisartan have a higher binding affinity to the hydrophobic pockets. Specifically, Telmisartan forms stable hydrogen bonds with SrtA, resulting in the highest binding energy. Our experiments prove that the efficiency of adhesion and invasion by S. aureus can be decreased without significantly affecting bacterial growth. Our work identifies Telmisartan as the most promising candidate for inhibiting SrtA, which can help combat S. aureus infection.

Annonaceous acetogenins as promising DNA methylation inhibitors to prevent and treat leukemogenesis - an in silico approach

Leukemia is a haematological malignancy affecting blood and bone marrow, ranking 10th among the other common cancers. DNA methylation is an epigenetic dysregulation that plays a critical role in leukemogenesis. DNA methyltransferases (DNMTs) such as DNMT1, DNMT3A and DNMT3B are the key enzymes catalysing DNA methylation. Inhibition of DNMT1 with secondary metabolites from medicinal plants helps reverse DNA methylation. The present study focuses on inhibiting DNMT1 protein (PDB ID: 3PTA) with annonaceous acetogenins through in-silico studies. The docking and molecular dynamic (MD) simulation study was carried out using Schrödinger Maestro and Desmond, respectively. These compounds' drug likeliness, ADMET properties and bioactivity scores were analysed. About 76 different acetogenins were chosen for this study, among which 17 showed the highest binding energy in the range of -8.312 to -10.266 kcal/mol. The compounds with the highest negative binding energy were found to be annohexocin (-10.266 kcal/mol), isoannonacinone (-10.209 kcal/mol) and annonacin (-9.839 kcal/mol). MD simulation results reveal that annonacin remains stable throughout the simulation time of 100 ns and also binds to the catalytic domain of DNMT1 protein. From the above results, it can be concluded that annonacin has the potential to inhibit the DNA methylation process and prevent leukemogenesis.Communicated by Ramaswamy H. Sarma.

Scipion-Chem: An Open Platform for Virtual Drug Screening

Virtual drug screening (VDS) tackles the problem of drug discovery by computationally reducing the number of potential pharmacological molecules that need to be tested experimentally to find a new drug. To do so, several approaches have been developed through the years, typically focusing on either the physicochemical characteristics of the receptor structure (structure-based virtual screening) or those of the potential ligands (ligand-based virtual screening). Scipion is a workflow engine well suited for structural studies of biological macromolecules. Here, we present Scipion-chem, a new branch oriented to VDS. A total of 11 plugins have already been integrated from the most common programs used in the field. They can be used through the Scipion graphical user interface to execute and analyze typical VDS tasks. In addition, we have developed several consensus protocols that combine results from the different integrated programs to generate more robust predictions. Backstage, Scipion also facilitates the interoperability of those different software packages while tracking all of the intermediate files, parameters, and user decisions. In summary, in this article, we present Scipion-chem. This accessible, interoperable, and traceable platform provides the user with all of the tools to carry out a successful VDS workflow. Scipion-chem is openly available at https://github.com/scipion-chem.

Improved prediction of protein-protein interactions by a modified strategy using three conventional docking software in combination

Proteins play a crucial role in many biological processes, where their interaction with other proteins are integral. Abnormal protein-protein interactions (PPIs) have been linked to various diseases including cancer, and thus targeting PPIs holds promise for drug development. However, experimental confirmation of the peculiarities of PPIs is challenging due to their dynamic and transient nature. As a complement to experimental technologies, multiple computational molecular docking (MD) methods have been developed to predict the structures of protein-protein complexes and their dynamics, still requiring further improvements in several issues. Here, we report an improved MD method, namely three-software docking (3SD), by employing three popular protein-peptide docking software (CABS-dock, HPEPDOCK, and HADDOCK) in combination to ensure constant quality for most targets. We validated our 3SD performance in known protein-peptide interactions (PpIs). We also enhanced MD performance in proteins having intrinsically disordered regions (IDRs) by applying the modified 3SD strategy, the three-software docking after removing random coiled IDR (3SD-RR), to the comparable crystal PpI structures. At the end, we applied 3SD-RR to the AlphaFold2-predicted receptors, yielding an efficient prediction of PpI pose with high relevance to the experimental data regardless of the presence of IDRs or the availability of receptor structures. Our study provides an improved solution to the challenges in studying PPIs through computational docking and has the potential to contribute to PPIs-targeted drug discovery. SIGNIFICANCE STATEMENT: Protein-protein interactions (PPIs) are integral to life, and abnormal PPIs are associated with diseases such as cancer. Studying protein-peptide interactions (PpIs) is challenging due to their dynamic and transient nature. Here we developed improved docking methods (3SD and 3SD-RR) to predict the PpI poses, ensuring constant quality in most targets and also addressing issues like intrinsically disordered regions (IDRs) and artificial intelligence-predicted structures. Our study provides an improved solution to the challenges in studying PpIs through computational docking and has the potential to contribute to PPIs-targeted drug discovery.

Exploring the anticancer properties of a new nicotinamide analogue: Investigations into in silico analysis, antiproliferative effects, selectivity, VEGFR-2 inhibition, apoptosis induction, and migration suppression

BACKGROUND

This study focuses on the development and evaluation of (E)-N-(3-(1-(2-(4-bromobenzoyl)hydrazono)ethyl)phenyl)nicotinamide (BHEPN) as a potential inhibitor of Vascular Endothelial Growth Factor Receptor-2 (VEGFR-2).

METHODS

Computational investigations as density function theory (DFT), docking, molecular dynamics (MD) simulations, and ADMET) in addition to in vitro (VEGFR-2 inhibition, cytotoxicity against HepG2 and MCF-7 cancer cell lines, selectivity index, cells cycle analysis, apoptosis investigation, and cells migration assay) studies were conducted.

RESULTS

DFT calculations determined the three-dimensional structure and indicated the reactivity of BHEPN. Molecular docking, and MD simulations analysis showed the BHEPN's binding affinity and its potential as a VEGFR-2 inhibitor. ADMET assessments predicted BHEPN's safety and drug-like characteristics. In vitro investigations confirmed the inhibition of VEGFR-2 with an IC50 value of 0.320 ± 0.012 µM. BHEPN also exhibited remarkable cytotoxic effects against HepG2 and MCF-7 cancer cell lines, with IC50 values of 0.19 ± 0.01 µM and 1.18 ± 0.01 µM, respectively, outperforming Sorafenib's IC50 values (2.24 ± 0.06 µM and 3.17 ± 0.01 µM), respectively. Notably, BHEPN displayed a higher IC50 value of 4.11 ± 0 µM against the non-carcinogenic Vero cell lines, indicating selectivity index values of 21.6 and 3.4 against the tested cancer cell lines, respectively. In a flow cytometry assay, BHEPN induced HepG2 cell cycle arrest at the G1/S phase. Moreover, BHEPN increased the incidence of early and late apoptosis in HepG2 cell lines (from 1.38% and 0.22%) in control cells to (4.11-26.02%) in the treated cells, respectively. Additionally, the percentage of necrosis raised to 13.39%, in contrast to 0.62% in control cells. Finally, BHEPN was able to reduce the migration and wound healing abilities in HepG2 cells to 38.89% compared to 87.92% in untreated cells after 48 h. These in vitro results aligned with the computational predictions, providing strong evidence of BHEPN's efficacy and safety in anticancer applications.

CONCLUSIONS

BHEPN is a promising candidate for the development of novel anticancer agents through further in vitro and in vivo investigations.

Computer-Assisted Drug Discovery of a Novel Theobromine Derivative as an EGFR Protein-Targeted Apoptosis Inducer

The overexpression of the Epidermal Growth Factor Receptor (EGFR) marks it as a pivotal target in cancer treatment, with the aim of reducing its proliferation and inducing apoptosis. This study aimed at the CADD of a new apoptotic EGFR inhibitor. The natural alkaloid, theobromine, was used as a starting point to obtain a new semisynthetic (di-ortho-chloro acetamide) derivative (T-1-DOCA). Firstly, T-1-DOCA's total electron density, energy gap, reactivity indices, and electrostatic surface potential were determined by DFT calculations, Then, molecular docking studies were carried out to predict the potential of T-1-DOCA against wild and mutant EGFR proteins. T-1-DOCA's correct binding was further confirmed by molecular dynamics (MD) over 100 ns, MM-GPSA, and PLIP experiments. In vitro, T-1-DOCA showed noticeable efficacy compared to erlotinib by suppressing EGFRWT and EGFRT790M with IC50 values of 56.94 and 269.01 nM, respectively. T-1-DOCA inhibited also the proliferation of H1975 and HCT-116 malignant cell lines, exhibiting IC50 values of 14.12 and 23.39 µM, with selectivity indices of 6.8 and 4.1, respectively, indicating its anticancer potential and general safety. The apoptotic effects of T-1-DOCA were indicated by flow cytometric analysis and were further confirmed through its potential to increase the levels of BAX, Casp3, and Casp9, and decrease Bcl-2 levels. In conclusion, T-1-DOCA, a new apoptotic EGFR inhibitor, was designed and evaluated both computationally and experimentally. The results suggest that T-1-DOCA is a promising candidate for further development as an anti-cancer drug.

Identification of Potential Therapeutic Targets for Myopic Choroidal Neovascularization via Discovery-Driven Data Mining

Purpose: Myopic choroidal neovascularization (mCNV) is a prevalent cause of vision loss. However, the development of effective therapeutic targets for mCNV has been hindered by the paucity of suitable animal models. Therefore, the aim of this study is to identify potential genes and pathways associated with mCNV and to unearth prospective therapeutic targets that can be utilized to devise efficacious treatments.Methods: Text data mining was used to identify genes linked to choroid, neovascularization, and myopia. g: Profiler was utilized to analyze the biological processes of gene ontology and the Reactome pathways. Protein interaction network analysis was performed using strings and visualized in Cytoscape. MCODE and cytoHubba were used for further screening.Results: Discovery-driven text data mining identified 55 potential genes related to choroid, neovascularization, and myopia. Gene enrichment analysis revealed 11 biological processes and seven Reactome pathways. A protein-protein interaction network with 47 nodes was constructed and analyzed using centrality ranking. Key clusters were identified through algorithm tools. Finally, 14 genes (IL6, FGF2, MMP9, IL10, TNF, MMP2, HGF, MMP3, IGF1, CCL2, CTNNB1, BDNF, NGF, and EDN1), in addition to VEGFA, were evaluated as targets with potential as future therapeutics.Conclusions: This study provides new potential therapeutic targets for mCNV, including IL6, FGF2, MMP9, IL10, TNF, MMP2, HGF, MMP3, IGF1, CCL2, CTNNB1, BDNF, NGF, and EDN1, which correspond to seven potential enriched pathways. These findings provide a basis for further research and offer new possibilities for developing therapeutic interventions for this condition.

Target-ligand binding affinity from single point enthalpy calculation and elemental composition

Reliable target-ligand binding thermodynamics data are essential for successful drug design and molecular engineering projects. Besides experimental methods, a number of theoretical approaches have been introduced for the generation of binding thermodynamics data. However, available approaches often neglect electronic effects or explicit water molecules influencing target-ligand interactions. To handle electronic effects within a reasonable time frame, we introduce a fast calculator QMH-L using a single target-ligand complex structure pre-optimized at the molecular mechanics level. QMH-L is composed of the semi-empirical quantum mechanics calculation of binding enthalpy with predicted explicit water molecules at the complex interface, and a simple descriptor based on the elemental composition of the ligand. QMH-L estimates the target-ligand binding free energy with a root mean square error (RMSE) of 0.94 kcal mol-1. The calculations also provide binding enthalpy values and they were compared with experimental binding thermodynamics data collected from the most reliable isothermal titration calorimetry studies of systems including various protein targets and challenging, large peptide ligands with a molecular weight of up to 2-3 thousand. The single point enthalpy calculations of QMH-L require modest computational resources and are based on short runs with open source and/or free software like Gromacs, Mopac, MobyWat, and Fragmenter. QMH-L can be applied for fast, automated scoring of drug candidates during a virtual screen, enthalpic engineering of new ligands or thermodynamic explanation of complex interactions.

Pharmacophore-Based Screening, Molecular Docking, and Dynamic Simulation of Fungal Metabolites as Inhibitors of Multi-Targets in Neurodegenerative Disorders

Neurodegenerative disorders, such as Alzheimer's disease (AD), negatively affect the economic and psychological system. For AD, there is still a lack of disease-altering treatments and promising cures due to its complex pathophysiology. In this study, we computationally screened the natural database of fungal metabolites against three known therapeutic target proteins of AD. Initially, a pharmacophore-based, drug-likeness category was employed for screening, and it filtered the 14 (A-N) best hits out of 17,544 fungal metabolites. The 14 best hits were docked individually against GSK-3β, the NMDA receptor, and BACE-1 to investigate the potential of finding a multitarget inhibitor. We found that compounds B, F, and L were immuno-toxic, whereas E, H, I, and J had a higher LD50 dose (5000 mg/kg). Among the examined metabolites, the Bisacremine-C (compound I) was found to be the most active molecule against GSK-3β (ΔG: -8.7 ± 0.2 Kcal/mol, Ki: 2.4 × 106 M-1), NMDA (ΔG: -9.5 ± 0.1 Kcal/mol, Ki: 9.2 × 106 M-1), and BACE-1 (ΔG: -9.1 ± 0.2 Kcal/mol, Ki: 4.7 × 106 M-1). It showed a 25-fold higher affinity with GSK-3β, 6.3-fold higher affinity with NMDA, and 9.04-fold higher affinity with BACE-1 than their native ligands, respectively. Molecular dynamic simulation parameters, such as RMSD, RMSF, Rg, and SASA, all confirmed that the overall structures of the targeted enzymes did not change significantly after binding with Bisacremine-C, and the ligand remained inside the binding cavity in a stable conformation for most of the simulation time. The most significant hydrophobic contacts for the GSK-3β-Bisacremine-C complex are with ILE62, VAL70, ALA83, and LEU188, whereas GLN185 is significant for H-bonds. In terms of hydrophobic contacts, TYR184 and PHE246 are the most important, while SER180 is vital for H-bonds in NMDA-Bisacremine-C. THR232 is the most crucial for H-bonds in BACE-1-Bisacremine-C and ILE110-produced hydrophobic contacts. This study laid a foundation for further experimental validation and clinical trials regarding the biopotency of Bisacremine-C.

A Quantitative Structure-Activity Relationship for Human Plasma Protein Binding: Prediction, Validation and Applicability Domain

Purpose

The purpose of this study was to develop a robust and externally predictive in silico QSAR-neural network model for predicting plasma protein binding of drugs. This model aims to enhance drug discovery processes by reducing the need for chemical synthesis and extensive laboratory testing.

Methods

A dataset of 277 drugs was used to develop the QSAR-neural network model. The model was constructed using a Filter method to select 55 molecular descriptors. The validation set's external accuracy was assessed through the predictive squared correlation coefficient Q2 and the root mean squared error (RMSE).

Results

The developed QSAR-neural network model demonstrated robustness and good applicability domain. The external accuracy of the validation set was high, with a predictive squared correlation coefficient Q2 of 0.966 and a root mean squared error (RMSE) of 0.063. Comparatively, this model outperformed previously published models in the literature.

Conclusion

The study successfully developed an advanced QSAR-neural network model capable of predicting plasma protein binding in human plasma for a diverse set of 277 drugs. This model's accuracy and robustness make it a valuable tool in drug discovery, potentially reducing the need for resource-intensive chemical synthesis and laboratory testing.

Peptidomimetics for CVD screened via TRADD-TRAF2 complex interface assessments

The main aim of this study is to screen and develop Peptidomimetics to treat atherosclerosis (AS) which is a Cardio Vascular Disease (CVD). Peptidomimetics were obtained from the protein-protein interaction interface of TRADD (Tumor necrosis factor receptor type 1-associated DEATH domain protein) and TRAF2 (TNF receptor-associated factor 2) complex. TRADD-TRAF2 interaction is critical in AS pathogenesis since it assists a series of signal transducers that activate NF-κB. Conceptually, the triggered NF-κB makes an extensive amount of nitric oxide (NO) synthesized by inducible nitric oxide synthase (iNOS), which boons the progress of AS. The examined TRADD-TRAF2 complex (PDB ID: 1F3V) and its interaction details revealed that the sequence range W11-G165 of TRADD highly interacts with TRAF2. The sequence range W11-G165 was selected for the design and preparation of the inhibitory peptide in silico. The selected sequence was mutated with the alanine scanning method to have a range of inhibitory peptides. With the help of different in silico tools, the top three, namely MIP11-25 L, MIP131-143 h, and MIP149-164 m peptides showed the best interaction with the critical residues of TRAF2. Thus, these three peptides were used for generating peptidomimetics using pepMMsMIMIC, a peptidomimetics virtual screening tool. Around 600 peptidomimetics were identified & and retrieved for further screening by employing molecular docking tools and MD (Molecular Dynamics) simulations. Density Functional Theory (DFT) and ADMET predictions were applied to validate the screened peptidomimetics druggability. In the results, peptidomimic compounds MMs03918858 and MMs03927281 with binding energy values of -9.6 kcal/mol and - 9.1 kcal/mol respectively were screened as the best and are proposed for further pre-clinical studies.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Using CADD tools to inhibit the overexpressed genes FAP, FN1, and MMP1 by repurposing ginsenoside C and Rg1 as a treatment for oral cancer

Oral cancer is one of the most common cancer types. Many factors can express certain genes that cause the proliferation of oral tissues. Overexpressed genes were detected in oral cancer patients; three were highly impacted. FAP, FN1, and MMP1 were the targeted genes that showed inhibition results in silico by ginsenoside C and Rg1. Approved drugs were retrieved from the DrugBank database. The docking scores show an excellent interaction between the ligands and the targeted macromolecules. Further molecular dynamics simulations showed the binding stability of the proposed natural products. This work recommends repurposing ginsenoside C and Rg1 as potential binders for the selected targets and endorses future experimental validation for the treatment of oral cancer.

Reformulating Reactivity Design for Data-Efficient Machine Learning

Machine learning (ML) can deliver rapid and accurate reaction barrier predictions for use in rational reactivity design. However, model training requires large data sets of typically thousands or tens of thousands of barriers that are very expensive to obtain computationally or experimentally. Furthermore, bespoke data sets are required for each region of interest in reaction space as models typically struggle to generalize. We have therefore reformulated the ML barrier prediction problem toward a much more data-efficient process: finding a reaction from a prespecified set with a desired target value. Our reformulation enables the rapid selection of reactions with purpose-specific activation barriers, for example, in the design of reactivity and selectivity in synthesis, catalyst design, toxicology, and covalent drug discovery, requiring just tens of accurately measured barriers. Importantly, our reformulation does not require generalization beyond the domain of the data set at hand, and we show excellent results for the highly toxicologically and synthetically relevant data sets of aza-Michael addition and transition-metal-catalyzed dihydrogen activation, typically requiring less than 20 accurately measured density functional theory (DFT) barriers. Even for incomplete data sets of E2 and SN2 reactions, with high numbers of missing barriers (74% and 56% respectively), our chosen ML search method still requires significantly fewer data points than the hundreds or thousands needed for more conventional uses of ML to predict activation barriers. Finally, we include a case study in which we use our process to guide the optimization of the dihydrogen activation catalyst. Our approach was able to identify a reaction within 1 kcal mol-1 of the target barrier by only having to run 12 DFT reaction barrier calculations, which illustrates the usage and real-world applicability of this reformulation for systems of high synthetic importance.

Role of Computer-Aided Drug Design in Drug Development

The introduction of computational techniques to pharmaceutical chemistry and molecular biology in the 20th century has changed the way people develop drugs [...].

Gargoyles: An Open Source Graph-Based Molecular Optimization Method Based on Deep Reinforcement Learning

Automatic optimization methods for compounds in the vast compound space are important for drug discovery and material design. Several machine learning-based molecular generative models for drug discovery have been proposed, but most of these methods generate compounds from scratch and are not suitable for exploring and optimizing user-defined compounds. In this study, we developed a compound optimization method based on molecular graphs using deep reinforcement learning. This method searches for compounds on a fragment-by-fragment basis and at high density by generating fragments to be added atom by atom. Experimental results confirmed that the quantum electrodynamics (QED), the optimization target set in this study, was enhanced by searching around the starting compound. As a use case, we successfully enhanced the activity of a compound by targeting dopamine receptor D2 (DRD2). This means that the generated compounds are not structurally dissimilar from the starting compounds, as well as increasing their activity, indicating that this method is suitable for optimizing molecules from a given compound. The source code is available at https://github.com/sekijima-lab/GARGOYLES.

Identification of potential novel therapeutic drug target against Elizabethkingia anophelis by integrative pan and subtractive genomic analysis: An in silico approach

Elizabethkingia anophelis is a human pathogen responsible for severe nosocomial infections in neonates and immunocompromised patients. The significantly higher mortality rate from E. anophelis infections and the lack of available regimens highlight the critical need to explore novel drug targets. The current study investigated effective novel drug targets by employing a comprehensive in silico subtractive genomic approach integrated with pangenomic analysis of E. anophelis strains. A total of 2809 core genomic proteins were found by pangenomic analysis of non-paralogous proteins. Subsequently, 156 pathogen-specific, 442 choke point, 202 virulence factor, 53 antibiotic resistant and 119 host-pathogen interacting proteins were identified in E. anophelis. By subtractive genomic approach, at first 791 proteins were found to be indispensable for the survival of E. anophelis. 558 and 315 proteins were detected as non-homologous to human and gut microflora respectively. Following that 245 cytoplasmic, 245 novel, and 23 broad-spectrum targets were selected and finally four proteins were considered as potential therapeutic targets of E. anophelis based on highest degree score in PPI network. Among those, three proteins were subjected to molecular docking and subsequent MD simulation as one protein did not contain a plausible binding pocket with sufficient surface area and volume. All the complexes were found to be stable and compact in 100 ns molecular dynamics simulation studies as measured by RMSD, RMSF, and Rg. These three short-listed targets identified in this study may lead to the development of novel antimicrobials capable of curing infections and pave the way to prevent and control the disease progression caused by the deadly agent E. anophelis.

Pharmacophore modeling, molecular docking, and molecular dynamics studies to identify new 5-HT2AR antagonists with the potential for design of new atypical antipsychotics

Some important atypical antipsychotic drugs target the serotonergic receptor 2A (5-HT2AR). Currently, new therapeutic strategies are needed to offer faster onset of action with fewer side effects and, therefore, greater efficacy in a substantial proportion of patients with neuropsychological disorders such as Autism and Parkinson. The main objective of this work was to use SBDD methods to identify new hit compounds potentially useful as precursors of novel and selective 5-HT2AR antagonists. A structure-based pharmacophore screening study based on a selective antagonist was carried out in ten databases. The set obtained was refined using molecular docking, and the five most promising compounds were subjected to molecular dynamics simulations. The most stable and promising hit occupied a side pocket present in the 5-HT2AR, a site that can be explored to obtain selective ligands. Simulations against 5-HT2CR and D2R showed that the best hit could not form stable complexes with these targets, strengthening the hypothesis that the hit presents selective binding by the receptor of interest. The selected hits showed some predicted toxicity risk or violated some drug-likeness property. However, it can be concluded that the identified hits are the most promising for performing in vitro assays. Once the presence of activity is confirmed, they could become precursors of optimized and selective antagonists of 5-HT2AR. An SBDD study was carried out to identify new selective 5-HT2AR ligands potentially useful for designing selective atypical antipsychotics.

Synthesis, in silico and in vitro antimicrobial efficacy of substituted arylidene-based quinazolin-4(3H)-one motifs

Introduction: Quinazolin-4(3H)-one derivatives have attracted considerable attention in the pharmacological profiling of therapeutic drug targets. The present article reveals the development of arylidene-based quinazolin-4(3H)-one motifs as potential antimicrobial drug candidates. Methods: The synthetic pathway was initiated through thermal cyclization of acetic anhydride on anthranilic acid to produce 2-methyl-4H-3,1-benzoxazan-4-one 1, which (upon condensation with hydrazine hydrate) gave 3-amino-2-methylquinazolin-4(3H)-one 2. The reaction of intermediate 2 at its amino side arm with various benzaldehyde derivatives furnished the final products, in the form of substituted benzylidene-based quinazolin-4(3H)-one motifs 3a-l, and with thiophene-2-carbaldehyde to afford 3 m. The purified targeted products 3a-m were effectively characterized for structural authentication using physicochemical parameters, microanalytical data, and spectroscopic methods, including IR, UV, and 1H- and 13C-NMR, as well as mass spectral data. The substituted arylidene-based quinazolin-4(3H)-one motifs 3a-m were screened for both in silico and in vitro antimicrobial properties against selected bacteria and fungi. The in silico studies carried out consisted of predicted ADMET screening, molecular docking, and molecular dynamics (MD) simulation studies. Furthermore, in vitro experimental validation was performed using the agar diffusion method, and the standard antibacterial and antifungal drugs used were gentamicin and ketoconazole, respectively. Results and discussion: Most of the compounds possessed good binding affinities according to the molecular docking studies, while MD simulation revealed their levels of structural stability in the protein-ligand complexes. 2-methyl-3-((thiophen-2-ylmethylene)amino) quinazolin-4(3H)-one 3 m emerged as both the most active antibacterial agent (with an minimum inhibitory concentration (MIC) value of 1.95 μg/mL) against Staphylococcus aureus and the most active antifungal agent (with an MIC value of 3.90 μg/mL) against Candida albicans, Aspergillus niger, and Rhizopus nigricans.

SubMDTA: drug target affinity prediction based on substructure extraction and multi-scale features

BACKGROUND

Drug-target affinity (DTA) prediction is a critical step in the field of drug discovery. In recent years, deep learning-based methods have emerged for DTA prediction. In order to solve the problem of fusion of substructure information of drug molecular graphs and utilize multi-scale information of protein, a self-supervised pre-training model based on substructure extraction and multi-scale features is proposed in this paper.

RESULTS

For drug molecules, the model obtains substructure information through the method of probability matrix, and the contrastive learning method is implemented on the graph-level representation and subgraph-level representation to pre-train the graph encoder for downstream tasks. For targets, a BiLSTM method that integrates multi-scale features is used to capture long-distance relationships in the amino acid sequence. The experimental results showed that our model achieved better performance for DTA prediction.

CONCLUSIONS

The proposed model improves the performance of the DTA prediction, which provides a novel strategy based on substructure extraction and multi-scale features.

Machine Learning-Based Approach to Developing Potent EGFR Inhibitors for Breast Cancer-Design, Synthesis, and In Vitro Evaluation

The epidermal growth factor receptor (EGFR) is vital for regulating cellular functions, including cell division, migration, survival, apoptosis, angiogenesis, and cancer. EGFR overexpression is an ideal target for anticancer drug development as it is absent from normal tissues, marking it as tumor-specific. Unfortunately, the development of medication resistance limits the therapeutic efficacy of the currently approved EGFR inhibitors, indicating the need for further development. Herein, a machine learning-based application that predicts the bioactivity of novel EGFR inhibitors is presented. Clustering of the EGFR small-molecule inhibitor (∼9000 compounds) library showed that N-substituted quinazolin-4-amine-based compounds made up the largest cluster of EGFR inhibitors (∼2500 compounds). Taking advantage of this finding, rational drug design was used to design a novel series of 4-anilinoquinazoline-based EGFR inhibitors, which were first tested by the developed artificial intelligence application, and only the compounds which were predicted to be active were then chosen to be synthesized. This led to the synthesis of 18 novel compounds, which were subsequently evaluated for cytotoxicity and EGFR inhibitory activity. Among the tested compounds, compound 9 demonstrated the most potent antiproliferative activity, with 2.50 and 1.96 μM activity over MCF-7 and MDA-MB-231 cancer cell lines, respectively. Moreover, compound 9 displayed an EGFR inhibitory activity of 2.53 nM and promising apoptotic results, marking it a potential candidate for breast cancer therapy.

Amphiphilic peptide Mastoparan-B induces conformational changes within the AdeB efflux pump, down-regulates adeB gene expression, and restores antibiotic susceptibility in an MDR strain of Acinetobacter baumannii

Mastoparan B (MP-B) is an amphiphilic peptide with a potent antimicrobial activity against most Gram-negative bacteria. However, there is little information available on the inhibition of the Acinetobacter baumannii resistance-nodulation-cell-division (RND) efflux pump using this antimicrobial peptide. Here, we carried out a series of in-silico experiments to find the mechanisms underlying the anti-efflux activity of MP-B using a multi-drug resistant (MDR) strain of A. baumannii (AB). According to our findings, MP-B demonstrated a potent antibacterial activity against an MDR-AB (minimum inhibitory concentration [MIC] = 1 μg/mL) followed by a 20-fold reduction in the adeB gene expression in the presence of sub-MIC of this peptide. Using Groningen Machine for Chemicals Simulation (GROMACS) via PyMOL Graphical User Interface (GUI), (we observed that, the AdeB transporter had conserved helix-turn-helix regions and a tight pore rich in Phe and Ala residues. To understand how inhibition of the AdeB is achieved, we generated 20 apo-MP-B poses using the InterPep and SiteMap tools. The high-quality model was created by homology modeling and used for docking via AutoDock/Vina to identify the MP-B binding sites. We established that the most apo-MP-B formed H-bonds to the backbone of five amino acids in the Helix-5. As a result, the dihedral angles of the involved amino acids shift by 9.0-9.6 Ǻ, causing a change in the conformation of the AdeB protein. This led to helix conformation stereoisomerization and block the AdeB activity. MP-B presumably has dual mechanisms. (1) It blocks the AdeB transporter by changing its conformation. (2) MP-B influences the adeB gene expression by binding to G-protein which laterally controls efflux regulators like MarA, RamA, SoxS, and Rob proteins.

Phytochemicals of Cocculus hirsutus deciphered SARS-CoV-2 inhibition by targeting main proteases in molecular docking, simulation, and pharmacological analyses

The COVID-19 pandemic is spreading rapidly due to the outbreak of novel coronavirus SARS-CoV-2 across the globe. Anti-COVID-19 drugs are urgently required in this situation. In this regard, the discovery of promising new anti-COVID-19 moieties is expected from traditional medicine. The study is aimed to discover phytochemicals of Cocculus hirsutus having anti-COVID-19 activity via inhibiting main proteases of SARS-CoV-2. Main proteases (Mpro) of SARS-CoV-2 serve as a protuberant target for anti-COVID-19 drug discovery because it is a key enzyme of coronaviruses and has a pivotal role in mediating viral replication and transcription that makes it an attractive drug target. Recent studies indicated the utility of C. hirsutus in the treatment of viral disorders like Dengue. Phytochemicals from C. hirsutus were docked against SARS-CoV-2 main proteases (6LU7, 5R7Y, 5R7Z, 5R80, 5R81, 5R82) using the PyRx virtual screen tool and discovery studio visualizer. Further, molecular dynamics simulations were performed (for 100 ns) to see conformational stability for all complexes. Pharmacokinetic properties and drug-likeness prediction of selected C. hirsutus phytoconstituents were also performed. Betulin, coclaurine, and quinic acid of C. hirsutus were found promising with significant binding affinity to SARS-CoV-2 Mpro in comparison to control. They have shown stable interactions with the amino acid residues present on the active site of most of the SARS-CoV-2 Mpro and were found as promising anti-COVID-19 candidates. These compounds could be potential leads for the development of target-specific anti-COVID-19 therapeutics while ethnomedicinal uses of this herb could further needed for its detailed antiviral therapeutic exploration.Communicated by Ramaswamy H. Sarma.

Discovery of novel RARα agonists using pharmacophore-based virtual screening, molecular docking, and molecular dynamics simulation studies

Nuclear retinoic acid receptors (RARs) are ligand-dependent transcription factors involved in various biological processes, such as embryogenesis, cell proliferation, differentiation, reproduction, and apoptosis. These receptors are regulated by retinoids, i.e., retinoic acid (RA) and its analogs, as receptor agonists. RAR agonists are promising therapeutic agents for the treatment of serious dermatological disorders, including some malignant conditions. By inducing apoptosis, they are able to inhibit the proliferation of diverse cancer cell lines. Also, RAR agonists have recently been identified as therapeutic options for some neurodegenerative diseases. These features make retinoids very attractive molecules for medical purposes. Synthetic selective RAR agonists have several advantages over endogenous ones, but they suffer poor pharmacokinetic properties. These compounds are normally lipophilic acids with unfavorable drug-like features such as poor oral bioavailability. Recently, highly selective, potent, and less toxic RAR agonists with proper lipophilicity, thus, good oral bioavailability have been developed for some therapeutic applications. In the present study, ligand and structure-based virtual screening technique was exploited to introduce some novel RARα agonists. Pharmacokinetic assessment was also performed in silico to suggest those compounds which have optimized drug-like features. Finally, two compounds with the best in silico pharmacological features are proposed as lead molecules for future development of RARα agonists.

Computational Guided Drug Targets Identification against Extended-Spectrum Beta-Lactamase-Producing Multi-Drug Resistant Uropathogenic Escherichia coli

Urinary tract infections (UTIs) are one of the most frequent bacterial infections in the world, both in the hospital and community settings. Uropathogenic Escherichia coli (UPEC) are the predominant etiological agents causing UTIs. Extended-spectrum beta-lactamase (ESBL) production is a prominent mechanism of resistance that hinders the antimicrobial treatment of UTIs caused by UPEC and poses a substantial danger to the arsenal of antibiotics now in use. As bacteria have several methods to counteract the effects of antibiotics, identifying new potential drug targets may help in the design of new antimicrobial agents, and in the control of the rising trend of antimicrobial resistance (AMR). The public availability of the entire genome sequences of humans and many disease-causing organisms has accelerated the hunt for viable therapeutic targets. Using a unique, hierarchical, in silico technique using computational tools, we discovered and described potential therapeutic drug targets against the ESBL-producing UPEC strain NA114. Three different sets of proteins (chokepoint, virulence, and resistance genes) were explored in phase 1. In phase 2, proteins shortlisted from phase 1 were analyzed for their essentiality, non-homology to the human genome, and gut flora. In phase 3, the further shortlisted putative drug targets were qualitatively characterized, including their subcellular location, broad-spectrum potential, and druggability evaluations. We found seven distinct targets for the pathogen that showed no similarity to the human proteome. Thus, possibilities for cross-reactivity between a target-specific antibacterial and human proteins were minimized. The subcellular locations of two targets, ECNA114_0085 and ECNA114_1060, were predicted as cytoplasmic and periplasmic, respectively. These proteins play an important role in bacterial peptidoglycan biosynthesis and inositol phosphate metabolism, and can be used in the design of drugs against these bacteria. Inhibition of these proteins will be helpful to combat infections caused by MDR UPEC.

Molecular Descriptors, Structure Generation, and Inverse QSAR/QSPR Based on SELFIES

For inverse QSAR/QSPR in conventional molecular design, several chemical structures must be generated and their molecular descriptors must be calculated. However, there is no one-to-one correspondence between the generated chemical structures and molecular descriptors. In this paper, molecular descriptors, structure generation, and inverse QSAR/QSPR based on self-referencing embedded strings (SELFIES), a 100% robust molecular string representation, are proposed. A one-hot vector is converted from SELFIES to SELFIES descriptors x, and an inverse analysis of the QSAR/QSPR model y = f(x) with the objective variable y and molecular descriptor x is conducted. Thus, x values that achieve a target y value are obtained. Based on these values, SELFIES strings or molecules are generated, meaning that inverse QSAR/QSPR is performed successfully. The SELFIES descriptors and SELFIES-based structure generation are verified using datasets of actual compounds. The successful construction of SELFIES-descriptor-based QSAR/QSPR models with predictive abilities comparable to those of models based on other fingerprints is confirmed. A large number of molecules with one-to-one relationships with the values of the SELFIES descriptors are generated. Furthermore, as a case study of inverse QSAR/QSPR, molecules with target y values are generated successfully. The Python code for the proposed method is available at https://github.com/hkaneko1985/dcekit.

Updates on drug designing approach through computational strategies: a review

The drug discovery and development (DDD) process in pursuit of novel drug candidates is a challenging procedure requiring lots of time and resources. Therefore, computer-aided drug design (CADD) methodologies are used extensively to promote proficiency in drug development in a systematic and time-effective manner. The point in reference is SARS-CoV-2 which has emerged as a global pandemic. In the absence of any confirmed drug moiety to treat the infection, the science fraternity adopted hit and trial methods to come up with a lead drug compound. This article is an overview of the virtual methodologies, which assist in finding novel hits and help in the progression of drug development in a short period with a specific medicinal solution.

Virtual Screening of Nrf2 Dietary-Derived Agonists and Safety by a New Deep-Learning Model and Verified In Vitro and In Vivo

Nuclear factor (erythroid-derived 2)-like 2 (Nrf2) is an essential regulatory target of antioxidants, but the lack of Nrf2 active site information has hindered discovery of new Nrf2 agonists from food-derived compounds by large-scale virtual screening. Two deep-learning models were separately trained to screen for Nrf2-agonists and safety. The trained models screened potentially active chemicals from approximately 70,000 dietary compounds within 5 min. Of the 169 potential Nrf2 agonists identified via deep-learning screening, 137 had not been reported before. Six compounds selected from the new Nrf2 agonists significantly increased (p < 0.05) the activity of Nrf2 on carbon tetrachloride (CCl₄)-intoxicated HepG2 cells (nicotiflorin (99.44 ± 18.5%), artemetin (97.91 ± 8.22%), daidzin (87.73 ± 3.77%), linonin (74.27 ± 5.73%), sinensetin (72.74 ± 10.41%), and tectoridin (77.78 ± 4.80%)), and their safety were demonstrated by an MTT assay. The safety and Nrf2 agonistic activity of nicotiflorin, artemetin, and daidzin were also reconfirm by a single-dose acute oral toxicity study and CCl₄-intoxicated rat assay.

Modes-of-action of antifungal compounds: Stressors and (target-site-specific) toxins, toxicants, or Toxin-stressors

Fungi and antifungal compounds are relevant to the United Nation's Sustainable Development Goals. However, the modes-of-action of antifungals-whether they are naturally occurring substances or anthropogenic fungicides-are often unknown or are misallocated in terms of their mechanistic category. Here, we consider the most effective approaches to identifying whether antifungal substances are cellular stressors, toxins/toxicants (that are target-site-specific), or have a hybrid mode-of-action as Toxin-stressors (that induce cellular stress yet are target-site-specific). This newly described 'toxin-stressor' category includes some photosensitisers that target the cell membrane and, once activated by light or ultraviolet radiation, cause oxidative damage. We provide a glossary of terms and a diagrammatic representation of diverse types of stressors, toxic substances, and Toxin-stressors, a classification that is pertinent to inhibitory substances not only for fungi but for all types of cellular life. A decision-tree approach can also be used to help differentiate toxic substances from cellular stressors (Curr Opin Biotechnol 2015 33: 228-259). For compounds that target specific sites in the cell, we evaluate the relative merits of using metabolite analyses, chemical genetics, chemoproteomics, transcriptomics, and the target-based drug-discovery approach (based on that used in pharmaceutical research), focusing on both ascomycete models and the less-studied basidiomycete fungi. Chemical genetic methods to elucidate modes-of-action currently have limited application for fungi where molecular tools are not yet available; we discuss ways to circumvent this bottleneck. We also discuss ecologically commonplace scenarios in which multiple substances act to limit the functionality of the fungal cell and a number of as-yet-unresolved questions about the modes-of-action of antifungal compounds pertaining to the Sustainable Development Goals.

Artificial Intelligence Technologies in Cardiology

As the world produces exabytes of data, there is a growing need to find new methods that are more suitable for dealing with complex datasets. Artificial intelligence (AI) has significant potential to impact the healthcare industry, which is already on the road to change with the digital transformation of vast quantities of information. The implementation of AI has already achieved success in the domains of molecular chemistry and drug discoveries. The reduction in costs and in the time needed for experiments to predict the pharmacological activities of new molecules is a milestone in science. These successful applications of AI algorithms provide hope for a revolution in healthcare systems. A significant part of artificial intelligence is machine learning (ML), of which there are three main types-supervised learning, unsupervised learning, and reinforcement learning. In this review, the full scope of the AI workflow is presented, with explanations of the most-often-used ML algorithms and descriptions of performance metrics for both regression and classification. A brief introduction to explainable artificial intelligence (XAI) is provided, with examples of technologies that have developed for XAI. We review important AI implementations in cardiology for supervised, unsupervised, and reinforcement learning and natural language processing, emphasizing the used algorithm. Finally, we discuss the need to establish legal, ethical, and methodical requirements for the deployment of AI models in medicine.

MSGNN-DTA: Multi-Scale Topological Feature Fusion Based on Graph Neural Networks for Drug-Target Binding Affinity Prediction

The accurate prediction of drug-target binding affinity (DTA) is an essential step in drug discovery and drug repositioning. Although deep learning methods have been widely adopted for DTA prediction, the complexity of extracting drug and target protein features hampers the accuracy of these predictions. In this study, we propose a novel model for DTA prediction named MSGNN-DTA, which leverages a fused multi-scale topological feature approach based on graph neural networks (GNNs). To address the challenge of accurately extracting drug and target protein features, we introduce a gated skip-connection mechanism during the feature learning process to fuse multi-scale topological features, resulting in information-rich representations of drugs and proteins. Our approach constructs drug atom graphs, motif graphs, and weighted protein graphs to fully extract topological information and provide a comprehensive understanding of underlying molecular interactions from multiple perspectives. Experimental results on two benchmark datasets demonstrate that MSGNN-DTA outperforms the state-of-the-art models in all evaluation metrics, showcasing the effectiveness of the proposed approach. Moreover, the study conducts a case study based on already FDA-approved drugs in the DrugBank dataset to highlight the potential of the MSGNN-DTA framework in identifying drug candidates for specific targets, which could accelerate the process of virtual screening and drug repositioning.

Galaxolide and Irgacure 369 are novel environmental androgens

Endocrine disruptors are environmental chemicals that can interfere with the endocrine system. However, research on endocrine disruptors that interfere with androgen's actions is still limited. The purpose of this study is to use in silico computation, i.e., molecular docking to facilitate the identification of environmental androgens. Computational docking was used to study the binding interactions of environmental/industrial compounds with the three dimensional structure of human androgen receptor (AR). Then reporter assay and cell proliferation assay using AR-expressing LNCaP prostate cancer cells were used to determine their in vitro androgenic activity. Animal studies using immature male rats were also carried out to test their in vivo androgenic activity. Two novel environmental androgens were identified. As a photoinitiator, 2-benzyl-2-(dimethylamino)-4'-morpholinobutyrophenone (Irgacure 369, abbreviated as IC-369) is widely used in the packaging and electronics industries. Galaxolide (HHCB) is widely used in the production of perfume, fabric softeners and detergents. We found that both IC-369 and HHCB could activate AR transcriptional activity and promote cell proliferation in AR-sensitive LNCaP cells. Furthermore, IC-369 and HHCB could induce cell proliferation and histological changes of seminal vesicles in immature rats. RNA sequencing and qPCR analysis showed that androgen-related genes in seminal vesicle tissue were up-regulated by IC-369 and HHCB. In conclusion, IC-369 and HHCB are new environmental androgens that bind AR and induce AR transcriptional activity, thereby exerting toxicological effects on the development of male reproductive organs.

Structure of the Spring Viraemia of Carp Virus Ribonucleoprotein Complex Reveals Its Assembly Mechanism and Application in Antiviral Drug Screening

Spring viremia of carp virus (SVCV) is a highly pathogenic Vesiculovirus infecting the common carp, yet neither a vaccine nor effective therapies are available to treat spring viremia of carp (SVC). Like all negative-sense viruses, SVCV contains an RNA genome that is encapsidated by the nucleoprotein (N) in the form of a ribonucleoprotein (RNP) complex, which serves as the template for viral replication and transcription. Here, the three-dimensional (3D) structure of SVCV RNP was resolved through cryo-electron microscopy (cryo-EM) at a resolution of 3.7 Å. RNP assembly was stabilized by N and C loops; RNA was wrapped in the groove between the N and C lobes with 9 nt nucleotide per protomer. Combined with mutational analysis, our results elucidated the mechanism of RNP formation. The RNA binding groove of SVCV N was used as a target for drug virtual screening, and it was found suramin had a good antiviral effect. This study provided insights into RNP assembly, and anti-SVCV drug screening was performed on the basis of this structure, providing a theoretical basis and efficient drug screening method for the prevention and treatment of SVC. IMPORTANCE Aquaculture accounts for about 70% of global aquatic products, and viral diseases severely harm the development of aquaculture industry. Spring viremia of carp virus (SVCV) is the pathogen causing highly contagious spring viremia of carp (SVC) disease in cyprinids, especially common carp (Cyprinus carpio), yet neither a vaccine nor effective therapies are available to treat this disease. In this study, we have elucidated the mechanism of SVCV ribonucleoprotein complex (RNP) formation by resolving the 3D structure of SVCV RNP and screened antiviral drugs based on the structure. It is found that suramin could competitively bind to the RNA binding groove and has good antiviral effects both in vivo and in vitro. Our study provides a template for rational drug discovery efforts to treat and prevent SVCV infections.

Antibiotic discovery against Piscirickettsia salmonis using a combined in silico and in vitro approach

Piscirickettsia salmonis is one of the main pathogens causing considerable economic losses in salmonid farming. The DNA gyrase of several pathogenic bacteria has been the target of choice for antibiotic design and discovery for years, due to its key function during DNA replication. In this study, we carried out a combined in silico and in vitro approach to antibiotic discovery targeting the GyrA subunit of Piscirickettsia salmonis. The in silico results of this work showed that flumequine (-6.6 kcal/mol), finafloxacin (-7.2 kcal/mol), rosoxacin (-6.6 kcal/mol), elvitegravir (-6.4 kcal/mol), sarafloxacin (-8.3 kcal/mol), orbifloxacin (-7.9 kcal/mol), and sparfloxacin (-7.2 kcal/mol) are docked with good affinities in the DNA binding domain of the Piscirickettsia salmonis GyrA subunit. In the in vitro inhibition assay, it was observed that most of these molecules inhibit the growth of Piscirickettsia salmonis, except for elvitegravir. We believe that this methodology could help to significantly reduce the time and cost of antibiotic discovery trials to combat Piscirickettsia salmonis within the salmonid farming industry.

Multitargeted Virtual Screening and Molecular Simulation of Natural Product-like Compounds against GSK3β, NMDA-Receptor, and BACE-1 for the Management of Alzheimer's Disease

The complexity of Alzheimer's disease (AD) and several side effects of currently available medication inclined us to search for a novel natural cure by targeting multiple key regulatory proteins. We initially virtually screened the natural product-like compounds against GSK3β, NMDA receptor, and BACE-1 and thereafter validated the best hit through molecular dynamics simulation (MDS). The results demonstrated that out of 2029 compounds, only 51 compounds exhibited better binding interactions than native ligands, with all three protein targets (NMDA, GSK3β, and BACE) considered multitarget inhibitors. Among them, F1094-0201 is the most potent inhibitor against multiple targets with binding energy -11.7, -10.6, and -12 kcal/mol, respectively. ADME-T analysis results showed that F1094-0201 was found to be suitable for CNS drug-likeness in addition to their other drug-likeness properties. The MDS results of RMSD, RMSF, Rg, SASA, SSE and residue interactions indicated the formation of a strong and stable association in the complex of ligands (F1094-0201) and proteins. These findings confirm the F1094-0201's ability to remain inside target proteins' binding pockets while forming a stable complex of protein-ligand. The free energies (MM/GBSA) of BACE-F1094-0201, GSK3β-F1094-0201, and NMDA-F1094-0201 complex formation were -73.78 ± 4.31 kcal mol^-1, -72.77 ± 3.43 kcal mol^-1, and -52.51 ± 2.85 kcal mol^-1, respectively. Amongst the target proteins, F1094-0201 have a more stable association with BACE, followed by NMDA and GSK3β. These attributes of F1094-0201 indicate it as a possible option for the management of pathophysiological pathways associated with AD.

Molecular mechanisms augmenting resistance to current therapies in clinics among cervical cancer patients

Cervical cancer (CC) is the fourth leading cause of cancer death (~ 324,000 deaths annually) among women internationally, with 85% of these deaths reported in developing regions, particularly sub-Saharan Africa and Southeast Asia. Human papillomavirus (HPV) is considered the major driver of CC, and with the availability of the prophylactic vaccine, HPV-associated CC is expected to be eliminated soon. However, female patients with advanced-stage cervical cancer demonstrated a high recurrence rate (50-70%) within two years of completing radiochemotherapy. Currently, 90% of failures in chemotherapy are during the invasion and metastasis of cancers related to drug resistance. Although molecular target therapies have shown promising results in the lab, they have had little success in patients due to the tumor heterogeneity fueling resistance to these therapies and bypass the targeted signaling pathway. The last two decades have seen the emergence of immunotherapy, especially immune checkpoint blockade (ICB) therapies, as an effective treatment against metastatic tumors. Unfortunately, only a small subgroup of patients (< 20%) have benefited from this approach, reflecting disease heterogeneity and manifestation with primary or acquired resistance over time. Thus, understanding the mechanisms driving drug resistance in CC could significantly improve the quality of medical care for cancer patients and steer them to accurate, individualized treatment. The rise of artificial intelligence and machine learning has also been a pivotal factor in cancer drug discovery. With the advancement in such technology, cervical cancer screening and diagnosis are expected to become easier. This review will systematically discuss the different tumor-intrinsic and extrinsic mechanisms CC cells to adapt to resist current treatments and scheme novel strategies to overcome cancer drug resistance.

QSAR and molecular docking studies of isatin and indole derivatives as SARS 3CLpro inhibitors

The 3C-like protease (3CL^pro), known as the main protease of SARS-COV, plays a vital role in the viral replication cycle and is a critical target for the development of SARS inhibitor. Comparative sequence analysis has shown that the 3CL^pro of two coronaviruses, SARS-CoV-2 and SARS-CoV, show high structural similarity, and several common features are shared among the substrates of 3CL^pro in different coronaviruses. The goal of this study is the development of validated QSAR models by CORAL software and Monte Carlo optimization to predict the inhibitory activity of 81 isatin and indole-based compounds against SARS CoV 3CL^pro. The models were built using a newer objective function optimization of this software, known as the index of ideality correlation (IIC), which provides favorable results. The entire set of molecules was randomly divided into four sets including: active training, passive training, calibration and validation sets. The optimal descriptors were selected from the hybrid model by combining SMILES and hydrogen suppressed graph (HSG) based on the objective function. According to the model interpretation results, eight synthesized compounds were extracted and introduced from the ChEMBL database as good SARS CoV 3CL^pro inhibitor. Also, the activity of the introduced molecules further was supported by docking studies using 3CL^pro of both SARS-COV-1 and SARS-COV-2. Based on the results of ADMET and OPE study, compounds CHEMBL4458417 and CHEMBL4565907 both containing an indole scaffold with the positive values of drug-likeness and the highest drug-score can be introduced as selected leads.

In-silico assay of a dosing vehicle based on chitosan-TiO2 and modified benzofuran-isatin molecules against Pseudomonas aeruginosa

A high priority of the World Health Organization (WHO) is the study of drugs against Pseudomonas aeruginosa, which has developed antibiotic resistance. In this order, recent research is analyzing biomaterials and metal oxide nanoparticles, such as chitosan (QT) and TiO2 (NT), which can transport molecules with biological activity against bacteria, to propose them as drug carrier candidates. In the present work, 10 modified benzofuran-isatin molecules were studied through computational simulation using density functional theory (DFT) and molecular docking assays against Hfq and LpxC (proteins of P. aeruginosa). The results show that the ligand efficiency of commercial drugs C-CP and C-AZI against Hfq is low compared with the best-designed molecule MOL-A. However, we highlight that the influence of NT promotes a better interaction of some molecules, where MOL-E generates a better interaction by 0.219 kcal/mol when NT is introduced in Hfq, forming the system Hfq-NT (Target-NT). Similar behavior is observed in the LpxC target, in which MOL-J is better at 0.072 kcal/mol. Finally, two pharmacophoric models for Hfq and LpxC implicate hydrophobic and aromatic-hydrophobic fragments.

A New and an Eco-Friendly Approach for the Construction of Multi-Functionalized Benzenes with Computational Studies

The new path chosen is more appropriate in the context of green chemistry. This research aims to construct 5,6,7,8-tetrahydronaphthalene-1,3-dicarbonitrile (THNDC) and 1,2,3,4-tetrahydroisoquinoline-6,8-dicarbonitrile (THIDC) derivatives via the cyclization of three easily obtainable reactants under an environmentally benign mortar and pestle grinding technique. Notably, the robust route offers an esteemed opportunity for the introduction of multi-substituted benzenes and ensures the good compatibility of bioactive molecules. Furthermore, the synthesized compounds are investigated using docking simulations with two representative drugs (6c and 6e) for target validation. The physicochemical, pharmacokinetic, drug-like properties (ADMET), and therapeutic friendliness characteristics of these synthesized compounds are computed.

An overview of compound properties, multiparameter optimization, and computational drug design methods for PARP-1 inhibitor drugs

Breast cancer treatment with PARP-1 inhibitors remains challenging due to emerging toxicities, drug resistance, and unaffordable costs of treatment options. How do we invent strategies to design better anti-cancer drugs? A part of the answer is in optimized compound properties, desirability functions, and modern computational drug design methods that drive selectivity and toxicity and have not been reviewed for PARP-1 inhibitors. Nonetheless, comparisons of these compound properties for PARP-1 inhibitors are not available in the literature. In this review, we analyze the physchem, PKPD space to identify inherent desirability functions characteristic of approved drugs that can be valuable for the design of better candidates. Recent literature utilizing ligand, structure-based drug design strategies and matched molecular pair analysis (MMPA) for the discovery of novel PARP-1 inhibitors are also reviewed. Thus, this perspective provides valuable insights into the medchem and multiparameter optimization of PARP-1 inhibitors that might be useful to other medicinal chemists.

Docking and Molecular Dynamics Identify Leads against 5 Alpha Reductase 2 for Benign Prostate Hyperplasia Treatment

Steroid 5 alpha-reductase 2 (5αR-2) is a membrane-embedded protein that together with other isoforms plays a key role in the metabolism of steroids. This enzyme catalyzes the reduction of testosterone to the more potent ligand, dihydrotestosterone (DHT) in the prostate. Androgens, testosterone, and DHT play important roles in prostate growth, development, and function. At the same time, both testosterone and DHT have been implicated in the pathogenesis of benign prostate hyperplasia (BPH). Inhibition of the DHT formation, therefore, provides a therapeutic strategy that offers the possibility of preventing, delaying, or treating BPH. Currently, two steroidal drugs that inhibit 5αR-2, dutasteride and finasteride, have been approved for clinical use. These two come at a high cost and also portray undesirable sexual side effects which necessitate the need to find new chemotherapeutic alternatives for the disease. Based on the aforementioned, finasteride and dutasteride were subjected to scaffold hopping, fragment-based de novo design, molecular docking, and molecular dynamics simulations employing databases like ChEMBL, DrugBank, PubChem, ChemSpider, and Zinc15 in the identification of potential hits targeting 5αR-2. Altogether, ten novel compounds targeting 5αR-2 were identified with binding energies lower or comparable to finasteride and dutasteride, the main inhibitors for this target. Molecular docking and molecular dynamics simulations studies identify amino acid residues Glu57, Phe219, Phe223, and Leu224 to be critical for ligand binding and complex stability. The physicochemical and pharmacological profiling suggests the potential of the hit compounds to be drug-like and orally active. Similarly, the quality parameter assessments revealed the hits possess LELP greater than 3 implying their promise as lead-like molecules. The compounds A5, A9, and A10 were, respectively, predicted to treat prostate disorders with Pa (0.188, 0.361, and 0.270) and Pi (0.176, 0.050, and 0.093), while A8 and A9 were found to be associated with BPH treatment with Pa (0.09 and 0.127) and Pi (0.077 and 0.033), respectively. Structural similarity searches via DrugBank identified the drugs faropenem, acemetacin, estradiol valerate, and yohimbine to be useful for BPH treatment suggesting the de novo designed ligands as potential chemotherapeutic agents for treating this disease.

Discovery of novel S6K1 inhibitors by an ensemble-based virtual screening method and molecular dynamics simulation

Ribosomal protein S6 kinase beta-1 (S6K1) is considered a potential target for the treatment of various diseases, such as obesity, type II diabetes, and cancer. Development of novel S6K1 inhibitors is an urgent and important task for the medicinal chemists. In this research, an effective ensemble-based virtual screening method, including common feature pharmacophore model, 3D-QSAR pharmacophore model, naïve Bayes classifier model, and molecular docking, was applied to discover potential S6K1 inhibitors from BioDiversity database with 29,158 compounds. Finally, 7 hits displayed considerable properties and considered as potential inhibitors against S6K1. Further, carefully analyzing the interactions between these 7 hits and key residues in the S6K1 active site, and comparing them with the reference compound PF-4708671, it was found that 2 hits exhibited better binding patterns. In order to further investigate the mechanism of the interactions between 2 hits and S6K1 at simulated physiological conditions, the molecular dynamics simulation was performed. The ΔG_bind energies for S6K1-Hit1 and S6K1-Hit2 were - 111.47 ± 1.29 and - 54.29 ± 1.19 kJ mol^-1, respectively. Furthermore, deep analysis of these results revealed that Hit1 was the most stable complex, which can stably bind to S6K1 active site, interact with all of the key residues, and induce H1, H2, and M-loop regions changes. Therefore, the identified Hit1 may be a promising lead compound for developing new S6K1 inhibitor for various metabolic diseases treatment.

Developing novel antifungals: lessons from G protein-coupled receptors

Up to 1.5 million people die yearly from fungal disease, but the repertoire of antifungal drug classes is minimal and the incidence of drug resistance is rising rapidly. This dilemma was recently declared by the World Health Organization as a global health emergency, but the discovery of new antifungal drug classes remains excruciatingly slow. This process could be accelerated by focusing on novel targets, such as G protein-coupled receptor (GPCR)-like proteins, that have a high likelihood of being druggable and have well-defined biology and roles in disease. We discuss recent successes in understanding the biology of virulence and in structure determination of yeast GPCRs, and highlight new approaches that might pay significant dividends in the urgent search for novel antifungal drugs.

Co-creation environment with cloud virtual reality and real-time artificial intelligence toward the design of molecular robots

This paper describes the design philosophy for our cloud-based virtual reality (VR) co-creation environment (CCE) for molecular modeling. Using interactive VR simulation can provide enhanced perspectives in molecular modeling for intuitive live demonstration and experimentation in the CCE. Then the use of the CCE can enhance knowledge creation by bringing people together to share and create ideas or knowledge that may not emerge otherwise. Our prototype CCE discussed here, which was developed to demonstrate our design philosophy, has already enabled multiple members to log in and touch virtual molecules running on a cloud server with no noticeable network latency via real-time artificial intelligence techniques. The CCE plays an essential role in the rational design of molecular robot parts, which consist of bio-molecules such as DNA and protein molecules.

Computational Chemistry for the Identification of Lead Compounds for Radiotracer Development

The use of computer-aided drug design (CADD) for the identification of lead compounds in radiotracer development is steadily increasing. Traditional CADD methods, such as structure-based and ligand-based virtual screening and optimization, have been successfully utilized in many drug discovery programs and are highlighted throughout this review. First, we discuss the use of virtual screening for hit identification at the beginning of drug discovery programs. This is followed by an analysis of how the hits derived from virtual screening can be filtered and culled to highly probable candidates to test in in vitro assays. We then illustrate how CADD can be used to optimize the potency of experimentally validated hit compounds from virtual screening for use in positron emission tomography (PET). Finally, we conclude with a survey of the newest techniques in CADD employing machine learning (ML).

An extensive survey on the use of supervised machine learning techniques in the past two decades for prediction of drug side effects

Approved drugs for sale must be effective and safe, implying that the drug's advantages outweigh its known harmful side effects. Side effects (SE) of drugs are one of the common reasons for drug failure that may halt the whole drug discovery pipeline. The side effects might vary from minor concerns like a runny nose to potentially life-threatening issues like liver damage, heart attack, and death. Therefore, predicting the side effects of the drug is vital in drug development, discovery, and design. Supervised machine learning-based side effects prediction task has recently received much attention since it reduces time, chemical waste, design complexity, risk of failure, and cost. The advancement of supervised learning approaches for predicting side effects have emerged as essential computational tools. Supervised machine learning technique provides early information on drug side effects to develop an effective drug based on drug properties. Still, there are several challenges to predicting drug side effects. Thus, a near-exhaustive survey is carried out in this paper on the use of supervised machine learning approaches employed in drug side effects prediction tasks in the past two decades. In addition, this paper also summarized the drug descriptor required for the side effects prediction task, commonly utilized drug properties sources, computational models, and their performances. Finally, the research gap, open problems, and challenges for the further supervised learning-based side effects prediction task have been discussed.