A statistical physics perspective on alignment-independent protein sequence comparison.

Motivation: Within bioinformatics, the textual alignment of amino acid sequences has long dominated the determination of similarity between proteins, with all that implies for shared structure, function, and evolutionary descent. Despite the relative success of modern-day sequence alignment algorithms, so-called alignment-free approaches offer a complementary means of determining and expressing similarity, with potential benefits in certain key applications, such as regression analysis of protein structure-function studies, where alignment-base similarity has performed poorly. Results: Here, we offer a fresh, statistical physics-based perspective focusing on the question of alignment-free comparison, in the process adapting results from “first passage probability distribution” to summarize statistics of ensemble averaged amino acid propensity values. In this paper, we introduce and elaborate this approach.

Sequence analysis of complementary DNAs encoding C3 in the nurse shark ({\it Ginglymostoma cirratum\/})

Mammalian C3 is a pivotal complement protein, encoded for by a single gene. In some vertebrate species multiple C3 isoforms are products of different C3 genes. The goal of this study was to determine whether multiple genes encode for shark C3. A protocol was developed for the isolation of mRNA from shark blood for the isolation of C3 cDNA clones. RT-PCR amplification of mRNA, using sense (GCGEQNM) and antisense (TWLTAYV) primers encoding conserved regions of human C3, yielded 21 clones. The C3-like clones isolated shared 97% similarity with each other and 40% similarity to human C3. RACE-PCR amplification of shark liver RNA, using gene specific primers, yielded products ranging from 1800bp to 3000bp. Deduced amino acid sequence, corresponding to 408bp of the 1800bp fragment, was obtained which showed 51% similarity to human C3. These results suggest that nurse shark C3 might be encoded for by more than one gene. ^

Cytokine analysis in Atlantic bottlenose dolphins: Molecular characterization of IL-4, IL-6, and IL-10

The health status of wild and captive Atlantic Bottlenose dolphins ( Tersiops truncatis) is difficult to ascertain. Mass strandings of these animals have been attributed to pollutants, as well as bacterial infections. Using human Enzyme Linked Immuno-Assays (ELISA) for immunological cytokines, I measured soluble cytokine levels with respect to their health status. In a retrospective analysis of dolphin sera, there was a trend of higher cytokine levels in “sick” animals. I cultured dolphin lymphocytes in the presence of a mitogen (PHA), a super antigen (Staph-A), Lipopolysaccharide (LPS), and a calcium flux inducer (PMA). Levels of messenger RNA, from these cultured cells, were assayed with Polymerase Chain Reaction (PCR) using primers for the human cytokines IL-2, IL-4, IL-6, IL-10, Tumor Necrosis Factor, and Interferon gamma. Only IL-4, IL-6, and IL-10 messages were obtained, inferring similar nucleotide homology to the human primer sequences. The PCR products were sequenced. Sixteen IL-4 sequences, twelve IL-6 sequences and seven IL-10 sequences were obtained and analyzed. Each cytokine exhibited the same nucleotide sequence in all dolphins examined. There was no difference in the cytokine profile in response to the various stimuli. The derived amino acid composition for each of the dolphin cytokines was used for molecular modeling, which showed that dolphin IL-4, IL-6, and IL-10 were structurally similar to the corresponding proteins of Perissodactyla. ^

The isolation and characterization of tescalcin, a novel gene differentially expressed in the mouse testis during the early stages of gonadal differentiation

Gonadal development is an ideal model to study organogenesis because a variety of developmental processes can be studied during the differentiation of the bipotential primordium into testis or ovary. To better understand this process, Representational Difference Analysis of cDNA was used to identify genes that are differentially expressed in mouse gonads at 13.5 days post-coitus. The analysis led to the identification of three testis specific genes and a sequence that was only expressed in the ovary. The male genes identified: renin, Col9a3, and a novel gene termed tescalcin had patterns of expression that suggested a role in testis determination. ^ Studies of the tescalcin gene revealed that it is organized into eight exons and seven introns. The gene was located at 64 cM in mouse chromosome 5, where it spans approximately 35 Kb. Three mRNA variants resulting from alternative splicing of intron 5 were identified in mouse tissues. Gel mobility shift assays demonstrated that Sp1 and Sp3 from Y-1, msc-1, and MIN-6 cells nuclear extracts bind the GC-boxes within the tescalcin proximal promoter. Bisulfite sequencing analysis of tescalcin CpG island revealed that it is differentially methylated in male and female mouse embryonic gonads, and that hypermethylation of this region represses expression of tescalcin in the β-TC3 cell line. ^ The major tescalcin mRNA encodes a protein with 214 amino acids that contains a consensus EF-hand Ca2+-binding domain and an N-myristoylation motif. The amino acid sequence of tescalcin is highly conserved among various species, and it showed the highest homology with calcineurin B homologous proteins 1 and 2, and calcineurin B. Western blot analysis using antibodies generated against the tescalcin protein confirmed its presence in specific mouse tissues and cell lines. Immunohistochemical analysis of mouse embryos confirmed the pattern of expression of tescalcin mRNA in fetal testis. Using pull-down assays, glyceraidehydes-3-phosphate dehydrogenase was identified as an interacting and potential functional partner of tescalcin. ^ The identification and characterization of tescalcin as a novel embryonic testicular marker will contribute to the elucidation of the genetic pathways involved in testis development and likely to the understanding of pathological conditions such as sex reversal and infertility. ^

(Table 2) Concentrations and D/L values for three amino acids in bulk foraminifera from ODP Sites 182-1126 and 182-1127

A number of studies have shown that methanogens are active in the presence of sulfate under some conditions. This phenomenon is especially exemplified in carbonate sediments of the southern Australian continental margin. Three sites cored during Ocean Drilling Program (ODP) Leg 182 in the Great Australian Bight have high concentrations of microbially-generated methane and hydrogen sulfide throughout almost 500 m of sediments. In these cores, the sulfate-reducing and methanogenic zones overlap completely; that is, the usual sulfate-methane transition zone is absent. Amino acid racemization data show that the gassy sediments consist of younger carbonates than the low-gas sites. High concentrations of the reduced gases also occur in two ODP sites on the margin of the Bahamas platform, both of which have similar sedimentary conditions to those of the high-gas sites of Leg 182. Co-generation of these reduced gases results from an unusual combination of conditions, including: (1) a thick Quaternary sequence of iron-poor carbonate sediments, (2) a sub-seafloor brine, and (3) moderate amounts of organic carbon. The probable explanation for the co-generation of hydrogen sulfide and methane in all these sites, as well as in other reported environments, is that methanogens are utilizing non-competitive substrates to produce methane within the sulfate-reducing zone. Taken together, these results form the basis of a new model for sulfate reduction and methanogenesis in marine sediments. The biogeochemical end-members of the model are: (1) minimal sulfate reduction, (2) complete sulfate reduction followed by methanogenesis, and (3) overlapping sulfate reduction and methanogenesis with no transition zone.

Granulometry, geochemistry, amino acids, radiocarbon ages and paleomagnetics of samples from Heidemann Valley

A Pliocene (2.6-3.5 Ma) age is determined from glacial sediments studied in a 20m long, 4 m deep trench excavated in Heidemann Valley, Vestfold Hills, East Antarctica. The age determination is based on a combined study of amino acid racemization, diatoms, foraminifera, and magnetic polarity, and supports earlier estimates of the age of the sedimentary section; all are beyond 14C range. Four till units are recognized and documented, and 16 subunits are identified. All are ascribed to deposition during a Late Pliocene glaciation that was probably the last time the entire Vestfold Hills was covered by an enlarged East Antarctic Ice Sheet (EAIS). Evidence for other more recent glacial events of the 'Vestfold Glaciation' may have been due to lateral expansion of the Sorsdal Glacier and limited expansion of the icesheet margin during the Last Glacial Maximum rather than a major expansion of the EAIS. The deposit appears to correlate with a marine deposition event recorded in Ocean Drilling Program Site 1166 in Prydz Bay, possibly with the Bardin Bluffs Formation of the Prince Charles Mountains and with part of the time represented in the ANDRILL AND-1B core in the Ross Sea.

The temperate Burkholderia phage AP3 of the Peduovirinae shows efficient antimicrobial activity against B. cenocepacia of the IIIA lineage

Burkholderia phage AP3 (vB_BceM_AP3) is a temperate virus of the Myoviridae and the Peduovirinae subfamily (P2likevirus genus). This phage specifically infects multidrug-resistant clinical Burkholderia cenocepacia lineage IIIA strains commonly isolated from cystic fibrosis patients. AP3 exhibits high pairwise nucleotide identity (61.7%) to Burkholderia phage KS5, specific to the same B. cenocepacia host, and has 46.7% - 49.5% identity to phages infecting other species of Burkholderia. The lysis cassette of these related phages has a similar organization (putative antiholin, putative holin, endolysin and spanins) and shows 29-98% homology between specific lysis genes, in contrast to Enterobacteria phage P2, the hallmark phage of this genus. The AP3 and KS5 lysis genes have conserved locations and high amino acid sequence similarity. The AP3 bacteriophage particles remain infective up to 5 h at pH 4-10 and are stable at 60°C for 30 min, but are sensitive to chloroform, with no remaining infective particles after 24 h of treatment. AP3 lysogeny can occur by stable genomic integration and by pseudo-lysogeny. The lysogenic bacterial mutants did not exhibit any significant changes in virulence compared to wild-type host strain when tested in the Galleria mellonella moth wax model. Moreover, AP3 treatment of larvae infected with B. cenocepacia revealed a significant increase (P < 0.0001) in larvae survival in comparison to AP3-untreated infected larvae. AP3 showed robust lytic activity, as evidenced by its broad host range, the absence of increased virulence in lysogenic isolates, the lack of bacterial gene disruption conditioned by bacterial tRNA downstream integration site, and the absence of detected toxin sequences. These data suggest the AP3 phage is a promising potent agent against bacteria belonging to most common B. cenocepacia IIIA lineage strains.

Molecular Characterization and Localization of the NAD(P)H Oxidase Components gp91-phox and p22-phox in Endothelial Cells

The production of reactive oxygen species (ROS) within endothelial cells may have several effects, including alterations in the activity of paracrine factors, gene expression, apoptosis, and cellular injury. Recent studies indicate that a phagocyte-type NAD(P)H oxidase is a major source of endothelial ROS. In contrast to the high-output phagocytic oxidase, the endothelial enzyme has much lower biochemical activity and a different substrate specificity (NADH.NADPH). In the present study, we (1) cloned and characterized the cDNA and predicted amino acid structures of the 2 major subunits of rat coronary microvascular endothelial cell NAD(P)H oxidase, gp91-phox and p22-phox; (2) undertook a detailed comparison with phagocytic NADPH oxidase sequences; and (3) studied the subcellular location of these subunits in endothelial cells. Although these studies revealed an overall high degree of homology (.90%) between the endothelial and phagocytic oxidase subunits, the endothelial gp91-phox sequence has potentially important differences in a putative NADPH-binding domain and in putative glycosylation sites. In addition, the subcellular location of the endothelial gp91-phox and p22-phox subunits is significantly different from that reported for the neutrophil oxidase, in that they are predominantly intracellular and collocated in the vicinity of the endoplasmic reticulum. This first detailed characterization of gp91-phox and p22-phox structure and location in endothelial cells provides new data that may account, in part, for the differences in function between the phagocytic and endothelial NAD(P)H oxidases.

An intimate link between antimicrobial peptide sequence diversity and binding to essential components of bacterial membranes

Antimicrobial peptides and proteins (AMPs) are widespread in the living kingdom. They are key effectors of defense reactions and mediators of competitions between organisms. They are often cationic and amphiphilic, which favors their interactions with the anionic membranes of microorganisms. Several AMP families do not directly alter membrane integrity but rather target conserved components of the bacterial membranes in a process that provides them with potent and specific antimicrobial activities. Thus, lipopolysaccharides (LPS), lipoteichoic acids (LTA) or the peptidoglycan precursor Lipid II are targeted by a broad series of AMPs. Studying the functional diversity of immune effectors tells us about the essential residues involved in AMP mechanism of action. Marine invertebrates have been found to produce a remarkable diversity of AMPs. Molluscan defensins and crustacean anti-LPS factors (ALF) are diverse in terms of amino acid sequence and show contrasted phenotypes in terms of antimicrobial activity. Their activity is directed essentially against Gram-positive or Gram-negative bacteria due their specific interactions with Lipid II or Lipid A, respectively. Through those interesting examples, we discuss here how sequence diversity generated throughout evolution informs us on residues required for essential molecular interaction at the bacterial membranes and subsequent antibacterial activity. Through the analysis of molecular variants having lost antibacterial activity or shaped novel functions, we also discuss the molecular bases of functional divergence in AMPs.

Glutathione Transferases as Detoxification Agents: Structural and Functional Studies

Glutathione transferases (GSTs) are a diverse family of enzymes that catalyze the glutathione-dependent detoxification of toxic compounds. GSTs are responsible for the conjugation of the tripeptide glutathione (GSH) to a wide range of electrophilic substrates. These include industrial pollutants, drugs, genotoxic carcinogen metabolites, antibiotics, insecticides and herbicides. In light of applications in biomedicine and biotechnology as cellular detoxification agents, detailed structural and functional studies of GSTs are required. Plant tau class GSTs play crucial catalytic and non-catalytic roles in cellular xenobiotic detoxification process in agronomically important crops. The abundant existence of GSTs in Glycine max and their ability to provide resistance to abiotic and biotic stresses such as herbicide tolerance is of great interest in agriculture because they provide effective and suitable tools for selective weed control. Structural and catalytic studies on tau class GST isoenzymes from Glycine max (GmGSTU10-10, GmGSTU chimeric clone 14 (Sh14), and GmGSTU2-2) were performed. Crystal structures of GmGSTU10-10 in complex with glutathione sulfenic acid (GSOH) and Sh14 in complex with S-(p-nitrobenzyl)-glutathione (Nb-GSH) were determined by molecular replacement at 1.6 Å and 1.75 Å, respectively. Major structural variations that affect substrate recognition and catalytic mechanism were revealed in the upper part of helix H4 and helix H9 of GmGSTU10-10. Structural analysis of Sh14 showed that the Trp114Cys point mutation is responsible for the enhanced catalytic activity of the enzyme. Furthermore, two salt bridges that trigger an allosteric effect between the H-sites were identified at the dimer interface between Glu66 and Lys104. The 3D structure of GmGSTU2-2 was predicted using homology modeling. Structural and phylogenetic analysis suggested GmGSTU2-2 shares residues that are crucial for the catalytic activity of other tau class GSTs–Phe10, Trp11, Ser13, Arg20, Tyr30, Leu37, Lys40, Lys53, Ile54, Glu66 and Ser67. This indicates that the catalytic and ligand binding site in GmGSTU2-2 are well-conserved. Nevertheless, at the ligandin binding site a significant variation was observed. Tyr32 is replaced by Ser32 in GmGSTU2-2 and thismay affect the ligand recognition and binding properties of GmGSTU2-2. Moreover, docking studies revealed important amino acid residues in the hydrophobic binding site that can affect the substrate specificity of the enzyme. Phe10, Pro12, Phe15, Leu37, Phe107, Trp114, Trp163, Phe208, Ile212, and Phe216 could form the hydrophobic ligand binding site and bind fluorodifen. Additionally, side chains of Arg111 and Lys215 could stabilize the binding through hydrogen bonds with the –NO2 groups of fluorodifen. GST gene family from the pathogenic soil bacterium Agrobacterium tumefaciens C58 was characterized and eight GST-like proteins in A. tumefaciens (AtuGSTs) were identified. Phylogenetic analysis revealed that four members of AtuGSTs belong to a previously recognized bacterial beta GST class and one member to theta class. Nevertheless, three AtuGSTs do not belong to any previously known GST classes. The 3D structures of AtuGSTs were predicted using homology modeling. Comparative structural and sequence analysis of the AtuGSTs showed local sequence and structural characteristics between different GST isoenzymes and classes. Interactions at the G-site are conserved, however, significant variations were seen at the active site and the H5b helix at the C-terminal domain. H5b contributes to the formation of the hydrophobic ligand binding site and is responsible for recognition of the electrophilic moiety of the xenobiotic. It is noted that the position of H5b varies among models, thus providing different specificities. Moreover, AtuGSTs appear to form functional dimers through diverse modes. AtuGST1, AtuGST3, AtuGST4 and AtuGST8 use hydrophobic ‘lock–and–key’-like motifs whereas the dimer interface of AtuGST2, AtuGST5, AtuGST6 and AtuGST7 is dominated by polar interactions. These results suggested that AtuGSTs could be involved in a broad range of biological functions including stress tolerance and detoxification of toxic compounds.

Structure-function relationship of new crotamine isoform from the Crotalus durissus cascavella

In this work we isolated a novel crotamine like protein from the Crotalus durissus cascavella venom by combination of molecular exclusion and analytical reverse phase HPLC. Its primary structure was:YKRCHKKGGHCFPKEKICLPPSSDLGKMDCRWKRK-CCKKGS GK. This protein showed a molecular mass of 4892.89 da that was determined by Matrix Assisted Laser Desorption Ionization Time-of-flight (MALDI-TOF) mass spectrometry. The approximately pI value of this protein was determined in 9.9 by two-dimensional electrophoresis. This crotamine-like protein isolated here and that named as Cro 2 produced skeletal muscle spasm and spastic paralysis in mice similarly to other crotamines like proteins. Cro 2 did not modify the insulin secretion at low glucose concentration (2.8 and 5.6 mM), but at high glucose concentration (16.7 mM) we observed an insulin secretion increasing of 2.7-3.0-fold than to control. The Na+ channel antagonist tetrodoxin (6 mM) decreased glucose and Cro 2-induced insulin secretion. These results suggested that Na+ channel are involved in the insulin secretion. In this article, we also purified some peptide fragment from the treatment of reduced and carboxymethylated Cro 2 (RC-Cro 2) with cyanogen bromide and protease V8 from Staphylococcus aureus. The isolated pancreatic beta-cells were then treated with peptides only at high glucose concentration (16.7 mM), in this condition only two peptides induced insulin secretion. The amino acid sequence homology analysis of the whole crotamine as well as the biologically-active peptide allowed determining the consensus region of the biologically-active crotamine responsible for insulin secretion was KGGHCFPKE and DCRWKWKCCKKGSG.

Characterization of a new platelet aggregating factor from crotoxin Crotalus durissus cascavella venom

In this article we investigated the platelet aggregating activity of whole crotoxin and its subunits isolated from Crotalus durissus cascavella venom. During the purification protocols of the venom, using HPLC molecular exclusion, we detected the presence of two different serine protease activities in the gyroxin fraction, and another in the crotoxin fraction, which induced strong and irreversible platelet aggregation, in addition to blood coagulation. From crotoxin, we isolated PLA(2), crotapotin (both fractions corresponding approximately 85% of whole crotoxin) and another minor fraction (F20) that exhibited serine protease activity. After a new fractionation on reverse phase HPLC chromatography, we obtained three other fractions named as F201, F202 and F203. F202 was obtained with high degree of molecular homogeneity with molecular mass of approximately 28 kDa and a high content of acidic amino residues, such as aspartic acid and glutamic acid. Other important amino acids were histidine, cysteine and lysine. This protein exhibited a high specificity for BApNA, a Michaelis-Menten behavior with Vmax estimated in 5.64 mu M/min and a Km value of 0.58 mM for this substrate. In this work, we investigated the ability of F202 to degrade fibrinogen and observed alpha and beta chain cleavage. Enzymatic as well as the platelet aggregation activities were strongly inhibited when incubated with TLCK and PMSF, specific inhibitors of serine protease. Also, F202 induced platelet aggregation in washed and platelet-rich plasma, and in both cases, TLCK inhibited its activity. The N-terminal amino acid sequence of F202 presented a high amino acid sequence homology with other thrombin-like proteins, but it was significantly different from gyroxin. These results showed that crotoxin is a highly heterogeneous protein composed of PLA(2), thrombin-like and other fractions that might explain the diversity of physiological and pharmacological activities of this protein.

Biological and structural characterization of a new PLA(2) from the Crotalus durissus collilineatus venom

In the present article we report on the biological characterization and amino acid sequence of a new basic Phospholipases A(2) (PLA(2)) isolated from the Crotalus durissus collilineatus venom (Cdcolli F6), which showed the presence of 122 amino acid residues with a pI value of 8.3, molecular mass of 14 kDa and revealed an amino acid sequence identity of 80% with crotalic PLA(2)s such as Mojave B, Cdt F15, and CROATOX. This homology, however, dropped to 50% if compared to other sources of PLA(2)s such as from the Bothrops snake venom. Also, this PLA(2) induced myonecrosis, although this effect was lower than that of BthTx-I or whole crotoxin and it was able to induce a strong blockage effect on the chick biventer neuromuscular preparation, independently of the presence of the acid subunid (crotapotin). The neurotoxic effect was strongly reduced by pre-incubation with heparin or with anhydrous acetic acid and rho-BPB showed a similar reduction. The rho-BPB did not reduce significantly the myotoxic activity induced by the PLA(2), but the anhydrous acetic acid treatment and the pre-incu-bation of PLA(2) with heparin reduced significantly its effects. This protein showed a strong antimicrobial activity against Xanthomonas axonopodis passiflorae (Gram-negative), which was drastically reduced by incubation of this PLA(2) with rho-BPB, but this effect was marginally reduced after treatment with anhydrous acetic acid. Our findings here allow to speculate that basic amino acid residues on the C-terminal and molecular regions near catalytic site regions such as Calcium binding loop or rho-wing region may be involved in the binding of this PLA(2) to the molecular receptor to induce the neurotoxic effect. The bactericidal effect, however, was completely dependent on the enzymatic activity of this protein.

Purification and characterization of the biological effects of phospholipase A(2) from sea anemone Bunodosoma caissarum

Sea anemones contain a variety of biologically active substances. Bunodosoma caissarum is a sea anemone from the Cnidaria phylum, found only in Brazilian coastal waters. The aim of the present work was to study the biological effects of PLA(2) isolated from the sea anemone B. caissarum on the isolated perfused kidney, the arteriolar mesenteric bed and on insulin secretion. Specimens of B. caissarum were collected from the Sao Vicente Channel on the southern coast of the State of São Paulo, Brazil. Reverse phase HPLC analysis of the crude extract of B. caissarum detected three PLA(2) proteins (named BcPLA(2)1, BCPLA(2)2 and BcPLA(2)3) found to be active in B. caissarum extracts. MALDI-TOF mass spectrometry of BcPLA(2)1 showed one main peak at 14.7 kDa. The N-terminal amino acid sequence of BcPLA(2)1 showed high amino acid sequence identity with PLA(2) group III protein isolated from the Mexican lizard (PA23 HELSU, HELSU, PA22 HELSU) and with the honey bee Apis mellifera (PLA(2) and 1POC_A). In addition, BcPLA(2)1 also showed significant overall homology to bee PLA(2). The enzymatic activity induced by native BCPLA(2)1 (20 mu g/well) was reduced by chemical treatment with p-bromophenacyl bromide (p-BPB) and with morin. BcPLA(2)1 strongly induced insulin secretion in presence of high glucose concentration. In isolated kidney, the PLA(2) from B. caissarum increased the perfusion pressure, renal vascular resistance, urinary flow, glomerular filtration rate, and sodium, potassium and chloride levels of excretion. BcPLA(2)1, however, did not increase the perfusion pressure on the mesenteric vascular bed. In conclusion, PLA(2), a group III phospholipase isolated from the sea anemone B. caissarum, exerted effects on renal function and induced insulin secretion in conditions of high glucose concentration. (C) 2009 Elsevier Ltd. All rights reserved.

Molecular biology of prolactin and prolactin receptor in a marine teleost: the sea bream(Sparus aurata)

Tese de dout. em Biologia, especialidade de Biologia Molecular, Unidade de Ciências e Tecnologias dos Recursos Aquáticos, Univ. do Algarve