Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Functional Consequences of the Subdomain Organization of the Sulfs

Sulf-1 and Sulf-2 are novel extracellular sulfatases that act on internal glucosamine 6-O-sulfate modifications within heparan sulfate proteoglycans and regulate their interactions with various signaling molecules, including Wnt ligands. Although the Sulfs are multidomain proteins, there is limited information available about how the subdomains contribute to their enzymatic and signaling activities. In this study, we found that both human Sulfs were synthesized as prepro-enzymes and cleaved by a furin-type proteinase to form disulfide-bond linked heterodimers of 75- and 50-kDa subunits. The mature Sulfs were secreted into conditioned medium, as well as retained on the cell membrane. Although the catalytic center resides in the N-terminal 75-kDa subunit, the C-terminal 50-kDa subunit was indispensable for both arylsufatase and glucosamine 6-O-sulfate-endosulfatase activity. We found that the hydrophilic regions of the Sulfs were essential for endosulfatase activity but not for arylsulfatase activity. Using Edman sequencing, we identified furin-type proteinase cleavage sites in Sulf-1 and Sulf-2. Deletion of these sequences resulted in uncleavable forms of Sulfs. The uncleavable Sulfs retained enzymatic activity. However, they were unable to potentiate Wnt signaling, which may be due to their defective localization into lipid rafts on the plasma membrane.Heparan sulfate proteoglycans (HSPGs)² are major components of the extracellular matrix/cell surface and regulate a variety of biological phenomena, including cell proliferation, cell migration, and differentiation (1). These effects are mediated through the ability of HSPGs to bind to a diverse repertoire of protein ligands. Among these are morphogens, growth factors, chemokines, and other classes of molecules (2, 3).HSPGs consist of multiple heparan sulfate (HS) chains covalently linked to a limited set of core proteins. The HS chains contain repeating uronic acid and glucosamine disaccharide units. The binding functions of HSPGs depend on the fine structure of the attached heparan sulfate chains where sulfation modifications occur in four positions (N-, 3-O, and 6-O of glucosamine and 2-O of uronic acid) in highly variegated, yet highly regulated patterns (3, 4). 6-O-Sulfation of glucosamine is established to be critical for certain HSPG functions in organisms from Drosophila through mammals (5, 6).Several years ago, we cloned cDNAs encoding two novel extracellular sulfatases (Sulf-1 and Sulf-2) in human and mouse (7), initiated by the identification of the Sulf-1 ortholog in the quail embryo (QSulf-1) (8). We and others showed that both Sulfs are neutral pH endosulfatases, which remove glucosamine-6-O-sulfate from internal glucosamine residues of highly sulfated subregions within heparin/HSPGs (7, 9, 10). The ability of these enzymes to modulate the heparin/HSPG interactions of a number of growth factors, morphogens, and chemokines has been confirmed in direct binding assays (9, 11–13). In some cellular contexts, the Sulfs act to promote signaling pathways (Wnts, bone morphogenetic protein, and glial cell-derived neurotrophic factor) (9–11, 14), whereas in others the Sulfs are inhibitory (fibroblast growth factor-2 and transforming growth factor-β) (15–17). The importance of the Sulfs in development has been revealed by gene knockdown (8) and knock-out studies (11, 18–20). The phenotypes in single and double null mice include abnormalities in general growth, muscle innervation, muscle regeneration, skeletal tissue, and lung development. The Sulfs have been extensively investigated in cancer with some studies consistent with tumor suppression activity (15, 21, 22) and others with a pro-oncogenic role (23–25).As is the case for the prototypic QSulf-1 (8), HSulf-1 and HSulf-2 consist of four domains from the N to C terminus: a signal peptide, a catalytic domain of 374 amino acids, a basic hydrophilic domain of 346/366 amino acids, and a C-terminal domain of 109/127 amino acids (7, 8). In the 17-member sulfatase family (26), the Sulfs share the most extensive sequence homology with lysosomal glucosamine-6-sulfatase in the catalytic and C-terminal domains, although the centrally inserted hydrophilic domain is absent from this enzyme and other sulfatases. Limited information has been available about the proteolytic processing of the Sulfs during synthesis. In the present study, we show that the mature form of each human Sulf consists of a heterodimer of 75- and 50-kDa subunits, which is formed through the action of a furin-type proteinase on a proprotein of 125 kDa. We investigate the structural requirements for the enzymatic and signaling activities of these proteins.     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Vinyl Sulfones as Antiparasitic Agents and a Structural Basis for Drug Design

Cysteine proteases of the papain superfamily are implicated in a number of cellular processes and are important virulence factors in the pathogenesis of parasitic disease. These enzymes have therefore emerged as promising targets for antiparasitic drugs. We report the crystal structures of three major parasite cysteine proteases, cruzain, falcipain-3, and the first reported structure of rhodesain, in complex with a class of potent, small molecule, cysteine protease inhibitors, the vinyl sulfones. These data, in conjunction with comparative inhibition kinetics, provide insight into the molecular mechanisms that drive cysteine protease inhibition by vinyl sulfones, the binding specificity of these important proteases and the potential of vinyl sulfones as antiparasitic drugs.Sleeping sickness (African trypanosomiasis), caused by Trypanosoma brucei, and malaria, caused by Plasmodium falciparum, are significant, parasitic diseases of sub-Saharan Africa (1). Chagas'' disease (South American trypanosomiasis), caused by Trypanosoma cruzi, affects approximately, 16–18 million people in South and Central America. For all three of these protozoan diseases, resistance and toxicity to current therapies makes treatment increasingly problematic, and thus the development of new drugs is an important priority (2–4).T. cruzi, T. brucei, and P. falciparum produce an array of potential target enzymes implicated in pathogenesis and host cell invasion, including a number of essential and closely related papain-family cysteine proteases (5, 6). Inhibitors of cruzain and rhodesain, major cathepsin L-like papain-family cysteine proteases of T. cruzi and T. brucei rhodesiense (7–10) display considerable antitrypanosomal activity (11, 12), and some classes have been shown to cure T. cruzi infection in mouse models (11, 13, 14).In P. falciparum, the papain-family cysteine proteases falcipain-2 (FP-2)⁶ and falcipain-3 (FP-3) are known to catalyze the proteolysis of host hemoglobin, a process that is essential for the development of erythrocytic parasites (15–17). Specific inhibitors, targeted to both enzymes, display antiplasmodial activity (18). However, although the abnormal phenotype of FP-2 knock-outs is “rescued” during later stages of trophozoite development (17), FP-3 has proved recalcitrant to gene knock-out (16) suggesting a critical function for this enzyme and underscoring its potential as a drug target.Sequence analyses and substrate profiling identify cruzain, rhodesain, and FP-3 as cathepsin L-like, and several studies describe classes of small molecule inhibitors that target multiple cathepsin L-like cysteine proteases, some with overlapping antiparasitic activity (19–22). Among these small molecules, vinyl sulfones have been shown to be effective inhibitors of a number of papain family-like cysteine proteases (19, 23–27). Vinyl sulfones have many desirable attributes, including selectivity for cysteine proteases over serine proteases, stable inactivation of the target enzyme, and relative inertness in the absence of the protease target active site (25). This class has also been shown to have desirable pharmacokinetic and safety profiles in rodents, dogs, and primates (28, 29). We have determined the crystal structures of cruzain, rhodesain, and FP-3 bound to vinyl sulfone inhibitors and performed inhibition kinetics for each enzyme. Our results highlight key areas of interaction between proteases and inhibitors. These results help validate the vinyl sulfones as a class of antiparasitic drugs and provide structural insights to facilitate the design or modification of other small molecule inhibitor scaffolds.     

In Vivo Identification of Nuclear Factors Interacting with the Conserved Elements of Box C/D Small Nucleolar RNAs

The U16 small nucleolar RNA (snoRNA) is encoded by the third intron of the L1 (L4, according to the novel nomenclature) ribosomal protein gene of Xenopus laevis and originates from processing of the pre-mRNA in which it resides. The U16 snoRNA belongs to the box C/D snoRNA family, whose members are known to assemble in ribonucleoprotein particles (snoRNPs) containing the protein fibrillarin. We have utilized U16 snoRNA in order to characterize the factors that interact with the conserved elements common to the other members of the box C/D class. In this study, we have analyzed the in vivo assembly of U16 snoRNP particles in X. laevis oocytes and identified the proteins which interact with the RNA by label transfer after UV cross-linking. This analysis revealed two proteins, of 40- and 68-kDa apparent molecular size, which require intact boxes C and D together with the conserved 5′,3′-terminal stem for binding. Immunoprecipitation experiments showed that the p40 protein corresponds to fibrillarin, indicating that this protein is intimately associated with the RNA. We propose that fibrillarin and p68 represent the RNA-binding factors common to box C/D snoRNPs and that both proteins are essential for the assembly of snoRNP particles and the stabilization of the snoRNA.One of the most interesting recent findings related to ribosome biogenesis has been the identification of a large number of small RNAs localized in the nucleolus (snoRNAs). So far, more than 60 snoRNAs have been identified in vertebrates (17), and more than 30 have been identified in yeast (2). The total number of snoRNAs is not known, but it is likely to be close to 200 (33, 38). These snoRNAs, with the exception of the mitochondrial RNA processing (MRP) species (38), can be grouped into two major families on the basis of conserved structural and sequence elements. The first group includes molecules referred to as box C/D snoRNAs, whereas the second one comprises the species belonging to the box H/ACA family (2, 15).The two families differ in many aspects. The box C/D snoRNAs are functionally heterogeneous. Most of them function as antisense RNAs in site-specific ribose methylation of the pre-rRNA (1, 10, 17, 26); a minority have been shown to play a direct role in pre-rRNA processing in both yeast and metazoan cells (11, 21). The box C/D snoRNAs play their role by means of unusually long (up to 21 contiguous nucleotides) regions of complementarity to highly conserved sequences of 28S and 18S rRNAs (1). In contrast, several members of the H/ACA RNA family have been shown to direct site-specific isomerization of uridines into pseudouridines and to display shorter regions of complementarity to rRNA (14, 24). Mutational analysis suggests that H/ACA snoRNAs can also play a role as antisense RNAs by base pairing with complementary regions on rRNA (15, 24).Another difference between the two families can be seen by comparison of secondary structures. A Y-shaped motif, where a 5′,3′-terminal stem adjoins the C and D conserved elements, has been proposed for many box C/D snoRNAs (16, 26, 40, 42), whereas box H/ACA snoRNAs have been proposed to fold into two conserved hairpin structures connected by a single-stranded hinge region, followed by a short 3′ tail (15).Despite these differences, analogies have been found in the roles played by the conserved box elements. Mutational analysis and competition experiments indicated that C/D and H/ACA boxes are required both for processing and stable accumulation of the mature snoRNA, suggesting that they represent binding sites for specific trans-acting factors (2, 3, 8, 15, 16, 28, 36, 41).All snoRNAs are associated with proteins to form specific ribonucleoparticles (snoRNPs). The study of these particles began only recently, and so far, very few aspects of their structure and biosynthesis have been clarified. The only detailed analysis performed was on the mammalian U3 (19) and the yeast snR30 (20) snoRNPs. Of the identified components, a few appear to be more general factors: fibrillarin, which was shown to be associated with C/D snoRNPs (3, 4, 8, 13, 28, 31, 39), and the nucleolar protein GAR1, which was found associated with H/ACA snoRNAs in yeast (20). Just as the study of small nuclear RNP (snRNP) particles was crucial to the understanding of the splicing process, a detailed structural and functional analysis of snoRNP particles will be essential to elucidate the complex process of ribosome biosynthesis.In this study, we have analyzed the snoRNP assembly of wild-type and mutant U16 snoRNAs by following the kinetics of complex formation in the in vivo system of the Xenopus laevis oocyte. By a UV cross-linking technique, we have identified two proteins, of 40- and 68-kDa apparent molecular mass, which require intact boxes C and D together with the terminal stem for their binding. The 40-kDa species is specifically recognized by fibrillarin antibodies, indicating that this protein is intimately associated with the RNA.