一个组织型纤溶性溶解酶原激活剂(T-pa)非共价类抑制剂的药效模型

为了研究组织型纤溶性溶解酶原激活剂（T－pa)的非共价类抑制剂的三维定量构效关系，本文采用分子模拟软件Catalyst4.0(molecular simulation company)构建T-pa的非共价类抑制剂的三维药效模型。14个非共价类的抑制剂及其体外活性数据被用于构建此药效模型。此药效模型有3个疏水区特性，一个正离子化区特性和氢键供体区特性，且结构与活性相关系数为r=0.962。与T－pa的晶体结构比较，此药物模型在化学状态和疏水性上与其能很好的相匹配。且利用此是找到的抑制剂（bbzi14)的活性构象，与晶体结构中此抑制剂的结合构象基本一致。通过此药效集团模型反映的抑制剂与受体的相互作用模式，能为发现新型亲和配基和抑制剂提供有用的启示。     

Structure-based design of inhibitors of the rice blast fungal enzyme trihydroxynaphthalene reductase

Rice Blast Disease, caused by the fungus Pyricularia oryzae, is one of the most important diseases of rice. Several enzymes in the melanin biosynthetic pathway have proven to be valuable targets for development of rice blast fungicides. In particular, inhibitors of trihydroxynaphthalene reductase (3HNR), which catalyzes the conversion of trihydroxynaphthalene to vermelone, have yielded commercially useful rice fungicides. The X-ray structure of 3HNR has been published recently, presenting an opportunity to use this information in the de novo design of novel 3HNR inhibitors that may exhibit useful rice blast activity. We used the LeapFrog program to develop a docking model for interaction of ligands with the active site of THNR. The final model gave a good correlation between calculated binding energy and log Ki and was used to design novel ligands and score compounds for synthesis. Using this as a tool, we synthesized inhibitors in the nanomolar range and also developed several inhibitors that did not conform to the properties of the THNR active site. Leapfrog was able to locate a previously unrecognized binding pocket that could accommodate these otherwise anomalous regions of structure.     

一个组织型纤溶性溶解酶原激活剂(T-pa)非共价类抑制剂的药效模型

为了研究组织型纤溶性溶解酶原激活剂(T-pa)的非共价类抑制剂的三维定量构效关系,本文采用分子模拟软件Catalyst 4.0(molecular simulation company)构建T-pa的非共价类抑制剂的三维药效模型。14个非共价类的抑制剂及其体外活性数据被用于构建此药效模型。此药效模型有3个疏水区特性,一个正离子化区特性和氢键供体区特性,且结构与活性相关系数为r=0.962。与T-pa的晶体结构比较,此药物模型在化学状态和疏水性上与其能很好的相匹配。且利用此模型查找到的抑制剂(bbzi14)的活性构象,与晶体结构中此抑制剂的结合构象基本一致。通过此药效集团模型反映的抑制剂与受体的相互作用模式,能为发现新型亲和配基和抑制剂提供有用的启示。     

QSAR and molecular graphics analysis of N2-phenylguanines as inhibitors of herpes simplex virus thymidine kinases

A quantitative structure-activity relationship study of N2-(substituted)-phenylguanines (PHG) as inhibitors of herpes simplex virus thymidine kinase (HSV TK) was performed. The activity of a set of PHG derivatives were analyzed against the thymidine kinase of herpes simplex virus types 1 (HSV1 TK) and 2 (HSV2 TK). Classic and calculated physicochemical parameters were included in the analysis. The results showed that there is an important difference in the activity of the meta substituted PHG derivatives against HSV1 TK and HSV2 TK. The activity of the meta derivatives against HSV2 TK is influenced by a steric effect, which is not observed against HSV1 TK. The superposition of the three-dimensional structures of the active sites of HSV1 TK (crystal structure) and HSV2 TK (homology model) revealed that the amino acid Ile97 is located near the meta position in the HSV1 TK active site, whereas the amino acid Leu97 is located near the meta position in the HSV2 TK active site. This single difference in the active sites of both enzymes can explain the source of the steric effect and serves as an indication that our previously proposed binding mode for the PHG derivatives is plausible. However, another observed mutation in the active site region, Ala168 by Ser168, suggests that an alternative binding mode, similar to that of ganciclovir, could be possible.     

Exploring the P2 and P3 ligand binding features for hepatitis C virus NS3 protease using some 3D QSAR techniques

Several three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been constructed using the comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and catalyst pharmacophore feature building programs for a series of 26 truncated ketoacid inhibitors designed particularly for exploring the P2 and P3 binding pockets of HCV NS3 protease. The structures of these inhibitors were built from a structure template extracted from the crystal structure of HCV NS3 protease. The structures were aligned through docking each inhibitor into the NS3 active site using program GOLD. The best CoMSIA model was identified from the stepwise analysis results and the corresponding pharmacophore features derived were used for constructing a pharmacophore hypothesis by the catalyst program. Pharmacophore features obtained by CoMFA and CoMSIA are found to be in accord with each other and are both mapped onto the molecular 5K surface of NS3 active site. These pharmacophore features were also compared with those obtained by the catalyst program and mapped onto the same NS3 molecular surface. The pharmacophore building process was also performed for 20 boronic acid based NS3 inhibitors characterized by a long hydrophobic side chain attached at position P2. This latter pharmacophore hypothesis built by the catalyst program was also mapped onto the molecular surface of NS3 active site to define a second hydrophobic feature at position P2. The possibility of using the pharmacophore features mapped P2 and P3 binding pocket to design more potent depeptidized NS3 inhibitors was discussed.     

Evaluation of the combination mode of azoles antifungal inhibitors with CACYP51 and the influence of Site-directed mutation

14α-demethylase (CYP51) is an essential metabolic enzyme for fungal survival and has been considered as an interesting target for the development of new antifungal inhibitors. Azoles antifungal inhibitors in the treatment of fungal diseases are good candidates via the interaction with the target enzyme CYP51 of fungus. In the study, we constructed the homology model for Candida albicans CYP51 (CACYP51) and analyzed the active site. In order to better understand the structural characteristics of azoles inhibitors and combination mode, the common feature pharmacophore model and the molecular docking were performed. The results suggest that the azoles inhibitors consist of three chemical features: the aromatic groups, phenyl groups and the azoles groups. The aromatic groups of inhibitors occupy the upper of active pocket, the phenyl groups and azoles groups occupy the bottom of active pocket. Further validation studies found these amino acid residues Tyr118, His310 and Ser378 play an important role in the substrate binding, and these amino acid residues with site-directed mutation will weaken the combining ability of the inhibitors.     

Pharmacophore modeling,virtual screening,docking and in silico ADMET analysis of protein kinase B (PKB β) inhibitors

Protein kinase B (PKB) is a key mediator of proliferation and survival pathways that are critical for cancer growth. Therefore, inhibitors of PKB are useful agents for the treatment of cancer. Herein, we describe pharmacophore-based virtual screening combined with docking study as a rational strategy for identification of novel hits or leads. Pharmacophore models of PKB β inhibitors were established using the DISCOtech and refined with GASP from compounds with IC₅₀ values ranging from 2.2 to 246 nM. The best pharmacophore model consists of one hydrogen bond acceptor (HBA), one hydrogen bond donor (HBD) site and two hydrophobic (HY) features. The pharmacophore models were validated through receiver operating characteristic (ROC) and Güner-Henry (GH) scoring methods indicated that the model-3 was statistically valuable and reliable in identifying PKB β inhibitors. Pharmacophore model as a 3D search query was searched against NCI database. Several compounds with different structures (scaffolds) were retrieved as hits. Molecules with a Q_fit value of more than 95 and three other known inhibitors were docked in the active site of PKB to further explore the binding mode of these compounds. Finally in silico pharmacokinetic and toxicities were predicted for active hit molecules. The hits reported here showed good potential to be PKB β inhibitors.     

Docking of sulfonamides to carbonic anhydrase II and IV

Starting with a known active site of a protein and a database of compounds, one would like to quickly identify a few compounds that "dock" into the active site and obtain "good" binding free energies. The main goal of current automated docking procedures is to predict the "best" substrate-enzyme complex while other programs such as UHBD and DelPhi can be used to compute binding free energies. In this paper, we will focus on the application of docking methods and parameters to study substrate-enzyme interactions of a metalloenzyme system. Specifically, we report on the docking of sulfonamides to carbonic anhydrase II and IV, which are of interest due to their application in glaucoma therapy. Using a standard docking protocol, it is possible to correctly predict not only the orientation of inhibitors to a specific isozyme, but also determine the qualitative affinity for a group of inhibitor for an isozyme.     

Structure-based design of inhibitors of NS3 serine protease of hepatitis C virus

We have designed small focused combinatorial library of hexapeptide inhibitors of NS3 serine protease of the hepatitis C virus (HCV) by structure-based molecular design complemented by combinatorial optimisation of the individual residues. Rational residue substitutions were guided by the structure and properties of the binding pockets of the enzyme's active site. The inhibitors were derived from peptides known to inhibit the NS3 serine protease by using unusual amino acids and alpha-ketocysteine or difluoroaminobutyric acid, which are known to bind to the S1 pocket of the catalytic site. Inhibition constants (Ki) of the designed library of inhibitors were predicted from a QSAR model that correlated experimental Ki of known peptidic inhibitors of NS3 with the enthalpies of enzyme-inhibitor interaction computed via molecular mechanics and the solvent effect contribution to the binding affinity derived from the continuum model of solvation. The library of the optimised inhibitors contains promising drug candidates-water-soluble anionic hexapeptides with predicted Ki* in the picomolar range.     

氰基丙烯酸酯类化合物的3D-QSAR研究

3D-Quantitative Structure—Activity Relationships of a series of cyanoacrylate inhibitors which are known to act by blocking photosynthetic electron transport close to photosystem II reaction center in the thylakoid membranes of plant chloroplasts have been investigated using comparative molecular field analysis(CoMFA). The active conformations of the title compounds were obtained by constrained systematic search program. The resulting CoMFA model with considerable predictive ability shows that the contribution of steric effects for activities of cyanoacrylate inhibitors is more important than electrostatic effects.     

Potent inhibitors precise to S1′ loop of MMP-13, a crucial target for osteoarthritis

Matrix metalloproteinase-13 (MMP-13) is the primary MMP involved in cartilage degradation through its particular ability to cleave type-II collagen. This protein is expressed by chondrocytes and synovial cells in human osteoarthritis and rheumatoid arthritis; hence, it is an attractive target for the treatment of arthritic diseases. Currently available inhibitors lack specificity for metalloproteinase because of a common Zn binding site in MMPs; thus, there is a need to identify selective MMP-13 inhibitors for osteoarthritis therapy. Because selectivity is the major concern, both ligand-based and protein-based pharmacophore methodologies were used to identity potent and selective MMP-13 inhibitors. Different hypotheses were validated, and the best hypothesis was used to screen Zinc (natural and chemical) databases to seek novel scaffolds as MMP-13 inhibitors. The identified hits were validated using different strategies, such as Glide Standard precision, extra precision, E-model energies and receiver operating curve (ROC). In addition, potent inhibitors were selected based on two criteria: a similar binding mode as that of the active site PB3 crystal ligand and crucial amino acid interactions that are catalytically important for the function of MMP-13. The candidate potent inhibitors ZINC 02535232, ZINC 08399795, ZINC 12419118 and ZINC 00624580 nearly reproduced the H-bond interactions formed in the crystal structure of 1XUC with reasonable RMSD values exhibiting a novel interaction pattern that was not previously observed in MMP-13 inhibitors. The identified potent hits with diverse chemical scaffolds may be useful in designing new MMP-13 inhibitors.     

新型二芳基取代-1,2,4-三唑类化合物的抗炎构效关系和分子对接研究

目的:建立新型二芳基取代-1,2,4-三唑类化合物的定量构效关系,设计新型COX-2抑制剂.方法:采用PM3半经验量子化学法全优化18种二芳基取代-1,2,4-三唑类选择性环氧合酶(COX-2)抑制剂的结构,从数据中搜寻或计算它们的226种参数,利用逐步回归法,建立经典结构-活性关系(2D-QSAR);用Autodock对接软件研究二芳基取代-1,2,4-三唑类化合物与环氧合酶(COX-2)的对接,分析该类化合物与环氧合酶在复合物的立体结构以及分子对接自由能与抑制活性的关系.结果:建立合理二芳基取代-1,2,4-三唑类化合物COX-2抑制剂定量构效关系,表明活性二芳基取代-1,2,4-三唑类选择性环氧合酶(COX-2)抑制剂具有类似塞来昔布等三环类环氧合酶-2抑制剂的立体结构,并且对接自由能与抑制剂活性有较好的相关性.结论:所得的模型可以解释已有的构效关系,而且预测同类化合物能力较好,可指导设计新抑制剂.     

Modeling and molecular dynamics of glutamine transaminase K/cysteine conjugate beta-lyase

The homodimeric, pyridoxal 5'-phosphate (PLP)-dependent enzyme glutamine transaminase K/cysteine conjugate beta-lyase (GTK/beta-lyase) has been implicated in the bioactivation of chemopreventive compounds. This paper describes the first homology model of rat renal GTK/beta-lyase and its active site residues, deduced from molecular dynamics (MD) simulations of the binding mode of 13 structurally diverse cysteine S-conjugates and amino acids after Amber-parametrization of PLP. Comparison with Thermus thermophilus aspartate aminotransferase (tAAT) and Trypanosoma cruzi tyrosine aminotransferase (tTAT), used as templates for modeling GTK/beta-lyase, showed that the PLP-binding site of GTK/beta-lyase is highly conserved. Binding of the ligand alpha-carboxylate-group occurred via the conserved residues Arg(432) and Asn(219), and Asn(50) and Gly(70). Two pockets accommodated the various ligand side chains. A small pocket, located directly above PLP, was of a highly hydrophobic and aromatic character. A larger pocket, formed partly by the substrate access channel, was more hydrophilic and notably involved the salt bridge partners Glu(54) and Arg(99*) (* denotes the other subunit). Ligand-binding residues included Leu(51), Phe(71), Tyr(135), Phe(373) and Phe(312*), and pi-stacking interactions were often observed. Tyr(135) and Asn(50) were prominent in hydrogen bonding with the sulfur-atom of cysteine S-conjugates.The observed binding mode of the ligands corresponded well with their experimentally determined inhibitory potency toward GTK/beta-lyase. The current homology model thus provides a starting point for further validation of the role of active site residues in ligand-binding by means of mutagenesis studies. Ultimately, insight in the binding of ligands to GTK/beta-lyase may result in the rational design of new ligands and selective inhibitors.     

Assessment of a mechanism for reactive inhibition of carboxypeptidase A with QM/MM methods

The mechanism for inhibition of carboxypeptidase A (CPA) by the two enantiomers of a reactive inhibitor, N-(2-chloroethyl)-N-methylphenylalanine, has been investigated using computational methods. Quantum mechanical and molecular mechanical (QM/MM) methods have been employed to find likely enzyme binding conformations by comparison with the observed rates of inactivation of the enzyme. The study has shown that the enzyme active site appears to be flexible enough to allow the nucleophilic deactivation reactions of both the (R) and (S) forms of a model of the inhibitor to be catalysed by the Zn(II) cofactor of CPA.     

Development of efficient docking strategies and structure-activity relationship study of the c-Met type II inhibitors

c-Met is a transmembrane receptor tyrosine kinase and an important therapeutic target for anticancer drugs. In the present study, we systematically investigated the influence of a range of parameters on the correlation between experimental and calculated binding free energies of type II c-Met inhibitors. We especially focused on evaluating the impact of different force fields, binding energy calculation methods, docking protocols, conformation sampling strategies, and conformations of the binding site captured in several crystallographic structures. Our results suggest that the force fields, the protein flexibility, and the selected conformation of the binding site substantially influence the correlation coefficient, while the sampling strategies and ensemble docking only mildly affect the prediction accuracy. Structure-activity relationship study suggests that the structural determinants to the high binding affinity of the type II inhibitors originate from its overall linear shape, hydrophobicity, and two conserved hydrogen bonds. Results from this study will form the basis for establishing an efficient computational docking approach for c-Met type II inhibitors design.     

Combining docking, scoring and molecular field analyses to probe influenza neuraminidase-ligand interactions

In this project, several docking conditions, scoring functions and corresponding protein-aligned molecular field analysis (CoMFA) models were evaluated for a diverse set of neuraminidase (NA) inhibitors. To this end, a group of inhibitors were docked into the active site of NA. The docked structures were utilized to construct a corresponding protein-aligned CoMFA models by employing probe-based (H+, OH, CH3) energy grids and genetic partial least squares (G/PLS) statistical analysis. A total of 16 different docking configurations were evaluated, of which some succeeded in producing self-consistent and predictive CoMFA models. However, the best model coincided with docking the ionized ligands into the hydrated form of the binding site via PLP1 scoring function (r2LOO=0.735, r2PRESS against 24 test compounds=0.828). The highest-ranking CoMFA models were employed to probe NA-ligand interactions. Further validation by comparison with a co-crystallized ligand-NA crystallographic structure was performed. This combination of docking/scoring/CoMFA modeling provided interesting insights into the binding of different NA inhibitors.     

Conformational analysis and docking study of potent factor XIIIa inhibitors having a cyclopropenone ring

A conformational analysis and docking study of potent factor XIIIa inhibitors having a cyclopropenone ring were carried out in an attempt to obtain structural insight into the inhibition mechanism. First, stable conformers of the inhibitors alone were obtained from the conformational analysis by systematic search and molecular dynamics. Next, a binding form model of factor XIIIa was built based on an X-ray crystal structure of the enzyme. Finally, the docking study of the inhibitors into the model’s binding site was performed. From the resulting stable complex structures, it was found that the cyclopropenone ring fits the active site located at the base of the binding cavity with high complementarity. The carbonyl oxygen of the cyclopropenone ring formed a hydrogen bond to the indole NH group of Trp279 and the terminal carbon atom of the reactive C=C double bond was in close proximity to the sulfur atom of the catalytic residue, Cys314. This binding mode suggests a possible inhibition mechanism, whereby the cysteine residue reacts with the cyclopropenone ring of the inhibitor, forming an enzyme-ligand adduct. In addition, the higher interaction energies between factor XIIIa and the inhibitors alluded to the probable binding sites of the ligand side chain.