In silico screening for non-nucleoside HIV-1 reverse transcriptase inhibitors using physicochemical filters and high-throughput docking followed by in vitro evaluation

Strategies in the Design and Development of Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs)

AIDS (acquired immunodeficiency syndrome) is a potentially life-threatening infectious disease caused by human immunodeficiency virus (HIV). To date, thousands of people have lost their lives annually due to HIV infection, and it continues to be a big public health issue globally. Since the discovery of the first drug, Zidovudine (AZT), a nucleoside reverse transcriptase inhibitor (NRTI), to date, 30 drugs have been approved by the FDA, primarily targeting reverse transcriptase, integrase, and/or protease enzymes. The majority of these drugs target the catalytic and allosteric sites of the HIV enzyme reverse transcriptase. Compared to the NRTI family of drugs, the diverse chemical class of non-nucleoside reverse transcriptase inhibitors (NNRTIs) has special anti-HIV activity with high specificity and low toxicity. However, current clinical usage of NRTI and NNRTI drugs has limited therapeutic value due to their adverse drug reactions and the emergence of multidrug-resistant (MDR) strains. To overcome drug resistance and efficacy issues, combination therapy is widely prescribed for HIV patients. Combination antiretroviral therapy (cART) includes more than one antiretroviral agent targeting two or more enzymes in the life cycle of the virus. Medicinal chemistry researchers apply different optimization strategies including structure- and fragment-based drug design, prodrug approach, scaffold hopping, molecular/fragment hybridization, bioisosterism, high-throughput screening, covalent-binding, targeting highly hydrophobic channel, targeting dual site, and multi-target-directed ligand to identify and develop novel NNRTIs with high antiviral activity against wild-type (WT) and mutant strains. The formulation experts design various delivery systems with single or combination therapies and long-acting regimens of NNRTIs to improve pharmacokinetic profiles and provide sustained therapeutic effects.

Relative assessment of different statistical instruments and measures for the prediction of promising outcomes using docking, virtual screening and ADMET analysis against HIV-RT

Reverse transcriptase is the most therapeutic target for the discovery of novel, potent, and non-toxic new anti-retroviral drugs. In the present work, various docking software such as Sybyl Surflex-Dock, OpenEye FRED, and Hermes GOLD were evaluated for their efficiency to reproduce known cognate inhibitors' conformations. Three metrics were used and compared to assess the performance of the applied scoring functions, i.e. enrichment factor, receiver operating characteristic (ROC) curves, and Bedroc analysis. Twelve different scoring functions of three softwares were used to assess their ability to rank the cognate ligand within the active site of its proteins. The extensive virtual screening task was performed on eight crystal structures, and the performance of docking and scoring was assessed by their ability to efficiently detect known active compounds enriched in the top-ranked of the list among a randomly selected dataset of the ten thousand compounds of the NCI database. The effectiveness of post-docking relaxation in Surflex was also evaluated. The top 20, 50, and 100 compounds were selected based on consensus scoring functions from all 48 proteins with different ligand complexes. Further, the shortlisted leads were subjected to ADMET via using Discovery Studio. The results further implicate the importance of various statistical tools that should be followed before large-scale virtual screening for the drug discovery process. In silico results demonstrating the experiment was successful. The study of the research covers the combinatorial in silico techniques such as benchmarking of the softwares and scoring functions, statistical tools applied for screening and different conformations of HIV-RT crystal structures for virtual screening with rigid and flexible molecular docking and molecular dynamics simulation approach. This study reveals a clear roadmap to identify novel scaffolds against HIV-RT for antiretroviral therapy, thus providing the remedial solutions of HIV related infections and other diseases caused by malfunctioning of the target protein.Communicated by Ramaswamy H. Sarma.

Docking Methodologies and Recent Advances

Docking, a molecular modelling method, has wide applications in identification and optimization in modern drug discovery. This chapter addresses the recent advances in the docking methodologies like fragment docking, covalent docking, inverse docking, post processing, hybrid techniques, homology modeling etc. and its protocol like searching and scoring functions. Advances in scoring functions for e.g. consensus scoring, quantum mechanics methods, clustering and entropy based methods, fingerprinting, etc. are used to overcome the limitations of the commonly used force-field, empirical and knowledge based scoring functions. It will cover crucial necessities and different algorithms of docking and scoring. Further different aspects like protein flexibility, ligand sampling and flexibility, and the performance of scoring function will be discussed.

Promising Inhibitory Effects of Anthraquinones, Naphthopyrone, and Naphthalene Glycosides, from Cassia obtusifolia on α-Glucosidase and Human Protein Tyrosine Phosphatases 1B

The present work aims to evaluate the anti-diabetic potentials of 16 anthraquinones, two naphthopyrone glycosides, and one naphthalene glycoside from Cassia obtusifolia via inhibition against the protein tyrosine phosphatases 1B (PTP1B) and α-glucosidase. Among them, anthraquinones emodin and alaternin exhibited the highest inhibitory activities on PTP1B and α-glucosidase, respectively. Moreover, we examined the effects of alaternin and emodin on stimulation of glucose uptake by insulin-resistant human HepG2 cells. The results showed that alaternin and emodin significantly increased the insulin-provoked glucose uptake. In addition, our kinetic study revealed that alaternin competitively inhibited PTP1B, and showed mixed-type inhibition against α-glucosidase. In order to confirm enzyme inhibition, we predicted the 3D structure of PTP1B using Autodock 4.2 to simulate the binding of alaternin. The docking simulation results demonstrated that four residues of PTP1B (Gly183, Arg221, Ile219, Gly220) interact with three hydroxyl groups of alaternin and that the binding energy was negative (−6.30 kcal/mol), indicating that the four hydrogen bonds stabilize the open form of the enzyme and potentiate tight binding of the active site of PTP1B, resulting in more effective PTP1B inhibition. The results of the present study clearly demonstrate that C. obtusifolia and its constituents have potential anti-diabetic activity and can be used as a functional food for the treatment of diabetes and associated complications.

Design, discovery, modelling, synthesis, and biological evaluation of novel and small, low toxicity s-triazine derivatives as HIV-1 non-nucleoside reverse transcriptase inhibitors

A set of top-ranked compounds from a multi-objective in silico screen was experimentally tested for toxicity and the ability to inhibit the activity of HIV-1 reverse transcriptase (RT) in cell-free assay and in cell-based assay using HIV-1 based virus-like particles. Detailed analysis of a commercial sample that indicated specific inhibition of HIV-1 reverse transcription revealed that a minor component that was structurally similar to that of the main compound was responsible for the strongest inhibition. As a result, novel s-triazine derivatives were proposed, modelled, discovered, and synthesised, and their antiviral activity and cellular toxicity were tested. Compounds 18a and 18b were found to be efficient HIV-1 RT inhibitors, with an IC50 of 5.6±1.1μM and 0.16±0.05μM in a cell-based assay using infectious HIV-1, respectively. Compound 18b also had no detectable toxicity for different human cell lines. Their binding mode and interactions with the RT suggest that there was strong and adaptable binding in a tight (NNRTI) hydrophobic pocket. In summary, this iterative study produced structural clues and led to a group of non-toxic, novel compounds to inhibit HIV-RT with up to nanomolar potency.

Kinetics and molecular docking studies of fucosterol and fucoxanthin, BACE1 inhibitors from brown algae Undaria pinnatifida and Ecklonia stolonifera

Since the action of β-site amyloid precursor protein cleaving enzyme 1 (BACE1) is strongly correlated with the onset of Alzheimer's disease (AD), the development of BACE1 inhibitors as therapeutic agents is being vigorously pursued. In our ongoing research aimed at identifying anti-AD remedies derived from maritime plants, we evaluated the BACE1 inhibitory activities of fucosterol and fucoxanthin from Ecklonia stolonifera and Undaria pinnatifida. In vitro anti-AD activities were performed via BACE1 inhibition assays, as well as enzyme kinetic and molecular docking predictions. Based on enzyme-based assays, fucosterol and fucoxanthin showed noncompetitive and mixed-type inhibition, respectively, against BACE1. In addition, docking simulation results demonstrated that the Lys224 residue of BACE1 interacted with one hydroxyl group of fucosterol, while two additional BACE1 residues (Gly11 and Ala127) interacted with two hydroxyl groups of fucoxanthin. Moreover, the binding energy of fucosterol and fucoxanthin was negative (-10.1 and -7.0 kcal/mol), indicating that hydrogen bonding may stabilize the open form of the enzyme and potentiate tight binding of the active site of BACE1, resulting in more effective BACE1 inhibition. The results suggest that fucosterol and fucoxanthin may be used beneficially in the treatment of AD and provide potential guidelines for the design of new BACE1 inhibitors.

Anti-HIV Drug Discovery and Development: Current Innovations and Future Trends

The early effectiveness of combinatorial antiretroviral therapy (cART) in the treatment of HIV infection has been compromised to some extent by rapid development of multidrug-resistant HIV strains, poor bioavailability, and cumulative toxicities, and so there is a need for alternative strategies of antiretroviral drug discovery and additional therapeutic agents with novel action modes or targets. From this perspective, we first review current strategies of antiretroviral drug discovery and optimization, with the aid of selected examples from the recent literature. We highlight the development of phosphate ester-based prodrugs as a means to improve the aqueous solubility of HIV inhibitors, and the introduction of the substrate envelope hypothesis as a new approach for overcoming HIV drug resistance. Finally, we discuss future directions for research, including opportunities for exploitation of novel antiretroviral targets, and the strategy of activation of latent HIV reservoirs as a means to eradicate the virus.

The (5Z)-5-Pentacosenoic and 5-Pentacosynoic Acids Inhibit the HIV-1 Reverse Transcriptase

The natural fatty acids (5Z)-5-pentacosenoic and (9Z)-9-pentacosenoic acids were synthesized for the first time in eight steps starting from either 4-bromo-1-butanol or 8-bromo-1-butanol and in 20-58% overall yields, while the novel fatty acids 5-pentacosynoic and 9-pentacosynoic acids were also synthesized in six steps and in 34-43% overall yields. The ∆(5) acids displayed the best IC50's (24-38 µM) against the HIV-1 reverse transcriptase (RT) enzyme, comparable to nervonic acid (IC50 = 12 µM). The ∆(9) acids were not as effective towards HIV-RT with the (9Z)-9-pentacosenoic acid displaying an IC50 = 54 µM and the 9-pentacosynoic acid not inhibiting the enzyme at all. Fatty acid chain length and position of the unsaturation was important for the observed inhibition. None of the synthesized fatty acids were toxic (IC50 > 500 µM) towards peripheral blood mononuclear cells. Molecular modeling studies indicated the structural determinants underlying the biological activity of the most potent compounds. These results provide new insights into the structural requirements that must be present in fatty acids so as to enhance their inhibitory potential towards HIV-RT.

Protein tyrosine phosphatase 1B inhibitory activity of alkaloids from Rhizoma Coptidis and their molecular docking studies

ETHNOPHARMACOLOGIC RELEVANCE

Rhizoma Coptidis (the rhizome of Coptis chinensis Franch) has commonly been used for treatment of diabetes mellitus in traditional Chinese medicine due to its blood sugar-lowering properties and therapeutic benefits which highly related to the alkaloids therein. However, a limited number of studies focused on the Coptis alkaloids other than berberine.

MATERIALS AND METHODS

In the present study, we investigated the anti-diabetic potential of Coptis alkaloids, including berberine (1), epiberberine (2), magnoflorine (3), and coptisine (4), by evaluating the ability of these compounds to inhibit protein tyrosine phosphatase 1B (PTP1B), and ONOO(-)-mediated protein tyrosine nitration. We scrutinized the potentials of Coptis alkaloids as PTP1B inhibitors via enzyme kinetics and molecular docking simulation.

RESULTS

The Coptis alkaloids 1-4 exhibited remarkable inhibitory activities against PTP1B with the IC50 values of 16.43, 24.19, 28.14, and 51.04 μM, respectively, when compared to the positive control ursolic acid. These alkaloids also suppressed ONOO(-)-mediated tyrosine nitration effectively in a dose dependent manner. In addition, our kinetic study using the Lineweaver-Burk and Dixon plots revealed that 1 and 2 showed a mixed-type inhibition against PTP1B, while 3 and 4 noncompetitively inhibited PTP1B. Moreover, molecular docking simulation of these compounds demonstrated negative binding energies (Autodock 4.0=-6.7 to -7.8 kcal/mol; Fred 2.0=-59.4 to -68.2 kcal/mol) and a high proximity to PTP1B residues, including Phe182 and Asp181 in the WPD loop, Cys215 in the active sites and Tyr46, Arg47, Asp48, Val49, Ser216, Ala217, Gly218, Ile219, Gly220, Arg221 and Gln262 in the pocket site, indicating a higher affinity and tighter binding capacity of these alkaloids for the active site of the enzyme.

CONCLUSION

Our results clearly indicate the promising anti-diabetic potential of Coptis alkaloids as inhibitors on PTP1B as well as suppressors of ONOO(-)-mediated protein tyrosine nitration, and thus hold promise as therapeutic agents for the treatment of diabetes and related disease.

Supercomputer investigation of the protein-ligand system low-energy minima

The accuracy ofthe protein-ligand binding energy calculations andligand positioning isstrongly influenced by the choice of the docking target function. This work demonstrates the evaluation of the five different target functions used in docking: functions based on MMFF94 force field and functions based on PM7 quantum-chemical method accounting orwithout accounting the implicit solvent model (PCM, COSMO or SGB). For these purposes the ligand positions corresponding to the minima of the target function and the experimentally known ligand positions in the protein active site (crystal ligand positions) were compared. Each function was examined on the same test-set of 16 protein-ligand complexes. The new parallelized docking program FLM based on Monte Carlo search algorithm was developed to perform the comprehensive low-energy minima search and to calculate the protein-ligand binding energy. This study demonstrates that the docking target function based on the MMFF94 force field can be used to detect the crystal or near crystal positions of the ligand by the finding the low-energy local minima spectrum of the target function. The importance of solvent accounting in the docking process for the accurate ligand positioning is also shown. The accuracy of the ligand positioning as well as the correlation between the calculated and experimentally determined protein-ligand binding energies are improved when the MMFF94 force field is substituted by the new PM7 method with implicit solvent accounting.

Kinetics and molecular docking studies of kaempferol and its prenylated derivatives as aldose reductase inhibitors

Aldose reductase inhibitors (ARIs) suppressing the hyperglycemia-induced polyol pathway have been provided as potential therapeutic candidates in the treatment and prevention of diabetic complications. Based upon structure-activity relationships of desmethylanhydroicaritin (1) and sophoflavescenol (2) as promising ARIs, 3,4'-dihydroxy flavonols with a prenyl or lavandulyl group at the C-8 position and a hydroxyl or methoxy group at the C-5 position are important for aldose reductase (AR) inhibition. In order to prove the above results, a combination of computational prediction and enzyme kinetics has begun to emerge as an effective screening technique for the potential. In the present study, we predicted the 3D structure of AR in rat and human using a docking algorithm to simulate binding between AR and prenylated flavonoids (1 and 2) and kaempferol (3) and scrutinized the reversible inhibition of AR by these ARIs. Docking simulation results of 1-3 demonstrated negative binding energies (Autodock 4.0=-9.11 to -7.64 kcal/mol; Fred 2.0=-79.54 to -51.84 kcal/mol) and an additional hydrogen bond through Phe122 and Trp219, in addition to the previously proposed interaction of AR and phenolics through Trp20, Tyr48, His110, and Trp111 residues, indicating that the presence of 8-prenyl and 5-methyl groups might potentiate tighter binding to the active site of the enzyme and more effective AR inhibitors. Moreover, types of AR inhibition were different depending on the presence or absence of the 8-prenyl group, in that 1 and 2 are mixed inhibitors with respective Ki values of 0.69 μM and 0.94 μM, while 3 showed noncompetitive inhibition with a Ki value of 4.65 μM. The present study suggests that an effective strategy for screening potential ARIs could be established by predicting 3D structural conformation of prenylated flavonoids and the orientation within the enzyme as well as by simultaneously determining the mode of enzyme inhibition.

In silico virtual screening approaches for anti-viral drug discovery

Despite the considerable advances in medical and pharmaceutical research during the past years, diseases caused by viruses have remained a major burden to public health. Virtual in silico screening has repeatedly proven to be useful to meet the special challenges of antiviral drug discovery. Large virtual compound libraries are filtered by different computational screening methods such as docking, ligand-based similarity searches or pharmacophore-based screening, reducing the number of candidate molecules to a smaller set of promising candidates that are then tested biologically. This rational approach makes the drug discovery process more goal-oriented and saves resources in terms of time and money. In this review we discuss how different virtual screening techniques can be applied to antiviral drug discovery, present recent success stories in this field and finally address the main differences between the methods.:

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.

Discovery and development of anti-HBV agents and their resistance

Hepatitis B virus (HBV) infection is a prime cause of liver diseases such as hepatitis, cirrhosis and hepatocellular carcinoma. The current drugs clinically available are nucleot(s)ide analogues that inhibit viral reverse transcriptase activity. Most drugs of this class are reported to have viral resistance with breakthrough. Recent advances in methods for in silico virtual screening of chemical libraries, together with a better understanding of the resistance mechanisms of existing drugs have expedited the discovery and development of novel anti-viral drugs. This review summarizes the current status of knowledge about and viral resistance of HBV drugs, approaches for the development of novel drugs as well as new viral and host targets for future drugs.