Naturally occurring Phe151Leu substitution near a conserved folding module lowers stability of glutathione transferase P1–1

Glutathione transferases (GSTs) are a family of enzymes that detoxify electrophilic compounds, such as carcinogens or drugs, by conjugating them to glutathione. The enzymes have contributed to the understanding of protein structure, due to large differences in amino acid sequence within the family, yet similar architecture and folding. Our objective was to conduct a systematic survey of GSTP1 polymorphisms and their function. Nearly all variants detected were known polymorphisms: IVS4+13C>A; Ile105Val; Ala114Val; and g.2596T>C (Ser185Ser). However, we also found a novel Phe151Leu substitution in an African-American subject (1 out of 111). Kinetic parameters for the conjugation reaction with 1-chloro-2,4-dinitrobenzene (CDNB) were determined for the novel variant enzyme purified via heterologous expression in Escherichia coli. Five substrates were used for measurement of specific activities, including isothiocyanate compounds that occur in cruciferous vegetables (benzylisothiocyanate, phenethylisothiocyanate, and sulforaphane). Such isothiocyanate substrates are potential cancer chemopreventive agents that are conjugated by GSTs. No major change in kinetic parameters was observed. However, the half-life at 50 °C of the Leu 151 enzyme was reduced to 12 min, as compared to 28 min for the Phe 151 enzyme. Residue 151 is located at the N-terminus of helix α6 in GST motif II, surrounded by hydrophobic residues, and near the conserved “hydrophobic staple” and N-capping box motifs. These local structural elements aid in formation of helix α6 and promote proper folding and protein stability. Analysis of the three-dimensional structure showed that substitution of Phe 151 with Leu produces a hydrophobic cavity in the GSTP1 core, thereby destabilizing its structure. Phe151Leu represents one of the first-described allelic variations in a protein folding motif.     

A conserved "hydrophobic staple motif" plays a crucial role in the refolding of human glutathione transferase P1-1

The specific (i, i+5) hydrophobic staple interaction involving a helix residue and a second residue located in the turn preceding the helix is a recurrent motif at the N terminus of alpha-helices. This motif is strictly conserved in the core of all soluble glutathione transferases (GSTs) as well as in other protein structures. Human GSTP1-1 variants mutated in amino acid Ile(149) and Tyr(154) of the hydrophobic staple motif of the alpha6-helix were analyzed. In particular, a double mutant cycle analysis has been performed to evaluate the role of the hydrophobic staple motif in the refolding process. The results show that this local interaction, by restricting the number of conformations of the alpha6-helix relative to the alpha1-helix, favors the formation of essential interdomain interactions and thereby accelerates the folding process. Thus, for the first time it is shown that the hydrophobic staple interaction has a role in the folding process of an intact protein. In P(i) class GSTs, Tyr(154) appears to be of particular structural importance, since it interacts with conserved residues Leu(21), Asp(24), and Gln(25) of the adjacent alpha1-helix which contributes to the active site. Human GSTP1-1 variants L21A and Y154F have also been analyzed in order to distinguish the role of interdomain interactions from that of the hydrophobic staple. The experimental results reported here suggest that the strict conservation of the hydrophobic staple motif reflects an evolutionary pressure for proteins to fold rapidly.     

Structures of thermolabile mutants of human glutathione transferase P1-1

An N-capping box motif (Ser/Thr-Xaa-Xaa-Asp) is strictly conserved at the beginning of helix alpha6 in the core of virtually all glutathione transferases (GST) and GST-related proteins. It has been demonstrated that this local motif is important in determining the alpha-helical propensity of the isolated alpha6-peptide and plays a crucial role in the folding and stability of GSTs. Its removal by site-directed mutagenesis generated temperature-sensitive folding mutants unable to refold at physiological temperature (37 degrees C). In the present work, variants of human GSTP1-1 (S150A and D153A), in which the capping residues have been substituted by alanine, have been generated and purified for structural analysis. Thus, for the first time, temperature-sensitive folding mutants of an enzyme, expressed at a permissive temperature, have been crystallized and their three-dimensional structures determined by X-ray crystallography. The crystal structures of human pi class GST temperature-sensitive mutants provide a basis for understanding the structural origin of the dramatic effects observed on the overall stability of the enzyme at higher temperatures upon single substitution of a capping residue.     

Identification of coding polymorphisms in human circadian rhythm genes PER1, PER2, PER3, CLOCK, ARNTL, CRY1, CRY2 and TIMELESS in a multi-ethnic screening panel.

STUDY OBJECTIVE: In this study, the exonic regions of the circadian rhythm genes PER1, PER2, PER3, CLOCK, ARNTL, CRY1, CRY2 and TIMELESS were re-sequenced and coding changes identified in a panel of 95 individuals varying in ethnicity. STUDY PARTICIPANTS: DNA screening panel consisting of 95 DNA samples (17 American Caucasians, 17 African Americans, 8 Ashkenazi Jews, 8 Chinese, 8 Japanese, 5 Mexican Indians, 8 Mexicans, 8 Northern Europeans, 8 Puerto Ricans, and 8 South Americans) selected from the Coriell Institute Human Variation Panel. RESULTS: In addition to coding changes already identified in the database dbSNP, novel coding changes were identified, including PER1: Pro37Ser, Pro351Ser, Gln988Pro, Ala998Thr; PER2: Leu83Arg, Leu157Leu, Thre174Ile, Phe400Phe, Pro822Pro, Ala828Thr, Ala861Val, Phe876Leu, Val883Met, Val903Ile, Ala923Pro; PER3: Pro67Pro, Val90Ile, His638His, Ala820Ala, Leu929Leu; ARNTL: Arg166Gln, Ser459Phe; CLOCK: Ala34Ala, Ser208Cys, Phe233Phe, Ser632Thr, Ser816Ser; TIMELESS: Met870Val and CRY2: His35His. No coding polymorphisms were identified in CRY1. CONCLUSIONS: Considerable genetic variation occurs within the coding region of the genes regulating circadian rhythm. Many of the non-synonymous coding polymorphisms could affect protein structure/function with the potential to affect molecular regulation of the sleep/wake cycle. Many of the potential functional effects could be ethnic group specific.     

Solution Structure of the N-terminal RNP Domain of U1A Protein: The Role of C-terminal Residues in Structure Stability and RNA Binding

The solution structure of a fragment of the human U1A spliceosomal protein containing residues 2 to 117 (U1A117) determined using multi-dimensional heteronuclear NMR is presented. The C-terminal region of the molecule is considerably more ordered in the free protein than thought previously and its conformation is different from that seen in the crystal structure of the complex with U1 RNA hairpin II. The residues between Asp90 and Lys98 form an α-helix that lies across the β-sheet, with residues Ile93, Ile94 and Met97 making contacts with Leu44, Phe56 and Ile58. This interaction prevents solvent exposure of hydrophobic residues on the surface of the β-sheet, thereby stabilising the protein. Upon RNA binding, helix C moves away from this position, changing its orientation by 135° to allow Tyr13, Phe56 and Gln54 to stack with bases of the RNA, and also allowing Leu44 to contact the RNA. The new position of helix C in the complex with RNA is stabilised by hydrophobic interactions from Ile93 and Ile94 to Ile58, Leu 41, Val62 and His10, as well as a hydrogen bond between Ser91 and Thr11. The movement of helix C mainly involves changes in the main-chain torsion angles of Thr89, Asp90 and Ser91, the helix thereby acting as a "lid" over the RNA binding surface.     

Polymorphisms of glutathione S-transferases M1, T1, P1 and A1 genes in the Tunisian population: an intra and interethnic comparative approach

Genetic polymorphisms in glutathione S-transferases (GSTs) genes might influence the detoxification activities of the enzymes predisposing individuals to cancer risk. Owing to the presence of these genetic variants, inter-individual and ethnic differences in GSTs detoxification capacity have been observed in various populations. Therefore, the present study was performed to determine the prevalence GSTM1 0/0, GSTT1 0/0, GSTP1 Ile(105)Val, and GSTA1 A/B polymorphisms in 154 healthy individuals from South Tunisia, and to compare them with those observed in North and Centre Tunisian populations and other ethnic groups. GSTM1 and GSTT1 polymorphisms were analyzed by a Multiplex-PCR approach, whereas GSTP1 and GSTA1 polymorphisms were examined by PCR-RFLP. The frequencies of GSTM10/0 and GSTT1 0/0 genotypes were 53.9% and 27.9%, respectively. The genotype distribution of GSTP1 was 47.4% (Ile/Ile), 40.9% (Ile/Val), and 11.7% (Val/Val). For GSTA1, the genotype distribution was 24.7% (A/A), 53.9% (A/B), and 21.4% (B/B). The combined genotypes distribution of GSTM1, GSTT1, GSTP1 and GSTA1 polymorphisms showed that thirty one of the 36 possible genotypes were present in our population; eight of them have a frequency greater than 5%. To the best of our knowledge, this is the first report of GSTs polymorphisms in South Tunisian population. Our findings demonstrate the impact of ethnicity and reveal a characteristic pattern for Tunisian population. The molecular studies in these enzymes provide basis for further epidemiological investigations in the population where these functional polymorphisms alter therapeutic response and act as susceptibility markers for various clinical conditions.     

Highly restricted distributions of hydrophobic and charged amino acids in longitudinal quadrants of alpha-helices

Helix formation in folding proteins is stabilized by binding of recurrent hydrophobic side chains in one longitudinal quadrant against the locally most hydrophobic region of the protein. To test this hypothesis, we fitted sequences of 247 alpha-helices of 55 proteins to the circular (infinite) template (symbol; see text) to maximize the strip-of-helix hydrophobicity index (the mean hydrophobicity of residues in (symbol; see text) positions). These template-predicted configurations closely matched crystallographic structures in 87% of four- or five-turn helices compared. We determined the longitudinal quadrant distributions of amino acids in the template-fitted, sheet projections of alpha-helices with respect to the best longitudinal, hydrophobic strip on each helix and to the N and C termini, interiors, and entire helices. Amino acids Leu, Ile, Val, and Phe were concentrated in one longitudinal quadrant (p less than 0.001). Lys, Arg, Asp, and Glu were not in the quadrant of Leu, Ile, Val, and Phe (p less than 0.001). Significant quadrant distributions for other amino acids and for termini of the helices were also found.     

Substrate specificity of the DnaK chaperone determined by screening cellulose-bound peptide libraries.

Hsp70 chaperones assist protein folding by ATP-dependent association with linear peptide segments of a large variety of folding intermediates. The molecular basis for this ability to differentiate between native and non-native conformers was investigated for the DnaK homolog of Escherichia coli. We identified binding sites and the recognition motif in substrates by screening 4360 cellulose-bound peptides scanning the sequences of 37 biologically relevant proteins. DnaK binding sites in protein sequences occurred statistically every 36 residues. In the folded proteins these sites are mostly buried and in the majority found in beta-sheet elements. The binding motif consists of a hydrophobic core of four to five residues enriched particularly in Leu, but also in Ile, Val, Phe and Tyr, and two flanking regions enriched in basic residues. Acidic residues are excluded from the core and disfavored in flanking regions. The energetic contribution of all 20 amino acids for DnaK binding was determined. On the basis of these data an algorithm was established that predicts DnaK binding sites in protein sequences with high accuracy.     

Modulation of the affinity of the single-stranded DNA-binding protein of Escherichia coli (E. coli SSB) to poly(dT) by site-directed mutagenesis

A vector for site-directed mutagenesis and overproduction of the Escherichia coli single-stranded-DNA-binding protein (E. coli SSB) was constructed. An E. coli strain carrying this vector produces up to 400 mg pure protein from 25 g wet cells. The vector was used to mutate specifically the Phe60 residue of E. coli SSB. Phe60 had been proposed to be located near the single-stranded-DNA-binding site. Substitution of the Phe60 residue by Val, Ser, Leu, His, Tyr and Trp gave proteins with no or only minor conformational changes, as detected by NMR spectroscopy. The affinity of the mutant E. coli SSB proteins for single-stranded DNA decreased in the order Trp greater than Phe (wild-type) greater than Tyr greater than Leu greater than His greater than Val greater than Ser, leading to the conclusion that position 60 is a site of hydrophobic interaction of the protein with DNA.     

Importance of secondary structural specificity determinants in protein folding: insertion of a native beta-sheet sequence into an alpha-helical coiled-coil

To examine how a short secondary structural element derived from a native protein folds when in a different protein environment, we inserted an 11-residue beta-sheet segment (cassette) from human immunoglobulin fold, Fab new, into an alpha-helical coiled-coil host protein (cassette holder). This de novo design protein model, the structural cassette mutagenesis (SCM) model, allows us to study protein folding principles involving both short- and long-range interactions that affect secondary structure stability and conformation. In this study, we address whether the insertion of this beta-sheet cassette into the alpha-helical coiled-coil protein would result in conformational change nucleated by the long-range tertiary stabilization of the coiled-coil, therefore overriding the local propensity of the cassette to form beta-sheet, observed in its native immunoglobulin fold. The results showed that not only did the nucleating helices of the coiled-coil on either end of the cassette fail to nucleate the beta-sheet cassette to fold with an alpha-helical conformation, but also the entire chimeric protein became a random coil. We identified two determinants in this cassette that prevented coiled-coil formation: (1) a tandem dipeptide NN motif at the N-terminal of the beta-sheet cassette, and (2) the hydrophilic Ser residue, which would be buried in the hydrophobic core if the coiled-coil structure were to fold. By amino acid substitution of these helix disruptive residues, that is, either the replacement of the NN motif with high helical propensity Ala residues or the substitution of Ser with Leu to enhance hydrophobicity, we were able to convert the random coil chimeric protein into a fully folded alpha-helical coiled-coil. We hypothesized that this NN motif is a "secondary structural specificity determinant" which is very selective for one type of secondary structure and may prevent neighboring residues from adopting an alternate protein fold. These sequences with secondary structural specificity determinants have very strong local propensity to fold into a specific secondary structure and may affect overall protein folding by acting as a folding initiation site.     

The effect of changing the hydrophobic S1' subsite of thermolysin-like proteases on substrate specificity.

The hydrophobic S1' subsite is one of the major determinants of the substrate specificity of thermolysin and related M4 family proteases. In the thermolysin-like protease (TLP) produced by Bacillus stearothermophilus (TLP-ste), the hydrophobic S1' subsite is mainly formed by Phe130, Phe133, Val139 and Leu202. In the present study, we have examined the effects of replacing Leu202 by smaller (Gly, Ala, Val) and larger (Phe, Tyr) hydrophobic residues. The mutational effects showed that the wild-type S1' pocket is optimal for binding leucine side chains. Reduction of the size of residue 202 resulted in a higher efficiency towards substrates with Phe in the P1' position. Rather unexpectedly, the Leu202-->Phe and Leu202-->Tyr mutations, which were expected to decrease the size of the S1' subsite, resulted in a large increase in activity towards dipeptide substrates with Phe in the P1' position. This is probably due to the fact that 202Phe and 202Tyr adopt a second possible rotamer that opens up the subsite compared to Leu202, and also favours interactions with the substrate. To validate these results, we constructed variants of thermolysin with changes in the S1' subsite. Thermolysin and TLP-ste variants with identical S1' subsites were highly similar in terms of their preference for Phe vs. Leu in the P1' position.     

Structural basis for the N‐degron specificity of ClpS1 from Arabidopsis thaliana

The N‐degron pathway determines the half‐life of proteins in both prokaryotes and eukaryotes by precisely recognizing the N‐terminal residue (N‐degron) of substrates. ClpS proteins from bacteria bind to substrates containing hydrophobic N‐degrons (Leu, Phe, Tyr, and Trp) and deliver them to the caseinolytic protease system ClpAP. This mechanism is preserved in organelles such as mitochondria and chloroplasts. Bacterial ClpS adaptors bind preferentially to Leu and Phe N‐degrons; however, ClpS1 from Arabidopsis thaliana (AtClpS1) shows a difference in that it binds strongly to Phe and Trp N‐degrons and only weakly to Leu. This difference in behavior cannot be explained without structural information due to the high sequence homology between bacterial and plant ClpS proteins. Here, we report the structure of AtClpS1 at 2.0 Å resolution in the presence of a bound N‐degron. The key determinants for α‐amino group recognition are conserved among all ClpS proteins, but the α3‐helix of eukaryotic AtClpS1 is significantly shortened, and consequently, a loop forming a pocket for the N‐degron is moved slightly outward to enlarge the pocket. In addition, amino acid replacement from Val to Ala causes a reduction in hydrophobic interactions with Leu N‐degron. A combination of the fine‐tuned hydrophobic residues in the pocket and the basic gatekeeper at the entrance of the pocket controls the N‐degron selectivity of the plant ClpS protein.     

Helix-forming tendencies of nonpolar amino acids predicted by Monte Carlo simulated annealing

Monte Carlo simulated annealing is applied to the study of the α-helix-forming tendencies of seven nonpolar amino acids, Ala, Leu, Met, Phe, Ile, Val, and Gly. Homooligomers of 10 amino acids are used and the helix tendency is calculated by folding α-helicies from completely random initial conformations. The results of the simulation imply that Met, Ala, and Leu are helix formers and that Val, Ile, and Gly are helix breakers, while Phe comes in between the two groups. The differences between helix formers and breakers turned out to be large in agreement with the recent experiments with short peptides. It is argued from the energy distributions of the obtained conformations that the helix tendency is small for the helix breakers because of steric hindrance of side chains. Homoglycine is shown to favor a random coil conformation. The β-strand tendencies of the same homooligomers are also considered, and they are shown to agree with the frequencies of amino acids in β-sheet from the protein data base. © 1994 Wiley-Liss, Inc.     

Amino acid sequence of a protease inhibitor isolated from Sarcophaga bullata determined by mass spectrometry.

The amino acid sequence of a protease inhibitor isolated from the hemolymph of Sarcophaga bullata larvae was determined by tandem mass spectrometry. Homology considerations with respect to other protease inhibitors with known primary structures assisted in the choice of the procedure followed in the sequence determination and in the alignment of the various peptides obtained from specific chemical cleavage at cysteines and enzyme digests of the S. bullata protease inhibitor. The resulting sequence of 57 residues is as follows: Val Asp Lys Ser Ala Cys Leu Gln Pro Lys Glu Val Gly Pro Cys Arg Lys Ser Asp Phe Val Phe Phe Tyr Asn Ala Asp Thr Lys Ala Cys Glu Glu Phe Leu Tyr Gly Gly Cys Arg Gly Asn Asp Asn Arg Phe Asn Thr Lys Glu Glu Cys Glu Lys Leu Cys Leu.     

Contribution of glycine 146 to a conserved folding module affecting stability and refolding of human glutathione transferase p1-1

In human glutathione transferase P1-1 (hGSTP1-1) position 146 is occupied by a glycine residue, which is located in a bend of a long loop that together with the alpha6-helix forms a substructure (GST motif II) maintained in all soluble GSTs. In the present study G146A and G146V mutants were generated by site-directed mutagenesis in order to investigate the function played by this conserved residue in folding and stability of hGSTP1-1. Crystallographic analysis of the G146V variant, expressed at the permissive temperature of 25 degrees C, indicates that the mutation causes a substantial change of the backbone conformation because of steric hindrance. Stability measurements indicate that this mutant is inactivated at a temperature as low as 32 degrees C. The structure of the G146A mutant is identical to that of the wild type with the mutated residue having main-chain bond angles in a high energy region of the Ramachandran plot. However even this Gly --> Ala substitution inactivates the enzyme at 37 degrees C. Thermodynamic analysis of all variants confirms, together with previous findings, the critical role played by GST motif II for overall protein stability. Analysis of reactivation in vitro indicates that any mutation of Gly-146 alters the folding pathway by favoring aggregation at 37 degrees C. It is hypothesized that the GST motif II is involved in the nucleation mechanism of the protein and that the substitution of Gly-146 alters this transient substructure. Gly-146 is part of the buried local sequence GXXh(T/S)XXDh (X is any residue and h is a hydrophobic residue), conserved in all GSTs and related proteins that seems to behave as a characteristic structural module important for protein folding and stability.     

Functional role of the lock and key motif at the subunit interface of glutathione transferase p1-1

The glutathione transferases (GSTs) represent a superfamily of dimeric proteins. Each subunit has an active site, but there is no evidence for the existence of catalytically active monomers. The lock and key motif is responsible for a highly conserved hydrophobic interaction in the subunit interface of pi, mu, and alpha class glutathione transferases. The key residue, which is either Phe or Tyr (Tyr(50) in human GSTP1-1) in one subunit, is wedged into a hydrophobic pocket of the other subunit. To study how an essentially inactive subunit influences the activity of the neighboring subunit, we have generated the heterodimer composed of subunits from the fully active human wild-type GSTP1-1 and the nearly inactive mutant Y50A obtained by mutation of the key residue Tyr(50) to Ala. Although the key residue is located far from the catalytic center, the k(cat) value of mutant Y50A decreased about 1300-fold in comparison with the wild-type enzyme. The decrease of the k(cat) value of the heterodimer by about 27-fold rather than the expected 2-fold in comparison with the wild-type enzyme indicates that the two active sites of the dimeric enzyme work synergistically. Further evidence for cooperativity was found in the nonhyperbolic GSH saturation curves. A network of hydrogen-bonded water molecules, found in crystal structures of GSTP1-1, connects the two active sites and the main chain carbonyl group of Tyr(50), thereby offering a mechanism for communication between the two active sites. It is concluded that a subunit becomes catalytically competent by positioning the key residue of one subunit into the lock pocket of the other subunit, thereby stabilizing the loop following the helix alpha2, which interacts directly with GSH.     

Transferase S-glutathione class pi gene (GSTP1) polymorphism in thyroid cancer patients

INTRODUCTION: One of the potential genes which can increase the risk of cancer is GSTP1 gene. It encodes enzyme called glutathione S-transferase pi class, which is involved in the detoxification of a variety of potential carcinogenic compounds. Polymorphism in this gene can cause the amino acid substitution. This substitution, close to the substrate binding site, changes the enzymatic activity for particular substrates and subsequently increases the risk of carcinogenesis. The aim of this study was to evaluate the function of GSTP1 polymorphism in thyroid cancer and possible association between GSTP1 polymorphism and age at diagnosis. MATERIAL AND METHODS: 103 patients with differentiated thyroid cancer (DTC) and 53 individuals from control group were examined using PCR-RFLP. RESULTS: Statistically insignificant association of studied polymorphisms with thyroid cancer was observed. Comparison of allele frequency between cases and control groups revealed the presence of risk alleles. For the first polymorphism Ile OR = 1.257; 95% CI [0.792-1.997] (p = 0.332), and for the second one Val OR = 1.283; 95% CI [0.6260-2.631] (p = 0.495). The presence of Val/Val (c.313A>G) led to a significant earlier age of onset as compared with other genotypes (p < 0.05). Mean age at diagnosis for Val/Val genotype was 41.1 +/- 15.2, and for Ile/Val + Ile/Ile reached 48.9 +/- 13.2. There was no association between age and genotype for c.341C>T polymorphism. CONCLUSIONS: Statistically insignificant association of GSTP1 gene polymorphism with thyroid cancer was observed in studied group of patients. The Val/Val genotype for c.313A>G polymorphism led to earlier age of tumour diagnosis as compared with other genotypes.     

Actin-inhibition and folding of vertebrate deoxyribonuclease I are affected by mutations at residues 67 and 114

Amino acid (aa) residues (Val-67 and Ala-114) have been suggested as being mainly responsible for actin-binding in human and bovine deoxyribonucleases I (DNase I). This study presents evidence of these two aa mutational mechanisms, not only for actin-binding but also for folding of DNase I in mammals, reptiles and amphibians. Human and viper snake (Agkistrodon blomhoffii) enzymes are inhibited by actin, whereas porcine, rat snake (Elaphe quadrivirgata), and African clawed frog (Xenopus laevis) enzymes are not. To investigate the role of aa at 67, mutants of rat snake (Ile67Val) and viper snake (Val67Ile) enzymes were constructed. After substitution, the rat snake was inhibited by actin, while the viper snake was not. For the role of aa at 114, mutants of viper snake (Phe114Ala), rat snake (Phe114Ala), African clawed frog (Phe114Ala), and porcine (Ser114Ala/Ser114Phe) enzymes were constructed. Strikingly, the substitute mutants for viper snake, rat snake and African clawed frog expressed no protein. The porcine (Ser114Ala) enzyme was inhibited by actin, but not the porcine (Ser114Phe) enzyme. These results suggest that Val-67 may be essential for actin-binding, that Phe-114 may be related to the folding of DNase I in reptiles and amphibians, and that Ala-114 may be indispensable for actin-binding in mammals.     

Structure prediction and R115866 binding study of human CYP26A1: homology modelling,fold recognition,molecular docking and MD simulations

Two three-dimensional (3D) models of human cytochrome P450 26A1 (CYP26A1) were constructed using the programs Modeller and Sybyl-GeneFold, respectively. After refinement by molecular mechanics and molecular dynamics (MD) simulations, the two models were validated by structure analysis-validation online server. Subsequently, a flexible docking study was performed on the model constructed by GeneFold with the potent and specific inhibitor R115866 to examine the enzyme–inhibitor interactions. From the docking results, we can see R115866 interacts with amino acid residues at the active site by multiple hydrophobic interactions including the side chains of His111, Trp112, Ser115, Val116, Leu125, Ser126, Leu221, Phe222, Glu296, Phe299, Gly300, Glu303, Thr304, Pro371 and the cofactor heme. Trp112 and Thr304 form hydrogen bonds with R115866 and play important roles in stabilising the complex. This constructed CYP26A1 model may provide an opportunity to understand the action mode of the enzyme and could be useful in designing novel retinoic acid metabolism blocking agents (RAMBAs).     

Halide-stabilizing residues of haloalkane dehalogenases studied by quantum mechanic calculations and site-directed mutagenesis

Haloalkane dehalogenases catalyze cleavage of the carbon-halogen bond in halogenated aliphatic compounds, resulting in the formation of an alcohol, a halide, and a proton as the reaction products. Three structural features of haloalkane dehalogenases are essential for their catalytic performance: (i) a catalytic triad, (ii) an oxyanion hole, and (iii) the halide-stabilizing residues. Halide-stabilizing residues are not structurally conserved among different haloalkane dehalogenases. The level of stabilization of the transition state structure of S(N)2 reaction and halide ion provided by each of the active site residues in the enzymes DhlA, LinB, and DhaA was quantified by quantum mechanic calculations. The residues that significantly stabilize the halide ion were assigned as the primary (essential) or the secondary (less important) halide-stabilizing residues. Site-directed mutagenesis was conducted with LinB enzyme to confirm location of its primary halide-stabilizing residues. Asn38Asp, Asn38Glu, Asn38Phe, Asn38Gln, Trp109Leu, Phe151Leu, Phe151Trp, Phe151Tyr, and Phe169Leu mutants of LinB were constructed, purified, and kinetically characterized. The following active site residues were classified as the primary halide-stabilizing residues: Trp125 and Trp175 of DhlA; Asn38 and Trp109 of LinB; and Asn41 and Trp107 of DhaA. All these residues make a hydrogen bond with the halide ion released from the substrate molecule, and their substitution results in enzymes with significantly modified catalytic properties. The following active site residues were classified as the secondary halide-stabilizing residues: Phe172, Pro223, and Val226 of DhlA; Trp207, Pro208, and Ile211 of LinB; and Phe205, Pro206, and Ile209 of DhaA. The differences in the halide stabilizing residues of three haloalkane dehalogenases are discussed in the light of molecular adaptation of these enzymes to their substrates.