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.

Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.

A comparison of computer-based methods for the determination of onset of muscle contraction using electromyography

Little consensus exists in the literature regarding methods for determination of the onset of electromyographic (EMG) activity. The aim of this study was to compare the relative accuracy of a range of computer-based techniques with respect to EMG onset determined visually by an experienced examiner. Twenty-seven methods were compared which varied in terms of EMG processing (low pass filtering at 10, 50 and 500 Hz), threshold value (1, 2 and 3 SD beyond mean of baseline activity) and the number of samples for which the mean must exceed the defined threshold (20, 50 and 100 ms). Three hundred randomly selected trials of a postural task were evaluated using each technique. The visual determination of EMG onset was found to be highly repeatable between days. Linear regression equations were calculated for the values selected by each computer method which indicated that the onset values selected by the majority of the parameter combinations deviated significantly from the visually derived onset values. Several methods accurately selected the time of onset of EMG activity and are recommended for future use. Copyright (C) 1996 Elsevier Science Ireland Ltd.

Carotid artery intimal medial thickness, brachial artery flow-mediated vasodilation and cardiovascular risk factors in diabetic and non-diabetic indigenous Australians

Background: Indigenous Australians are at high risk for cardiovascular disease and type 2 diabetes. Carotid artery intimal medial thickness (CIMT) and brachial artery flow-mediated vasodilation (FMD) are ultrasound imaging based surrogate markers of cardiovascular risk. This study examines the relative contributions of traditional cardiovascular risk factors on CIMT and FMD in adult Indigenous Australians with and without type 2 diabetes mellitus. Method: One hundred and nineteen Indigenous Australians were recruited. Physical and biochemical markers of cardiovascular risk, together with CIMT and FMD were meausred for all subjects. Results: Fifty-three Indigenous Australians subjects (45%) had type 2 diabetes mellitus. There was a significantly greater mean CIMT in diabetic versus non-diabetic subjects (p = 0.049). In the non-diabetic group with non-parametric analyses, there were significant correlations between CIMT and: age (r = 0.64, p < 0.001), systolic blood pressure (r = 0.47, p < 0.001) and non-smokers (r = -0.30, p = 0.018). In the diabetic group, non-parametric analysis showed correlations between CIMT, age (r = 0.36, p = 0.009) and duration of diabetes (r = 0.30, p = 0.035) only. Adjusting forage, sex, smoking and history of cardiovascular disease, Hb(A1c) became the sole significant correlate of CIMT (r = 0.35,p = 0.01) in the diabetic group. In non-parametric analysis, age was the sole significant correlate of FMD (r = -0.31,p = 0.013), and only in non-diabetic subjects. Linear regression analysis showed significant associations between CIMT and age (t = 4.6,p < 0.001), systolic blood pressure (t = 2.6, p = 0.010) and Hb(A1c) (t = 2.6, p = 0.012), smoking (t = 2.1, p = 0.04) and fasting LDL-cholesterol (t = 2.1, p = 0.04). There were no significant associations between FMD and examined cardiovascular risk factors with linear regression analysis Conclusions: CIMT appears to be a useful surrogate marker of cardiovascular risk in this sample of Indigenous Australian subjects, correlating better than FMD with established cardiovascular risk factors. A lifestyle intervention programme may alleviate the burden of cardiovascular disease in Indigenous Australians by reducing central obesity, lowering blood pressure, correcting dyslipidaemia and improving glycaemic control. CIMT may prove to be a useful tool to assess efficacy of such an intervention programme. (c) 2004 Elsevier Ireland Ltd. All rights reserved.

Effects of rock properties on specific cutting energy in linear cutting of sandstones by picks

Specific cutting energy (SE) has been widely used to assess the rock cuttability for mechanical excavation purposes. Some prediction models were developed for SE through correlating rock properties with SE values. However, some of the textural and compositional rock parameters i.e. texture coefficient and feldspar, mafic, and felsic mineral contents were not considered. The present study is to investigate the effects of previously ignored rock parameters along with engineering rock properties on SE. Mineralogical and petrographic analyses, rock mechanics, and linear rock cutting tests were performed on sandstone samples taken from sites around Ankara, Turkey. Relationships between SE and rock properties were evaluated using bivariate correlation and linear regression analyses. The tests and subsequent analyses revealed that the texture coefficient and feldspar content of sandstones affected rock cuttability, evidenced by significant correlations between these parameters and SE at a 90% confidence level. Felsic and mafic mineral contents of sandstones did not exhibit any statistically significant correlation against SE. Cementation coefficient, effective porosity, and pore volume had good correlations against SE. Poisson's ratio, Brazilian tensile strength, Shore scleroscope hardness, Schmidt hammer hardness, dry density, and point load strength index showed very strong linear correlations against SE at confidence levels of 95% and above, all of which were also found suitable to be used in predicting SE individually, depending on the results of regression analysis, ANOVA, Student's t-tests, and R2 values. Poisson's ratio exhibited the highest correlation with SE and seemed to be the most reliable SE prediction tool in sandstones.

Effects of rock properties on specific cutting energy in linear cutting of sandstones by picks

Specific cutting energy (SE) has been widely used to assess the rock cuttability for mechanical excavation purposes. Some prediction models were developed for SE through correlating rock properties with SE values. However, some of the textural and compositional rock parameters i.e. texture coefficient and feldspar, mafic, and felsic mineral contents were not considered. The present study is to investigate the effects of previously ignored rock parameters along with engineering rock properties on SE. Mineralogical and petrographic analyses, rock mechanics, and linear rock cutting tests were performed on sandstone samples taken from sites around Ankara, Turkey. Relationships between SE and rock properties were evaluated using bivariate correlation and linear regression analyses. The tests and subsequent analyses revealed that the texture coefficient and feldspar content of sandstones affected rock cuttability, evidenced by significant correlations between these parameters and SE at a 90% confidence level. Felsic and mafic mineral contents of sandstones did not exhibit any statistically significant correlation against SE. Cementation coefficient, effective porosity, and pore volume had good correlations against SE. Poisson's ratio, Brazilian tensile strength, Shore scleroscope hardness, Schmidt hammer hardness, dry density, and point load strength index showed very strong linear correlations against SE at confidence levels of 95% and above, all of which were also found suitable to be used in predicting SE individually, depending on the results of regression analysis, ANOVA, Student's t-tests, and R-2 values. Poisson's ratio exhibited the highest correlation with SE and seemed to be the most reliable SE prediction tool in sandstones.

Predicting the solvent accessibility of transmembrane residues from protein sequence

In this study, we propose a novel method to predict the solvent accessible surface areas of transmembrane residues. For both transmembrane alpha-helix and beta-barrel residues, the correlation coefficients between the predicted and observed accessible surface areas are around 0.65. On the basis of predicted accessible surface areas, residues exposed to the lipid environment or buried inside a protein can be identified by using certain cutoff thresholds. We have extensively examined our approach based on different definitions of accessible surface areas and a variety of sets of control parameters. Given that experimentally determining the structures of membrane proteins is very difficult and membrane proteins are actually abundant in nature, our approach is useful for theoretically modeling membrane protein tertiary structures, particularly for modeling the assembly of transmembrane domains. This approach can be used to annotate the membrane proteins in proteomes to provide extra structural and functional information.

The water balance on sloping land in rainfed lowland rice ecosystem

Quantifying water losses in paddy fields assists estimation of water availability in rainfed lowland rice ecosystem. There is, however, no definite method for determining the water losses, and little information is available on water balance in different toposequence positions of a sloped rainfed lowland. Therefore, the aim of this work was to quantify percolation and the lateral water flow with special reference to the toposequential variation. Data used for the analysis was collected in Laos and northeast Thailand. Percolation and water tables were measured on a daily basis. The percolator is a steel cylindrical tube with a lid to prevent water loss from evapotranspiration. The water table meter is a short PVC tube for determination of perched water table and a long PVC tube for groundwater table, and the side is perforated with 5-mm diameter holes at 20-mm distance. Percolation rate was determined using linear regression analysis of cumulative percolation. Assuming that the total amount of evaporation and transpiration was equivalent to potential evapotranspiration, the lateral water flow was estimated using the water balance equation. Our results are in agreement with the previously reported findings, and the methodology of estimating water balance components appears reasonably acceptable. With regard to the toposequential variation, the higher position in the topoesquence, the greater potential of the water losses because of higher percolation and lateral flow rates.

Dependence of public support for survival of wildlife species on their likeability

We surveyed a sample of 204 individuals selected from the public in Brisbane, Australia, to ascertain the extent to which they like or dislike 24 species of wildlife present in tropical Australia. The species belong to three classes: mammals, birds and reptiles. We calculated likeability indices for each of these species. We also asked respondents if they favoured the survival of each of these species and so the percentage of respondents favouring survival of each of these species could be calculated. Thus, using linear regression analysis, the percentage of respondents favouring survival of each of the species was related to their indices of likeability. In addition, the data enables the average likeability of species in the three classes (mammals, birds and reptiles) to be compared with the average support for survival of species in each of these three classes. As a result, we are able to assess how important stated likeability seems to be for preferences for survival of species, and to reconsider the hypothesis in the literature that there is likely to be more public support for the survival of mammals than for birds than for reptiles.

Increased hepatic iron concentration in nonalcoholic steatohepatitis is associated with increased fibrosis

Background & Aims: Nonalcoholic steatohepatitis (NASH) is a chronic liver disease that occasionally progresses to cirrhosis but usually has a benign course. The aim of this study was to investigate the role of the hemochromatosis mutation Cys282Tyr in development of the mild hepatic iron overload found in some patients with NASH and its association with hepatic damage in these patients. Methods: Fifty-one patients with NASH were studied. The presence of the Cys282Tyr mutation was tested in all patients, and the data were analyzed with respect to the histological grade of steatosis, inflammation, Perls' staining, hepatic iron concentration (HIC), and serum iron indices. Results: Thirty-one percent of patients with NASH were either homozygous or heterozygous for the Cys282Tyr mutation. This mutation was significantly associated with Perls' stain grade (P < 0.005), HIC (P < 0.005), and transferrin saturation percentage (P < 0.005) but not with serum ferritin levels. Linear regression analysis showed that increased hepatic iron (Perls' stain or HIC) had the greatest association with the severity of fibrosis (P < 0.0001). Conclusions: The Cys282Tyr mutation is responsible for most of the mild iron overload found in NASH and thus has a significant association with hepatic damage in these patients. Heterozygosity for the hemochromatosis gene mutation therefore cannot always be considered benign.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

Relative dispersions of intra-albumin transit times across rat and elasmobranch perfused livers, and implications for intra- and inter-species scaling of hepatic clearance using microsomal data

It is recognized that vascular dispersion in the liver is a determinant of high first-pass extraction of solutes by that organ. Such dispersion is also required for translation of in-vitro microsomal activity into in-vivo predictions of hepatic extraction for any solute. We therefore investigated the relative dispersion of albumin transit times (CV2) in the livers of adult and weanling rats and in elasmobranch livers. The mean and normalized variance of the hepatic transit time distribution of albumin was estimated using parametric non-linear regression (with a correction for catheter influence) after an impulse (bolus) input of labelled albumin into a single-pass liver perfusion. The mean +/- s.e. of CV2 for albumin determined in each of the liver groups were 0.85 +/- 0.20 (n = 12), 1.48 +/- 0.33 (n = 7) and 0.90 +/- 0.18 (n = 4) for the livers of adult and weanling rats and elasmobranch livers, respectively. These CV2 are comparable with that reported previously for the dog and suggest that the CV2 Of the liver is of a similar order of magnitude irrespective of the age and morphological development of the species. It might, therefore, be justified, in the absence of other information, to predict the hepatic clearances and availabilities of highly extracted solutes by scaling within and between species livers using hepatic elimination models such as the dispersion model with a CV2 of approximately unity.

Initial results of ultrasound screening for aneurysm of the abdominal aorta in Western Australia: relevance for endoluminal treatment of aneurysm disease

Background. Increased life expectancy in men during the last thirty years is largely due to the decrease in mortality from cardiovascular disease in the age group 29-69 yr. This change has resulted in a change in the disease profile of the population with conditions such as aneurysm of the abdominal aorta (AAA) becoming more prevalent. The advent of endoluminal treatment for AAA has encouraged prophylactic intervention and fuelled the argument to screen for the disease. The feasibility of inserting an endoluminal graft is dependent on the morphology and growth characteristics of the aneurysm. This study used data from a randomized controlled trial of ultrasound screening for AAA in men aged 65-83 yr in Western Australia for the purpose of determining the norms of the living anatomy in the pressurized infrarenal aorta. Aims. To examine (1) the diameters of the infra-renal aorta in aneurysmal and non-aneurysmal cases, (2) the implications for treatment modalities, with particular reference to endoluminal grafting, which is most dependent on normal and aneurysmal morphology, and (3) any evidence to support the notion that northern Europeans are predisposed to aneurysmal disease. Methods. Using ultrasound, a randomized control trial was established in Western Australia to assess the value of a screening program in males aged 65-83 yr, The infra-renal aorta was defined as aneurysmal if the maximum diameter was 30 mm or more. Aortic diameter was modelled both as a continuous tin mm) and as a binary outcome variable, for those men who had an infra-renal diameter of 30 mm or more. ANOVA and linear regression were used for modelling aortic diameter as a continuum, while chi-square analysis and logistic regression were used in comparing men with and without the diagnosis of AAA. Findings. By December 1998, of 19.583 men had been invited to undergo ultrasound screening for AAA, 12.203 accepted the invitation (corrected response fraction 70.8%). The prevalence of AAA increased with age from 4.8% at 65 yr to 10.8% at 80 yr (chi (2) = 77.9, df = 3, P<0.001). The median (IQR) diameter for the non-aneurysmal group was 21.4 mm (3.3 mm) and there was an increase (<chi>(2) = 76.0, df = 1, P<0.001) in the diameter of the infra-renal aorta with age. Since 27 mm is the 95th centile for the non-aneurysmal infra-renal aorta, a diameter of 30 mm or more is justified as defining an aneurysm. The risk of AAA was higher in men of Australian (OR = 1.0) and northern European origin (OR = 1.0, 95%CL: 0.9. 1.2) compared with those of Mediterranean origin (OR = 0.5, 99%CL: 0.4, 0.7). Conclusion. Although screening has not yet been shown to reduce mortality from AAA. these population-based data assist the understanding of aneurysmal disease and the further development and use of endoluminal grafts for this condition. (C) 2001 Published by Elsevier Science Ltd on behalf of The International Society for Cardiovascular Surgery.

Stochastic structure and individual-tree growth models

The majority of past and current individual-tree growth modelling methodologies have failed to characterise and incorporate structured stochastic components. Rather, they have relied on deterministic predictions or have added an unstructured random component to predictions. In particular, spatial stochastic structure has been neglected, despite being present in most applications of individual-tree growth models. Spatial stochastic structure (also called spatial dependence or spatial autocorrelation) eventuates when spatial influences such as competition and micro-site effects are not fully captured in models. Temporal stochastic structure (also called temporal dependence or temporal autocorrelation) eventuates when a sequence of measurements is taken on an individual-tree over time, and variables explaining temporal variation in these measurements are not included in the model. Nested stochastic structure eventuates when measurements are combined across sampling units and differences among the sampling units are not fully captured in the model. This review examines spatial, temporal, and nested stochastic structure and instances where each has been characterised in the forest biometry and statistical literature. Methodologies for incorporating stochastic structure in growth model estimation and prediction are described. Benefits from incorporation of stochastic structure include valid statistical inference, improved estimation efficiency, and more realistic and theoretically sound predictions. It is proposed in this review that individual-tree modelling methodologies need to characterise and include structured stochasticity. Possibilities for future research are discussed. (C) 2001 Elsevier Science B.V. All rights reserved.

Weight loss during prolonged lactation in rural Bangladeshi mothers

To determine the duration of lactation which is associated with weight loss in rural Bangladeshi mothers and also to determine the relationship with consumption patterns of principal food items, a cross-sectional study was carried out among 791 lactating rural Bangladeshi mothers aged 18-40 years. Results were compared with 333 non-pregnant and non-lactating mothers of a similar age group. The duration of lactation was up to 60 months. The mean difference in body-weight and body mass index (BMI) of lactating mothers who breastfed their children up to 24 months was significantly lower compared to non-lactating mothers of the same age group, but no differences were observed for those who breastfed beyond 24 months. The frequency of consumption of principal food items was comparable between the non-lactating and the lactating mothers who breastfed beyond 24 months. Results of multiple linear regression analysis showed that body-weight of mothers was negatively correlated with 1-12 month(s) and 13-24 months of lactation after controlling for height, education, and food consumption (slope -1.04, p < 0.05 and slope -1.23, p < 0.05 respectively). Height and consumption of meat and milk were significantly positively correlated with body-weight (slope 0.53, p < 0.001; slope 1.44, p < 0.001; and slope 0.75, p < 0.05 respectively). The study concluded that Bangladeshi women who breastfed up to 24 months were of lower weight than non-lactating mothers, most likely due to the effect of lactation. These mothers were not taking any additional foods during their lactating period. Based on the findings of the study, it is recommended that mothers consume additional energy-rich foods during the first 24 months of lactation to prevent weight loss.