Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Molecular cloning of the nahG gene encoding salicylate hydroxylase from Pseudomonas fluorescens

A gene encoding the salicylate hydroxylase was cloned from the genomic DNA of Pseudomonas fluorescens SME11. The DNA fragment containing the nahG gene for the salicylate hydroxylase was mapped with restriction endonucleases and sequenced. The DNA fragment contained an ORF of 1,305 bp encoding a polypeptide of 434 amino acid residues. The nucleotide and amino acid sequences of the salicylate hydroxylase revealed several conserved regions with those of the enzyme encoded in P. putida PpG7: The homology of the nucleotide sequence is 83% and that of amino acid sequence is 72%. We found large conserved regions of the amino acid sequence at FAD and NADH binding regions. The FAD binding site is located at the amino terminal region and a lysine residue functions as a NADH-binding site.     

Involvement of Glycine and Aspartate Residues in the Binding Capacity of FAD in the NADH Dehydrogenase from an Alkaliphilic <Emphasis Type="Italic">Bacillus</Emphasis>

There is a region exhibiting a similarity of amino acid sequence near the carboxyl-terminal segment of each FAD-containing oxidoreductase. In this region, four amino acid residues—Thr, Ala, Gly, and Asp—are highly conserved. To determine the involvement of the four amino acid residues (Thr-469, Ala-476, Gly-478, and Asp-479) in the activity of NADH dehydrogenase of an alkaliphilic Bacillus, mutations of these amino acid residues were conducted. In spite of high conservation, mutations of Thr-469 and Ala-476 to Ala and Ser, respectively, did not lead to a critical loss of enzyme activity. However, mutations of Gly-478 and Asp-479 to Ala caused a complete loss of the activity, which appears to result from the loss of binding capacity of FAD. Received: 3 July 2002 / Accepted: 29 July 2002     

Involvement of glycine and aspartate residues in the binding capacity of FAD in the NADH dehydrogenase from an alkaliphilic Bacillus

There is a region exhibiting a similarity of amino acid sequence near the carboxyl-terminal segment of each FAD-containing oxidoreductase. In this region, four amino acid residues-Thr, Ala, Gly, and Asp-are highly conserved. To determine the involvement of the four amino acid residues (Thr-469, Ala-476, Gly-478, and Asp-479) in the activity of NADH dehydrogenase of an alkaliphilic Bacillus, mutations of these amino acid residues were conducted. In spite of high conservation, mutations of Thr-469 and Ala-476 to Ala and Ser, respectively, did not lead to a critical loss of enzyme activity. However, mutations of Gly-478 and Asp-479 to Ala caused a complete loss of the activity, which appears to result from the loss of binding capacity of FAD.     

Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.     

Multiple forms of liver microsomal flavin-containing monooxygenases: complete covalent structure of form 2

Hepatic flavin-containing monooxygenases catalyze NADPH-dependent oxygenation of a wide variety of drugs that possess a nucleophilic heteroatom. Two forms of these enzymes (form 1 and 2) have been isolated from rabbit liver microsomes and partially characterized (Ozols, J., 1989, Biochem. Biophys. Res. Commun. 163, 49-55). The complete amino acid sequence of form 2 is presented here. Sequence determination was achieved by pulsed liquid-phase and solid-phase sequencing of 40 peptides generated by chemical and enzymatic cleavages, including CNBr cleavage of tryptophanyl residues. Form 2 monooxygenase contains 533 amino acid residues and has a molecular weight of 60,089. The COOH terminus of this enzyme is very hydrophobic and presumably functions to anchor the protein to the membrane. Form 2 is readily degraded, since a form lacking residues 1 to 278 and a form without the COOH-terminal segment were also isolated from solubilized membrane preparations. The amino acid sequence of form 2 is 52% identical to that of form 1 and shows 55% identity to the sequence of rabbit lung monooxygenase derived from the cDNA data. The putative FAD and NADP binding segments around residues 9 and 190 are conserved in all three forms. Three variable segments can also be identified in these isoforms. These are residues 308 to 321, residues 408 to 421, and the membrane binding domain, residues 505 to 533. A comparison of the presently limited amino acid sequence data of flavin-containing monooxygenases (FMOs) implies that a particular FMO in different mammalian species may be very similar, but isozymes within a species may exhibit more extensive variability with respect to homology and catalytic activity. This study documents the structural diversity of a second hepatic FMO from rabbit liver and establishes this class of drug-metabolizing enzymes as a family of related proteins.     

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.     

Role for a conserved structural motif in assembly of a class I aminoacyl-tRNA synthetase active site

The catalytic domains of class I aminoacyl-tRNA synthetases are built around a conserved Rossmann nucleotide binding fold, with additional polypeptide domains responsible for tRNA binding or hydrolytic editing of misacylated substrates. Structural comparisons identified a conserved motif bridging the catalytic and anticodon binding domains of class Ia and Ib enzymes. This stem contact fold (SCF) has been proposed to globally orient each enzyme's cognate tRNA by interacting with the inner corner of the L-shaped tRNA. Despite the structural similarity of the SCF among class Ia/Ib enzymes, the sequence conservation is low. We replaced amino acids of the MetRS SCF with portions of the structurally similar glutaminyl-tRNA synthetase (GlnRS) motif or with alanine residues. Chimeric variants retained significant tRNA methionylation activity, indicating that structural integrity of the helix-turn-strand-helix motif contributes more to tRNA aminoacylation than does amino acid identity. In contrast, chimeras were significantly reduced in methionyl adenylate synthesis, suggesting a role for the SCF in formation of a structured active site domain. A highly conserved aspartic acid within the MetRS SCF is proposed to make an electrostatic interaction with an active site lysine; these residues were replaced with alanines or conservative substitutions. Both methionyl adenylate formation and methionine transfer were impaired, and activity was not significantly recovered by making the compensatory double substitution.     

Identification of a novel conserved sequence motif in flavoprotein hydroxylases with a putative dual function in FAD/NAD(P)H binding.

A novel conserved sequence motif has been located among the flavoprotein hydroxylases. Based on the crystal structure and site-directed mutagenesis studies of p-hydroxybenzoate hydroxylase (PHBH) from Pseudomonas fluorescens, this amino acid fingerprint sequence is proposed to play a dual function in both FAD and NAD(P)H binding. In PHBH, the novel sequence motif (residues 153-166) includes strand A4 and the N-terminal part of helix H7. The conserved amino acids Asp 159, Gly 160, and Arg 166 are necessary for maintaining the structure. The backbone oxygen of Cys 158 and backbone nitrogens of Gly 160 and Phe 161 interact indirectly with the pyrophosphate moiety of FAD, whereas it is known from mutagenesis studies that the side chain of the moderately conserved His 162 is involved in NADPH binding.     

Isolation and characterization of the gene from Pseudomonas syringae pv. phaseolicola encoding the phaseolotoxin-insensitive ornithine carbamoyltransferase

Summary The gene coding for the phaseolotoxin-insensitive ornithine carbamoyltransferase (OCTase) fromPseudomonas syringae pv.phaseolicola has been cloned and sequenced. The gene has a deduced coding capacity for a polypeptide with a calculated M, of 36520 daltons. Comparison of the amino acid sequence of the OCTase enzymes encoded by theP. aeruginosa argF and theEscherichia coli argI andargF genes with the deduced sequence of the newly identified gene shows that 79 amino acid residues are strictly conserved in all four polypeptides; among these 7 out of 9 residues are involved in enzyme function. Of three amino acid regions that have been implicated in substrate binding or catalysis, two are strictly conserved, and the third involved in carbamoylphosphate binding differs. This correlates well with published data showing that phaseolotoxin competes for the carbamoylphosphate binding site in the phaseolotoxin-sensitive OCTases. We propose that the gene be namedargK.     

Binding of the Covalent Flavin Assembly Factor to the Flavoprotein Subunit of Complex II

Escherichia coli harbors two highly conserved homologs of the essential mitochondrial respiratory complex II (succinate:ubiquinone oxidoreductase). Aerobically the bacterium synthesizes succinate:quinone reductase as part of its respiratory chain, whereas under microaerophilic conditions, the quinol:fumarate reductase can be utilized. All complex II enzymes harbor a covalently bound FAD co-factor that is essential for their ability to oxidize succinate. In eukaryotes and many bacteria, assembly of the covalent flavin linkage is facilitated by a small protein assembly factor, termed SdhE in E. coli. How SdhE assists with formation of the covalent flavin bond and how it binds the flavoprotein subunit of complex II remain unknown. Using photo-cross-linking, we report the interaction site between the flavoprotein of complex II and the SdhE assembly factor. These data indicate that SdhE binds to the flavoprotein between two independently folded domains and that this binding mode likely influences the interdomain orientation. In so doing, SdhE likely orients amino acid residues near the dicarboxylate and FAD binding site, which facilitates formation of the covalent flavin linkage. These studies identify how the conserved SdhE assembly factor and its homologs participate in complex II maturation.     

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 (379)NXX[S/T]HX[S/T]QDXXXT(391). 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.     

Identification of ADP in the iron-sulfur flavoprotein trimethylamine dehydrogenase

Analysis of the 2.4-A resolution electron density map of trimethylamine dehydrogenase has revealed the unexpected presence of one molecule of ADP/subunit. This binding has been confirmed chemically. The binding site is located at the analogous position of the ADP moiety of FAD in glutathione reductase, the FAD and NADPH binding domains of which resemble two of the domains of trimethylamine dehydrogenase. Comparison of the environments of the ADP moieties in the two proteins indicates that 32 residues in 6 peptides are in equivalent positions with a root mean square deviation for C alpha positions of 1.11 A. Twelve of these amino acids are identical, based on the electron density-derived "x-ray" sequence of trimethylamine dehydrogenase. Detailed analysis of the environment of the ADP moiety indicates that most of the conserved residues are not in direct contact with the cofactor. Some of them probably represent the "fingerprint" of the beta alpha beta binding fold found in dinucleotide binding proteins, but the remaining conserved residues may indicate a closer evolutionary relationship between these two proteins.     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

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.     

Evolutionary information hidden in a single protein structure

The knowledge of conserved sequences in proteins is valuable in identifying functionally or structurally important residues. Generating the conservation profile of a sequence requires aligning families of homologous sequences and having knowledge of their evolutionary relationships. Here, we report that the conservation profile at the residue level can be quantitatively derived from a single protein structure with only backbone information. We found that the reciprocal packing density profiles of protein structures closely resemble their sequence conservation profiles. For a set of 554 nonhomologous enzymes, 74% (408/554) of the proteins have a correlation coefficient > 0.5 between these two profiles. Our results indicate that the three-dimensional structure, instead of being a mere scaffold for positioning amino acid residues, exerts such strong evolutionary constraints on the residues of the protein that its profile of sequence conservation essentially reflects that of its structural characteristics.     

Functional properties of the alternative NADH:ubiquinone oxidoreductase from E. coli through comparative 3-D modelling

The alternative NADH:ubiquinone oxidoreductase (NDH-2) from Escherichia coli is a membrane protein playing a prominent role in respiration by linking the reduction of NADH to the quinone pool. Remote sequence similarity reveals an evolutionary relation between alternative NADH:quinone oxidoreductases and the SCOP-family "FAD/NAD-linked reductases". We have created a structural model for NDH-2 from E. coli through comparative modelling onto a template from this family. Combined analysis of our model and sequence conservation allowed us to include the cofactor FAD and the substrate NADH in atomic detail. Furthermore, we propose the most plausible orientation of NDH-2 relative to the membrane and specify a region of the protein potentially involved in ubiquinone binding.     

Sequence homology around the biotin-binding site of human propionyl-CoA carboxylase and pyruvate carboxylase

Biotin-dependent carboxylases require covalently bound biotin for enzymatic activity. The biotin is attached through a lysine residue, which in a number of bacterial, avian, and mammalian carboxylases, is found within the conserved sequence Ala-Met-Lys-Met. We have determined the partial nucleotide sequence of cDNA clones for human propionyl-CoA carboxylase and pyruvate carboxylase. The predicted amino acid sequence of both these proteins contains the conserved tetrapeptide 35 residues from the carboxy terminus. In addition, both proteins contain the tripeptide, Pro-Met-Pro, 26 residues toward the amino terminus from the biotin attachment site. The overall amino acid homology through this region is 43%. Similar findings have been made for the biotin-containing polypeptides of transcarboxylase of Propionibacterium shermanii and acetyl-CoA carboxylase of Escherichia coli (W. L. Maloy, B. U. Bowien, G. K. Zwolinski, K. G. Kumar, and H. G. Wood (1979) J. Biol. Chem. 254, 11615-11622). The implications of this sequence conservation with regard to the function and evolution of biotin-dependent carboxylases is discussed. We propose that the 60 amino acids surrounding the biotin site are bounded by a proline "hinge" and the carboxy terminus has remained conserved as a result of constraints imposed by biotinylation of the enzyme.     

Primary structure and eubacterial relationships of the pyruvate:Ferredoxin oxidoreductase of the amitochondriate eukaryoteTrichomonas vaginalis

In the eukaryotic unicellular organismTrichomonas vaginalis a key step of energy metabolism, the oxidative decarboxylation of pyruvate with the formation of acetyl-CoA, is catalyzed by the iron-sulfur protein pyruvate:ferredoxin oxidoreductase (PFO) and not by the almost-ubiquitous pyruvate dehydrogenase multienzyme complex. This enzyme is localized in the hydrogenosome, an organelle bounded by a double membrane. PFO and its closely related homolog, pyruvate: flavodoxin oxidoreductase, are enzymes found in a number of archaebacteria and eubacteria. The presence of these enzymes in eukaryotes is restricted, however, to a few amitochondriate groups. To gain more insight into the evolutionary relationships ofT. vaginalis PFO we determined the primary structure of its two genes (pfoA andpfoB). The deduced amino acid sequences showed 95% positional identity. Motifs implicated in related enzymes in liganding the Fe-S centers and thiamine pyrophosphate were well conserved. TheT. vaginalis PFOs were found to be homologous to eubacterial pyruvate: flavodoxin oxidoreductases and showed about 40% amino acid identity to these enzymes over their entire length. Lack of eubacterial PFO sequences precluded a comparison.pfoA andpfoB revealed a greater distance from related enzymes of Archaebacteria. The conceptual translation of the nucleotide sequences predicted an amino-terminal pentapeptide not present in the mature protein. This processed leader sequence was similar to but shorter than leader sequences noted in other hydrogenosomal proteins. These sequences are assumed to be involved in organellar targeting and import. The results underscore the unusual characteristics ofT. vaginalis metabolism and of their hydrogenosomes. They also suggest that in its energy metabolismT. vaginalis is closer to eubacteria than archaebacteria.Abbreviations PCR DNA polymerase chain reaction - PDH pyruvate dehydrogenase - PFO pyruvate:ferredoxin oxidoreductase - TPP thiamine pyrophosphate Correspondence to: M. Müller