Predicting residue-wise contact orders in proteins by support vector regression

BACKGROUND: The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. RESULTS: We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. CONCLUSION: The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences

Identification of an ATP/P2X7/Mast Cell Pathway Mediating Ozone-Induced Bronchial Hyperresponsiveness

Ozone is a highly reactive environmental pollutant with well-recognized adverse effects on lung health. Bronchial hyperresponsiveness (BHR) is one consequence of ozone exposure, particularly for individuals with underlying lung disease. Our data demonstrated that ozone induced substantial ATP release from human airway epithelia in vitro and into the airways of mice in vivo and that ATP served as a potent inducer of mast cell degranulation and BHR, acting through P2X7 receptors on mast cells. Both mast cell-deficient and P2X7 receptor-deficient (P2X7-/-) mice demonstrated markedly attenuated BHR to ozone. Reconstitution of mast cell-deficient mice with WT mast cells and P2X7-/- mast cells restored ozone-induced BHR. Despite equal numbers of mast cells in reconstituted mouse lungs, mice reconstituted with P2X7-/- mast cells demonstrated significantly less robust BHR than mice reconstituted with WT mast cells. These results support a model where P2X7 on mast cells and other cell types contribute to ozone-induced BHR

Shared midgut binding sites for Cry1A.105, Cry1Aa, Cry1Ab, Cry1Ac and Cry1Fa proteins from Bacillus thuringiensis in two important corn pests, Ostrinia nubilalis and Spodoptera frugiperda

First generation of insect-protected transgenic corn (Bt-corn) was based on the expression of Cry1Ab or Cry1Fa proteins. Currently, the trend is the combination of two or more genes expressing proteins that bind to different targets. In addition to broadening the spectrum of action, this strategy helps to delay the evolution of resistance in exposed insect populations. One of such examples is the combination of Cry1A.105 with Cry1Fa and Cry2Ab to control O. nubilalis and S. frugiperda. Cry1A.105 is a chimeric protein with domains I and II and the C-terminal half of the protein from Cry1Ac, and domain III almost identical to Cry1Fa. The aim of the present study was to determine whether the chimeric Cry1A.105 has shared binding sites either with Cry1A proteins, with Cry1Fa, or with both, in O. nubilalis and in S. frugiperda. Brush-border membrane vesicles (BBMV) from last instar larval midguts were used in competition binding assays with 125I-labeled Cry1A.105, Cry1Ab, and Cry1Fa, and unlabeled Cry1A.105, Cry1Aa, Cry1Ab, Cry1Ac, Cry1Fa, Cry2Ab and Cry2Ae. The results showed that Cry1A.105, Cry1Ab, Cry1Ac and Cry1Fa competed with high affinity for the same binding sites in both insect species. However, Cry2Ab and Cry2Ae did not compete for the binding sites of Cry1 proteins. Therefore, according to our results, the development of cross-resistance among Cry1Ab/Ac, Cry1A.105, and Cry1Fa proteins is possible in these two insect species if the alteration of shared binding sites occurs. Conversely, cross-resistance between these proteins and Cry2A proteins is very unlikely in such case

Phenotypic characterization and 16S rDNA identification of culturable non-obligate halophilic bacterial communities from a hypersaline lake, La Sal del Rey, in extreme South Texas (USA)

Background: La Sal del Rey ( the King’s Salt”) is one of several naturally-occurring salt lakes in Hidalgo County, Texas and is part of the Lower Rio Grande Valley National Wildlife Refuge. The research objective was to isolate and characterize halophilic microorganisms from La Sal del Rey. Water samples were collected from the lake and a small creek that feeds into the lake. Soil samples were collected from land adjacent to the water sample locations. Sample salinity was determined using a refractometer. Samples were diluted and cultured on a synthetic saline medium to grow halophilic bacteria. The density of halophiles was estimated by viable plate counts. A collection of isolates was selected, gram-stained, tested for catalase, and characterized using API 20E® test strips. Isolates were putatively identified by sequencing the 16S rDNA. Carbon source utilization by the microbial community from each sample site was examined using EcoPlate™ assays and the carbon utilization total activity of the community was determined. Results: Results showed that salinity ranged from 4 parts per thousand (ppt) at the lake water source to 420 ppt in water samples taken just along the lake shore. The density of halophilic bacteria in water samples ranged from 1.2 × 102 - 5.2 × 103 colony forming units per ml (cfu ml-1) whereas the density in soil samples ranged from 4.0 × 105 - 2.5 × 106 colony forming units per gram (cfu g-1). In general, as salinity increased the density of the bacterial community decreased. Microbial communities from water and soil samples were able to utilize 12 - 31 carbon substrates. The greatest number of substrates utilized was by water-borne communities compared to soil-based communities, especially at lower salinities. The majority of bacteria isolated were gram-negative, catalase-positive, rods. Biochemical profiles constructed from API 20E® test strips showed that bacterial isolates from low-salinity water samples (4 ppt) showed the greatest phenotypic diversity with regards to the types and number of positive tests from the strip. Isolates taken from water samples at the highest salinity (420 ppt) tended to be less diverse and have only a limited number of positive tests. Sequencing of 16S DNA displayed the presence of members of bacterial genera Bacillus, Halomonas, Pseudomonas, Exiguobacterium and others. The genus Bacillus was most commonly identified. None of the isolates were members of the Archaea probably due to dilution of salts in the samples. Conclusions: The La Sal del Rey ecosystem supports a robust and diverse bacterial community despite the high salinity of the lake and soil. However, salinity does appear to a limiting factor with

Identification of gene expression changes associated with the initiation of diapause in the brain of the cotton bollworm, Helicoverpa armigera

<p>Abstract</p> <p>Background</p> <p>Diapause, a state of arrested development accompanied by a marked decrease of metabolic rate, helps insects to overcome unfavorable seasons. <it>Helicoverpa armigera </it>(Har) undergoes pupal diapause, but the molecular mechanism of diapause initiation is unclear. Using suppression subtractive hybridization (SSH), we investigated differentially expressed genes in diapause- and nondiapause-destined pupal brains at diapause initiation.</p> <p>Results</p> <p>We constructed two SSH libraries (forward, F and reverse, R) to isolate genes that are up-regulated or down-regulated at diapause initiation. We obtained 194 unique sequences in the F library and 115 unique sequences in the R library. Further, genes expression at the mRNA and protein level in diapause- and nondiapause-destined pupal brains were confirmed by RT-PCR, Northern blot or Western blot analysis. Finally, we classified the genes and predicted their possible roles at diapause initiation.</p> <p>Conclusion</p> <p>Differentially expressed genes at pupal diapause initiation are possibly involved in the regulation of metabolism, energy, stress resistance, signaling pathways, cell cycle, transcription and translation.</p

Urinary levels of N-nitroso compounds in relation to risk of gastric cancer: Findings from the Shanghai cohort study

Background: N-Nitroso compounds are thought to play a significant role in the development of gastric cancer. Epidemiological data, however, are sparse in examining the associations between biomarkers of exposure to N-nitroso compounds and the risk of gastric cancer. Methods: A nested case-control study within a prospective cohort of 18,244 middle-aged and older men in Shanghai, China, was conducted to examine the association between urinary level of N-nitroso compounds and risk of gastric cancer. Information on demographics, usual dietary intake, and use of alcohol and tobacco was collected through in-person interviews at enrollment. Urinary levels of nitrate, nitrite, N-nitroso-2-methylthiazolidine-4-carboxylic acid (NMTCA), N-nitrosoproline (NPRO), N-nitrososarcosine (NSAR), N-nitrosothiazolidine-4-carboxylic acid (NTCA), as well as serum H. pylori antibodies were quantified in 191 gastric cancer cases and 569 individually matched controls. Logistic regression method was used to assess the association between urinary levels of N-nitroso compounds and risk of gastric cancer. Results: Compared with controls, gastric cancer patients had overall comparable levels of urinary nitrate, nitrite, and N-nitroso compounds. Among individuals seronegative for antibodies to H. pylori, elevated levels of urinary nitrate were associated with increased risk of gastric cancer. The multivariate-adjusted odds ratios for the second and third tertiles of nitrate were 3.27 (95% confidence interval = 0.76-14.04) and 4.82 (95% confidence interval = 1.05-22.17), respectively, compared with the lowest tertile (P for trend = 0.042). There was no statistically significant association between urinary levels of nitrite or N-nitroso compounds and risk of gastric cancer. Urinary NMTCA level was significantly associated with consumption of alcohol and preserved meat and fish food items. Conclusion: The present study demonstrates that exposure to nitrate, a precursor of N-nitroso compounds, may increase the risk of gastric cancer among individuals without a history of H. pylori infection

Kinome-wide interaction modelling using alignment-based and alignment-independent approaches for kinase description and linear and non-linear data analysis techniques

<p>Abstract</p> <p>Background</p> <p>Protein kinases play crucial roles in cell growth, differentiation, and apoptosis. Abnormal function of protein kinases can lead to many serious diseases, such as cancer. Kinase inhibitors have potential for treatment of these diseases. However, current inhibitors interact with a broad variety of kinases and interfere with multiple vital cellular processes, which causes toxic effects. Bioinformatics approaches that can predict inhibitor-kinase interactions from the chemical properties of the inhibitors and the kinase macromolecules might aid in design of more selective therapeutic agents, that show better efficacy and lower toxicity.</p> <p>Results</p> <p>We applied proteochemometric modelling to correlate the properties of 317 wild-type and mutated kinases and 38 inhibitors (12,046 inhibitor-kinase combinations) to the respective combination's interaction dissociation constant (K_d). We compared six approaches for description of protein kinases and several linear and non-linear correlation methods. The best performing models encoded kinase sequences with amino acid physico-chemical z-scale descriptors and used support vector machines or partial least- squares projections to latent structures for the correlations. Modelling performance was estimated by double cross-validation. The best models showed high predictive ability; the squared correlation coefficient for new kinase-inhibitor pairs ranging P²= 0.67-0.73; for new kinases it ranged P²_kin= 0.65-0.70. Models could also separate interacting from non-interacting inhibitor-kinase pairs with high sensitivity and specificity; the areas under the ROC curves ranging AUC = 0.92-0.93. We also investigated the relationship between the number of protein kinases in the dataset and the modelling results. Using only 10% of all data still a valid model was obtained with P²= 0.47, P²_kin= 0.42 and AUC = 0.83.</p> <p>Conclusions</p> <p>Our results strongly support the applicability of proteochemometrics for kinome-wide interaction modelling. Proteochemometrics might be used to speed-up identification and optimization of protein kinase targeted and multi-targeted inhibitors.</p

Enhancement Effects of Martentoxin on Glioma BK Channel and BK Channel (α+β1) Subtypes

BACKGROUND: BK channels are usually activated by membrane depolarization and cytoplasmic Ca(2+). Especially,the activity of BK channel (α+β4) can be modulated by martentoxin, a 37 residues peptide, with Ca(2+)-dependent manner. gBK channel (glioma BK channel) and BK channel (α+β1) possessed higher Ca(2+) sensitivity than other known BK channel subtypes. METHODOLOGY AND PRINCIPAL FINDINGS: The present study investigated the modulatory characteristics of martentoxin on these two BK channel subtypes by electrophysiological recordings, cell proliferation and Ca(2+) imaging. In the presence of cytoplasmic Ca(2+), martentoxin could enhance the activities of both gBK and BK channel (α+β1) subtypes in dose-dependent manner with EC(50) of 46.7 nM and 495 nM respectively, while not shift the steady-state activation of these channels. The enhancement ratio of martentoxin on gBK and BK channel (α+β1) was unrelated to the quantitative change of cytoplasmic Ca(2+) concentrations though the interaction between martentoxin and BK channel (α+β1) was accelerated under higher cytoplasmic Ca(2+). The selective BK pore blocker iberiotoxin could fully abolish the enhancement of these two BK subtypes induced by martentoxin, suggesting that the auxiliary β subunit might contribute to the docking for martentoxin. However, in the absence of cytoplasmic Ca(2+), the activity of gBK channel would be surprisingly inhibited by martentoxin while BK channel (α+β1) couldn't be affected by the toxin. CONCLUSIONS AND SIGNIFICANCE: Thus, the results shown here provide the novel evidence that martentoxin could increase the two Ca(2+)-hypersensitive BK channel subtypes activities in a new manner and indicate that β subunit of these BK channels plays a vital role in this enhancement by martentoxin

Down Regulation of a Gene for Cadherin, but Not Alkaline Phosphatase, Associated with Cry1Ab Resistance in the Sugarcane Borer Diatraea saccharalis

The sugarcane borer, Diatraea saccharalis, is a major target pest of transgenic corn expressing Bacillus thuringiensis (Bt) proteins (i.e., Cry1Ab) in South America and the mid-southern region of the United States. Evolution of insecticide resistance in such target pests is a major threat to the durability of transgenic Bt crops. Understanding the pests' resistance mechanisms will facilitate development of effective strategies for delaying or countering resistance. Alterations in expression of cadherin- and alkaline phosphatase (ALP) have been associated with Bt resistance in several species of pest insects. In this study, neither the activity nor gene regulation of ALP was associated with Cry1Ab resistance in D. saccharalis. Total ALP enzymatic activity was similar between Cry1Ab-susceptible (Cry1Ab-SS) and -resistant (Cry1Ab-RR) strains of D. saccharalis. In addition, expression levels of three ALP genes were also similar between Cry1Ab-SS and -RR, and cDNA sequences did not differ between susceptible and resistant larvae. In contrast, altered expression of a midgut cadherin (DsCAD1) was associated with the Cry1Ab resistance. Whereas cDNA sequences of DsCAD1 were identical between the two strains, the transcript abundance of DsCAD1 was significantly lower in Cry1Ab-RR. To verify the involvement of DsCAD1 in susceptibility to Cry1Ab, RNA interference (RNAi) was employed to knock-down DsCAD1 expression in the susceptible larvae. Down-regulation of DsCAD1 expression by RNAi was functionally correlated with a decrease in Cry1Ab susceptibility. These results suggest that down-regulation of DsCAD1 is associated with resistance to Cry1Ab in D. saccharalis

Utilization patterns of Chinese medicine and Western medicine under the National Health Insurance Program in Taiwan, a population-based study from 1997 to 2003

<p>Abstract</p> <p>Background</p> <p>In 1995, Taiwan has launched a national health-care system (the National Health Insurance Program, NHI) covering the use of both Western medicine (WM) and Chinese medicine (CM). This population-based study was conducted to understand the role of CM in this dual medical system by determining the utilization patterns of CM and WM and to analyze the demographic characteristics and primary indications influencing the choice of the medical services for the development of strategies to enhance the appropriate use and reduce unnecessary use of CM.</p> <p>Methods</p> <p>This study used the NHI sample files from 1997 to 2003 consisting of comprehensive utilization and enrolment information for a random sample of 200,432 NHI beneficiaries of the total enrolees from 1995 to 2000. A total of 136,720 subjects with valid and complete enrolment and utilization data were included in this study. The logistic regression method was employed to estimate the odds ratios (ORs) for utilization of CM and WM. The usage, frequency of services, and primary indications for CM and WM were evaluated. A significance level of α = 0.05 was selected.</p> <p>Results</p> <p>Compared with WM, the odds of CM increased from 1997 to 2003. The odds of using CM (OR = 1.48; 95% CI: 1.45–1.50; p < 0.001) and WM (OR = 1.74; 95% CI: 1.72–1.77; p < 0.001) were higher in females and that of CM increased with age to a peak in the 45–54-year-group (OR = 1.75; 95% CI: 1.68–1.82; p < 0.001) and WM (OR = 1.09; 95% CI: 1.05–1.13; p < 0.001) in the elderly subjects (≥ 65 years). The odds of CM and WM were similar in all income groups. However, those of CM were higher in Central (OR = 1.65; 95% CI: 1.56–1.74; p < 0.001) and Southern Taiwan (OR = 1.18; 95% CI: 1.12–1.25; p < 0.001) and lower in the remote areas (OR = 0.57; 95% CI: 0.52–0.63; p < 0.001). Most of the patients had one ambulatory visit of both medical services annually. However, the utilization of WM predominated over CM. Over 90% of CM service was provided by clinics, whereas over 60% of WM service by hospitals. Diseases of the respiratory system was the most frequent primary indication in CM and WM. Herbal medication was the most commonly used form of CM (68.4–72.7%).</p> <p>Conclusion</p> <p>In recent years, there is an increasing trend in the utilization of CM in Taiwan. This increasing trend may be due to the covering of CM in the national health insurance system.</p