Molecular docking and 3D-QSAR studies of HIV-1 protease inhibitors

Environmentally benign synthesis of unsymmetrical ureas and their evaluation as potential HIV-1 protease inhibitors via a computational approach

The present work reports the cost-effective, high yielding and environmentally acceptable preparation of unsymmetrical ureas from thiocarbamate salts using sodium percarbonate as an oxidant. Efficacy of the unsymmetrical ureas as potential human immune deficiency virus (HIV-1) protease inhibitors has been evaluated via in silico approach. The results revealed interactions of the urea compounds at the active site of the enzyme with favorable binding affinities causing possible mutations hindering the functioning of the enzyme. Further computational assessment of IC50 using known references satisfactorily authenticated the inhibitory action of the selected compounds against HIV-1 protease. Added to the easy synthesis of the ureas following an environmentally benign protocol, this work may be a valuable addition to the ongoing search for drugs with better efficacy profiles and reduced toxicity against HIV.

DFT based QSAR study on quinolone-triazole derivatives as antibacterial agents

QSAR modeling was performed on 39 quinolone-triazole derivatives against gram-positive Staphylococcus aureus and gram-negative Pseudomonas aeruginosa bacteria. The molecular structures were optimized using the DFT/B3LYP method and 6-31 G basis set. Molecular descriptors were extracted using quantum mechanical calculations. The hierarchical cluster analysis was performed for a rational subset division. The initial dataset was divided into calibration and validation sets, and modeling was done by stepwise MLR method for each of the two bacteria. Internal and external validation methods confirmed the robustness and predictability of the obtained models. According to the obtained model for S. aureus (R² = 0.889, R²_ext = 0.938, Q²_LOO = 0.853), the four descriptors- partial atomic charges for the N1 atom in triazole and C7 of the quinolone nucleus, 4-carbonyl bond length, and ¹³C-NMR chemical shift of 3-carboxylic acid- were found to be the descriptors controlling the activity. According to the obtained model for P. aeruginosa (R² = 0.957, R²_ext = 0.923, Q²_LOO = 0.909), the O atom's partial charge in carbonyl, LUMO-HOMO energy gap, and logP were found to be the descriptors having the highest correlation with the antibacterial activity. Finally, some new compounds with higher activities were designed and proposed.

The recent application of 3D-QSAR and docking studies to novel HIV-protease inhibitor drug discovery

INTRODUCTION

Despite the availability of FDA approved inhibitors of HIV protease, numerous efforts are still ongoing to achieve 'near-perfect' drugs devoid of characteristic adverse side effects, toxicities, and mutational resistance. While experimental methods have been plagued with huge consumption of time and resources, there has been an incessant shift towards the use of computational simulations in HIV protease inhibitor drug discovery.

AREAS COVERED

Herein, the authors review the numerous applications of 3D-QSAR modeling methods over recent years relative to the design of new HIV protease inhibitors from a series of experimentally derived compounds. Also, the augmentative contributions of molecular docking are discussed.

EXPERT OPINION

Efforts to optimize 3D QSAR and molecular docking for HIV-1 drug discovery are ongoing, which could further incorporate inhibitor motions at the active site using molecular dynamics parameters. Also, highly predictive machine learning algorithms such as random forest, K-means, decision trees, linear regression, hierarchical clustering, and Bayesian classifiers could be employed.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

Target-Specific Prediction of Ligand Affinity with Structure-Based Interaction Fingerprints

Discovery and optimization of small molecule inhibitors as therapeutic drugs have immensely benefited from rational structure-based drug design. With recent advances in high-resolution structure determination, computational power, and machine learning methodology, it is becoming more tractable to elucidate the structural basis of drug potency. However, the applicability of machine learning models to drug design is limited by the interpretability of the resulting models in terms of feature importance. Here, we take advantage of the large number of available inhibitor-bound HIV-1 protease structures and associated potencies to evaluate inhibitor diversity and machine learning models to predict ligand affinity. First, using a hierarchical clustering approach, we grouped HIV-1 protease inhibitors and identified distinct core structures. Explicit features including protein-ligand interactions were extracted from high-resolution cocrystal structures as 3D-based fingerprints. We found that a gradient boosting machine learning model with this explicit feature attribution can predict binding affinity with high accuracy. Finally, Shapley values were derived to explain local feature importance. We found specific van der Waals (vdW) interactions of key protein residues are pivotal for the predicted potency. Protein-specific and interpretable prediction models can guide the optimization of many small molecule drugs for improved potency.

QSAR, docking studies of 1,3-thiazinan-3-yl isonicotinamide derivatives for antitubercular activity

The enzyme - enoyl acyl carrier protein reductase (enoyl ACP reductase) is a validated target for antitubercular activity. Inhibition of this enzyme interferes with mycolic acid synthesis which is crucial for Mycobacterium tuberculosis cell growth. In the present work 2D and 3D quantitative structure activity relationship (QSAR) studies were carried out on a series of thiazinan-Isoniazid pharmacophore to design newer analogues. For 2D QSAR, the best statistical model was generated using SA-MLR method (r²=0.958, q²=0.922) while 3D QSAR model was derived using the SA KNN method (q²=0.8498). These studies could guide the topological, electrostatic, steric, hydrophobic substitutions around the nucleus based on which the NCEs were designed. Furthermore, molecular docking was performed to gauze the binding affinity of the designed analogues for enoyl ACP reductase enzyme. Amongst all the designed analogues the binding energies of SKS 01 and SKS 05 were found to be -5.267kcal/mol and -5.237kcal/mol respectively which was comparable with the binding energy of the standard Isoniazid (-6.254kcal/mol).

Molecular modeling studies of dihydro-alkyloxy-benzyl-oxopyrimidines (DABOs) as non-nucleoside inhibitors of HIV-1 reverse transcriptase using 3D-QSAR, Topomer CoMFA and molecular docking simulations

The human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) is generally regarded as a target for the treatment of acquired immune deficiency syndrome (AIDS).

A combined 3D-QSAR and docking studies for the In-silico prediction of HIV-protease inhibitors

BACKGROUND

Tremendous research from last twenty years has been pursued to cure human life against HIV virus. A large number of HIV protease inhibitors are in clinical trials but still it is an interesting target for researchers due to the viral ability to get mutated. Mutated viral strains led the drug ineffective but still used to increase the life span of HIV patients.

RESULTS

In the present work, 3D-QSAR and docking studies were performed on a series of Danuravir derivatives, the most potent HIV- protease inhibitor known so far. Combined study of 3D-QSAR was applied for Danuravir derivatives using ligand-based and receptor-based protocols and generated models were compared. The results were in good agreement with the experimental results. Additionally, docking analysis of most active 32 and least active 46 compounds into wild type and mutated protein structures further verified our results. The 3D-QSAR and docking results revealed that compound 32 bind efficiently to the wild and mutated protein whereas, sufficient interactions were lost in compound 46.

CONCLUSION

The combination of two computational techniques would helped to make a clear decision that compound 32 with well inhibitory activity bind more efficiently within the binding pocket even in case of mutant virus whereas compound 46 lost its interactions on mutation and marked as least active compound of the series. This is all due to the presence or absence of substituents on core structure, evaluated by 3D-QSAR studies. This set of information could be used to design highly potent drug candidates for both wild and mutated form of viruses.

Analysis of Ah receptor binding affinities of polybrominated diphenyl ethers via in silico molecular docking and 3D-QSAR

Polybrominated diphenyl ethers (PBDEs) have become ubiquitous contaminations due to their use as flame retardants. The structural similarity of PBDE to some dioxin-like compounds suggested that they may share similar toxicological effects: they might activate the aryl hydrocarbon receptor (AhR) signal transduction pathway and thus might have adverse effects on wildlife and humans. In this study, in silico computational workflow combining molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) was performed to investigate the binding interactions between PBDEs and AhR and the structural features affecting the AhR binding affinity of PBDE. The molecular docking showed that hydrogen-bond and hydrophobic interactions were the major driving forces for the binding of ligands to AhR, and several key amino acid residues were also identified. The CoMSIA model was developed from the conformations obtained from molecular docking and exhibited satisfactory results as q (2) of 0.605 and r (2) of 0.996. Furthermore, the derived model had good robustness and statistical significance in both internal and external validations. The 3D contour maps generated from CoMSIA provided important structural features influence the binding affinity. The obtained results were beneficial to better understand the toxicological mechanism of PBDEs.

Recent advances in fragment-based QSAR and multi-dimensional QSAR methods

This paper provides an overview of recently developed two dimensional (2D) fragment-based QSAR methods as well as other multi-dimensional approaches. In particular, we present recent fragment-based QSAR methods such as fragment-similarity-based QSAR (FS-QSAR), fragment-based QSAR (FB-QSAR), Hologram QSAR (HQSAR), and top priority fragment QSAR in addition to 3D- and nD-QSAR methods such as comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA), Topomer CoMFA, self-organizing molecular field analysis (SOMFA), comparative molecular moment analysis (COMMA), autocorrelation of molecular surfaces properties (AMSP), weighted holistic invariant molecular (WHIM) descriptor-based QSAR (WHIM), grid-independent descriptors (GRIND)-based QSAR, 4D-QSAR, 5D-QSAR and 6D-QSAR methods.