Accurate prediction of protein secondary structure and solvent accessibility by consensus combiners of sequence and structure information

Background  

Structural properties of proteins such as secondary structure and solvent accessibility contribute to three-dimensional structure prediction, not only in the ab initio case but also when homology information to known structures is available. Structural properties are also routinely used in protein analysis even when homology is available, largely because homology modelling is lower throughput than, say, secondary structure prediction. Nonetheless, predictors of secondary structure and solvent accessibility are virtually always ab initio.     

Cross-adaptation and molecular modeling study of receptor mechanisms common to four taste stimuli in humans

Psychophysical cross-adaptation experiments were performed with two carbohydrates, sucrose (SUC) and fructose (FRU), and two sweeteners, acesulfame-K (MOD) and dulcin (DUL). Seven subjects were asked to match concentrations that elicited the same intensity as a sucrose reference (30 g/l). Cross-adaptation levels were calculated as the ratio of isointense concentrations measured for a given stimulus before and under adaptation. On average, cross-adaptation between SUC and FRU is low and apparently reciprocal. By contrast, cross-adaptation between SUC and MOD is clearly non-reciprocal: SUC adapts MOD significantly (24%, P < 0.005), but MOD fails to adapt SUC (2%, P < 0.79). Significant and reciprocal cross-enhancement is observed between DUL and MOD (approximately -20%, P < 0.03), and also between SUC and DUL (approximately -15%, P < 0.08). In parallel, molecular modeling of the four tastants was performed in order to look for the 12 common binding motifs that were isolated on 14 other tastants in a previous study. SUC and FRU each display 10 out of the 12 binding motifs, whereas DUL and MOD only display four and five distinct motifs respectively and do not have any motif in common. Experimental cross-adaptation levels seem to correlate well with the number of motifs that molecules have in common. FRU and SUC share a majority of binding motifs and correlatively show mutual cross-adaptation. Four motifs of MOD are found among the 10 motifs of SUC, which may explain why SUC cross-adapts MOD but not vice versa. By contrast, DUL and MOD do not share any motif and do not cross- adapt. The various molecular mechanisms that may be responsible for cross-adaptation and/or cross-enhancement are discussed in light of our results.      

Kinetics of thin filament activation probed by fluorescence of N-((2-(Iodoacetoxy)ethyl)-N-methyl)amino-7-nitrobenz-2-oxa-1, 3-diazole-labeled troponin I incorporated into skinned fibers of rabbit psoas muscle: implications for regulation of muscle contraction

Making use of troponin with fluorescently labeled troponin I subunit (N-((2-(iodoacetoxy)ethyl)-N-methyl)amino-7-nitrobenz-2-oxa-1, 3-diazole-troponin I, IANBD-TnI) that had previously been described in solution studies as a probe for thin filament activation (. Proc. Natl. Acad. Sci. 77:7209-7213), we present a new approach that allows the kinetics of thin filament activation to be studied in skinned muscle fibers. After the exchange of native troponin for fluorescently labeled troponin, the fluorescence intensity is sensitive to both changes in calcium concentration and actin attachment of cross-bridges in their strong binding states (. Biophys. J. 77:000-000). Imposing rapid changes in the fraction of strongly attached cross-bridges, e.g., by switching from isometric contraction to high-speed shortening, causes changes in thin filament activation at fixed Ca(2+) concentrations that can be followed by recording fluorescence intensity. Upon changing to high-speed shortening we observed small (<20%) changes in fluorescence that became faster at higher Ca(2+) concentrations. At all Ca(2+) concentrations, these changes are more than 10-fold faster than force redevelopment subsequent to the period of unloaded shortening. We interpret this as an indication that equilibration among different states of the thin filament is rapid and becomes faster as Ca(2+) is raised. Fast equilibration suggests that the rate constant of force redevelopment is not limited by changes in the activation level of thin filaments induced by the isotonic contraction before force redevelopment. Instead, our modeling shows that, in agreement with our previous proposal for the regulation of muscle contraction, a rapid and Ca(2+)-dependent equilibration among different states of the thin filament can fully account for the Ca(2+) dependence of force redevelopment and the fluorescence changes described in this study.     

Single cell studies of the cell cycle and some models

Analysis of growth and division often involves measurements made on cell populations, which tend to average data. The value of single cell analysis needs to be appreciated, and models based on findings from single cells should be taken into greater consideration in our understanding of the way in which cell size and division are co-ordinated. Examples are given of some single cell analyses in mammalian cells, yeast and other microorganisms. There is also a short discussion on how far the results are in accord with simple models.     

Two new species of rhabdochonids (Nematoda: Rhabdochonidae) from freshwater fishes in Mexico, with a description of a new genus

Two new nematode species, Beaninema nayaritense n. gen., n. sp. and Rhabdochona xiphophori n. sp., are described on the basis of the specimens recovered from the gall bladder and intestine of 2 fishes, Cichlasoma beani (Jordan) (Cichlidae, Perciformes) and Xiphophorus sp. (Poeciliidae, Cyprinodontiformes), respectively, from Mexico. The monotypic genus Beaninema differs from other rhabdochonid genera mainly in the presence of large conical teeth in the middle of the prostom. Rhabdochona xiphophori is characterized mainly by a unique structure of the prostom (presence of 6 anterior teeth; dorsal and ventral teeth unusually broad, with 2 lateral horns) and the shape of the deirids (bifurcated, with markedly long base).     

An approach to large scale identification of non-obvious structural similarities between proteins

Background

A new sequence independent bioinformatics approach allowing genome-wide search for proteins with similar three dimensional structures has been developed. By utilizing the numerical output of the sequence threading it establishes putative non-obvious structural similarities between proteins. When applied to the testing set of proteins with known three dimensional structures the developed approach was able to recognize structurally similar proteins with high accuracy.

Results

The method has been developed to identify pathogenic proteins with low sequence identity and high structural similarity to host analogues. Such protein structure relationships would be hypothesized to arise through convergent evolution or through ancient horizontal gene transfer events, now undetectable using current sequence alignment techniques. The pathogen proteins, which could mimic or interfere with host activities, would represent candidate virulence factors.The developed approach utilizes the numerical outputs from the sequence-structure threading. It identifies the potential structural similarity between a pair of proteins by correlating the threading scores of the corresponding two primary sequences against the library of the standard folds. This approach allowed up to 64% sensitivity and 99.9% specificity in distinguishing protein pairs with high structural similarity.

Conclusion

Preliminary results obtained by comparison of the genomes of Homo sapiens and several strains of Chlamydia trachomatis have demonstrated the potential usefulness of the method in the identification of bacterial proteins with known or potential roles in virulence.

     

A GWAS assessment of the contribution of genomic imprinting to the variation of body mass index in mice

Background

Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis.

Results

10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10⁻⁴ to 1.842 ×10⁻⁴ when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered.

Conclusions

When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).     

Defining window-boundaries for genomic analyses using smoothing spline techniques

Background

High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome.

Results

Simulations applying this method were performed to identify selection signatures from pooled sequencing F_ST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing F_ST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach.

Conclusions

We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on F_ST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.