Plasma and cerebrospinal fluid population pharmacokinetics of temozolomide in malignant glioma patients

RESUME : Objectif: Le glioblastome multiforme (GBM) est la tumeur cérébrale maligne la plus agressive qui conduit au décès de la majorité des patients moins d'une année après le diagnostic. La plupart des agents chimiothérapeutiques actuellement disponibles ne traversent pas la barrière hémato¬encéphalique et ne peuvent par conséquent pas être utilisés pour ce type de tumeur. Le Temozolomide (TMZ) est un nouvel agent alkylant récemment développé pour le traitement des gliomes malins. A ce jour, très peu d'informations sont disponibles sur la pénétration intra-cérébrale de cet agent. Au cours d'une étude pilote de phase II menée auprès de 64 patients atteints de GBM, l'administration précoce de TMZ combinée à une radiothérapie standard (RT) afin d'intervenir au plus tôt dans l'évolution de la maladie, a permis de prolonger la survie de ces patients, résultat qui pu être confirmé par la suite lors de l'étude randomisée de phase III. L'objectif de cette étude a été de déterminer les paramètres pharmacocinétique du TMZ dans le plasma et le liquide céphalo-rachidien (LCR), d'évaluer l'influence de certains facteurs individuels (âge, sexe, surface corporelle, fonction rénale/hépatique, co-médications, RT concomitante) sur ces différents paramètres, et enfin d'explorer la relation existant entre l'exposition au TMZ et certains marqueurs cliniques d'efficacité et de toxicité. Matériel et Méthode: Les concentrations de TMZ ont été mesurées par chromatographie liquide à haute performance (HPLC) dans le plasma et le LCR de 35 patients atteints de GBM nouvellement diagnostiqués (étude pilote) ou de gliomes malins en récidive (étude récidive). L'analyse pharmacocinétique de population a été réalisée à l'aide du programme NONMEM. L'exposition systémique et cérébrale, définie par les AUC (Area Under the time-concentration Curve) dans le plasma et le LCR, a été estimée pour chaque patient et corrélée à la toxicité, la survie ainsi que la survie sans progression tumorale. Résultats: Un modèle à 1 compartiment avec une cinétique d'absorption et de transfert Kplasma -> LCR de ordre a été retenu afin de décrire le profil pharmacocinétique du TMZ. Les valeurs moyennes de population ont été de 10 L/h pour la clairance, de 30.3 L pour le volume de distribution, de 2.1 h pour la 1/2 vie d'élimination, de 5.78 hE-1 pour la constante d'absorption, de 7.2 10E4 hE-1 pour Kplasma->LCR et de 0.76 hE-1 pour KLCR plasma. La surface corporelle a montré une influence significative sur la clairance et le volume de distribution, alors que le sexe influence la clairance uniquement. L'AUC mesurée dans le LCR représente ~20% de celle du plasma et une augmentation de 15% de Kplasma->LCR a été observée lors du traitement concomitant de radiochimiothérapie. Conclusions: Cette étude est la première analyse pharmacocinétique effectuée chez l'homme permettant de quantifier la pénétration intra-cérébrale du TMZ. Le rapport AUC LCR/AUC Plasma a été de 20%. Le degré d'exposition systémique et cérébral au TMZ ne semble pas être un meilleur facteur prédictif de la survie ou de la tolérance au produit que ne l'est la dose cumulée seule. ABSTRACT Purpose: Scarce information is available on the brain penetration of temozolomide (TMZ), although this novel methylating agent is mainly used for the treatment of ma¬lignant brain tumors. The purpose was to assess TNIZ phar¬macokinetics in plasma and cerebrospinal fluid (CSF) along with its inter-individual variability, to characterize covari¬ates and to explore relationships between systemic or cere¬bral drug exposure and clinical outcomes. Experimental Design: TMZ levels were measured by high-performance liquid chromatography in plasma and CSF samples from 35 patients with newly diagnosed or recurrent malignant gliomas. The population pharmacoki¬netic analysis was performed with nonlinear mixed-effect modeling software. Drug exposure, defined by the area un¬der the concentration-time curve (AUC) in plasma and CSF, was estimated for each patient and correlated with toxicity, survival, and progression-free survival. Results: A three-compartment model with first-order absorption and transfer rates between plasma and CSF described the data appropriately. Oral clearance was 10 liter/h; volume of distribution (VD), 30.3 liters; absorption constant rate, 5.8 hE-1; elimination half-time, 2.1 h; transfer rate from plasma to CSF (Kplasma->CSF), 7.2 x 10E-4hE-1 and the backwards rate, 0.76hE-1. Body surface area signifi¬cantly influenced both clearance and VD, and clearance was sex dependent. The AU CSF corresponded to 20% of the AUCplasma. A trend toward an increased K plasma->CSF of 15% was observed in case of concomitant radiochemo-therapy. No significant correlations between AUC in plasma or CSF and toxicity, survival, or progression-free survival were apparent after deduction of dose-effect. Conclusions: This is the first human pharmacokinetic study on TMZ to quantify CSF penetration. The AUC CSF/ AUC plasma ratio was 20%. Systemic or cerebral exposures are not better predictors than the cumulative dose alone for both efficacy and safety.

Structural alterations of the coronary arterial wall are associated with myocardial flow heterogeneity in type 2 diabetes mellitus.

PURPOSE: To determine the relationship between carotid intima-media thickness (IMT), coronary artery calcification (CAC), and myocardial blood flow (MBF) at rest and during vasomotor stress in type 2 diabetes mellitus (DM). METHODS: In 68 individuals, carotid IMT was measured using high-resolution vascular ultrasound, while the presence of CAC was determined with electron beam tomography (EBT). Global and regional MBF was determined in milliliters per gram per minute with (13)N-ammonia and positron emission tomography (PET) at rest, during cold pressor testing (CPT), and during adenosine (ADO) stimulation. RESULTS: There was neither a relationship between carotid IMT and CAC (r = 0.10, p = 0.32) nor between carotid IMT and coronary circulatory function in response to CPT and during ADO (r = -0.18, p = 0.25 and r = 0.10, p = 0.54, respectively). In 33 individuals, EBT detected CAC with a mean Agatston-derived calcium score of 44 +/- 18. There was a significant difference in regional MBFs between territories with and without CAC at rest and during ADO-stimulated hyperemia (0.69 +/- 0.24 vs. 0.74 +/- 0.23 and 1.82 +/- 0.50 vs. 1.95 +/- 0.51 ml/g/min; p < or = 0.05, respectively) and also during CPT in DM but less pronounced (0.81 +/- 0.24 vs. 0.83 +/- 0.23 ml/g/min; p = ns). The increase in CAC was paralleled with a progressive regional decrease in resting as well as in CPT- and ADO-related MBFs (r = -0.36, p < or = 0.014; r = -0.46, p < or = 0.007; and r = -0.33, p < or = 0.041, respectively). CONCLUSIONS: The absence of any correlation between carotid IMT and coronary circulatory function in type 2 DM suggests different features and stages of early atherosclerosis in the peripheral and coronary circulation. PET-measured MBF heterogeneity at rest and during vasomotor stress may reflect downstream fluid dynamic effects of coronary artery disease (CAD)-related early structural alterations of the arterial wall.

Vasopressin-dependent coupling between sodium transport and water flow in a mouse cortical collecting duct cell line.

Water balance is achieved through the ability of the kidney to control water reabsorption in the connecting tubule and the collecting duct. In a mouse cortical collecting duct cell line (mCCD(c11)), physiological concentrations of arginine vasopressin increased both electrogenic, amiloride-sensitive, epithelial sodium channel (ENaC)-mediated sodium transport measured by the short-circuit current (Isc) method and water flow (Jv apical to basal) measured by gravimetry with similar activation coefficient K(1/2) (6 and 12 pM, respectively). Jv increased linearly according to the osmotic gradient across the monolayer. A small but highly significant Jv was also measured under isoosmotic conditions. To test the coupling between sodium reabsorption and water flow, mCCD(c11) cells were treated for 24 h under isoosmotic condition with either diluent, amiloride, vasopressin or vasopressin and amiloride. Isc, Jv, and net chemical sodium fluxes were measured across the same monolayers. Around 30% of baseline and 50% of vasopressin-induced water flow is coupled to an amiloride-sensitive, ENaC-mediated, electrogenic sodium transport, whereas the remaining flow is coupled to an amiloride-insensitive, nonelectrogenic sodium transport mediated by an unknown electroneutral transporter. The mCCD(c11) cell line is a first example of a mammalian tight epithelium allowing quantitative study of the coupling between sodium and water transport. Our data are consistent with the 'near isoosmotic' fluid transport model.

Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

Comment on recent modeling studies of astrocyte-neuron metabolic interactions.

Recent years have seen a surge in mathematical modeling of the various aspects of neuron-astrocyte interactions, and the field of brain energy metabolism is no exception in that regard. Despite the advent of biophysical models in the field, the long-lasting debate on the role of lactate in brain energy metabolism is still unresolved. Quite the contrary, it has been ported to the world of differential equations. Here, we summarize the present state of this discussion from the modeler's point of view and bring some crucial points to the attention of the non-mathematically proficient reader.

Modeling complex wells with the multi-scale finite-volume method

In this paper, an extension of the multi-scale finite-volume (MSFV) method is devised, which allows to Simulate flow and transport in reservoirs with complex well configurations. The new framework fits nicely into the data Structure of the original MSFV method,and has the important property that large patches covering the whole well are not required. For each well. an additional degree of freedom is introduced. While the treatment of pressure-constraint wells is trivial (the well-bore reference pressure is explicitly specified), additional equations have to be solved to obtain the unknown well-bore pressure of rate-constraint wells. Numerical Simulations of test cases with multiple complex wells demonstrate the ability of the new algorithm to capture the interference between the various wells and the reservoir accurately. (c) 2008 Elsevier Inc. All rights reserved.

Carbon and oxygen isotope study of hydrothermal carbonates in the zinc-lead deposits of the San Vicente district, central Peru: A quantitative modeling on mixing processes and CO2 degassing

Carbon and oxygen isotope studies of the host and gangue carbonates of Mississippi Valley-type zinc-lead deposits in the San Vicente District hosted in the Upper Triassic to Lower Jurassic dolostones of the Pucara basin (central Peru) were used to constrain models of the ore formation. A mixing model between an incoming hot saline slightly acidic radiogenic (Pb, Sr) fluid and the native formation water explains the overall isotopic variation (delta(13)C = - 11.5 to + 2.5 parts per thousand relative to PDB and delta(18)O = + 18.0 to + 24.3 parts per thousand relative to SMOW) of the carbonate generations. The dolomites formed during the main ore stage show a narrower range (delta(13)C = - 0.1 to + 1.7 parts per thousand and delta(18)O = + 18.7 to + 23.4 parts per thousand) which is explained by exchange between the mineralizing fluids and the host carbonates combined with changes in temperature and pressure. This model of fluid-rock interaction explains the pervasive alteration of the host dolomite I and precipitation of sphalerite I. The open-space filling hydrothermal white sparry dolomite and the coexisting sphalerite II formed by prolonged fluid-host dolomite interaction and limited CO2 degassing. Late void-filling dolomite III (or calcite) and the associated sphalerite III formed as the consequence of CO2 degassing and concomitant pH increase of a slightly acidic ore fluid. Widespread brecciation is associated to CO2 outgassing. Consequently, pressure variability plays a major role in the ore precipitation during the late hydrothermal events in San Vicente. The presence of native sulfur associated with extremely carbon-light calcites replacing evaporitic sulfates (e.g., delta(13)C = - 11.5 parts per thousand), altered native organic matter and heavier hydrothermal bitumen (from - 27.0 to - 23.0 parts per thousand delta(13)C) points to thermochemical reduction of sulfate and/or thiosulfate. The delta(13)C- and delta(18)O-values of the altered host dolostone and hydrothermal carbonates, and the carbon isotope composition of the associated organic matter show a strong regional homogeneity. These results coupled with the strong mineralogical and petrographic similarities of the different MVT occurrences perhaps reflects the fact that the mineralizing processes were similar in the whole San Vicente belt, suggesting the existence of a common regional mineralizing hydrothermal system with interconnected plumbing.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Properties of debris flow deposits and source materials compared: implications for debris flow characterization

In the Alps, debris flow deposits generally contain < 5% clay-size particles, and the role of the surface-charged < 2 mu m particles is often neglected, although these particles may have a significant impact on the rheological properties of the interstitial fluid. The objective of this study was to compare debris flow deposits and parent materials from two neighbouring catchments of the Swiss Alps, with special emphasis on the colloidal constituents. The catchments are small in area (4 km(2)), 2.5 km long, similar in morphology, but different in geology. The average slopes are 35-40%. The catchments were monitored for debris flow events and mapped for surface aspect and erosion activity. Debris flow deposits and parent materials were sampled, the clay and silt fractions extracted and the bulk density, < 2 mm fraction bulk density, particle size distribution, chemical composition, cation exchange capacity (CEC) and mineralogy analysed. The results show that the deposits are similar to the parent screes in terms of chemical composition, but differ in terms of: (i) particle size distribution; and (ii) mineralogy, reactivity and density of the < 2 mm fraction. In this fraction, compared with the parent materials the deposits show dense materials enriched in coarse monocrystalline particles, of which the smallest and more reactive particles were leached. The results suggest that deposit samples should not be considered as representative of source or flow materials, particularly with respect to their physical properties.

MORM--a Petri net based model for assessing OH&S risks in industrial processes: modeling qualitative aspects.

Because of the increase in workplace automation and the diversification of industrial processes, workplaces have become more and more complex. The classical approaches used to address workplace hazard concerns, such as checklists or sequence models, are, therefore, of limited use in such complex systems. Moreover, because of the multifaceted nature of workplaces, the use of single-oriented methods, such as AEA (man oriented), FMEA (system oriented), or HAZOP (process oriented), is not satisfactory. The use of a dynamic modeling approach in order to allow multiple-oriented analyses may constitute an alternative to overcome this limitation. The qualitative modeling aspects of the MORM (man-machine occupational risk modeling) model are discussed in this article. The model, realized on an object-oriented Petri net tool (CO-OPN), has been developed to simulate and analyze industrial processes in an OH&S perspective. The industrial process is modeled as a set of interconnected subnets (state spaces), which describe its constitutive machines. Process-related factors are introduced, in an explicit way, through machine interconnections and flow properties. While man-machine interactions are modeled as triggering events for the state spaces of the machines, the CREAM cognitive behavior model is used in order to establish the relevant triggering events. In the CO-OPN formalism, the model is expressed as a set of interconnected CO-OPN objects defined over data types expressing the measure attached to the flow of entities transiting through the machines. Constraints on the measures assigned to these entities are used to determine the state changes in each machine. Interconnecting machines implies the composition of such flow and consequently the interconnection of the measure constraints. This is reflected by the construction of constraint enrichment hierarchies, which can be used for simulation and analysis optimization in a clear mathematical framework. The use of Petri nets to perform multiple-oriented analysis opens perspectives in the field of industrial risk management. It may significantly reduce the duration of the assessment process. But, most of all, it opens perspectives in the field of risk comparisons and integrated risk management. Moreover, because of the generic nature of the model and tool used, the same concepts and patterns may be used to model a wide range of systems and application fields.

Budgeting rockfall and modeling sedimentary delivery in torrent systems

This thesis is a compilation of projects to study sediment processes recharging debris flow channels. These works, conducted during my stay at the University of Lausanne, focus in the geological and morphological implications of torrent catchments to characterize debris supply, a fundamental element to predict debris flows. Other aspects of sediment dynamics are considered, e.g. the coupling headwaters - torrent, as well as the development of a modeling software that simulates sediment transfer in torrent systems. The sediment activity at Manival, an active torrent system of the northern French Alps, was investigated using terrestrial laser scanning and supplemented with geostructural investigations and a survey of sediment transferred in the main torrent. A full year of sediment flux could be observed, which coincided with two debris flows and several bedload transport events. This study revealed that both debris flows generated in the torrent and were preceded in time by recharge of material from the headwaters. Debris production occurred mostly during winter - early spring time and was caused by large slope failures. Sediment transfers were more puzzling, occurring almost exclusively in early spring subordinated to runoffconditions and in autumn during long rainfall. Intense rainstorms in summer did not affect debris storage that seems to rely on the stability of debris deposits. The morpho-geological implication in debris supply was evaluated using DEM and field surveys. A slope angle-based classification of topography could characterize the mode of debris production and transfer. A slope stability analysis derived from the structures in rock mass could assess susceptibility to failure. The modeled rockfall source areas included more than 97% of the recorded events and the sediment budgets appeared to be correlated to the density of potential slope failure. This work showed that the analysis of process-related terrain morphology and of susceptibility to slope failure document the sediment dynamics to quantitatively assess erosion zones leading to debris flow activity. The development of erosional landforms was evaluated by analyzing their geometry with the orientations of potential rock slope failure and with the direction of the maximum joint frequency. Structure in rock mass, but in particular wedge failure and the dominant discontinuities, appear as a first-order control of erosional mechanisms affecting bedrock- dominated catchment. They represent some weaknesses that are exploited primarily by mass wasting processes and erosion, promoting not only the initiation of rock couloirs and gullies, but also their propagation. Incorporating the geological control in geomorphic processes contributes to better understand the landscape evolution of active catchments. A sediment flux algorithm was implemented in a sediment cascade model that discretizes the torrent catchment in channel reaches and individual process-response systems. Each conceptual element includes in simple manner geomorphological and sediment flux information derived from GIS complemented with field mapping. This tool enables to simulate sediment transfers in channels considering evolving debris supply and conveyance, and helps reducing the uncertainty inherent to sediment budget prediction in torrent systems. Cette thèse est un recueil de projets d'études des processus de recharges sédimentaires des chenaux torrentiels. Ces travaux, réalisés lorsque j'étais employé à l'Université de Lausanne, se concentrent sur les implications géologiques et morphologiques des bassins dans l'apport de sédiments, élément fondamental dans la prédiction de laves torrentielles. D'autres aspects de dynamique sédimentaire ont été abordés, p. ex. le couplage torrent - bassin, ainsi qu'un modèle de simulation du transfert sédimentaire en milieu torrentiel. L'activité sédimentaire du Manival, un système torrentiel actif des Alpes françaises, a été étudiée par relevés au laser scanner terrestre et complétée par une étude géostructurale ainsi qu'un suivi du transfert en sédiments du torrent. Une année de flux sédimentaire a pu être observée, coïncidant avec deux laves torrentielles et plusieurs phénomènes de charriages. Cette étude a révélé que les laves s'étaient générées dans le torrent et étaient précédées par une recharge de débris depuis les versants. La production de débris s'est passée principalement en l'hiver - début du printemps, causée par de grandes ruptures de pentes. Le transfert était plus étrange, se produisant presque exclusivement au début du printemps subordonné aux conditions d'écoulement et en automne lors de longues pluies. Les orages d'été n'affectèrent guère les dépôts, qui semblent dépendre de leur stabilité. Les implications morpho-géologiques dans l'apport sédimentaire ont été évaluées à l'aide de MNT et études de terrain. Une classification de la topographie basée sur la pente a permis de charactériser le mode de production et transfert. Une analyse de stabilité de pente à partir des structures de roches a permis d'estimer la susceptibilité à la rupture. Les zones sources modélisées comprennent plus de 97% des chutes de blocs observées et les bilans sédimentaires sont corrélés à la densité de ruptures potentielles. Ce travail d'analyses des morphologies du terrain et de susceptibilité à la rupture documente la dynamique sédimentaire pour l'estimation quantitative des zones érosives induisant l'activité torrentielle. Le développement des formes d'érosion a été évalué par l'analyse de leur géométrie avec celle des ruptures potentielles et avec la direction de la fréquence maximale des joints. Les structures de roches, mais en particulier les dièdres et les discontinuités dominantes, semblent être très influents dans les mécanismes d'érosion affectant les bassins rocheux. Ils représentent des zones de faiblesse exploitées en priorité par les processus de démantèlement et d'érosion, encourageant l'initiation de ravines et couloirs, mais aussi leur propagation. L'incorporation du control géologique dans les processus de surface contribue à une meilleure compréhension de l'évolution topographique de bassins actifs. Un algorithme de flux sédimentaire a été implémenté dans un modèle en cascade, lequel divise le bassin en biefs et en systèmes individuels répondant aux processus. Chaque unité inclut de façon simple les informations géomorpologiques et celles du flux sédimentaire dérivées à partir de SIG et de cartographie de terrain. Cet outil permet la simulation des transferts de masse dans les chenaux, considérants la variabilité de l'apport et son transport, et aide à réduire l'incertitude liée à la prédiction de bilans sédimentaires torrentiels. Ce travail vise très humblement d'éclairer quelques aspects de la dynamique sédimentaire en milieu torrentiel.

Flow structures at an idealized bifurcation : a numerical experiment

River bifurcations are key nodes within braided river systems controlling the flow and sediment partitioning and therefore the dynamics of the river braiding process. Recent research has shown that certain geometrical configurations induce instabilities that lead to downstream mid-channel bar formation and the formation of bifurcations. However, we currently have a poor understanding of the flow division process within bifurcations and the flow dynamics in the downstream bifurcates, both of which are needed to understand bifurcation stability. This paper presents results of a numerical sensitivity experiment undertaken using computational fluid dynamics (CFD) with the purpose of understanding the flow dynamics of a series of idealized bifurcations. A geometric sensitivity analysis is undertaken for a range of channel slopes (0.005 to 0.03), bifurcation angles (22 degrees to 42 degrees) and a restricted set of inflow conditions based upon simulating flow through meander bends with different curvature on the flow field dynamics through the bifurcation. The results demonstrate that the overall slope of the bifurcation affects the velocity of flow through the bifurcation and when slope asymmetry is introduced, the flow structures in the bifurcation are modified. In terms of bifurcation evolution the most important observation appears to be that once slope asymmetry is greater than 0.2 the flow within the steep bifurcate shows potential instability and the potential for alternate channel bar formation. Bifurcation angle also defines the flow structures within the bifurcation with an increase in bifurcation angle increasing the flow velocity down both bifurcates. However, redistributive effects of secondary circulation caused by upstream curvature can very easily counter the effects of local bifurcation characteristics. Copyright (C) 2011 John Wiley & Sons, Ltd.

Measurements of glial metabolic fluxes with C-11-acetate using positron emission and an adapted NMR-based metabolic modeling approach

Objectives: Acetate brain metabolism has the particularity to occur specifically in glial cells. Labeling studies, using acetate labeled either with 13C (NMR) or 11C (PET), are governed by the same biochemical reactions and thus follow the same mathematical principles. In this study, the objective was to adapt an NMR acetate brain metabolism model to analyse [1-11C]acetate infusion in rats. Methods: Brain acetate infusion experiments were modeled using a two-compartment model approach used in NMR.1-3 The [1-11C]acetate labeling study was done using a beta scintillator.4 The measured radioactive signal represents the time evolution of the sum of all labeled metabolites in the brain. Using a coincidence counter in parallel, an arterial input curve was measured. The 11C at position C-1 of acetate is metabolized in the first turn of the TCA cycle to the position 5 of glutamate (Figure 1A). Through the neurotransmission process, it is further transported to the position 5 of glutamine and the position 5 of neuronal glutamate. After the second turn of the TCA cycle, tracer from [1-11C]acetate (and also a part from glial [5-11C]glutamate) is transferred to glial [1-11C]glutamate and further to [1-11C]glutamine and neuronal glutamate through the neurotransmission cycle. Brain poster session: oxidative mechanisms S460 Journal of Cerebral Blood Flow & Metabolism (2009) 29, S455-S466 Results: The standard acetate two-pool PET model describes the system by a plasma pool and a tissue pool linked by rate constants. Experimental data are not fully described with only one tissue compartment (Figure 1B). The modified NMR model was fitted successfully to tissue time-activity curves from 6 single animals, by varying the glial mitochondrial fluxes and the neurotransmission flux Vnt. A glial composite rate constant Kgtg=Vgtg/[Ace]plasma was extracted. Considering an average acetate concentration in plasma of 1 mmol/g5 and the negligible additional amount injected, we found an average Vgtg = 0.08±0.02 (n = 6), in agreement with previous NMR measurements.1 The tissue time-activity curve is dominated by glial glutamate and later by glutamine (Figure 1B). Labeling of neuronal pools has a low influence, at least for the 20 mins of beta-probe acquisition. Based on the high diffusivity of CO2 across the blood-brain barrier; 11CO2 is not predominant in the total tissue curve, even if the brain CO2 pool is big compared with other metabolites, due to its strong dilution through unlabeled CO2 from neuronal metabolism and diffusion from plasma. Conclusion: The two-compartment model presented here is also able to fit data of positron emission experiments and to extract specific glial metabolic fluxes. 11C-labeled acetate presents an alternative for faster measurements of glial oxidative metabolism compared to NMR, potentially applicable to human PET imaging. However, to quantify the relative value of the TCA cycle flux compared to the transmitochondrial flux, the chemical sensitivity of NMR is required. PET and NMR are thus complementary.

Brain lactate kinetics: Modeling evidence for neuronal lactate uptake upon activation.

A critical issue in brain energy metabolism is whether lactate produced within the brain by astrocytes is taken up and metabolized by neurons upon activation. Although there is ample evidence that neurons can efficiently use lactate as an energy substrate, at least in vitro, few experimental data exist to indicate that it is indeed the case in vivo. To address this question, we used a modeling approach to determine which mechanisms are necessary to explain typical brain lactate kinetics observed upon activation. On the basis of a previously validated model that takes into account the compartmentalization of energy metabolism, we developed a mathematical model of brain lactate kinetics, which was applied to published data describing the changes in extracellular lactate levels upon activation. Results show that the initial dip in the extracellular lactate concentration observed at the onset of stimulation can only be satisfactorily explained by a rapid uptake within an intraparenchymal cellular compartment. In contrast, neither blood flow increase, nor extracellular pH variation can be major causes of the lactate initial dip, whereas tissue lactate diffusion only tends to reduce its amplitude. The kinetic properties of monocarboxylate transporter isoforms strongly suggest that neurons represent the most likely compartment for activation-induced lactate uptake and that neuronal lactate utilization occurring early after activation onset is responsible for the initial dip in brain lactate levels observed in both animals and humans.