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]     

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]     

Arachidonic Acid Stimulates Cell Adhesion through a Novel p38 MAPK-RhoA Signaling Pathway That Involves Heat Shock Protein 27

Rho GTPases are critical components of cellular signal transduction pathways. Both hyperactivity and overexpression of these proteins have been observed in human cancers and have been implicated as important factors in metastasis. We previously showed that dietary n-6 fatty acids increase cancer cell adhesion to extracellular matrix proteins, such as type IV collagen. Here we report that in MDA-MB-435 human melanoma cells, arachidonic acid activates RhoA, and inhibition of RhoA signaling with either C3 exoenzyme or dominant negative Rho blocked arachidonic acid-induced cell adhesion. Inhibition of the Rho kinase (ROCK) with either small molecule inhibitors or ROCK II-specific small interfering RNA (siRNA) blocked the fatty acid-induced adhesion. However, unlike other systems, inhibition of ROCK did not block the activation of p38 mitogen-activated protein kinase (MAPK); instead, Rho activation depended on p38 MAPK activity and the presence of heat shock protein 27 (HSP27), which is phosphorylated downstream of p38 after arachidonic acid treatment. HSP27 associated with p115RhoGEF in fatty acid-treated cells, and this association was blocked when p38 was inhibited. Furthermore, siRNA knockdown of HSP27 blocked the fatty acid-stimulated Rho activity. Expression of dominant negative p115-RhoGEF or p115RhoGEF-specific siRNA inhibited both RhoA activation and adhesion on type IV collagen, whereas a constitutively active p115RhoGEF restored the arachidonic acid stimulation in cells in which the p38 MAPK had been inhibited. These data suggest that n-6 dietary fatty acids stimulate a set of interactions that regulates cell adhesion through RhoA and ROCK II via a p38 MAPK-dependent association of HSP27 and p115RhoGEF.The ability of tumor cells to metastasize to secondary sites is a hallmark of neoplastic disease. Unfortunately, this propensity to spread is the primary cause of morbidity and death in cancer patients (1). Metastasis is clearly a highly regulated, multistep process that occurs in a spatiotemporal manner (2–4). To escape the restrictive compartment boundaries characteristic of adult tissue, separate intravasation and extravasation steps requiring alterations in co-adhesion, adhesion, invasion, and migration must occur. Execution of these biological processes, involving multiple proteins and cellular organelles, require highly coordinated cell signaling mechanisms.The Rho family of small GTPases regulates many facets of cytoskeletal rearrangements that facilitate cell attachment and migration (5–7). Rho GTPases act as molecular switches by changing from an inactive GDP-bound conformation to an active GTP-bound conformation, thereby regulating a signaling pathway. These proteins are directly regulated by Rho guanine nucleotide exchange factors (GEFs),² Rho GTPase activating proteins, and Rho GDP-dissociation inhibitors (8–12). RhoGEFs bind to the GTPase to catalyze the dissociation of GDP, allowing the binding of GTP and thereby promoting Rho activation (8). The RGS (regulators of G protein signaling) domain-containing RhoGEFs are a recently described family of GEFs. Currently, there are three members of this family, PDZ-RhoGEF, LARG, and p115RhoGEF (13–15), in which the RGS domains function as a heterotrimeric GTPase-activating domain (13, 15, 16). The RGS family of RhoGEFs has been shown to regulate Rho during several processes including cytoskeletal rearrangements, cell adhesion, and cancer progression (17–21).There is significant interplay between the activity of small GTPases and signaling derived from fatty acid metabolism (22–28). Linoleic acid, which is metabolized to arachidonic acid, is an n-6 polyunsaturated fatty acid that is present at high levels in most western diets (29). In animal models, diets high in n-6 polyunsaturated fatty acids have been shown to enhance tumor progression and metastasis (30, 31). Additionally, arachidonic acid is stored in cell membranes and is made available by phospholipases under conditions of increased inflammatory response (32). Arachidonic acid is further metabolized by cyclooxygenases (COX), lipoxygenases (LOX), and cytochrome P450 monooxygenases to yield bioactive products that have myriad effects on cells, and altered metabolism of arachidonic acid by COX, LOX, and P450 has been implicated in cancer progression (31, 33–36).We have studied mechanisms of cell adhesion using the MDA-MB-435 cells as a model of a highly metastatic human cancer cell line (37). These cells have been extensively studied for their ability to recapitulate the metastatic cascade in vivo and in vitro, although recent work indicates that the cells currently in use are most likely a human melanoma line (38). We initially observed that arachidonic acid (AA) enhanced adhesion of MDA-MB-435 cells to type IV collagen through specific integrin-mediated pathways (37). Exogenous AA led to the activation of mitogen-activated protein kinase (MAPK)-activated protein kinase 2 and the phosphorylation of heat shock protein 27 (HSP27) via a p38 MAPK-dependent process (39). Inhibition of p38 MAPK activation blocked cell adhesion as did function-blocking antibodies specific for subunits of the collagen receptor (40). More recently, we identified the key metabolite of AA (15-(S)- hydroxyeicosatetraenoic acid) and the upstream kinases (TAK1 and MKK6) that are responsible for activation of p38 MAPK in this system (41).In this study we investigated the role of Rho activation in the MDA-MB-435 cells after exposure to arachidonic acid. Several aspects of the regulation of Rho signaling in these cells provide insights into the cross-talk between important signaling pathways.     

Differential Regulation of Cell Type-specific Apoptosis by Stromelysin-3: A POTENTIAL MECHANISM VIA THE CLEAVAGE OF THE LAMININ RECEPTOR DURING TAIL RESORPTION IN XENOPUS LAEVIS*

Matrix metalloproteinases (MMPs) have been extensively studied because of their functional attributes in development and diseases. However, relatively few in vivo functional studies have been reported on the roles of MMPs in postembryonic organ development. Amphibian metamorphosis is a unique model for studying MMP function during vertebrate development because of its dependence on thyroid hormone (T3) and the ability to easily manipulate this process with exogenous T3. The MMP stromelysin-3 (ST3) is induced by T3, and its expression correlates with cell death during metamorphosis. We have previously shown that ST3 is both necessary and sufficient for larval epithelial cell death in the remodeling intestine. To investigate the roles of ST3 in other organs and especially on different cell types, we have analyzed the effect of transgenic overexpression of ST3 in the tail of premetamorphic tadpoles. We report for the first time that ST3 expression, in the absence of T3, caused significant muscle cell death in the tail of premetamorphic transgenic tadpoles. On the other hand, only relatively low levels of epidermal cell death were induced by precocious ST3 expression in the tail, contrasting what takes place during natural and T3-induced metamorphosis when ST3 expression is high. This cell type-specific apoptotic response to ST3 in the tail suggests distinct mechanisms regulating cell death in different tissues. Furthermore, our analyses of laminin receptor, an in vivo substrate of ST3 in the intestine, suggest that laminin receptor cleavage may be an underlying mechanism for the cell type-specific effects of ST3.The extracellular matrix (ECM),³ the dynamic milieu of the cell microenvironment, plays a critical role in dictating the fate of the cell. The cross-talk between the cell and ECM and the timely catabolism of the ECM are crucial for tissue remodeling during development (1). Matrix metalloproteinases (MMPs), extrinsic proteolytic regulators of the ECM, mediate this process to a large extent. MMPs are a large family of Zn²⁺-dependent endopeptidases potentially capable of cleaving the extracellular as well as nonextracellular proteins (2–9). The MMP superfamily includes collagenases, gelatinases, stromelysins, and membrane-type MMPs based on substrate specificity and domain organization (2–4). MMPs have been implicated to influence a wide range of physiological and pathological processes (10–13). The roles of MMPs appear to be very complex. For example, MMPs have been suggested to play roles in both tumor promotion and suppression (13–19). Unfortunately, relatively few functional studies have been carried out in vivo, especially in relation to the mechanisms involved during vertebrate development.Amphibian metamorphosis presents a fascinating experimental model to study MMP function during postembryonic development. A unique and salient feature of the metamorphic process is the absolute dependence on the signaling of thyroid hormone (20–23). This makes it possible to prevent metamorphosis by simply inhibiting the synthesis of endogenous T3 or to induce precocious metamorphosis by merely adding physiological levels of T3 in the rearing water of premetamorphic tadpoles. Gene expression screens have identified the MMP stromelysin-3 (ST3) as a direct T3 response gene (24–27). Expression studies have revealed a distinct spatial and temporal ST3 expression profile in correlation with metamorphic event, especially cell death (25, 28–31). Organ culture studies on intestinal remodeling have directly substantiated an essential role of ST3 in larval epithelial cell death and ECM remodeling (32). Furthermore, precocious expression of ST3 alone in premetamorphic tadpoles through transgenesis is sufficient to induce ECM remodeling and larval epithelial apoptosis in the tadpole intestine (33). Thus, ST3 appears to be necessary and sufficient for intestinal epithelial cell death during metamorphosis.ST3 was first isolated as a breast cancer-associated gene (34), and unlike most other MMPs, ST3 is secreted as an active protease through a furin-dependent intracellular activation mechanism (35). Like many other MMPs, ST3 is expressed in a number of pathological processes, including most human carcinomas (11, 36–40), as well as in many developmental processes in mammals (10, 34, 41–43), although the physiological and pathological roles of ST3 in vivo are largely unknown in mammals. Interestingly, compared with other MMPs, ST3 has only weak activities toward ECM proteins in vitro but stronger activities against non-ECM proteins like α1 proteinase inhibitor and IGFBP-1 (44–46). Although ST3 may cleave ECM proteins strongly in the in vivo environment, these findings suggest that the cleavage of non-ECM proteins is likely important for its biological roles. Consistently, we have recently identified a cell surface receptor, laminin receptor (LR) as an in vivo substrate of ST3 in the tadpole intestine during metamorphosis (47–49). Analyses of LR expression and cleavage suggest that LR cleavage by ST3 is likely an important mechanism by which ST3 regulates the interaction between the larval epithelial cells and the ECM to induce cell death during intestinal remodeling (47, 48).Here, to investigate the role of ST3 in the apoptosis in other tissues during metamorphosis and whether LR cleavage serves as a mechanism for ST3 to regulate the fate of different cell types, we have analyzed the effects of precocious expression of ST3 in premetamorphic tadpole tail. The tail offers an opportunity to examine the effects of ST3 on different cell types. The epidermis, the fast and slow muscles, and the connective tissue underlying the epidermis in the myotendinous junctions and surrounding the notochord constitute the major tissue types in tail (50). Even though death is the destiny of all these cell types, it is not clear whether they all die through similar or different mechanisms. Microscopic and histochemical analyses have shown that at least the muscle and epidermal cells undergo T3-dependent apoptosis during metamorphosis (23, 29, 51, 52). To study whether ST3 regulates apoptosis of these two cell types, we have made use of the transgenic animals that express a transgenic ST3 under the control of a heat shock-inducible promoter (33). We show that whereas extensive apoptosis is present in both the epidermis and muscles during natural as well as T3-induced metamorphosis, transgenic expression of ST3 induces cell death predominantly in the muscles. Furthermore, we show that LR is expressed in the epidermis and connective tissue but not in muscles of the tadpole tail. More importantly, LR cleavage products are present in the tail during natural metamorphosis but not in transgenic tadpoles overexpressing ST3. These results suggest that ST3 has distinct effects on the epidermis and muscles in the tail, possibly because of the tissue-specific expression and function of LR.     

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 Mutational Analysis of YtvA from Bacillus subtilis: MECHANISM OF LIGHT ACTIVATION OF THE GENERAL STRESS RESPONSE

The general stress response of Bacillus subtilis can be activated by stimuli such as the addition of salt or ethanol and with blue light. In the latter response, YtvA activates σ^B through a cascade of Rsb proteins, organized in stressosomes. YtvA is composed of an N-terminal LOV (light, oxygen, and voltage) domain and a C-terminal STAS (sulfate transporter and anti-sigma factor) domain and shows light-modulated GTP binding in vitro. Here, we examine the mechanism of YtvA-mediated activation of σ^B in vivo with site-directed mutagenesis. Constitutive off and constitutive on mutations have been identified. Disruption of GTP binding in the STAS domain eliminates light activation of σ^B. In contrast, modification of sites relevant for phosphorylation of STAS domains does not affect the stress response significantly. The data obtained are integrated into a model for the structure of full-length YtvA, which presumably functions as a dimer.LOV² domains (1), members of the superfamily of PAS domains (2, 3), are abundant in all domains of life and were first identified in plant phototropins (4). These photoreceptors regulate stomatal opening, phototropism, etc. and contain two N-terminal LOV domains that confer light regulation on the C-terminal Ser/Thr kinase domain (4). They also occur in bacteria, in which YtvA from Bacillus subtilis has been best characterized (for a review, see e.g. Ref. 5). Its N-terminal LOV domain binds FMN and shows the typical LOV photochemistry (6, 7): covalent adduct formation between a cysteine and the FMN chromophore. A linker helix, denoted Jα (7), connects the LOV domain to a STAS domain. The latter domain is present in many regulators of the general stress response of B. subtilis (8, 9). Stress via the addition of salt or ethanol (for a review, see Ref. 10) and blue light (11, 12) activates the general stress response via the environmental pathway, which integrates various signals via a large multiprotein complex, called the stressosome (13, 14). YtvA, which mediates light activation of σ^B (11, 12, 15), co-purifies with other STAS domain proteins in the stressosomes (16).When cells are stressed, STAS domains of several stressosome proteins (e.g. RsbS and RsbR) are phosphorylated by another intrinsic stressosome component, the serine/threonine kinase RsbT (9, 14, 17, 18). Next, RsbT is released from the complex to trigger RsbU, a protein phosphatase, thus (indirectly) activating σ^B (19). Phosphorylation of YtvA, however, has never been detected. Rather, it has been demonstrated in vitro that YtvA shows light-dependent GTP binding, presumably at its NTP-binding site in the STAS domain (20).Little is known about the mechanism of signal transmission in and by YtvA, except that in the C62A mutant, photochemistry in vitro (12) and light activation of σ^B in vivo (12, 15) are abolished. More detailed information is available for LOV domains of phototropins. A conserved glutamine, which is in hydrogen-bonding contact with the isoalloxazine ring of FMN, rotates its side chain by 180° upon covalent adduct formation (21). Replacement of this residue by leucine in the LOV2 domain of Phy3 from Adiantum results in a considerable reduction of the light-induced structural change (22). The corresponding mutation in phototropin 1 from Arabidopsis impairs autophosphorylation activity (23). The signal generated in the LOV2 domain is transmitted to the downstream kinase domain of phototropin 1 of Avena sativa through disruption of the interaction between its central β-sheet and the C-terminal linker region, the Jα-helix (24).Here, we study the mechanism of activation of YtvA in vivo, i.e. light-induced activation of the σ^B response, with site-directed mutagenesis. We focus on three regions of the protein, the flavin-binding pocket, the β-sheet of the LOV domain, and the GTP-binding site, and on potential phosphorylation sites of the STAS domain. We demonstrate that light-activated GTP binding is crucial for functional YtvA. A computational approach was used to model the structure of full-length YtvA. The model suggests that light modulates accessibility of the GTP-binding site of the STAS domain of YtvA.