Exploration of histone deacetylase 8 inhibitors through classification QSAR study: Part II

Exploring molecular fragments for fraction unbound in human plasma of chemicals: a fragment-based cheminformatics approach

Fraction unbound in plasma (fu,p) of drugs is an significant factor for drug delivery and other biological incidences related to the pharmacokinetic behaviours of drugs. Exploration of different molecular fragments for fu,p of different small molecules/agents can facilitate in identification of suitable candidates in the preliminary stage of drug discovery. Different researchers have implemented strategies to build several prediction models for fu,p of different drugs. However, these studies did not focus on the identification of responsible molecular fragments to determine the fraction unbound in plasma. In the current work, we tried to focus on the development of robust classification-based QSAR models and evaluated these models with multiple statistical metrics to identify essential molecular fragments/structural attributes for fractions unbound in plasma. The study unequivocally suggests various N-containing aromatic rings and aliphatic groups have positive influences and sulphur-containing thiadiazole rings have negative influences for the fu,p values. The molecular fragments may help for the assessment of the fu,p values of different small molecules/drugs in a speedy way in comparison to experiment-based in vivo and in vitro studies.

Investigation of bacterial DNA gyrase Inhibitor classification models and structural requirements utilizing multiple machine learning methods

Infections from multidrug-resistant (MDR) bacteria have emerged as a paramount global health concern, and the therapeutic effectiveness of current treatments is swiftly diminishing. An urgent need exists to explore innovative strategies for countering drug-resistant bacteria. Bacterial DNA gyrase, functioning as an ATP-dependent enzyme, plays a pivotal role in the intricate processes of transcription, replication, and chromosome segregation within bacterial DNA. This renders it a prime target for the development of innovative antibacterial agents. However, the experimental identification of bacterial DNA gyrase inhibitors faces multifaceted challenges due to current methodological constraints. Recognizing its significance, this study developed 56 computational models designed for predicting bacterial DNA gyrase inhibitors. These models employed seven distinct molecular fingerprints and eight machine learning algorithms. Among these models, Model_2D, created using KlekotaRoth fingerprints and the SVM algorithm, stands out as the most robust performer (ACC = 0.86, MCC = 0.63, G-mean = 0.82). Moreover, given the limited exploration of structural fragments required for DNA Gyrase B inhibitors, crucial structural fingerprints influencing DNA Gyrase B inhibitors were identified through Bayesian classification. Subsequently, we conducted molecular docking to reveal the binding modes between these crucial structural fingerprints and the active site of DNA gyrase B. In conclusion, the present study aimed to develop the optimal classification model for bacterial DNA gyrase inhibitors, offering invaluable support to medicinal chemists creating innovative DNA gyrase inhibitors.

Unveiling critical structural features for effective HDAC8 inhibition: a comprehensive study using quantitative read-across structure-activity relationship (q-RASAR) and pharmacophore modeling

Histone deacetylases constitute a group of enzymes that participate in several biological processes. Notably, inhibiting HDAC8 has become a therapeutic strategy for various diseases. The current inhibitors for HDAC8 lack selectivity and target multiple HDACs. Consequently, there is a growing recognition of the need for selective HDAC8 inhibitors to enhance the effectiveness of therapeutic interventions. In our current study, we have utilized a multi-faceted approach, including Quantitative Structure-Activity Relationship (QSAR) combined with Quantitative Read-Across Structure-Activity Relationship (q-RASAR) modeling, pharmacophore mapping, molecular docking, and molecular dynamics (MD) simulations. The developed q-RASAR model has a high statistical significance and predictive ability (Q2F1:0.778, Q2F2:0.775). The contributions of important descriptors are discussed in detail to gain insight into the crucial structural features in HDAC8 inhibition. The best pharmacophore hypothesis exhibits a high regression coefficient (0.969) and a low root mean square deviation (0.944), highlighting the importance of correctly orienting hydrogen bond acceptor (HBA), ring aromatic (RA), and zinc-binding group (ZBG) features in designing potent HDAC8 inhibitors. To confirm the results of q-RASAR and pharmacophore mapping, molecular docking analysis of the five potent compounds (44, 54, 82, 102, and 118) was performed to gain further insights into these structural features crucial for interaction with the HDAC8 enzyme. Lastly, MD simulation studies of the most active compound (54, mapped correctly with the pharmacophore hypothesis) and the least active compound (34, mapped poorly with the pharmacophore hypothesis) were carried out to validate the observations of the studies above. This study not only refines our understanding of essential structural features for HDAC8 inhibition but also provides a robust framework for the rational design of novel selective HDAC8 inhibitors which may offer insights to medicinal chemists and researchers engaged in the development of HDAC8-targeted therapeutics.

The Importance of Structural Water in HDAC8 for Correct Binding Pose Applied for Drug Design of Anticancer Molecules

AIMS

Validating the docking procedure and maintaining the structural water molecules at HDAC8 catalytic site.

BACKGROUND

Molecular docking simulations play a significant role in Computer-Aided Drug Design, contributing to the development of new molecules. To ensure the reliability of these simulations, a validation process called "self-docking or re-docking" is employed, focusing on the binding mode of a ligand co-crystallized with the protein of interest.

OBJECTIVE

In this study, several molecular docking studies were conducted using five X-ray structures of HDAC8-ligand complexes from the PDB.

METHODS

Ligands initially complexed with HDAC8 were removed and re-docked onto the free protein, revealing a poor reproduction of the expected binding mode. In response to this, we observed that most HDAC8-ligand complexes contained one to two water molecules in the catalytic site, which were crucial for maintaining the cocrystallized ligand.

RESULTS

These water molecules enhance the binding mode of the co-crystallized ligand by stabilizing the proteinligand complex through hydrogen bond interactions between ligand and water molecules. Notably, these interactions are lost if water molecules are removed, as is often done in classical docking methodologies. Considering this, molecular docking simulations were repeated, both with and without one or two conserved water molecules near Zn+2 in the catalytic cavity. Simulations indicated that replicating the native binding pose of co-crystallized ligands on free HDAC8 without these water molecules was challenging, showing greater coordinate displacements (RMSD) compared to those including conserved water molecules from crystals.

CONCLUSION

The study highlighted the importance of conserved water molecules within the active site, as their presence significantly influenced the successful reproduction of the ligands' native binding modes. The results suggest an optimal molecular docking procedure for validating methods suitable for filtering new HDAC8 inhibitors for future experimental assays.

Exploring structural requirements of HDAC10 inhibitors through comparative machine learning approaches

Histone deacetylase (HDAC) inhibitors are in the limelight of anticancer drug development and research. HDAC10 is one of the class-IIb HDACs, responsible for cancer progression. The search for potent and effective HDAC10 selective inhibitors is going on. However, the absence of human HDAC10 crystal/NMR structure hampers the structure-based drug design of HDAC10 inhibitors. Different ligand-based modeling techniques are the only hope to speed up the inhibitor design. In this study, we applied different ligand-based modeling techniques on a diverse set of HDAC10 inhibitors (n = 484). Machine learning (ML) models were developed that could be used to screen unknown compounds as HDAC10 inhibitors from a large chemical database. Moreover, Bayesian classification and Recursive partitioning models were used to identify the structural fingerprints regulating the HDAC10 inhibitory activity. Additionally, a molecular docking study was performed to understand the binding pattern of the identified structural fingerprints towards the active site of HDAC10. Overall, the modeling insight might offer helpful information for medicinal chemists to design and develop efficient HDAC10 inhibitors.

A comparative quantitative structural assessment of benzothiazine-derived HDAC8 inhibitors by predictive ligand-based drug designing approaches

Histone deacetylase 8 (HDAC8) is a verified biomolecular target associated with diverse diseases including cancer. Though several HDAC inhibitors emerged effective against such diseases, no selective HDAC8 inhibitor is approved to date. Therefore, the development of potent HDAC8-selective inhibitors is inevitable to combat such diseases. Here, some benzothiazine-derived HDAC8 inhibitors were considered for a comparative QSAR analysis which may elucidate the prime structural components responsible for modulating their efficacy. Several outcomes from these diverse modelling techniques justified one another and thus validated each other. The ligand-based pharmacophore modelling study identified ring aromatic, positive ionizable, and hydrophobic features as essential structural attributes for HDAC8 inhibition. Besides, MLR, HQSAR and field-based 3D-QSAR studies signified the utility of the positive ionizable and hydrophobic features for potent HDAC8 inhibition. Again, the field-based 3D-QSAR study provided useful insight regarding the substitution in the fused phenyl ring. Moreover, the current observations also validated the previously reported molecular docking observations. Based on the outcomes, some new molecules were designed and predicted. Therefore, this comparative structural analysis of these HDAC8 inhibitors will surely assist in the development of potent HDAC8 inhibitors as promising anticancer therapeutics in the future.

A critical analysis of urea transporter B inhibitors: molecular fingerprints, pharmacophore features for the development of next-generation diuretics

Urea transporter is a membrane transport protein. It is involved in the transferring of urea across the cell membrane in humans. Along with urea transporter A, urea transporter B (UT-B) is also responsible for the management of urea concentration and blood pressure of human. The inhibitors of urea transporters have already generated a huge attention to be developed as alternate safe class of diuretic. Unlike conventional diuretics, these inhibitors are suitable for long-term therapy without hampering the precious electrolyte imbalance in the human body. In this study, UT-B inhibitors were analysed by using multi-chemometric modelling approaches. The possible pharmacophore features along with favourable and unfavourable sub-structural fingerprints for UT-B inhibition are extracted. This information will guide the medicinal chemist to design potent UT-B inhibitors in future.

Applying comparative molecular modelling techniques on diverse hydroxamate-based HDAC2 inhibitors: an attempt to identify promising structural features for potent HDAC2 inhibition

Histone deacetylase 2 (HDAC2) has been implicated in a variety of cardiovascular and neurodegenerative disorders as well as in cancers. Thus, HDAC2 has become an exclusive target for anticancer drug development. Therefore, the development of newer HDAC2 inhibitors in disease conditions is a prime goal to restrain such a scenario. Although a handful of HDAC inhibitors was accepted for the treatment of HDAC-related disease conditions, the non-selective nature of these entities is one of the major setbacks in the treatment of specific HDAC isoform-related pathophysiology. In this framework, the analyses of pre-existing molecules are essential to identify the important structural features that can fulfil the requirements for the cap and linker moieties to obtain potent and effective HDAC2 inhibition. Thus, in this study, the implementation of a combined comparative 2D and 3D molecular modelling techniques was done on a group of 92 diverse hydroxamate derivatives having a wide range of HDAC2 inhibitory potency. Besides other crucial features, this study upheld the importance of groups like triazole and benzyl moieties along with the molecular fields that are crucial for regulating HDAC2 inhibition. The outcomes of this study may be employed for the designing of HDAC2 inhibitors in future.

Ranked binding energies of residues and data fusion to identify the active and selective pyrimidine-based Janus kinases 3 (JAK3) inhibitors

The idea of using ranked binding energies of residues and data fusion are presented here for the first time as a valuable tool to classify active and selective inhibitors. Selective inhibitors of JAK3 can inhibit inflammatory cytokine while preventing targeting other subtypes of JAK1 and JAK2. Herein, we report a novel way to identify both active JAK3 and selective JAK1/JAK3 and JAK2/JAK3 inhibitors using the effective activity and selectivity classifications. The most important residues (top 10) responsible for the inhibition mechanism are sorted from high to low energies, which are considered as variables in the classification process. In addition, the ranked energies of ligands' heteroatoms (top 5), ranked energies of hydrogen bonds (top 5) and important molecular descriptors (top 10) were used to construct different data fusion possibilities. It is shown that the proposed data fusion strategy can increase the accuracy of the activity classification to 100% and the selectivity classification to 96.4%. The proposed strategies represented in this paper can help medicinal or pharmaceutical chemist in evaluation of both active and selective inhibitors before synthesizing new pharmaceuticals.

First molecular modelling report on tri-substituted pyrazolines as phosphodiesterase 5 (PDE5) inhibitors through classical and machine learning based multi-QSAR analysis

Phosphodiesterase 5 (PDE5) falls under a broad category of metallohydrolase enzymes responsible for the catalysis of the phosphodiesterase bond, and thus it can terminate the action of cyclic guanosine monophosphate (cGMP). Overexpression of this enzyme leads to development of a number of pathological conditions. Thus, targeting the enzyme to develop inhibitors could be useful for the treatment of erectile dysfunction as well as pulmonary hypertension. In the current study, several molecular modelling techniques were utilized including Bayesian classification, single tree and forest tree recursive partitioning, and genetic function approximation to identify crucial structural fingerprints important for optimization of tri-substituted pyrazoline derivatives as PDE5 inhibitors. Later, various machine learning models were also developed that could be utilized to predict and screen PDE5 inhibitors in the future.

Chemical-informatics approach to COVID-19 drug discovery: Exploration of important fragments and data mining based prediction of some hits from natural origins as main protease (Mpro) inhibitors

As the world struggles against current global pandemic of novel coronavirus disease (COVID-19), it is challenging to trigger drug discovery efforts to search broad-spectrum antiviral agents. Thus, there is a need of strong and sustainable global collaborative works especially in terms of new and existing data analysis and sharing which will join the dots of knowledge gap. Our present chemical-informatics based data analysis approach is an attempt of application of previous activity data of SARS-CoV main protease (Mpro) inhibitors to accelerate the search of present SARS-CoV-2 Mpro inhibitors. The study design was composed of three major aspects: (1) classification QSAR based data mining of diverse SARS-CoV Mpro inhibitors, (2) identification of favourable and/or unfavourable molecular features/fingerprints/substructures regulating the Mpro inhibitory properties, (3) data mining based prediction to validate recently reported virtual hits from natural origin against SARS-CoV-2 Mpro enzyme. Our Structural and physico-chemical interpretation (SPCI) analysis suggested that heterocyclic nucleus like diazole, furan and pyridine have clear positive contribution while, thiophen, thiazole and pyrimidine may exhibit negative contribution to the SARS-CoV Mpro inhibition. Several Monte Carlo optimization based QSAR models were developed and the best model was used for screening of some natural product hits from recent publications. The resulted active molecules were analysed further from the aspects of fragment analysis. This approach set a stage for fragment exploration and QSAR based screening of active molecules against putative SARS-CoV-2 Mpro enzyme. We believe the future in vitro and in vivo studies would provide more perspectives for anti-SARS-CoV-2 agents.

Exploring indole derivatives as myeloid cell leukaemia-1 (Mcl-1) inhibitors with multi-QSAR approach: a novel hope in anti-cancer drug discovery

In humans, the over-expression of Mcl-1 protein causes different cancers and it is also responsible for cancer resistance to different cytotoxic agents.