Identification and analysis of chromodomain-containing proteins encoded in the mouse transcriptome

The chromodomain is 40-50 amino acids in length and is conserved in a wide range of chromatic and regulatory proteins involved in chromatin remodeling. Chromodomain-containing proteins can be classified into families based on their broader characteristics, in particular the presence of other types of domains, and which correlate with different subclasses of the chromodomains themselves. Hidden Markov model (HMM)-generated profiles of different subclasses of chromodomains were used here to identify sequences encoding chromodomain-containing proteins in the mouse transcriptome and genome. A total of 36 different loci encoding proteins containing chromodomains, including 17 novel loci, were identified. Six of these loci (including three apparent pseudogenes, a novel HP1 ortholog, and two novel Msl-3 transcription factor-like proteins) are not present in the human genome, whereas the human genome contains four loci (two CDY orthologs and two apparent CDY pseuclogenes) that are not present in mouse. A number of these loci exhibit alternative splicing to produce different isoforms, including 43 novel variants, some of which lack the chromodomain. The likely functions of these proteins are discussed in relation to the known functions of other chromodomain-containing proteins within the same family.

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.

Cost-effectiveness of combination peginterferon alpha-2a and ribavirin compared with interferon alpha-2b and ribavirin in patients with chronic hepatitis C

BACKGROUND: Sustained virological response (SVR) is the primary objective in the treatment of chronic hepatitis C (CHC). Results from a recent clinical trial of patients with previously untreated CHC demonstrate that the combination of peginterferon alpha-2a and ribavirin produces a greater SVR than interferon alpha-2b and ribavirin combination therapy. However, the cost-effectiveness of peginterferon alpha-2a plus ribavirin in the U.S. setting has not been investigated. METHODS: A Markov model was developed to investigate cost-effectiveness in patients with CHC using genotype to guide treatment duration. SVR and disease progression parameters were derived from the clinical trials and epidemiologic studies. The impact of treatment on life expectancy and costs were projected for a lifetime. Patients who had an SVR were assumed to remain virus-free for the rest of their lives. In genotype 1 patients, the SVRs were 46% for peginterferon alpha-2a plus ribavirin and 36% for interferon alpha-2b plus ribavirin. In genotype 2/3 patients, the SVRs were 76% for peginterferon alpha-2a plus ribavirin and 61% for interferon alpha-2b plus ribavirin. Quality of life and costs were based on estimates from the literature. All costs were based on published U.S. medical care costs and were adjusted to 2003 U.S. dollars. Costs and benefits beyond the first year were discounted at 3%. RESULTS: In genotype 1, peginterferon alpha-2a plus ribavirin increases quality-adjusted life expectancy (QALY) by 0.70 yr compared to interferon alpha-2b plus ribavirin, producing a cost-effectiveness ratio of $2,600 per QALY gained. In genotype 2/3 patients, peginterferon alpha-2a plus ribavirin increases QALY by 1.05 yr in comparison to interferon alpha-2b plus ribavirin. Peginterferon alpha-2a combination therapy in patients with HCV genotype 2 or 3 is dominant (more effective and cost saving) compared to interferon alpha-2b plus ribavirin. Results weighted by genotype prevalence (75% genotype 1; 25% genotype 2 or 3) also show that peginterferon alpha-2a plus ribavirin is dominant. Peginterferon alpha-2a and ribavirin remained cost-effective (below $16,500 per QALY gained) under sensitivity analyses on key clinical and cost parameters. CONCLUSION: Peginterferon alpha-2a in combination with ribavirin with duration of therapy based on genotype, is cost-effective compared with conventional interferon alpha-2b in combination with ribavirin when given to treatment-naive adults with CHC.

A cost-effectiveness analysis of fluvastatin in patients with diabetes after successful percutaneous coronary intervention

Background: The Lescol Intervention Prevention Study (LIPS) was a multinational randomized controlled trial that showed a 47% reduction in the relative risk of cardiac death and a 22% reduction in major adverse cardiac events (MACEs) from the routine use of fluvastatin, compared with controls, in patients undergoing percutaneous coronary intervention (PCI, defined as angioplasty with or without stents). In this study, MACEs included cardiac death, nonfatal myocardial infarction, and subsequent PCI and coronary artery bypass graft. Diabetes was the greatest risk factor for MACEs. Objective: This study estimated the cost-effectiveness of fluvastatin when used for secondary prevention of MACEs after PCI in people with diabetes. Methods: A post hoc subgroup analysis of patients with diabetes from the LIPS was used to estimate the effectiveness of fluvastatin in reducing myocardial infarction, revascularization, and cardiac death. A probabilistic Markov model was developed using United Kingdom resource and cost data to estimate the additional costs and quality-adjusted life-years (QALYs) gained over 10 years from the perspective of the British National Health Service. The model contained 6 health states, and the transition probabilities were derived from the LIPS data. Crossover from fluvastatin to other lipid-lowering drugs, withdrawal from fluvastatin, and the use of lipid-lowering drugs in the control group were included. Results: In the subgroup of 202 patients with diabetes in the LIPS trial, 18 (15.0%) of 120 fluvastatin patients and 21 (25.6%) of 82 control participants were insulin dependent (P = NS). Compared with the control group, patients treated with fluvastatin can expect to gain an additional mean (SD) of 0.196 (0.139) QALY per patient over 10 years (P < 0.001) and will cost the health service an additional mean (SD) of 10 (E448) (P = NS) (mean [SD] US $16 [$689]). The additional cost per QALY gained was;(51 (US $78). The key determinants of cost-effectiveness included the probabilities of repeat interventions, cardiac death, the cost of fluvastatin, and the time horizon used for the evaluation. Conclusion: Fluvastatin was an economically efficient treatment to prevent MACEs in these patients with diabetes undergoing PCI.

Quantification of lean bodyweight

Background: Lean bodyweight (LBW) has been recommended for scaling drug doses. However, the current methods for predicting LBW are inconsistent at extremes of size and could be misleading with respect to interpreting weight-based regimens. Objective: The objective of the present study was to develop a semi-mechanistic model to predict fat-free mass (FFM) from subject characteristics in a population that includes extremes of size. FFM is considered to closely approximate LBW. There are several reference methods for assessing FFM, whereas there are no reference standards for LBW. Patients and methods: A total of 373 patients (168 male, 205 female) were included in the study. These data arose from two populations. Population A (index dataset) contained anthropometric characteristics, FFM estimated by dual-energy x-ray absorptiometry (DXA - a reference method) and bioelectrical impedance analysis (BIA) data. Population B (test dataset) contained the same anthropometric measures and FFM data as population A, but excluded BIA data. The patients in population A had a wide range of age (18-82 years), bodyweight (40.7-216.5kg) and BMI values (17.1-69.9 kg/m(2)). Patients in population B had BMI values of 18.7-38.4 kg/m(2). A two-stage semi-mechanistic model to predict FFM was developed from the demographics from population A. For stage 1 a model was developed to predict impedance and for stage 2 a model that incorporated predicted impedance was used to predict FFM. These two models were combined to provide an overall model to predict FFM from patient characteristics. The developed model for FFM was externally evaluated by predicting into population B. Results: The semi-mechanistic model to predict impedance incorporated sex, height and bodyweight. The developed model provides a good predictor of impedance for both males and females (r(2) = 0.78, mean error [ME] = 2.30 x 10(-3), root mean square error [RMSE] = 51.56 [approximately 10% of mean]). The final model for FFM incorporated sex, height and bodyweight. The developed model for FFM provided good predictive performance for both males and females (r(2) = 0.93, ME = -0.77, RMSE = 3.33 [approximately 6% of mean]). In addition, the model accurately predicted the FFM of subjects in population B (r(2) = 0.85, ME -0.04, RMSE = 4.39 [approximately 7% of mean]). Conclusions: A semi-mechanistic model has been developed to predict FFM (and therefore LBW) from easily accessible patient characteristics. This model has been prospectively evaluated and shown to have good predictive performance.

Optimizing Presence–Absence Surveys For Detecting Population Trends

Presence-absence surveys are a commonly used method for monitoring broad-scale changes in wildlife distributions. However, the lack of power of these surveys for detecting population trends is problematic for their application in wildlife management. Options for improving power include increasing the sampling effort or arbitrarily relaxing the type I error rate. We present an alternative, whereby targeted sampling of particular habitats in the landscape using information from a habitat model increases power. The advantage of this approach is that it does not require a trade-off with either cost or the Pr(type I error) to achieve greater power. We use a demographic model of koala (Phascolarctos cinereus) population dynamics and simulations of the monitoring process to estimate the power to detect a trend in occupancy for a range of strategies, thereby demonstrating that targeting particular habitat qualities can improve power substantially. If the objective is to detect a decline in occupancy, the optimal strategy is to sample high-quality habitats. Alternatively, if the objective is to detect an increase in occupancy, the optimal strategy is to sample intermediate-quality habitats. The strategies with the highest power remained the same under a range of parameter assumptions, although observation error had a strong influence on the optimal strategy. Our approach specifically applies to monitoring for detecting long-term trends in occupancy or abundance. This is a common and important monitoring objective for wildlife managers, and we provide guidelines for more effectively achieving it.

Speeding up the EM algorithm for mixture model-based segmentation of magnetic resonance images

Mixture models implemented via the expectation-maximization (EM) algorithm are being increasingly used in a wide range of problems in pattern recognition such as image segmentation. However, the EM algorithm requires considerable computational time in its application to huge data sets such as a three-dimensional magnetic resonance (MR) image of over 10 million voxels. Recently, it was shown that a sparse, incremental version of the EM algorithm could improve its rate of convergence. In this paper, we show how this modified EM algorithm can be speeded up further by adopting a multiresolution kd-tree structure in performing the E-step. The proposed algorithm outperforms some other variants of the EM algorithm for segmenting MR images of the human brain. (C) 2004 Pattern Recognition Society. Published by Elsevier Ltd. All rights reserved.

A simulation model of cereal-legume intercropping systems for semi-arid regions I. Model development

Cereal-legume intercropping plays an important role in subsistence food production in developing countries, especially in situations of limited water resources. Crop simulation can be used to assess risk for intercrop productivity over time and space. In this study, a simple model for intercropping was developed for cereal and legume growth and yield, under semi-arid conditions. The model is based on radiation interception and use, and incorporates a water stress factor. Total dry matter and yield are functions of photosynthetically active radiation (PAR), the fraction of radiation intercepted and radiation use efficiency (RUE). One of two PAR sub-models was used to estimate PAR from solar radiation; either PAR is 50% of solar radiation or the ratio of PAR to solar radiation (PAR/SR) is a function of the clearness index (K-T). The fraction of radiation intercepted was calculated either based on Beer's Law with crop extinction coefficients (K) from field experiments or from previous reports. RUE was calculated as a function of available soil water to a depth of 900 mm (ASW). Either the soil water balance method or the decay curve approach was used to determine ASW. Thus, two alternatives for each of three factors, i.e., PAR/SR, K and ASW, were considered, giving eight possible models (2 methods x 3 factors). The model calibration and validation were carried out with maize-bean intercropping systems using data collected in a semi-arid region (Bloemfontein, Free State, South Africa) during seven growing seasons (1996/1997-2002/2003). The combination of PAR estimated from the clearness index, a crop extinction coefficient from the field experiment and the decay curve model gave the most reasonable and acceptable result. The intercrop model developed in this study is simple, so this modelling approach can be employed to develop other cereal-legume intercrop models for semi-arid regions. (c) 2004 Elsevier B.V. All rights reserved.

A simulation model of cereal-legume intercropping systems for semi-arid regions II. Model application

Smallholder farmers in Africa practice traditional cropping techniques such as intercropping. Intercropping is thought to offer higher productivity and resource milisation than sole cropping. In this study, risk associated with maize-bean intercropping was evaluated by quantifying long-term yield in both intercropping and sole cropping in a semi-arid region of South Africa (Bloemfontein, Free State) with reference to rainfall variability. The crop simulation model was run with different cultural practices (planting date and plant density) for 52 summer crop growing seasons (1950/1951-2001/2002). Eighty-one scenarios, consisted of three levels of initial soil water, planting date, maize population, and bean population, were simulated. From the simulation outputs, the total land equivalent ratio (LER) was greater than one. The intercrop (equivalent to sole maize) had greater energy value (EV) than sole beans, and the intercrop (equivalent to sole beans) had greater monetary value (MV) than sole maize. From these results, it can be concluded that maize-bean intercropping is advantageous for this semi-arid region. Soil water at planting was the most important factor of all scenario factors, followed by planting date. Irrigation application at planting, November/December planting and high plant density of maize for EV and beans for MV can be one of the most effective cultural practices in the study region. With regard to rainfall variability, seasonal (October-April) rainfall positively affected EV and MV, but not LER. There was more intercrop production in La Nina years than in El Nino years. Thus, better cultural practices may be selected to maximize maize-bean intercrop yields for specific seasons in the semi-arid region based on the global seasonal outlook. (c) 2004 Elsevier B.V. All rights reserved.