A wearable system for multi-segment foot kinetics measurement.

This study aims to design a wearable system for kinetics measurement of multi-segment foot joints in long-distance walking and to investigate its suitability for clinical evaluations. The wearable system consisted of inertial sensors (3D gyroscopes and 3D accelerometers) on toes, forefoot, hindfoot, and shank, and a plantar pressure insole. After calibration in a laboratory, 10 healthy elderly subjects and 12 patients with ankle osteoarthritis walked 50m twice wearing this system. Using inverse dynamics, 3D forces, moments, and power were calculated in the joint sections among toes, forefoot, hindfoot, and shank. Compared to those we previously estimated for a one-segment foot model, the sagittal and transverse moments and power in the ankle joint, as measured via multi-segment foot model, showed a normalized RMS difference of less than 11%, 14%, and 13%, respectively, for healthy subjects, and 13%, 15%, and 14%, for patients. Similar to our previous study, the coronal moments were not analyzed. Maxima-minima values of anterior-posterior and vertical force, sagittal moment, and power in shank-hindfoot and hindfoot-forefoot joints were significantly different between patients and healthy subjects. Except for power, the inter-subject repeatability of these parameters was CMC>0.90 for healthy subjects and CMC>0.70 for patients. Repeatability of these parameters was lower for the forefoot-toes joint. The proposed measurement system estimated multi-segment foot joints kinetics with acceptable repeatability but showed difference, compared to those previously estimated for the one-segment foot model. These parameters also could distinguish patients from healthy subjects. Thus, this system is suggested for outcome evaluations of foot treatments.

Effect of age toxicokinetics among human volunteers exposed to propylene glycol monomethyl ether (PGME)

Rationale: Aging adults represent the fastest growing population segment in many countries. Physiological and metabolic changes in the aging process may alter how aging adults respond to exposures compared to younger workers. Current preventive workplace exposure measures may therefore not be sufficiently protective for the aging workforce. In a controlled human toxicokinetic study (exposure chamber; 12m3), the volunteers (n=11) were men and women over the age of 58 years and exposed to a commonly used, low neurotoxic glycol ether; PGME (CAS no. 107-98- 2) (50 ppm, 6 hours). Oxidative metabolism (Michaelis-Menten) is the major pathway and conjugation the minor in humans. Metabolites, conjugated and free PGME are eliminated through the kidneys, and the elimination kinetics is dose-dependent (0 order). Scope: (1) compare the toxicokinetic profile of PGME obtained in the aging volunteers (58- 62 years) to young volunteers (20-25 years) from a previous study; (2) Test the predictive power of an existing PGME toxicokinetic compartment model for aging persons against urinary PGME concentrations found in volunteers from our experimental study. Experimental procedure: Urine samples were collected before, every 2-hour during exposures for six hours, and ad-lib for additional 20 hours. Urinary analysis of free and total PGME was performed using capillary GC/FID. The toxicokinetic model (Berkley Madonna software) was ageadjusted. Results. Urinary free and total PGME concentration rose rapidly, and did not reach an apparent plateau level during exposure. Less conjugation was observed in the older group. The predictive model developed for the young group predicted well total PGME in the aging group but not free PGME. The age adjusted toxicokinetic model's Vmax1 had to be changed for the aging group, implying slower enzymatic pathway. Conclusion: The toxicokinetic model did not predict well if only the physiological parameters were adjusted for aging adults (existing model); a substance specific metabolic rate parameter was also needed.

Interplay between T cell receptor binding kinetics and the level of cognate peptide presented by major histocompatibility complexes governs CD8+ T cell responsiveness.

Through a rational design approach, we generated a panel of HLA-A*0201/NY-ESO-1(157-165)-specific T cell receptors (TCR) with increasing affinities of up to 150-fold from the wild-type TCR. Using these TCR variants which extend just beyond the natural affinity range, along with an extreme supraphysiologic one having 1400-fold enhanced affinity, and a low-binding one, we sought to determine the effect of TCR binding properties along with cognate peptide concentration on CD8(+) T cell responsiveness. Major histocompatibility complexes (MHC) expressed on the surface of various antigen presenting cells were peptide-pulsed and used to stimulate human CD8(+) T cells expressing the different TCR via lentiviral transduction. At intermediate peptide concentration we measured maximum cytokine/chemokine secretion, cytotoxicity, and Ca(2+) flux for CD8(+) T cells expressing TCR within a dissociation constant (K(D)) range of ∼1-5 μM. Under these same conditions there was a gradual attenuation in activity for supraphysiologic affinity TCR with K(D) < ∼1 μM, irrespective of CD8 co-engagement and of half-life (t(1/2) = ln 2/k(off)) values. With increased peptide concentration, however, the activity levels of CD8(+) T cells expressing supraphysiologic affinity TCR were gradually restored. Together our data support the productive hit rate model of T cell activation arguing that it is not the absolute number of TCR/pMHC complexes formed at equilibrium, but rather their productive turnover, that controls levels of biological activity. Our findings have important implications for various immunotherapies under development such as adoptive cell transfer of TCR-engineered CD8(+) T cells, as well as for peptide vaccination strategies.

Impact of biological and environmental variabilities on biological monitoring: an approach using toxicokinetic models

Biological monitoring of occupational exposure is characterized by important variability, due both to variability in the environment and to biological differences between workers. A quantitative description and understanding of this variability is important for a dependable application of biological monitoring. This work describes this variability,using a toxicokinetic model, for a large range of chemicals for which reference biological reference values exist. A toxicokinetic compartmental model describing both the parent compound and its metabolites was used. For each chemical, compartments were given physiological meaning. Models were elaborated based on physiological, physicochemical, and biochemical data when available, and on half-lives and central compartment concentrations when not available. Fourteen chemicals were studied (arsenic, cadmium, carbon monoxide, chromium, cobalt, ethylbenzene, ethyleneglycol monomethylether, fluorides, lead, mercury, methyl isobutyl ketone, penthachlorophenol, phenol, and toluene), representing 20 biological indicators. Occupational exposures were simulated using Monte Carlo techniques with realistic distributions of both individual physiological parameters and exposure conditions. Resulting biological indicator levels were then analyzed to identify the contribution of environmental and biological variability to total variability. Comparison of predicted biological indicator levels with biological exposure limits showed a high correlation with the model for 19 out of 20 indicators. Variability associated with changes in exposure levels (GSD of 1.5 and 2.0) is shown to be mainly influenced by the kinetics of the biological indicator. Thus, with regard to variability, we can conclude that, for the 14 chemicals modeled, biological monitoring would be preferable to air monitoring. For short half-lives (less than 7 hr), this is very similar to the environmental variability. However, for longer half-lives, estimated variability decreased. [Supplementary materials are available for this article. Go to the publisher's online edition of Journal of Occupational and Environmental Hygiene for the following free supplemental resource: tables detailing the CBTK models for all 14 chemicals and the symbol nomenclature that was used.] [Authors]

High activity of Fosfomycin and Rifampin against methicillin-resistant staphylococcus aureus biofilm in vitro and in an experimental foreign-body infection model.

Increasing antimicrobial resistance reduces treatment options for implant-associated infections caused by methicillin-resistant Staphylococcus aureus (MRSA). We evaluated the activity of fosfomycin alone and in combination with vancomycin, daptomycin, rifampin, and tigecycline against MRSA (ATCC 43300) in a foreign-body (implantable cage) infection model. The MICs of the individual agents were as follows: fosfomycin, 1 μg/ml; daptomycin, 0.125 μg/ml; vancomycin, 1 μg/ml; rifampin, 0.04 μg/ml; and tigecycline, 0.125 μg/ml. Microcalorimetry showed synergistic activity of fosfomycin and rifampin at subinhibitory concentrations against planktonic and biofilm MRSA. In time-kill curves, fosfomycin exhibited time-dependent activity against MRSA with a reduction of 2.5 log10 CFU/ml at 128 × the MIC. In the animal model, planktonic bacteria in cage fluid were reduced by <1 log10 CFU/ml with fosfomycin and tigecycline, 1.7 log10 with daptomycin, 2.2 log10 with fosfomycin-tigecycline and fosfomycin-vancomycin, 3.8 log10 with fosfomycin-daptomycin, and >6.0 log10 with daptomycin-rifampin and fosfomycin-rifampin. Daptomycin-rifampin cured 67% of cage-associated infections and fosfomycin-rifampin cured 83%, whereas all single drugs (fosfomycin, daptomycin, and tigecycline) and rifampin-free fosfomycin combinations showed no cure of MRSA cage-associated infections. No emergence of fosfomycin resistance was observed in animals; however, a 4-fold increase in fosfomycin MIC (from 2 to 16 μg/ml) occurred in the fosfomycin-vancomycin group. In summary, the highest eradication of MRSA cage-associated infections was achieved with fosfomycin in combination with rifampin (83%). Fosfomycin may be used in combination with rifampin against MRSA implant-associated infections, but it cannot replace rifampin as an antibiofilm agent.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

A coherent neurobiological framework for functional neuroimaging provided by a model integrating compartmentalized energy metabolism.

Functional neuroimaging has undergone spectacular developments in recent years. Paradoxically, its neurobiological bases have remained elusive, resulting in an intense debate around the cellular mechanisms taking place upon activation that could contribute to the signals measured. Taking advantage of a modeling approach, we propose here a coherent neurobiological framework that not only explains several in vitro and in vivo observations but also provides a physiological basis to interpret imaging signals. First, based on a model of compartmentalized energy metabolism, we show that complex kinetics of NADH changes observed in vitro can be accounted for by distinct metabolic responses in two cell populations reminiscent of neurons and astrocytes. Second, extended application of the model to an in vivo situation allowed us to reproduce the evolution of intraparenchymal oxygen levels upon activation as measured experimentally without substantially altering the initial parameter values. Finally, applying the same model to functional neuroimaging in humans, we were able to determine that the early negative component of the blood oxygenation level-dependent response recorded with functional MRI, known as the initial dip, critically depends on the oxidative response of neurons, whereas the late aspects of the signal correspond to a combination of responses from cell types with two distinct metabolic profiles that could be neurons and astrocytes. In summary, our results, obtained with such a modeling approach, support the concept that both neuronal and glial metabolic responses form essential components of neuroimaging signals.

Fat oxidation kinetics : effect of exercise

ABSTRACT Fat oxidation kinetics: effect of exercise. During graded exercise, absolute whole body fat oxidation rates increase from low to moderate intensities, and then markedly decline at high intensities, implying an exercise intensity (Fatmax) at which the fat oxidation rate is maximal (MFO). The main aim of the present work was to examine the effect of exercise on whole body fat oxidation kinetics. For this purpose, a sinusoidal mathematical model (SIN) has been developped in the first study to provide an accurate description of the shape of fat oxidation kinetics during graded exercise, represented as a function of exercise intensity, and to determine Fatmax and MFO. The SIN model incorporates three independent variables (i.e., dilatation, symmetry, and translation) that correspond to main expected modulations of the basic fat oxidation curve because of factors such as mode of exercise or training status. The results of study 1 showed that the SIN model was a valuable tool to determine Fatmax and MFO, and to precisely characterize and quantify the different shape of fat oxidation kinetics through its three variables. The effectiveness of the SIN model to detect differences in fat oxidation kinetics induced by a specific factor was then confirmed in the second study, which quantitatively described and compared fat oxidation kinetics in two different popular modes of exercise: running and cycling. It was found that the mean fat oxidation kinetics during running was characterized by a greater dilatation and a rightward asymmetry compared with the symmetric parabolic curve in cycling. In the two subsequent studies, the effect of a prior endurance exercise of different intensities and durations on whole body fat oxidation kinetics was examined. Study 3 determined the impact of a 1-h continuous exercise bout at an exercise intensity corresponding to Fatmax on fat oxidation kinetics during a subsequent graded test, while study 4 investigated the effect of an exercise leading to a more pronounced muscle glycogen depletion. The results of these two latter studies showed that fat oxidation rates, MFO, and Fatmax were enhanced following endurance exercise, but were increased to a greater extent with a more severe mucle glycogen depletion, inducing therefore modifications in the postexercise fat oxidation kinetics (i.e., greater dilatation and rightward asymmetry). In perspective, further studies have been suggested 1) to assess physiological meaning of the three independent variables of the SIN model; and 2) to compare the effect of two different training programs on fat oxidation kinetics in obese subjects.

Drug-eluting beads loaded with antiangiogenic agents for chemoembolization: in vitro sunitinib loading and release and in vivo pharmacokinetics in an animal model.

PURPOSE: The combination of embolic beads with a multitargeted tyrosine kinase inhibitor that inhibits tumor vessel growth is suggested as an alternative and improvement to the current standard doxorubicin-eluting beads for use in transarterial chemoembolization. This study demonstrates the in vitro loading and release kinetics of sunitinib using commercially available embolization microspheres and evaluates the in vitro biologic efficacy on cell cultures and the resulting in vivo pharmacokinetics profiles in an animal model. MATERIALS AND METHODS: DC Bead microspheres, 70-150 µm and 100-300 µm (Biocompatibles Ltd., Farnham, United Kingdom), were loaded by immersion in sunitinib solution. Drug release was measured in saline in a USP-approved flow-through apparatus and quantified by spectrophotometry. Activity after release was confirmed in cell culture. For pharmacokinetics and in vivo toxicity evaluation, New Zealand white rabbits received sunitinib either by intraarterial injection of 100-300 µm sized beads or per os. Plasma and liver tissue drug concentrations were assessed by liquid chromatography-tandem mass spectroscopy. RESULTS: Sunitinib loading on beads was close to complete and homogeneous. A total release of 80% in saline was measured, with similar fast-release profiles for both sphere sizes. After embolization, drug plasma levels remained below the therapeutic threshold (< 50 ng/mL), but high concentrations at 6 hours (14.9 µg/g) and 24 hours (3.4 µg/g) were found in the liver tissue. CONCLUSIONS: DC Bead microspheres of two sizes were efficiently loaded with sunitinib and displayed a fast and almost complete release in saline. High liver drug concentrations and low systemic levels indicated the potential of sunitinib-eluting beads for use in embolization.

hMMP9 as predictive factor for response and progression free survival in breast cancer patients treated with bevacizumab and pegylated liposomal doxorubicin (PLD)

Background: The anti-angiogenic drug, bevacizumab (Bv), is currently used in the treatment of different malignancies including breast cancer. Many angiogenesis-associated molecules are found in the circulation of cancer patients. Until now, there are no prognostic or predictive factors identified in breast cancer patients treated with Bv. We present here the first results of the prospective monitoring of 6 angiogenesis-related molecules in the peripheral blood of breast cancer patients treated with a combination of Bv and PLD in the phase II trial, SAKK 24/06. Methods: Patients were treated with PLD (20 mg/m2) and Bv (10 mg/kg) on days 1 and 15 of each 4-week cycle for a maximum of 6 cycles, followed by Bv monotherapy maintenance (10 mg/m2 q2 weeks) until progression or severe toxicity. Plasma and serum samples were collected at baseline, after 2 months of therapy, then every 3 months and at treatment discontinuation. Enzyme-linked immunosorbent assays (Quantikine, R&D Systems and Reliatech) were used to measure the expression levels of human vascular endothelial growth factor (hVEGF), placental growth factor (hPlGF), matrix metalloproteinase 9 (hMMP9) and soluble VEGF receptors hsVEGFR-1, hsVEGFR-2 and hsVEGFR-3. The log-transformed data (to reduce the skewness) for each marker was analyzed using an analysis of variance (ANOVA) model to determine if there was a difference between the mean of the subgroups of interest (where α = 0.05). The untransformed data was also analyzed in the same manner as a "sensitivity" check. Results: 132 blood samples were collected in 41 out of 43 enrolled patients. Baseline levels of the molecules were compared to disease status according to RECIST. There was a statistically significant difference in the mean of the log-transformed levels of hMMP9 between responders [CR+PR] versus the mean in patients with PD (p-value=0.0004, log fold change=0.7536), and between patients with disease control [CR+PR+SD] and those with PD (p-value=<0.0001, log fold change=0.81559), with the log-transformed level of hMMP9 being higher for the responder group. The mean of the log-transformed levels of hsVEGFR-1 was statistically significantly different between patients with disease control [CR+PR+SD] and those with PD (p-value=0.0068, log fold change=-0.6089), where the log-transformed level of hsVEGFR-1 was lower for the responder group. The log-transformed level of hMMP9 at baseline was identified as a significant prognostic factor in terms of progression free survival (PFS): p-value=0.0417, hazard ratio (HR)=0.574 with a corresponding 95% confidence interval (0.336 - 0.979)). No strong correlation was shown either between the log-transformed levels of hsVEGF, hPlGF, hsVEGFR-2 or hsVEGFR-3 and clinical response or the occurrence of severe toxicity, or between the levels of the different molecules. Conclusions: Our results suggest that baseline plasma level of the matrix metalloproteinase, hMMP9, could predict tumor response and PFS in patients treated with a combination of Bv and PLD. These data justify further investigation in breast cancer patients treated with anti-angiogenic therapy.

Preparation, maintenance, and use of serum-free aggregating brain cell cultures.

Serum-free aggregating brain cell cultures are free-floating three-dimensional primary cell cultures able to reconstitute spontaneously a histotypic brain architecture to reproduce critical steps of brain development and to reach a high level of structural and functional maturity. This culture system offers, therefore, a unique model for neurotoxicity testing both during the development and at advanced cellular differentiation, and the high number of aggregates available combined with the excellent reproducibility of the cultures facilitates routine test procedures. This chapter presents a detailed description of the preparation, maintenance, and use of these cultures for neurotoxicity studies and a comparison of the developmental characteristics between cultures derived from the telencephalon and cultures derived from the whole brain. For culture preparation, mechanically dissociated embryonic brain tissue is used. The initial cell suspension, composed of neural stem cells, neural progenitor cells, immature postmitotic neurons, glioblasts, and microglial cells, is kept in a serum-free, chemically defined medium under continuous gyratory agitation. Spherical aggregates form spontaneously and are maintained in suspension culture for several weeks. Within the aggregates, the cells rearrange and mature, reproducing critical morphogenic events, such as migration, proliferation, differentiation, synaptogenesis, and myelination. For experimentation, replicate cultures are prepared by the randomization of aggregates from several original flasks. The high yield and reproducibility of the cultures enable multiparametric endpoint analyses, including "omics" approaches.

Comparative localization of glioma-reactive monoclonal antibodies in vivo in an athymic mouse human glioma xenograft model.

Radioiodinated murine monoclonal antibodies (Mabs) 81C6, Me 1-14, C12, D12, and E9, made against or reactive with human gliomas but not normal brain, and Mab UJ13A, a pan-neuroectodermal Mab reactive with normal human glial and neural cells, were evaluated in paired label studies in the D-54 MG subcutaneous human glioma xenograft model system in nude mice. Following intravenous injection in the tail vein of mice bearing 200-400 mm3 tumors, specific localization of Mabs to tumor over time (6 h-9 days) was evaluated by tissue counting; each Mab demonstrated a unique localization profile. The comparison of localization indices (LI), determined as a ratio of tissue level of Mab to control immunoglobulin with simultaneous correction for blood levels of each, showed Mabs 81C6 and Me 1-14 to steadily accumulate in glioma xenografts, maintaining LI from 5-20 at 7-9 days after Mab injection. Mab UJ13A peaked at day 1, maintaining this level through day 2, and declining thereafter. Mabs D12 and C12 peaked at days 3 and 4, respectively, and E9 maintained an LI of greater than 3 from days 3-9. Percent injected dose localized/g of tumor varied from a peak high of 16% (81C6) to a low of 5% (Me 1-14 and UJ13A). Immunoperoxidase histochemistry, performed with each Mab on a battery of primary human brain neoplasms, revealed that Mabs 81C6 and E9, which demonstrated the highest levels of percent injected dose localized/g of tumor over time, reacted with antigens expressed in the extracellular matrix. This finding suggests that extracellular matrix localization of antigen represents a biologically significant factor affecting localization and/or binding in the xenograft model used. The demonstration of significant localization, varied kinetics and patterns of localization of this localizing Mab panel warrants their continued investigation as potential imaging and therapeutic agents for human trials.

Effect of a 1-hour single bout of moderate-intensity exercise on fat oxidation kinetics.

The present study aimed to examine the effects of a prior 1-hour continuous exercise bout (CONT) at an intensity (Fat(max)) that elicits the maximal fat oxidation (MFO) on the fat oxidation kinetics during a subsequent submaximal incremental test (IncrC). Twenty moderately trained subjects (9 men and 11 women) performed a graded test on a treadmill (Incr), with 3-minute stages and 1-km.h(-1) increments. Fat oxidation was measured using indirect calorimetry and plotted as a function of exercise intensity. A mathematical model (SIN) including 3 independent variables (dilatation, symmetry, and translation) was used to characterize the shape of fat oxidation kinetics and to determine Fat(max) and MFO. On a second visit, the subjects performed CONT at Fat(max) followed by IncrC. After CONT performed at 57% +/- 3% (means +/- SE) maximal oxygen uptake (Vo(2max)), the respiratory exchange ratio during IncrC was lower at every stage compared with Incr (P < .05). Fat(max) (56.4% +/- 2.3% vs 51.5% +/- 2.4% Vo(2max), P = .013), MFO (0.50 +/- 0.03 vs 0.40 +/- 0.03 g.min(-1), P < .001), and fat oxidation rates from 35% to 70% Vo(2max) (P < .05) were significantly greater during IncrC compared with Incr. However, dilatation and translation were not significantly different (P > .05), whereas symmetry tended to be greater in IncrC (P = .096). This study showed that the prior 1-hour continuous moderate-intensity exercise bout increased Fat(max), MFO, and fat oxidation rates over a wide range of intensities during the postexercise incremental test. Moreover, the shape of the postexercise fat oxidation kinetics tended to have a rightward asymmetry.

Determination of four sequential stages during microautophagy in vitro.

Microautophagy is the transfer of cytosolic components into the lysosome by direct invagination of the lysosomal membrane and subsequent budding of vesicles into the lysosomal lumen. This process is topologically equivalent to membrane invagination during multivesicular body formation and to the budding of enveloped viruses. Vacuoles are lysosomal compartments of yeasts. Vacuolar membrane invagination can be reconstituted in vitro with purified yeast vacuoles, serving as a model system for budding of vesicles into the lumen of an organelle. Using this in vitro system, we defined different reaction states. We identified inhibitors of microautophagy in vitro and used them as tools for kinetic analysis. This allowed us to characterize four biochemically distinguishable steps of the reaction. We propose that these correspond to sequential stages of vacuole invagination and vesicle scission. Formation of vacuolar invaginations was slow and temperature-dependent, whereas the final scission of the vesicle from a preformed invagination was fast and proceeded even on ice. Our observations suggest that the formation of invaginations rather than the scission of vesicles is the rate-limiting step of the overall reaction.

Prediction of resting energy expenditure from fat-free mass and fat mass.

On the basis of literature values, the relationship between fat-free mass (FFM), fat mass (FM), and resting energy expenditure [REE (kJ/24 h)] was determined for 213 adults (86 males, 127 females). The objectives were to develop a mathematical model to predict REE based on body composition and to evaluate the contribution of FFM and FM to REE. The following regression equations were derived: 1) REE = 1265 + (93.3 x FFM) (r2 = 0.727, P < 0.001); 2) REE = 1114 + (90.4 x FFM) + (13.2 x FM) (R2 = 0.743, P < 0.001); and 3) REE = (108 x FFM) + (16.9 x FM) (R2 = 0.986, P < 0.001). FM explained only a small part of the variation remaining after FFM was accounted for. The models that include both FFM and FM are useful in examination of the changes in REE that occur with a change in both the FFM and FM. To account for more of the variability in REE, FFM will have to be divided into organ mass and skeletal muscle mass in future analyses.