Particular interaction between pyrimethamine derivatives and quadruple mutant type dihydrofolate reductase of Plasmodium falciparum: CoMFA and quantum chemical calculations studies

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 = 0.702, S_press = 0.608, = 0.980, s = 0.156, and = 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.     

CoMFA of artemisinin derivatives: effect of location and size of lattice.

A CoMFA study of artemisinin derivatives with changes of the location and the number of lattice points was performed. The location of probe atoms in a large lattice has practically no effect on the cross-validated r(2) value (r(2)(cv)). The selection of only 18 probe atoms around the peroxide bond, considering the action mechanism of artemisinin, provided a less time-demanding and more reliable CoMFA model, which forecasts better than the large lattice model despite the lower r(2)(cv) value. Only 1 A displacement of the small lattice caused a reduction of cross-validated r(2) value of more than 50%, which indicates the lattice location played an important role in this small lattice model.     

CoMFA and HQSAR of acylhydrazide cruzain inhibitors

An approach combining CoMFA and HQSAR methods was used to describe QSAR models for a series of cruzain inhibitors having the acylhydrazide framework. A CoMFA study using two alignment orientations (I and II), three different probe atoms and changes of the lattice spacing (1 and 2 A) was performed. Alignment II and an sp3 probe carbon atom yielded good cross-validation (q2=0.688) employing lattice spacing of 1 A. The best HQSAR model was generated using atoms, bond, and connectivity as fragment distinction and fragment size default (4-5) showing similar cross-validated value of CoMFA (q2=0.689). Based upon the information derived from CoMFA and HQSAR, we have identified some key features that may be used to design new acylhydrazide derivatives that may be more potent cruzain inhibitors.     

Computer-aided molecular design of 1H-imidazole-2,4-diamine derivatives as potential inhibitors of <Emphasis Type="Italic">Plasmodium falciparum</Emphasis> DHFR enzyme

Design and discovery of new potential inhibitors of Plasmodium falciparum dihydrofolate reductase (PfDHFR), equally active against both the wild-type and mutant strains, is urgently needed. In this study, a computer-aided molecular design approach that involved ab initio molecular orbital and density functional theory calculations, along with molecular electrostatic potential analysis, and molecular docking studies was employed to design 15 1H-imidazole-2,4-diamine derivatives as potential inhibitors of PfDHFR enzyme. Visual inspection of the binding modes of the compounds demonstrated that they all interact, via H-bond interactions, with key amino acid residues (Asp54, Ileu/Leu164, Asn/Ser108 and Ile14) similar to those of WR99210 (3) in the active site of the enzymes used in the study. These interactions are known to be essential for enzyme inhibition. These compounds showed better or comparable binding affinities to that of the bound ligand (WR99210). In silico toxicity predictions, carried out using TOPKAT software, also indicated that the compounds are non-toxic.     

3-D-QSAR CoMFA and CoMSIA studies on tetrahydrofuroyl-L-phenylalanine derivatives as VLA-4 antagonists

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on tetrahydrofuroyl-L-phenylalanine derivatives as VLA-4 antagonists. The best CoMFA and CoMSIA models that were generated using atom based alignment from a training set of twenty five tetrahydrofuroyl-L-phenylalanine derivatives, are six-component models with good statistics; CoMFA: r(2)(cv)=0.366, r(2)=0.983, s=0.099, F=172.661 and PRESS=4.435; CoMSIA: r(2)(cv)=0.528, r(2)=0.995, s=0.054, F=577.87 and PRESS=3.563. Both of these 3-D-QSAR models were validated using a test set of eleven compounds, whose predicted pIC(50) values fall within one log unit of the actual pIC(50). The contour diagrams obtained for the various CoMFA and CoMSIA field contributions can be mapped back onto structural features to explain the activity trends of the molecules analysed. Based on the spatial arrangement of the various field contributions, novel molecules with improved activity can be designed.     

Structure-Guided Modification of Rhizomucor miehei Lipase for Production of Structured Lipids

To improve the performance of yeast surface-displayed Rhizomucor miehei lipase (RML) in the production of human milk fat substitute (HMFS), we mutated amino acids in the lipase substrate-binding pocket based on protein hydrophobicity, to improve esterification activity. Five mutants: Asn87Ile, Asn87Ile/Asp91Val, His108Leu/Lys109Ile, Asp256Ile/His257Leu, and His108Leu/Lys109Ile/Asp256Ile/His257Leu were obtained and their hydrolytic and esterification activities were assayed. Using Discovery Studio 3.1 to build models and calculate the binding energy between lipase and substrates, compared to wild-type, the mutant Asp256Ile/His257Leu was found to have significantly lower energy when oleic acid (3.97 KJ/mol decrease) and tripalmitin (7.55 KJ/mol decrease) were substrates. This result was in accordance with the esterification activity of Asp256Ile/His257Leu (2.37-fold of wild-type). The four mutants were also evaluated for the production of HMFS in organic solvent and in a solvent-free system. Asp256Ile/His257Leu had an oleic acid incorporation of 28.27% for catalyzing tripalmitin and oleic acid, and 53.18% for the reaction of palm oil with oleic acid. The efficiency of Asp256Ile/His257Leu was 1.82-fold and 1.65-fold that of the wild-type enzyme for the two reactions. The oleic acid incorporation of Asp256Ile/His257Leu was similar to commercial Lipozyme RM IM for palm oil acidolysis with oleic acid. Yeast surface-displayed RML mutant Asp256Ile/His257Leu is a potential, economically feasible catalyst for the production of structured lipids.     

Role of hydrogen bond networks and dynamics in positive and negative cooperative stabilization of a protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.     

CoMFA and HQSAR studies on 6,7-dimethoxy-4-pyrrolidylquinazoline derivatives as phosphodiesterase10A inhibitors

Phosphodiesterase10A (PDE10A) is an important enzyme with diverse biological actions in intracellular signaling systems, making it an emerging target for diseases such as schizophrenia, Huntington's disease, and diabetes mellitus. The objective of the current 3D QSAR study is to uncover some of the structural parameters which govern PDE10A inhibitory activity over PDE3A/B. Thus, comparative molecular field analysis (CoMFA) and hologram quantitative structure-activity relationship (HQSAR) studies were carried out on recently reported 6,7-dimethoxy-4-pyrrolidylquinazoline derivatives as PDE10A inhibitors. The best CoMFA model using atom-fit alignment approach with the bound conformation of compound 21 as the template yielded the steric parameter as a major contributor (nearly 70%) to the observed variations in biological activity. The best CoMFA model produced statistically significant results, with the cross-validated (r(cv)(2)) and conventional correlation (r(ncv)(2)) coefficients being 0.557 and 0.991, respectively, for the 21 training set compounds. Validation of the model by external set of six compounds yielded a high (0.919) predictive value. The CoMFA models of PDE10A and PDE3A/B activity were compared in order to address the selectivity issue of these inhibitors. The best HQSAR model for PDE10A was obtained with an r(cv)(2) of 0.704 and r(ncv)(2) of 0.902 using atoms, bonds, connections, chirality, donor, and acceptor as fragment distinction and default fragment size of 4-7 with three components for the 21 compounds. The HQSAR model predicted the external test-set of compounds well since a good agreement between the experimental and predicted values was verified. Taken together, the present QSAR models were found to accurately predict the PDE10A inhibitory activity of the test-set compounds and to yield reliable clues for further optimization of the quinazoline derivatives in the dataset.     

Modeling study of Rusticyanin-Cytochrome C(4) complex: an insight to possible H-bond mediated recognition and electron--transfer process

Rusticyanin (RCy) mediated transfer of electron to Cytochrome C(4) (Cytc(4)) from the extracellular Fe(+2) ion is primarily involved in the Thiobacillus ferrooxidans induced bio-leaching of pyrite ore and also in the metabolism of this acidophilic bacteria. The modeling studies have revealed the two possible mode of RCy-Cytc(4) complexation involving nearly the same stabilization energy approximately -15 x 10(3) kJ/mol, one through N-terminal Asp 15 and another -C terminal Glu 121 of Cytc(4) with the Cu-bonded His 143 of RCy. The Asp 15:His 143 associated complex (DH) of Cytc(4)-RCy was stabilized by the intermolecular H-bonds of the carboxyl oxygen atoms O(delta1) and O(delta2) of Asp 15 with the Nepsilon-atom of His 143 and O(b) atoms of Ala 8 and Asp 5 (of Cytc(4)) with the Thr 146 and Phe 51 (of RCy). But the other Glu 121:His 143 associated complex (EH) of Cytc(4)-RCy was stabilized by the H-bonding interaction of the oxygen atoms O(epsilon1) and O(epsilon2) of Glu 121 with the Nepsilon and Ogamma atoms of His 143 and Thr 146 of RCy. The six water molecules were present in the binding region of the two proteins in the energy minimized autosolvated DH and EH-complexes. The MD studies also revealed the presence of six interacting water molecules at the binding region between the two proteins in both the complexes. Several residues Gly 82 and 84, His 143 (RCy) were participated through the water mediated (W 389, W 430, W 413, W 431, W 373, and W 478) interaction with the Asp 15, Ile 82, and 62, Tyr 63 (Cytc(4)) in DH complex, whereas in EH complex the Phe 51, Asn 80, Tyr 146 (RCy) residues were observed to interact with Asn 108, Met 120, Glu 121 (of Cytc(4)) through the water molecules W 507, W 445, W 401, W 446, and W 440. The direct water mediated (W 478) interaction of His 143 (RCy) to Asp 15 (of Cytc(4)) was observed only in the DH complex but not in EH. These direct and water mediated H-bonding between the two respective proteins and the binding free energy with higher interacting buried surface area of the DH complex compare to other EH complex have indicated an alternative possibility of the electron transfer route through the interaction of His 143 of RCy and the N-terminal Asp 15 of Cytc(4).     

3-D QSAR studies on new dibenzyltin(IV) anticancer agents by comparative molecular field analysis (CoMFA).

Dibenzyltin(IV)dichloride and dibenzyltin(IV)diisothiocyanate derivatives with N,S-donor ligands show significant cytotoxic activities against human cancer cell lines, and are well compared to analogous dialkyltin(IV) derivatives. CoMFA models were generated for the first time for these organotin derivatives using the cytotoxic activities (against two human tumor cell lines, MCF-7, a mammary carcinoma and WiDr, a colon carcinoma) of 21 complexes. High r(2) and r(2)(cv) values for both CoMFA models indicate good predictive power for the models.     

Docking-based CoMFA and CoMSIA analyses of tetrahydro-β-carboline derivatives as type-5 phosphodiesterase inhibitors

Tetrahydro-β-carboline derivatives (THBCs) have been identified as a class of potent Type-5 Phosphodiesterase (PDE5) inhibitors, showing benefits for the treatment of erectile dysfunction and also bearing anticancer properties. A computational strategy based on molecular docking studies, followed by docking-based Comparative Molecular Fields Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA), has been used to elucidate the atomic details of the PDE5/THBC interactions and to identify the most important features impacting the THBC PDE5 inhibitory activity. The final CoMSIA model resulted to be the more predictive, showing r(ncv)(2) = 0.96, r(cv)(2) = 0.688, SEE = 0.248, F = 104.800, and r(2)(pred) = 0.78. The results allowed us to obtain useful information for the design of new THBC analogues, potentially acting as PDE5 inhibitors, and to predict their potency prior to synthesis.     

The role of tryptophan-48 in catalysis and binding of inhibitors of Plasmodium falciparum dihydrofolate reductase

Dihydrofolate reductases (DHFRs) from Plasmodium falciparum (Pf) and various species of both prokaryotic and eukaryotic organisms have a conserved tryptophan (Trp) at position 48 in the active site. The role in catalysis and binding of inhibitors of the conserved Trp48 of PfDHFR has been analysed by site-specific mutagenesis, enzyme kinetics and use of a bacterial surrogate system. All 19 mutant enzymes showed undetectable or very low specific activities, with the highest value of k(cat)/K(m) from the Tyr48 (W48Y) mutant (0.12 versus 11.94M(-1)s(-1)), of about 1% of the wild-type enzyme. The inhibition constants for pyrimethamine, cycloguanil and WR99210 of the W48Y mutants are 2.5-5.3 times those of the wild-type enzyme. All mutants, except W48Y, failed to support the growth of Escherichia coli transformed with the parasite gene in the presence of trimethoprim, indicating the loss of functional activity of the parasite enzyme. Hence, Trp48 plays a crucial role in catalysis and inhibitor binding of PfDHFR. Interestingly, W48Y with an additional mutation at Asn188Tyr (N188Y) was found to promote bacterial growth and yielded a higher amount of purified enzyme. However, the kinetic parameters of the purified W48Y+N188Y enzyme were comparable with W48Y and the binding affinities for DHFR inhibitors were also similar to the wild-type enzyme. Due to its conserved nature, Trp48 of PfDHFR is a potential site for interaction with antimalarial inhibitors which would not be compromised by its mutations.     

3D-QSAR studies of farnesyltransferase inhibitors: a comparative molecular field analysis approach

3D-QSAR analysis has been performed on a series of previously synthesized benzonitrile derivatives, which were screened as farnesyltransferase inhibitors, using comparative molecular field analysis (CoMFA) with partial least-square fit to predict the steric and electrostatic molecular field interactions for the activity. The CoMFA study was carried out using a training set of 34 compounds. The predictive ability of the model developed was assessed using a test set of eight compounds (r(pred)(2) as high as 0.770). The analyzed 3D-QSAR CoMFA model has demonstrated a good fit, having r(2) value of 0.991 and cross-validated coefficient q(2) value as 0.619. The analysis of CoMFA contour maps provided insight into the possible modification of the molecules for better activity.     

Production of an activated form of Bacillus stearothermophilus L-2-hydroxyacid dehydrogenase by directed evolution

Bacillus stearothermophillus lactate dehydrogenase (bsLDH) is activated in the presence of fructose 1,6 bisphosphate (FBP). The activator is expensive and representative of the sort of co-factor complications that are undesirable in industrial processes. Three rounds of random mutagenesis and screening produced a mutant (6A) which is almost fully activated in the absence of FBP. Wild-type bsLDH has a K(pyr)(M) of 5 mM in the absence of FBP but when activated (+FBP) the K(pyr)(M) drops to 0.05 mM. The mutant 6A has a K(pyr)(M) of 0.07 mM in the absence of FBP. 6A has three amino acid substitutions-R118C, Q203L and N307S-resulting in a 70-fold activation, none of the mutations are near the active site. The activation of wild type bsLDH is due to an FBP induced tetramerization of dimeric bsLDH bringing about a structural rearrangement of key active site residues. The most likely explanation for the activation of 6A is derived from the position of Q203L, which is at the dimer-dimer interface. The suggestion is that the hydrophilic to hydrophobic change has altered the dimer-tetramer equilibrium position towards that of the tetramer. What is significant is the activation of bsLDH by a subtle long range event produced by the 'blind' directed evolution approach.     

Site-directed mutagenesis of highly conserved amino acids in the first cytoplasmic loop of Drosophila Rh1 opsin blocks rhodopsin synthesis in the nascent state.

The cytoplasmic surface of Drosophila melanogaster Rh1 rhodopsin (ninaE) harbours amino acids which are highly conserved among G-protein-coupled receptors. Site-directed mutations which cause Leu81Gln or Asn86Ile amino acid substitutions in the first cytoplasmic loop of the Rh1 opsin protein, are shown to block rhodopsin synthesis in the nascent, glycosylated state from which the mutant opsin is degraded rapidly. In mutants Leu81Gln and Asn86Ile, only 20-30% and <2% respectively, of functional rhodopsins are synthesized and transported to the photoreceptive membrane. Thus, conserved amino acids in opsin's cytoplasmic surface are a critical factor in the interaction of opsin with proteins of the rhodopsin processing machinery. Photoreceptor cells expressing mutant rhodopsins undergo age-dependent degeneration in a recessive manner.     

Pharmacophore modeling, virtual screening and 3D-QSAR studies of 5-tetrahydroquinolinylidine aminoguanidine derivatives as sodium hydrogen exchanger inhibitors

Sodium hydrogen exchanger (SHE) inhibitor is one of the most important targets in treatment of myocardial ischemia. In the course of our research into new types of non-acylguanidine, SHE inhibitory activities of 5-tetrahydroquinolinylidine aminoguanidine derivatives were used to build pharmacophore and 3D-QSAR models. Genetic Algorithm Similarity Program (GASP) was used to derive a 3D pharmacophore model which was used in effective alignment of data set. Eight molecules were selected on the basis of structure diversity to build 10 different pharmacophore models. Model 1 was considered as the best model as it has highest fitness score compared to other nine models. The obtained model contained two acceptor sites, two donor atoms and one hydrophobic region. Pharmacophore modeling was followed by substructure searching and virtual screening. The best CoMFA model, representing steric and electrostatic fields, obtained for 30 training set molecules was statistically significant with cross-validated coefficient (q(2)) of 0.673 and conventional coefficient (r(2)) of 0.988. In addition to steric and electrostatic fields observed in CoMFA, CoMSIA also represents hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. CoMSIA model was also significant with cross-validated coefficient (q(2)) and conventional coefficient (r(2)) of 0.636 and 0.986, respectively. Both models were validated by an external test set of eight compounds and gave satisfactory prediction (r(pred)(2)) of 0.772 and 0.701 for CoMFA and CoMSIA models, respectively. This pharmacophore based 3D-QSAR approach provides significant insights that can be used to design novel, potent and selective SHE inhibitors.     

New potent 5-substituted benzofuroxans as inhibitors of Trypanosoma cruzi growth: quantitative structure-activity relationship studies

Benzofuroxan derivatives have been shown to inhibit the growth of Trypanosoma cruzi, the etiological agent of Chagas' disease. Therefore, 2D- and 3D-QSAR models of their in vitro antichagasic activity were developed. Six new derivatives were synthesized to complete a final set of 26 structurally diverse benzofuroxans. The 2D-QSAR model (r = 0.939, r(adj)(2) = 0.849) was generated using multiple regression analysis of tabulated substituents' physicochemical properties and indicator variables. In addition, a 3D-QSAR model (r(2) = 0.997, q(2) = 0.802) was obtained using a comparative molecular field analysis (CoMFA). Due to the well-known benzofuroxan tautomerism, in both approaches (2D- and 3D-QSAR) it was necessary to include an indicator variable to consider the N-oxide position (I(6)). This parameter was established using low-temperature NMR experiments. Both QSAR models identified the electrophilic character of the substituent alpha-atom as a requirement for activity. Further support was found using a density functional theory (DFT) approach.     

过氧化氢酶Ⅰ结构延伸突变改善酶热稳定性的初步研究

嗜热脂肪芽孢杆菌过氧化氢酶Ⅰ及其取代突变酶８８（２５）的延伸突变是指在酶蛋白的Ｃ末端上连接一段随机肽链,从而改变了酶蛋白的结构．研究结果表明,这种延伸突变的方法非常有效地提高了酶的热稳定性,并且随机肽链的疏水性与其对应的延伸突变体酶的热稳定性呈现一定的负相关性,即随机肽链的疏水性越高,对应的突变体酶热稳定性越低