Search

Background: JAK2 inhibitors have been proposed as a new therapeutic option for thalassemia therapy. The objective of this study was to discover the key structural features for improving 2-aminopyrimidine derivatives as potential JAK2 inhibitors. Materials & methods: Quantitative structure–activity relationship (QSAR) approaches (hologram QSAR and comparative molecular similarity indices analysis), molecular dynamics simulations, binding energy calculations and pharmacokinetic predictions were employed. Results: Reliable QSAR models, binding mode and binding interactions of JAK2 inhibitors were obtained and these obtained results were used as the key information for rational design of highly potent JAK2 inhibitors. Conclusion: The concept of new potential JAK2 inhibitors integrated from the obtained results was proved, producing two newly designed compounds, D01 and D02, with potential for use as JAK2 inhibitors.

Aim: In silico screening approaches were performed to discover novel InhA inhibitors. Methods: Candidate InhA inhibitors were obtained from the combination of virtual screening and pharmacokinetic prediction. In addition, molecular mechanics Poisson–Boltzmann surface area, molecular mechanics Generalized Born surface area and WaterSwap methods were performed to investigate the binding interactions and binding energy of candidate compounds. Results: Four candidate compounds with suitable physicochemical, pharmacokinetic and antibacterial properties are proposed. The crucial interactions of the candidate compounds were H-bond, pi–pi and sigma–pi interactions observed in the InhA binding site. The binding affinity of these compounds was improved by hydrophobic interactions with hydrophobic side chains in the InhA pocket. Conclusion: The four newly identified InhA inhibitors reported in this study could serve as promising hit compounds against Mycobacterium tuberculosis and may be considered for further experimental studies.

The enoyl-acyl carrier protein reductase InhA of Mycobacterium tuberculosis is an attractive, validated target for antituberculosis drug development. Moreover, direct inhibitors of InhA remain effective against InhA variants with mutations associated with isoniazid resistance, offering the potential for activity against MDR isolates. Here, structure-based virtual screening supported by biological assays was applied to identify novel InhA inhibitors as potential antituberculosis agents. High-speed Glide SP docking was initially performed against two conformations of InhA differing in the orientation of the active site Tyr158. The resulting hits were filtered for drug-likeness based on Lipinski’s rule and avoidance of PAINS-like properties and finally subjected to Glide XP docking to improve accuracy. Sixteen compounds were identified and selected for in vitro biological assays, of which two (compounds 1 and 7) showed MIC of 12.5 and 25 μg/mL against M. tuberculosis H37Rv, respectively. Inhibition assays against purified recombinant InhA determined IC50 values for these compounds of 0.38 and 0.22 μM, respectively. A crystal structure of the most potent compound, compound 7, bound to InhA revealed the inhibitor to occupy a hydrophobic pocket implicated in binding the aliphatic portions of InhA substrates but distant from the NADH cofactor, i.e., in a site distinct from those occupied by the great majority of known InhA inhibitors. This compound provides an attractive starting template for ligand optimization aimed at discovery of new and effective compounds against M. tuberculosis that act by targeting InhA.

The aims of this study were to investigate the potential benefits of antioxidant, anti-inflammatory, anti-hepatotoxic, and anti-tyrosinase activities of a methanolic extract of fresh tea leaves (FTE) (Camellia sinensis L.). The antioxidant capacity was investigated using three different methods at different temperatures. The anti-inflammatory activity was studied in vitro by the inhibition of 5-lipoxygenase assay. The anti-hepatotoxic effect was investigated in CCl4-induced liver injury in rats. The anti-tyrosinase activities of the FTE and its principal phenolic compounds were investigated in l-3,4-dihydroxyphenylalanine (l-DOPA) oxidation by a mushroom tyrosinase. A molecular docking study was conducted to determine how the FTE’s principal catechins interact with the tyrosinase. The FTE exhibited the best shelf life at low temperatures and demonstrated concentration-dependent antioxidant, anti-inflammatory, anti-hepatotoxic, and anti-tyrosinase effects compared to positive references. Treatment of rats with the FTE at 2000 mg/kg/day for 28 consecutive days reversed CCl4-induced oxidative damage in hepatic tissues by lowering the levels of alanine aminotransferase by 69% and malondialdehyde by 90%. Our findings suggest that the FTE has the capacity to scavenge free radicals and can protect against oxidative stress induced by CCl4 intoxication. The docking results were consistent with our in vitro data, indicating the anti-tyrosinase potency of the principal catechins.

Azanaphthoquinone annelated pyrrole derivatives have been developed and synthesized with a continuous attempt to develop novel DNA intercalating agents as anti-cancer compounds with lower organ toxicity. With the remarkable antiproliferative activity of synthesized azanaphthoquinone annelated pyrrole derivatives, a structurally novel scaffold of these compounds is appropriated for further development of novel anti-cancer agents. Therefore, in the present study, 3D QSAR study (CoMSIA) was applied on 28 azanaphthoquinone annelated pyrrole derivatives to evaluate the structural requirement of these compounds. The resulting CoMSIA model is satisfied with r(2) of 0.99 and q(2) of 0.65. The interpretation of CoMSIA contours reveals the significant importance of steric, electrostatic, hydrophobic and hydrogen acceptor descriptors on the activities of azanaphthoquinone annelated pyrrole derivatives. Remarkably, the structural requirement of six substituent positions on the azanaphthoquinone annelated pyrrole scaffold was elucidated here. This result is the useful concept for design of new and more active azanaphthoquinone annelated pyrrole derivatives. Moreover, MD simulations using AMBER program were performed to model the binding of azanaphthoquinone annelated pyrrole derivatives in the intercalation site of the DNA duplex. Based on MD simulations, the information in terms of ligand-DNA interaction, complex structure and binding free energy was provided in this work. Therefore, the integrated results are informative for further modification of azanaphthoquinone annelated pyrrole scaffold leading to gain novel azanaphthoquinone annelated pyrrole derivatives possessing better antiproliferative activity.

Reverse transcriptase (RT), an essential enzyme for HIV-1 (human immunodeficiency virus type-1) life cycle, is a key target in drug discovery efforts against HIV-1 infection. Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are very specific to HIV-1 RT and have relatively less toxicity than the nucleoside reverse transcriptase inhibitors (NRTIs). However, the rapid emergence of drug-resistant viral strains has limited the therapeutic efficacy of these inhibitors. In this review, recent advances in computer-aided drug design (CADD) for design of new compounds against HIV-1 RT based on ligand-based drug design (LBDD) using 2D-and 3D-QSAR approaches, structure-based drug design (SBDD) with the combination of molecular docking, virtual screening and de novo drug design, molecular simulations and particular interaction calculated from quantum chemical calculations are discussed. Their successful applications are also highlighted.

Snakebite envenomations cause severe local tissue necrosis and the venom metalloproteinases are thought to be the key toxins involved. In this study, the ethanolic extract from seed kernels of Thai mango (Mangifera indica L. cv. 'Fahlun') (Anacardiaceae) and its major phenolic principle (pentagalloylglucopyranose) exhibited potent and dose-dependent inhibitory effects on the caseinolytic and fibrinogenolytic activities of Malayan pit viper and Thai cobra venoms in in vitro tests. molecular docking studies revealed that the binding orientations of the phenolic principles were in the binding pockets of snake venom metalloproteinases (SVMPs). The phenolic principles could form hydrogen bonds with the three histidine residues in the conserved zinc-binding motif and could chelate the Zn(2+) atom of the SVMPs, which could potentially result in inhibition of the venom enzymatic activities and thereby inhibit tissue necrosis.

The ethanolic extract from seed kernels of Thai mango (MSKE) ( Mangifera indica L. cv. ‘Fahlun’) (Anacardiaceae) and its major phenolic principle (pentagalloyl glucopyranose) exhibited dose-dependent inhibitory effects on enzymatic activities of phospholipase A 2 (PLA 2 ), hyaluronidase and L -amino acid oxidase (LAAO) of Calloselasma rhodostoma (CR) and Naja naja kaouthia (NK) venoms by in vitro tests. The anti-hemorrhagic and anti-dermonecrotic activities of MSKE against both venoms were clearly supported by in vivo tests. Molecular docking studies indicated that the phenolic molecules of the MSKE could selectively bind to the active sites or their proximity, or modify conserved residues that are critical for the catalysis of PLA 2 , and selectively bind to the LAAO binding pocket of both CR and NK venoms and thereby inhibit their enzymatic activities. The results imply a potential use of MSKE against snake venoms. Keywords: Calloselasma rhodostoma ; Mangifera indica L.; Molecular docking study

Comparative molecular field analysis (CoMFA) was performed on twenty-three pyrimethamine (pyr) derivatives active against quadruple mutant type (Asn51Ile, Cys59Arg, Ser108Asn, Ile164Leu) dihydrofolate reductase of Plasmodium falcipaarum (PfDHFR). The represented CoMFA models were evaluated based on the various three different probe atoms, C(sp3) (+1), O(sp3) (-1) and H (+1), resulting in the best model with combined three types of probe atoms. The statistical results were r(2)(cv) = 0.702, S(press) = 0.608, r(2)(nv) = 0.980, s = 0.156, and r(2)(test-set) = 0.698 which can explain steric contribution of about 50%. In addition, an understanding of particular interaction energy between inhibitor and surrounding residues in the binding pocket was performed by using MP2/6-31G(d,p) quantum chemical calculations. The obtained results clearly demonstrate that Asn108 is the cause of pyr resistance with the highest repulsive interaction energy. Therefore, CoMFA and particular interaction energy analyses can be useful for identifying the structural features of potent pyr derivatives active against quadruple mutant type PfDHFR.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous contaminants of the terrestrial environment that have been designated as Environmental Protection Agency (EPA) Priority Pollutants. In this study, molecular modeling was used to examine the physical and chemical characteristics of soil organic matter (SOM), fulvic acid (FA) and humic acid (HA), as well as their binding interactions with PAHs. The molecular structures of 18 PAHs were built by using the SYBYL 7.0 program and then fully optimized by a semiempirical (AM1) method. A molecular docking program, AutoDock 3.05, was used to calculate the binding interactions between the PAHs, and three molecular structure models including FA (Buffle’s model), HA (Stevenson’s model) and SOM (Schulten and Schnitzer’s model). The À-À interactions and H-bonding interactions were found to play an important role in the intermolecular bonding of the SOM/PAHs complexes. In addition, significant correlations between two chemical properties, boiling point (bp) and octanol/water partition coefficient (Log Kow) and final docking energies were observed. The preliminary docking results provided knowledge of the important binding modes to FA, HA and SOM, and thereby to predict the sorption behavior of PAHs and other pollutants.

A series of scopoletin derivatives incorporated with the pyridinium moiety was synthesized and evaluated for their acetylcholinesterase (AChE) inhibitory activity by the colorimetric Ellman's method. A 2-fluorobenzylpyridinium derivative was the most potent among the tested compounds, with an IC50 value of 0.215±0.015μM, which was greatly improved from that of scopoletin. Docking studies revealed that the scopoletin portion of the mentioned compound was bound to the peripheral anionic site of the AChE, whereas the N-benzylpyridinium residue to the catalytic anionic site.

The Mycobacterium tuberculosis protein kinase B (PknB) is critical for growth and survival of M. tuberculosis within the host. The series of aminopyrimidine derivatives show impressive activity against PknB (IC50 .5 μM). However, most of them show weak or no cellular activity against M. tuberculosis (MIC 63 μM). Consequently, the key structural features related to activity against of both PknB and M. tuberculosis need to be investigated. Here, two- and three-dimensional quantitative structure-activity relationship (2D and 3D QSAR) analyses combined with molecular dynamics (MD) simulations were employed with the aim to evaluate these key structural features of aminopyrimidine derivatives. Hologram quantitative structure-activity relationship (HQSAR) and CoMSIA models constructed from IC50 and MIC values of aminopyrimidine compounds could establish the structural requirements for better activity against of both PknB and M. tuberculosis. The NH linker and the R1 substituent of the template compound are not only crucial for the biological activity against PknB but also for the biological activity against M. tuberculosis. Moreover, the results obtained from MD simulations show that these moieties are the key fragments for binding of aminopyrimidine compounds in PknB. The combination of QSAR analysis and MD simulations helps us to provide a structural concept that could guide future design of PknB inhibitors with improved potency against both the purified enzyme and whole M. tuberculosis cells.

Ligand- and structure-based design approaches have been applied to an extended series of 74 efavirenz compounds effectively inhibiting wild type (WT) and mutant type (K103N) HIV-1 reverse transcriptase (RT). For ligand-based approach, three dimensional quantitative structure-activity relationship (3D-QSAR) methods, comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA), were performed. The starting geometry of efavirenz was obtained from X-ray crystallographic data. The efavirenz derivatives were constructed and fully optimized by ab-initio molecular orbital method at HF/3-21G level. Reliable QSAR models for high predictive abilities were developed. Regarding WT and K103N inhibitions, CoMFA models with r2/cv = 0.651 and 0.678 and CoMSIA models with r2/cv = 0.662 and 0.743 were derived, respectively. The interpretation obtained from the models highlights different structural requirements for inhibition of WT and K103N HIV-1 RT. To elucidate potential binding modes of efavirenz derivatives in the binding pocket of WT and K103N HIV-1 RT, structure-based approach based on computational docking studies of selected efavirenz compounds were performed by using GOLD and FlexX programs. The results derived from docking analysis give additional information and further probe the inhibitor-enzyme interactions. The correlation of the results obtained from 3D QSAR and docking models validate each other and lead to better understanding of the structural requirements for the activity. Therefore, these integrated results are informative to provide key features and a helpful guideline for novel compound design active against HIV-1 RT.

Nevirapine (Viramune) belongs to the first generation of non-nucleoside reverse transcriptase inhibitors (NNRTIs). Its efficiency is limited by drug resistant mutations, such as K103N and Y181C, so, the aim of this work was to design novel nevirapine analogues insensitive to the K103N and Y181C HIV-1 RT. 360 Nevirapine derivatives were designed using a combinatorial library design approach and these compounds were docked into the binding pocket of mutant HIV-1 RT enzyme structures, using the GOLD program. 124 Compounds having a GoldScore higher than that of nevirapine (55.00 and 52.00 for K103N and Y181C mutants, respectively) were first retrieved and submitted to a topological analysis with the SILVER program. Consequently, 31 compounds presenting a significant percentage of the surfaces buried upon binding (>80%) and exhibiting hydrogen bonds to either N103 or C181 residues of the HIV-RT were selected. To ensure that these compounds had hydrogen bonding interaction to either N103 or C181 residues, their interaction energies were estimated by quantum chemical calculations (QCCs). Finally, QCCs represent an alternative method for performing post docking procedure.

Substituted aminopyrimidine inhibitors have recently been introduced as antituberculosis agents. These inhibitors show impressive activity against protein kinase B, a Ser/Thr protein kinase that is essential for cell growth of M. tuberculosis. However, up to now, X-ray structures of the protein kinase B enzyme complexes with the substituted aminopyrimidine inhibitors are currently unavailable. Consequently, structural details of their binding modes are questionable, prohibiting the structural-based design of more potent protein kinase B inhibitors in the future. Here, molecular dynamics simulations, in conjunction with molecular mechanics/Poisson–Boltzmann surface area binding free-energy analysis, were employed to gain insight into the complex structures of the protein kinase B inhibitors and their binding energetics. The complex structures obtained by the molecular dynamics simulations show binding free energies in good agreement with experiment. The detailed analysis of molecular dynamics results shows that Glu93, Val95, and Leu17 are key residues responsible to the binding of the protein kinase B inhibitors. The aminopyrazole group and the pyrimidine core are the crucial moieties of substituted aminopyrimidine inhibitors for interaction with the key residues. Our results provide a structural concept that can be used as a guide for the future design of protein kinase B inhibitors with highly increased antagonistic activity.

A series of diphenyl ether derivatives were developed and showed promising potency for inhibiting InhA, an essential enoyl acyl carrier protein reductase involved in mycolic acid biosynthesis, leading to the lysis of Mycobacterium tuberculosis. To understand the structural basis of diphenyl ether derivatives for designing more potent inhibitors, molecular dynamics (MD) simulations were performed. Based on the obtained results, the dynamic behaviour in terms of flexibility, binding free energy, binding energy decomposition, conformation, and the inhibitor-enzyme interaction of diphenyl ether inhibitors were elucidated. Phe149, Tyr158, Met161, Met199, Val203 and NAD+ are the key residues for binding of diphenyl ether inhibitors in the InhA binding pocket. Our results could provide the structural concept to design new diphenyl ether inhibitors with better enzyme inhibitory activity against M. tuberculosis InhA. The present work facilitates the design of new and potentially more effective anti-tuberculosis agents.

Arylamides have been identified as direct InhA inhibitors which overcome the drug-resistance problem of isoniazid, the first-line drug for tuberculosis treatment. However, arylamide properties are not yet optimal against Mycobacterium tuberculosis . Arylamides show high potency in InhA enzyme assay, but they fail in antimycobacterial assay. To achieve the structural basis to improve antimycobacterial activity, the dynamic behavior of arylamide inhibitors and a substrate, trans -2-hexadecenoyl-( N -acetylcysteamine)-thioester, were carried out by molecular dynamics (MD) simulations. Arylamide inhibitors and a substrate are positioned at the same site which indicates the competitive inhibitor function of arylamides. Based on our findings, the amide carbonyl oxygen causes the selectivity of arylamide inhibitors for InhA inhibition. Moreover, this moiety is crucial for the affinity of the arylamide–InhA interactions with Tyr158 and NADH to form hydrogen bonds. It is possible to enhance the selectivity of arylamide inhibitors to reach the InhA target by introducing a hydrophilic substituent into the aryl ring A. In order to increase the membrane permeability of arylamide inhibitors, more lipophilic properties should be incorporated into the substituent B. Therefore, based on the obtained results, the correct balance between the selectivity and the membrane permeability of arylamide inhibitors should improve their inhibitory activity against M. tuberculosis strain.

The enoyl-ACP reductase enzyme (InhA) from M. tuberculosis is recognized as the primary target of isoniazid (INH), a first-line antibiotic for tuberculosis treatment. To identify the specific interactions of INH-NAD adduct and its derivative adducts in InhA binding pocket, molecular docking calculations and quantum chemical calculations were performed on a set of INH derivative adducts. Reliable binding modes of INH derivative adducts in the InhA pocket were established using the Autodock 3.05 program, which shows a good ability to reproduce the X-ray bound conformation with rmsd of less than 1.0 Å. The interaction energies of the INH-NAD adduct and its derivative adducts with individual amino acids in the InhA binding pocket were computed based on quantum chemical calculations at the MP2/6-31G (d) level. The molecular docking and quantum chemical calculation results reveal that hydrogen bond interactions are the main interactions for adduct binding. To clearly delineate the linear relationship between structure and activity of these adducts, CoMFA and CoMSIA models were set up based on molecular docking alignment. The resulting CoMFA and CoMSIA models are in conformity with the best statistical qualities, in which r2cv is 0.67 and 0.74, respectively. Structural requirements of isoniazid derivatives that can be incorporated into the isoniazid framework to improve the activity have been identified through CoMFA and CoMSIA steric and electrostatic contour maps. The integrated results from structure-based, ligand-based design approaches and quantum chemical calculations provide useful structural information facilitating the design of new and more potentially effective antitubercular agents as follow: the R substituents of isoniazid derivatives should contain a large plane and both sides of the plane should contain an electropositive group. Moreover, the steric and electrostatic fields of the 4-pyridyl ring are optimal for greater potency.

The alcoholic extract from seed kernels of Thai mango (Mangifera indica L. cv. ‘Fahlun’) (Anacardiaceae) and its major phenolic principle (pentagalloylglucopyranose) exhibited potent, dose-dependent inhibitory effects on tyrosinase with respect to L-DOPA. Molecular docking studies revealed that the binding orientations of the phenolic principles were in the tyrosinase binding pocket and their orientations were located in the hydrophobic binding pocket surrounding the binuclear copper active site. The results indicated a possible mechanism for their anti-tyrosinase activity which may involve an ability to chelate the copper atoms which are required for the catalytic activity of tyrosinase.