Docking mode of delvardine and its analogues into the p66 domain of HIV-1 reverse transcriptase: screening using molecular mechanics-generalized born/surface area and absorption, distribution, metabolism and excretion properties

Pharmacophore based virtual screening, molecular docking and molecular dynamic simulation studies for finding ROS1 kinase inhibitors as potential drug molecules

Proto-oncogene receptor tyrosine kinase ROS-1 is one of the clinically important biomarker and plays a crucial role in regulation of a number of cellular functions including cell proliferation, migration and angiogenesis. Recently, inhibition of ROS1 kinase has proven to be a promising target of anticancer drugs for non-small cell lung cancer (NSCLC). The very few compounds have been used as potent drug molecules so far and the selective ROS1 inhibitors are relatively rare. Besides the currently available drugs such as Crizotinib and PF-06463922 are becoming sensitive due to mutations in the ROS1 protein. To curtail the problem of the resistant, present study was designed to identify the potent inhibitors against ROS1. Three different screening approaches such as structure based, Atom-based and pharmacophore based screening were carried out against commercially available databases and the retrieved best hits were further evaluated by Lipinski's filter. Thereafter the lead molecule was subjected to pocket specific docking with ROS1. The results show that, total of 9 molecules (3 from each screening) has good docking score (with range of -9.288 to -12.49 Kcal/Mol) and binding interactions within the active site of ROS1. In order to analyze the stability of the ligand- protein complexes, molecular dynamics simulation was performed. Thus, these identified potential lead molecules with good binding score and binding affinity with ROS1 may act as the potent ROS1 inhibitor, and that are worth considering for further experimental studies.Communicated by Ramaswamy H. Sarma.

Pharmacoinformatics approaches to identify potential hits against tetraacyldisaccharide 4'-kinase (LpxK) of Pseudomonas aeruginosa

Pseudomonas aeruginosa infection can cause pneumonia and urinary tract infection and the management of Pseudomonas aeruginosa infection is critical in multidrug resistance, hospital-acquired bacteremia and ventilator-associated pneumonia. The key enzymes of lipid A biosynthesis in Pseudomonas aeruginosa are promising drug targets. However, the enzyme tetraacyldisaccharide 4'-kinase (LpxK) has not been explored as a drug target so far. Several pharmacoinformatics tools such as comparative metabolic pathway analysis (Metacyc), data mining from a database of essential genes (DEG), homology modeling, molecular docking, pharmacophore based virtual screening, ADMET prediction and molecular dynamics simulation were used in identifying novel lead compounds against this target. The top virtual hits STOCK6S-33288, 43621, 39892, 37164 and 35740 may serve as the templates for the design and synthesis of potent LpxK inhibitors in the management of serious Pseudomonas aeruginosa infection.

Identifying dual leucine zipper kinase (DLK) inhibitors using e-pharamacophore screening and molecular docking

Alzheimer's is a neural disorder causing gradual loss in structure and function of nerve cell. To treat such disorders, c-Jun N-terminal Kinase (JNK) Pathway inhibitors were developed by representing chemical compounds that were used to inhibit the JNK signaling pathways. DLK is the stress sensor and implicating as regulatory factor in JNK pathway. Therefore, in the present investigation, pharmacophore screening was tried to identify the chemical compounds that involving inhibition of DLK proteins. To explore the pharmacophore region and mode of binding with DLK protein, N- (I H-pyrazol-3-y l) pyridin-2-aminer inhibitors were docked with DLK. Results reveal the information on the interaction mechanism of protein and ligand with chemical characteristics required to inhibit DLK protein. Such predicted information (AAAARH) was used as query to find out potential novel lead compounds sourced from public database. As an outcome of 65 compounds were listed based on the fitness score (2≥), and were subjected to glide HTVS.SP and XP. Best performing 5 lead compounds were shortlisted for dynamic simulations. This exhibited a constant RMSD over 20 ns of timescale.

Identification of potential inhibitors for Penicillinbinding protein (PBP) from Staphylococcus aureus

Staphylococcus aureus is an infectious agent that causes severe skin and soft tissue infection in hospitalized patients. Therefore, it is of interest to develop potent inhibitors for S. aureus. Penicillin Binding protein (PBP) is a known drug target for inhibition of cell wall biosynthesis in S. aureus. Hence, PBP was screened with compounds from six databases using virtual screening approaches. Results shows that the screened lead compound produced higher docking score (-9.87 kcal/mol) compared to resistant drugs. Antimicrobial activity using screened lead compounds and resistant drugs showed maximum activity in potential screened compounds compared to resistant compounds.

In Silico study of the active site of Helicobacter pylori urease and its inhibition by hydroxamic acids

Hydroxamic acids have emerged as the most promising candidates from among the different classes of inhibitors of urease. In order to understand the mechanism of their action, we have studied in detail using quantum mechanics the active site of Helicobacter pylori urease complexed with acetohydroxamic acid. A diverse library of ligands having the hydroxamate moiety has been prepared and docked into the active site of urease using the QM/MM methodology. It is found that hydroxamic acids with hydrophobic groups attached to them are more potent inhibitors of urease because they can easily penetrate the hydrophobic environment surrounding the active site. The -CONHO^- moiety of the hydroxamic acid is also found to be absolutely necessary for chelation and inhibition of urease. In order to determine the roles of residues His 221 and Ala 365, which are not part of the active site, but are nevertheless involved in hydrogen bonding with the ligand, we have performed Molecular Dynamics simulations, both on the wild urease and also on its mutated counterpart, with the two residues substituted, respectively, by alanine and glycine.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.

Molecular Docking and Molecular Dynamics to Identify a Novel Human Immunodeficiency Virus Inhibitor from Alkaloids of Toddalia asiatica

Background:

Acquired immunodeficiency syndrome caused by human immunodeficiency virus (HIV) is an immunosuppressive disease. Over the past decades, it has plagued human health due to the grave consequences in its harness.

Objective:

For this reason, anti-HIV agents are imperative, and the search for the same from natural resources would assure the safety.

Materials and Methods:

In this investigation we have performed molecular docking, molecular property prediction, drug-likeness score, and molecular dynamics (MD) simulation to develop a novel anti-HIV drug. We have screened 12 alkaloids from a medicinal plant Toddalia asiatica for its probabilistic binding with the active site of the HIV-1-reverse transcriptase (HIV-1-RT) domain (the major contributor to the onset of the disease).

Results:

The docking results were evaluated based on free energies of binding (ΔG), and the results suggested toddanol, toddanone, and toddalenone to be potent inhibitors of HIV-1-RT. In addition, the alkaloids were subjected to molecular property prediction analysis. Toddanol and toddanone with more rotatable bonds were found to have a drug-likeness score of 0.23 and 0.11, respectively. These scores were comparable with the standard anti-HIV drug zidovudine with a model score 0.28. Finally, two characteristic protein-ligand complexes were exposed to MD simulation to determine the stability of the predicted conformations.

Conclusion:

The toddanol-RT complex showed higher stability and stronger H-bonds than toddanone-RT complex. Based on these observations, we firmly believe that the alkaloid toddanol could aid in efficient HIV-1 drug discovery.

SUMMARY

In the present study, the molecular docking and MD simulations are performed to explore the possible binding mode of HIV 1 RT with 12 alkaloids of T. asiatica. Molecular docking by AutoDock4 revealed three alkaloids toddanol, toddanone, and toddalenone with highest binding affinity towards HIV 1 RT. The drug likeness model score revealed a positive score for toddanol and toddanone which is comparable to the drug likeness score of the standard anti HIV drug zidovudine. Results from simulation analysis revealed that toddanol RT complex is more stable than toddanone RT complex inferring toddanol as a potential anti HIV drug molecule.

Abbreviations used: HIV: Human immunodeficiency virus, HIV 1 RT: HIV 1 reverse transcriptase, RNase H: Ribonuclease H, MD: Molecular dynamics, PDB: Protein databank, RMSD: Root mean square deviation, RMSF: Root mean square fluctuation.

Structural insights on identification of potential lead compounds targeting WbpP in Vibrio vulnificus through structure-based approaches

WbpP encoding UDP-GlcNAC C4 epimerase is responsible for the activation of virulence factor in marine pathogen Vibrio vulnificus (V. vulnificus) and it is linked to many aquatic diseases, thus making it a potential therapeutic target. There are few reported compounds that include several natural products and synthetic compounds targeting Vibrio sp, but specific inhibitor targeting WbpP are unavailable. Here, we performed structure-based virtual screening using chemical libraries such as Binding, TOSLab and Maybridge to identify small molecule inhibitors of WbpP with better drug-like properties. Deficient structural information forced to model the structure and the stable protein structure was obtained through 30 ns of MD simulations. Druggability regions are focused for new lead compounds and our screening protocol provides fast docking of entire small molecule library with screening criteria of ADME/Lipinski filter/Docking followed by re-docking of top hits using a method that incorporates both ligand and protein flexibility. Docking conformations of lead molecules interface displays strong H-bond interactions with the key residues Gly101, Ser102, Val195, Tyr165, Arg298, Val209, Ser142, Arg233 and Gln200. Subsequently, the top-ranking compounds were prioritized using the molecular dynamics simulation-based conformation and stability studies. Our study suggests that the proposed compounds may aid as a starting point for the rational design of novel therapeutic agents.

Aspartic acid functions as carbonyl trapper to inhibit the formation of advanced glycation end products by chemical chaperone activity

Advanced glycation end products (AGEs) were implicated in pathology of numerous diseases. In this study, we present the bioactivity of aspartic acid (Asp) to inhibit the AGEs. Hemoglobin and bovine serum albumin (BSA) were glycated with glucose, fructose, and ribose in the presence and absence of Asp (100-200 μM). HbA1c inhibition was investigated using human blood and characterized by micro-column ion exchange chromatography. The effect of methyl glyoxal (MG) on hemoglobin and BSA was evaluated by fluorescence spectroscopy and gel electrophoresis. The effect of MG on red blood cells morphology was characterized by scanning electron micrographs. Molecular docking was performed on BSA with Asp. Asp is capable of inhibiting the formation of fluorescent AGEs by reacting with the reducing sugars. The presence of Asp as supplement in whole blood reduced the HbA1c% from 8.8 to 6.1. The presence of MG showed an increase in fluorescence and the presence of Asp inhibited the glycation thereby the fluorescence was quenched. MG also affected the electrophoretic mobility of hemoglobin and BSA by forming high molecular weight aggregates. Normal RBCs showed typical biconcave shape. MG modified RBCs showed twisted and elongated shape whereas the presence of ASP tends to protect RBC from twisting. Asp interacted with arginine residues of bovine serum albumin particularly ARG 194, ARG 198, and ARG 217 thereby stabilized the protein complex. We conclude that Asp has dual functions as a chemical chaperone to stabilize protein and as a dicarbonyl trapper, and thereby it can prevent the complications caused by glycation.

Advances in computationally modeling human oral bioavailability

Although significant progress has been made in experimental high throughput screening (HTS) of ADME (absorption, distribution, metabolism, excretion) and pharmacokinetic properties, the ADME and Toxicity (ADME-Tox) in silico modeling is still indispensable in drug discovery as it can guide us to wisely select drug candidates prior to expensive ADME screenings and clinical trials. Compared to other ADME-Tox properties, human oral bioavailability (HOBA) is particularly important but extremely difficult to predict. In this paper, the advances in human oral bioavailability modeling will be reviewed. Moreover, our deep insight on how to construct more accurate and reliable HOBA QSAR and classification models will also discussed.

Mechanistic insights of SrtA–LPXTG blockers targeting the transpeptidase mechanism in Streptococcus mutans

The SrtA–LPXTG interaction plays a key role in transpeptidation reaction, cell wall and biofilm formations. This study explains the blocking of LEU interactions with SrtA will results as SrtA inhibitors through MD simulation and energy calculations methods.

Bioassay-guided isolation, identification and molecular ligand-target insight of lipoxygenase inhibitors from leaves of Anisomeles malabarica R.Br

Background:

Anisomeles malabarica R. Br. (Lamiaceae) is extensively used in traditional medicine in major parts of India for several medicinal purposes, including their use in rheumatism.

Materials and Methods:

The air-dried leaves of A. malabarica were extracted with ethanol, defatted with n-hexane and then successively partitioned into chloroform and n-butanol fractions. Bioassay-guided fractionation and purification of chloroform fraction from A. malabarica lead to the isolation of lipoxygenase (LOX) inhibitors. The structures of isolated compounds were elucidated by ultraviolet, infrared, ¹H nuclear magnetic resonance (NMR), ¹³C NMR and mass spectrometry spectroscopic techniques and assessed further by in vitro soybean lipoxygenase (sLOX) assay. In addition, the enzyme type inhibition was evaluated through molecular docking technique as a part of computational study.

Results:

The bioactive compounds 3, 4 dihydroxy benzoic acid (1) and 4’, 5, 7-trihydroxyflavone (2) were isolated from chloroform fraction of A. malabarica, whose bioactivity was observed to be dose-dependent compared to n-butanol fraction. Among the compounds, 3, 4 dihydroxy benzoic acid showed significant sLOX inhibitory activity with 74.04% ±2.6% followed by 4’, 5, 7-trihydroxyflavone (34.68% ±1.9%). The computational analysis of compounds showed their molecular interaction with important amino acid residues and nonheme iron atom in the catalytic site of LOX by enlightening their potential binding mode at molecular level.

Conclusions:

The LOX inhibitory constituents were identified from A. malabarica by means of bioassay-guided fractionation process. The results derived from in vitro and computational experiments confirm the potential of the isolated compounds and provide additional evidence for its traditional use in inflammatory disorders.

Exploring the selectivity of a ligand complex with CDK2/CDK1: a molecular dynamics simulation approach

Cyclin-dependent kinases (CDKs) are core components of the cell cycle machinery that govern the transition between phases during cell cycle progression. Abnormalities in CDKs activity and regulation are common features of cancer, making CDK family members attractive targets for the development of anticancer drugs. Their inhibitors have entered in clinical trials to treat cancer. Very recently, Heathcote et al. (J. Med. Chem. 2010, 53:8508-8522) have found a ligand BS194 that has a high affinity with CDK2 (IC(50) = 3 nM) but shows low affinity with CDK1 (IC(50) = 30 nM). To understand the selectivity, we used homology modeling, molecular docking, molecular dynamics, and free-energy calculation to analyze the interactions. A rational three-dimensional model of the CDK1/BS194 complex is built. We found that Leu83 is a key residue that recognizes BS194 more effectively with CDK2 with good binding free energies rather than CDK1. Energetic analysis reveals that van der Waals interaction and non-polar contributions to solvent are favorable in the formation of complexes and amine group of the ligand, which plays a crucial role for binding selectivity between CDK2 and CDK1.

Benzoxazinones as Human Peroxisome Proliferator Activated Receptor Gamma (PPARγ) Agonists: A Docking Study Using Glide

The purpose of the present study is to undertake a docking study of some benzoxazinone derivatives on human peroxisome proliferator activated receptor co-crystallized with an alpha-aryloxyphenylacetic acid agonist using Glide 4.5. The QikProp program was used to obtain the absorption, distribution, metabolism and excretion properties of the analogues. The intermolecular hydrogen bonding interaction of the best-fit ligands were found to be associated with Tyr473, Ser289, Hie 449, Hip 323, Ser 342 and Gly 284 amino acid residue at the receptor active site. Among all the observed interaction with similar binding pattern, the presence of methyl carboxypentyl side chain (Lig. No. 21) showed additional interaction with Ser 342 and the affinity was increased by carboxyl oxygen (as hydrogen bond acceptor) with a best Glide score of -14.54 as compared to the co-crystallized aryloxyphenyl acetic acid which achieved a glide score of -12.50.

Synthesis, characterization and antiproliferative activity of 1,2-naphthoquinone and its derivatives

In the present study substituted 1,2-naphthoquinones were synthesized, purified and characterized by spectroscopic studies (UV, FT-IR, ¹H NMR, ¹³ C NMR and elemental analysis). These compounds were evaluated for cytotoxicity against a panel of human cancer cell lines (Hep-G₂ for liver sarcoma, MG-63 for osteosarcoma and MCF-7 for human breast cancer). The cells were dosed with these ortho-naphthoquinone derivatives at varying concentrations, and cell viability was measured by a 3-(4,5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide (MTT) assay with doxorubicin as positive control. Significant anticancer activities were observed in vitro for some members of the series, and compounds 1,2-naphthoquinone 2-thiosemicarbazone, 1,2-naphthoquinone-2-semicarbazone, 4-amino-1,2-naphthoquinone 2-thiosemicarbazone and 4-amino-1,2-naphthoquinone-2-semicarbazone are active cytotoxic agents against different cancer cell lines with IC₅₀ values in the range of 5.73-17.67 μM. The obtained data suggested that better anticancer activity was linked with introduction of thiosemicarbazone and semicarbazone moiety in 1,2-naphthoquinone ring system. Outcomes of experimentation also reveal that incorporation of amino group in 1,2-naphthoquinone moiety contributes positively for cytotoxic action of compounds. Docking experiments showed a good correlation between their calculated interaction energies with the topoisomerase-II and the observed IC₅₀ values of all these compounds.

HIV-1 NNRTIs: structural diversity, pharmacophore similarity, and implications for drug design

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) nowadays represent very potent and most promising anti-AIDS agents that specifically target the HIV-1 reverse transcriptase (RT). However, the effectiveness of NNRTI drugs can be hampered by rapid emergence of drug-resistant viruses and severe side effects upon long-term use. Therefore, there is an urgent need to develop novel, highly potent NNRTIs with broad spectrum antiviral activity and improved pharmacokinetic properties, and more efficient strategies that facilitate and shorten the drug discovery process would be extremely beneficial. Fortunately, the structural diversity of NNRTIs provided a wide space for novel lead discovery, and the pharmacophore similarity of NNRTIs gave valuable hints for lead discovery and optimization. More importantly, with the continued efforts in the development of computational tools and increased crystallographic information on RT/NNRTI complexes, structure-based approaches using a combination of traditional medicinal chemistry, structural biology, and computational chemistry are being used increasingly in the design of NNRTIs. First, this review covers two decades of research and development for various NNRTI families based on their chemical scaffolds, and then describes the structural similarity of NNRTIs. We have attempted to assemble a comprehensive overview of the general approaches in NNRTI lead discovery and optimization reported in the literature during the last decade. The successful applications of medicinal chemistry strategies, crystallography, and computational tools for designing novel NNRTIs are highlighted. Future directions for research are also outlined.