Identification of active sites in amidase: Evolutionary relationship between amide bond- and peptide bond-cleaving enzymes

Mainly based on various inhibitor studies previously performed, amidases came to be regarded as sulfhydryl enzymes. Not completely satisfied with this generally accepted interpretation, we performed a series of site-directed mutagenesis studies on one particular amidase of Rhodococcus rhodochrous J1 that was involved in its nitrile metabolism. For these experiments, the recombinant amidase was produced as the inclusion body in Escherichia coli to greatly facilitate its recovery and subsequent purification. With regard to the presumptive active site residue Cys203, a Cys203 → Ala mutant enzyme still retained 11.5% of the original specific activity. In sharp contrast, substitutions in certain other positions in the neighborhood of Cys203 had a far more dramatic effect on the amidase. Glutamic acid substitution of Asp191 reduced the specific activity of the mutant enzyme to 1.33% of the wild-type activity. Furthermore, Asp191 → Asn substitution as well as Ser195 → Ala substitution completely abolished the specific activity. It would thus appear that, among various conserved residues residing within the so-called signature sequence common to all amidases, the real active site residues are Asp191 and Ser195 rather than Cys203. Inasmuch as an amide bond (CO-NH2) in the amide substrate is not too far structurally removed from a peptide bond (CO-NH-), the signature sequences of various amidases were compared with the active site sequences of various types of proteases. It was found that aspartic acid and serine residues corresponding to Asp191 and Ser195 of the Rhodococcus amidase are present within the active site sequences of aspartic proteinases, thus suggesting the evolutionary relationship between the two.

Cellulose-binding domains promote hydrolysis of different sites on crystalline cellulose

The cohesin-dockerin interaction in Clostridium thermocellum cellulosome mediates the tight binding of cellulolytic enzymes to the cellulosome-integrating protein CipA. Here, this interaction was used to study the effect of different cellulose-binding domains (CBDs) on the enzymatic activity of C. thermocellum endoglucanase CelD (1,4-β-d endoglucanase, EC3.2.1.4) toward various cellulosic substrates. The seventh cohesin domain of CipA was fused to CBDs originating from the Trichoderma reesei cellobiohydrolases I and II (CBDCBH1 and CBDCBH2) (1,4-β-d glucan-cellobiohydrolase, EC3.2.1.91), from the Cellulomonas fimi xylanase/exoglucanase Cex (CBDCex) (β-1,4-d glucanase, EC3.2.1.8), and from C. thermocellum CipA (CBDCipA). The CBD-cohesin hybrids interacted with the dockerin domain of CelD, leading to the formation of CelD-CBD complexes. Each of the CBDs increased the fraction of cellulose accessible to hydrolysis by CelD in the order CBDCBH1 < CBDCBH2 ≈ CBDCex < CBDCipA. In all cases, the extent of hydrolysis was limited by the disappearance of sites accessible to CelD. Addition of a batch of fresh cellulose after completion of the reaction resulted in a new burst of activity, proving the reversible binding of the intact complexes despite the apparent binding irreversibility of some CBDs. Furthermore, burst of activity also was observed upon adding new batches of CelD–CBD complexes that contained a CBD differing from the first one. This complementation between different CBDs suggests that the sites made available for hydrolysis by each of the CBDs are at least partially nonoverlapping. The only exception was CBDCipA, whose sites appeared to overlap all of the other sites.

Induction and repression of DAN1 and the family of anaerobic mannoprotein genes in Saccharomyces cerevisiae occurs through a complex array of regulatory sites

The DAN/TIR mannoprotein genes of Saccharomyces cerevisiae (DAN1, DAN2, DAN3, DAN4, TIR1, TIR2, TIR3 and TIR4) are expressed in anaerobic cells while the predominant cell wall proteins Cwp1 and Cwp2 are down-regulated. Elements involved in activation and repression of the DAN/TIR genes were defined in this study, using the DAN1 promoter as a model. Nested deletions in a DAN1/lacZ reporter pinpointed regions carrying activation and repression elements. Inspection revealed two consensus sequences subsequently shown to be independent anaerobic response elements (AR1, consensus TCGTTYAG; AR2, consensus AAAAATTGTTGA). AR1 is found in all of the DAN/TIR promoters; AR2 is found in DAN1, DAN2 and DAN3. A 120 bp segment carrying two copies of AR1 preferentially activated transcription of lacZ under anaerobic conditions. A fusion of three synthetic copies of AR1 to MEL1 was also expressed anaerobically. Mutations in either AR1 site within the 120 bp segment caused a drastic loss of expression, indicating that both are necessary for activation and implying cooperativity between adjacent transcriptional activation complexes. A single AR2 site carried on a 46 bp fragment from the DAN1 promoter activated lacZ transcription under anaerobic conditions, as did a 26 bp synthetic AR2 fragment fused to MEL1. Nucleotide substitutions within the AR2 sequence eliminated the activity of the 46 bp segment. Ablation of the AR2 sequences in the full promoter caused a partial reduction of expression. The presence of the ATTGTT core (recognized by HMG proteins) in the AR2 sequence suggests that an HMG protein may activate through AR2. One region was implicated in aerobic repression of DAN1. It contains sites for the heme-induced Mot3 and Rox1 repressors.

Mapping of contact sites in complex formation between transducin and light-activated rhodopsin by covalent crosslinking: Use of a photoactivatable reagent

Interaction of light-activated rhodopsin with transducin (T) is the first event in visual signal transduction. We use covalent crosslinking approaches to map the contact sites in interaction between the two proteins. Here we use a photoactivatable reagent, N-[(2-pyridyldithio)-ethyl], 4-azido salicylamide. The reagent is attached to the SH group of cytoplasmic monocysteine rhodopsin mutants by a disulfide-exchange reaction with the pyridylthio group, and the derivatized rhodopsin then is complexed with T by illumination at λ >495 nm. Subsequent irradiation of the complex at λ310 nm generates covalent crosslinks between the two proteins. Crosslinking was demonstrated between T and a number of single cysteine rhodopsin mutants. However, sites of crosslinks were investigated in detail only between T and the rhodopsin mutant S240C (cytoplasmic loop V-VI). Crosslinking occurred predominantly with Tα. For identification of the sites of crosslinks in Tα, the strategy used involved: (i) derivatization of all of the free cysteines in the crosslinked proteins with N-ethylmaleimide; (ii) reduction of the disulfide bond linking the two proteins and isolation of all of the Tα species carrying the crosslinked moiety with a free SH group; (iii) adduct formation of the latter with the N-maleimide moiety of the reagent, maleimido-butyryl-biocytin, containing a biotinyl group; (iv) trypsin degradation of the resulting Tα derivatives and isolation of Tα peptides carrying maleimido-butyryl-biocytin by avidin-agarose chromatography; and (v) identification of the isolated peptides by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. We found that crosslinking occurred mainly to two C-terminal peptides in Tα containing the amino acid sequences 310–313 and 342–345.

Mapping of contact sites in complex formation between light-activated rhodopsin and transducin by covalent crosslinking: Use of a chemically preactivated reagent

Contact sites in interaction between light-activated rhodopsin and transducin (T) have been investigated by using a chemically preactivated crosslinking reagent, N-succinimidyl 3-(2-pyridyldithio)propionate. The 3 propionyl-N-succinimidyl group in the reagent was attached by a disulfide exchange reaction to rhodopsin mutants containing single reactive cysteine groups in the cytoplasmic loops. Complex formation between the derivatized rhodopsin mutants and T was carried out by illumination at λ > 495 nm. Subsequent increase in pH (from 6 to 7.5 or higher) of the complex resulted in crosslinking of rhodopsin to the Tα subunit. Crosslinking to Tα was demonstrated for the rhodopsin mutants K141C, S240C, and K248C, and the crosslinked sites in Tα were identified for the rhodopsin mutant S240C. The peptides carrying the crosslinking moiety were isolated from the trypsin-digested peptide mixture, and their identification was carried out by matrix-assisted laser desorption ionization–time of flight mass spectrometry. The main site of crosslinking is within the peptide sequence, Leu-19–Arg-28 at the N-terminal region of Tα. The total results show that both the N and the C termini of Tα are in close vicinity to the third cytoplasmic loop of rhodopsin in the complex between rhodopsin and T.

Random mutagenesis of Thermus aquaticus DNA polymerase I: concordance of immutable sites in vivo with the crystal structure.

Expression of Thermus aquaticus (Taq) DNA polymerase I (pol I) in Escherichia, coli complements the growth defect caused by a temperature-sensitive mutation in the host pol I. We replaced the nucleotide sequence encoding amino acids 659-671 of the O-helix of Taq DNA pol I, corresponding to the substrate binding site, with an oligonucleotide containing random nucleotides. Functional Taq pol I mutants were selected based on colony formation at the nonpermissive temperature. By using a library with 9% random substitutions at each of 39 positions, we identified 61 active Taq pol I mutants, each of which contained from one to four amino acid substitutions. Some amino acids, such as alanine-661 and threonine-664, were tolerant of several or even many diverse replacements. In contrast, no replacements or only conservative replacements were identified at arginine-659, lysine-663, and tyrosine-671. By using a library with totally random nucleotides at five different codons (arginine-659, arginine-660, lysine-663, phenylalanine-667, and glycine-668), we confirmed that arginine-659 and lysine-663 were immutable, and observed that only tyrosine substituted for phenylalanine-667. The two immutable residues and the two residues that tolerate only highly conservative replacements lie on the side of O-helix facing the incoming deoxynucleoside triphosphate, as determined by x-ray analysis. Thus, we offer a new approach to assess concordance of the active conformation of an enzyme, as interpreted from the crystal structure, with the active conformation inferred from in vivo function.

Negative electrostatic surface potential of protein sites specific for anionic ligands.

Determination of the crystal structure of an "open" unliganded active mutant (T141D) form of the Escherichia coli phosphate receptor for active transport has allowed calculation of the electrostatic surface potential for it and two other comparably modeled receptor structures (wild type and D137N). A discovery of considerable implication is the intensely negative potential of the phosphate-binding cleft. We report similar findings for a sulfate transport receptor, a DNA-binding protein, and, even more dramatically, redox proteins. Evidently, for proteins such as these, which rely almost exclusively on hydrogen bonding for anion interactions and electrostatic balance, a noncomplementary surface potential is not a barrier to binding. Moreover, experimental results show that the exquisite specificity and high affinity of the phosphate and sulfate receptors for unions are insensitive to modulations of charge potential, but extremely sensitive to conditions that leave a hydrogen bond donor or acceptor unpaired.

Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis.

Competing hypotheses seek to explain the evolution of oxygenic and anoxygenic processes of photosynthesis. Since chlorophyll is less reduced and precedes bacteriochlorophyll on the modern biosynthetic pathway, it has been proposed that chlorophyll preceded bacteriochlorophyll in its evolution. However, recent analyses of nucleotide sequences that encode chlorophyll and bacteriochlorophyll biosynthetic enzymes appear to provide support for an alternative hypothesis. This is that the evolution of bacteriochlorophyll occurred earlier than the evolution of chlorophyll. Here we demonstrate that the presence of invariant sites in sequence datasets leads to inconsistency in tree building (including maximum-likelihood methods). Homologous sequences with different biological functions often share invariant sites at the same nucleotide positions. However, different constraints can also result in additional invariant sites unique to the genes, which have specific and different biological functions. Consequently, the distribution of these sites can be uneven between the different types of homologous genes. The presence of invariant sites, shared by related biosynthetic genes as well as those unique to only some of these genes, has misled the recent evolutionary analysis of oxygenic and anoxygenic photosynthetic pigments. We evaluate an alternative scheme for the evolution of chlorophyll and bacteriochlorophyll.

Mutagenesis of palmitoylation sites in endothelial nitric oxide synthase identifies a novel motif for dual acylation and subcellular targeting.

The endothelial nitric oxide synthase (ec-NOS) plays a key role in the transduction of signals from the bloodstream to the underlying smooth muscle. ecNOS undergoes a complex series of covalent modifications, including myristoylation and palmitoylation, which appear to play a role in ecNOS membrane association. Mutagenesis of the myristoylation site, which prevents both myristoylation and palmitoylation, blocks ecNOS targeting to cell membranes. Further, as described for some G-protein alpha subunits, both membrane association and palmitoylation of ecNOS are dynamically regulated: in response to agonists, the enzyme undergoes partial redistribution to the cell cytosol concomitant with depalmitoylation. To clarify the role of palmitoylation in determining ecNOS subcellular localization, we have constructed palmitoylation-deficient mutants of ecNOS. Serine was substituted for cysteine at two potential palmitoylation sites (Cys-15 and Cys-26) by site-directed mutagenesis. Immunoprecipitation of ecNOS mutants following cDNA transfection and biosynthetic labeling with [3H]palmitate revealed that mutagenesis of either cysteine residue attenuated palmitoylation, whereas replacement of both residues completely eliminated palmitoylation. Analysis of N-terminal deletion mutations of ecNOS demonstrated that the region containing these two cysteine residues is both necessary and sufficient for enzyme palmitoylation. The cysteines thus identified as the palmitoylation sites for ecNOS are separated by an unusual (Gly-Leu)5 sequence and appear to define a sequence motif for dual acylation. We analyzed the subcellular distribution of ecNOS mutants by differential ultracentrifugation and found that mutagenesis of the ecNOS palmitoylation sites markedly reduced membrane association of the enzyme. These results document that ecNOS palmitoylation is an important determinant for the subcellular distribution of ecNOS and identify a new motif for the reversible palmitoylation of signaling proteins.

Relating aromatic hydrocarbon-induced DNA adducts and c-H-ras mutations in mouse skin papillomas: the role of apurinic sites.

Mouse skin tumors contain activated c-H-ras oncogenes, often caused by point mutations at codons 12 and 13 in exon 1 and codons 59 and 61 in exon 2. Mutagenesis by the noncoding apurinic sites can produce G-->T and A-->T transversions by DNA misreplication with more frequent insertion of deoxyadenosine opposite the apurinic site. Papillomas were induced in mouse skin by several aromatic hydrocarbons, and mutations in the c-H-ras gene were determined to elucidate the relationship among DNA adducts, apurinic sites, and ras oncogene mutations. Dibenzo[a,l]pyrene (DB[a,l]P), DB[a,l]P-11,12-dihydrodiol, anti-DB[a,l]P-11,12-diol-13,14-epoxide, DB[a,l]P-8,9-dihydrodiol, 7,12-dimethylbenz[a]anthracene (DMBA), and 1,2,3,4-tetrahydro-DMBA consistently induced a CAA-->CTA mutation in codon 61 of the c-H-ras oncogene. Benzo[a]pyrene induced a GGC-->GTC mutation in codon 13 in 54% of tumors and a CAA-->CTA mutation in codon 61 in 15%. The pattern of mutations induced by each hydrocarbon correlated with its profile of DNA adducts. For example, both DB[a,l]P and DMBA primarily form DNA adducts at the N-3 and/or N-7 of deoxyadenosine that are lost from the DNA by depurination, generating apurinic sites. Thus, these results support the hypothesis that misreplication of unrepaired apurinic sites generated by loss of hydrocarbon-DNA adducts is responsible for transforming mutations leading to papillomas in mouse skin.

Phytosulfokine-α, a sulfated pentapeptide, stimulates the proliferation of rice cells by means of specific high- and low-affinity binding sites

Peptide growth factors were isolated from conditioned medium derived from rice (Oryza sativa L.) suspension cultures and identified to be a sulfated pentapeptide [H-Tyr(SO3H)-Ile-Tyr(SO3H)-Thr-Gln-OH] and its C-terminal-truncated tetrapeptide [H-Tyr(SO3H)-Ile-Tyr(SO3H)-Thr-OH]. These structures were identical to the phytosulfokines originally found in asparagus (Asparagus officinalis L.) mesophyll cultures. The pentapeptide [phytosulfokine-α (PSK-α)] very strongly stimulated colony formation of rice protoplasts at concentrations above 10−8 M, indicating a similar mode of action in rice of phytosulfokines. Binding assays using 35S-labeled PSK-α demonstrated the existence of both high- and low-affinity specific saturable binding sites on the surface of rice cells in suspension. Analysis of [35S]PSK-α binding in differential centrifugation fractions suggested association of the binding with a plasma membrane-enriched fraction. The apparent Kd values for [35S]PSK-α binding were found to be 1 × 10−9 M for the high-affinity type and 1 × 10−7 M for the low-affinity type, with maximal numbers of binding sites of 1 × 104 sites per cell and 1 × 105 sites per cell, respectively. Competition studies with [35S]PSK-α and several synthetic PSK-α analogs demonstrated that only peptides that possesses mitogenic activity can effectively displace the radioligand. These results suggest that a signal transduction pathway mediated by peptide factors is involved in plant cell proliferation.

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22 815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT–AG junctions (22 199 entries) and 0.56% have non-canonical GC–AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC–AG pairs (of which one was an error that corrected to GC–AG), 61 errors corrected to GT–AG canonical pairs, six AT–AC pairs (of which two were errors corrected to AT–AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac.uk/spldb/SpliceDB.html and at http://www.softberry.com/spldb/SpliceDB.html.

Identification of three physically and functionally distinct binding sites for C3b in human complement factor H by deletion mutagenesis.

Human complement factor H controls spontaneous activation of complement in plasma and appears to play a role in distinguishing host cells from activators of the alternative pathway of complement. In both mice and humans, the protein is composed of 20 homologous short consensus repeat (SCR) domains. The size of the protein suggests that portions of the structure outside the known C3b binding site (SCR 1-4) possess a significant biological role. We have expressed the full-length cDNA of factor H in the baculovirus system and have shown the recombinant protein to be fully active. Mutants of this full-length protein have now been prepared, purified, and examined for cofactor activity and binding to C3b and heparin. The results demonstrate (i) that factor H has at least three sites that bind C3b, (ii) that one of these sites is located in SCR domains 1-4, as has been shown by others, (iii) that a second site exists in the domain 6-10 region, (iv) that a third site resides in the SCR 16-20 region, and (v) that two heparin binding sites exist in factor H, one near SCR 13 and another in the SCR 6-10 region. Functional assays demonstrated that only the first C3b site located in SCR 1-4 expresses factor I cofactor activity. Mutant proteins lacking any one of the three C3b binding sites exhibited 6- to 8-fold reductions in affinity for C3b on sheep erythrocytes, indicating that all three sites contribute to the control of complement activation on erythrocytes. The identification of multiple functionally distinct sites on factor H clarifies many of the heretofore unexplainable behaviors of this protein, including the heterogeneous binding of factor H to surface-bound C3b, the effects of trypsin cleavage, and the differential control of complement activation on activators and nonactivators of the alternative pathway of complement.

Presynaptic association of Rad51 protein with selected sites in meiotic chromatin.

Eukaryotic homologs of Escherichia coli Rec-A protein have been shown to form nucleoprotein filaments with single-stranded DNA that recognize homologous sequences in duplex DNA. Several recent reports in four widely diverse species have demonstrated the association of RecA homologs with meiotic prophase chromatin. The current immunocytological study on mouse spermatocytes and oocytes shows that a eukaryotic homolog, Rad5l, associates with a subset of chromatin sites as early as premeiotic S phase, hours before either the appearance of precursors of synaptonemal complexes or the initiation of synapsis. When homologous chromosomes do begin to pair, the Rad5l-associated sequences are sites of initial contact between homologues and of localized DNA synthesis. Distribution of Rad5l foci on the chromatin of fully synapsed bivalents at early pachynema corresponds to an R-band pattern of mitotic chromosomes. R-bands are known to be preferred sites of both synaptic initiation and recombination. The time course of appearance of Rad51 association with chromatin, its distribution, and its interaction with other Rad5l-associated sequences suggests that it plays an important role preselection of sequences and synaptic initiation.

Different cleavage sites are aligned differently in the active site of M1 RNA, the catalytic subunit of Escherichia coli RNase P.

We have studied RNase P RNA (M1 RNA) cleavage of model tRNA precursors that are cleaved at two independent positions. Here we present data demonstrating that cleavage at both sites depends on the 2'-OH immediately 5' of the respective cleavage site. However, we show that the 2-amino group of a guanosine at the cleavage site plays a significant role in cleavage at one of these sites but not at the other. These data suggest that these two cleavage sites are handled differently by the ribozyme. This theory is supported by our finding that the cross-linking pattern between Ml RNA and tRNA precursors carrying 4-thioU showed distinct differences, depending on the location of the 4-thioU relative to the respective cleavage site. These findings lead us to suggest that different cleavage sites are aligned differently in the active site, possibly as a result of different binding modes of a substrate to M1 RNA. We discuss a model in which the interaction between the 3'-terminal "RCCA" motif (first three residues interact) of a tRNA precursor and M1 RNA plays a significant role in this process.