Sequence fingerprint and structural analysis of the SCOR enzyme A3DFK9 from Clostridium thermocellum

We have identified a highly conserved fingerprint of 40 residues in the TGYK subfamily of the short‐chain oxidoreductase enzymes. The TGYK subfamily is defined by the presence of an N‐terminal TGxxxGxG motif and a catalytic YxxxK motif. This subfamily contains more than 12,000 members, with individual members displaying unique substrate specificities. The 40 fingerprint residues are critical to catalysis, cofactor binding, protein folding, and oligomerization but are substrate independent. Their conservation provides critical insight into evolution of the folding and function of TGYK enzymes. Substrate specificity is determined by distinct combinations of residues in three flexible loops that make up the substrate‐binding pocket. Here, we report the structure determinations of the TGYK enzyme A3DFK9 from Clostridium thermocellum in its apo form and with bound NAD⁺ cofactor. The function of this protein is unknown, but our analysis of the substrate‐binding loops putatively identifies A3DFK9 as a carbohydrate or polyalcohol metabolizing enzyme. C. thermocellum has potential commercial applications because of its ability to convert biomaterial into ethanol. A3DFK9 contains 31 of the 40 TGYK subfamily fingerprint residues. The most significant variations are the substitution of a cysteine (Cys84) for a highly conserved glycine within a characteristic VNNAG motif, and the substitution of a glycine (Gly106) for a highly conserved asparagine residue at a helical kink. Both of these variations occur at positions typically participating in the formation of a catalytically important proton transfer network. An alternate means of stabilizing this proton wire was observed in the A3DFK9 crystal structures. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Specificity and phenetic relationships of iron- and manganese-containing superoxide dismutases on the basis of structure and sequence comparisons

The iron- and manganese-containing superoxide dismutases (Fe/Mn-SOD) share the same chemical function and spatial structure but can be distinguished according to their modes of oligomerization and their metal ion specificity. They appear as homodimers or homotetramers and usually require a specific metal for activity. On the basis of 261 aligned SOD sequences and 12 superimposed x-ray structures, two phenetic trees were constructed, one sequence-based and the other structure-based. Their comparison reveals the imperfect correlation of sequence and structural changes; hyperthermophilicity requires the largest sequence alterations, whereas dimer/tetramer and manganese/iron specificities are induced by the most sizable structural differences within the monomers. A systematic investigation of sequence and structure characteristics conserved in all aligned SOD sequences or in subsets sharing common oligomeric and/or metal specificities was performed. Several residues were identified as guaranteeing the common function and dimeric conformation, others as determining the tetramer formation, and yet others as potentially responsible for metal specificity. Some form cation-pi interactions between an aromatic ring and a fully or partially positively charged group, suggesting that these interactions play a significant role in the structure and function of SOD enzymes. Dimer/tetramer- and iron/manganese-specific fingerprints were derived from the set of conserved residues; they can be used to propose selected residue substitutions in view of the experimental validation of our in silico derived hypotheses.     

Prion protein region 23–32 interacts with tubulin and inhibits microtubule assembly

In previous studies we have demonstrated that prion protein (PrP) binds directly to tubulin and this interaction leads to the inhibition of microtubule formation by inducement of tubulin oligomerization. This report is aimed at mapping the regions of PrP and tubulin involved in the interaction and identification of PrP domains responsible for tubulin oligomerization. Preliminary studies focused our attention to the N‐terminal flexible part of PrP encompassing residues 23–110. Using a panel of deletion mutants of PrP, we identified two microtubule‐binding motifs at both ends of this part of the molecule. We found that residues 23–32 constitute a major site of interaction, whereas residues 101–110 represent a weak binding site. The crucial role of the 23–32 sequence in the interaction with tubulin was confirmed employing chymotryptic fragments of PrP. Surprisingly, the octarepeat region linking the above motifs plays only a supporting role in the interaction. The binding of Cu²⁺ to PrP did not affect the interaction. We also demonstrate that PrP deletion mutants lacking residues 23–32 exhibit very low efficiency in the inducement of tubulin oligomerization. Moreover, a synthetic peptide corresponding to this sequence, but not that identical with fragment 101–110, mimics the effects of the full‐length protein on tubulin oligomerization and microtubule assembly. At the cellular level, peptide composed of the PrP motive 23–30 and signal sequence (1–22) disrupted the microtubular cytoskeleton. Using tryptic and chymotryptic fragments of α‐ and β‐tubulin, we mapped the docking sites for PrP within the C‐terminal domains constituting the outer surface of microtubule. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Orientation and oligomerization specificity of the Bcr coiled-coil oligomerization domain

The Bcr oligomerization domain, from the Bcr-Abl oncoprotein, is an attractive therapeutic target for treating leukemias because it is required for cellular transformation. The domain homodimerizes via an antiparallel coiled coil with an adjacent short, helical swap domain. Inspection of the coiled-coil sequence does not reveal obvious determinants of helix-orientation specificity, raising the possibility that the antiparallel orientation preference and/or the dimeric oligomerization state are due to interactions of the swap domains. To better understand how structural specificity is encoded in Bcr, coiled-coil constructs containing either an N- or C-terminal cysteine were synthesized without the swap domain. When cross-linked to adopt exclusively parallel or antiparallel orientations, these showed similar circular dichroism spectra. Both constructs formed coiled-coil dimers, but the antiparallel construct was approximately 16 degrees C more stable than the parallel to thermal denaturation. Equilibrium disulfide-exchange studies confirmed that the isolated coiled-coil homodimer shows a very strong preference for the antiparallel orientation. We conclude that the orientation and oligomerization preferences of Bcr are not caused by the presence of the swap domains, but rather are directly encoded in the coiled-coil sequence. We further explored possible determinants of structural specificity by mutating residues in the d position of the coiled-coil core. Some of the mutations caused a change in orientation specificity, and all of the mutations led to the formation of higher-order oligomers. In the absence of the swap domain, these residues play an important role in disfavoring alternate states and are especially important for encoding dimeric oligomerization specificity.     

Oligomerization properties of GCN4 leucine zipper e and g position mutants.

Putative intersubunit electrostatic interactions between charged amino acids on the surfaces of the dimer interfaces of leucine zippers (g-e'' ion pairs) have been implicated as determinants of dimerization specificity. To evaluate the importance of these ionic interactions in determining the specificity of dimer formation, we constructed a pool of > 65,000 GCN4 leucine zipper mutants in which all the e and g positions are occupied by different combinations of alanine, glutamic acid, lysine, or threonine. The oligomerization properties of these mutants were evaluated based on the phenotypes of cells expressing lambda repressor-leucine zipper fusion proteins. About 90% of the mutants do not form stable homooligomers. Surprisingly, approximately 8% of the mutant sequences have phenotypes consistent with the formation of higher-order (> dimer) oligomers, which can be classified into three types based on sequence features. The oligomerization states of mutants from two of these types were determined by characterizing purified fusion proteins. The Type I mutant behaved as a tetramer under all tested conditions, whereas the Type III mutant formed a variety of higher-order oligomers, depending on the solution conditions. Stable homodimers comprise less than 3% of the pool; several g-e'' positions in these mutants could form attractive ion pairs. Putative repulsive ion pairs are not found among the homodimeric mutants. However, patterns of charged residues at the e and g positions do not seem to be sufficient to predict either homodimer or heterodimer formation among the mutants.     

Role of sequence evolution and conformational dynamics in the substrate specificity and oligomerization mode of thymidylate kinases

Thymidylate kinase (TMK) is a key enzyme for the synthesis of DNA, making it an important target for the development of anticancer, antibacterial, and antiparasitic drugs. TMK homologs exhibit significant variations in sequence, residue conformation, substrate specificity, and oligomerization mode. However, the influence of sequence evolution and conformational dynamics on its quaternary structure and function has not been studied before. Based on extensive sequence and structure analyses, our study detected several non-conserved residues which are linked by co-evolution and are implicated in the observed variations in flexibility, oligomeric assembly, and substrate specificity among the homologs. These lead to differences in the pattern of interactions at the active site in TMKs of different specificity. The method was further tested on TMK from Sulfolobus tokodaii (StTMK) which has substantial differences in sequence and structure compared to other TMKs. Our analyses pointed to a more flexible dTMP-binding site in StTMK compared to the other homologs. Binding assays proved that the protein can accommodate both purine and pyrimidine nucleotides at the dTMP binding site with comparable affinity. Additionally, the residues responsible for the narrow specificity of Brugia malayi TMK, whose three-dimensional structure is unavailable, were detected. Our study provides a residue-level understanding of the differences observed among TMK homologs in previous experiments. It also illustrates the correlation among sequence evolution, conformational dynamics, oligomerization mode, and substrate recognition in TMKs and detects co-evolving residues that affect binding, which should be taken into account while designing novel inhibitors.     

Relaxing the nicotinamide cofactor specificity of phosphite dehydrogenase by rational design

Homology modeling was used to identify two particular residues, Glu175 and Ala176, in Pseudomonas stutzeri phosphite dehydrogenase (PTDH) as the principal determinants of nicotinamide cofactor (NAD(+) and NADP(+)) specificity. Replacement of these two residues by site-directed mutagenesis with Ala175 and Arg176 both separately and in combination resulted in PTDH mutants with relaxed cofactor specificity. All three mutants exhibited significantly better catalytic efficiency for both cofactors, with the best kinetic parameters displayed by the double mutant, which had a 3.6-fold higher catalytic efficiency for NAD(+) and a 1000-fold higher efficiency for NADP(+). The cofactor specificity was changed from 100-fold in favor of NAD(+) for the wild-type enzyme to 3-fold in favor of NADP(+) for the double mutant. Isoelectric focusing of the proteins in a nondenaturing gel showed that the replacement with more basic residues indeed changed the effective pI of the protein. HPLC analysis of the enzymatic products of the double mutant verified that the reaction proceeded to completion using either substrate and produced only the corresponding reduced cofactor and phosphate. Thermal inactivation studies showed that the double mutant was protected from thermal inactivation by both cofactors, while the wild-type enzyme was protected by only NAD(+). The combined results provide clear evidence that Glu175 and Ala176 are both critical for nicotinamide cofactor specificity. The rationally designed double mutant might be useful for the development of an efficient in vitro NAD(P)H regeneration system for reductive biocatalysis.     

Roles of low specificity and cofactor interaction sites on thrombin during factor XIII activation. Competition for cofactor sites on thrombin determines its fate

Factor XIII is activated by thrombin, and this reaction is enhanced by the presence of fibrin(ogen). Using a substrate-based screening assay for factor XIII activity complemented by kinetic analysis of activation peptide cleavage, we show by using thrombin mutants of surface-exposed residues that Arg-178, Arg-180, Asp-183, Glu-229, Arg-233, and Trp-50 of thrombin are necessary for direct activation of factor XIII. These residues define a low specificity site known to be important also for both protein C activation and for inhibition of thrombin by antithrombin. The enhancing effect of fibrinogen occurs as a consequence of its conversion to fibrin and subsequent polymerization. Surface residues of thrombin further involved in high specificity fibrin-enhanced factor XIII activation were identified as His-66, Tyr-71, and Asn-74. These residues represent a distinct interaction site on thrombin (within exosite I) also employed by thrombomodulin in its cofactor-enhanced activation of protein C. In competition experiments, thrombomodulin inhibited fibrin-enhanced factor XIII activation. Based upon these and prior published results, we propose that the polymerization process forms a fibrin cofactor that acts to approximate thrombin and factor XIII bound to separate and complementary domains of fibrinogen. This enables enhanced factor XIII activation to be localized around the fibrin clot. We also conclude that proximity to and competition for cofactor interaction sites primarily directs the fate of thrombin.     

Roles of individual domains and conserved motifs of the AAA+ chaperone ClpB in oligomerization,ATP hydrolysis,and chaperone activity

ClpB of Escherichia coli is an ATP-dependent ring-forming chaperone that mediates the resolubilization of aggregated proteins in cooperation with the DnaK chaperone system. ClpB belongs to the Hsp100/Clp subfamily of AAA+ proteins and is composed of an N-terminal domain and two AAA-domains that are separated by a "linker" region. Here we present a detailed structure-function analysis of ClpB, dissecting the individual roles of ClpB domains and conserved motifs in oligomerization, ATP hydrolysis, and chaperone activity. Our results show that ClpB oligomerization is strictly dependent on the presence of the C-terminal domain of the second AAA-domain, while ATP binding to the first AAA-domains stabilized the ClpB oligomer. Analysis of mutants of conserved residues in Walker A and B and sensor 2 motifs revealed that both AAA-domains contribute to the basal ATPase activity of ClpB and communicate in a complex manner. Chaperone activity strictly depends on ClpB oligomerization and the presence of a residual ATPase activity. The N-domain is dispensable for oligomerization and for the disaggregating activity in vitro and in vivo. In contrast the presence of the linker region, although not involved in oligomerization, is essential for ClpB chaperone activity.     

The three-dimensional solution structure of Aesculus hippocastanum antimicrobial protein 1 determined by 1H nuclear magnetic resonance

Aesculus hippocastanum antimicrobial protein 1 (Ah-AMP1) is a plant defensin isolated from horse chestnuts. The plant defensins have been divided in several subfamilies according to their amino acid sequence homology. Ah-AMP1, belonging to subfamily A2, inhibits growth of a broad range of fungi. So far, a three-dimensional structure has been determined only for members of subfamilies A3 and B2. In order to understand activity and specificity of these plant defensins, the structure of a protein belonging to subfamily A2 is needed. We report the three-dimensional solution structure of Ah-AMP1 as determined from two-dimensional 1H nuclear magnetic resonance data. The structure features all the characteristics of the "cysteine-stabilized alpha beta-motif." A comparison of the structure, the electrostatic potential surface and regions important for interaction with the fungal receptor, is made with Rs-AFP1 (plant defensin of subfamily A3). Thus, residues important for activity and specificity have been assigned.     

The insect immune protein scolexin is a novel serine proteinase homolog

Scolexin is a coagulation-provoking plasma protein induced in response to bacterial or viral infection of larval Manduca sexta, a large lepidopterous insect. Here we report the isolation and sequencing of two cDNA clones that code for scolexin isoforms sharing 80% sequence identity. The scolexin sequences have low but recognizable sequence similarity to members of the chymotrypsin family and represent a new subfamily of chymotrypsin-like serine proteinases. Comparison with known structures reveals the conservation of key catalytic residues and a possible specificity for small nonpolar residues. Most remarkable is the absence of a canonical activation peptide cleavage site. This suggests that the regulation of scolexin activity will involve a novel activation mechanism.     

Acidic residues at the carboxyl terminus of p60c-src are required for regulation of tyrosine kinase activity and transformation.

Protein phosphorylation sites act to transduce signals into changes in enzymatic activity, representing a point of interaction within a regulatory pathway. The amino acid sequence surrounding a phosphorylation site may well have several functions, including recognition by an appropriate kinase. By generating random mutations in its immediate vicinity, we have examined the sequence requirements of a regulatory tyrosine phosphorylation site, Tyr527, in the proto-oncogene product, p60c-src. The transforming and kinase activities of p60c-src are repressed by phosphorylation of Tyr527. Mutations were made around Tyr527 without changing Tyr527 or the kinase domain. Twenty-nine mutants were sequenced and classified as transforming or nontransforming for Rat-2 cells. Nontransforming mutants contained a surprising variety of COOH-terminal mutations, although acidic residues were present at positions 518 and 524 in all nontransforming mutants. Transforming mutants that contained single-residue changes at Asp518 and Ser522 demonstrated the importance of these residues. Other transforming mutants contained two or more substitutions, but the results are most simply explained if residues Glu524 and Thr523 are also important for normal regulation. Transforming mutations reduced the phosphorylation of Tyr527. We conclude that only a few of the residues in the COOH terminus other than Tyr527 are required to ensure normal phosphorylation and repression of activity in fibroblasts. Other residues may have been conserved during evolution to permit normal function and regulation in other cell types.     

The conserved sequence NXX[S/T]HX[S/T]QDXXXT of the lactate/pyruvate:H+ symporter subfamily defines the function of the substrate translocation pathway

In Saccharomyces cerevisiae Jen1p is a lactate/proton symporter belonging to the lactate/pyruvate:H⁺ symporter subfamily (TC#2.A.1.12.2) of the Major Facilitator Superfamily. We investigated structure-function relationships of Jen1p using a rational mutational analysis based on the identification of conserved amino acid residues. In particular, we studied the conserved sequence ³⁷⁹NXX[S/T]HX[S/T]QDXXXT³⁹¹. Substitution of amino acid residues N379, H383 or D387, even with very similar amino acids, resulted in a dramatic reduction of lactate and pyruvate uptake, but conserved measurable acetate transport. Acetate transport inhibition assays showed that these mutants conserve the ability to bind, but do not transport, lactate and pyruvate. More interestingly, the double mutation H383D/D387H, while behaving as a total loss-of-function allele for lactate and pyruvate uptake, can fully restore the kinetic parameters of Jen1p for acetate transport. Thus, residues N379, H383 or D387 affect both the transport capacity and the specificity of Jen1p. Substitutions of Q386 and T391 resulted in no or moderate changes in Jen1p transport capacities for lactate, pyruvate and acetate. On the other hand, Q386N reduces the binding affinities for all Jen1p substrates, while Q386A increases the affinity specifically for pyruvate. We also tested Jen1p specificity for a range of monocarboxylates. Several of the mutants studied showed altered inhibition constants for these acids. These results and 3D in silico modelling by homology threading suggest that the conserved motif analyzed is part of the substrate translocation pathway in the lactate/pyruvate:H⁺ symporter subfamily.     

A detailed thermodynamic analysis of ras/effector complex interfaces

Many cellular functions are based on the interaction and crosstalk of various signaling proteins. Among these, members of the Ras family of small GTP-binding proteins are important for communicating signals into different pathways. In order to answer the question of how binding affinity and specificity is achieved, we analyzed binding energetics on the molecular level, with reference to the available structural data. The interaction of two members of the Ras subfamily with two different effector proteins, namely Raf and RalGDS, were investigated using isothermal titration calorimetry and a fluorescence-based method. Experiments with alanine mutants, located in the complex interfaces, yielded an energy map for the contact areas of the Ras/effector complexes, which could be differentiated into enthalpy and entropy contributions. In addition, by using double mutant cycle analysis, we probed the energetic contribution of selected pairs of amino acid residues. The resulting energy landscapes of the Ras/effector interface areas show a highly different topology when comparing the two effectors, Raf and RalGDS, demonstrating the specificity of the respective interactions. Particularly, we observe a high degree of compensating effects between enthalpy and entropy; differences between these components are much greater than the overall free energy differences. This is observed also when using the software FOLD-X to predict the effect of point mutations on the crystal structures of the different complexes. Prediction of the free energy changes shows a very good correlation with the experimentally observed energies. Furthermore, in line with experimental data, energy decomposition indicates that many different components of large magnitude counteract each other to produce a smaller change in overall free energy, illustrating the importance of long-range electrostatic forces in complex formation.     

Inhibition of the thioredoxin-dependent activation of the NADP-malate dehydrogenase and cofactor specificity

The chloroplastic NADP-malate dehydrogenase is activated by reduction of its N- and C-terminal disulfides by reduced thioredoxin. The activation is inhibited by NADP(+), the oxidized form of the cofactor. Previous studies suggested that the C-terminal disulfide was involved in this process. Recent structural data pointed toward a possible direct interaction between the C terminus of the oxidized enzyme and the cofactor. In the present study, the relationship between the cofactor specificity for catalysis and for inhibition of activation has been investigated by changing the cofactor specificity of the enzyme by substitution of selected residues of the cofactor-binding site. An NAD-specific thiol-regulated MDH was engineered. Its activation was inhibited by NAD(+) but no longer by NADP(+). These results demonstrate that the oxidized cofactor is bound at the same site as the reduced cofactor and support the idea of a direct interaction between the negatively charged C-terminal end of the enzyme and the positively charged nicotinamide ring of the cofactor, in agreement with the structural data. The structural requirements for cofactor specificity are modeled and discussed.     

Molecular determinants of the cofactor specificity of ribitol dehydrogenase, a short-chain dehydrogenase/reductase

Ribitol dehydrogenase from Zymomonas mobilis (ZmRDH) catalyzes the conversion of ribitol to d-ribulose and concomitantly reduces NAD(P)(+) to NAD(P)H. A systematic approach involving an initial sequence alignment-based residue screening, followed by a homology model-based screening and site-directed mutagenesis of the screened residues, was used to study the molecular determinants of the cofactor specificity of ZmRDH. A homologous conserved amino acid, Ser156, in the substrate-binding pocket of the wild-type ZmRDH was identified as an important residue affecting the cofactor specificity of ZmRDH. Further insights into the function of the Ser156 residue were obtained by substituting it with other hydrophobic nonpolar or polar amino acids. Substituting Ser156 with the negatively charged amino acids (Asp and Glu) altered the cofactor specificity of ZmRDH toward NAD(+) (S156D, [k(cat)/K(m)(,NAD)]/[k(cat)/K(m)(,NADP)] = 10.9, where K(m)(,NAD) is the K(m) for NAD(+) and K(m)(,NADP) is the K(m) for NADP(+)). In contrast, the mutants containing positively charged amino acids (His, Lys, or Arg) at position 156 showed a higher efficiency with NADP(+) as the cofactor (S156H, [k(cat)/K(m)(,NAD)]/[k(cat)/K(m)(,NADP)] = 0.11). These data, in addition to those of molecular dynamics and isothermal titration calorimetry studies, suggest that the cofactor specificity of ZmRDH can be modulated by manipulating the amino acid residue at position 156.     

Insights into the role of d-amino acid oxidase mutations in amyotrophic lateral sclerosis

Missense mutations in the coding region of d -amino acid oxidase (DAO) have been found in patients suffering from amyotrophic lateral sclerosis (ALS). Mutations primarily impair the enzymatic activity of DAO and cause neurodegeneration due to an abnormal accumulation of d -serine in the spinal cord. However, the structural and dynamic changes that lead to impaired enzymatic activity are not fully understood. We present here extensive molecular dynamics simulations of wild-type, and all reported ALS-associated DAO mutants to elucidate the plausible mechanisms of impaired enzymatic activity, a critical function needed for neuroprotection. Simulation results show that DAO mutations disrupt several key interactions with the active site residues and decrease the conformational flexibility of active site loop comprising 216 to 228 residues, necessary for substrate binding and product release. This conformational restriction of the active site loop in the mutants is mainly due to the distortion of critical salt bridge and hydrogen bond interactions compared with wild-type. Furthermore, binding free energy calculations show that DAO mutants have a lower binding affinity toward cofactor flavin adenine dinucleotide and substrate imino-serine than the wild-type. A closer look at the cofactor and substrate interaction profiles further show that DAO mutants have lost several critical interactions with the neighboring residues as compared with wild-type. Taken together, this study provides first-hand explanation of crucial structural features that lead to the loss of enzymatic function in DAO mutants and highlights the need of further genomic scans of patients with ALS to map the association of novel DAO variants in ALS pathophysiology.     

Contribution of basic residues of the autolysis loop to the substrate and inhibitor specificity of factor IXa

The autolysis loop (residues 143-154 in chymotrypsinogen numbering) plays a pivotal role in determining the macromolecular substrate and inhibitor specificity of coagulation proteases. This loop in factor IXa (FIXa) has 3 basic residues (Arg143, Lys147, and Arg150) whose contribution to the protease specificity of factor IXa has not been studied. Here, we substituted these residues individually with Ala in Gla-domainless forms of recombinant factor IX expressed in mammalian cells. All mutants exhibited normal amidolytic activities toward a FIXa-specific chromogenic substrate. However, Arg143 and Lys147 mutants showed a approximately 3- to 6-fold impairment in FX activation, whereas the Arg150 mutant activated factor X normally both in the absence and presence of factor VIIIa. By contrast, Arg143 and Lys147 mutants reacted normally with antithrombin (AT) in both the absence and presence of the cofactor, heparin. However, the reactivity of the Arg150 mutant with AT was impaired 6.6-fold in the absence of heparin and 33- to 70-fold in the presence of pentasaccharide and full-length heparins. These results suggest that Arg143 and Lys147 of the autolysis loop are recognition sites for FX independent of factor VIIIa, and Arg150 is a specific recognition site for AT that can effectively interact with AT only if the serpin is in the heparin-activated conformation.