Prediction of coenzyme specificity in dehydrogenases/reductases. A hidden Markov model-based method and its application on complete genomes

Dehydrogenases and reductases are enzymes of fundamental metabolic importance that often adopt a specific structure known as the Rossmann fold. This fold, consisting of a six-stranded beta-sheet surrounded by alpha-helices, is responsible for coenzyme binding. We have developed a method to identify Rossmann folds and predict their coenzyme specificity (NAD, NADP or FAD) using only the amino acid sequence as input. The method is based upon hidden Markov models and sequence pattern analysis. The prediction sensitivity is 79% and the selectivity close to 100%. The method was applied on a set of 68 genomes, representing the three kingdoms archaea, bacteria and eukaryota. In prokaryotes, 3% of the genes were found to code for Rossmann-fold proteins, while the corresponding ratio in eukaryotes is only around 1%. In all genomes, NAD is the most preferred cofactor (41-49%), followed by NADP with 30-38%, while FAD is the least preferred cofactor (21%). However, the NAD preponderance over NADP is most pronounced in archaea, and least in eukaryotes. In all three kingdoms, only 3-8% of the Rossmann proteins are predicted to have more than one membrane-spanning segment, which is much lower than the frequency of membrane proteins in general. Analysis of the major protein types in eukaryotes reveals that the most common type (26%) of the Rossmann proteins are short-chain dehydrogenases/reductases. In addition, the identified Rossmann proteins were analyzed with respect to further protein types, enzyme classes and redundancy. The described method is available at http://www.ifm.liu.se/bioinfo, where the preferred coenzyme and its binding region are predicted given an amino acid sequence as input.     

An Ancient Fingerprint Indicates the Common Ancestry of Rossmann-Fold Enzymes Utilizing Different Ribose-Based Cofactors

Nucleoside-based cofactors are presumed to have preceded proteins. The Rossmann fold is one of the most ancient and functionally diverse protein folds, and most Rossmann enzymes utilize nucleoside-based cofactors. We analyzed an omnipresent Rossmann ribose-binding interaction: a carboxylate side chain at the tip of the second β-strand (β2-Asp/Glu). We identified a canonical motif, defined by the β2-topology and unique geometry. The latter relates to the interaction being bidentate (both ribose hydroxyls interacting with the carboxylate oxygens), to the angle between the carboxylate and the ribose, and to the ribose’s ring configuration. We found that this canonical motif exhibits hallmarks of divergence rather than convergence. It is uniquely found in Rossmann enzymes that use different cofactors, primarily SAM (S-adenosyl methionine), NAD (nicotinamide adenine dinucleotide), and FAD (flavin adenine dinucleotide). Ribose-carboxylate bidentate interactions in other folds are not only rare but also have a different topology and geometry. We further show that the canonical geometry is not dictated by a physical constraint—geometries found in noncanonical interactions have similar calculated bond energies. Overall, these data indicate the divergence of several major Rossmann-fold enzyme classes, with different cofactors and catalytic chemistries, from a common pre-LUCA (last universal common ancestor) ancestor that possessed the β2-Asp/Glu motif.     

Flavogenomics--a genomic and structural view of flavin-dependent proteins

Riboflavin (vitamin B(2)) serves as the precursor for FMN and FAD in almost all organisms that utilize the redox-active isoalloxazine ring system as a coenzyme in enzymatic reactions. The role of flavin, however, is not limited to redox processes, as ～ 10% of flavin-dependent enzymes catalyze nonredox reactions. Moreover, the flavin cofactor is also widely used as a signaling and sensing molecule in biological processes such as phototropism and nitrogen fixation. Here, we present a study of 374 flavin-dependent proteins analyzed with regard to their function, structure and distribution among 22 archaeal, eubacterial, protozoan and eukaryotic genomes. More than 90% of flavin-dependent enzymes are oxidoreductases, and the remaining enzymes are classified as transferases (4.3%), lyases (2.9%), isomerases (1.4%) and ligases (0.4%). The majority of enzymes utilize FAD (75%) rather than FMN (25%), and bind the cofactor noncovalently (90%). High-resolution structures are available for about half of the flavoproteins. FAD-containing proteins predominantly bind the cofactor in a Rossmann fold (～ 50%), whereas FMN-containing proteins preferably adopt a (βα)(8)-(TIM)-barrel-like or flavodoxin-like fold. The number of genes encoding flavin-dependent proteins varies greatly in the genomes analyzed, and covers a range from ～ 0.1% to 3.5% of the predicted genes. It appears that some species depend heavily on flavin-dependent oxidoreductases for degradation or biosynthesis, whereas others have minimized their flavoprotein arsenal. An understanding of 'flavin-intensive' lifestyles, such as in the human pathogen Mycobacterium tuberculosis, may result in valuable new intervention strategies that target either riboflavin biosynthesis or uptake.     

A structurally conserved water molecule in Rossmann dinucleotide-binding domains

A computational comparison of 102 high-resolution (相似文献   

A shared binding site for NAD+ and coenzyme A in an acetaldehyde dehydrogenase involved in bacterial degradation of aromatic compounds

The meta-cleavage pathway for catechol is a central pathway for the bacterial dissimilation of a wide variety of aromatic compounds, including phenols, methylphenols, naphthalenes, and biphenyls. The last enzyme of the pathway is a bifunctional aldolase/dehydrogenase that converts 4-hydroxy-2-ketovalerate to pyruvate and acetyl-CoA via acetaldehyde. The structure of the NAD (+)/CoASH-dependent aldehyde dehydrogenase subunit is similar to that of glyceraldehyde-3-phosphate dehydrogenase, with a Rossmann fold-based NAD (+) binding site observed in the NAD (+)-enzyme complex [Manjasetty, B. A., et al. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 6992-6997]. However, the location of the CoASH binding site was not determined. In this study, hydrogen-deuterium exchange experiments, coupled with peptic digest and mass spectrometry, were used to examine cofactor binding. The pattern of hydrogen-deuterium exchange in the presence of CoASH was almost identical to that observed with NAD (+), consistent with the two cofactors sharing a binding site. This is further supported by the observations that either CoASH or NAD (+) is able to elute the enzyme from an NAD (+) affinity column and that preincubation of the enzyme with NAD (+) protects against inactivation by CoASH. Consistent with these data, models of the CoASH complex generated using AUTODOCK showed that the docked conformation of CoASH can fully occupy the cavity containing the enzyme active site, superimposing with the NAD (+) cofactor observed in the X-ray crystal structure. Although CoASH binding Rossmann folds have been described previously, this is the first reported example of a Rossmann fold that can alternately bind CoASH or NAD (+) cofactors required for enzymatic catalysis.     

Knock down of cytosolic phospholipase A2: an antisense oligonucleotide having a nuclear localization binds a C-terminal motif of glyceraldehyde-3-phosphate dehydrogenase

We have previously shown that an antisense, effective in the knock down of cytosolic phospholipase A2 (cPLA2), localizes mainly in the nucleus of human endothelial cells and monocytes and that glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is involved in its nuclear localization. In this study, we clarify how GAPDH participates in the nuclear localization of this antisense oligodeoxynucleotide (ODN) directed against cPLA2 mRNA. A central TAAAT motif providing specificity and high affinity binding was assumed to interact with the enzyme Rossmann fold region on the basis of competition to this site by NAD+. To asses whether the TAAAT motif interacts directly with the enzyme Rossmann fold region, we evaluated the binding to GAPDH of different oligonucleotides and the effect of competitors such as NAD+, NADH, mononucleotides, DNA, polyribonucleic acids and polyanions. We found that the dissociation constant for TAAAT containing oligonucleotides was three--to fivefold higher with respect to oligo not containing this motif. By covalently linking 32P-labeled cPLA2p(N)16 to GAPDH and after executing hydrolysis with hydroxylamine, the labeling was exclusively found in the C-terminal domain (aa 286-334). These results indicate that the antisense oligonucleotide interacts with a site not having a defined function but which can be negatively allosterically regulated when NAD+ or polynucleotides are bound to Rossmann fold.     

The yeast F(1)F(0)-ATP synthase: analysis of the molecular organization of subunit g and the importance of a conserved GXXXG motif

The F(1)F(0)-ATP synthase enzyme is located in the inner mitochondrial membrane, where it forms dimeric complexes. Dimerization of the ATP synthase involves the physical association of the neighboring membrane-embedded F(0)-sectors. In yeast, the F(0)-sector subunits g and e (Su g and Su e, respectively) play a key role in supporting the formation of ATP synthase dimers. In this study we have focused on Su g to gain a better understanding of the function and the molecular organization of this subunit within the ATP synthase complex. Su g proteins contain a GXXXG motif (G is glycine, and X is any amino acid) in their single transmembrane segment. GXXXG can be a dimerization motif that supports helix-helix interactions between neighboring transmembrane segments. We demonstrate here that the GXXXG motif is important for the function and in particular for the stability of Su g within the ATP synthase. Using site-directed mutagenesis and cross-linking approaches, we demonstrate that Su g and Su e interact, and our findings emphasize the importance of the membrane anchor regions of these proteins for their interaction. Su e also contains a conserved GXXXG motif in its membrane anchor. However, data presented here would suggest that an intact GXXXG motif in Su g is not essential for the Su g-Su e interaction. We suggest that the GXXXG motif may not be the sole basis for a Su g-Su e interaction, and possibly these dimerization motifs may enable both Su g and Su e to interact with another mitochondrial protein.     

Sequence-structure analysis of FAD-containing proteins

We have analyzed structure-sequence relationships in 32 families of flavin adenine dinucleotide (FAD)-binding proteins, to prepare for genomic-scale analyses of this family. Four different FAD-family folds were identified, each containing at least two or more protein families. Three of these families, exemplified by glutathione reductase (GR), ferredoxin reductase (FR), and p-cresol methylhydroxylase (PCMH) were previously defined, and a family represented by pyruvate oxidase (PO) is newly defined. For each of the families, several conserved sequence motifs have been characterized. Several newly recognized sequence motifs are reported here for the PO, GR, and PCMH families. Each FAD fold can be uniquely identified by the presence of distinctive conserved sequence motifs. We also analyzed cofactor properties, some of which are conserved within a family fold while others display variability. Among the conserved properties is cofactor directionality: in some FAD-structural families, the adenine ring of the FAD points toward the FAD-binding domain, whereas in others the isoalloxazine ring points toward this domain. In contrast, the FAD conformation and orientation are conserved in some families while in others it displays some variability. Nevertheless, there are clear correlations among the FAD-family fold, the shape of the pocket, and the FAD conformation. Our general findings are as follows: (a) no single protein 'pharmacophore' exists for binding FAD; (b) in every FAD-binding family, the pyrophosphate moiety binds to the most strongly conserved sequence motif, suggesting that pyrophosphate binding is a significant component of molecular recognition; and (c) sequence motifs can identify proteins that bind phosphate-containing ligands.     

Identification of residues important for NAD+ binding by the Thermotoga maritima alpha-glucosidase AglA,a member of glycoside hydrolase family 4

The NAD+-requiring enzymes of glycoside hydrolase family 4 (GHF4) contain a region with a conserved Gly-XXX-Gly-Ser (GXGS) motif near their N-termini that is reminiscent of the fingerprint region of the Rossmann fold, a conserved structural motif of classical nicotinamide nucleotide-binding proteins. The function of this putative NAD+-binding motif in the alpha-glucosidase AglA of Thermotoga maritima was probed by directed mutagenesis. The K(d) for NAD+ of the AglA mutants G10A, G12A and S13A was increased by about 300-, 5-, and 9-fold, respectively, while their K(m) for p-nitrophenyl-alpha-glucopyranoside was not seriously affected. The results indicate that the GXGS motif is indeed important for NAD+ binding by the glycosidases of GHF4.     

Cytochrome b5 reductase: role of the si-face residues, proline 92 and tyrosine 93, in structure and catalysis

The conserved sequence motif "RxY(T)(S)xx(S)(N)" coordinates flavin binding in NADH:cytochrome b(5) reductase (cb(5)r) and other members of the flavin transhydrogenase superfamily of oxidoreductases. To investigate the roles of Y93, the third and only aromatic residue of the "RxY(T)(S)xx(S)(N)" motif, that stacks against the si-face of the flavin isoalloxazine ring, and P92, the second residue in the motif that is also in close proximity to the FAD moiety, a series of rat cb(5)r variants were produced with substitutions at either P92 or Y93, respectively. The proline mutants P92A, G, and S together with the tyrosine mutants Y93A, D, F, H, S, and W were recombinantly expressed in E. coli and purified to homogeneity. Each mutant protein was found to bind FAD in a 1:1 cofactor:protein stoichiometry while UV CD spectra suggested similar secondary structure organization among all nine variants. The tyrosine variants Y93A, D, F, H, and S exhibited varying degrees of blue-shift in the flavin visible absorption maxima while visible CD spectra of the Y93A, D, H, S, and W mutants exhibited similar blue-shifted maxima together with changes in absorption intensity. Intrinsic flavin fluorescence was quenched in the wild type, P92S and A, and Y93H and W mutants while Y93A, D, F, and S mutants exhibited increased fluorescence when compared to free FAD. The tyrosine variants Y93A, D, F, and S also exhibited greater thermolability of FAD binding. The specificity constant (k(cat)/K(m)(NADH)) for NADH:FR activity decreased in the order wild type > P92S > P92A > P92G > Y93F > Y93S > Y93A > Y93D > Y93H > Y93W with the Y93W variant retaining only 0.5% of wild-type efficiency. Both K(s)(H4NAD) and K(s)(NAD+) values suggested that Y93A, F, and W mutants had compromised NADH and NAD(+) binding. Thermodynamic measurements of the midpoint potential (E degrees ', n = 2) of the FAD/FADH(2) redox couple revealed that the potentials of the Y93A and S variants were approximately 30 mV more positive than that of wild-type cb(5)r (E degrees ' = -268 mV) while that of Y93H was approximately 30 mV more negative. These results indicate that neither P92 nor Y93 are critical for flavin incorporation in cb(5)r and that an aromatic side chain is not essential at position 93, but they demonstrate that Y93 forms contacts with the FAD that effectively modulate the spectroscopic, catalytic, and thermodynamic properties of the bound cofactor.     

A path from primary protein sequence to ligand recognition

A novel method to organize protein structural information based solely on sequence is presented. The method clusters proteins into families that correlate with the three-dimensional protein structure and the conformation of the bound ligands. This procedure was applied to nicotinamide adenine dinucleotide [NAD(P)]-utilizing enzymes to identify a total of 94 sequence families, 53 of which are structurally characterized. Each of the structurally characterized proteins within a sequence family correlates to a single protein fold and to a common bound conformation of NAD(P). A wide range of structural folds is identified that recognize NAD(P), including Rossmann folds and beta/alpha barrels. The defined sequence families can be used to identify the type and prevalence of NAD(P)-utilizing enzymes in the proteomes of sequenced organisms. The proteome of Mycobacterium tuberculosis was mined to generate a proteome-wide profile of NAD(P)-utilizing enzymes coded by this organism. This enzyme family comprises approximately 6% of the open reading frames, with the largest subgroup being the Rossmann fold, short-chain dehydrogenases. The preponderance of short-chain dehydrogenases correlates strongly with the phenotype of M. tuberculosis, which is characterized as having one of the most complex prokaryotic cell walls.     

The crystal structure of NADPH:ferredoxin reductase from Azotobacter vinelandii.

NADPH:ferredoxin reductase (AvFPR) is involved in the response to oxidative stress in Azotobacter vinelandii. The crystal structure of AvFPR has been determined at 2.0 A resolution. The polypeptide fold is homologous with six other oxidoreductases whose structures have been solved including Escherichia coli flavodoxin reductase (EcFldR) and spinach, and Anabaena ferredoxin:NADP+ reductases (FNR). AvFPR is overall most homologous to EcFldR. The structure is comprised of a N-terminal six-stranded antiparallel beta-barrel domain, which binds FAD, and a C-terminal five-stranded parallel beta-sheet domain, which binds NADPH/NADP+ and has a classical nucleotide binding fold. The two domains associate to form a deep cleft where the NADPH and FAD binding sites are juxtaposed. The structure displays sequence conserved motifs in the region surrounding the two dinucleotide binding sites, which are characteristic of the homologous enzymes. The folded over conformation of FAD in AvFPR is similar to that in EcFldR due to stacking of Phe255 on the adenine ring of FAD, but it differs from that in the FNR enzymes, which lack a homologous aromatic residue. The structure of AvFPR displays three unique features in the environment of the bound FAD. Two features may affect the rate of reduction of FAD: the absence of an aromatic residue stacked on the isoalloxazine ring in the NADPH binding site; and the interaction of a carbonyl group with N10 of the flavin. Both of these features are due to the substitution of a conserved C-terminal tyrosine residue with alanine (Ala254) in AvFPR. An additional unique feature may affect the interaction of AvFPR with its redox partner ferredoxin I (FdI). This is the extension of the C-terminus by three residues relative to EcFldR and by four residues relative to FNR. The C-terminal residue, Lys258, interacts with the AMP phosphate of FAD. Consequently, both phosphate groups are paired with a basic group due to the simultaneous interaction of the FMN phosphate with Arg51 in a conserved FAD binding motif. The fourth feature, common to homologous oxidoreductases, is a concentration of 10 basic residues on the face of the protein surrounding the active site, in addition to Arg51 and Lys258.     

GXXXG and AXXXA: common alpha-helical interaction motifs in proteins, particularly in extremophiles

The GXXXG motif is a frequently occurring sequence of residues that is known to favor helix-helix interactions in membrane proteins. Here we show that the GXXXG motif is also prevalent in soluble proteins whose structures have been determined. Some 152 proteins from a non-redundant PDB set contain at least one alpha-helix with the GXXXG motif, 41 +/- 9% more than expected if glycine residues were uniformly distributed in those alpha-helices. More than 50% of the GXXXG-containing alpha-helices participate in helix-helix interactions. In fact, 26 of those helix-helix interactions are structurally similar to the helix-helix interaction of the glycophorin A dimer, where two transmembrane helices associate to form a dimer stabilized by the GXXXG motif. As for the glycophorin A structure, we find backbone-to-backbone atomic contacts of the C alpha-H...O type in each of these 26 helix-helix interactions that display the stereochemical hallmarks of hydrogen bond formation. These glycophorin A-like helix-helix interactions are enriched in the general set of helix-helix interactions containing the GXXXG motif, suggesting that the inferred C alpha-H...O hydrogen bonds stabilize the helix-helix interactions. In addition to the GXXXG motif, some 808 proteins from the non-redundant PDB set contain at least one alpha-helix with the AXXXA motif (30 +/- 3% greater than expected). Both the GXXXG and AXXXA motifs occur frequently in predicted alpha-helices from 24 fully sequenced genomes. Occurrence of the AXXXA motif is enhanced to a greater extent in thermophiles than in mesophiles, suggesting that helical interaction based on the AXXXA motif may be a common mechanism of thermostability in protein structures. We conclude that the GXXXG sequence motif stabilizes helix-helix interactions in proteins, and that the AXXXA sequence motif also stabilizes the folded state of proteins.     

The crystal structure of NAD(P)H oxidase from Lactobacillus sanfranciscensis: insights into the conversion of O2 into two water molecules by the flavoenzyme

The FAD-dependent NAD(P)H oxidase from Lactobacillus sanfrancisensis (L.san-Nox2) catalyzes the oxidation of 2 equivalents of either NADH or NADPH and reduces 1 equivalent of O(2) to yield 2 equivalents of water. During steady-state turnover only 0.5% of the reducing equivalents are detected in solution as hydrogen peroxide, suggesting that it is not released from the enzyme after the oxidation of the first equivalent of NAD(P)H and reaction with O(2). Here we report the crystal structure of L.san-Nox2 to 1.8 A resolution. The enzyme crystallizes as a dimer with each monomer consisting of a FAD binding domain (residues 1-120), a NAD(P)H binding domain (residues 150-250), and a dimerization domain (residues 325-451). The electron density for the redox-active Cys42 residue located adjacent to the si-face FAD is consistent with oxidation to the sulfenic acid (Cys-SOH) state. The side chain of Cys42 is also observed in two conformations; in one the sulfenic acid is hydrogen bonded to His10 and in the other it hydrogen bonds with the FAD O2' atom. Surprisingly, the NAD(P)H binding domains each contain an ADP ligand as established by electron density maps and MALDI-TOF analysis of the ligands released from heat-denatured enzyme. The ADP ligand copurifies with the enzyme, and its presence does not inhibit enzyme activity. Consequently, we hypothesize that either NADPH or NADH substrates bind via a long channel that extends from the enzyme exterior and terminates at the FAD re-face. A homology model of the NADH oxidase from Lactococcus lactis (L.lac-Nox2) was also generated using the crystal structure of L.san-Nox2, which reveals several important similarities and differences between the two enzymes. HPLC analysis of ligands released from denatured L.lac-Nox2 indicates that it does not bind ADP, which correlates with the specificity of the enzyme for oxidation of NADH.     

Cofactory: Sequence‐based prediction of cofactor specificity of Rossmann folds

Obtaining optimal cofactor balance to drive production is a challenge in metabolically engineered microbial production strains. To facilitate identification of heterologous enzymes with desirable altered cofactor requirements from native content, we have developed Cofactory, a method for prediction of enzyme cofactor specificity using only primary amino acid sequence information. The algorithm identifies potential cofactor binding Rossmann folds and predicts the specificity for the cofactors FAD(H₂), NAD(H), and NADP(H). The Rossmann fold sequence search is carried out using hidden Markov models whereas artificial neural networks are used for specificity prediction. Training was carried out using experimental data from protein–cofactor structure complexes. The overall performance was benchmarked against an independent evaluation set obtaining Matthews correlation coefficients of 0.94, 0.79, and 0.65 for FAD(H₂), NAD(H), and NADP(H), respectively. The Cofactory method is made publicly available at http://www.cbs.dtu.dk/services/Cofactory . Proteins 2014; 82:1819–1828. © 2014 Wiley Periodicals, Inc.     

Interaction Signatures Stabilizing the NAD(P)-Binding Rossmann Fold: A Structure Network Approach

The fidelity of the folding pathways being encoded in the amino acid sequence is met with challenge in instances where proteins with no sequence homology, performing different functions and no apparent evolutionary linkage, adopt a similar fold. The problem stated otherwise is that a limited fold space is available to a repertoire of diverse sequences. The key question is what factors lead to the formation of a fold from diverse sequences. Here, with the NAD(P)-binding Rossmann fold domains as a case study and using the concepts of network theory, we have unveiled the consensus structural features that drive the formation of this fold. We have proposed a graph theoretic formalism to capture the structural details in terms of the conserved atomic interactions in global milieu, and hence extract the essential topological features from diverse sequences. A unified mathematical representation of the different structures together with a judicious concoction of several network parameters enabled us to probe into the structural features driving the adoption of the NAD(P)-binding Rossmann fold. The atomic interactions at key positions seem to be better conserved in proteins, as compared to the residues participating in these interactions. We propose a “spatial motif” and several “fold specific hot spots” that form the signature structural blueprints of the NAD(P)-binding Rossmann fold domain. Excellent agreement of our data with previous experimental and theoretical studies validates the robustness and validity of the approach. Additionally, comparison of our results with statistical coupling analysis (SCA) provides further support. The methodology proposed here is general and can be applied to similar problems of interest.     

Solution structure of the cysteine-rich domain of the Escherichia coli chaperone protein DnaJ

The solution structure of the cysteine-rich (CR) domain of Escherichia coli DnaJ has been solved by NMR methods. The structure of a 79 residue CR domain construct shows a novel fold with an overall V-shaped extended beta-hairpin topology. The CR domain is characterized by four C-X-X-C-X-G-X-G sequence motifs that bind two zinc ions. Residues in these two zinc modules show strong similarities in the grouping of resonances in the (15)N-(1)H HSQC spectrum and display pseudo-symmetry of the motifs in the calculated structures. The conformation of the cysteine residues coordinated to the zinc ion resembles that of the rubredoxin-knuckle, but there are significant differences in hydrogen bonding patterns in the two motifs. Zinc (15)N-(1)H HSQC titrations indicate that the fold of the isolated DnaJ CR domain is zinc-dependent and that one zinc module folds before the other. The C-X-X-C-X-G-X-G sequence motif is highly conserved in CR domains from a wide variety of species. The three-dimensional structure of the E. coli CR domain indicates that this sequence conservation is likely to result in a conserved structural motif.     

New sequence motifs in flavoproteins: evidence for common ancestry and tools to predict structure

We describe two new sequence motifs, present in several families of flavoproteins. The "GG motif" (RxGGRxxS/T) is found shortly after the betaalphabetadinucleotide-binding motif (DBM) in L-amino acid oxidases, achacin and aplysianin-A, monoamine oxidases, corticosteroid-binding proteins, and tryptophan 2-monooxygenases. Other disperse sequence similarities between these families suggest a common origin. A GG motif is also found in protoporphyrinogen oxidase and carotenoid desaturases and, reduced to the central GG doublet, in the THI4 protein, dTDP-4-dehydrorhamnose reductase, soluble fumarate reductase, steroid dehydrogenases, Rab GDP-dissociation inhibitor, and in most flavoproteins with two dinucleotide-binding domains (glutathione reductase, glutamate synthase, flavin-containing monooxygenase, trimethylamine dehydrogenase...). In the latter families, an "ATG motif" (oxhhhATG) is found in both the FAD- and NAD(P)H-binding domains, forming the fourth beta-strand of the Rossman fold and the connecting loop. On the basis of these and previously described motifs, we present a classification of dinucleotide-binding proteins that could also serve as an evolutionary scheme. Like the DBM, the ATG motif appears to predate the divergence of NAD(P)H- and FAD-binding proteins. We propose that flavoproteins have evolved from a well-differentiated NAD(P)H-binding protein. The bulk of the substrate-binding domain was formed by an insertion after the fourth beta-strand, either of a closely related NAD(P)H-binding domain or of a domain of completely different origin.     

ATG8 Family Proteins Act as Scaffolds for Assembly of the ULK Complex: SEQUENCE REQUIREMENTS FOR LC3-INTERACTING REGION (LIR) MOTIFS*

Autophagy is a lysosome-dependent degradation system conserved among eukaryotes. The mammalian Atg1 homologues, Unc-51 like kinase (ULK) 1 and 2, are multifunctional proteins with roles in autophagy, neurite outgrowth, and vesicle transport. The mammalian ULK complex involved in autophagy consists of ULK1, ULK2, ATG13, FIP200, and ATG101. We have used pulldown and peptide array overlay assays to study interactions between the ULK complex and six different ATG8 family proteins. Strikingly, in addition to ULK1 and ULK2, ATG13 and FIP200 interacted with human ATG8 proteins, all with strong preference for the GABARAP subfamily. Similarly, yeast and Drosophila Atg1 interacted with their respective Atg8 proteins, demonstrating the evolutionary conservation of the interaction. Use of peptide arrays allowed precise mapping of the functional LIR motifs, and two-dimensional scans of the ULK1 and ATG13 LIR motifs revealed which substitutions that were tolerated. This information, combined with an analysis of known LIR motifs, provides us with a clearer picture of sequence requirements for LIR motifs. In addition to the known requirements of the aromatic and hydrophobic residues of the core motif, we found the interactions to depend strongly on acidic residues surrounding the central core LIR motifs. A preference for either a hydrophobic residue or an acidic residue following the aromatic residue in the LIR motif is also evident. Importantly, the LIR motif is required for starvation-induced association of ULK1 with autophagosomes. Our data suggest that ATG8 proteins act as scaffolds for assembly of the ULK complex at the phagophore.     

Mutational analysis of the role of GXXXG motif in the function of human organic anion transporter 1 (hOAT1)

Human organic anion transporter hOAT1 plays a critical role in the body disposition of environmental toxins and clinically important drugs including anti-HIV therapeutics, anti-tumor drugs, antibiotics, anti-hypertensives, and anti-inflammatories. hOAT1 has two GXXXG motifs in its transmembrane domains 2 and 5, a motif linked to the protein processing and oligomerization of other proteins. In the current study, we substituted glycine of these GXXXG motifs with alanine and evaluated the effect of such mutations on the expression and function of hOAT1. Mutations of GXXXG motif in the transmembrane domain 2 resulted in mutants G144A and G148A, both of which had no transport activity due to complete loss in the surface and total cell expression of the transporter protein. Treatment of G144A- and G148A-expressing cells with proteasomal inhibitor resulted in the recovery of ER-resident immature form of hOAT1, but not its surface-resident mature form, whereas treatment of these cells with lysosomal inhibitor had no effect on the expression of the mutant transporters. Mutations of GXXXG motif in the transmembrane domain 5 resulted in mutants G223A and G227A, among which only G227 had dramatic reduction of transport activity due to dramatic loss in the surface and total cell expression of the transporter. The reduction in the surface expression of G227 was consistent with the decrease in maximum transport velocity Vmax. Treatment of G227A-expressing cells with proteasomal inhibitor or lysosomal inhibitor resulted in partial recovery of both the immature form and the mature form of hOAT1 in the total cell extracts. However, such partial recovery of the mature form in total cell extracts did not lead to the partial recovery of surface expression and function of the transporter. Our data suggest that the GXXXG motifs in transmembrane domains 2 and 5 play critical roles in the stability of hOAT1.