Transcript annotation in FANTOM3: Mouse gene catalog based on physical cDNAs

T he international FANTOM consortium aims to produce a comprehensive picture of the mammalian transcriptome, based upon an extensive cDNA collection and functional annotation of full-length enriched cDNAs. The previous dataset, FANTOM(2), comprised 60,770 full- length enriched cDNAs. Functional annotation revealed that this cDNA dataset contained only about half of the estimated number of mouse protein- coding genes, indicating that a number of cDNAs still remained to be collected and identified. To pursue the complete gene catalog that covers all predicted mouse genes, cloning and sequencing of full- length enriched cDNAs has been continued since FANTOM2. In FANTOM3, 42,031 newly isolated cDNAs were subjected to functional annotation, and the annotation of 4,347 FANTOM2 cDNAs was updated. To accomplish accurate functional annotation, we improved our automated annotation pipeline by introducing new coding sequence prediction programs and developed a Web- based annotation interface for simplifying the annotation procedures to reduce manual annotation errors. Automated coding sequence and function prediction was followed with manual curation and review by expert curators. A total of 102,801 full- length enriched mouse cDNAs were annotated. Out of 102,801 transcripts, 56,722 were functionally annotated as protein coding ( including partial or truncated transcripts), providing to our knowledge the greatest current coverage of the mouse proteome by full- length cDNAs. The total number of distinct non- protein- coding transcripts increased to 34,030. The FANTOM3 annotation system, consisting of automated computational prediction, manual curation, and. nal expert curation, facilitated the comprehensive characterization of the mouse transcriptome, and could be applied to the transcriptomes of other species.

Use of tissue Doppler imaging to facilitate the prediction of events in patients with abnormal left ventricular function dobutamine echocardiography

The extent of abnormality in patients with positive do-butamine echocardiography (DE) is predictive of risk, but the wall motion score (WMS) has low concordance among observers. We sought whether quantifying the extent of abnormal wall motion using tissue Doppler (TD) could guide risk assessment in patients with abnormal DE in 576 patients with known or suspected coronary artery disease; standard DE was combined with color TD imaging at peak dose. WMS was assessed by an expert observer and studies were identified as abnormal in the presence of 2:1 segments with resting or stress-induced wall motion abnormalities. Patients with abnormal DE had peak systolic velocity measured in each segment. Tissue tracking was used to measure myocardial displacement. Follow-up for death or infarction was per-formed after. 16 +/- 12 months. Of 251 patients with abnormal DE, 22 patients died (20 from cardiac causes) and 7 had nonfatal myocardial infarctionis. The average WMS in patients with events was 1.8 +/- 0.5, compared with 1.7 +/- 0.5 in patients without events (p = NS). The average systolic velocity in patients with events was 4.9 +/- 1.7 cm/s and 6.4 +/- 6.5 cm/s in the patients without events (p <0.001). The average tissue tracking in patients with events was 4.5 +/- 1.5 mm and was significant. (5.7 +/- 3.1 mm),in those,without events (p <0.001). Thus, TD is an alternative to WMS for quantifying the total extent of abnormal left ventricular function-at DE, and appears to be superior for predicting adverse outcomes. (C) 2004 by Excerpta Medica, Inc.

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Improved slope adjustment functions for soil erosion prediction

Numerous studies in the last 60 years have investigated the relationship between land slope and soil erosion rates. However, relatively few of these have investigated slope gradient responses: ( a) for steep slopes, (b) for specific erosion processes, and ( c) as a function of soil properties. Simulated rainfall was applied in the laboratory on 16 soils and 16 overburdens at 100 mm/h to 3 replicates of unconsolidated flume plots 3 m long by 0.8 m wide and 0.15 m deep at slopes of 20, 5, 10, 15, and 30% slope in that order. Sediment delivery at each slope was measured to determine the relationship between slope steepness and erosion rate. Data from this study were evaluated alongside data and existing slope adjustment functions from more than 55 other studies from the literature. Data and the literature strongly support a logistic slope adjustment function of the form S = A + B/[1 + exp (C - D sin theta)] where S is the slope adjustment factor and A, B, C, and D are coefficients that depend on the dominant detachment and transport processes. Average coefficient values when interill-only processes are active are A - 1.50, B 6.51, C 0.94, and D 5.30 (r(2) = 0.99). When rill erosion is also potentially active, the average slope response is greater and coefficient values are A - 1.12, B 16.05, C 2.61, and D 8.32 (r(2) = 0.93). The interill-only function predicts increases in sediment delivery rates from 5 to 30% slope that are approximately double the predictions based on existing published interill functions. The rill + interill function is similar to a previously reported value. The above relationships represent a mean slope response for all soils, yet the response of individual soils varied substantially from a 2.5-fold to a 50-fold increase over the range of slopes studied. The magnitude of the slope response was found to be inversely related ( log - log linear) to the dispersed silt and clay content of the soil, and 3 slope adjustment equations are proposed that provide a better estimate of slope response when this soil property is known. Evaluation of the slope adjustment equations proposed in this paper using independent datasets showed that the new equations can improve soil erosion predictions.

Quantitative myocardial contrast echocardiography for prediction of thrombolysis in myocardial infarction flow in acute myocardial infarction

Clinical evaluation of arterial potency in acute ST-elevation myocardial infarction (STEMI) is unreliable. We sought to identify infarction and predict infarct-related artery potency measured by the Thrombolysis In Myocardial Infarction (TIMI) score with qualitative and quantitative intravenous myocardial contrast echocardiography (MCE). Thirty-four patients with suspected STEMI underwent MCE before emergency angiography and planned angioplasty. MCE was performed with harmonic imaging and variable triggering intervals during intravenous administration of Optison. Myocardial perfusion was quantified offline, fitting an exponential function to contrast intensity at various pulsing intervals. Plateau myocardial contrast intensity (A), rate of rise (beta), and myocardial flow (Q = A x beta) were assessed in 6 segments. Qualitative assessment of perfusion defects was sensitive for the diagnosis of infarction (sensitivity 93%) and did not differ between anterior and inferior infarctions. However, qualitative assessment had only moderate specificity (50%), and perfusion defects were unrelated to TIMI flow. In patients with STEMI, quantitatively derived myocardial blood flow Q (A x beta) was significantly lower in territories subtended by an artery with impaired (TIMI 0 to 2) flow than those territories supplied by a reperfused artery with TIMI 3 flow (10.2 +/- 9.1 vs 44.3 +/- 50.4, p = 0.03). Quantitative flow was also lower in segments with impaired flow in the subtending artery compared with normal patients with TIMI 3 flow (42.8 +/- 36.6, p = 0.006) and all segments with TIMI 3 flow (35.3 +/- 32.9, p = 0.018). An receiver-operator characteristic curve derived cut-off Q value of

Wurst: a protein threading server with a structural scoring function, sequence profiles and optimized substitution matrices

Wurst is a protein threading program with an emphasis on high quality sequence to structure alignments (http://www.zbh.uni-hamburg.de/wurst). Submitted sequences are aligned to each of about 3000 templates with a conventional dynamic programming algorithm, but using a score function with sophisticated structure and sequence terms. The structure terms are a log-odds probability of sequence to structure fragment compatibility, obtained from a Bayesian classification procedure. A simplex optimization was used to optimize the sequence-based terms for the goal of alignment and model quality and to balance the sequence and structural contributions against each other. Both sequence and structural terms operate with sequence profiles.

Statistical tests of load-unload response ratio signals by lattice solid model: Implication to tidal triggering and earthquake prediction

Statistical tests of Load-Unload Response Ratio (LURR) signals are carried in order to verify statistical robustness of the previous studies using the Lattice Solid Model (MORA et al., 2002b). In each case 24 groups of samples with the same macroscopic parameters (tidal perturbation amplitude A, period T and tectonic loading rate k) but different particle arrangements are employed. Results of uni-axial compression experiments show that before the normalized time of catastrophic failure, the ensemble average LURR value rises significantly, in agreement with the observations of high LURR prior to the large earthquakes. In shearing tests, two parameters are found to control the correlation between earthquake occurrence and tidal stress. One is, A/(kT) controlling the phase shift between the peak seismicity rate and the peak amplitude of the perturbation stress. With an increase of this parameter, the phase shift is found to decrease. Another parameter, AT/k, controls the height of the probability density function (Pdf) of modeled seismicity. As this parameter increases, the Pdf becomes sharper and narrower, indicating a strong triggering. Statistical studies of LURR signals in shearing tests also suggest that except in strong triggering cases, where LURR cannot be calculated due to poor data in unloading cycles, the larger events are more likely to occur in higher LURR periods than the smaller ones, supporting the LURR hypothesis.

Determinants of functional capacity in patients with chronic heart failure: Role of filling pressure and systolic and diastolic function

Background Previous work suggesting a better correlation of diastolic than systolic function with exercise capacity in heart failure may reflect the -relative insensitivity and load-dependence of ejection fraction (EF). We sought the correlation of new and more sensitive methods of quantifying systolic and diastolic function and filling pressure with functional capacity. Methods We studied 155 consecutive exercise tests on 95 patients with congestive heart failure (81 male, aged 62 +/- 10 years), who underwent resting 2-climensional echocardiography and tissue Doppler imaging before and after measurement of maximum oxygen uptake (peak VO2)Results The resting EF was 3 1 % 10% and a peak VO(2)was 13 +/- 5 mL/kg/min; the majority of these patients (80%) had an ischemic cardiornyopathy. Resting EF (r 0.14, P =.09) correlated poorly with peak VO2 and mean systolic (r = 0.23, P =.004) and diastolic tissue velocities (r 0.18, P =.02). Peak EF was weakly correlated with the mean systolic (r = 0.18, P =.02) and diastolic velocities (r = 0.16, P <.04). The mean sum of systolic and diastolic velocities in both annuli (r = 0.30, P <.001) and E/Ea ratio (r 0.31, P <.001) were better correlated with peak VO2 Prediction of peak VO2 was similar with models based on models of filling pressure (R = 0.61), systolic factors (R = 0.63), and diastolic factors (R 0.59), although a composite model of filling pressure, systolic and diastolic function was a superior predictor of peak VO2 (R 0.69; all P<.001). Conclusions The reported association of diastolic rather than systolic function with functional capacity may have reflected the limitations of EF. Functional capacity appears related not only to diastolic function, but also to systolic function and filling pressure, and is most closely associated with a combination of these factors.

Better prediction of protein contact number using a support vector regression analysis of amino acid sequence

Background: Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of C-beta atoms in other residues within a sphere around the C-beta atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence. Results: We predict contact number from protein sequence using a novel support vector regression algorithm. Using protein local sequences with multiple sequence alignments (PSI-BLAST profiles), we demonstrate a correlation coefficient between predicted and observed contact numbers of 0.70, which outperforms previously achieved accuracies. Including additional information about sequence weight and amino acid composition further improves prediction accuracies significantly with the correlation coefficient reaching 0.73. If residues are classified as being either contacted or non-contacted, the prediction accuracies are all greater than 77%, regardless of the choice of classification thresholds. Conclusion: The successful application of support vector regression to the prediction of protein contact number reported here, together with previous applications of this approach to the prediction of protein accessible surface area and B-factor profile, suggests that a support vector regression approach may be very useful for determining the structure-function relation between primary sequence and higher order consecutive protein structural and functional properties.

Accurate prediction of scorpion toxin functional properties from primary structures

Scorpion toxins are common experimental tools for studies of biochemical and pharmacological properties of ion channels. The number of functionally annotated scorpion toxins is steadily growing, but the number of identified toxin sequences is increasing at much faster pace. With an estimated 100,000 different variants, bioinformatic analysis of scorpion toxins is becoming a necessary tool for their systematic functional analysis. Here, we report a bioinformatics-driven system involving scorpion toxin structural classification, functional annotation, database technology, sequence comparison, nearest neighbour analysis, and decision rules which produces highly accurate predictions of scorpion toxin functional properties. (c) 2005 Elsevier Inc. All rights reserved.

Benefits of all work and no play: The relationship between neuroticism and performance as a function of resource allocation

The authors evaluate a model suggesting that the performance of highly neurotic individuals, relative to their stable counterparts, is more strongly influenced by factors relating to the allocation of attentional resources. First, an air traffic control simulation was used to examine the interaction between effort intensity and scores on the Anxiety subscale of Eysenck Personality Profiler Neuroticism in the prediction of task performance. Overall effort intensity enhanced performance for highly anxious individuals more so than for individuals with low anxiety. Second, a longitudinal field study was used to examine the interaction between office busyness and Eysenck Personality Inventory Neuroticism in the prediction of telesales performance. Changes in office busyness were associated with greater performance improvements for highly neurotic individuals compared with less neurotic individuals. These studies suggest that highly neurotic individuals outperform their stable counterparts in a busy work environment or if they are expending a high level of effort.

Prediction of vapor-moisture equilibriums for a wood-moisture system using a modified UNIQUAC model

A modified UNIQUAC model has been extended to describe and predict the equilibrium relative humidity and moisture content for wood. The method is validated over a range of moisture content from oven-dried state to fiber saturation point, and over a temperature range of 20-70 degrees C. Adjustable parameters and binary interaction parameters of the UNIQUAC model were estimated from experimental data for Caribbean pine and Hoop pine as well as data available in the literature. The two group-interaction parameters for the wood-moisture system were consistent with using function group contributions for H2O, -OH and -CHO. The result reconfirms that the main contributors to water adsorption in cell walls are the hydroxyl groups of the carbohydrates in cellulose and hemicelluloses. This provides some physical insight into the intermolecular force and energy between bound water and the wood material. (c) 2006 Elsevier Ltd. All rights reserved.

Prediction of HLA-DQ3.2 beta ligands: evidence of multiple registers in class II binding peptides

Motivation: While processing of MHC class II antigens for presentation to helper T-cells is essential for normal immune response, it is also implicated in the pathogenesis of autoimmune disorders and hypersensitivity reactions. Sequence-based computational techniques for predicting HLA-DQ binding peptides have encountered limited success, with few prediction techniques developed using three-dimensional models. Methods: We describe a structure-based prediction model for modeling peptide-DQ3.2 beta complexes. We have developed a rapid and accurate protocol for docking candidate peptides into the DQ3.2 beta receptor and a scoring function to discriminate binders from the background. The scoring function was rigorously trained, tested and validated using experimentally verified DQ3.2 beta binding and non-binding peptides obtained from biochemical and functional studies. Results: Our model predicts DQ3.2 beta binding peptides with high accuracy [area under the receiver operating characteristic (ROC) curve A(ROC) > 0.90], compared with experimental data. We investigated the binding patterns of DQ3.2 beta peptides and illustrate that several registers exist within a candidate binding peptide. Further analysis reveals that peptides with multiple registers occur predominantly for high-affinity binders.

SCORPION2: A database for structure-function analysis of scorpion toxins

Scorpion toxins are important experimental tools for characterization of vast array of ion channels and serve as scaffolds for drug design. General public database entries contain limited annotation whereby rich structure-function information from mutation studies is typically not available. SCORPION2 contains more than 800 records of native and mutant toxin sequences enriched with binding affinity and toxicity information, 624 three-dimensional structures and some 500 references. SCORPION2 has a set of search and prediction tools that allow users to extract and perform specific queries: text searches of scorpion toxin records, sequence similarity search, extraction of sequences, visualization of scorpion toxin structures, analysis of toxic activity, and functional annotation of previously uncharacterized scorpion toxins. The SCORPION2 database is available at http://sdmc.i2r.a-star.edu.sg/scorpion/. (c) 2006 Elsevier Ltd. All rights reserved.