Testing statistical significance scores of sequence comparison methods with structure similarity

Background  

In the past years the Smith-Waterman sequence comparison algorithm has gained popularity due to improved implementations and rapidly increasing computing power. However, the quality and sensitivity of a database search is not only determined by the algorithm but also by the statistical significance testing for an alignment. The e-value is the most commonly used statistical validation method for sequence database searching. The CluSTr database and the Protein World database have been created using an alternative statistical significance test: a Z-score based on Monte-Carlo statistics. Several papers have described the superiority of the Z-score as compared to the e-value, using simulated data. We were interested if this could be validated when applied to existing, evolutionary related protein sequences.     

Metal ion reconstitution studies of yeast copper-zinc superoxide dismutase: the "phantom" subunit and the possible role of Lys7p

 Using a corrected molar extinction coefficient for yeast apo copper-zinc superoxide dismutase (CuZnSOD), we have confirmed that the metal binding properties of this protein in vitro differ greatly from those of the bovine and human CuZnSOD enzymes. Thus yeast apo CuZnSOD was found to bind only one Co²⁺ per protein dimer under the conditions in which the bovine and human CuZnSOD apoenzymes readily bind two per dimer. The spectroscopic properties characteristic of the two Cu²⁺ plus two Co²⁺ per dimer or four Cu²⁺ per dimer metal-substituted bovine apo CuZnSOD derivatives were obtained for the yeast apoprotein but by the addition of only half of the appropriate metals, i.e., one Cu²⁺ plus one Co²⁺ per dimer or two Cu²⁺ per dimer. This half-metallated yeast CuZnSOD has been characterized by UV-visible and EPR spectroscopy as well as by native polyacrylamide gel electrophoresis. We conclude that yeast apo CuZnSOD, unlike the bovine and human apoproteins, cannot be reconstituted fully with metal ions under the same conditions. Instead, only one subunit of the homodimer, the "normal" subunit, can be remetalled in a fashion reminiscent of the well-characterized bovine protein. The other "phantom" subunit is not competent to bind metals in this fashion. Furthermore, we have shown that CuZnSOD protein isolated from Saccharomyces cerevisiae that lacks the gene coding for the copper chaperone, Lys7p, contains only one metal ion, Zn²⁺, per protein dimer. The possibility that yeast CuZnSOD can exist in multiple conformational states may represent an increased propensity of the yeast protein to undergo changes that can occur in all CuZnSODs, and may have implications for amyotrophic lateral sclerosis. Received: 8 June 1998 / Accepted: 9 September 1998     

Direct immobilization of gangliosides onto gold-carboxymethyldextran sensor surfaces by hydrophobic interaction: applications to antibody characterization

We describe a novel immobilization technique to investigate interactions between immobilized gangliosides (GD3, GM1, and GM2) and their respective antibodies, antibody fragments, or binding partners using an optical biosensor. Immobilization was performed by direct injection onto a carboxymethyldextran sensor chip and did not require derivatization of the sensor surface or the ganglioside. The ganglioside appeared to bind to the sensor surface by hydrophobic interaction, leaving the carbohydrate epitope available for antibody or, in the case of GM1, cholera toxin binding. The carboxyl group of the dextran chains on the sensor surface did not appear to be involved in the immobilization as evidenced by equivalent levels of immobilization following conversion of the carboxyl groups into acyl amino esters, but rather the dextran layer provided a hydrophilic coverage of the sensor chip which was essential to prevent nonspecific binding. This technique gave better reactivity and specificity for anti- ganglioside monoclonal antibodies (anti-GD3: KM871, KM641, R24; and anti-GM2: KM966) than immobilization by hydrophobic interaction onto a gold sensor surface or photoactivated cross-linking onto carboxymethydextran. This rapid immobilization procedure has facilitated detailed kinetic analysis of ganglioside/antibody interactions, with the surface remaining viable for a large number of cycles (>125). Kinetic constants were determined from the biosensor data using linear regression, nonlinear least squares and equilibrium analysis. The values of kd, ka, and KAobtained by nonlinear analysis (KAKM871 = 1.05, KM641 = 1.66, R24 = 0.14, and KM966 = 0.65 x 10(7) M- 1) were essentially independent of concentration and showed good agreement with data obtained by other analytical methods.      

Uclacyanins,stellacyanins, and plantacyanins are distinct subfamilies of phytocyanins: plant-specific mononuclear blue copper proteins.

The cDNAs encoding plantacyanin from spinach were isolated and characterized. In addition, four new cDNA sequences from Arabidopsis ESTs were identified that encode polypeptides resembling phytocyanins, plant-specific proteins constituting a distinct family of mononuclear blue copper proteins. One of them encodes plantacyanin from Arabidopsis, while three others, designated as uclacyanin 1, 2, and 3, encode protein precursors that are closely related to precursors of stellacyanins and a blue copper protein from pea pods. Comparative analyses with known phytocyanins allow further classification of these proteins into three distinct subfamilies designated as uclacyanins, stellacyanins, and plantacyanins. This specification is based on (1) their spectroscopic properties, (2) their glycosylation state, (3) the domain organization of their precursors, and (4) their copper-binding amino acids. The recombinant copper binding domain of Arabidopsis uclacyanin 1 was expressed, purified, and shown to bind a copper atom in a fashion known as "blue" or type 1. The mutant of cucumber stellacyanin in which the glutamine axial ligand was substituted by a methionine (Q99M) was purified and shown to possess spectroscopic properties similar to uclacyanin 1 rather than to plantacyanins. Its redox potential was determined by cyclic voltammetry to be +420 mV, a value that is significantly higher than that determined for the wild-type protein (+260 mV). The available structural data suggest that stellacyanins (and possibly other phytocyanins) might not be diffusible electron-transfer proteins participating in long-range electron-transfer processes. Conceivably, they are involved in redox reactions occurring during primary defense responses in plants and/or in lignin formation.     

Influence of a natural and a synthetic inhibitor of factor XIIIa on fibrin clot rheology

We investigated the origins of greater clot rigidity associated with FXIIIa-dependent cross-linking. Fibrin clots were examined in which cross-linking was controlled through the use of two inhibitors: a highly specific active-center-directed synthetic inhibitor of FXIIIa, 1,3-dimethyl-4,5-diphenyl-2[2(oxopropyl)thio]imidazolium trifluoromethylsulfonate, and a patient-derived immunoglobulin directed mainly against the thrombin-activated catalytic A subunits of thrombin-activated FXIII. Cross-linked fibrin chains were identified and quantified by one- and two-dimensional gel electrophoresis and immunostaining with antibodies specific for the alpha- and gamma-chains of fibrin. Gamma-dimers, gamma-multimers, alpha(n)-polymers, and alpha(p)gamma(q)-hybrids were detected. The synthetic inhibitor was highly effective in preventing the production of all cross-linked species. In contrast, the autoimmune antibody of the patient caused primarily an inhibition of alpha-chain cross-linking. Clot rigidities (storage moduli, G') were measured with a cone and plate rheometer and correlated with the distributions of the various cross-linked species found in the clots. Our findings indicate that the FXIIIa-induced dimeric cross-linking of gamma-chains by itself is not sufficient to stiffen the fibrin networks. Instead, the augmentation of clot rigidity was more strongly correlated with the formation of gamma-multimers, alpha(n)-polymers, and alpha(p)gamma(q)-hybrid cross-links. A mechanism is proposed to explain how these cross-linked species may enhance clot rigidity.     

Loss of in vitro metal ion binding specificity in mutant copper-zinc superoxide dismutases associated with familial amyotrophic lateral sclerosis

The presence of the copper ion at the active site of human wild type copper-zinc superoxide dismutase (CuZnSOD) is essential to its ability to catalyze the disproportionation of superoxide into dioxygen and hydrogen peroxide. Wild type CuZnSOD and several of the mutants associated with familial amyotrophic lateral sclerosis (FALS) (Ala(4) --> Val, Gly(93) --> Ala, and Leu(38) --> Val) were expressed in Saccharomyces cerevisiae. Purified metal-free (apoproteins) and various remetallated derivatives were analyzed by metal titrations monitored by UV-visible spectroscopy, histidine modification studies using diethylpyrocarbonate, and enzymatic activity measurements using pulse radiolysis. From these studies it was concluded that the FALS mutant CuZnSOD apoproteins, in direct contrast to the human wild type apoprotein, have lost their ability to partition and bind copper and zinc ions in their proper locations in vitro. Similar studies of the wild type and FALS mutant CuZnSOD holoenzymes in the "as isolated" metallation state showed abnormally low copper-to-zinc ratios, although all of the copper acquired was located at the native copper binding sites. Thus, the copper ions are properly directed to their native binding sites in vivo, presumably as a result of the action of the yeast copper chaperone Lys7p (yeast CCS). The loss of metal ion binding specificity of FALS mutant CuZnSODs in vitro may be related to their role in ALS.     

Marsupials and monotremes sort genome treasures from junk

A recent landmark paper demonstrates the unique contribution of marsupials and monotremes to comparative genome analysis, filling an evolutionary gap between the eutherian mammals (including humans) and more distant vertebrate species.     

Individual colour patches as multicomponent signals

Colour patches are complex traits, the components of which may evolve independently through a variety of mechanisms. Although usually treated as simple, two-dimensional characters and classified as either structural or pigmentary, in reality colour patches are complicated, three-dimensional structures that often contain multiple pigment types and structural features. The basic dermal chromatophore unit of fishes, reptiles and amphibians consists of three contiguous cell layers. Xanthophores and erythrophores in the outermost layer contain carotenoid and pteridine pigments that absorb short-wave light; iridophores in the middle layer contain crystalline platelets that reflect light back through the xanthophores; and melanophores in the basal layer contain melanins that absorb light across the spectrum. Changes in any one component of a chromatophore unit can drastically alter the reflectance spectrum produced, and for any given adaptive outcome (e.g. an increase in visibility), there may be multiple biochemical or cellular routes that evolution could take, allowing for divergent responses by different populations or species to similar selection regimes. All of the mechanisms of signal evolution that previously have been applied to single ornaments (including whole colour patches) could potentially be applied to the individual components of colour patches. To reach a complete understanding of colour patch evolution, however, it may be necessary to take an explicitly multi-trait approach. Here, we review multiple trait evolution theory and the basic mechanisms of colour production in fishes, reptiles and amphibians, and use a combination of computer simulations and empirical examples to show how multiple trait evolution theory can be applied to the components of single colour patches. This integrative perspective on animal colouration opens up a host of new questions and hypotheses. We offer specific, testable functional hypotheses for the most common pigmentary (carotenoid, pteridine and melanin) and structural components of vertebrate colour patches.     

Mapping the Naked Neck (NA) and Polydactyly (PO) mutants of the chicken with microsatellite molecular markers

The bulked segregant analysis methodology has been used to map, with microsatellite markers, two morphological mutations in the chicken: polydactyly (PO) and naked neck (NA). These autosomal mutations show partial dominance for NA, and dominance with incomplete penetrance for PO. They were mapped previously to different linkage groups of the classical map, PO to the linkage group IV and NA being linked to the erythrocyte antigen CPPP. An informative family of 70 offspring was produced by mating a sire, heterozygous for each of the mutations, to 7 dams homozygous recessive for each locus. Three DNA pools were prepared, pool PO included 20 chicks exhibiting at least one extra-toe, pool NA included 20 non-polydactyly chicks showing the typical phenotype associated with heterozygosity for the naked neck mutation, and pool NP included 20 chicks exhibiting neither of the mutant phenotypes. Typings were done on an ABI-373 automatic sequencer with 147 microsatellite markers covering most of the genome. An unbalanced distribution of sire marker alleles were detected between pool PO, and pools NA and NP, for two markers of chromosome 2p, MCW0082 and MCW0247. A linkage analysis taking into account the incomplete penetrance of polydactyly (80%) was performed with additional markers of this region and showed that the closest marker to the PO locus was MCW0071 (5 cM, lod score = 9). MCW0071 lies within the engrailed gene EN2 in the chicken. In the mouse, the homologous gene maps on chromosome 5, close to the hemimelic extra-toes mutation Hx. In the case of the NA locus, markers of chromosome 3 were selected because CPPP was mapped on this chromosome. Analysis of individual typings showed a linkage of 5.7 cM (lod score = 13) between the NA locus and ADL0237 in the distal region of chromosome 3q. These results contribute to connecting the former classical map to the molecular genetic map of the chicken, and open the way to the identification of the molecular nature of two developmental mutations of the chicken that are known to occur in many breeds of chickens.