Convergent evolution within the V3 loop domain of human immunodeficiency virus type 1 in association with disease progression.

Phylogenetic analysis was used to study in vivo genetic variation of the V3 region of human immunodeficiency virus type 1 in relation to disease progression in six infants with vertically acquired human immunodeficiency virus type 1 infection. Nucleotide sequences from each infant formed a monophyletic group with similar average branch lengths separating the sets of sequences. In contrast to the star-shaped phylogeny characteristic of interinfant viral evolution, the shape of the phylogeny formed by sequences from the infants who developed AIDS tended to be linear. A computer program, DISTRATE, was written to analyze changes in DNA distance values over time. For the six infants, the rate of divergence from the initial variant was inversely correlated with CD4 cell counts averaged over the first 11 to 15 months of life (r = -0.87, P = 0.024). To uncover evolutionary relationships that might be dictated by protein structure and function, tree-building methods were applied to inferred amino acid sequences. Trees constructed from the full-length protein fragment (92 amino acids) showed that viruses from each infant formed a monophyletic group. Unexpectedly, V3 loop protein sequences (35 amino acids) that were found at later time points from the two infants who developed AIDS clustered together. Furthermore, these sequences uniquely shared amino acids that have been shown to confer a T-cell line tropic phenotype. The evolutionary pattern suggests that viruses from these infants with AIDS acquired similar and possibly more virulent phenotypes.     

Characterisation of alpha-1 giardin: an immunodominant Giardia lamblia annexin with glycosaminoglycan-binding activity

Alpha-1 giardin is an immunodominant protein in the intestinal protozoan parasite Giardia lamblia. The Triage((R)) parasite panel, used to detect copro-antigens in stool from giardiasis patients, reacts with an epitope between amino acids 160 and 200 in alpha-1 giardin. This region of the protein is also highly immunogenic during human infections. Alpha-1 giardin is related to annexins and like many other annexins it was shown to be plasma membrane associated. Immunoelectron and immunofluorescence microscopy revealed that some alpha-1 giardin are displayed on the surface of recently excysted cells. Recombinant alpha-1 giardin displayed a Ca(2+)-dependent binding to glycosaminoglycans (GAGs), in particular heparan sulphate, a common GAG in the intestinal tract. Recombinant alpha-1 giardin bound to thin sections of human small intestine, a binding which could be inhibited by adding increasing concentrations of sulphated sugars. A surface associated trypsin activated Giardia lectin (taglin) has been suggested to be important for G. lamblia attachment. In this study we show that a monoclonal antibody that inhibits taglin recognises alpha-1 and alpha-2 giardin. Thus, alpha-1 giardin is a highly immunoreactive GAG-binding protein, which may play a key role in the parasite-host interaction. Our results further show a conserved function of annexins from lower to higher eukaryotes.     

Evolutionary relationships of the glucokinase from the amitochondriate protist, Trichomonas vaginalis

Two genes coding for Trichomonas vaginalis glucokinase were isolated and sequenced. The putative translation products have molecular masses of 41,584 and 41,772 Da, corresponding to 375 and 377 amino acids, respectively. These values agree with data determined by sodium dodecylsulfate-polyacrylamide gel electrophoresis (SDS-PAGE) for the enzyme purified from the organism. The two sequences showed 78% amino acid identity. The sequences and their phylogenetic reconstruction show that they are members of a glucokinase/fructokinase protein family found in eubacteria and also in the eukaryote Giardia lamblia and are only distantly related to typical eukaryotic hexokinases. The results indicate that the evolutionary past of this enzyme, catalyzing the first step of glycolysis in T. vaginalis, is different from that of the enzyme performing this key role in almost all other eukaryotes.     

Positive and negative design in stability and thermal adaptation of natural proteins

The aim of this work is to elucidate how physical principles of protein design are reflected in natural sequences that evolved in response to the thermal conditions of the environment. Using an exactly solvable lattice model, we design sequences with selected thermal properties. Compositional analysis of designed model sequences and natural proteomes reveals a specific trend in amino acid compositions in response to the requirement of stability at elevated environmental temperature: the increase of fractions of hydrophobic and charged amino acid residues at the expense of polar ones. We show that this “from both ends of the hydrophobicity scale” trend is due to positive (to stabilize the native state) and negative (to destabilize misfolded states) components of protein design. Negative design strengthens specific repulsive non-native interactions that appear in misfolded structures. A pressure to preserve specific repulsive interactions in non-native conformations may result in correlated mutations between amino acids that are far apart in the native state but may be in contact in misfolded conformations. Such correlated mutations are indeed found in TIM barrel and other proteins.     

Variation of enzyme activities and metabolite levels in 24 Arabidopsis accessions growing in carbon-limited conditions

Our understanding of the interaction of carbon (C) metabolism with nitrogen (N) metabolism and growth is based mainly on studies of responses to environmental treatments, and studies of mutants and transformants. Here, we investigate which metabolic parameters vary and which parameters change in a coordinated manner in 24 genetically diverse Arabidopsis (Arabidopsis thaliana) accessions, grown in C-limited conditions. The accessions were grown in short days, moderate light, and high nitrate, and analyzed for rosette biomass, levels of structural components (protein, chlorophyll), total phenols and major metabolic intermediates (sugars, starch, nitrate, amino acids), and the activities of seven representative enzymes from central C and N metabolism. The largest variation was found for plant weight, reducing sugars, starch at the end of the night, and several enzyme activities. High levels of one sugar correlated with high levels of other sugars and starch, and a trend to increased amino acids, slightly lower nitrate, and higher protein. The activities of enzymes at the interface of C and N metabolism correlated with each other, but were unrelated to carbohydrates, amino acid levels, and total protein. Rosette weight was unrelated or showed a weak negative trend to sugar and amino acid contents at the end of the day in most of the accessions, and was negatively correlated with starch at the end of the night. Rosette weight was positively correlated with several enzyme activities. We propose that growth is not related to the absolute levels of starch, sugars, and amino acids; instead, it is related to flux, which is indicated by the enzymatic capacity to use these central resources.     

In silico reconstruction of the amino acid metabolic pathways of Trypanosoma cruzi

Trypanosoma cruzi is the epidemiological agent of Chagas' disease, affecting most of Central and South America, constituting a significant health and socio-economic problem. The parasite has a metabolism largely based on the consumption of amino acids, which participate in a diversity of metabolic pathways, leading to many crucial compounds for the survival of this parasite. Study of its enzymes has the potential to disclose new therapeutic targets and foster the development of new drugs. In this study, we employed computational approaches to reconstruct in silico the amino acid metabolic pathways of T. cruzi, aiming to link genomic information with functional information. For that, protein sequences from 570 EC classes belonging to 25 different pathways in general amino acid metabolism were downloaded from KEGG. A subset of 471 EC classes had at least one sequence deposited. Clustering of the proteins belonging to each EC class was performed using a similarity-based approach implemented in the tool AnEnPi. Reconstruction of the metabolic pathways comprising the amino acid metabolism of T. cruzi was performed by analyzing the output of BLASTP, using as query the dataset of predicted proteins of T. cruzi against all sequences of each individual cluster. This approach allowed us to identify 764 T. cruzi proteins probably involved in the metabolism of amino acids as well as the identification of several putative cases of analogy. Furthermore, we were able to identify several enzymatic activities of T. cruzi that were not previously included in KEGG.     

Primary Structure and Phylogenetic Relationships of Glyceraldehyde-3-Phosphate Dehydrogenase Genes of Free-Living and Parasitic Diplomonad Flagellates

ABSTRACT. Complete nucleotide sequences have been established for two genes (gap1 and gap2) coding for glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC 1.2.1.12) homologs in the diplomonad Giardia lamblia. In addition, almost complete sequences of the GAPDH open reading frames were obtained from PCR products for two free-living diplomonad species, Trepomonas agilis and Hexamita inflata, and a parasite of Atlantic salmon, an as yet unnamed species with morphological affinities to Spironucleus. Giardia lamblia gap 1 and the genes from the three other diplomonad species show high similarity to each other and to other glycolytic GAPDH genes. All amino-acyl residues known to be highly conserved in this enzyme are also conserved in these sequences. Giardia lamblia gap2 gene is more divergent and its putative translation reveals the presence of a cysteine and serine-rich insertion resembling a metal binding finger. This motif has not yet been noted in other GAPDH molecules. All sequences contain an S-loop signature with characteristics close to those of eukaryotes. In phylogenetic reconstructions based on the derived amino acid sequences with neighborjoining, parsimony and maximum-likelihood methods the four typical GAPDH sequences of diplomonads cluster into a single clade. Within this clade, G. lamblia gap1 shares a common ancestor with the rest of the genes. The latter are more closely related to each other, indicating an early separation of the lineage leading to the genus Giardia from the lineage encompassing the morphologically less differentiated genera, Trepomonas, Hexamita and that of the unnamed species. This result is discordant with the orthogonal evolution of diplomonads suggested on the basis of comparative morphology. In neighbor-joining reconstructions G. lamblia gap2 occupies a variable position, due to its great divergence. In parsimony and maximum likelihood analysis however, it shares a most recent common ancestor with the typical G. lamblia gap1 gene, suggesting that it diverged after the separation of the Giardia lineage. The position of the diplomonad clade in broader phylogenetic reconstructions is firmly within the typical cytosolic glycolytic representatives of GAPDH of eukaryotes.     

Functional dissection of the C-terminal domain of type II DNA topoisomerase from the kinetoplastid hemoflagellate Leishmania donovani

The amino acid sequences of the C-terminal domain (CTD) of the type II DNA topoisomerases are divergent and species specific as compared with the highly conserved N-terminal and central domains. A set of C-terminal deletion mutants of Leishmania donovani topoisomerase II was constructed. Removal of more than 178 amino acids out of 1236 amino acid residues from the C-terminus inactivates the enzyme, whereas removal of 118 amino acids or less has no apparent effect on the ability of the parasite enzyme to complement a temperature-sensitive mutation of the Saccharomyces cerevisiae topoisomerase II gene. Deletion analysis revealed a potent nuclear localization signal (NLS) within the amino acid residues 998–1058. Immunomicroscopy results suggest that the removal of an NLS in the CTD is likely to contribute to the physiological dysfunction of these proteins. Modeling of the LdTOP2 based on the crystal structure of the yeast type II DNA topoisomerase showed that the parasite protein assumes a structure similar to its yeast counterpart harboring all the conserved residues in a structurally similar position. However, a marked difference in electrostatic potential was found in a span of 60 amino acid residues (998–1058), which also do not have any homology with topoisomerase II sequences. Such significant differences can be exploited by the structure-based design of selective inhibitors using the structure of the Leishmania enzyme as a template.     

一种新的蛋白质序列的向量表示方法及其应用

为了更多地挖掘隐藏在蛋白质序列中的信息,本研究将20种氨基酸均匀地排列在单位圆周上,得到每种氨基酸对应的二维坐标,再与氨基酸的6个理化指标结合起来,最终用一个八维向量来刻画蛋白质序列。为避免数据极差对分析结果造成的影响,本研究对蛋白质序列所对应的八维向量作归一化处理。基于归一化后的蛋白质序列的向量表示,运用神经网络对蛋白质序列进行分类,并根据向量之间的欧式距离来量化序列之间的相似性。最后,以9个不同物种的ND5蛋白质序列以及8个不同物种的ND6蛋白质序列为例,Clustal W序列比对方法为基准,对本研究的方法与5-字母方法进行验证和比较,结果表明本研的方法是有效的。     

Ultrastructural localization of giardins to the edges of disk microribbons of Giarida lamblia and the nucleotide and deduced protein sequence of alpha giardin

The giardins are a group of 29-38-kD proteins in the ventral disk of the protozoan parasite Giardia lamblia. The disk attaches the parasite to the host's intestinal epithelium and is composed of parallel, coiled microtubules that are adjacent to the ventral plasma membrane and from which processes called microribbons extend into the cytoplasm; the microribbons are connected by crossbridges. G. lamblia cytoskeletons, consisting of disks and attached flagella, were isolated and used to show that the 29-38-kD proteins separate into five bands by one-dimensional electrophoresis and into 23 species by two-dimensional analysis. Rabbit antibodies raised against a 33-kD protein band, purified by one-dimensional gel electrophoresis and shown to contain three proteins by two-dimensional electrophoresis, recognized 17 proteins by two-dimensional immunoblot analysis. By immunofluorescence these antibodies reacted with the ventral disk but not with the flagella in isolated cytoskeletons. Electron microscopy revealed that the anti-giardin antibodies bound to the edges of the microribbons but not to the microtubules, crossbridges, or other, nondisk structures. Antibodies to tubulin reacted with both the disk and flagella in isolated cytoskeletons but bound only to the microtubules in these structures. The amino-terminal sequence of the 33-kD immunogen was determined and used to construct a DNA oligomer, and the oligomer was used to isolate the alpha giardin gene. The gene was used to hybrid select RNA, and the in vitro translation product from this RNA was precipitated by the antibodies against the 33-kD immunogen. The gene sequence was a single open reading frame of 885 nucleotides that predicted a protein of 33.8 kD. The protein sequence is unique, having no significant homology to two other giardin sequences or to any sequences within the Protein Identification Resource. It is predicted to be 82% alpha helical. The downstream sequence of the gene indicates that the sequence AGT-PuAA is located six to nine nucleotides beyond the stop codon in all protein-encoding genes of G. lamblia that have been sequenced and reported to date.     

Looking at structure, stability, and evolution of proteins through the principal eigenvector of contact matrices and hydrophobicity profiles

We review and further develop an analytical model that describes how thermodynamic constraints on the stability of the native state influence protein evolution in a site-specific manner. To this end, we represent both protein sequences and protein structures as vectors: structures are represented by the principal eigenvector (PE) of the protein contact matrix, a quantity that resembles closely the effective connectivity of each site; sequences are represented through the "interactivity" of each amino acid type, using novel parameters that are correlated with hydropathy scales. These interactivity parameters are more strongly correlated than the other hydropathy scales that we examine with: (1) the change upon mutations of the unfolding free energy of proteins with two-states thermodynamics; (2) genomic properties as the genome-size and the genome-wide GC content; (3) the main eigenvectors of the substitution matrices. The evolutionary average of the interactivity vector correlates very strongly with the PE of a protein structure. Using this result, we derive an analytic expression for site-specific distributions of amino acids across protein families in the form of Boltzmann distributions whose "inverse temperature" is a function of the PE component. We show that our predictions are in agreement with site-specific amino acid distributions obtained from the Protein Data Bank, and we determine the mutational model that best fits the observed site-specific amino acid distributions. Interestingly, the optimal model almost minimizes the rate at which deleterious mutations are eliminated by natural selection.     

An immunogenic Mr 23,000 integral membrane protein of Schistosoma mansoni worms that closely resembles a human tumor-associated antigen

A Mr 23,000 Ag of the human trematode parasite, Schistosoma mansoni, has been identified by immunoscreening an adult worm cDNA library with antibody affinity purified on the Mr 23,000 to 25,000 integral membrane protein fraction of the parasite. This Ag is immunogenic in infected humans as well as in rabbits exposed to S. mansoni. The protein sequence of the Ag as deduced from cloned DNA sequences is 218 amino acids long and contains four putative transmembrane regions. Of particular significance, the Ag is strikingly similar, with respect to both amino acid sequence (36% identity) and putative domain structure to ME491, a human stage-specific melanoma-associated Ag.     

Molecular characterisation of a predominant antigenic region of Giardia lamblia variant surface protein H7

During infection, the intestinal protozoan parasite Giardia lamblia undergoes continuous antigenic variation which is determined by diversification of the parasite's major surface antigen, named VSP (variant surface protein). One member from this protein family, VSP H7, is expressed by G. lamblia clone GS/M-83-H7. In the present study, we characterised a highly antigenic portion of VSP H7 which is positioned inside a 130 amino acid C-terminal region of the protein. This region overlaps with a cysteine-rich motif that is rather conserved within the VSP family. Detailed molecular dissection of the antigenic portion monitored a 12 amino acid peptidyl structure which constitutes a non-conformational epitope of VSP H7. In the murine host, this epitope is recognised relatively early (before day 10 p.i.) during infection and stimulates a strong intestinal immunoglobulin A response. At late infective stages (after day 10 p.i.) this immune reaction is progressively complemented by reactions against 'late' antigenic epitopes which are also located inside the 130 amino acid antigenic portion but in closer proximity to the C-terminal end of VSP H7 than the 12 amino acid epitope. Both the high antigenicity and the conserved character suggest that the 12 amino acid epitope is a key factor within the immunological interplay between G. lamblia and the experimental murine host.     

The Uptake and Metabolism of Cysteine by Giardia lamblia Trophozoites

ABSTRACT. The cysteine, cystine, methionine and sulfate uptake and cysteine metabolism of Giardia lamblia was studied. Initial experiments indicated that bathocuproine sulphonate (20 μM) added to Keister's modified TYI-S-33 medium supported the growth of G. lamblia at low L-cysteine concentration. This allowed the use of high specific activity radiolabeled L-cysteine for further studies. The analyses of L-cysteine uptake by G. lamblia indicate the presence of at least two different transport systems. The total cysteine uptake was non saturable, with a capacity of 3.7 pmoles per 10⁶ cells per min per μM of cysteine, and probably represent passive diffusion. However, cysteine transport was partially inhibited by L-methionine, D-cysteine and DL-homocysteine. indicating that another system specific for SH-containing amino acids is also present. Cysteine uptake was markedly decreased in medium without serum. In contrast to cysteine, the uptake of L-methionine and sulfate were carried out by saiurable systems with apparent K_m, of 71 and 72 μM, respectively, but the Vmax of the uptake of sulfate was six orders of magnitude lower than the Vmax of methionine uptake. Cystine was not incorporated into trophozoites. [³⁵S]-labeled L-cysteine and L-methionine, but not [³⁵S]sulfate, were incorporated into Giardia proteins, indicating that the parasite lacks the capacity to synthesize cysteine or methionine from sulfate. Neither cystathionine γ lyase nor crystathionine γ synthase activities was detected in homogenates of Giardia lamblia , suggesting that the transsulfuration pathway is not active and there is no conversion of methionine to cysteine. Our data indicate that cysteine is essential for Giardia because the parasite: a) cannot take up cystine, and b) cannot synthesize cysteine de novo.     

Codon Usage in Giardia lamblia

ABSTRACT A codon usage table for the intestinal parasite Giardia lamblia was generated by analysis of the nucleotide sequences of eight genes comprising 3,135 codons. Codon usage revealed a biased use of synonomous codons with a preference for NNC codons (42.1%). The codon usage of G. lamblia more closely resembles that of the prokaryote Halobacterium halobium (correlation coefficient r = 0.73) rather than that of other eukaryotic protozoans, i.e. Trypanosoma brucei ( r = 0.434) and Plasmodium falciparum ( r =–0.31). These observations are consistent with the view that G. lamblia represents the first line of descent from the ancestral cells that first took on eukaryotic features.     

Codon usage in Giardia lamblia.

A codon usage table for the intestinal parasite Giardia lamblia was generated by analysis of the nucleotide sequences of eight genes comprising 3,135 codons. Codon usage revealed a biased use of synonymous codons with a preference for NNC codons (42.1%). The codon usage of G. lamblia more closely resembles that of the prokaryote Halobacterium halobium (correlation coefficient r = 0.73) rather than that of other eukaryotic protozoans, i.e. Trypanosoma brucei (r = 0.434) and Plasmodium falciparum (r = -0.31). These observations are consistent with the view that G. lamblia represents the first line of descent from the ancestral cells that first took on eukaryotic features.     

Core histones of the amitochondriate protist, Giardia lamblia

Genes coding for the core histones H2a, H2b, H3, and H4 of Giardia lamblia were sequenced. A conserved organism- and gene-specific element, GRGCGCAGATTTVGG, was found upstream of the coding region in all core histone genes. The derived amino acid sequences of all four histones were similar to their homologs in other eukaryotes, although they were among the most divergent members of this protein family. Comparative protein structure modeling combined with energy evaluation of the resulting models indicated that the G. lamblia core histones individually and together can assume the same three-dimensional structures that were established by X-ray crystallography for Xenopus laevis histones and the nucleosome core particle. Since G. lamblia represents one of the earliest-diverging eukaryotes in many different molecular trees, the structure of its histones is potentially of relevance to understanding histone evolution. The G. lamblia proteins do not represent an intermediate stage between archaeal and eukaryotic histones.     

Nucleotide sequences of the mRNA''s encoding the vesicular stomatitis virus G and M proteins determined from cDNA clones containing the complete coding regions.

The complete nucleotide sequences of the vesicular stomatitis virus mRNA's encoding the glycoprotein (G) and the matrix protein (M) have been determined from cDNA clones that contain the complete coding sequences from each mRNA. The G protein mRNA is 1,665 nucleotides long, excluding polyadenylic acid, and encodes a protein of 511 amino acids including a signal peptide of 16 amino acids. G protein contains two large hydrophobic domains, one in the signal peptide and the other in the transmembrane segment near the COOH terminus. Two sites of glycosylation are predicted at amino acid residues 178 and 335. The close correspondence of the positions of these sites with the reported timing of the addition of the two oligosaccharides during synthesis of G suggests that glycosylation occurs as soon as the appropriate asparagine residues traverse the membrane of the rough endoplasmic reticulum. The mRNA encoding the vesicular stomatitis virus M protein is 831 nucleotides long, excluding polyadenylic acid, and encodes a protein of 229 amino acids. The predicted M protein sequence does not contain any long hydrophobic or nonpolar domains that might promote membrane association. The protein is rich in basic amino acids and contains a highly basic amino terminal domain. Details of construction of the nearly full-length cDNA clones are presented.     

Rapid assessment of correlated amino acids from pair-to-pair (P2P) substitution matrices

Identification of correlated amino acids in proteins has been a topic of broad interest in view of its functional implications and importance in protein design. A new set of pair-to-pair (P2P) substitution matrices for amino acids was recently introduced as a useful tool for inferring information on such correlated sites. We present a website developed for automated application of these matrices for analysis of query sequences. The site offers options for graphical analysis of correlations, as well as visualization of correlated amino acids on representative, structurally characterized, members of the examined family of sequences. Availability: http://www.ccbb.pitt.edu/p2p.