1001 optimal PDB structure alignments: integer programming methods for finding the maximum contact map overlap.

Protein structure comparison is a fundamental problem for structural genomics, with applications to drug design, fold prediction, protein clustering, and evolutionary studies. Despite its importance, there are very few rigorous methods and widely accepted similarity measures known for this problem. In this paper we describe the last few years of developments on the study of an emerging measure, the contact map overlap (CMO), for protein structure comparison. A contact map is a list of pairs of residues which lie in three-dimensional proximity in the protein's native fold. Although this measure is in principle computationally hard to optimize, we show how it can in fact be computed with great accuracy for related proteins by integer linear programming techniques. These methods have the advantage of providing certificates of near-optimality by means of upper bounds to the optimal alignment value. We also illustrate effective heuristics, such as local search and genetic algorithms. We were able to obtain for the first time optimal alignments for large similar proteins (about 1,000 residues and 2,000 contacts) and used the CMO measure to cluster proteins in families. The clusters obtained were compared to SCOP classification in order to validate the measure. Extensive computational experiments showed that alignments which are off by at most 10% from the optimal value can be computed in a short time. Further experiments showed how this measure reacts to the choice of the threshold defining a contact and how to choose this threshold in a sensible way.     

On the biotopology of protein biosynthesis

The present paper is an attempt to outline a possible approach to the study of concrete cellular systems in terms of relational biology as developed by Rashevsky and Rosen. The basic ideas and the formalism of Rosen’s (M,R)-systems, proposed as a model of abstract biological systems, are used in order to represent the cellular protein biosynthesis. A diagram corresponding to the activation of amino acids and synthesis of amino-acyl-transfer RNA, the attachment of _t RNA to a specific codon of messenger RNA and peptide bond synthesis with the release of a protein molecule, is constructed. The systemM thus obtained for the synthesis of a proteinp _k receives a set of environmental inputs, that is, the twently naturally occurring amino acids and emits a single output, thep _k protein. The problem of noncontractibility of inputs in the system is then analyzed. In our context, it is found that the noncontractibility is not associated with the whole amino acid setS _pk but with an “essential amino acid set” , so that and represent the set of amino acids which can be replaced or absent. According to our considerations, the biochemical concept of “essential amino acid” acquires a new significance, that is, what seems “essential” is linked with the ability to form a giventRNA _t ∼a _i complex in a suitable augmented dependent set essential for the biosynthesis of a functional protein. Eventually the discussion of re-establishability leads to some important biological implications concerning the existence of ambiguous codons and the degeneracy phenomenon in the genetic code, as anecessary biochemical tool involved in adaptive processes.     

Induction of vacuolar Ca2+-ATPase and H+/Ca2+ exchange activity in yeast mutants lacking Pmr1, the Golgi Ca2+-ATPase.

We have analyzed Ca2+ transport activity in defined subcellular fractions of an isogenic set of wild-type and mutant yeast. The results, together with measurements of polypeptide expression levels and promoter::reporter gene activity, show that the Golgi Ca2+-ATPase, Pmr1, is the major Ca2+ pump under normal growth conditions. In the absence of Pmr1, we show a massive, calcineurin-dependent compensatory induction of the vacuolar Ca2+-ATPase, Pmc1. In addition, H+/Ca2+ exchange activity, that may be distinct from the vacuolar exchanger Vcx1, is also increased.     

Diffusion of transcription factors can drastically enhance the noise in gene expression

We study by Green's Function Reaction Dynamics the effect of the diffusive motion of repressor molecules on the noise in mRNA and protein levels for a gene that is under the control of a repressor. We find that spatial fluctuations due to diffusion can drastically enhance the noise in gene expression. After dissociation from the operator, a repressor can rapidly rebind to the DNA. Our results show that the rebinding trajectories are so short that, on this timescale, the RNA polymerase (RNAP) cannot effectively compete with the repressor for binding to the promoter. As a result, a dissociated repressor molecule will on average rebind many times, before it eventually diffuses away. These rebindings thus lower the effective dissociation rate, and this increases the noise in gene expression. Another consequence of the timescale separation between repressor rebinding and RNAP association is that the effect of spatial fluctuations can be described by a well-stirred, zero-dimensional, model by renormalizing the reaction rates for repressor-DNA (un) binding. Our results thus support the use of well-stirred, zero-dimensional models for describing noise in gene expression. We also show that for a fixed repressor strength, the noise due to diffusion can be minimized by increasing the number of repressors or by decreasing the rate of the open complex formation. Lastly, our results emphasize that power spectra are a highly useful tool for studying the propagation of noise through the different stages of gene expression.     

Quantitative real-time PCR study of the influence of probiotic culture soluble fraction on the expression of Pseudomonas aeruginosa quorum sensing genes

The opportunistic pathogen Pseudomonas (Ps.) aeruginosa causes severe infections, particularly in immunocompromised individuals and patients with cystic fibrosis (CF). A serious side effect of antibiotic therapy in Ps. aeruginosa infections is the development of resistance to antibiotics. During the infection process Ps. aeruginosa forms biofilms, rendering bacterial cells more resistant to disinfectants, antibiotics and the action of host immune defense effectors. Pseudomonas aeruginosa employs the intercellular communication system, known as quorum sensing (QS) to coordinate the expression of tissue-damaging factors. Since the QS systems controls the production of different virulence factors, it is possible that the inhibition of its regulatory activity to severely compromise the ability of Ps. aeruginosa to cause infections in humans. Many studies have shown that some probiotic strains exhibit inhibitory activity on different virulence properties of pathogenic bacteria (adherence to cellular or inert substrate, soluble virulence factors expression). The aim of the present study was to investigate by real-time RT-qPCR the influence of probiotic culture soluble factors on the QS genes expression in 30 Ps. aeruginosa strains isolated from patients hospitalized in the National Institute for Cardiovascular Infections, Prof. C.C. Iliescu Fundeni Hospital, Bucharest. The results of the real time RT-qPCR have shown that in all Ps. aeruginosa strains grown in the presence of probiotic culture sterile filtrates, the level of QS genes expression was reduced comparatively with those from control cultures. In conclusion, these results proved that the inhibition of virulence factors regulation mechanisms by soluble molecules secreted by probiotics could represent an interesting way pathogenicity and virulence attenuation in Ps. aeruginosa nosocomial strains.     

Babel's tower revisited: a universal resource for cross-referencing across annotation databases

MOTIVATION: Annotation databases are widely used as public repositories of biological knowledge. However, most of these resources have been developed by independent groups which used different designs and different identifiers for the same biological entities. As we show in this article, incoherent name spaces between various databases represent a serious impediment to using the existing annotations at their full potential. Navigating between various such name spaces by mapping IDs from one database to another is a very important issue which is not properly addressed at the moment. RESULTS: We have developed a web-based resource, Onto-Translate (OT), which effectively addresses this problem. OT is able to map onto each other different types of biological entities from the following annotation databases: Swiss-Prot, TrEMBL, NREF, PIR, Gene Ontology, KEGG, Entrez Gene, GenBank, GenPept, IMAGE, RefSeq, UniGene, OMIM, PDB, Eukaryotic Promoter Database, HUGO Gene Nomenclature Committee and NetAffx. Currently, OT is able to perform 462 types of mappings between 29 different types of IDs from 17 databases concerning 53 organisms. Among these, over 300 types of translations and 15 types of IDs are not currently supported by any other tool or resource. On average, OT is able to correctly map between 96 and 99% of the biological entities provided as input. In terms of speed, sets of approximately 20 000 IDs can be translated in <30 s, in most cases. AVAILABILITY: OT is a part of Onto-Tools, which is freely available at http://vortex.cs.wayne.edu/Projects.html     

Global functional profiling of gene expression

The typical result of a microarray experiment is a list of tens or hundreds of genes found to be differentially regulated in the condition under study. Independent of the methods used to select these genes, the common task faced by any researcher is to translate these lists of genes into a better understanding of the biological phenomena involved. Currently, this is done through a tedious combination of searches through the literature and a number of public databases. We developed Onto-Express (OE) as a novel tool able to automatically translate such lists of differentially regulated genes into functional profiles characterizing the impact of the condition studied. OE constructs functional profiles (using Gene Ontology terms) for the following categories: biochemical function, biological process, cellular role, cellular component, molecular function, and chromosome location. Statistical significance values are calculated for each category. We demonstrate the validity and the utility of this comprehensive global analysis of gene function by analyzing two breast cancer datasets from two separate laboratories. OE was able to identify correctly all biological processes postulated by the original authors, as well as discover novel relevant mechanisms.     

Solid state NMR sequential resonance assignments and conformational analysis of the 2×10.4 kDa dimeric form of the Bacillus subtilis protein Crh

Solid state NMR sample preparation and resonance assignments of the U-[¹³C,¹⁵N] 2×10.4 kDa dimeric form of the regulatory protein Crh in microcrystalline, PEG precipitated form are presented. Intra– and interresidue correlations using dipolar polarization transfer methods led to nearly complete sequential assignments of the protein, and to 88% of all ¹⁵N, ¹³C chemical shifts. For several residues, the resonance assignments differ significantly from those reported for the monomeric form analyzed by solution state NMR. Dihedral angles obtained from a TALOS-based statistical analysis suggest that the microcrystalline arrangement of Crh must be similar to the domain-swapped dimeric structure of a single crystal form recently solved using X-ray crystallography. For a limited number of protein residues, a remarkable doubling of the observed NMR resonances is observed indicative of local static or dynamic conformational disorder. Our study reports resonance assignments for the largest protein investigated by solid state NMR so far and describes the conformational dimeric variant of Crh with previously unknown chemical shifts.