In vitro and in vivo analyses of a Phe/Tyr switch controlling product specificity of histone lysine methyltransferases

The functional significance of mono-, di-, and tri-methylation of lysine residues within histone proteins is under investigation. Evidence from several model organisms suggests that different methylated states of H3 Lys(9) (H3K9) are generated by specific histone methyltransferases (MTases) to mark distinct types of silent chromatin. Sequence alignment of all histone lysine MTases with known product specificity suggested that a key residue in the active site determines how many methyl groups they add. We examined this possibility both in vitro and in vivo and found that a Phe at the position equivalent to Phe(281) of Neurospora crassa DIM-5 or Phe(1205) of human G9a allows the enzyme to perform di and tri-methylation, whereas a Tyr at this position is restrictive, inhibiting tri-methylation and thus yielding a mono- or di-MTase. Phe to Tyr mutants of both DIM-5 and G9a restrict product specificity in vitro and in vivo without compromising overall catalysis. These mutants were employed to probe the biological significance of mono-, di-, and tri-methylation of H3K9 in both mouse embryonic stem cells and N. crassa. G9a F1205Y, when expressed in G9a (-/-) embryonic stem cells, rescued only H3K9 mono-methylation, but not di-methylation, to wild-type levels yet silenced Mage-a gene expression. When expressed in dim-5 strains, DIM-5 F281Y generated significant levels of mono- and di-H3K9 methylation (which are not observed in wild type Neurospora) as well as tri-methyl H3K9. The altered DIM-5 rescued the growth defect characteristic of dim-5 N. crassa but did not fully rescue the gross DNA hypomethylation of dim-5 strains.     

Identification of Tyr(290(6.58)) of the human gonadotropin-releasing hormone (GnRH) receptor as a contact residue for both GnRH I and GnRH II: importance for high-affinity binding and receptor activation

Molecular modeling showed interactions of Tyr (290(6.58)) in transmembrane domain 6 of the GnRH receptor with Tyr (5) of GnRH I, and His (5) of GnRH II. The wild-type receptor exhibited high affinity for [Phe (5)]GnRH I and [Tyr (5)]GnRH II, but 127- and 177-fold decreased affinity for [Ala (5)]GnRH I and [Ala (5)]GnRH II, indicating that the aromatic ring in position 5 is crucial for receptor binding. The receptor mutation Y290F decreased affinity for GnRH I, [Phe (5)]GnRH I, GnRH II and [Tyr (5)]GnRH II, while Y290A and Y290L caused larger decreases, suggesting that both the para-OH and aromatic ring of Tyr (290(6.58)) are important for binding of ligands with aromatic residues in position 5. Mutating Tyr (290(6.58)) to Gln increased affinity for Tyr (5)-containing GnRH analogues 3-12-fold compared with the Y290A and Y290L mutants, suggesting a hydrogen-bond between Gln of the Y290Q mutant and Tyr (5) of GnRH analogues. All mutations had small effects on affinity of GnRH analogues that lack an aromatic residue in position 5. These results support direct interactions of the Tyr (290(6.58)) side chain with Tyr (5) of GnRH I and His (5) of GnRH II. Tyr (290(6.58)) mutations, except for Y290F, caused larger decreases in GnRH potency than affinity, indicating that an aromatic ring is important for the agonist-induced receptor conformational switch.     

Autoreduction of ferryl myoglobin: discrimination among the three tyrosine and two tryptophan residues as electron donors

Ferric myoglobin undergoes a two-electron oxidation in its reaction with H(2)O(2). One oxidation equivalent is used to oxidize Fe(III) to the Fe(IV) ferryl species, while the second is associated with a protein radical but is rapidly dissipated. The ferryl species is then slowly reduced back to the ferric state by unknown mechanisms. To clarify this process, the formation and stability of the ferryl forms of the Tyr --> Phe and Trp --> Phe mutants of recombinant sperm whale myoglobin (SwMb) were investigated. Kinetic studies showed that all the mutants react normally with H(2)O(2) to give the ferryl species. However, the rapid phase of ferryl autoreduction typical of wild-type SwMb was absent in the triple Tyr --> Phe mutant and considerably reduced in the Y103F and Y151F mutants, strongly implicating these two residues as intramolecular electron donors. Replacement of Tyr146, Trp7, or Trp14 did not significantly alter the autoreduction, indicating that these residues do not contribute to ferryl reduction despite the fact that Tyr146 is closer to the iron than Tyr151 or Tyr103. Furthermore, analysis of the fast phase of autoreduction in the dimer versus recovered monomer of the Tyr --> Phe mutant K102Q/Y103F/Y146F indicates that the Tyr151-Tyr151 cross-link is a particularly effective electron donor. The presence of an additional, slow phase of reduction in the triple Tyr --> Phe mutant indicates that alternative but normally minor electron-transfer pathways exist in SwMb. These results demonstrate that internal electron transfer is governed as much by the tyrosine pK(a) and oxidation potential as by its distance from the electron accepting iron atom.     

Mutational analysis of cytochrome b at the ubiquinol oxidation site of yeast complex III

The cytochrome bc1 complex is a dimeric enzyme of the inner mitochondrial membrane that links electron transfer from ubiquinol to cytochrome c by a protonmotive Q cycle mechanism in which ubiquinol is oxidized at one center in the enzyme, referred to as center P, and ubiquinone is rereduced at a second center, referred to as center N. To better understand the mechanism of ubiquinol oxidation, we have examined catalytic activities and pre-steady-state reduction kinetics of yeast cytochrome bc1 complexes with mutations in cytochrome b that we expected would affect oxidation of ubiquinol. We mutated two residues thought to be involved in proton conduction linked to ubiquinol oxidation, Tyr132 and Glu272, and two residues proposed to be involved in docking ubiquinol into the center P pocket, Phe129 and Tyr279. Substitution of Phe129 by lysine or arginine yielded a respiration-deficient phenotype and lipid-dependent catalytic activity. Increased bypass reactions were detectable for both variants, with F129K showing the more severe effects. Substitution with lysine leads to a disturbed coordination of a b heme as deduced from changes in the midpoint potential and the EPR signature. Removal of the aromatic side chain in position Tyr279 lowers the catalytic activity accompanied by a low level of bypass reactions. Pre-steady-state kinetics of the enzymes modified at Glu272 and Tyr132 confirmed the importance of their functional groups for electron transfer. Altered center N kinetics and activation of ubiquinol oxidation by binding of cytochrome c in the Y132F and E272D enzymes indicate long range effects of these mutations.     

Pressure-dependent ionization of Tyr 9 in glutathione S-transferase A1-1: contribution of the C-terminal helix to a "soft" active site.

The glutathione S-transferase (GST) isozyme A1-1 contains at its active site a catalytic tyrosine, Tyr9, which hydrogen bonds to, and stabilizes, the thiolate form of glutathione, GS-. In the substrate-free GST A1-1, the Tyr 9 has an unusually low pKa, approximately 8.2, for which the ionization to tyrosinate is monitored conveniently by UV and fluorescence spectroscopy in the tryptophan-free mutant, W21F. In addition, a short alpha-helix, residues 208-222, provides part of the GSH and hydrophobic ligand binding sites, and the helix becomes "disordered" in the absence of ligands. Here, hydrostatic pressure has been used to probe the conformational dynamics of the C-terminal helix, which are apparently linked to Tyr 9 ionization. The extent of ionization of Tyr 9 at pH 7.6 is increased dramatically at low pressures (p1/2 = 0.52 kbar), based on fluorescence titration of Tyr 9. The mutant protein W21F:Y9F exhibits no changes in tyrosine fluorescence up to 1.2 kbar; pressure specifically ionizes Tyr 9. The volume change, delta V, for the pressure-dependent ionization of Tyr 9 at pH 7.6, 19 degrees C, was -33 +/- 3 mL/mol. In contrast, N-acetyl tyrosine exhibits a delta V for deprotonation of -11 +/- 1 mL/mol, beginning from the same extent of initial ionization, pH 9.5. The pressure-dependent ionization is completely reversible for both Tyr 9 and N-acetyl tyrosine. Addition of S-methyl GSH converted the "soft" active site to a noncompressible site that exhibited negligible pressure-dependent ionization of Tyr 9 below 0.8 kbar. In addition, Phe 220 forms part of an "aromatic cluster" with Tyr 9 and Phe 10, and interactions among these residues were hypothesized to control the order of the C-terminal helix. The amino acid substitutions F220Y, F2201, and F220L afford proteins that undergo pressure-dependent ionization of Tyr 9 with delta V values of 31 +/- 2 mL/mol, 43 +/- 3 mL/mol, and 29 +/- 2 mL/mol, respectively. The p1/2 values for Tyr 9 ionization were 0.61 kbar, 0.41 kbar, and 0.46 kbar for F220Y, F220I, and F220L, respectively. Together, the results suggest that the C-terminal helix is conformationally heterogeneous in the absence of ligands. The conformations differ little in free energy, but they are significantly different in volume, and mutations at Phe 220 control the conformational distribution.     

Does "butyrylization" of acetylcholinesterase through substitution of the six divergent aromatic amino acids in the active center gorge generate an enzyme mimic of butyrylcholinesterase?

The active center gorge of human acetylcholinesterase (HuAChE) is lined by 14 aromatic residues, whereas in the closely related human butyrylcholinesterase (HuBChE) 3 of the aromatic active center residues (Phe295, Phe297, Tyr337) as well as 3 of the residues at the gorge entrance (Tyr72, Tyr124, Trp286) are replaced by aliphatic amino acids. To investigate whether this structural variability can account for the reactivity differences between the two enzymes, gradual replacement of up to all of the 6 aromatic residues in HuAChE by the corresponding residues in HuBChE was carried out. The affinities of the hexamutant (Y72N/Y124Q/W286A/F295L/F297V/Y337A) toward tacrine, decamethonium, edrophonium, huperzine A, or BW284C51 differed by about 5-, 80-, 170-, 25000-, and 17000-fold, respectively, from those of the wild-type HuAChE. For most of these prototypical noncovalent active center and peripheral site ligands, the hexamutant HuAChE displayed a reactivity phenotype closely resembling that of HuBChE. These results support the accepted view that the active center architectures of AChE and BChE differ mainly by the presence of a larger void space in BChE. Nevertheless, reactivity of the hexamutant HuAChE toward the substrates acetylthiocholine and butyrylthiocholine, or covalent ligands such as phosphonates and the transition state analogue m-(N,N,N-trimethylammonio)trifluoroacetophenone (TMTFA), is about 45-170-fold lower than that of HuBChE. Most of this reduction in reactivity can be related to the combined replacements of the three aromatic residues at the active center, Phe295, Phe297, and Tyr337. We propose that the hexamutant HuAChE, unlike BChE, is impaired in its capacity to accommodate certain tetrahedral species in the active center. This impairment may be related to the enhanced mobility of the catalytic histidine His447, which is observed in molecular dynamics simulations of the hexamutant and the F295L/F297V/Y337A HuAChE enzymes but not in the wild-type HuAChE.     

Role of Phe286 in the recognition mechanism of cyclomaltooligosaccharides (cyclodextrins) by Thermoactinomyces vulgaris R-47 alpha-amylase 2 (TVAII). X-ray structures of the mutant TVAIIs, F286A and F286Y, and kinetic analyses of the Phe286-replaced mutant TVAIIs

Phe286 located in the center of the active site of alpha-amylase 2 from Thermoactinomyces vulgaris R-47 (TVAII) plays an important role in the substrate recognition for cyclomaltooligosaccharides (cyclodextrins). The X-ray structures of mutant TVAIIs with the replacement of Phe286 by Ala (F286A) and Tyr (F286Y) were determined at 3.2 A resolution. Their structures have no significant differences from that of the wild-type enzyme. The kinetic analyses of Phe286-replaced variants showed that the variants with non-aromatic residues, Ala (F286A) and Leu (F286L), have lower enzymatic activities than those with aromatic residues, Tyr (F286Y) and Trp (F286W), and the replacement of Phe286 affects enzymatic activities for CDs more than those for starch.     

SET7/9 Catalytic Mutants Reveal the Role of Active Site Water Molecules in Lysine Multiple Methylation

SET domain lysine methyltransferases (KMTs) methylate specific lysine residues in histone and non-histone substrates. These enzymes also display product specificity by catalyzing distinct degrees of methylation of the lysine ϵ-amino group. To elucidate the molecular mechanism underlying this specificity, we have characterized the Y245A and Y305F mutants of the human KMT SET7/9 (also known as KMT7) that alter its product specificity from a monomethyltransferase to a di- and a trimethyltransferase, respectively. Crystal structures of these mutants in complex with peptides bearing unmodified, mono-, di-, and trimethylated lysines illustrate the roles of active site water molecules in aligning the lysine ϵ-amino group for methyl transfer with S-adenosylmethionine. Displacement or dissociation of these solvent molecules enlarges the diameter of the active site, accommodating the increasing size of the methylated ϵ-amino group during successive methyl transfer reactions. Together, these results furnish new insights into the roles of active site water molecules in modulating lysine multiple methylation by SET domain KMTs and provide the first molecular snapshots of the mono-, di-, and trimethyl transfer reactions catalyzed by these enzymes.     

The S2 subsites of cathepsins K and L and their contribution to collagen degradation

The exchange of residues 67 and 205 of the S2 pocket of human cysteine cathepsins K and L induces a permutation of their substrate specificity toward fluorogenic peptide substrates. While the cathepsin L-like cathepsin K (Tyr67Leu/Leu205Ala) mutant has a marked preference for Phe, the Leu67Tyr/Ala205Leu cathepsin L variant shows an effective cathepsin K-like preference for Leu and Pro. A similar turnaround of inhibition was observed by using specific inhibitors of cathepsin K [1-(N-Benzyloxycarbonyl-leucyl)-5-(N-Boc-phenylalanyl-leucyl)carbohydrazide] and cathepsin L [N-(4-biphenylacetyl)-S-methylcysteine-(D)-Arg-Phe-beta-phenethylamide]. Molecular modeling studies indicated that mutations alter the character of both S2 and S3 subsites, while docking calculations were consistent with kinetics data. The cathepsin K-like cathepsin L was unable to mimic the collagen-degrading activity of cathepsin K against collagens I and II, DQ-collagens I and IV, and elastin-Congo Red. In summary, double mutations of the S2 pocket of cathepsins K (Y67L/L205A) and L (L67Y/A205L) induce a switch of their enzymatic specificity toward small selective inhibitors and peptidyl substrates, confirming the key role of residues 67 and 205. However, mutations in the S2 subsite pocket of cathepsin L alone without engineering of binding sites to chondroitin sulfate are not sufficient to generate a cathepsin K-like collagenase, emphasizing the pivotal role of the complex formation between glycosaminoglycans and cathepsin K for its unique collagenolytic activity.     

Structurally conserved aromaticity of Tyr249 and Phe264 in helix 7 is important for toxicity of the Bacillus thuringiensis Cry4Ba toxin

Functional elements of the conserved helix 7 in the poreforming domain of the Bacillus thuringiensis Cry delta- endotoxins have not yet been clearly identified. Here, we initially performed alanine substitutions of four highly conserved aromatic residues, Trp(243), Phe(246), Tyr(249) and Phe(264), in helix 7 of the Cry4Ba mosquito-larvicidal protein. All mutant toxins were overexpressed in Escherichia coli as 130-kDa protoxins at levels comparable to the wild-type. Bioassays against Stegomyia aegypti mosquito larvae revealed that only W243A, Y249A or F264A mutant toxins displayed a dramatic decrease in toxicity. Further mutagenic analysis showed that replacements with an aromatic residue particularly at Tyr(249) and Phe(264) still retained the high-level toxin activity. In addition, a nearly complete loss in larvicidal activity was found for Y249L/F264L or F264A/ Y249A double mutants, confirming the involvement in toxicity of both aromatic residues which face towards the same direction. Furthermore, the Y249L/F264L mutant was found to be structurally stable upon toxin solubilisation and trypsin digestion, albeit a small change in the circular dichroism spectrum. Altogether, the present study provides for the first time an insight into the highly conserved aromaticity of Tyr(249) and Phe(264) within helix 7 playing an important role in larvicidal activity of the Cry4Ba toxin.     

Determinants of the substrate specificity of multidrug resistance protein 1: role of amino acid residues with hydrogen bonding potential in predicted transmembrane helix 17

Human multidrug resistance protein 1 (MRP1) confers resistance to many natural product chemotherapeutic agents and actively transports structurally diverse organic anion conjugates. We previously demonstrated that two hydrogen-bonding amino acid residues in the predicted transmembrane 17 (TM17) of MRP1, Thr(1242) and Trp(1246), were important for drug resistance and 17beta-estradiol 17-(beta-d-glucuronide) (E(2)17betaG) transport. To determine whether other residues with hydrogen bonding potential within TM17 influence substrate specificity, we replaced Ser(1233), Ser(1235), Ser(1237), Gln(1239), Thr(1241), and Asn(1245) with Ala and Tyr(1236) and Tyr(1243) with Phe. Mutations S1233A, S1235A, S1237A, and Q1239A had no effect on any substrate tested. In contrast, mutations Y1236F and T1241A decreased resistance to vincristine but not to VP-16, doxorubicin, and epirubicin. Mutation Y1243F reduced resistance to all drugs tested by 2-3-fold. Replacement of Asn(1245) with Ala also decreased resistance to VP-16, doxorubicin, and epirubicin but increased resistance to vincristine. This mutation also decreased E(2)17betaG transport approximately 5-fold. Only mutation Y1243F altered the ability of MRP1 to transport both leukotriene 4 and E(2)17betaG. Together with our previous results, these findings suggest that residues with side chain hydrogen bonding potential, clustered in the cytoplasmic half of TM17, participate in the formation of a substrate binding site.     

Characterization of oligomeric intermediates in alpha-synuclein fibrillation: FRET studies of Y125W/Y133F/Y136F alpha-synuclein

The aggregation of alpha-synuclein is believed to be a critical step in the etiology of Parkinson's disease. A variety of biophysical techniques were used to investigate the aggregation and fibrillation of alpha-synuclein in which one of the four intrinsic Tyr residues was replaced by Trp, and two others by Phe, in order to permit fluorescence resonance energy transfer (FRET) between residues 39 (Tyr) and 125 (Trp). The mutant Y125W/Y133F/Y136F alpha-synuclein (one Tyr, one Trp) showed fibrillation kinetics similar to that of the wild-type, as did the Y125F/Y133F/Y136F (one Tyr, no Trp) and Y39F/Y125W/Y133F/Y136F (no Tyr, one Trp) mutants. Time-dependent changes in FRET, Fourier transform infrared, Trp fluorescence, dynamic light-scattering and other probes, indicate the existence of a transient oligomer, whose population reaches a maximum at the end of the lag time. This oligomer, in which the alpha-synuclein is in a partially folded conformation, is subsequently converted into fibrils, and has physical properties that are distinct from those of the monomer and fibrils. In addition, another population of soluble oligomers was observed to coexist with fibrils at completion of the reaction. The average distance between Tyr39 and Trp125 decreases from 24.9A in the monomer to 21.9A in the early oligomer and 18.8A in the late oligomer. Trp125 remains solvent-exposed in both the oligomers and fibrils, indicating that the C-terminal domain is not part of the fibril core. No FRET was observed in the fibrils, due to quenching of Tyr39 fluorescence in the fibril core. Thus, aggregation of alpha-synuclein involves multiple oligomeric intermediates and competing pathways.     

Molecular determinants of xenobiotic metabolism: QM/MM simulation of the conversion of 1-chloro-2,4-dinitrobenzene catalyzed by M1-1 glutathione S-transferase

Modeling methods allow the identification and analysis of determinants of reactivity and specificity in enzymes. The reaction between glutathione and 1-chloro-2,4-dinitrobenzene (CDNB) is widely used as a standard activity assay for glutathione S-transferases (GSTs). It is important to understand the causes of differences between catalytic GST isoenzymes and the effects of mutations and genetic polymorphisms. Quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations have been performed here to investigate the addition of the glutathione anion to CDNB in the wild-type M1-1 GST isoenzyme from rat and in three single point mutant (Tyr6Phe, Tyr115Phe, and Met108Ala) M1-1 GST enzymes. We have developed a specifically parameterized QM/MM method (AM1-SRP/CHARMM22) to model this reaction by fitting to experimental heats of formation and ionization potentials. Free energy profiles were obtained from molecular dynamics simulations of the reaction using umbrella sampling and weighted histogram analysis techniques. The reaction in solution has also been simulated and is compared to the enzymatic reaction. The free energies are in excellent agreement with experimental results. Overall the results of the present study show that QM/MM reaction pathway analysis provides detailed insight into the chemistry of GST and can be used to obtain mechanistic insight into the effects of specific mutations on this catalytic process.     

A model of structure and catalysis for ketoreductase domains in modular polyketide synthases

A putative catalytic triad consisting of tyrosine, serine, and lysine residues was identified in the ketoreductase (KR) domains of modular polyketide synthases (PKSs) based on homology modeling to the short chain dehydrogenase/reductase (SDR) superfamily of enzymes. This was tested by constructing point mutations for each of these three amino acid residues in the KR domain of module 6 of the 6-deoxyerythronolide B synthase (DEBS) and determining the effect on ketoreduction. Experiments conducted in vitro with the truncated DEBS Module 6+TE (M6+TE) enzyme purified from Escherichia coli indicated that any of three mutations, Tyr --> Phe, Ser --> Ala, and Lys --> Glu, abolish KR activity in formation of the triketide lactone product from a diketide substrate. The same mutations were also introduced in module 6 of the full DEBS gene set and expressed in Streptomyces lividans for in vivo analysis. In this case, the Tyr --> Phe mutation appeared to completely eliminate KR6 activity, leading to the 3-keto derivative of 6-deoxyerythronolide B, whereas the other two mutations, Ser --> Ala and Lys --> Glu, result in a mixture of both reduced and unreduced compounds at the C-3 position. The results support a model analogous to SDRs in which the conserved tyrosine serves as a proton donating catalytic residue. In contrast to deletion of the entire KR6 domain of DEBS, which causes a loss in substrate specificity of the adjacent acyltransferase (AT) domain in module 6, these mutations do not affect the AT6 specificity and offer a potentially superior approach to KR inactivation for engineered biosynthesis of novel polyketides. The homology modeling studies also led to identification of amino acid residues predictive of the stereochemical nature of KR domains. Finally, a method is described for the rapid purification of engineered PKS modules that consists of a biotin recognition sequence C-terminal to the thioesterase domain and adsorption of the biotinylated module from crude extracts to immobilized streptavidin. Immobilized M6+TE obtained by this method was over 95% pure and as catalytically effective as M6+TE in solution.     

Tyr26 and Phe73 are essential for full biological activity of the Fd Gene 5 protein

Abstract Using site-saturation mutagenesis, we have established all possible amino acid substitutions at Tyr26 and Phe73 that are compatible with biological activity of the gene 5 protein in vivo. No substitutions were found at either site that gave rise to a fully functional gene 5 protein, indicating that these two amino acid residues are crucial. However, partial activity was found if either residue was replaced by another aromatic amino acid (Y26F, Y26W, F73Y, F73W). The results suggest that both Tyr26 and Phe73 are important for base stacking in the nucleoprotein complex. The functional consequences of the removal of the hydroxyl group from Tyr26 argue that this residue may, in addition, be involved in hydrogen bond formation to confer greater stability on the complex. In contrast, the addition of such a group to Phe73 reduces activity.     

Identification of tyrosine residues involved in ligand recognition by the phosphatidylinositol 3-kinase Src homology 3 domain: circular dichroism and UV resonance Raman studies.

Src homology 3 (SH3) domains are small noncatalytic protein modules capable of mediating protein-protein interactions. We previously demonstrated that the association of a ligand peptide RLP1 (RKLPPRPSK) causes environmental and structural changes of Trp55 and some of seven Tyr residues in the phosphatidylinositol 3-kinase (PI3K) SH3 domain by circular dichroism (CD) and 235-nm excited UV resonance Raman (UVRR) spectroscopies [Okishio, N., et al. (2000) Biopolymers 57, 208-217]. In this work, the affected Tyr residues were identified as Tyr12, Tyr14, and Tyr73 by the CD analysis of a series of mutants, in which every single Tyr residue was replaced by a Phe residue. Among these three residues, Tyr14 was found to be a main contributor to the UVRR spectral change upon the RLP1 binding. Interestingly, CD and UVRR analyses revealed that RLP1 associates with the Y14F and Y14H mutants in different ways. These results suggest that Tyr14 plays a crucial role in the ligand recognition, and the amino acid substitution at Tyr14 affects the mode of PI3K SH3-ligand interaction. Our findings give an insight into how SH3 domains can produce diversity and specificity to transduce signaling within cells.     

Novel urushiol derivatives as HDAC8 inhibitors: rational design,virtual screening,molecular docking and molecular dynamics studies

Three series of novel urushiol derivatives were designed by introducing a hydroxamic acid moiety into the tail of an alkyl side chain and substituents with differing electronic properties or steric bulk onto the benzene ring and alkyl side chain. The compounds’ binding affinity toward HDAC8 was screened by Glide docking. The highest-scoring compounds were processed further with molecular docking, MD simulations, and binding free energy studies to analyze the binding modes and mechanisms. Ten compounds had Glide scores of ?8.2 to ?10.2, which revealed that introducing hydroxy, carbonyl, amino, or methyl ether groups into the alkyl side chain or addition of –F, –Cl, sulfonamide, benzamido, amino, or hydroxy substituents on the benzene ring could significantly increase binding affinity. Molecular docking studies revealed that zinc ion coordination, hydrogen bonding, and hydrophobic interactions contributed to the high calculated binding affinities of these compounds toward HDAC8. MD simulations and binding free energy studies showed that all complexes possessed good stability, as characterized by low RMSDs, low RMSFs of residues, moderate hydrogen bonding and zinc ion coordination and low values of binding free energies. Hie147, Tyr121, Phe175, Hip110, Phe119, Tyr273, Lys21, Gly118, Gln230, Leu122, Gly269, and Gly107 contributed favorably to the binding; and Van der Waals and electrostatic interactions provided major contributions to the stability of these complexes. These results show the potential of urushiol derivatives as HDAC8 binding lead compounds, which have great therapeutic potential in the treatment of various malignancies, neurological disorders, and human parasitic diseases.     

Importance of exposed aromatic residues in chitinase B from Serratia marcescens 2170 for crystalline chitin hydrolysis

Chitinase B (ChiB) of S. marcescens has five exposed aromatic residues linearly aligned toward the catalytic cleft, Tyr481 and Trp479 in the C-terminal domain, and Trp252, Tyr240 and Phe190 in the catalytic domain. To determine the contribution of these residues to the hydrolysis of crystalline beta-chitin, site-directed mutagenesis, to replace them by alanine, was carried out. The Y481A, W479A, W252A, and Y240A mutations all decreased the binding activity and hydrolyzing activity toward beta-chitin microfibrils. Substitution of Trp residues affected the binding activity more severely than that of Tyr residues. The F190A mutation decreased neither the binding activity nor the hydrolyzing activity. None of the mutations decreased the hydrolyzing activity toward soluble substrates. These results suggest that ChiB hydrolyzes crystalline beta-chitin via a mechanism in which four exposed aromatic residues play important roles, similar to the mechanism of hydrolysis by ChiA of this bacterium, although the directions of hydrolysis of the two chitinases are opposite.     

Quaternary ammonium compounds as water channel blockers. Specificity, potency, and site of action

Excessive water uptake through Aquaporins (AQP) can be life-threatening and reversible AQP inhibitors are needed. Here, we determined the specificity, potency, and binding site of tetraethylammonium (TEA) to block Aquaporin water permeability. Using oocytes, externally applied TEA blocked AQP1/AQP2/AQP4 with IC50 values of 1.4, 6.2, and 9.8 microM, respectively. Related tetraammonium compounds yielded some (propyl) or no (methyl, butyl, or pentyl) inhibition. TEA inhibition was lost upon a Tyr to Phe amino acid switch in the external water pore of AQP1/AQP2/AQP4, whereas the water permeability of AQP3 and AQP5, which lack a corresponding Tyr, was not blocked by TEA. Consistent with experimental data, multi-nanosecond molecular dynamics simulations showed one stable binding site for TEA, but not tetramethyl (TMA), in AQP1, resulting in a nearly 50% water permeability inhibition, which was reduced in AQP1-Y186F due to effects on the TEA inhibitory binding region. Moreover, in the simulation TEA interacted with charged residues in the C (Asp128) and E (Asp185) loop, and the A(Tyr37-Asn42-Thr44) loop of the neighboring monomer, but not directly with Tyr186. The loss of TEA inhibition in oocytes expressing properly folded AQP1-N42A or -T44A is in line with the computationally predicted binding mode. Our data reveal that the molecular interaction of TEA with AQP1 differs and is about 1000-fold more effective on AQPs than on potassium channels. Moreover, the observed experimental and simulated similarities open the way for rational design and virtual screening for AQP-specific inhibitors, with quaternary ammonium compounds in general, and TEA in particular as a lead compound.     

Conservative substitutions in the hydrophobic core of Rhodobacter sphaeroides thioredoxin produce distinct functional effects.

The internal residue Phe 25 in Rhodobacter sphaeroides thioredoxin was changed to five amino acids (Ala, Val, Leu, Ile, Tyr) by site-directed mutagenesis, and the mutant proteins were characterized in vitro and in vivo using the mutant trxA genes in an Escherichia coli TrxA- background. The substitution F25A severely impaired the functional properties of the enzyme. Strains expressing all other mutations can grow on methionine sulfoxide with growth efficiencies of 45-60% that of the wild type at 37 degrees, and essentially identical at 42 degrees. At both temperatures, however, strains harboring the substitutions F25V and F25Y had lower growth rates and formed smaller colonies. In another in vivo assay, only the wild type and the F25I substitution allowed growth of phage T3/7 at 37 degrees, demonstrating that subtle modifications of the protein interior at position 25 Ile/Leu or Phe/Tyr) can produce significant biological effects. All F25 mutants were good substrates for E. coli thioredoxin reductase. Although turnover rates and apparent Km values were significantly lower for all mutants compared to the wild type, catalytic efficiency of thioredoxin reductase was similar for all substrates. Determination of the free energy of unfolding showed that the aliphatic substitutions (Val, Leu, Ile) significantly destabilized the protein, whereas the F25Y substitution did not affect protein stability. Thus, thermodynamic stability of R. sphaeroides thioredoxin variants is not correlated with the distinct functional effects observed both in vivo and in vitro.