Amino Acid Substitutions in Cold-Adapted Proteins from Halorubrum lacusprofundi,an Extremely Halophilic Microbe from Antarctica

The halophilic Archaeon Halorubrum lacusprofundi, isolated from the perennially cold and hypersaline Deep Lake in Antarctica, was recently sequenced and compared to 12 Haloarchaea from temperate climates by comparative genomics. Amino acid substitutions for 604 H. lacusprofundi proteins belonging to conserved haloarchaeal orthologous groups (cHOGs) were determined and found to occur at 7.85% of positions invariant in proteins from mesophilic Haloarchaea. The following substitutions were observed most frequently: (a) glutamic acid with aspartic acid or alanine; (b) small polar residues with other small polar or non-polar amino acids; (c) small non-polar residues with other small non-polar residues; (d) aromatic residues, especially tryptophan, with other aromatic residues; and (e) some larger polar residues with other similar residues. Amino acid substitutions for a cold-active H. lacusprofundi β-galactosidase were then examined in the context of a homology modeled structure at residues invariant in homologous enzymes from mesophilic Haloarchaea. Similar substitutions were observed as in the genome-wide approach, with the surface accessible regions of β-galactosidase displaying reduced acidity and increased hydrophobicity, and internal regions displaying mainly subtle changes among smaller non-polar and polar residues. These findings are consistent with H. lacusprofundi proteins displaying amino acid substitutions that increase structural flexibility and protein function at low temperature. We discuss the likely mechanisms of protein adaptation to a cold, hypersaline environment on Earth, with possible relevance to life elsewhere.     

Evidence for Production of a New Lantibiotic (Butyrivibriocin OR79A) by the Ruminal Anaerobe Butyrivibrio fibrisolvens OR79: Characterization of the Structural Gene Encoding Butyrivibriocin OR79A

The ruminal anaerobe Butyrivibrio fibrisolvens OR79 produces a bacteriocin-like activity demonstrating a very broad spectrum of activity. An inhibitor was isolated from spent culture fluid by a combination of ammonium sulfate and acidic precipitations, reverse-phase chromatography, and high-resolution gel filtration. N-terminal analysis of the isolated inhibitor yielded a 15-amino-acid sequence (G-N/Q-G/P-V-I-L-X-I-X-H-E-X-S-M-N). Two different amino acid residues were detected in the second and third positions from the N terminus, indicating the presence of two distinct peptides. A gene with significant homology to one combination of the determined N-terminal sequence was cloned, and expression of the gene was confirmed by Northern blotting. The gene (bvi79A) encoded a prepeptide of 47 amino acids and a mature peptide, butyrivibriocin OR79A, of 25 amino acids. Significant sequence homology was found between this peptide and previously reported lantibiotics containing the double-glycine leader peptidase processing site. Immediately downstream of bvi79A was a second, partial open reading frame encoding a peptide with significant homology to proteins which are believed to be involved in the synthesis of lanthionine residues. These findings indicate that the isolated inhibitory peptides represent new lantibiotics. Results from both total and N-terminal amino acid sequencing indicated that the second peptide was identical to butyrivibriocin OR79A except for amino acid substitutions in positions 2 and 3 of the mature lantibiotic. Only a single coding region was detected when restriction enzyme digests of total DNA were probed either with an oligonucleotide based on the 5′ region of bvi79A or with degenerate oligonucleotides based on the predicted sequence of the second peptide.     

Engineered Ssp DnaX inteins for protein splicing with flanking proline residues

Inteins are internal protein sequences capable of catalyzing a protein splicing reaction by self-excising from a precursor protein and simultaneously joining the flanking sequences with a peptide bond. Split inteins have separate pieces (N-intein and C-intein) that reassemble non-covalently to catalyze a protein trans-splicing reaction joining two polypeptides. Protein splicing has become increasingly useful tools in many fields of biological research and biotechnology. However, natural and engineered inteins have failed previously to function when being flanked by proline residue at the −1 or +2 positions, which limits general uses of inteins. In this study, different engineered inteins were tested. We found that engineered Ssp DnaX mini-intein and split inteins could carry out protein splicing with proline at the +2 positions or at both −1 and +2 positions. Under in vivo conditions in E. coli cells, the mini-intein, S1 split intein, and S11 split intein spliced efficiently, whereas the S0 split intein did not splice with proline at both −1 and +2 positions. The S1 and S11 split inteins also trans-spliced efficiently in vitro with proline at the +2 positions or at both −1 and +2 positions, but the S0 split intein trans-spliced inefficiently with proline at the +2 position and did not trans-splice with proline at both −1 and +2 positions. These findings contribute significantly to the toolbox of intein-based technologies by allowing the use of inteins in proteins having proline at the splicing point.     

Two genes encoding gas vacuole proteins in Halobacterium halobium

Summary The archaebacterium Halobacterium halobium contains two related gas vacuole protein-encoding genes (vac). One of these genes encodes a protein of 76 amino acids and resides on the major plasmid. The second gene is located on the chromosome in a (G+C)-rich DNA fraction and encodes a slightly larger but highly homologous protein consisting of 79 amino acids. The plasmid encoded vac gene is transcribed constitutively throughout the growth cycle while the chromosomal vac gene is expressed during the stationary phase of growth. Comparison of the nucleotide sequences of the two genes indicates differences in the putative promoter regions as well as 35 single base-pair exchanges within the coding regions of the two genes. The majority of the nucleotide exchanges in the coding region occur in the third position of a codon triplet generating the codon synonym. The only differences between the two encoded proteins are the exchange of 2 amino acids (positions 8 and 29) and a deletion of 3 amino acids near the carboxy-terminus of the plasmid encoded vac protein. The genomic DNAs from other halobacterial isolates (Halobacterium sp. SB3, GN101 and YC819-9) were found to contain only a chromosomal vac gene copy. There is a high conservation of the chromosomal vac gene and the genomic region surrounding it among the halobacterial strains investigated.     

Structural position correlation analysis (SPCA) for protein family

Background

The proteins in a family, which perform the similar biological functions, may have very different amino acid composition, but they must share the similar 3D structures, and keep a stable central region. In the conservative structure region similar biological functions are performed by two or three catalytic residues with the collaboration of several functional residues at key positions. Communication signals are conducted in a position network, adjusting the biological functions in the protein family.

Methodology

A computational approach, namely structural position correlation analysis (SPCA), is developed to analyze the correlation relationship between structural segments (or positions). The basic hypothesis of SPCA is that in a protein family the structural conservation is more important than the sequence conservation, and the local structural changes may contain information of biology functional evolution. A standard protein P⁽⁰⁾ is defined in a protein family, which consists of the most-frequent amino acids and takes the average structure of the protein family. The foundational variables of SPCA is the structural position displacements between the standard protein P⁽⁰⁾ and individual proteins P_i of the family. The structural positions are organized as segments, which are the stable units in structural displacements of the protein family. The biological function differences of protein members are determined by the position structural displacements of individual protein P_i to the standard protein P⁽⁰⁾. Correlation analysis is used to analyze the communication network among segments.

Conclusions

The structural position correlation analysis (SPCA) is able to find the correlation relationship among the structural segments (or positions) in a protein family, which cannot be detected by the amino acid sequence and frequency-based methods. The functional communication network among the structural segments (or positions) in protein family, revealed by SPCA approach, well illustrate the distantly allosteric interactions, and contains valuable information for protein engineering study.     

Genes in S. cerevisiae encoding proteins with domains homologous to the mammalian ras proteins

The ras genes, which were first identified by their presence in RNA tumor viruses and which belong to a highly conserved gene family in vertebrates, have two close homologs in yeast, detectable by Southern blotting. We have cloned both genes (RAS1 and RAS2) from plasmid libraries and determined the complete nucleotide sequence of their coding regions. They encode proteins with nearly 90% homology to the first 80 positions of the mammalian ras proteins, and nearly 50% homology to the next 80 amino acids. Yeast RAS1 and RAS2 proteins are more homologous to each other, with about 90% homology for the first 180 positions. After this, at nearly the same position that the mammalian ras proteins begin to diverge from each other, the two yeast ras proteins diverge radically. The yeast ras proteins, like the proteins encoded by the mammalian genes, terminate with the sequence cysAAX, where A is an aliphatic amino acid. Thus the yeast ras proteins have the same overall structure and interrelationship as the family of mammalian ras proteins. The domains of divergence may correspond to functional domains of the ras proteins. Monoclonal antibody directed against mammalian ras proteins immunoprecipitates protein in yeast cells containing high copy numbers of the yeast RAS2 gene.     

The drosophila ras oncogenes: Structure and nucleotide sequence

Three Drosophila genes homologous to the Ha-ras probe were isolated and mapped to positions 85D, 64B, and 62B on chromosome 3. Two of these genes (termed Dras1 and Dras2) were sequenced. In the case of Dras1, which contains multiple introns, a cDNA clone was isolated and sequenced. In the case of Dras2, the nucleotide sequence of the genomic clone was determined. Each gene codes for a protein with a predicted molecular weight of 21.6 kd. Alignment of the amino acid sequence of Dras1 with the vertebrate Ha-ras protein shows that at the amino terminus and central portion (residues 1–121 and 137–164) the two proteins are remarkably similar, and have an overall homology of 75%. The Dras2 gene lacks significant homology to the vertebrate counterpart at the extreme amino terminus and is homologous only between positions 28–120 and 139–161 (overall homology of 50%). This result suggests that the N terminus of p21 forms a distinct regulatory or functional domain. At the carboxy terminus, the major region of variability among the vertebrate ras proteins, the two Drosophila sequences also display considerable variability. However, both appear to be more similar to exon 4B of the Ki-ras gene.     

Structural proteins of Amsacta moorei,Euxoa auxiliaris, and Melanoplus sanguinipes entomopoxviruses

The structural proteins of Amsacta moorei, Euxoa auxiliaris, and Melanoplus sanguinipes entomopoxviruses (EPVs) were separated by electrophoresis on sodium dodecyl sulfate (SDS)-polyacrylamide gels. More than 35 structural proteins were detected in each virus. Based on the distribution and the variation in the molecular weights of the virus structural proteins little homology was detected between the EPVs and vaccinia virus. The molecular weight of Amsacta EPV occlusion body matrix protein (110,000) was determined by SDS-acrylamide gel electrophoresis. The occlusion body matrix protein of Amsacta EPV occluded virus isolated from infected E. acrea larvae was rapidly degraded at pH 10.6 to peptides of approximately 94,000 and 60,000 daltons. After 2 hr incubation at alkaline pH, Amsacta EPV occlusion body protein was degraded to approximately 56,000 daltons. Proteolysis of occlusion body protein was inhibited by SDS. No proteolytic degradation was detected in occlusion body matrix protein isolated from Amsacta EPV infected BTI-EAA cells. Amino acid analysis indicates that entomopoxvirus occlusion body matrix protein consists of approximately 20% acidic amino acids and 9% of the sulfur-containing amino acids cysteine and methionine.     

SpaK/SpaR Two-component System Characterized by a Structure-driven Domain-fusion Method and in Vitro Phosphorylation Studies

Here we introduce a quantitative structure-driven computational domain-fusion method, which we used to predict the structures of proteins believed to be involved in regulation of the subtilin pathway in Bacillus subtilis, and used to predict a protein-protein complex formed by interaction between the proteins. Homology modeling of SpaK and SpaR yielded preliminary structural models based on a best template for SpaK comprising a dimer of a histidine kinase, and for SpaR a response regulator protein. Our LGA code was used to identify multi-domain proteins with structure homology to both modeled structures, yielding a set of domain-fusion templates then used to model a hypothetical SpaK/SpaR complex. The models were used to identify putative functional residues and residues at the protein-protein interface, and bioinformatics was used to compare functionally and structurally relevant residues in corresponding positions among proteins with structural homology to the templates. Models of the complex were evaluated in light of known properties of the functional residues within two-component systems involving His-Asp phosphorelays. Based on this analysis, a phosphotransferase complexed with a beryllofluoride was selected as the optimal template for modeling a SpaK/SpaR complex conformation. In vitro phosphorylation studies performed using wild type and site-directed SpaK mutant proteins validated the predictions derived from application of the structure-driven domain-fusion method: SpaK was phosphorylated in the presence of ³²P-ATP and the phosphate moiety was subsequently transferred to SpaR, supporting the hypothesis that SpaK and SpaR function as sensor and response regulator, respectively, in a two-component signal transduction system, and furthermore suggesting that the structure-driven domain-fusion approach correctly predicted a physical interaction between SpaK and SpaR. Our domain-fusion algorithm leverages quantitative structure information and provides a tool for generation of hypotheses regarding protein function, which can then be tested using empirical methods.     

Uptake,incorporation and calcium-ionophore-stimulated mobilization of arachidonic,eicosapentaenoic and docosahexaenoic acids by leucocytes of the rainbow trout,Salmo gairdneri

Rainbow trout leucocytes contain high levels of neutral lipid (about 70% of total lipid on a wt% basis) consisting of mostly triacylglycerol, free sterols and sterol esters (25%, 15% and 52% of neutral lipid, respectively). The phospholipids, separated by thin-layer chromatography, consisted predominantly of phosphatidylcholine, phosphatidylethanolamine and phosphatidylserine, each present at about 30% of the total phospholipid. Radiolabelling of the leucocytes for 1 h with 1 μCi (approx. 6 μM) [1−¹⁴C]20:4(n−6), [1−¹⁴C]20:5(n−3) or [1−¹⁴C]22:6(n−3) each gave similar uptake values (approx. 1 · 10⁵ cpm/10⁷ leucocytes). The incorporation into total phospholipids was highest for 22:6(n−3) and lowest for 20:4(n−6). A higher percentage of radiolabel from [1−¹⁴C]22:6(n − 3) was found incorporated into phosphatidylcholine and phosphatidylethanolamine as compared to that from [1−¹⁴C]20:4(n − 6) and [1−^14C]20:5(n−3), while the reverse situation was found with phosphatidylinositol and phosphatidylserine. The relative rates of incorporation into the different phospholipid classes for all three fatty acids were in the order phosphatidylinositol > sphingomyelin > diphosphatidylglycerol > phosphatidylcholine > phosphatidylethanolamine > phosphatidylserine. Calcium ionophore-challenge did not significantly alter the pattern of phospholipid radiolabel. Ionophore-challenge released large amounts of radiolabel, much of which was recovered after high-performance liquid chromatographic separation as free fatty acid/monohydroxy fatty acids, although only approx. 0.3% was recovered in leukotriene B₄ and leukotriene B₅ for the [1−¹⁴C]20:4(n−6) and [1−¹⁴C]20:5(n−3) labelled leucocytes, respectively. Other lipoxygenase products were also radiolabelled and tentatively identified as 20-carboxy-LTB₄, 20-hydroxy-LTB₄, 6-trans-LTB₄, 6-trans-12-epi-LTB₄, 6-trans-8-cis-12-epi-LTB₄ and the corresponding LTB₅ structures. No ‘6-series’ leukotrienes were produced from [1−¹⁴C]22:6(n−3), nor was there any evidence for the synthesis of ‘5-series’ leukotrienes via retroconversion of 22:6(n−3) to 20:5(n−3). This latter finding shows that, despite the preponderance of 22:6(n−3) in the membranes of trout leucocytes, this fatty acid is not a substrate for leukotriene generation.     

VITAL NMR: using chemical shift derived secondary structure information for a limited set of amino acids to assess homology model accuracy

Homology modeling is a powerful tool for predicting protein structures, whose success depends on obtaining a reasonable alignment between a given structural template and the protein sequence being analyzed. In order to leverage greater predictive power for proteins with few structural templates, we have developed a method to rank homology models based upon their compliance to secondary structure derived from experimental solid-state NMR (SSNMR) data. Such data is obtainable in a rapid manner by simple SSNMR experiments (e.g., ¹³C–¹³C 2D correlation spectra). To test our homology model scoring procedure for various amino acid labeling schemes, we generated a library of 7,474 homology models for 22 protein targets culled from the TALOS+/SPARTA+ training set of protein structures. Using subsets of amino acids that are plausibly assigned by SSNMR, we discovered that pairs of the residues Val, Ile, Thr, Ala and Leu (VITAL) emulate an ideal dataset where all residues are site specifically assigned. Scoring the models with a predicted VITAL site-specific dataset and calculating secondary structure with the Chemical Shift Index resulted in a Pearson correlation coefficient (−0.75) commensurate to the control (−0.77), where secondary structure was scored site specifically for all amino acids (ALL 20) using STRIDE. This method promises to accelerate structure procurement by SSNMR for proteins with unknown folds through guiding the selection of remotely homologous protein templates and assessing model quality.     

Prediction of operator-binding protein by discriminant analysis

A number of operator-binding proteins contain similar sequence features to Cro and cI repressors of bacteriophage and CAP protein of Escherichia coli, such as conserved amino acids at constant positions. However, these sequence patterns also occur in proteins that are not operator-binding. We use sequence analogy information in conjunction with a pattern recognition algorithm. The functional and structural properties, e.g., distributions of hydrophobicity, hydrophilicity, charged amino acids, electrostatic free energy, and helical structures of protein are also considered. Within the framework of discriminant analysis, we calculate the above variables and search for a better combination of variables. To assess the discriminatory power of these variables, we allocated additional sequences and predict DNA-binding regions of regulatory proteins not included in the training set.     

Improved thermostability of a Bacillus subtilis esterase by domain exchange

A moderately thermostable esterase from Geobacillus stearothermophilus (BsteE) and its homolog from Bacillus subtilis (BsubE) show a high structural similarity with more than 95 % homology and 74 % amino acid identity. Interestingly, their thermal stability differs significantly by 30 °C in their melting temperature. In order to identify the positions that are responsible for this difference, most of the flexible amino acids assumed to confer instability were found to be in the cap region. For this reason, a 30 amino acid long cap domain fragment containing ten differing positions derived from BsteE was incorporated into the homologous gene encoding for the more labile BsubE by spliced overlap-extension PCR. The melting temperature of the two wild-type esterases and the mutant was evaluated by circular dichroism spectroscopy, while the kinetic parameters and the stability were determined with a photometric assay. The cap domain mutant maintained its activity, with a catalytic efficiency more similar to BsteE, while it exhibited an increase of the melting temperature by 4 °C compared to BsubE. Additional point mutations based on the differences of the parent enzymes gave a further increase of the thermostability up to 11 °C compared to BsubE; however, a significant reduction in activity was observed.     

The complete amino acid sequences of the β1- and β2-subunits of the isolectins LoL1 and LoL11 from seeds of Lathyrus ochrus (L.) DC

The complete amino acid sequences of the β₁- and β₂-subunits of the isolectins (LoL1 and LoL11) from seeds of Lathyrus ochrus were determined by analysis of peptides derived from the proteins by digestion with trypsin, chymotrypsin, pepsin and the S. aureus V₈ protease, as well as fragments produced by cleavage with iodosobenzoic acid. Both β-subunits consisted of singlepolypeptide chains of 181 amino acids, which differed from one another in only 3 positions. The homology of the Lathyrus ochrus isolectins with the other two-chain lectins of the tribe Vicieae, and the single-chain lectins of other tribes of the Leguminosae is discussed.     

Changes of fatty acid composition of phopholipids and lipid structural order in rat liver mitochondrial membrane subsequent to galactosamine intoxication. Effect of clofibrate

Since it has been earlier reported that d-galactosamine induces an inhibition of palmitoylcarnitine transferase I and a depletion of mitochondrial phospholipids which were both prevented by clofibrate, an evaluation of the effects of these drugs on mitochondrial fatty acid composition was made. Galactosamine does not alter the fatty acid pattern of these fatty acids whereas clofibrate induces a 2-fold increase in monounsaturated / saturated fatty acids ratio and a 10-fold decrease of the 20:4 (n − 6)/20:3 (n − 6) ratio in phosphatidylcholine. These alterations suggest an increase of Δ⁹-desaturation and decrease of Δ⁵-desaturation. To determine whether the drug-induced changes in mitochondrial phospholipids has an effect on the physical properties of the membrane, the lipid structural order of mitochondrial preparations was studied using the lipophilic probes DPH and TMA-DPH. Mitochondria isolated either from galactosamine- or clofibrate-treated rats showed a decrease in fluorescence polarization, indicating an overall decrease in lipid structural order. This alteration is more drastic when both drugs are administered. This phenomenon suggests drastic changes in the bulk phase of inner mitochondrial membrane lipids after treatments and could explain the altered kinetic properties of palmitoylcarnitine transferase I.     

Purification and characterization of folate binding proteins from rat placenta

Rat placenta contains virtually no unsaturated (i.e., apo-form) folate binding protein. However, by lowering the pH of a solubilized membrane preparation of this tissue to 3.5, the endogenous bound folate was dissociated from the protein and adsorbed to charcoal. The apo-form of the folate binding protein thus obtained was purified by affinity chromatography using pteroylglutamic acid covalently coupled to Sepharose 4B. A single protein band with an apparent M_r of 36 000 was observed by SDS-polyacrylamide gel electrophoresis of the eluate from the affinity matrix. Western blot of this preparation using a rabbit antiserum raised with the affinity eluate also identified a single 36 kDa protein band. However, peptide sequencing of the N-terminal region of the proteins in the affinity eluate established that it contained two homologous proteins. Computer alignment of the first 22 N-terminal amino acids of each rat placental protein with human, bovine milk and mouse folate binding proteins showed 50–64% identical homology and 27% homology when the eight proteins were aligned together. The affinity of both rat proteins is highest for pteroylglutamic acid (K_a = 1.6 − 109 l/mol) lower for N⁵-methyltetrahydrofolate and substantially lower for N⁵-formyltetrahydrofolate. In the dose-response range studied there was no apparent affinity for methotrexate. The folate binding proteins could be released from a preparation of placental membranes using phospholipase C indicating that these proteins belong to the class of proteins anchored to the plasma membrane by a glycosyl phosphatidylinositol adduct.     

Arabidopsis thaliana sequence analysis confirms the presenceof cyt b-561 in plants: Evidence for a novel protein family

Recent advances in the Arabidopsis sequencing project has elucidated the presence of two genes Atb561-A and Atb561-B that show limited homology to the DNA sequence encoding for the mammalian chromaffin granule cytochrome b-561 (cyt b-561). Detailed analysis of the structural features and conserved residues reveals, however, that the structural homology between the presumptive Arabidopsis proteins and the animal proteins is very high. All proteins are hydrophobic and show highly conserved transmembrane helices. The presumably heme-binding histidine residues in the plant and animal sequences as well as the suggested binding site for the electron acceptor, monodehydroascorbate, are strictly conserved. In contrast, the suggested electron donor (ascorbate) binding site is not very well conserved between the plant and animal sequences questioning the function of this motif. Sequence analysis of the Atb561-B gene demonstrates a different splicing than that initially predicted in silico resulting in a protein with nine extra amino acids and a significantly higher homology to the other cyt b-561 sequences. The homology between the plant and animal sequences is further supported by the strong similarity between a number of biochemical properties of the chromaffin cyt b-561 and the cyt b-561 isolated from bean hook plasma membranes. Since the mammalian cyt b-561 is considered specific to neuroadrenergic tissues, the identification of a closely related homologue in an aneural organism demonstrates that these proteins constitute a new class of widely occurring membrane proteins. Both the plant and animal cyt b-561 are involved in transmembrane electron transport using ascorbate as an electron donor. The similarity between these proteins therefore suggests, for the first time, that this transport supports a number of different cell physiological processes. An evolutionary relationship between the plant and animal proteins is presented.     

Cluster analysis identifies aminoacid compositional features that indicate Toxoplasma gondii adhesin proteins

Toxoplasma gondii invade host cells using a multi-step process that depends on the regulated secretion of adhesions. To identify key primary sequence features of adhesins in this parasite, we analyze the relative frequency of individual amino acids, their dipeptide frequencies, and the polarity, polarizability and Van der Waals volume of the individual amino acids by using cluster analysis. This method identified cysteine as a key amino acid in the Toxoplasma adhesin group. The best vector algorithm of non-concatenated features was for 2 attributes: the single amino acid relative frequency and the dipeptide frequency. Polarity, polarizability and Van der Waals volume were not good classificatory attributes. Single amino acid attributes clustered unambiguously 67 apicomplexan hypothetical adhesins. This algorithm was also useful for clustering hypothetical Toxoplasma target host receptors. All of the cluster performances had over 70% sensitivity and 80% specificity. Compositional aminoacid data can be useful for improving machine learning-based prediction software when homology and structural data are not sufficient.     

The rules of variation: Amino acid exchange according to the rotating circular genetic code

General guidelines for the molecular basis of functional variation are presented while focused on the rotating circular genetic code and allowable exchanges that make it resistant to genetic diseases under normal conditions. The rules of variation, bioinformatics aids for preventative medicine, are: (1) same position in the four quadrants for hydrophobic codons, (2) same or contiguous position in two quadrants for synonymous or related codons, and (3) same quadrant for equivalent codons. To preserve protein function, amino acid exchange according to the first rule takes into account the positional homology of essential hydrophobic amino acids with every codon with a central uracil in the four quadrants, the second rule includes codons for identical, acidic, or their amidic amino acids present in two quadrants, and the third rule, the smaller, aromatic, stop codons, and basic amino acids, each in proximity within a 90 degree angle. I also define codifying genes and palindromati, CTCGTGCCGAATTCGGCACGAG.     

The nucleotide sequence of the β-glucosidase gene from Cellvibrio gilvus

A gene encoding β-glucosidase from Cellvibrio gilvus, a cellobiose-producing bacterium, was cloned into Escherichia coli and sequenced. The structural gene consisted of 2565 bp encoding 854 amino acid residues with a characteristic signal peptide. A typical promoter sequence and SD region were located upstream of the initiation ATG codon. A sequence (180 amino acids) having high homology with those of β-glucosidases from several microorganisms was found in the deduced amino acid sequence of C. gilvus β-glucosidase. This sequence contains the aspartic acid residue which was found to be an active site residue in Aspergillus wentii β-glucosidase A3. The β-glucosidase gene of C. gilvus contains a high amount (69.4%) of G+C. These bases are localized not in the 3rd position of the codon, as is usually observed in G+C-rich genes, but rather in the 1st position. This result in a peptide which contains an extremely high amount (48%) of four amino acids (Pro, Ala, Arg, Gly) coded by CCN, GCN, CGN, and GGN.