Drug Design by Pharmacophore and Virtual Screening Approach

Discovery of selective HDAC6 inhibitors based on a multi-layer virtual screening strategy

The abnormal enhancement of histone deacetylase 6 (HDAC6) has been demonstrated to be closely related to the occurrence and development of various malignant tumors, attracting extensive attention as a promising target for cancer therapy. Currently, only limited selective HDAC6 inhibitors have entered clinical trials, making the rapid discovery of selective HDAC6 inhibitors with safety profiles particularly urgent. In this study, a multi-layer virtual screening workflow was established, and the representative compounds screened were biologically evaluated in combination with enzyme inhibitory and anti-tumor cell proliferation experiments. The experimental results showed that the screened compounds L-25, L-32, L-45 and L-81 exhibited nanomolar inhibitory activity against HDAC6, and exerted a certain degree of anti-proliferative activities against tumor cells, especially the cytotoxicity of L-45 to A375 (IC₅₀ = 11.23 ± 1.27 μM) and the cytotoxicity of L-81 against HCT-116 (IC₅₀ = 12.25 ± 1.13 μM). Additionally, the molecular mechanisms underlying the subtype selective inhibitory activities of the selected compounds were further elucidated using computational approaches, and the hotspot residues on HDAC6 contributing to the ligands' binding were identified. In summary, this study established a multi-layer screening scheme to quickly and effectively screen out hit compounds with enzyme inhibitory activity and anti-tumor cell proliferation, providing novel scaffolds for the subsequent anti-tumor drug design based on HDAC6 target.

Developing new PI3Kγ inhibitors by combining pharmacophore modeling, molecular dynamic simulation, molecular docking, fragment-based drug design, and virtual screening

Since dysregulation of the phosphatidylinositol 3-kinase gamma (PI3Kγ) signaling pathway is associated with the pathogenesis of cancer, inflammation, and autoimmunity, PI3Kγ has emerged as an attractive target for drug development. IPI-549 is the only selective PI3Kγ inhibitor that has advanced to clinical trials, thus, IPI-549 could serve as a promising template for designing novel PI3Kγ inhibitors. In this present study, a modeling strategy consisting of common feature pharmacophore modeling, receptor-ligand pharmacophore modeling, and molecular dynamics simulation was utilized to identify the key pharmacodynamic characteristic elements of the target compound and the key residue information of the PI3Kγ interaction with the inhibitors. Then, 10 molecules were designed based on the structure-activity relationships, and some of them exhibited satisfactory predicted binding affinities to PI3Kγ. Finally, a hierarchical multistage virtual screening method, involving the developed common feature and receptor-ligand pharmacophore model and molecular docking, was constructed for screening the potential PI3Kγ inhibitors. Overall, we hope these findings would provide some guidance for the development of novel PI3Kγ inhibitors.

Structure-based identification of potential substrate antagonists for isethionate sulfite-lyase enzyme of Bilophila Wadsworthia: Towards novel therapeutic intervention to curb gut-associated illness

Bilophila wadsworthia is one of the prominent sources of hydrogen sulfide (H₂S) production in appendices, excessive levels of which can result in a weaker colonic mucus barrier, inflammatory bowel disease, and colorectal cancer. Isethionate sulfite-lyase (IslA) enzyme catalyzes H2S production by cleaving CS bond in isethionate, producing acetaldehyde and sulfite. In this study, we aimed to identify potential substrate antagonists for IsIA using a structure-based drug design. Initially, pharmacophore-based computational screening of the ZINC20 database yielded 66 hits that were subjected to molecular docking targeting the isethionate binding site of IsIA. Based on striking docking scores, nine compounds showed strong interaction with critical IsIA residues (Arg189, Gln193, Glu470, Cys468, and Arg678), drug-like features, appropriate adsorption, metabolism, excretion, and excretion profile with non-toxicity. Molecular dynamics simulations uncovered the significant impact of binding the compounds on protein conformational dynamics. Finally, binding free energies revealed substantial binding affinity (ranging from -35.23 to -53.88 kcal/mol) of compounds (ZINC913876497, ZINC913856647, ZINC914263733, ZINC914137795, ZINC915757996, ZINC914357083, ZINC913934833, ZINC9143362047, and ZINC913854740) for IsIA. The compounds proposed herein through a multi-faceted computational strategy can be experimentally validated as potential substrate antagonists of B. wadsworthia's IsIA for developing new medications to curb gut-associated illness in the future.

In Silico Methodologies to Improve Antioxidants' Characterization from Marine Organisms

Marine organisms have been reported to be valuable sources of bioactive molecules that have found applications in different industrial fields. From organism sampling to the identification and bioactivity characterization of a specific compound, different steps are necessary, which are time- and cost-consuming. Thanks to the advent of the -omic era, numerous genome, metagenome, transcriptome, metatranscriptome, proteome and microbiome data have been reported and deposited in public databases. These advancements have been fundamental for the development of in silico strategies for basic and applied research. In silico studies represent a convenient and efficient approach to the bioactivity prediction of known and newly identified marine molecules, reducing the time and costs of "wet-lab" experiments. This review focuses on in silico approaches applied to bioactive molecule discoveries from marine organisms. When available, validation studies reporting a bioactivity assay to confirm the presence of an antioxidant molecule or enzyme are reported, as well. Overall, this review suggests that in silico approaches can offer a valuable alternative to most expensive approaches and proposes them as a little explored field in which to invest.

Advances in the discovery of new chemotypes through ultra-large library docking

INTRODUCTION

The size and complexity of virtual screening libraries in drug discovery have skyrocketed in recent years, reaching up to multiple billions of accessible compounds. However, virtual screening of such ultra-large libraries poses several challenges associated with preparing the libraries, sampling, and pre-selection of suitable compounds. The utilization of artificial intelligence (AI)-assisted screening approaches, such as deep learning, poses a promising countermeasure to deal with this rapidly expanding chemical space. For example, various AI-driven methods were recently successfully used to identify novel small molecule inhibitors of the SARS-CoV-2 main protease (M^pro).

AREAS COVERED

This review focuses on presenting various kinds of virtual screening methods suitable for dealing with ultra-large libraries. Challenges associated with these computational methodologies are discussed, and recent advances are highlighted in the example of the discovery of novel M^pro inhibitors targeting the SARS-CoV-2 virus.

EXPERT OPINION

With the rapid expansion of the virtual chemical space, the methodologies for docking and screening such quantities of molecules need to keep pace. Employment of AI-driven screening compounds has already been shown to be effective in a range from a few thousand to multiple billion compounds, furthered by de novo generation of drug-like molecules without human interference.

Computational Approaches to the Rational Design of Tubulin-Targeting Agents

Microtubules are highly dynamic polymers of α,β-tubulin dimers which play an essential role in numerous cellular processes such as cell proliferation and intracellular transport, making them an attractive target for cancer and neurodegeneration research. To date, a large number of known tubulin binders were derived from natural products, while only one was developed by rational structure-based drug design. Several of these tubulin binders show promising in vitro profiles while presenting unacceptable off-target effects when tested in patients. Therefore, there is a continuing demand for the discovery of safer and more efficient tubulin-targeting agents. Since tubulin structural data is readily available, the employment of computer-aided design techniques can be a key element to focus on the relevant chemical space and guide the design process. Due to the high diversity and quantity of structural data available, we compiled here a guide to the accessible tubulin-ligand structures. Furthermore, we review different ligand and structure-based methods recently used for the successful selection and design of new tubulin-targeting agents.

Computer-Aided Drug Design towards New Psychotropic and Neurological Drugs

Central nervous system (CNS) disorders are a therapeutic area in drug discovery where demand for new treatments greatly exceeds approved treatment options. This is complicated by the high failure rate in late-stage clinical trials, resulting in exorbitant costs associated with bringing new CNS drugs to market. Computer-aided drug design (CADD) techniques minimise the time and cost burdens associated with drug research and development by ensuring an advantageous starting point for pre-clinical and clinical assessments. The key elements of CADD are divided into ligand-based and structure-based methods. Ligand-based methods encompass techniques including pharmacophore modelling and quantitative structure activity relationships (QSARs), which use the relationship between biological activity and chemical structure to ascertain suitable lead molecules. In contrast, structure-based methods use information about the binding site architecture from an established protein structure to select suitable molecules for further investigation. In recent years, deep learning techniques have been applied in drug design and present an exciting addition to CADD workflows. Despite the difficulties associated with CNS drug discovery, advances towards new pharmaceutical treatments continue to be made, and CADD has supported these findings. This review explores various CADD techniques and discusses applications in CNS drug discovery from 2018 to November 2022.

Computational design of prospective molecular targets for Burkholderia cepacia complex by molecular docking and dynamic simulation studies

The study aimed to screen prospective molecular targets of BCC and potential natural lead candidates as effective binders by computational modeling, molecular docking, and dynamic (MD) simulation studies. Based on the virulent functions, tRNA 5-methylaminomethyl-2-thiouridine biosynthesis bifunctional protein (mnmC) and pyrimidine/purine nucleoside phosphorylase (ppnP) were selected as the prospective molecular targets. In the absence of experimental data, the three-dimensional (3D) structures of these targets were computationally predicted. After a thorough literature survey and database search, the drug-likeness, and pharmacokinetic properties of 70 natural molecules were computationally predicted and the effectual binding of the best lead molecules against both the targets was predicted by molecular docking. The stabilities of the best-docked complexes were validated by MD simulation and the binding energy calculations were carried out by MM-GBSA approaches. The present study revealed that the hypothetical models of mnmC and ppnP showed stereochemical accuracy. The study also showed that among 70 natural compounds subjected to computational screening, Honokiol (3',5-Di(prop-2-en-1-yl) [1,1'-biphenyl]-2,4'-diol) present in Magnolia showed ideal drug-likeness, pharmacokinetic features and showed effectual binding with mnmC and ppnP (binding energies -7.3 kcal/mol and -6.6 kcal/mol, respectively). The MD simulation and GBSA calculation studies showed that the ligand-protein complexes stabilized throughout tMD simulation. The present study suggests that Honokiol can be used as a potential lead molecule against mnmC and ppnP targets of BCC and this study provides insight into further experimental validation for alternative lead development against drug resistant BCC.

An in-silico pharmacophore-based molecular docking study to evaluate the inhibitory potentials of novel fungal triterpenoid Astrakurkurone analogues against a hypothetical mutated main protease of SARS-CoV-2 virus

BACKGROUND

The main protease is an important structural protein of SARS-CoV-2, essential for its survivability inside a human host. Considering current vaccines' limitations and the absence of approved therapeutic targets, Mpro may be regarded as the potential candidate drug target. Novel fungal phytocompound Astrakurkurone may be studied as the potential Mpro inhibitor, considering its medicinal properties reported elsewhere.

METHODS

In silico molecular docking was performed with Astrakurkurone and its twenty pharmacophore-based analogues against the native Mpro protein. A hypothetical Mpro was also constructed with seven mutations and targeted by Astrakurkurone and its analogues. Furthermore, multiple parameters such as statistical analysis (Principal Component Analysis), pharmacophore alignment, and drug likeness evaluation were performed to understand the mechanism of protein-ligand molecular interaction. Finally, molecular dynamic simulation was done for the top-ranking ligands to validate the result.

RESULT

We identified twenty Astrakurkurone analogues through pharmacophore screening methodology. Among these twenty compounds, two analogues namely, ZINC89341287 and ZINC12128321 showed the highest inhibitory potentials against native and our hypothetical mutant Mpro, respectively (-7.7 and -7.3 kcal mol-1) when compared with the control drug Telaprevir (-5.9 and -6.0 kcal mol-1). Finally, we observed that functional groups of ligands namely two aromatic and one acceptor groups were responsible for the residual interaction with the target proteins. The molecular dynamic simulation further revealed that these compounds could make a stable complex with their respective protein targets in the near-native physiological condition.

CONCLUSION

To conclude, Astrakurkurone analogues ZINC89341287 and ZINC12128321 can be potential therapeutic agents against the highly infectious SARS-CoV-2 virus.

Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: II. regression and dielectric constant

End-point free energy calculations as a powerful tool have been widely applied in protein-ligand and protein-protein interactions. It is often recognized that these end-point techniques serve as an option of intermediate accuracy and computational cost compared with more rigorous statistical mechanic models (e.g., alchemical transformation) and coarser molecular docking. However, it is observed that this intermediate level of accuracy does not hold in relatively simple and prototypical host-guest systems. Specifically, in our previous work investigating a set of carboxylated-pillar[6]arene host-guest complexes, end-point methods provide free energy estimates deviating significantly from the experimental reference, and the rank of binding affinities is also incorrectly computed. These observations suggest the unsuitability and inapplicability of standard end-point free energy techniques in host-guest systems, and alteration and development are required to make them practically usable. In this work, we consider two ways to improve the performance of end-point techniques. The first one is the PBSA_E regression that varies the weights of different free energy terms in the end-point calculation procedure, while the second one is considering the interior dielectric constant as an additional variable in the end-point equation. By detailed investigation of the calculation procedure and the simulation outcome, we prove that these two treatments (i.e., regression and dielectric constant) are manipulating the end-point equation in a somehow similar way, i.e., weakening the electrostatic contribution and strengthening the non-polar terms, although there are still many detailed differences between these two methods. With the trained end-point scheme, the RMSE of the computed affinities is improved from the standard ~ 12 kcal/mol to ~ 2.4 kcal/mol, which is comparable to another altered end-point method (ELIE) trained with system-specific data. By tuning PBSA_E weighting factors with the host-specific data, it is possible to further decrease the prediction error to ~ 2.1 kcal/mol. These observations along with the extremely efficient optimized-structure computation procedure suggest the regression (i.e., PBSA_E as well as its GBSA_E extension) as a practically applicable solution that brings end-point methods back into the library of usable tools for host-guest binding. However, the dielectric-constant-variable scheme cannot effectively minimize the experiment-calculation discrepancy for absolute binding affinities, but is able to improve the calculation of affinity ranks. This phenomenon is somehow different from the protein-ligand case and suggests the difference between host-guest and biomacromolecular (protein-ligand and protein-protein) systems. Therefore, the spectrum of tools usable for protein-ligand complexes could be unsuitable for host-guest binding, and numerical validations are necessary to screen out really workable solutions in these 'prototypical' situations.

Identification of Active Compounds against Melanoma Growth by Virtual Screening for Non-Classical Human DHFR Inhibitors

Antifolates such as methotrexate (MTX) have been largely known as anticancer agents because of their role in blocking nucleic acid synthesis and cell proliferation. Their mechanism of action lies in their ability to inhibit enzymes involved in the folic acid cycle, especially human dihydrofolate reductase (hDHFR). However, most of them have a classical structure that has proven ineffective against melanoma, and, therefore, inhibitors with a non-classical lipophilic structure are increasingly becoming an attractive alternative to circumvent this clinical resistance. In this study, we conducted a protocol combining virtual screening (VS) and cell-based assays to identify new potential non-classical hDHFR inhibitors. Among 173 hit compounds identified (average logP = 3.68; average MW = 378.34 Da), two-herein, called C1 and C2-exhibited activity against melanoma cell lines B16 and A375 by MTT and Trypan-Blue assays. C1 showed cell growth arrest (39% and 56%) and C2 showed potent cytotoxic activity (77% and 51%) in a dose-dependent manner. The effects of C2 on A375 cell viability were greater than MTX (98% vs 60%) at equivalent concentrations and times. Our results indicate that the integrated in silico/in vitro approach provided a benchmark to identify novel promising non-classical DHFR inhibitors showing activity against melanoma cells.

Discovery of Novel HIV Protease Inhibitors Using Modern Computational Techniques

The human immunodeficiency virus type 1 (HIV-1) has continued to be a global concern. With the new HIV incidence, the emergence of multi-drug resistance and the untoward side effects of currently used anti-HIV drugs, there is an urgent need to discover more efficient anti-HIV drugs. Modern computational tools have played vital roles in facilitating the drug discovery process. This research focuses on a pharmacophore-based similarity search to screen 111,566,735 unique compounds in the PubChem database to discover novel HIV-1 protease inhibitors (PIs). We used an in silico approach involving a 3D-similarity search, physicochemical and ADMET evaluations, HIV protease-inhibitor prediction (IC50/percent inhibition), rigid receptor-molecular docking studies, binding free energy calculations and molecular dynamics (MD) simulations. The 10 FDA-approved HIV PIs (saquinavir, lopinavir, ritonavir, amprenavir, fosamprenavir, atazanavir, nelfinavir, darunavir, tipranavir and indinavir) were used as reference. The in silico analysis revealed that fourteen out of the twenty-eight selected optimized hit molecules were within the acceptable range of all the parameters investigated. The hit molecules demonstrated significant binding affinity to the HIV protease (PR) when compared to the reference drugs. The important amino acid residues involved in hydrogen bonding and п-п stacked interactions include ASP25, GLY27, ASP29, ASP30 and ILE50. These interactions help to stabilize the optimized hit molecules in the active binding site of the HIV-1 PR (PDB ID: 2Q5K). HPS/002 and HPS/004 have been found to be most promising in terms of IC50/percent inhibition (90.15%) of HIV-1 PR, in addition to their drug metabolism and safety profile. These hit candidates should be investigated further as possible HIV-1 PIs with improved efficacy and low toxicity through in vitro experiments and clinical trial investigations.

Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: I. Standard procedure

Despite the massive application of end-point free energy methods in protein-ligand and protein-protein interactions, computational understandings about their performance in relatively simple and prototypical host-guest systems are limited. In this work, we present a comprehensive benchmark calculation with standard end-point free energy techniques in a recent host-guest dataset containing 13 host-guest pairs involving the carboxylated-pillar[6]arene host. We first assess the charge schemes for solutes by comparing the charge-produced electrostatics with many ab initio references, in order to obtain a preliminary albeit detailed view of the charge quality. Then, we focus on four modelling details of end-point free energy calculations, including the docking procedure for the generation of initial condition, the charge scheme for host and guest molecules, the water model used in explicit-solvent sampling, and the end-point methods for free energy estimation. The binding thermodynamics obtained with different modelling schemes are compared with experimental references, and some practical guidelines on maximizing the performance of end-point methods in practical host-guest systems are summarized. Further, we compare our simulation outcome with predictions in the grand challenge and discuss further developments to improve the prediction quality of end-point free energy methods. Overall, unlike the widely acknowledged applicability in protein-ligand binding, the standard end-point calculations cannot produce useful outcomes in host-guest binding and thus are not recommended unless alterations are performed.

Small Molecules as Toll-like Receptor 4 Modulators Drug and In-House Computational Repurposing

The innate immunity toll-like receptor 4 (TLR4) system is a receptor of paramount importance as a therapeutic target. Virtual screening following a “computer-aided drug repurposing” approach was applied to the discovery of novel TLR4 modulators with a non-lipopolysaccharide-like structure. We screened almost 29,000 approved drugs and drug-like molecules from commercial, public, and in-house academia chemical libraries and, after biological assays, identified several compounds with TLR4 antagonist activity. Our computational protocol showed to be a robust approach for the identification of hits with drug-like scaffolds as possible inhibitors of the TLR4 innate immune pathways. Our collaborative work broadens the chemical diversity for inspiration of new classes of TLR4 modulators.

Application of Mathematical Modeling and Computational Tools in the Modern Drug Design and Development Process

The conventional drug discovery approach is an expensive and time-consuming process, but its limitations have been overcome with the help of mathematical modeling and computational drug design approaches. Previously, finding a small molecular candidate as a drug against a disease was very costly and required a long time to screen a compound against a specific target. The development of novel targets and small molecular candidates against different diseases including emerging and reemerging diseases remains a major concern and necessitates the development of novel therapeutic targets as well as drug candidates as early as possible. In this regard, computational and mathematical modeling approaches for drug development are advantageous due to their fastest predictive ability and cost-effectiveness features. Computer-aided drug design (CADD) techniques utilize different computer programs as well as mathematics formulas to comprehend the interaction of a target and drugs. Traditional methods to determine small-molecule candidates as a drug have several limitations, but CADD utilizes novel methods that require little time and accurately predict a compound against a specific disease with minimal cost. Therefore, this review aims to provide a brief insight into the mathematical modeling and computational approaches for identifying a novel target and small molecular candidates for curing a specific disease. The comprehensive review mainly focuses on biological target prediction, structure-based and ligand-based drug design methods, molecular docking, virtual screening, pharmacophore modeling, quantitative structure-activity relationship (QSAR) models, molecular dynamics simulation, and MM-GBSA/MM-PBSA approaches along with valuable database resources and tools for identifying novel targets and therapeutics against a disease. This review will help researchers in a way that may open the road for the development of effective drugs and preventative measures against a disease in the future as early as possible.