BAliBASE (Benchmark Alignment dataBASE): enhancements for repeats, transmembrane sequences and circular permutations

BAliBASE is specifically designed to serve as an evaluation resource to address all the problems encountered when aligning complete sequences. The database contains high quality, manually constructed multiple sequence alignments together with detailed annotations. The alignments are all based on three-dimensional structural superpositions, with the exception of the transmembrane sequences. The first release provided sets of reference alignments dealing with the problems of high variability, unequal repartition and large N/C-terminal extensions and internal insertions. Here we describe version 2.0 of the database, which incorporates three new reference sets of alignments containing structural repeats, trans­membrane sequences and circular permutations to evaluate the accuracy of detection/prediction and alignment of these complex sequences. BAliBASE can be viewed at the web site http://www-igbmc.u-strasbg.fr/BioInfo/BAliBASE2/index.html or can be downloaded from ftp://ftp-igbmc.u-strasbg.fr/pub/BAliBASE2/.

Random circular permutation of genes and expressed polypeptide chains: application of the method to the catalytic chains of aspartate transcarbamoylase.

Recent studies on proteins whose N and C termini are in close proximity have demonstrated that folding of polypeptide chains and assembly of oligomers can be accomplished with circularly permuted chains. As yet no methodical study has been conducted to determine how extensively new termini can be introduced and where such termini cannot be tolerated. We have devised a procedure to generate random circular permutations of the catalytic chains of Escherichia coli aspartate transcarbamoylase (ATCase; EC 2.1.3.2) and to select clones that produce active or stable holoenzyme containing permuted chains. A tandem gene construct was made, based on the desired linkage between amino acid residues in the C- and N-terminal regions of the polypeptide chain, and this DNA was treated with a suitable restriction enzyme to yield a fragment containing the rearranged coding sequence for the chain. Circularization achieved with DNA ligase, followed by linearization at random with DNase I, and incorporation of the linearized, repaired, blunt-ended, rearranged genes into a suitable plasmid permitted the expression of randomly permuted polypeptide chains. The plasmid with appropriate stop codons also contained pyrI, the gene encoding the regulatory chain of ATCase. Colonies expressing detectable amounts of ATCase-like molecules containing permuted catalytic chains were identified by an immunoblot technique or by their ability to grow in the absence of pyrimidines in the growth medium. Sequencing of positive clones revealed a variety of novel circular permutations. Some had N and C termini within helices of the wild-type enzyme as well as deletions and insertions. Permutations were concentrated in the C-terminal domain and only few were detected in the N-terminal domain. The technique, which is adaptable generally to proteins whose N and C termini are near each other, can be of value in relating in vivo folding of nascent, growing polypeptide chains to in vitro renaturation of complete chains and determining the role of protein sequence in folding kinetics.

Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: correlation with exon skipping.

The cytochrome P450 2C24 gene is characterized by the capability to generate, in rat kidney, a transcript containing exons 2 and 4 spliced at correct sites but having the donor site of exon 4 directly joined to the acceptor site of exon 2 (exon scrambling). By reverse transcriptase-PCR analysis, it is now shown that the only exons present in the scrambled transcript are exons 2, 3, and 4 and that this molecule lacks a poly(A)+ tail. Furthermore, the use of PCR primers in both orientations of either exon 2 or exon 4 revealed that the orders of the exons in the scrambled transcript are 2-3-4-2 and 4-2-3-4, respectively. These results, combined with the observation that P450 2C24 is a single-copy gene, with no duplication of the exon 2 to exon 4 segment, suggest that the scrambled transcript has properties consistent with that of a circular molecule. In line with this is the observation of an increased resistance of the transcript to phosphodiesterase I, a 3'-exonuclease. Moreover, an alternatively processed cytochrome P450 2C24 mRNA, lacking the three scrambled exons and having exon 1 directly joined to exon 5, has been identified in kidney and liver, tissues that express the scrambled transcript. This complete identity of the exons that are absent in the alternatively processed mRNA but present in the scrambled transcript is interpreted as indicative of the possibility that exon scrambling and exon skipping might be interrelated phenomena. It is therefore proposed that alternative pre-mRNA processing has the potential to generate not only mRNAs lacking one or more exons but also circular RNA molecules.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

Origin and evolution of circular waves and spirals in Dictyostelium discoideum territories.

Randomly distributed Dictyostelium discoideum cells form cooperative territories by signaling to each other with cAMP. Cells initiate the process by sending out pulsatile signals, which propagate as waves. With time, circular and spiral patterns form. We show that by adding spatial and temporal noise to the levels of an important regulator of external cAMP levels, the cAMP phosphodiesterase inhibitor, we can explain the natural progression of the system from randomly firing cells to circular waves whose symmetries break to form double- and single- or multi-armed spirals. When phosphodiesterase inhibitor is increased with time, mimicking experimental data, the wavelength of the spirals shortens, and a proportion of them evolve into pairs of connected spirals. We compare these results to recent experiments, finding that the temporal and spatial correspondence between experiment and model is very close.

Mutagenesis in the C-terminal region of human interleukin 5 reveals a central patch for receptor alpha chain recognition.

Cassette mutagenesis was used to identify side chains in human interleukin 5 (hIL-5) that mediate binding to hIL-5 receptor alpha chain (hIL-5R alpha). A series of single alanine substitutions was introduced into a stretch of residues in the C-terminal region, including helix D, which previously had been implicated in receptor alpha chain recognition and which is aligned on the IL-5 surface so as to allow the topography of receptor binding residues to be examined. hIL-5 and single site mutants were expressed in COS cells, their interactions with hIL-5R alpha were measured by a sandwich surface plasmon resonance biosensor method, and their biological activities were measured by an IL-5-dependent cell proliferation assay. A pattern of mutagenesis effects was observed, with greatest impact near the interface between the two four-helix bundles of IL-5, in particular at residues Glu-110 and Trp-111, and least at the distal ends of the D helices. This pattern suggests the possibility that residues near the interface of the two four-helix bundles in hIL-5 comprise a central patch or hot spot, which constitutes an energetically important alpha chain recognition site. This hypothesis suggests a structural explanation for the 1:1 stoichiometry observed for the complex of hIL-5 with hIL-5R alpha.

Real-time 3D semi-local surface patch extraction using GPGPU

Feature vectors can be anything from simple surface normals to more complex feature descriptors. Feature extraction is important to solve various computer vision problems: e.g. registration, object recognition and scene understanding. Most of these techniques cannot be computed online due to their complexity and the context where they are applied. Therefore, computing these features in real-time for many points in the scene is impossible. In this work, a hardware-based implementation of 3D feature extraction and 3D object recognition is proposed to accelerate these methods and therefore the entire pipeline of RGBD based computer vision systems where such features are typically used. The use of a GPU as a general purpose processor can achieve considerable speed-ups compared with a CPU implementation. In this work, advantageous results are obtained using the GPU to accelerate the computation of a 3D descriptor based on the calculation of 3D semi-local surface patches of partial views. This allows descriptor computation at several points of a scene in real-time. Benefits of the accelerated descriptor have been demonstrated in object recognition tasks. Source code will be made publicly available as contribution to the Open Source Point Cloud Library.

Data analysis using circular causality in networks

Complex systems in causal relationships are known to be circular rather than linear; this means that a particular result is not produced by a single cause, but rather that both positive and negative feedback processes are involved. However, although interpreting systemic interrelationships requires a language formed by circles, this has only been developed at the diagram level, and not from an axiomatic point of view. The first difficulty encountered when analysing any complex system is that usually the only data available relate to the various variables, so the first objective was to transform these data into cause-and-effect relationships. Once this initial step was taken, our discrete chaos theory could be applied by finding the causal circles that will form part of the system attractor and allow their behavior to be interpreted. As an application of the technique presented, we analyzed the system associated with the transcription factors of inflammatory diseases.

Shaking-table tests of flat-bottom circular silos containing grain-like material

According to Eurocode 8, the seismic design of flat-bottom circular silos containing grain-like material is based on a rough estimate of the inertial force imposed on the structure by the ensiled content during an earthquake: 80% of the mass of the content multiplied by the peak ground acceleration. A recent analytical consideration of the horizontal shear force mobilised within the ensiled material during an earthquake proposed by some of the authors has resulted in a radically reduced estimate of this load suggesting that, in practice, the effective mass of the content is significantly less than that specified. This paper describes a series of laboratory tests that featured shaking table and a silo model, which were conducted in order to obtain some experimental data to verify the proposed theoretical formulations and to compare with the established code provisions. Several tests have been performed with different heights of ensiled material – about 0.5 mm diameter Ballotini glass – and different magnitudes of grain–wall friction. The results indicate that in all cases, the effective mass is indeed lower than the Eurocode specification, suggesting that the specification is overly conservative, and that the wall–grain friction coefficient strongly affects the overturning moment at the silo base. At peak ground accelerations up to around 0.35 g, the proposed analytical formulation provides an improved estimate of the inertial force imposed on such structures by their contents.

The deterioration of Circular Mausoleum, Roman Necropolis of Carmona, Spain

The Circular Mausoleum tomb in the Roman Necropolis of Carmona was carved on a calcarenite sequence in an ancient quarry located in the town of Carmona, Southern Spain. This rock-cut tomb, representative of Roman burial practices, currently suffers from serious deterioration. A detailed survey over several years permitted the identification of the main tomb's pathologies and damaging processes, which include loss of material (scaling, flaking, granular disintegration), surface modifications (efflorescences, crusts and deposits) and extensive biological colonization. The results obtained in this study indicated that anthropogenic changes were largely responsible and enhanced the main alteration mechanisms observed in the Circular Mausoleum. Based on the deterioration diagnosis, effective corrective actions were proposed. This study shows that any conservative intervention in the interior of the tomb should be preceded by accurate in situ measurements and laboratory analyses to ascribe the source of the deterioration damages and thus designing effective treatments.