Crystal contacts as nature's docking solutions

The assumption that crystal contacts reflect natural macromolecular interactions makes a basis for many studies in structural biology. However, the crystal state may correspond to a global minimum of free energy where biologically relevant interactions are sacrificed in favor to unspecific contacts. A large‐scale docking experiment was performed to assess the extent of misrepresentation of natural (in‐solvent) protein dimers by crystal packing. As found, the failure rate of docking may be quantitatively interpreted if both calculation errors and misrepresentation effects are taken into account. The failure rate analysis is based on the assumption that crystal structures reflect thermodynamic equilibrium between different dimeric configurations. The analysis gives an estimate of misrepresentation probability, which suggests that weakly bound complexes with K_D ≥ 100 μM (some 20% of all dimers in the PDB) have higher than 50% chances to be misrepresented by crystals. The developed theoretical framework is applicable in other studies, where experimental results may be viewed as snapshots of systems in thermodynamic equilibrium. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function,efficient optimization,and multithreading

AutoDock Vina, a new program for molecular docking and virtual screening, is presented. AutoDock Vina achieves an approximately two orders of magnitude speed‐up compared with the molecular docking software previously developed in our lab (AutoDock 4), while also significantly improving the accuracy of the binding mode predictions, judging by our tests on the training set used in AutoDock 4 development. Further speed‐up is achieved from parallelism, by using multithreading on multicore machines. AutoDock Vina automatically calculates the grid maps and clusters the results in a way transparent to the user. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Molecular docking studies of selected tricyclic and quinone derivatives on trypanothione reductase of Leishmania infantum

Visceral leishmaniasis, most lethal form of Leishmaniasis, is caused by Leishmania infantum in the Old world. Current therapeutics for the disease is associated with a risk of high toxicity and development of drug resistant strains. Thiol‐redox metabolism involving trypanothione and trypanothione reductase, key for survival of Leishmania, is a validated target for rational drug design. Recently published structure of trypanothione reductase (TryR) from L. infantum, in oxidized and reduced form along with Sb(III), provides vital clues on active site of the enzyme. In continuation with our attempts to identify potent inhibitors of TryR, we have modeled binding modes of selected tricyclic compounds and quinone derivatives, using AutoDock4. Here, we report a unique binding mode for quinone derivatives and 9‐aminoacridine derivatives, at the FAD binding domain. A conserved hydrogen bonding pattern was observed in all these compounds with residues Thr335, Lys60, His461. With the fact that these residues aid in the orientation of FAD towards the active site forming the core of the FAD binding domain, designing selective and potent compounds that could replace FAD in vivo during the synthesis of Trypanothione reductase can be deployed as an effective strategy in designing new drugs towards Leishmaniasis. We also report the binding of Phenothiazine and 9‐aminoacridine derivatives at the Z site of the protein. The biological significance and possible mode of inhibition by quinone derivatives, which binds to FAD binding domain, along with other compounds are discussed. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

A combined molecular modeling study on gelatinases and their potent inhibitors

Zinc‐dependent matrix metalloproteinase (MMP) family is considered to be an attractive target because of its important role in many physiological and pathological processes. In the present work, a molecular modeling study combining protein‐, ligand‐ and complex‐based computational methods was performed to analyze a new series of β‐N‐biaryl ether sulfonamide hydroxamates as potent inhibitors of gelatinase A (MMP‐2) and gelatinase B (MMP‐9). Firstly, the similarities and differences between the binding sites of MMP‐2 and MMP‐9 were analyzed through sequence alignment and structural superimposition. Secondly, in order to extract structural features influencing the activities of these inhibitors, quantitative structure‐activity relationship (QSAR) models using genetic algorithm‐multiple linear regression (GA‐MLR), comparative molecular field (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were developed. The proposed QSAR models could give good predictive ability for the studied inhibitors. Thirdly, docking study was employed to further explore the binding mode between the ligand and protein. The results from all the above analyses could provide the information about the similarities and differences of the binding mode between the MMP‐2, MMP‐9 and their potent inhibitors. The obtained results can provide very useful information for the design of new potential inhibitors. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

A new Lamarckian genetic algorithm for flexible ligand‐receptor docking

We present a Lamarckian genetic algorithm (LGA) variant for flexible ligand‐receptor docking which allows to handle a large number of degrees of freedom. Our hybrid method combines a multi‐deme LGA with a recently published gradient‐based method for local optimization of molecular complexes. We compared the performance of our new hybrid method to two non gradient‐based search heuristics on the Astex diverse set for flexible ligand‐receptor docking. Our results show that the novel approach is clearly superior to other LGAs employing a stochastic optimization method. The new algorithm features a shorter run time and gives substantially better results, especially with increasing complexity of the ligands. Thus, it may be used to dock ligands with many rotatable bonds with high efficiency. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Comparison of PBDE congener profiles and concentration levels in human specimens from China and the US and identification of human exposure sources

In an effort to investigate the status of human exposure to PBDEs in China,available monitoring data in human specimens(including breast milk,serums,and blood) was collected from the general population as well as specific groups that are occupationally exposed.PBDEs exposure profiles and concentration levels were compared with their counterparts in the United States of America.It was found that PBDE burdens in general Chinese population are one order lower and have different congener profiles from that in t...     

Indole- and Pyrazole-Glycyrrhetinic Acid Derivatives as PTP1B Inhibitors: Synthesis,In Vitro and In Silico Studies

Regulating insulin and leptin levels using a protein tyrosine phosphatase 1B (PTP1B) inhibitor is an attractive strategy to treat diabetes and obesity. Glycyrrhetinic acid (GA), a triterpenoid, may weakly inhibit this enzyme. Nonetheless, semisynthetic derivatives of GA have not been developed as PTP1B inhibitors to date. Herein we describe the synthesis and evaluation of two series of indole- and N-phenylpyrazole-GA derivatives (4a–f and 5a–f). We measured their inhibitory activity and enzyme kinetics against PTP1B using p-nitrophenylphosphate (pNPP) assay. GA derivatives bearing substituted indoles or N-phenylpyrazoles fused to their A-ring showed a 50% inhibitory concentration for PTP1B in a range from 2.5 to 10.1 µM. The trifluoromethyl derivative of indole-GA (4f) exhibited non-competitive inhibition of PTP1B as well as higher potency (IC₅₀ = 2.5 µM) than that of positive controls ursolic acid (IC₅₀ = 5.6 µM), claramine (IC₅₀ = 13.7 µM) and suramin (IC₅₀ = 4.1 µM). Finally, docking and molecular dynamics simulations provided the theoretical basis for the favorable activity of the designed compounds.     

Chromenol Derivatives as Novel Antifungal Agents: Synthesis,In Silico and In Vitro Evaluation

Herein we report the synthesis of some new 1H-1,2,4-triazole functionalized chromenols (3a–3n) via tandem reactions of 1-(alkyl/aryl)-2-(1H-1,2,4-triazole-1-yl) with salicylic aldehydes and the evaluation of their antifungal activity. In silico prediction of biological activity with computer program PASS indicate that the compounds have a high novelty compared to the known antifungal agents. We did not find any close analog among the over 580,000 pharmaceutical agents in the Cortellis Drug Discovery Intelligence database at the similarity cutoff of 70%. The evaluation of antifungal activity in vitro revealed that the highest activity was exhibited by compound 3k, followed by 3n. Their MIC values for different fungi were 22.1–184.2 and 71.3–199.8 µM, respectively. Twelve from fourteen tested compounds were more active than the reference drugs ketoconazole and bifonazole. The most sensitive fungus appeared to be Trichoderma viride, while Aspergillus fumigatus was the most resistant one. It was found that the presence of the 2-(tert-butyl)-2H-chromen-2-ol substituent on the 4th position of the triazole ring is very beneficial for antifungal activity. Molecular docking studies on C. albicans sterol 14α-demethylase (CYP51) and DNA topoisomerase IV were used to predict the mechanism of antifungal activities. According to the docking results, the inhibition of CYP51 is a putative mechanism of antifungal activity of the novel chromenol derivatives. We also showed that most active compounds have a low cytotoxicity, which allows us to consider them promising antifungal agents for the subsequent testing activity in in vivo assays.     

Application of the two-dimensional differential transform method to heat conduction problem for heat transfer in longitudinal rectangular and convex parabolic fins

In this article, approximate analytical (series) solutions for the temperature distribution in a longitudinal rectangular and convex parabolic fins with temperature dependent thermal conductivity and heat transfer coefficient are derived. The transient heat conduction problem is solved for the first time using the two-dimensional differential transform method (2D DTM). The effects of some physical parameters such as the thermo-geometric parameter, exponent and thermal conductivity gradient on temperature distribution are studied. Furthermore, we study the temperature profile at the fin tip.     

FT-midIR determination of fatty acid profiles, including trans fatty acids, in bakery products after focused microwave-assisted Soxhlet extraction

A study of the feasibility of Fourier transform medium infrared spectroscopy (FT-midIR) for analytical determination of fatty acid profiles, including trans fatty acids, is presented. The training and validation sets—75% (102 samples) and 25% (36 samples) of the samples once the spectral outliers have been removed—to develop FT-midIR general equations, were built with samples from 140 commercial and home-made bakery products. The concentration of the analytes in the samples used for this study is within the typical range found in these kinds of products. Both sets were independent; thus, the validation set was only used for testing the equations. The criterion used for the selection of the validation set was samples with the highest number of neighbours and the most separation between them (H<0.6). Partial least squares regression and cross validation were used for multivariate calibration. The FT-midIR method does not require post-extraction manipulation and gives information about the fatty acid profile in two min. The 14:0, 16:0, 18:0, 18:1 and 18:2 fatty acids can be determined with excellent precision and other fatty acids with good precision according to the Shenk criteria, R ²≥0.90, SEP=1–1.5 SEL and R ²=0.70-0.89, SEP=2–3 SEL, respectively. The results obtained with the proposed method were compared with those provided by the conventional method based on GC-MS. At 95% significance level, the differences between the values obtained for the different fatty acids were within the experimental error. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.