Structure of a cDNA for the pro alpha 2 chain of human type I procollagen. Comparison with chick cDNA for pro alpha 2(I) identifies structurally conserved features of the protein and the gene

Nucleotide sequences were determined for cloned cDNAs encoding for more than half of the pro alpha 2 chain of type I procollagen from man. Comparisons with previously published data on homologous cDNAs from chick embryos made it possible to examine evolution of the gene in two species which have diverged for 250-300 million years. The amino acid sequence of the alpha-chain domain supported previous indications that there is a strong selective pressure to maintain glycine as every third amino acid and to maintain a prescribed distribution of charged amino acids. However, there is little apparent selective pressure on other amino acids. The amino acid sequence of the C-propeptide domain showed less divergence than the alpha-chain domain. The 5' end or N terminus of the human C-propeptide, however, contained an insert of 12 bases coding for 4 amino acids not found in the chick C-propeptide. About 100 amino acid residues from the N terminus, two residues found in the chick sequence were missing from the human. In the second half of the C-propeptide, there was complete conservation of a 37 amino acid sequence and conservation of 50 out of 51 amino acids in the same region, an observation which suggested that the region serves some special purpose such as directing the association of one pro alpha 2(I) C-propeptide with two pro alpha 1(I) C-propeptides so as to produce the heteropolymeric structure of type I procollagen. In addition, comparison of human and chick DNAs for pro alpha 2(I) revealed three different classes of conservation of nucleotide sequence which have no apparent effect on the structure of the protein: a preference for U on the third base position of codons for glycine, proline, and alanine; a high degree of nucleotide conservation in the 51 amino acid highly conserved region of the C-propeptide; a high degree of nucleotide conservation in the 3'-noncoding region. These three classes of nucleotide conservation may reflect unusual features of collagen genes, such as their high GC content or their highly repetitive coding sequences.     

Structural-based analysis of dihydrofolate reductase evolution

The evolution of dihydrofolate reductase (DHFR) was studied through a comprehensive structural-based analysis. An amino acid sequence alignment was generated from a superposition of experimentally determined X-ray crystal structures of wild-type (wt) DHFR from the Protein Data Bank (PDB). Using this structure-based alignment of DHFR, a metric was generated for the degree of conservation at each alignment site - not only in terms of amino acid residue, but also secondary structure, and residue class. A phylogenetic tree was generated using the alignment that compared favorably with the canonical phylogeny. This structure-based alignment was used to confirm that the degree of conservation of active-site residues in terms of both sequence as well as structure was significantly greater than non-active site residues. These results can be used in helping to understand the likely future evolution of DHFR in response to novel therapies.     

胡杨锌指蛋白基因克隆及其结构分析

锌指蛋白属于核转录因子家族,在原核生物与真核生物基因转录调控中发挥作用。分析了耐盐锌指蛋白Alfin-1基因在苜蓿与拟南芥中的保守性后,设计了一对引物。以胡杨水培叶片为材料,从总RNA中通过RT-PCR分离得到一个锌指蛋白基因,其cDNA长924bp。分析其氨基酸序列表明,存在一个典型的Cys2/His2锌指结构,从第556位开始有一个富含G的启动子结合位点GTGGGG。由于具有相同功能的转录因子在结构和DNA结合区的氨基酸序列上具有保守性,因此,从结构分析上可以推测该基因与Alfin-1在功能上是有一定的相关性。     

In vitro and in silico cloning of Xenopus laevis SOD2 cDNA and its phylogenetic analysis

By using the methodology of both wet and dry biology (i.e., RT-PCR and cycle sequencing, and biocomputational technology, respectively) and the data obtained through the Genome Projects, we have cloned Xenopus laevis SOD2 (MnSOD) cDNA and determined its nucleotide sequence. These data and the deduced protein primary structure were compared with all the other SOD2 nucleotide and amino acid sequences from eukaryotes and prokaryotes, published in public databases. The analysis was performed by using both Clustal W, a well known and widely used program for sequence analysis, and AntiClustAl, a new algorithm recently created and implemented by our group. Our results demonstrate a very high conservation of the enzyme amino acid sequence during evolution, which proves a close structure-function relationship. This is to be expected for very ancient molecules endowed with critical biological functions, performed through a specific structural organization. The nucleotide sequence conservation is less pronounced: this too was foreseeable, due to neutral mutations and to the species-specific codon usage. The data obtained by using AntiClustAl are comparable with those produced with Clustal W, which validates this algorithm as an important new tool for biocomputational analysis. Finally, it is noteworthy that evolutionary trees, drawn by using all the available data on SOD2 nucleotide sequences and amino acid and either Clustal W or AntiClustAl, are comparable to those obtained through phylogenetic analysis based on fossil records.     

Determination of the base recognition positions of zinc fingers from sequence analysis.

The CC/HH zinc finger is a small independently folded DNA recognition motif found in many eukaryotic proteins, which ligates zinc through two cysteine and two histidine ligands. A database of 1340 zinc fingers from 221 proteins has been constructed and a program for analysis of aligned sequences written. This paper describes sequence analysis aimed at determining the amino acid positions that recognize the DNA bases, by comparing two types of sequence variation. Using the idea that long runs of adjacent zinc fingers have arisen from internal gene duplication, the conservation of each position of the finger within the runs was calculated. The conservation of each position of the finger between homologous proteins from different species was also noted. A correlation of the two types of conservation showed clusters of related amino acids. One cluster of three positions was found to be especially variable within long runs, but highly conserved between corresponding fingers of homologous proteins; these positions are predicted to be the base contact positions. They match the amino acid positions that contact the bases in the co-crystal structure determined by Pavletich and Pabo [Science, 240, 809-817 (1991)]. An adjacent cluster of four positions on the plot may also be associated with DNA binding. This analysis shows that the base recognition positions can be identified even in the absence of a known structure for a zinc finger. These results are applicable to zinc fingers where the structure of the complex is unknown, in particular suggesting that the individual finger--DNA interaction seen in the Zif268--DNA structure has been conserved in many zinc finger--DNA interactions.     

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.     

Improvement in protein functional site prediction by distinguishing structural and functional constraints on protein family evolution using computational design

The prediction of functional sites in newly solved protein structures is a challenge for computational structural biology. Most methods for approaching this problem use evolutionary conservation as the primary indicator of the location of functional sites. However, sequence conservation reflects not only evolutionary selection at functional sites to maintain protein function, but also selection throughout the protein to maintain the stability of the folded state. To disentangle sequence conservation due to protein functional constraints from sequence conservation due to protein structural constraints, we use all atom computational protein design methodology to predict sequence profiles expected under solely structural constraints, and to compute the free energy difference between the naturally occurring amino acid and the lowest free energy amino acid at each position. We show that functional sites are more likely than non-functional sites to have computed sequence profiles which differ significantly from the naturally occurring sequence profiles and to have residues with sub-optimal free energies, and that incorporation of these two measures improves sequence based prediction of protein functional sites. The combined sequence and structure based functional site prediction method has been implemented in a publicly available web server.     

The primary structure of the alpha subunit of human elongation factor 1. Structural aspects of guanine-nucleotide-binding sites

The primary structure of the alpha subunit of elongation factor 1 (EF-1 alpha) from human MOLT 4 cells was determined by cDNA sequencing. The data show that the conservation of the amino acid sequence is more than 80% when compared with yeast and Artemia EF-1 alpha. An inventory of amino acid sequences around the guanine-nucleotide-binding site in elongation factor Tu from Escherichia coli and homologous amino acid sequences in G proteins, initiation and elongation factors and proteins from the RAS family shows two regions containing conserved sequence elements. Region I has the sequence apolar-Xaa-Xaa-Xaa-Gly-Xaa-Xaa-Yaa-Xaa-Gly-LYs-Thr(Ser)- -Xaa-Xaa-Xaa-Xaa-X-apolar. Except for RAS proteins, Yaa is always an acidic amino acid residue. Region II is characterized by the invariant sequence apolar-apolar-Xaa-Xaa-Asn-Lys-Xaa-Asp. In order to facilitate sequence comparison we have used a graphic display, which is based on the hydrophilicity values of individual amino acids in a sequence.     

Gene recognition in eukaryotic DNA by comparison of genomic sequences.

MOTIVATION: Sequencing of complete eukaryotic genomes and large syntenic fragments of genomes makes it possible to apply genomic comparison for gene recognition. RESULTS: This paper describes a spliced alignment algorithm that aligns candidate exon chains of two homologous genomic sequence fragments from different species. The algorithm is implemented in Pro-Gen software. Unlike other algorithms, Pro-Gen does not assume conservation of the exon-intron structure. Amino acid sequences obtained by the formal translation of candidate exons are aligned instead of nucleotide sequences, which allows for distant comparisons. The algorithm was tested on a sample of human-mammal (mouse), human-vertebrate (Xenopus ) and human-invertebrate (Drosophila ) gene pairs. Surprisingly, the best results, 97-98% correlation between the actual and predicted genes, were obtained for more distant comparisons, whereas the correlation on the human-mouse sample was only 93%. The latter value increases to 95% if conservation of the exon-intron structure is assumed. This is caused by a large amount of sequence conservation in non-coding regions of the human and mouse genes probably due to regulatory elements. AVAILABILITY: Pro-Gen v. 3.0 is available to academic researchers free of charge at http://www.anchorgen.com/pro_gen/pro_gen.html.     

Sequence and genomic structure of the human adult skeletal muscle sodium channel alpha subunit gene on 17q.

The amino acid sequence of the sodium channel alpha subunit from adult human skeletal muscle has been deduced by cross-species PCR-mediated cloning and sequencing of the cDNA. The protein consists of 1836 amino acid residues. The amino acid sequence shows 93% identity to the alpha subunit from rat adult skeletal muscle and 70% identity to the alpha subunit from other mammalian tissues. A 500 kb YAC clone containing the complete coding sequence and two overlapping lambda clones covering 68% of the cDNA were used to estimate the gene size at 35 kb. The YAC clone proved crucial for gene structure studies as the high conservation between ion channel genes made hybridization studies with total genomic DNA difficult. Our results provide valuable information for the study of periodic paralysis and paramyotonia congenita, two inherited neurological disorders which are caused by point mutations within this gene.     

Human liver fatty acid binding protein. Isolation of a full length cDNA and comparative sequence analyses of orthologous and paralogous proteins

We have determined the primary structure of human liver fatty acid binding protein from an analysis of a full length cDNA. This 127-residue 14,178-Da protein exhibits a high degree of sequence conservation when compared to its orthologous homologue, rat liver fatty acid binding protein. It appears likely that this polypeptide arose from two intragenic duplication events. Using a variety of computational techniques, we were unable to find any evidence of amphipathic alpha helical domains in this protein nor any sequence similarities to apolipoproteins and serum albumins. A family of paralogous proteins was defined, whose members share a remarkable degree of sequence homology with share a remarkable degree of sequence homology with human liver fatty acid binding protein. These include rat intestinal fatty acid binding protein, the cellular the P2 protein of myelin. It appears that the small cytosolic fatty acid binding proteins have evolved structural features necessary for lipid-protein interaction which are different from those present in some familiar and better studied extracellular sequences.     

The ConSurf-HSSP database: the mapping of evolutionary conservation among homologs onto PDB structures

The HSSP (Homology-Derived Secondary Structure of Proteins) database provides multiple sequence alignments (MSAs) for proteins of known three-dimensional (3D) structure in the Protein Data Bank (PDB). The database also contains an estimate of the degree of evolutionary conservation at each amino acid position. This estimate, which is based on the relative entropy, correlates with the functional importance of the position; evolutionarily conserved positions (i.e., positions with limited variability and low entropy) are occasionally important to maintain the 3D structure and biological function(s) of the protein. We recently developed the Rate4Site algorithm for scoring amino acid conservation based on their calculated evolutionary rate. This algorithm takes into account the phylogenetic relationships between the homologs and the stochastic nature of the evolutionary process. Here we present the ConSurf-HSSP database of Rate4Site estimates of the evolutionary rates of the amino acid positions, calculated using HSSP's MSAs. The database provides precalculated evolutionary rates for nearly all of the PDB. These rates are projected, using a color code, onto the protein structure, and can be viewed online using the ConSurf server interface. To exemplify the database, we analyzed in detail the conservation pattern obtained for pyruvate kinase and compared the results with those observed using the relative entropy scores of the HSSP database. It is reassuring to know that the main functional region of the enzyme is detectable using both conservation scores. Interestingly, the ConSurf-HSSP calculations mapped additional functionally important regions, which are moderately conserved and were overlooked by the original HSSP estimate. The ConSurf-HSSP database is available online (http://consurf-hssp.tau.ac.il).     

Accurate Prediction of Protein Catalytic Residues by Side Chain Orientation and Residue Contact Density

Prediction of protein catalytic residues provides useful information for the studies of protein functions. Most of the existing methods combine both structure and sequence information but heavily rely on sequence conservation from multiple sequence alignments. The contribution of structure information is usually less than that of sequence conservation in existing methods. We found a novel structure feature, residue side chain orientation, which is the first structure-based feature that achieves prediction results comparable to that of evolutionary sequence conservation. We developed a structure-based method, Enzyme Catalytic residue SIde-chain Arrangement (EXIA), which is based on residue side chain orientations and backbone flexibility of protein structure. The prediction that uses EXIA outperforms existing structure-based features. The prediction quality of combing EXIA and sequence conservation exceeds that of the state-of-the-art prediction methods. EXIA is designed to predict catalytic residues from single protein structure without needing sequence or structure alignments. It provides invaluable information when there is no sufficient or reliable homology information for target protein. We found that catalytic residues have very special side chain orientation and designed the EXIA method based on the newly discovered feature. It was also found that EXIA performs well for a dataset of enzymes without any bounded ligand in their crystallographic structures.     

Determination of peptide regions on the surface of the eubacterial and archaebacterial ribosome by limited proteolytic digestion.

Limited proteolysis was used in combination with two-dimensional gel electrophoresis, blotting, and amino acid sequence analysis to investigate the surface of intact ribosomal subunits at the peptide and amino acid level. Surface sites of 14 ribosomal proteins from Escherichia coli 50S subunits were determined using proteases with different specificities. To assess the evolutionary conservation of ribosomal topography among eubacteria, large subunits from Bacillus stearothermophilus were also subjected to limited proteolysis. The results obtained indicate a conservation of the three-dimensional ribosomal structure at the peptide level. The data for the eubacterial ribosomes are in full agreement with the model of the 50S protein topography derived from immunological data. Furthermore, peptide surface regions of archaebacterial ribosomes have been investigated. The results presented in this work prove that limited proteolysis can successfully be applied to halophilic and thermophilic ribosomes from archaebacteria.     

Dynamic Behavior of an Intrinsically Unstructured Linker Domain Is Conserved in the Face of Negligible Amino Acid Sequence Conservation

Proteins or regions of proteins that do not form compact globular structures are classified as intrinsically unstructured proteins (IUPs). IUPs are common in nature and have essential molecular functions, but even a limited understanding of the evolution of their dynamic behavior is lacking. The primary objective of this work was to test the evolutionary conservation of dynamic behavior for a particular class of IUPs that form intrinsically unstructured linker domains (IULD) that tether flanking folded domains. This objective was accomplished by measuring the backbone flexibility of several IULD homologues using nuclear magnetic resonance (NMR) spectroscopy. The backbone flexibility of five IULDs, representing three kingdoms, was measured and analyzed. Two IULDs from animals, one IULD from fungi, and two IULDs from plants showed similar levels of backbone flexibility that were consistent with the absence of a compact globular structure. In contrast, the amino acid sequences of the IULDs from these three taxa showed no significant similarity. To investigate how the dynamic behavior of the IULDs could be conserved in the absence of detectable sequence conservation, evolutionary rate studies were performed on a set of nine mammalian IULDs. The results of this analysis showed that many sites in the IULD are evolving neutrally, suggesting that dynamic behavior can be maintained in the absence of natural selection. This work represents the first experimental test of the evolutionary conservation of dynamic behavior and demonstrates that amino acid sequence conservation is not required for the conservation of dynamic behavior and presumably molecular function.     

The complete amino acid sequence of proapolipoprotein A-I of chicken high density lipoproteins

The complete amino acid sequence of proapolipoprotein (proapo) A-I of chicken high density lipoproteins was determined by sequencing overlapping peptides produced by trypsin, S. aureus V8 protease, and cyanogen bromide cleavage. There are 240 amino acid residues in mature chicken apoA-I. By direct sequence analysis of a cyanogen bromide peptide, we also determined the sequence of a 6-amino-acid prosegment which is present at approx. 10% the molar amount of the mature peptide in chicken plasma. Sequence comparison among apoA-I from chicken, human, rabbit, dog and rat, and secondary structure analysis indicate that while the degree of sequence homology is only moderate (less than 50% between chicken and man), there is good conservation of apoA-I secondary structure, especially in the N-terminal two-thirds of the protein in these widely separated species.     

SARS冠状病毒RNA依赖的RNA聚合酶(RdRp)蛋白保守氨基酸基序的鉴定

严重急性呼吸综合片冠状病毒(SARS病毒)的高危害性,使得研究其分子机制并开发有效的治疗药物成为当前生物学家面临的紧迫任务.     

Nucleotide and amino acid sequence of a Cucurbita phytochrome cDNA clone: identification of conserved features by comparison with Avena phytochrome

The amino acid (aa) sequence of Cucurbita phytochrome has been deduced from the nucleotide (nt) sequence of a cDNA clone which was initially identified by hybridization to an Avena phytochrome cDNA clone. Cucurbita, a dicot, and Avena, a monocot, represent evolutionarily divergent groups of plants. The Cucurbita phytochrome polypeptide is 1123 aa in length, corresponding to 125 kDa. Overall, the Cucurbita and Avena phytochrome sequences are 65% homologous at both the nt and aa levels but this sequence conservation is not evenly distributed. Most of the N-terminal two-thirds of the aligned polypeptide chains exhibits localized regions of high conservation, while the extreme N terminus and the C-terminal one-third are less homologous. Comparison of the predicted hydropathic properties of these polypeptides also indicates conservation of domains of phytochrome structure. The possible correlation of these conserved structural features with previously identified functional domains of phytochrome is discussed.     

喜盐鸢尾Na~+/H~+逆转运蛋白基因IhNHX1的克隆及序列分析

采用同源克隆和RACE法克隆获得喜盐鸢尾(Iris halophila Pall.)Na+/H+逆转运蛋白基因IhNHX1的全长序列,该基因序列的全长为1 946 bp,包含1个长度为1 611 bp的开放阅读框(ORF),编码537个氨基酸。序列对比及系统树分析结果表明:IhNHX1基因编码的氨基酸序列与另外11种植物NHX1基因编码的氨基酸序列的一致性高达96.2%,相同序列占61.7%,表明该氨基酸序列保守性较高;在系统树上喜盐鸢尾与其他植物的分支长均大于1.2,表明它们的亲缘关系均较远;IhNHX1基因编码的氨基酸序列含有2个保守结构域,即氨氯吡嗪结合位点和CaM结合结构域,分别是NHX1蛋白的标志性结合位点和重要调节区域。该蛋白质的二级结构和跨膜结构域分析结果表明:在IhNHX1基因编码的蛋白质的二级结构中,α螺旋占48.60%、不规则卷曲占32.03%、延伸链占14.71%、氢键转角占4.66%;该蛋白质含有10个跨膜结构域。此外,对5’RACE方法中5’端正向引物的优化设计步骤进行了归纳,以提高PCR扩增的特异性。