Respiratory variability in mechanically ventilated patients. Critical Care 2011

Introduction: Increased respiratory pattern variability is associated with improved oxygenation. Pressure support (PS) is a widely used partial-assist mechanical ventilation (MV) mode, in which each breathing cycle is initiated by flow or pressure variation at the airway due to patient inspiratory effort. Neurally adjusted ventilatory assist (NAVA) is relatively new and uses the electrical activity of the diaphragm (Eadi) to deliver ventilatory support proportional to the patient's inspiratory demand. We hypothesize that respiratory variability should be greater with NAVA compared with PS.Methods: Twenty-two patients underwent 20 minutes of PS followed by 20 minutes of NAVA. Flow and Eadi curves were used to obtain tidal volume (Vt) and ∫Eadi for 300 to 400 breaths in each patient. Patient-specific cumulative distribution functions (CDF) show the percentage Vt and ∫Eadi within a clinically defined (±10%) variability band for each patient. Values are normalized to patient-specific medians for direct comparison. Variability in Vt (outcome) is thus expressed in terms of variability in ∫Eadi (demand) on the same plot.Results: Variability in Vt relative to variability in ∫Eadi is significantly greater for NAVA than PS (P = 0.00012). Hence, greater variability in outcome Vt is obtained for a given demand in ∫Eadi, under NAVA, as illustrated in Figure 1 for a typical patient. A Fisher 2 × 2 contingency analysis showed that 45% of patients under NAVA had a Vt variability in equal proportion to ∫Eadi variability, versus 0% for PS (P < 0.05).Conclusions: NAVA yields greater variability in tidal volume, relative to ∫Eadi demand, and a better match between Vt and ∫Eadi. These results indicate that NAVA could achieve improved oxygenation compared with PS when sufficient underlying variability in ∫Eadi is present, due to its ability to achieve higher tidal volume variability from a given variability in ∫Eadi.

Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Comparison of the glomerular filtration rate in children by the new revised Schwartz formula and a new generalized formula.

The most widely used formula for estimating glomerular filtration rate (eGFR) in children is the Schwartz formula. It was revised in 2009 using iohexol clearances with measured GFR (mGFR) ranging between 15 and 75 ml/min × 1.73 m(2). Here we assessed the accuracy of the Schwartz formula using the inulin clearance (iGFR) method to evaluate its accuracy for children with less renal impairment comparing 551 iGFRs of 392 children with their Schwartz eGFRs. Serum creatinine was measured using the compensated Jaffe method. In order to find the best relationship between iGFR and eGFR, a linear quadratic regression model was fitted and a more accurate formula was derived. This quadratic formula was: 0.68 × (Height (cm)/serum creatinine (mg/dl))-0.0008 × (height (cm)/serum creatinine (mg/dl))(2)+0.48 × age (years)-(21.53 in males or 25.68 in females). This formula was validated using a split-half cross-validation technique and also externally validated with a new cohort of 127 children. Results show that the Schwartz formula is accurate until a height (Ht)/serum creatinine value of 251, corresponding to an iGFR of 103 ml/min × 1.73 m(2), but significantly unreliable for higher values. For an accuracy of 20 percent, the quadratic formula was significantly better than the Schwartz formula for all patients and for patients with a Ht/serum creatinine of 251 or greater. Thus, the new quadratic formula could replace the revised Schwartz formula, which is accurate for children with moderate renal failure but not for those with less renal impairment or hyperfiltration.

Obesity in Switzerland: a critical assessment of prevalence in children and adults.

The objective was to analyze the situation in Switzerland regarding the prevalence of overweight or obesity in children, adolescents and adults. The data were compared with France, an adjacent much larger country. The results showed that there is a definitive lack of objective information in Switzerland on the prevalence of obesity at different ages. As in other European studies, the fact that many national surveys are classically based on subject interviews (self-reported weights and heights rather than measured values) implies that the overweight/obesity prevalence is largely underestimated in adulthood. For example, in a recent Swiss epidemiological study, the prevalence of obesity (BMI greater than 30 kg/m(2)) averaged 6-7% in young men and women (25-34 y), the prevalence being underestimated by a factor of two to three when body weight was self-reported rather than measured. This phenomenon has already been observed in previous European studies. It is concluded that National Surveys based on telephone interviews generally produce biased obesity prevalence results, although the direction of the changes in prevalence of obesity and its evolution with repeated surveys using strict standardized methodology may be evaluated correctly. Therefore, these surveys should be complemented by large-scale epidemiological studies (based on measured anthropomeric variables rather than declared) covering the different linguistic areas of Switzerland. An epidemiological body weight (BMI) monitoring surveillance system, using a harmonized methodology among European countries, would help to accurately assess differences in obesity prevalence across Europe without methodological bias. It will permit monitoring of the dynamic evolution of obesity prevalence as well as the development of appropriate strategies (taking into account the specificity of each country) for obesity prevention and treatment.

Carnitine deficiency in chronic critical illness.

PURPOSE OF REVIEW: New insight in mitochondrial physiology has highlighted the importance of mitochondrial dysfunction in the metabolic and neuroendocrine changes observed in patients presenting with chronic critical illness. This review highlights specifically the importance of carnitine status in this particular patient population and its impact on beta-oxidation and mitochondrial function. RECENT FINDINGS: The main function of carnitine is long chain fatty acid esterification and transport through the mitochondrial membrane. Carnitine depletion should be suspected in critically ill patients with risk factors such as prolonged continuous renal replacement therapy or chronic parenteral nutrition, and evidence of beta-oxidation impairments such as inappropriate hypertriglyceridemia or hyperlactatemia. When fatty acid oxidation is impaired, acyl-CoAs accumulate and deplete the CoA intramitochondrial pool, hence causing a generalized mitochondrial dysfunction and multiorgan failure, with clinical consequences such as muscle weakness, rhabdomyolysis, cardiomyopathy, arrhythmia or sudden death. In such situations, carnitine plasma levels should be measured along with a complete assessment of plasma amino acid, plasma acylcarnitines and urinary organic acid analysis. Supplementation should be initiated if below normal levels (20 μmol/l) of carnitine are observed. In the absence of current guidelines, we recommend an initial supplementation of 0.5-1 g/day. SUMMARY: Metabolic modifications associated with chronic critical illness are just being explored. Carnitine deficiency in critically ill patients is one aspect of these profound and complex changes associated with prolonged stay in ICU. It is readily measurable in the plasma and can easily be substituted if needed, although guidelines are currently missing.

Using toxicokinetic modeling approach to determine biological reference values (BRVs) and to assess human exposure to pesticides

Exposure to various pesticides has been characterized in workers and the general population, but interpretation and assessment of biomonitoring data from a health risk perspective remains an issue. For workers, a Biological Exposure Index (BEI®) has been proposed for some substances, but most BEIs are based on urinary biomarker concentrations at Threshold Limit Value - Time Weighted Average (TLV-TWA) airborne exposure while occupational exposure can potentially occurs through multiple routes, particularly by skin contact (i.e.captan, chlorpyrifos, malathion). Similarly, several biomonitoring studies have been conducted to assess environmental exposure to pesticides in different populations, but dose estimates or health risks related to these environmental exposures (mainly through the diet), were rarely characterized. Recently, biological reference values (BRVs) in the form of urinary pesticide metabolites have been proposed for both occupationally exposed workers and children. These BRVs were established using toxicokinetic models developed for each substance, and correspond to safe levels of absorption in humans, regardless of the exposure scenario. The purpose of this chapter is to present a review of a toxicokinetic modeling approach used to determine biological reference values. These are then used to facilitate health risk assessments and decision-making on occupational and environmental pesticide exposures. Such models have the ability to link absorbed dose of the parent compound to exposure biomarkers and critical biological effects. To obtain the safest BRVs for the studied population, simulations of exposure scenarios were performed using a conservative reference dose such as a no-observed-effect level (NOEL). The various examples discussed in this chapter show the importance of knowledge on urine collections (i.e. spot samples and complete 8-h, 12-h or 24-h collections), sampling strategies, metabolism, relative proportions of the different metabolites in urine, absorption fraction, route of exposure and background contribution of prior exposures. They also show that relying on urinary measurements of specific metabolites appears more accurate when applying this approach to the case of occupational exposures. Conversely, relying on semi-specific metabolites (metabolites common to a category of pesticides) appears more accurate for the health risk assessment of environmental exposures given that the precise pesticides to which subjects are exposed are often unknown. In conclusion, the modeling approach to define BRVs for the relevant pesticides may be useful for public health authorities for managing issues related to health risks resulting from environmental and occupational exposures to pesticides.

A critical evaluation of secondary cancer risk models applied to Monte Carlo dose distributions of 2-dimensional, 3-dimensional conformal and hybrid intensity-modulated radiation therapy for breast cancer.

The comparison of radiotherapy techniques regarding secondary cancer risk has yielded contradictory results possibly stemming from the many different approaches used to estimate risk. The purpose of this study was to make a comprehensive evaluation of different available risk models applied to detailed whole-body dose distributions computed by Monte Carlo for various breast radiotherapy techniques including conventional open tangents, 3D conformal wedged tangents and hybrid intensity modulated radiation therapy (IMRT). First, organ-specific linear risk models developed by the International Commission on Radiological Protection (ICRP) and the Biological Effects of Ionizing Radiation (BEIR) VII committee were applied to mean doses for remote organs only and all solid organs. Then, different general non-linear risk models were applied to the whole body dose distribution. Finally, organ-specific non-linear risk models for the lung and breast were used to assess the secondary cancer risk for these two specific organs. A total of 32 different calculated absolute risks resulted in a broad range of values (between 0.1% and 48.5%) underlying the large uncertainties in absolute risk calculation. The ratio of risk between two techniques has often been proposed as a more robust assessment of risk than the absolute risk. We found that the ratio of risk between two techniques could also vary substantially considering the different approaches to risk estimation. Sometimes the ratio of risk between two techniques would range between values smaller and larger than one, which then translates into inconsistent results on the potential higher risk of one technique compared to another. We found however that the hybrid IMRT technique resulted in a systematic reduction of risk compared to the other techniques investigated even though the magnitude of this reduction varied substantially with the different approaches investigated. Based on the epidemiological data available, a reasonable approach to risk estimation would be to use organ-specific non-linear risk models applied to the dose distributions of organs within or near the treatment fields (lungs and contralateral breast in the case of breast radiotherapy) as the majority of radiation-induced secondary cancers are found in the beam-bordering regions.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Répercussions de l'agression sur le métabolisme et la composition corporelle chez l'individu obèse : Nutrition de l'obèse agressé. [Repercussions of injury on metabolic and body composition factors in obese individuals : Nutrition and obese critical patients]

La réponse métabolique de l'obèse apparemment « sainen situation d'agression aiguë (polytraumatisés, traumatisés crâniens, patients chirurgicaux, grands brûlés, opérations électives) ne se distingue pas ou peu de celle de l'individu non-obèse. Cependant, les complications médicales liées à l'agression (insuffisances respiratoire et cardiaque, bronchopneumonie, infections de plaies, thrombophlébites et embolies) demeurent plus importantes chez l'obèse morbide que chez l'individu de poids normal. Grâce à l'inflation de ses réserves énergétiques, l'obèse apparemment sain est avantagé, par rapport au sujet mince, au cours d'une agression nutritionnelle chronique telle que le jeûne prolongé. Le facteur fonctionnel limitant la survie dépend avant tout de la composition corporelle initiale et du degré d'adaptation métabolique (et comportementale) en particulier du degré de conservation de la masse maigre par rapport à la masse grasse. La mobilisation accrue de la masse grasse associée à la perte de poids chez l'obèse (par rapport à son homologue non-obèse) est favorable à une prolongation de la vie, car, en brûlant davantage de graisse corporelle, la part des protéines corporelles endogènes utilisée à des fins énergétiques est plus faible. Il s'ensuit chez l'obèse qu'un niveau de masse maigre critique pour la survie n'est atteint qu'après une réduction très marquée de ses réserves énergétiques. En revanche, le sujet mince perd davantage de masse maigre lors de l'amaigrissement et, par conséquent, son métabolisme de repos diminuera plus rapidement que celui du sujet obèse. Cela peut constituer un avantage énergétique évident en termes d'économie d'énergie consécutive à l'adaptation métabolique, mais un inconvénient majeur quant à la durée de la survie. The metabolic response of « apparently healthyobese individuals following acute injury (multiple trauma, head injury and surgical patients, extended burns, elective surgery) is not dramatically different from that of a non-obese individuals. However, the medical complications following the injury (respiratory and cardiac insufficiency, broncho-pneumonia, infections of wounds, trombophlebitis and embolism) are more prevalent in morbid obese patients than in individuals of normal body weight. Because of a large increase in their individuals energy store, "apparently healthy" obese individuals have an advantage over very lean subjects when exposed to a chronic nutritional aggression such as total fasting. The functional limiting factor for survival depends primarily on initial body composition and the magnitude of metabolic adaptation (including behavioral adaptation). The key factor is the extent to which the fat-free mass is maintained (versus to the fat mass) during weight loss. The increased proportion of body fat mobilized during weight loss in obese patients, compared with their non-obese counterparts, favors prolonged survival, because more adipose tissue is burned off, the fraction of body protein endogenously utilized for energy purpose individuals, is smaller. This implies that obese individuals do not reach a fat-free mass "critical" for their survival until their energy stores reach very low values. In contrast, lean subject tend to lose more fat-free mass during weight loss than obese subjects and, as a result, their energy expenditure drops more rapidly. This may offer a potential advantage in terms of energy economy (more energy saving) but a major disadvantage in terms of duration of survival.

PRPF mutations are associated with generalized defects in spliceosome formation and pre-mRNA splicing in patients with retinitis pigmentosa.

Proteins PRPF31, PRPF3 and PRPF8 (RP-PRPFs) are ubiquitously expressed components of the spliceosome, a macromolecular complex that processes nearly all pre-mRNAs. Although these spliceosomal proteins are conserved in eukaryotes and are essential for survival, heterozygous mutations in human RP-PRPF genes lead to retinitis pigmentosa, a hereditary disease restricted to the eye. Using cells from patients with 10 different mutations, we show that all clinically relevant RP-PRPF defects affect the stoichiometry of spliceosomal small nuclear RNAs (snRNAs), the protein composition of tri-small nuclear ribonucleoproteins and the kinetics of spliceosome assembly. These mutations cause inefficient splicing in vitro and affect constitutive splicing ex-vivo by impairing the removal of at least 9% of endogenously expressed introns. Alternative splicing choices are also affected when RP-PRPF defects are present. Furthermore, we show that the steady-state levels of snRNAs and processed pre-mRNAs are highest in the retina, indicating a particularly elevated splicing activity. Our results suggest a role for PRPFs defects in the etiology of PRPF-linked retinitis pigmentosa, which appears to be a truly systemic splicing disease. Although these mutations cause widespread and important splicing defects, they are likely tolerated by the majority of human tissues but are critical for retinal cell survival.

Teaching clinical anatomy of the female pelvic floor to undergraduate students: A critical review of nevralgic points

Pelvic floor anatomy is complex and its three-dimensional organization is often difficult to understand for both undergrad- uate and postgraduate students. Here, we focused on several critical points that need to be considered when teaching the perineum. We have to deal with a mixed population of students and with a variety of interest. Yet, a perfect knowledge of the pelvic floor is the basis for any gynecologist and for any surgical intervention. Our objectives are several-fold; i) to estab- lish the objectives and the best way of teaching, ii) to identify and localize areas in the female pelvic floor that are suscepti- ble to generate problems in understanding the three-dimensional organization, iii) to create novel approaches by respecting the anatomical surroundings, and iv) prospectively, to identify elements that may create problems during surgery i.e. to have a closer look at nerve trajectories and on compression sites that may cause neuralgia or postoperative pain. A feedback from students concludes that they have difficulties to assimilate this much information, especially the different imaging tech- niques. Eventually, this will lead to a severe selection of what has to be taught and included in lectures or practicals. Another consequence is that more time to study prosected pelves needs to be given.

Heritability, determinants and reference values of renal length: a family-based population study.

OBJECTIVES: In this population-based study, reference values were generated for renal length, and the heritability and factors associated with kidney length were assessed. METHODS: Anthropometric parameters and renal ultrasound measurements were assessed in randomly selected nuclear families of European ancestry (Switzerland). The adjusted narrow sense heritability of kidney size parameters was estimated by maximum likelihood assuming multivariate normality after power transformation. Gender-specific reference centiles were generated for renal length according to body height in the subset of non-diabetic non-obese participants with normal renal function. RESULTS: We included 374 men and 419 women (mean ± SD, age 47 ± 18 and 48 ± 17 years, BMI 26.2 ± 4 and 24.5 ± 5 kg/m(2), respectively) from 205 families. Renal length was 11.4 ± 0.8 cm in men and 10.7 ± 0.8 cm in women; there was no difference between right and left renal length. Body height, weight and estimated glomerular filtration rate (eGFR) were positively associated with renal length, kidney function negatively, age quadratically, whereas gender and hypertension were not. The adjusted heritability estimates of renal length and volume were 47.3 ± 8.5 % and 45.5 ± 8.8 %, respectively (P < 0.001). CONCLUSION: The significant heritability of renal length and volume highlights the familial aggregation of this trait, independently of age and body size. Population-based references for renal length provide a useful guide for clinicians. KEY POINTS: • Renal length and volume are heritable traits, independent of age and size. • Based on a European population, gender-specific reference values/percentiles are provided for renal length. • Renal length correlates positively with body length and weight. • There was no difference between right and left renal lengths in this study. • This negates general teaching that the left kidney is larger and longer.

Critical evaluation of outcome scales to assess outcome after lateral ankle ligament repair

Introduction: Several scores are commonly used to evaluate patients' postoperative satisfaction after lateral ankle ligament repair, including: AOFAS, FAAM, CAIT and CAIS. Comparing published studies in the literature is difficult, as the same patient can have markedly different results depending on which scoring system is used. The current study aims to address this gap in the literature by developing a system to compare these tests, to allow better analysis and comparison of published studies. Patients and methods: This is a retrospective cohort study of 47 patients following lateral ankle ligament repair using a modified Broström-Gould technique. All patients were operated between 2005 and 2010 by a single surgeon and followed the same post operative rehabilitation protocol. Six patients were excluded from the study because of concomitant surgery. Patients were assessed by an independent observer. We used the Pearson correlation coefficient to analyse the concordance of the scores, as well as scatter plots to assess the linear relationship between them. Results: A linear distribution between the scores was found when the results were analysed using scatter plots. We were thus able to use the Pearson correlation coefficient to evaluate the relationship between each of the different postoperative scores. The correlation was found to be above 0.5 in all cases except for the comparison between the CAIT and the FAAM for the activities of daily living (0.39). We were, therefore, able to compare the results obtained and assess the relative concordance of the scoring systems. The results showed that the more specific the scale is, the worst the score is and inversely. So the CAIT and the CAIS appeared to be more severe than the AOFAS and the FAAM measuring the activities of daily living. The sports subscale of the FAAM demonstrated intermediate results. Conclusion: This study outlines a system to compare different postoperative scores commonly used to evaluate outcome after ankle stabilization surgery. The impact of this study is that it makes comparison of published studies easier, even though they use a variety of different clinical scores, thus facilitating better outcome analysis of operative techniques.

Pharmacokinetics and safety of panobacumab: specific adjunctive immunotherapy in critical patients with nosocomial Pseudomonas aeruginosa O11 pneumonia.

Objectives Nosocomial Pseudomonas aeruginosa pneumonia remains a major concern in critically ill patients. We explored the potential impact of microorganism-targeted adjunctive immunotherapy in such patients. Patients and methods This multicentre, open pilot Phase 2a clinical trial (NCT00851435) prospectively evaluated the safety, pharmacokinetics and potential efficacy of three doses of 1.2 mg/kg panobacumab, a fully human monoclonal anti-lipopolysaccharide IgM, given every 72 h in 18 patients developing nosocomial P. aeruginosa (serotype O11) pneumonia. Results Seventeen out of 18 patients were included in the pharmacokinetic analysis. In 13 patients receiving three doses, the maximal concentration after the third infusion was 33.9 ± 8.0 μg/mL, total area under the serum concentration-time curve was 5397 ± 1993 μg h/mL and elimination half-life was 102.3 ± 47.8 h. Panobacumab was well tolerated, induced no immunogenicity and was detected in respiratory samples. In contrast to Acute Physiology and Chronic Health Evaluation II (APACHE II) prediction, all 13 patients receiving three doses survived, with a mean clinical resolution in 9.0 ± 2.7 days. Two patients suffered a recurrence at days 17 and 20. Conclusions These data suggest that panobacumab is safe, with a pharmacokinetic profile similar to that in healthy volunteers. It was associated with high clinical cure and survival rates in patients developing nosocomial P. aeruginosa O11 pneumonia. We concluded that these promising results warrant further trials.

Development and current use of parenteral nutrition in critical care - an opinion paper.

Critically ill patients depend on artificial nutrition for the maintenance of their metabolic functions and lean body mass, as well as for limiting underfeeding-related complications. Current guidelines recommend enteral nutrition (EN), possibly within the first 48 hours, as the best way to provide the nutrients and prevent infections. EN may be difficult to realize or may be contraindicated in some patients, such as those presenting anatomic intestinal continuity problems or splanchnic ischemia. A series of contradictory trials regarding the best route and timing for feeding have left the medical community with great uncertainty regarding the place of parenteral nutrition (PN) in critically ill patients. Many of the deleterious effects attributed to PN result from inadequate indications, or from overfeeding. The latter is due firstly to the easier delivery of nutrients by PN compared with EN increasing the risk of overfeeding, and secondly to the use of approximate energy targets, generally based on predictive equations: these equations are static and inaccurate in about 70% of patients. Such high uncertainty about requirements compromises attempts at conducting nutrition trials without indirect calorimetry support because the results cannot be trusted; indeed, both underfeeding and overfeeding are equally deleterious. An individualized therapy is required. A pragmatic approach to feeding is proposed: at first to attempt EN whenever and as early as possible, then to use indirect calorimetry if available, and to monitor delivery and response to feeding, and finally to consider the option of combining EN with PN in case of insufficient EN from day 4 onwards.