Predicting reliable regions in protein alignments from sequence profiles

For applications such as comparative modelling one major issue is the reliability of sequence alignments. Reliable regions in alignments can be predicted using sub-optimal alignments of the same pair of sequences. Here we show that reliable regions in alignments can also be predicted from multiple sequence profile information alone.Alignments were created for a set of remotely related pairs of proteins using five different test methods. Structural alignments were used to assess the quality of the alignments and the aligned positions were scored using information from the observed frequencies of amino acid residues in sequence profiles pre-generated for each template structure. High-scoring regions of these profile-derived alignment scores were a good predictor of reliably aligned regions.These profile-derived alignment scores are easy to obtain and are applicable to any alignment method. They can be used to detect those regions of alignments that are reliably aligned and to help predict the quality of an alignment. For those residues within secondary structure elements, the regions predicted as reliably aligned agreed with the structural alignments for between 92% and 97.4% of the residues. In loop regions just under 92% of the residues predicted to be reliable agreed with the structural alignments. The percentage of residues predicted as reliable ranged from 32.1% for helix residues to 52.8% for strand residues.This information could also be used to help predict conserved binding sites from sequence alignments. Residues in the template that were identified as binding sites, that aligned to an identical amino acid residue and where the sequence alignment agreed with the structural alignment were in highly conserved, high scoring regions over 80% of the time. This suggests that many binding sites that are present in both target and template sequences are in sequence-conserved regions and that there is the possibility of translating reliability to binding site prediction.     

Combining protein sequence,structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence‐structure‐dynamics‐function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence‐conserved residues and build phylogenetic tree. Three‐dimensional structure alignment was also applied to obtain structure‐conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.     

茶树肌动蛋白基因（CsActin1）全长cDNA克隆与生物信息学分析

以茶树（Camellia sinensis）萌动芽为材料,根据茶树萌动芽芽抑制消减杂交文库中分离得到的肌动蛋白（actin）基因的5′-片段设计引物,利用3′-RACE技术克隆了其cDNA全长序列,该基因cDNA全长1 470 bp,命名为CsActin1（GenBank登录号HQ235647）。序列分析表明,CsActin1开放阅读框长1 134 bp,编码377个氨基酸,5′非编码区100 bp,3′非编码区236 bp。推测的蛋白质分子量为41.70 kD,等电点约为5.31,具有肌动蛋白家族的特征信号序列（YVGDEAQs.KRG和WIAKaEYDE）和肌动蛋白相关蛋白的特征信号序列（LLTEApLNPkaNR）。CsActin1与GenBank中注册的其它植物肌动蛋白核苷酸序列的相似性在80%以上,氨基酸序列相似性在95%以上。与其它植物肌动蛋白的进化树分析结果表明,茶树肌动蛋白与杨树的两个肌动蛋白间的亲缘关系最为密切。并对推导的蛋白结构进行了分析。     

Comparison of the three-dimensional structure of two human rhinoviruses (HRV2 and HRV14)

An attempt has been made to build a model of human rhinovirus 2 (HRV2) based on the known human rhinovirus 14 (HRV14) structure. HRV2 was selected because its amino acid sequence is known and because it belongs to the minor rhinovirus receptor class as compared to HRV14, which belongs to the major class. Initial alignment of HRV2 with HRV14 based on the primary sequence and the knowledge of the three-dimensional structure of HRV14 showed that the most probable position of the majority of insertions and deletions occurred in the vicinity of the neutralizing immunogenic sites (NIm). Out of a total of 855 amino acids present in one copy of each of the capsid proteins VP1 through VP4 of HRV14, 411 are different between the two viruses. There are also 6 amino acid residues inserted and 14 residues deleted in HRV2 relative to HRV14. Examination of amino acid interactions showed several cases of conservation of function, e.g., salt bridges or the filling of restricted space. The largest variation amongst the residues lining the canyon, the putative receptor binding site, was in the carboxy-terminal residues of VP1.     

Comparative sequence analysis of acid sensitive/resistance proteins in Escherichia coli and Shigella flexneri

The molecular basis for the survival of bacteria under extreme conditions in which growth is inhibited is a question of great current interest. A preliminary study was carried out to determine residue pattern conservation among the antiporters of enteric bacteria, responsible for extreme acid sensitivity especially in Escherichia coli and Shigella flexneri. Here we found the molecular evidence that proved the relationship between E. coli and S. flexneri. Multiple sequence alignment of the gadC coded acid sensitive antiporter showed many conserved residue patterns at regular intervals at the N-terminal region. It was observed that as the alignment approaches towards the C-terminal, the number of conserved residues decreases, indicating that the N-terminal region of this protein has much active role when compared to the carboxyl terminal. The motif, FHLVFFLLLGG, is well conserved within the entire gadC coded protein at the amino terminal. The motif is also partially conserved among other antiporters (which are not coded by gadC) but involved in acid sensitive/resistance mechanism. Phylogenetic cluster analysis proves the relationship of Escherichia coli and Shigella flexneri. The gadC coded proteins are converged as a clade and diverged from other antiporters belongs to the amino acid-polyamine-organocation (APC) superfamily.     

一个鼻咽癌相关EST的鉴定及其全长cDNA序列分析

鼻咽癌是我国南方及东南亚地区常见的恶性肿瘤之一.通过对鼻咽癌染色体高频率杂合性丢失区域3p21的表达序列标签(expressedsequencetag,EST)进行同源性比较分析,运用逆转录聚合酶链式反应的方法,筛选到一个在41.18%(14/34)的鼻咽癌活检组织及20.0%(1/5)的鼻咽癌细胞系中表达下调的ESTBG772301;并用Northern杂交方法,检测了该EST在多种正常成人组织中的表达状况及其所代表基因的转录本大小.在此基础上,对该EST来源的cDNA克隆(IMAGE:4839190)进行直接测序,获得了一个全长为2377bp的新cDNA序列;经生物信息学分析,发现它与已知基因序列无明显同源性,属于一个新基因,定位于染色体3p21.3,被命名为鼻咽癌表达下调基因(NPCEDRG,GenBank登录号:AF538150).其编码的蛋白质含169个氨基酸,与一个已报道的在进化上相对保守、功能未知的人类蛋白Nicolin1(简称NICN1)N端170个氨基酸残基的序列同源性为97%,但缺少NICN1蛋白C端43个氨基酸残基,可能是nicolin1基因不同剪接本的编码产物.     

来源于瘤胃菌Ruminococcus sp.的D-阿洛酮糖3-差向异构酶的底物结合位点分析

【目的】研究来源于瘤胃菌Ruminococcus sp.的D-阿洛酮糖3-差向异构酶的底物结合机制。【方法】通过同源模拟和同源序列比对,筛选与其底物结合相关的关键位点,进而通过定点突变构建突变体并对其动力学性质进行研究。【结果】筛选得到关键位点Y6和A109,构建了突变体Y6F、Y6I、A109P及A109L。【结论】Y6既与底物结合又与催化能力相关,其-OH只与底物结合相关,芳香环则与催化能力和结合能力均相关;而A109则只是底物结合的位点。该研究结果为D-阿洛酮糖3-差向异构酶的催化机理研究及分子改造提供了借鉴。     

Molecular analysis of the Corynebacterium glutamicum transketolase gene

Transketolase is important in production of the aromatic amino acids in Corynebacterium glutamicum. The complete nucleotide sequence of the C. glutamicum transketolase gene has been identified. The DNA-derived protein sequence is highly similar to the transketolase of Mycobacterium tuberculosis, taxonomically related to C. glutamicum. The alignment of the N-terminus regions between both transketolases showed TTG to be the most probable start codon. Potential ribosomal binding and promoter regions were situated upstream from the TTG. The deduced amino acid sequence consists of 700 residues with a calculated molecular mass of 75 kDa, and contains all amino acid residues involved in cofactor and substrate binding in the well-characterized yeast transketolase sequence.     

Alpha-bungarotoxin and the competing antibody WF6 interact with different amino acids within the same cholinergic subsite

In the nicotinic acetylcholine receptors (AChRs), the sequence segment surrounding two invariant vicinal cysteinyl residues at positions 192 and 193 of the alpha subunit contains important structural component(s) of the binding site for acetylcholine and high molecular weight cholinergic antagonists, like snake alpha-neurotoxins. At least a second sequence region contributes to the formation of the cholinergic site. Studying the binding of alpha-bungarotoxin and three different monoclonal antibodies, able to compete with alpha-neurotoxins and cholinergic ligands, to a panel of synthetic peptides as representative structural elements of the AChR from Torpedo, we recently identified the sequence segments alpha 181-200 and alpha 55-74 as contributing to form the cholinergic site (Conti-Tronconi et al., 1990). As a first attempt to elucidate the structural requirements for ligand binding to the subsite formed by the sequence alpha 181-200, we have now studied the binding of alpha-bungarotoxin and of antibody WF6 to the synthetic peptide alpha 181-200, and to a panel of peptide analogues differing from the parental sequence alpha 181-200 by substitution of a single amino acid residue. CD spectral analysis of the synthetic peptide analogues indicated that they all have comparable structures in solution, and they can therefore be used to analyze the influence of single amino acid residues on ligand binding. Distinct clusters of amino acid residues, discontinuously positioned along the sequence 181-200, seem to serve as attachment points for the two ligands studied, and the residues necessary for binding of alpha-bungarotoxin are different from those crucial for binding of antibody WF6. In particular, residues at positions 188-190 (VYY) and 192-194 (CCP) were necessary for binding of alpha-bungarotoxin, while residues W187, T191, and Y198 and the three residues at positions 193-195 (CPD) were necessary for binding of WF6. Comparison of the CD spectra of the toxin/peptide complexes, and those obtained for the same peptides and alpha-bungarotoxin in solution, indicates that structural changes of the ligand(s) occur upon binding, with a net increase of the beta-structure component. The cholinergic binding site is therefore a complex surface area, formed by discontinuous clusters of amino acid residues from different sequence regions. Such complex structural arrangement is similar to the "discontinuous epitopes" observed by X-ray diffraction studies of antibody/antigen complexes [reviewed in Davies et al. (1988)]. Within this relatively large structure, cholinergic ligands bind with multiple points of attachment, and ligand-specific patterns of the attachment points exist.(ABSTRACT TRUNCATED AT 400 WORDS)     

A family-based approach reveals the function of residues in the nuclear receptor ligand-binding domain

Literature studies, 3D structure data, and a series of sequence analysis techniques were combined to reveal important residues in the structure and function of the ligand-binding domain of nuclear hormone receptors. A structure-based multiple sequence alignment allowed for the seamless combination of data from many different studies on different receptors into one single functional model. It was recently shown that a combined analysis of sequence entropy and variability can divide residues in five classes; (1) the main function or active site, (2) support for the main function, (3) signal transduction, (4) modulator or ligand binding and (5) the rest. Mutation data extracted from the literature and intermolecular contacts observed in nuclear receptor structures were analyzed in view of this classification and showed that the main function or active site residues of the nuclear receptor ligand-binding domain are involved in cofactor recruitment. Furthermore, the sequence entropy-variability analysis identified the presence of signal transduction residues that are located between the ligand, cofactor and dimer sites, suggesting communication between these regulatory binding sites. Experimental and computational results agreed well for most residues for which mutation data and intermolecular contact data were available. This allows us to predict the role of the residues for which no functional data is available yet. This study illustrates the power of family-based approaches towards the analysis of protein function, and it points out the problems and possibilities presented by the massive amounts of data that are becoming available in the "omics era". The results shed light on the nuclear receptor family that is involved in processes ranging from cancer to infertility, and that is one of the more important targets in the pharmaceutical industry.     

Ferric hydroxamate binding protein FhuD from Escherichia coli: mutants in conserved and non-conserved regions

Uptake of iron complexes into the Gram-negative bacterial cell requires highly specific outer membrane receptors and specific ATP-dependent (ATP-Binding-Cassette (ABC)) transport systems located in the inner membrane. The latter type of import system is characterized by a periplasmic binding protein (BP), integral membrane proteins, and membrane-associated ATP-hydrolyzing proteins. In Gram-positive bacteria lacking the periplasmic space, the binding proteins are lipoproteins tethered to the cytoplasmic membrane. To date, there is little structural information about the components of ABC transport systems involved in iron complex transport. The recently determined structure of the Escherichia coli periplasmic ferric siderophore binding protein FhuD is unique for an ABC transport system (Clarke et al. 2000). Unlike other BP's, FhuD has two domains connected by a long -helix. The ligand binds in a shallow pocket between the two domains. In vivo and in vitro analysis of single amino acid mutants of FhuD identified several residues that are important for proper functioning of the protein. In this study, the mutated residues were mapped to the protein structure to define special areas and specific amino acid residues in E. coli FhuD that are vital for correct protein function. A number of these important residues were localized in conserved regions according to a multiple sequence alignment of E. coli FhuD with other BP's that transport siderophores, heme, and vitamin B₁₂. The alignment and structure prediction of these polypeptides indicate that they form a distinct family of periplasmic binding proteins.     

Identification of Amino Acids in the N-Terminal Domain of Corticotropin-Releasing Factor Receptor 1 that Are Important Determinants of High-Affinity Ligand Binding

Abstract : The aim of the present study was to identify the N-terminal regions of human corticotropin-releasing factor (CRF) receptor type 1 (hCRF-R1) that are crucial for ligand binding. Mutant receptors were constructed by replacing specific residues in hCRF-R1 with amino acids from the corresponding position in the N-terminal region of the human vasoactive intestinal peptide receptor type 2 (hVIP-R2). In cyclic AMP stimulation and CRF binding assays, it was established that two regions within the N-terminal domain were crucial for the binding of CRF receptor agonists and antagonists : one region mapping to amino acids 43-50 and a second amino acid sequence extending from position 76 to 84 of hCRF-R1. Recently, it was found that the latter sequence plays a very important role in determining the high ligand selectivity of the Xenopus CRF-R1 (xCRF-R1). Replacement of amino acids 76-84 of hCRF-R1 with residues from the same segment of the hVIP-R2 N terminus markedly reduced the binding affinity of CRF ligands. Mutation of Arg⁷⁶ or Asn⁸¹ but not Gly⁸³ of hCRF-R1 to the corresponding amino acids of xCRF-R1 or hVIP-R2 resulted in 100-1,000-fold lower affinities for human/rat CRF, rat urocortin, and astressin. These data underline the importance of the N-terminal domain of CRF-R1 in high-affinity ligand binding.     

Prediction of residues involved in inhibitor specificity in the dihydrofolate reductase family

Dihydrofolate reductase (DHFR) is of significant recent interest as a target for drugs against parasitic and opportunistic infections. Understanding factors which influence DHFR homolog inhibitor specificity is critical for the design of compounds that selectively target DHFRs from pathogenic organisms over the human homolog. This paper presents a novel approach for predicting residues involved in ligand discrimination in a protein family using DHFR as a model system. In this approach, the relationship between inhibitor specificity and amino acid composition for sets of protein homolog pairs is examined. Similar inhibitor specificity profiles correlate with increased sequence homology at specific alignment positions. Residue positions that exhibit the strongest correlations are predicted as specificity determinants. Correlation analysis requires a quantitative measure of similarity in inhibitor specificity (S(lig)) for a pair of homologs. To this end, a method of calculating S(lig) values using K(I) values for the two homologs against a set of inhibitors as input was developed. Correlation analysis of S(lig) values to amino acid sequence similarity scores - obtained via multiple sequence alignments - was performed for individual residue alignment positions and sets of residues on 13 DHFRs. Eighteen alignment positions were identified with a strong correlation of S(lig) to sequence similarity. Of these, three lie in the active site; four are located proximal to the active site, four are clustered together in the adenosine binding domain and five on the βFβG loop. The validity of the method is supported by agreement between experimental findings and current predictions involving active site residues.     

Characterization and sequencing of cDNA clone encoding the phloem protein PP2 of Cucurbita pepo

Direct N-terminal amino acid sequencing of the phloem protein 2 (PP2) from 3-month old Cucurbita pepo L. (pumpkin), purified by SDS-PAGE and blotted onto PVDF membrane, showed that the protein had a blocked N-terminus. However, after in situ cleavage of the polypeptide in a gel slice by cyanogen bromide, 75 residues of sequence on two cyanogen bromide fragments were determined. An oligonucle-otide probe based on this amino acid sequence was used to screen a cDNA library, constructed from mRNA of 3–5-day old seedling hypocotyls, in ZAP II. A cDNA clone (p11A) predicted an amino acid sequence of 218 residues, in full agreement with the sequences determined for two CNBr fragments of PP2, and suggests that the N-terminus of the protein is a blocked methionine residue which is cleaved off by CNBr. Two additional cDNA clones were sequenced but no heterogeneity in the PP2 sequence was found. The deduced amino acid sequence of C. pepo differs in nine residues from the recently published sequence of Cucurbita maxima (Bostwick et al., Plant Cell 4 (1992) 1539–1548). Southern blot showed that PP2 is encoded by a gene family with a relatively large number of members (estimated as 7–15 per haploid genome).     

Complete amino acid sequence of Bombyx egg-specific protein deduced from cDNA clone

Complementary (c)DNA coding for an insect yolk protein, the egg-specific protein of the silkworm Bombyx mori was cloned and the nucleotide sequence determined. The sequence covers the entire coding region of 1,677 base pairs with 5′ and 3′ noncoding regions (21 and 115 base pairs, respectively). The deduced amino acid sequence of the egg-specific protein consists of 559 amino acid residues. The NH₂-terminal 18 amino acid sequence is enriched in hydrophobic amino acids and assumed to be a signal peptide. A sequence, Asn-X-Thr, a potential N-linked glycosylation site, is found at positions 191 to 193. A serine-rich domain is localized in the region from 63 to 90, in which phosphorylation takes place. Cys His motif in 405 to 415 is analogous to a proposed metal binding sequence. Lys¹³²-Asn¹³³ and Arg²²⁸-Asp²²⁹ are probably the sites cleaved by the egg-specific protein protease that appears during embryogenesis. The derived amino acid sequence has no appreciable homology to other sequenced proteins.     

Molecular Analysis of the Corynebacterium glutamicum Transketolase Gene

Transketolase is important in production of the aromatic amino acids in Corynebacterium glutamicum. The complete nucleotide sequence of the C. glutamicum transketolase gene has been identified. The DNA-derived protein sequence is highly similar to the transketolase of Mycobacterium tuberculosis, taxonomically related to C. glutamicum. The alignment of the N-terminus regions between both transketolases showed TTG to be the most probable start codon. Potential ribosomal binding and promoter regions were situated upstream from the TTG. The deduced amino acid sequence consists of 700 residues with a calculated molecular mass of 75 kDa, and contains all amino acid residues involved in cofactor and substrate binding in the well-characterized yeast transketolase sequence.     

甜瓜STK类R基因同源序列的克隆与分析

以植物丝氨酸/苏氨酸蛋白激酶类( serine-threonine kinase,STK)抗病基因产物催化结构域I和Ⅸ的保守氨基酸序列( FGK/V/L/SVYK/RG,DY/IYSF/YGV/I/M)设计简并引物,对甜瓜(Cucumis melo L.)基因组DNA进行PCR扩增,得到大约500 bp的目的条带,通过重组质粒克隆并经PCR检测后得到12条不同的DNA序列,命名为tg1～tg12,其中tg2、tg5、tg9和tg12(Genbank登录号为JN646853 ～JN646856)可以编码完整的氨基酸序列.Blast分析结果显示:4条序列均具有ATP结合部位、底物结合部位和激酶结构域的活化环(A-loop)等,属于典型的蛋白激酶基因家族,可能是STK类R基因的同源序列片段;4条序列与蓖麻(Ricinus communisL.)的STK同源性均较高.氨基酸序列比对结果显示tg2、tg5、tg9和tg12均具有R基因的9个保守结构域,为STK类候选抗病基因类序列.分子系统树显示tg2、tg5、tg9和tg12与已知的R基因(Pto、Lr10和Lectin)在氨基酸水平上的相似性仅为33.5％ ～53.4％,且4个甜瓜同源序列的氨基酸相似性也较低,表明甜瓜RGAs标记可能具有较高的特异性.     

The unique ligand-binding pocket for the human prostacyclin receptor. Site-directed mutagenesis and molecular modeling

The human prostacyclin receptor is a seven-transmembrane alpha-helical G-protein coupled receptor, which plays important roles in both vascular smooth muscle relaxation as well as prevention of blood coagulation. The position of the native ligand-binding pocket for prostacyclin as well as other derivatives of the 20-carbon eicosanoid, arachidonic acid, has yet to be determined. Through the use of prostanoid receptor sequence alignments, site-directed mutagenesis, and the 2.8-A x-ray crystallographic structure of bovine rhodopsin, we have developed a three-dimensional model of the agonist-binding pocket within the seven-transmembrane (TM) domains of the human prostacyclin receptor. Upon mutation to alanine, 11 of 29 candidate residues within TM domains II, III, IV, V, and VII exhibited a marked decrease in agonist binding. Of this group, four amino acids, Arg-279 (TMVII), Phe-278 (TMVII), Tyr-75 (TMII), and Phe-95 (TMIII), were identified (via receptor amino acid sequence alignment, ligand structural comparison, and computer-assisted homology modeling) as having direct molecular interactions with ligand side-chain constituents. This binding pocket is distinct from that of the biogenic amine receptors and rhodopsin where the native ligands (also composed of a carbon ring and a carbon chain) are accommodated in an opposing direction. These findings should assist in the development of novel and highly specific ligands including selective antagonists for further molecular pharmacogenetic studies of the human prostacyclin receptor.     

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.