Protein and lipid binding parameters in rainbow trout (Oncorhynchus mykiss) blood and liver fractions to extrapolate from an in vitro metabolic degradation assay to in vivo bioaccumulation potential of hydrophobic organic chemicals

Binding of hydrophobic chemicals to colloids such as proteins or lipids is difficult to measure using classical microdialysis methods due to low aqueous concentrations, adsorption to dialysis membranes and test vessels, and slow kinetics of equilibration. Here, we employed a three-phase partitioning system where silicone (polydimethylsiloxane, PDMS) serves as a third phase to determine partitioning between water and colloids and acts at the same time as a dosing device for hydrophobic chemicals. The applicability of this method was demonstrated with bovine serum albumin (BSA). Measured binding constants (K(BSAw)) for chlorpyrifos, methoxychlor, nonylphenol, and pyrene were in good agreement with an established quantitative structure-activity relationship (QSAR). A fifth compound, fluoxypyr-methyl-heptyl ester, was excluded from the analysis because of apparent abiotic degradation. The PDMS depletion method was then used to determine partition coefficients for test chemicals in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (K(S9w)) and blood plasma (K(bloodw)). Measured K(S9w) and K(bloodw) values were consistent with predictions obtained using a mass-balance model that employs the octanol-water partition coefficient (K(ow)) as a surrogate for lipid partitioning and K(BSAw) to represent protein binding. For each compound, K(bloodw) was substantially greater than K(S9w), primarily because blood contains more lipid than liver S9 fractions (1.84% of wet weight vs 0.051%). Measured liver S9 and blood plasma binding parameters were subsequently implemented in an in vitro to in vivo extrapolation model to link the in vitro liver S9 metabolic degradation assay to in vivo metabolism in fish. Apparent volumes of distribution (V(d)) calculated from the experimental data were similar to literature estimates. However, the calculated binding ratios (f(u)) used to relate in vitro metabolic clearance to clearance by the intact liver were 10 to 100 times lower than values used in previous modeling efforts. Bioconcentration factors (BCF) predicted using the experimental binding data were substantially higher than the predicted values obtained in earlier studies and correlated poorly with measured BCF values in fish. One possible explanation for this finding is that chemicals bound to proteins can desorb rapidly and thus contribute to metabolic turnover of the chemicals. This hypothesis remains to be investigated in future studies, ideally with chemicals of higher hydrophobicity.

Metabolomics reveals the metabolic map of procainamide in humans and mice

Procainamide, a type I antiarrhythmic agent, is used to treat a variety of atrial and ventricular dysrhythmias. It was reported that long-term therapy with procainamide may cause lupus erythematosus in 25-30% of patients. Interestingly, procainamide does not induce lupus erythematosus in mouse models. To explore the differences in this side-effect of procainamide between humans and mouse models, metabolomic analysis using ultra-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC-ESI-QTOFMS) was conducted on urine samples from procainamide-treated humans, CYP2D6-humanized mice, and wild-type mice. Thirteen urinary procainamide metabolites, including nine novel metabolites, derived from P450-dependent, FMO-dependent oxidations and acylation reactions, were identified and structurally elucidated. In vivo metabolism of procainamide in CYP2D6-humanized mice as well as in vitro incubations with microsomes and recombinant P450s suggested that human CYP2D6 plays a major role in procainamide metabolism. Significant differences in N-acylation and N-oxidation of the drug between humans and mice largely account for the interspecies differences in procainamide metabolism. Significant levels of the novel N-oxide metabolites produced by FMO1 and FMO3 in humans might be associated with the development of procainamide-induced systemic lupus erythematosus. Observations based on this metabolomic study offer clues to understanding procainamide-induced lupus in humans and the effect of P450s and FMOs on procainamide N-oxidation.

Metabolic effects of dialyzate glucose in chronic hemodialysis: results from a prospective, randomized crossover trial

There is no agreement concerning dialyzate glucose concentration in hemodialysis (HD) and 100 and 200 mg/dL (G100 and G200) are frequently used. G200 may result in diffusive glucose flux into the patient, with consequent hyperglycemia and hyperinsulinism, and electrolyte alterations, in particular potassium (K) and phosphorus (P). This trial compared metabolic effects of G100 versus G200.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Establishment of a novel intervertebral disc/endplate culture model: analysis of an ex vivo in vitro whole-organ rabbit culture system

STUDY DESIGN: Ex vivo in vitro study evaluating a novel intervertebral disc/endplate culture system. OBJECTIVES: To establish a whole-organ intervertebral disc culture model for the study of disc degeneration in vitro, including the characterization of basic cell and organ function. SUMMARY OF BACKGROUND DATA: With current in vivo models for the study of disc and endplate degeneration, it remains difficult to investigate the complex disc metabolism and signaling cascades. In contrast, more controlled but simplified in vitro systems using isolated cells or disc fragments are difficult to culture due to the unconstrained conditions, with often-observed cell death or cell dedifferentiation. Therefore, there is a demand for a controlled culture model with preserved cell function that offers the possibility to investigate disc and endplate pathologies in a structurally intact organ. METHODS: Naturally constrained intervertebral disc/endplate units from rabbits were cultured in multi-well plates. Cell viability, metabolic activity, matrix composition, and matrix gene expression profile were monitored using the Live/Dead cell viability test (Invitrogen, Basel, Switzerland), tetrazolium salt reduction (WST-8), proteoglycan and deoxyribonucleic acid quantification assays, and quantitative polymerase chain reaction. RESULTS: Viability and organ integrity were preserved for at least 4 weeks, while proteoglycan and deoxyribonucleic acid content decreased slightly, and matrix genes exhibited a degenerative profile with up-regulation of type I collagen and suppression of collagen type II and aggrecan genes. Additionally, cell metabolic activity was reduced to one third of the initial value. CONCLUSIONS: Naturally constrained intervertebral rabbit discs could be cultured for several weeks without losing cell viability. Structural integrity and matrix composition were retained. However, the organ responded to the artificial environment with a degenerative gene expression pattern and decreased metabolic rate. Therefore, the described system serves as a promising in vitro model to study disc degeneration in a whole organ.

Late-onset manifestation of antenatal Bartter syndrome as a result of residual function of the mutated renal Na+-K+-2Cl- co-transporter

Genetic defects of the Na+-K+-2Cl- (NKCC2) sodium potassium chloride co-transporter result in severe, prenatal-onset renal salt wasting accompanied by polyhydramnios, prematurity, and life-threatening hypovolemia of the neonate (antenatal Bartter syndrome or hyperprostaglandin E syndrome). Herein are described two brothers who presented with hyperuricemia, mild metabolic alkalosis, low serum potassium levels, and bilateral medullary nephrocalcinosis at the ages of 13 and 15 yr. Impaired function of sodium chloride reabsorption along the thick ascending limb of Henle's loop was deduced from a reduced increase in diuresis and urinary chloride excretion upon application of furosemide. Molecular genetic analysis revealed that the brothers were compound heterozygotes for mutations in the SLC12A1 gene coding for the NKCC2 co-transporter. Functional analysis of the mutated rat NKCC2 protein by tracer-flux assays after heterologous expression in Xenopus oocytes revealed significant residual transport activity of the NKCC2 p.F177Y mutant construct in contrast to no activity of the NKCC2-D918fs frameshift mutant construct. However, coexpression of the two mutants was not significantly different from that of NKCC2-F177Y alone or wild type. Membrane expression of NKCC2-F177Y as determined by luminometric surface quantification was not significantly different from wild-type protein, pointing to an intrinsic partial transport defect caused by the p.F177Y mutation. The partial function of NKCC2-F177Y, which is not negatively affected by NKCC2-D918fs, therefore explains a mild and late-onset phenotype and for the first time establishes a mild phenotype-associated SLC12A1 gene mutation.

Risk of adverse effects of intensified treatment in insulin-dependent diabetes mellitus: a meta-analysis

While the benefits of intensified insulin treatment in insulin-dependent (Type 1) diabetes mellitus (IDDM) are well recognized, the risks have not been comprehensively characterized. We examined the risk of severe hypoglycaemia, ketoacidosis, and death in a meta-analysis of randomized controlled trials. The MEDLINE database, reference lists, and specialist journals were searched electronically or by hand to identify relevant studies with at least 6 months of follow-up and the monitoring of glycaemia by glycosylated haemoglobin measurements. Logistic regression was used for calculation of combined odds ratios and 95% confidence intervals (95% CI). The influence of covariates was examined by including covariate-by-treatment interaction terms. Methodological study quality was assessed and sensitivity analyses were performed. Fourteen trials were identified. These contributed 16 comparisons with 1028 patients allocated to intensified and 1039 allocated to conventional treatment. A total of 846 patients suffered at least one episode of severe hypoglycaemia, 175 patients experienced ketoacidosis and 26 patients died. The combined odds ratio (95% CI) for hypoglycaemia was 2.99 (2.45-3.64), for ketoacidosis 1.74 (1.27-2.38) and for death from all causes 1.40 (0.65-3.01). The risk of severe hypoglycaemia was determined by the degree of normalization of glycaemia achieved (p=0.005 for interaction term), with the results from the Diabetes Control and Complications Trial (DCCT) in line with the other trials. Ketoacidosis risk depended on the type of intensified treatment used. Odds ratios (95% CI) were 7.20 (2.95-17.58) for exclusive use of pumps, 1.13 (0.15-8.35) for multiple daily injections and 1.28 (0.90-1.83) for trials offering a choice between the two (p = 0.004 for interaction). Mortality was significantly (p = 0.007) increased for causes potentially associated with acute complications (7 vs 0 deaths, 5 deaths attributed to ketoacidosis, and 2 sudden deaths), and non-significantly (p = 0.16) decreased for macrovascular causes (3 vs 8 deaths). We conclude that there is a substantial risk of severe adverse effects associated with intensified insulin treatment. Mortality from acute metabolic causes is increased; however, this is largely counterbalanced by a reduction in cardiovascular mortality. The excess of severe hypoglycemia in the DCCT is not exceptional. Multiple daily injection schemes may be safer than treatment with insulin pumps.

Liver repopulation with bone marrow derived cells improves the metabolic disorder in the Gunn rat

BACKGROUND: Reversible ischaemia/reperfusion (I/R) liver injury has been used to induce engraftment and hepatic parenchymal differentiation of exogenous beta2-microglubulin(-)/Thy1(+) bone marrow derived cells. AIM: To test the ability of this method of hepatic parenchymal repopulation, theoretically applicable to clinical practice, to correct the metabolic disorder in a rat model of congenital hyperbilirubinaemia. METHODS AND RESULTS: Analysis by confocal laser microscopy of fluorescence labelled cells and by immunohistochemistry for beta2-microglubulin, 72 hours after intraportal delivery, showed engraftment of infused cells in liver parenchyma of rats with I/R, but not in control animals with non-injured liver. Transplantation of bone marrow derived cells obtained from GFP-transgenic rats into Lewis rats resulted in the presence of up to 20% of GFP positive hepatocytes in I/R liver lobes after one month. The repopulation rate was proportional to the number of transplanted cells. Infusion of GFP negative bone marrow derived cells into GFP positive transgenic rats resulted in the appearance of GFP negative hepatocytes, suggesting that the main mechanism underlying parenchymal repopulation was differentiation rather than cell fusion. Transplantation of wild type bone marrow derived cells into hyperbilirubinaemic Gunn rats with deficient bilirubin conjugation after I/R damage resulted in 30% decrease in serum bilirubin, the appearance of bilirubin conjugates in bile, and the expression of normal UDP-glucuronyltransferase enzyme evaluated by polymerase chain reaction. CONCLUSIONS: I/R injury induced hepatic parenchymal engraftment and differentiation into hepatocyte-like cells of bone marrow derived cells. Transplantation of bone marrow derived cells from non-affected animals resulted in the partial correction of hyperbilirubinaemia in the Gunn rat.

Depression, the metabolic syndrome and cardiovascular risk

BACKGROUND: The relationship between depression and the metabolic syndrome is unclear, and whether metabolic syndrome explains the association between depression and cardiovascular disease (CVD) risk is unknown. METHODS: We studied 652 women who received coronary angiography as part of the Women's Ischemia Syndrome Evaluation (WISE) study and completed the Beck Depression Inventory (BDI). Women who had both elevated depressive symptoms (BDI > or =10) and a previous diagnosis of depression were considered at highest risk, whereas those with one of the two conditions represented an intermediate group. The metabolic syndrome was defined according to the ATP-III criteria. The main outcome was incidence of adverse CVD events (hospitalizations for myocardial infarction, stroke, congestive heart failure, and CVD-related mortality) over a median follow-up of 5.9 years. RESULTS: After adjusting for demographic factors, lifestyle and functional status, both depression categories were associated with about 60% increased odds for metabolic syndrome compared with no depression (p = .03). The number of metabolic syndrome risk factors increased gradually across the three depression categories (p = .003). During follow-up, 104 women (15.9%) experienced CVD events. In multivariable analysis, women with both elevated symptoms and a previous diagnosis of depression had 2.6 times higher risk of CVD. When metabolic syndrome was added to the model, the risk associated with depression only decreased by 7%, and both depression and metabolic syndrome remained significant predictors of CVD. CONCLUSIONS: In women with suspected coronary artery disease, the metabolic syndrome is independently associated with depression but explains only a small portion of the association between depression and incident CVD.

Human and embodied energy analysis applied to water source protection and household water treatment interventions used in Mali, West Africa

As water quality interventions are scaled up to meet the Millennium Development Goal of halving the proportion of the population without access to safe drinking water by 2015 there has been much discussion on the merits of household- and source-level interventions. This study furthers the discussion by examining specific interventions through the use of embodied human and material energy. Embodied energy quantifies the total energy required to produce and use an intervention, including all upstream energy transactions. This model uses material quantities and prices to calculate embodied energy using national economic input/output-based models from China, the United States and Mali. Embodied energy is a measure of aggregate environmental impacts of the interventions. Human energy quantifies the caloric expenditure associated with the installation and operation of an intervention is calculated using the physical activity ratios (PARs) and basal metabolic rates (BMRs). Human energy is a measure of aggregate social impacts of an intervention. A total of four household treatment interventions – biosand filtration, chlorination, ceramic filtration and boiling – and four water source-level interventions – an improved well, a rope pump, a hand pump and a solar pump – are evaluated in the context of Mali, West Africa. Source-level interventions slightly out-perform household-level interventions in terms of having less total embodied energy. Human energy, typically assumed to be a negligible portion of total embodied energy, is shown to be significant to all eight interventions, and contributing over half of total embodied energy in four of the interventions. Traditional gender roles in Mali dictate the types of work performed by men and women. When the human energy is disaggregated by gender, it is seen that women perform over 99% of the work associated with seven of the eight interventions. This has profound implications for gender equality in the context of water quality interventions, and may justify investment in interventions that reduce human energy burdens.

Analysis of in-cylinder pressure transducer data quality utilizing a SIDI turbocharged engine

In-cylinder pressure transducers have been used for decades to record combustion pressure inside a running engine. However, due to the extreme operating environment, transducer design and installation must be considered in order to minimize measurement error. One such error is caused by thermal shock, where the pressure transducer experiences a high heat flux that can distort the pressure transducer diaphragm and also change the crystal sensitivity. This research focused on investigating the effects of thermal shock on in-cylinder pressure transducer data quality using a 2.0L, four-cylinder, spark-ignited, direct-injected, turbo-charged GM engine. Cylinder four was modified with five ports to accommodate pressure transducers of different manufacturers. They included an AVL GH14D, an AVL GH15D, a Kistler 6125C, and a Kistler 6054AR. The GH14D, GH15D, and 6054AR were M5 size transducers. The 6125C was a larger, 6.2mm transducer. Note that both of the AVL pressure transducers utilized a PH03 flame arrestor. Sweeps of ignition timing (spark sweep), engine speed, and engine load were performed to study the effects of thermal shock on each pressure transducer. The project consisted of two distinct phases which included experimental engine testing as well as simulation using a commercially available software package. A comparison was performed to characterize the quality of the data between the actual cylinder pressure and the simulated results. This comparison was valuable because the simulation results did not include thermal shock effects. All three sets of tests showed the peak cylinder pressure was basically unaffected by thermal shock. Comparison of the experimental data with the simulated results showed very good correlation. The spark sweep was performed at 1300 RPM and 3.3 bar NMEP and showed that the differences between the simulated results (no thermal shock) and the experimental data for the indicated mean effective pressure (IMEP) and the pumping mean effective pressure (PMEP) were significantly less than the published accuracies. All transducers had an IMEP percent difference less than 0.038% and less than 0.32% for PMEP. Kistler and AVL publish that the accuracy of their pressure transducers are within plus or minus 1% for the IMEP (AVL 2011; Kistler 2011). In addition, the difference in average exhaust absolute pressure between the simulated results and experimental data was the greatest for the two Kistler pressure transducers. The location and lack of flame arrestor are believed to be the cause of the increased error. For the engine speed sweep, the torque output was held constant at 203 Nm (150 ft-lbf) from 1500 to 4000 RPM. The difference in IMEP was less than 0.01% and the PMEP was less than 1%, except for the AVL GH14D which was 5% and the AVL GH15DK which was 2.25%. A noticeable error in PMEP appeared as the load increased during the engine speed sweeps, as expected. The load sweep was conducted at 2000 RPM over a range of NMEP from 1.1 to 14 bar. The difference in IMEP values were less 0.08% while the PMEP values were below 1% except for the AVL GH14D which was 1.8% and the AVL GH15DK which was at 1.25%. In-cylinder pressure transducer data quality was effectively analyzed using a combination of experimental data and simulation results. Several criteria can be used to investigate the impact of thermal shock on data quality as well as determine the best location and thermal protection for various transducers.

DENITRIFICATION IN SOILS: FROM GENES TO ENVIRONMENTAL OUTCOMES

Denitrification is an important process of global nitrogen cycle as it removes reactive nitrogen from the biosphere, and acts as the primary source of nitrous oxide (N2O). This thesis seeks to gain better understanding of the biogeochemistry of denitrification by investigating the process from four different aspects: genetic basis, enzymatic kinetics, environmental interactions, and environmental consequences. Laboratory and field experiments were combined with modeling efforts to unravel the complexity of denitrification process under microbiological and environmental controls. Dynamics of denitrification products observed in laboratory experiments revealed an important role of constitutive denitrification enzymes, whose presence were further confirmed with quantitative analysis of functional genes encoding nitrite reductase and nitrous oxide reductase. A metabolic model of denitrification developed with explicit denitrification enzyme kinetics and representation of constitutive enzymes successfully reproduced the dynamics of N2O and N2 accumulation observed in the incubation experiments, revealing important regulatory effect of denitrification enzyme kinetics on the accumulation of denitrification products. Field studies demonstrated complex interaction of belowground N2O production, consumption and transport, resulting in two pulse pattern in the surface flux. Coupled soil gas diffusion/denitrification model showed great potential in simulating the dynamics of N2O below ground, with explicit representation of the activity of constitutive denitrification enzymes. A complete survey of environmental variables showed distinct regulation regimes on the denitrification activity from constitutive enzymes and new synthesized enzymes. Uncertainties in N2O estimation with current biogeochemical models may be reduced as accurate simulation of the dynamics of N2O in soil and surface fluxes is possible with a coupled diffusion/denitrification model that includes explicit representation of denitrification enzyme kinetics. In conclusion, denitrification is a complex ecological function regulated at cellular level. To assess the environmental consequences of denitrification and develop useful tools to mitigate N2O emissions require a comprehensive understanding of the regulatory network of denitrification with respect to microbial physiology and environmental interactions.

Association of arterial blood pressure and vasopressor load with septic shock mortality: a post hoc analysis of a multicenter trial

INTRODUCTION: It is unclear to which level mean arterial blood pressure (MAP) should be increased during septic shock in order to improve outcome. In this study we investigated the association between MAP values of 70 mmHg or higher, vasopressor load, 28-day mortality and disease-related events in septic shock. METHODS: This is a post hoc analysis of data of the control group of a multicenter trial and includes 290 septic shock patients in whom a mean MAP > or = 70 mmHg could be maintained during shock. Demographic and clinical data, MAP, vasopressor requirements during the shock period, disease-related events and 28-day mortality were documented. Logistic regression models adjusted for the geographic region of the study center, age, presence of chronic arterial hypertension, simplified acute physiology score (SAPS) II and the mean vasopressor load during the shock period was calculated to investigate the association between MAP or MAP quartiles > or = 70 mmHg and mortality or the frequency and occurrence of disease-related events. RESULTS: There was no association between MAP or MAP quartiles and mortality or the occurrence of disease-related events. These associations were not influenced by age or pre-existent arterial hypertension (all P > 0.05). The mean vasopressor load was associated with mortality (relative risk (RR), 1.83; confidence interval (CI) 95%, 1.4-2.38; P < 0.001), the number of disease-related events (P < 0.001) and the occurrence of acute circulatory failure (RR, 1.64; CI 95%, 1.28-2.11; P < 0.001), metabolic acidosis (RR, 1.79; CI 95%, 1.38-2.32; P < 0.001), renal failure (RR, 1.49; CI 95%, 1.17-1.89; P = 0.001) and thrombocytopenia (RR, 1.33; CI 95%, 1.06-1.68; P = 0.01). CONCLUSIONS: MAP levels of 70 mmHg or higher do not appear to be associated with improved survival in septic shock. Elevating MAP >70 mmHg by augmenting vasopressor dosages may increase mortality. Future trials are needed to identify the lowest acceptable MAP level to ensure tissue perfusion and avoid unnecessary high catecholamine infusions.

In vivo analysis of efavirenz metabolism in individuals with impaired CYP2A6 function

INTRODUCTION: The antiretroviral drug efavirenz (EFV) is extensively metabolized into three primary metabolites: 8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV. There is a wide interindividual variability in EFV plasma exposure, explained to a great extent by cytochrome P450 2B6 (CYP2B6), the main isoenzyme responsible for EFV metabolism and involved in the major metabolic pathway (8-hydroxylation) and to a lesser extent in 7-hydroxylation. When CYP2B6 function is impaired, the relevance of CYP2A6, the main isoenzyme responsible for 7-hydroxylation may increase. We hypothesize that genetic variability in this gene may contribute to the particularly high, unexplained variability in EFV exposure in individuals with limited CYP2B6 function. METHODS: This study characterized CYP2A6 variation (14 alleles) in individuals (N=169) previously characterized for functional variants in CYP2B6 (18 alleles). Plasma concentrations of EFV and its primary metabolites (8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV) were measured in different genetic backgrounds in vivo. RESULTS: The accessory metabolic pathway CYP2A6 has a critical role in limiting drug accumulation in individuals characterized as CYP2B6 slow metabolizers. CONCLUSION: Dual CYP2B6 and CYP2A6 slow metabolism occurs at significant frequency in various human populations, leading to extremely high EFV exposure.

Hierarchical modeling gave plausible estimates of associations between metabolic syndrome and components of antiretroviral therapy

OBJECTIVE: Hierarchical modeling has been proposed as a solution to the multiple exposure problem. We estimate associations between metabolic syndrome and different components of antiretroviral therapy using both conventional and hierarchical models. STUDY DESIGN AND SETTING: We use discrete time survival analysis to estimate the association between metabolic syndrome and cumulative exposure to 16 antiretrovirals from four drug classes. We fit a hierarchical model where the drug class provides a prior model of the association between metabolic syndrome and exposure to each antiretroviral. RESULTS: One thousand two hundred and eighteen patients were followed for a median of 27 months, with 242 cases of metabolic syndrome (20%) at a rate of 7.5 cases per 100 patient years. Metabolic syndrome was more likely to develop in patients exposed to stavudine, but was less likely to develop in those exposed to atazanavir. The estimate for exposure to atazanavir increased from hazard ratio of 0.06 per 6 months' use in the conventional model to 0.37 in the hierarchical model (or from 0.57 to 0.81 when using spline-based covariate adjustment). CONCLUSION: These results are consistent with trials that show the disadvantage of stavudine and advantage of atazanavir relative to other drugs in their respective classes. The hierarchical model gave more plausible results than the equivalent conventional model.