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]     

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]     

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.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[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]     

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]     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[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]     

Evidence that Antibody-Mediated Neutralization of Human Immunodeficiency Virus Type 1 by Sera from Infected Individuals Is Independent of Coreceptor Usage

Human immunodeficiency virus type 1 (HIV-1) uses a variety of chemokine receptors as coreceptors for virus entry, and the ability of the virus to be neutralized by antibody may depend on which coreceptors are used. In particular, laboratory-adapted variants of the virus that use CXCR4 as a coreceptor are highly sensitive to neutralization by sera from HIV-1-infected individuals, whereas primary isolates that use CCR5 instead of, or in addition to, CXCR4 are neutralized poorly. To determine whether this dichotomy in neutralization sensitivity could be explained by differential coreceptor usage, virus neutralization by serum samples from HIV-1-infected individuals was assessed in MT-2 cells, which express CXCR4 but not CCR5, and in mitogen-stimulated human peripheral blood mononuclear cells (PBMC), where multiple coreceptors including CXCR4 and CCR5 are available for use. Our results showed that three of four primary isolates with a syncytium-inducing (SI) phenotype and that use CXCR4 and CCR5 were neutralized poorly in both MT-2 cells and PBMC. The fourth isolate, designated 89.6, was more sensitive to neutralization in MT-2 cells than in PBMC. We showed that the neutralization of 89.6 in PBMC was not improved when CCR5 was blocked by having RANTES, MIP-1α, and MIP-1β in the culture medium, indicating that CCR5 usage was not responsible for the decreased sensitivity to neutralization in PBMC. Consistent with this finding, a laboratory-adapted strain of virus (IIIB) was significantly more sensitive to neutralization in CCR5-deficient PBMC (homozygous Δ32-CCR5 allele) than were two of two SI primary isolates tested. The results indicate that the ability of HIV-1 to be neutralized by sera from infected individuals depends on factors other than coreceptor usage.Human immunodeficiency virus type 1 (HIV-1), the etiologic agent of AIDS, utilizes the HLA class II receptor, CD4, as its primary receptor to gain entry into cells (17, 30). Entry is initiated by a high-affinity interaction between CD4 and the surface gp120 of the virus (32). Subsequent to this interaction, conformational changes that permit fusion of the viral membrane with cellular membranes occur within the viral transmembrane gp41 (9, 58, 59). In addition to CD4, one or more recently described viral coreceptors are needed for fusion to take place. These coreceptors belong to a family of seven-transmembrane G-protein-coupled proteins and include the CXC chemokine receptor CXCR4 (3, 4, 24, 44), the CC chemokine receptors CCR5 (1, 12, 13, 18, 21, 23, 45) and, less commonly, CCR3 and CCR2b (12, 21), and two related orphan receptors termed BONZO/STRL33 and BOB (19, 34). Coreceptor usage by HIV-1 can be blocked by naturally occurring ligands, including SDF-1 for CXCR4 (4, 44), RANTES, MIP-1α, and MIP-1β in the case of CCR5 (13, 45), and eotaxin for CCR3 (12).The selective cellular tropisms of different strains of HIV-1 may be determined in part by coreceptor usage. For example, all culturable HIV-1 variants replicate initially in mitogen-stimulated human peripheral blood mononuclear cells (PBMC), but only a minor fraction are able to infect established CD4⁺ T-cell lines (43). This differential tropism is explained by the expression of CXCR4 together with CCR5 and other CC chemokine coreceptors on PBMC and the lack of expression of CCR5 on most T-cell lines (5, 10, 19, 35, 39, 50, 53). Indeed, low-passage field strains (i.e., primary isolates) of HIV-1 that fail to replicate in T-cell lines use CCR5 as their major coreceptor and are unable to use CXCR4 (1, 12, 18, 21, 23, 28). Because these isolates rarely produce syncytia in PBMC and fail to infect MT-2 cells, they are often classified as having a non-syncytium-inducing (NSI) phenotype. Primary isolates with a syncytium-inducing (SI) phenotype are able to use CXCR4 alone or, more usually, in addition to CCR5 (16, 20, 51). HIV-1 variants that have been passaged multiple times in CD4⁺ T-cell lines, and therefore considered to be laboratory adapted, exhibit a pattern of coreceptor usage that resembles that of SI primary isolates. Most studies have shown that the laboratory-adapted strain IIIB uses CXCR4 alone (3, 13, 20, 24, 51) and that MN and SF-2 use CXCR4 primarily and CCR5 to a lesser degree (11, 13). Sequences within the V3 loop of gp120 have been shown to be important, either directly or indirectly, for the interaction of HIV-1 with both CXCR4 (52) and CCR5 (12, 14, 54, 60). This region of gp120 contains multiple determinants of cellular tropism (43) and is a major target for neutralizing antibodies to laboratory-adapted HIV-1 but not to primary isolates (29, 46, 57).It has been known for some time that the ability of sera from HIV-1-infected individuals to neutralize laboratory-adapted strains of HIV-1 does not predict their ability to neutralize primary isolates in vitro (7). In general, the former viruses are highly sensitive to neutralization whereas the latter viruses are neutralized poorly by antibodies induced in response to HIV-1 infection (7, 43). Importantly, neutralizing antibodies generated by candidate HIV-1 subunit vaccines have been highly specific for laboratory-adapted viruses (26, 37, 38). In principle, the dichotomy in neutralization sensitivity between these two categories of virus could be related to coreceptor usage. To test this, we investigated whether the use of CXCR4 in the absence of CCR5 would render SI primary isolates highly sensitive to neutralization in vitro by sera from HIV-1-infected individuals. Two similar studies using human monoclonal antibodies and soluble CD4 have been reported (31a, 55).     

Question Processing and Clustering in INDOC: A Biomedical Question Answering System

The exponential growth in the volume of publications in the biomedical domain has made it impossible for an individual to keep pace with the advances. Even though evidence-based medicine has gained wide acceptance, the physicians are unable to access the relevant information in the required time, leaving most of the questions unanswered. This accentuates the need for fast and accurate biomedical question answering systems. In this paper we introduce INDOC—a biomedical question answering system based on novel ideas of indexing and extracting the answer to the questions posed. INDOC displays the results in clusters to help the user arrive the most relevant set of documents quickly. Evaluation was done against the standard OHSUMED test collection. Our system achieves high accuracy and minimizes user effort.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]     

Differential Regulation of Elastic Fiber Formation by Fibulin-4 and -5

Fibulin-4 and -5 are extracellular glycoproteins with essential non-compensatory roles in elastic fiber assembly. We have determined how they interact with tropoelastin, lysyl oxidase, and fibrillin-1, thereby revealing how they differentially regulate assembly. Strong binding between fibulin-4 and lysyl oxidase enhanced the interaction of fibulin-4 with tropoelastin, forming ternary complexes that may direct elastin cross-linking. In contrast, fibulin-5 did not bind lysyl oxidase strongly but bound tropoelastin in terminal and central regions and could concurrently bind fibulin-4. Both fibulins differentially bound N-terminal fibrillin-1, which strongly inhibited their binding to lysyl oxidase and tropoelastin. Knockdown experiments revealed that fibulin-5 controlled elastin deposition on microfibrils, although fibulin-4 can also bind fibrillin-1. These experiments provide a molecular account of the distinct roles of fibulin-4 and -5 in elastic fiber assembly and how they act in concert to chaperone cross-linked elastin onto microfibrils.Fibulins are a family of extracellular glycoproteins containing contiguous calcium-binding epidermal growth factor-like domains (cbEGFs)³ and a characteristic C-terminal fibulin (FC) domain (1–3). Recent studies have revealed that fibulin-4 and -5 are both essential for elastic fiber formation (4–7). Fibulin-4 is widely expressed from early embryogenesis and is necessary for normal vascular, lung, and skin development, since mice that lack fibulin-4 do not form elastic fibers and die perinatally (5). Furthermore, mice with reduced fibulin-4 expression develop aneurysms (8). Fibulin-5 is abundant in the aorta and large arteries during embryogenesis and following vascular injury (9, 10). Lack of fibulin-5 causes a less severe phenotype, with viable homozygous mice, but the elastic fibers in skin, lungs, and aorta are irregular and fragmented (6, 7), and there is altered vascular remodeling (11). These mice models also highlight that fibulin-4 and -5 have non-compensatory roles in elastic fiber formation. Mutations in both molecules can cause cutis laxa, a heritable disorder associated with elastic fiber degeneration leading to sagging skin, vascular tortuosity, and emphysematous lungs (12–15). A third isoform, fibulin-3, may play a minor role in elastic fiber formation, since its deficiency disrupts elastic fibers in Bruch''s membrane of the eye (16) and vaginal tissues (17).Elastic fiber formation is a complex multistep process (18–20). Initial pericellular microassembly of tropoelastin, which may involve the 67-kDa elastin-binding protein receptor, generates elastin globules that are stabilized by desmosine cross-links catalyzed mainly by lysyl oxidase (LOX) but also by LOXL1 (LOX-like 1). These globules are deposited on a fibrillin microfibril template, where they coalesce and undergo further cross-linking to form the elastin core of mature fibers. The ability of fibulin-4 and -5 to bind tropoelastin and fibrillin-1, the major structural component of microfibrils, supports a model in which these fibulins direct elastin deposition on microfibrils (4–7, 21–25). This model does not delineate the unique molecular contributions of fibulin-4 and -5 to elastic fiber formation, but some molecular differences have emerged. Tropoelastin was bound more strongly by fibulin-5 than by fibulin-4, whereas fibulin-5 was at the microfibril-elastin interface, but perichondrial fibulin-4 localized mainly to microfibrils (4).Fibulin-4 null mice offer tantalizing clues to how fibulin-4 contributes to elastic fiber formation (5). They had dramatically reduced (94%) desmosine cross-links despite no change in elastin or LOX expression levels, and electron-dense rodlike structures were prominent within elastin aggregates. Morphologically similar structures seen after chemically inhibiting LOX were previously identified as glycosaminoglycans, which can bind charged free ϵ-amino groups on lysines in tropoelastin (26). However, fibulin-4^+/− mice showed ∼20% increase in desmosine (5). LOX-null mice have phenotypic features similar to those of fibulin-4 null mice, dying perinatally with 60% reduced desmosine cross-links and major abnormalities in vascular and other elastic tissues (27, 28). In contrast, LOXL1-null mice are viable but have reduced desmosine (29), whereas fibulin-5 null mice have a 16% reduction in desmosine cross-links and survive well into adulthood (7). Detection of the LOXL1 pro-domain in fibulin-5 null mice skin but not wild-type skin implicates fibulin-5 in activation of LOXL1 (30).We and others have shown that fibrillin-1 and the microfibrillar protein MAGP-1 can both directly bind tropoelastin (31–34). However, the fibulin-null mice show that the fibrillin-1 interaction with tropoelastin is insufficient to support elastic fiber formation in vivo. Fibulin-5 has been reported to facilitate tropoelastin binding to the N-terminal half of fibrillin-1 (21). A study of elastin polypeptide self-assembly through coacervation and maturation phases showed that, although the N-terminal half of fibrillin-1 increased maturation velocity and droplet clustering, fibulin-4 and -5 both slowed maturation and limited globule growth (35). These studies imply that fibulins and fibrillin-1 act together to regulate elastin accretion on microfibrils.To gain further insights into the contributions of fibulin-4 and -5 to elastic fiber formation, we have delineated how they interact with tropoelastin, LOX, and fibrillin-1. Novel findings are that fibulin-4 directly binds LOX, and this interaction enhances fibulin-4 binding to tropoelastin, thus forming a ternary complex that may be critical for elastin cross-linking. Fibulin-5 can concurrently bind fibulin-4 and tropoelastin, but the interaction of both fibulins with fibrillin-1 strongly inhibits their binding to tropoelastin. These interactions indicate the molecular basis of how fibulins act as chaperones for deposition of elastin onto microfibrils. Our study thus provides a molecular account of the differential roles of fibulins-4 and -5 in elastic fiber formation.     

Photoprotective Energy Dissipation in Higher Plants Involves Alteration of the Excited State Energy of the Emitting Chlorophyll(s) in the Light Harvesting Antenna II (LHCII)

Non-photochemical quenching (NPQ), a mechanism of energy dissipation in higher plants protects photosystem II (PSII) reaction centers from damage by excess light. NPQ involves a reduction in the chlorophyll excited state lifetime in the PSII harvesting antenna (LHCII) by a quencher. Yet, little is known about the effect of the quencher on chlorophyll excited state energy and dynamics. Application of picosecond time-resolved fluorescence spectroscopy demonstrated that NPQ involves a red-shift (60 ± 5 cm⁻¹) and slight enhancement of the vibronic satellite of the main PSII lifetime component present in intact chloroplasts. Whereas this fluorescence red-shift was enhanced by the presence of zeaxanthin, it was not dependent upon it. The red-shifted fluorescence of intact chloroplasts in the NPQ state was accompanied by red-shifted chlorophyll a absorption. Nearly identical absorption and fluorescence changes were observed in isolated LHCII complexes quenched in a low detergent media, suggesting that the mechanism of quenching is the same in both systems. In both cases, the extent of the fluorescence red-shift was shown to correlate with the lifetime of a component. The alteration in the energy of the emitting chlorophyll(s) in intact chloroplasts and isolated LHCII was also accompanied by changes in lutein 1 observed in their 77K fluorescence excitation spectra. We suggest that the characteristic red-shifted fluorescence emission reflects an altered environment of the emitting chlorophyll(s) in LHCII brought about by their closer interaction with lutein 1 in the quenching locus.The photosynthetic apparatus of plants and algae has evolved a process of feedback control of energy input into the photosynthetic reaction centers, revealed in observations of the non-photochemical quenching of chlorophyll fluorescence (NPQ)² (1–3). Exposure to high light intensity leads to closure of a large fraction of photosystem II (PSII) reaction centers and the buildup of the proton gradient across the photosynthetic membrane (ΔpH). The latter prompts the transition of the PSII light harvesting antenna (LHCII) into the photoprotective state (4–6), where a large proportion of the absorbed energy is rapidly converted into heat, preventing photoinhibitory damage of reaction centers. The molecular mechanism of NPQ is currently under scrupulous multidisciplinary investigation (7–14).Many of these studies have centered on the search for pigments associated with light harvesting complexes that may act as energy traps/quenchers. In recent years certain evidence has accumulated that xanthophylls may play a central role in photoprotective energy dissipation. A carotenoid radical cation, suggested to be zeaxanthin, was found to correlate with the level of NPQ and was proposed to act as a direct quencher of chlorophyll excitation (7, 10). Other studies indicated the involvement of energy transfer to the S₁ state of LHCII-bound lutein 1 (9, 15). Evidence of these putative quenchers occurring in vivo was provided (7, 9), yet little is known about how ΔpH activates them. ΔpH is believed to protonate luminal residues on the PSII subunit PsbS and LHC complexes causing certain conformational changes to occur, which are necessary for NPQ (5, 6). These conformational changes are not well understood but can be traced by changes in LHCII-bound pigments detected by absorption and resonance Raman spectroscopy (9, 15–20).Both the proposed zeaxanthin- and lutein-quenching mechanisms assume very close (van der Waals) contact between the quenching xanthophyll and chlorophyll. Such an interaction would, most likely, affect the energy levels of the pigments involved, which may be required to convert the xanthophyll into an energy quencher. Indeed, a subpopulation of zeaxanthin molecules in PSII was reported to undergo an absorption red-shift in the NPQ state (18–20). Recently, lutein absorption was also found to be decreased in the dissipative state (20), interestingly this loss of molar extinction was enhanced by the presence of zeaxanthin.Whereas many studies have probed the role of xanthophylls, relatively little is known about the state of chlorophylls upon formation of NPQ. Some research suggests that chlorophyll associates may be formed in the NPQ state, which potentially could act as quenchers (21–23). It has been proposed that such associates have red-shifted absorption and fluorescence spectra. Low temperature (77K) fluorescence spectra of quenched LHCII aggregates possess a 20-nm red-shifted band, F700 (22), which was also detected in leaves and thylakoids frozen in the NPQ state (21, 24). The F700 emitting chlorophylls appeared to possess very little or no absorption, indicating they may arise from chlorophyll-chlorophyll excitonic interactions (22). Later, low temperature (77K) time-resolved fluorescence spectra demonstrated F700 was one of several red-shifted bands present in aggregated LHCII, however, because they all possessed much longer lifetimes than the main F680 band it was considered unlikely that they could act as quenchers (25, 26). Indeed, the most red-shifted of these bands, F700, was absent in non-aggregated LHCII quenched in a gel media (15). F700 was therefore suggested to reflect intersubunit LHCII chlorophyll interactions rather than being directly related the quencher (15). While the latter study did not exclude the event of LHCII aggregation during NPQ formation, it did demonstrate that an intrinsic change in pigment interactions rather than protein aggregation per se is the cause of quenching. Nonetheless, a small red-shift of both the chlorophyll a absorption and fluorescence emission at room temperature was still observed and found to be in good correlation with the strength of quenching (15).To investigate the effect of NPQ on chlorophyll excited state energy and dynamics we have undertaken a combined time-resolved and steady state spectroscopic study on isolated intact chloroplasts. Common quenching-induced features were found in the time-resolved and excitation fluorescence spectra in both chloroplasts and LHCII, suggesting a common mechanism of quenching involving changes in the environment of the emitting chlorophyll(s) and lutein 1.