Identification of pharmaceutical tablets by Raman spectroscopy and chemometrics

Raman spectroscopy has become an attractive tool for the analysis of pharmaceutical solid dosage forms. In the present study it is used to ensure the identity of tablets. The two main applications of this method are release of final products in quality control and detection of counterfeits. Twenty-five product families of tablets have been included in the spectral library and a non-linear classification method, the Support Vector Machines (SVMs), has been employed. Two calibrations have been developed in cascade: the first one identifies the product family while the second one specifies the formulation. A product family comprises different formulations that have the same active pharmaceutical ingredient (API) but in a different amount. Once the tablets have been classified by the SVM model, API peaks detection and correlation are applied in order to have a specific method for the identification and allow in the future to discriminate counterfeits from genuine products. This calibration strategy enables the identification of 25 product families without error and in the absence of prior information about the sample. Raman spectroscopy coupled with chemometrics is therefore a fast and accurate tool for the identification of pharmaceutical tablets.

Population pharmacokinetic study of gentamicin in a large cohort of premature and term neonates.

AIM: This study aims to investigate the clinical and demographic factors influencing gentamicin pharmacokinetics in a large cohort of unselected premature and term newborns and to evaluate optimal regimens in this population. METHODS: All gentamicin concentration data, along with clinical and demographic characteristics, were retrieved from medical charts in a Neonatal Intensive Care Unit over 5 years within the frame of a routine therapeutic drug monitoring programme. Data were described using non-linear mixed-effects regression analysis ( nonmem®). RESULTS: A total of 3039 gentamicin concentrations collected in 994 preterm and 455 term newborns were included in the analysis. A two compartment model best characterized gentamicin disposition. The average parameter estimates, for a median body weight of 2170 g, were clearance (CL) 0.089 l h(-1) (CV 28%), central volume of distribution (Vc ) 0.908 l (CV 18%), intercompartmental clearance (Q) 0.157 l h(-1) and peripheral volume of distribution (Vp ) 0.560 l. Body weight, gestational age and post-natal age positively influenced CL. Dopamine co-administration had a significant negative effect on CL, whereas the influence of indomethacin and furosemide was not significant. Both body weight and gestational age significantly influenced Vc . Model-based simulations confirmed that, compared with term neonates, preterm infants need higher doses, superior to 4 mg kg(-1) , at extended intervals to achieve adequate concentrations. CONCLUSIONS: This observational study conducted in a large cohort of newborns confirms the importance of body weight and gestational age for dosage adjustment. The model will serve to set up dosing recommendations and elaborate a Bayesian tool for dosage individualization based on concentration monitoring.

3-d residual eddy current field characterisation: applied to diffusion weighted magnetic resonance imaging.

Clinical use of the Stejskal-Tanner diffusion weighted images is hampered by the geometric distortions that result from the large residual 3-D eddy current field induced. In this work, we aimed to predict, using linear response theory, the residual 3-D eddy current field required for geometric distortion correction based on phantom eddy current field measurements. The predicted 3-D eddy current field induced by the diffusion-weighting gradients was able to reduce the root mean square error of the residual eddy current field to ~1 Hz. The model's performance was tested on diffusion weighted images of four normal volunteers, following distortion correction, the quality of the Stejskal-Tanner diffusion-weighted images was found to have comparable quality to image registration based corrections (FSL) at low b-values. Unlike registration techniques the correction was not hindered by low SNR at high b-values, and results in improved image quality relative to FSL. Characterization of the 3-D eddy current field with linear response theory enables the prediction of the 3-D eddy current field required to correct eddy current induced geometric distortions for a wide range of clinical and high b-value protocols.

Ganciclovir exposure under a 450mg daily dosage of valganciclovir for cytomegalovirus prevention in kidney transplantation: a prospective study.

Manuel O, Pascual M, Perrottet N, Lamoth F, Venetz J-P, Decosterd LA, Buclin T, Meylan PR. Ganciclovir exposure under a 450 mg daily dosage of valganciclovir for cytomegalovirus prevention in kidney transplantation: a prospective study. 
Clin Transplant 2010: 24: 794-800. Abstract:  This prospective study aimed at determining the ganciclovir exposure observed under a daily dosage of 450 mg valganciclovir routinely applied to kidney transplant recipients with a GFR above 25 mL/min at risk for cytomegalovirus (CMV) disease. Ganciclovir levels at trough (C(trough) ) and at peak (C(3 h) ) were measured monthly. Ganciclovir exposure (area under the curve [AUC(0-24) ]) was estimated using Bayesian non-linear mixed-effect modeling (NONMEM). Thirty-six patients received 450 mg of valganciclovir daily for three months. Median ganciclovir C(3 h) was 3.9 mg/L (range: 1.3-7.1), and C(trough) was 0.4 mg/L (range 0.1-2.7). Median AUC(0-24) of ganciclovir was 59.3 mg h/L (39.0-85.3) in patients with GFR(MDRD) 26-39 mL/min, 35.8 mg h/L (24.9-55.8) in patients with GFR(MDRD) 40-59 mL/min, and 29.6 mg h/L (22.0-43.2) in patients with GFR(MDRD)  ≥ 60 mL/min. No major differences in adverse events according to ganciclovir exposure were observed. CMV viremia was not detected during prophylaxis. After discontinuing prophylaxis, CMV viremia was seen in 8/36 patients (22%), and 4/36 patients (11%) developed CMV disease. Ganciclovir exposure after administration of valganciclovir 450 mg daily in recipients with GFR ≥60 mL/min was comparable to those previously reported with oral ganciclovir. A routine daily dose of 450 mg valganciclovir appears to be acceptable for CMV prophylaxis in most kidney transplant recipients.

Imatinib plasma concentrations variability in oncologic patients

Imatinib (Glivec®) has transformed the treatment and short-term prognosis of chronic myeloid leukemia (CML) and gastrointestinal stromal tumor (GIST). However, the treatment must be taken indefinitely, it is not devoid of inconvenience and toxicity. Moreover, resistance or escape from disease control occurs in a significant number of patients. Imatinib is a substrate of the cytochromes P450 CYP3A4/5 and of the multidrug transporter P-glycoprotein (product of the MDR1 gene). Considering the large inter-individual differences in the expression and function of those systems, the disposition and clinical activity of imatinib can be expected to vary widely among patients, calling for dosage individualization. The aim of this exploratory study was to determine the average pharmacokinetic parameters characterizing the disposition of imatinib in the target population, to assess their inter-individual variability, and to identify influential factors affecting them. A total of 321 plasma concentrations, taken at various sampling times after the latest dose, were measured in 59 patients receiving Glivec at diverse regimens, using a validated HPLC-UV method developed for this study. The results were analyzed by non-linear mixed effect modeling (NONMEM). A one-compartment model with first-order absorption appeared appropriate to describe the data, with an average apparent clearance of 12.4 l/h, a distribution volume of 268 l and an absorption constant of 0.47 h-1. The clearance was affected by body weight, age and sex. No influences of interacting drugs were found. DNA samples were used for pharmacogenetic explorations. At present, only the MDR1 polymorphism has been assessed and seems to affect the pharmacokinetic parameters of imatinib. Large inter-individual variability remained unexplained by the demographic covariates considered, both on clearance (40 %) and distribution volume (71 %). Together with intra-patient variability (34 %), this translates into an 8-fold width of the 90 %-prediction interval of plasma concentrations expected under a fixed dosing regimen. This is a strong argument to further investigate the possible usefulness of a therapeutic drug monitoring program for imatinib. It may help to individualize the dosing regimen before overt disease progression or observation of treatment toxicity, thus improving both the long-term therapeutic effectiveness and tolerability of this drug.

Population pharmacokinetics of imatinib in CML and GIST patients under long-term treatment

Imatinib (Glivec®) has transformed the treatment and short-term prognosis of chronic myeloid leukaemia (CML) and gastro-intestinal stromal tumour (GIST). However, the treatment must be taken indefinitely, it is not devoid of inconvenience and toxicity. Moreover, resistance or escape from disease control occur in a significant number of patients. Imatinib is a substrate of the cytochromes P450 CYP3A4/5 and of the multidrug transporter P glycoprotein (product of the MDR1 gene). Considering the large inter-individual differences in the expression and function of those systems, the disposition and clinical activity of imatinib can be expected to vary widely among patients, calling for dosage individualisation. The aim of this exploratory study was to determine the average pharmacokinetic parameters characterizing the disposition of imatinib in the target population, to assess their inter-individual variability, and to identify influential factors affecting them. A total of 321 plasma concentrations, taken at various sampling times after latest dose, were measured in 59 patients receiving Glivec® at diverse regimens, using a validated chromatographic method (HPLC-UV) developed for this study. The results were analysed by non-linear mixed effect modelling (NONMEM). A one- compartment model with first-order absorption appeared appropriate to describe the data, with an average apparent clearance of 12.4 l/h, a distribution volume of 268 l and an absorption constant of 0.47 h-1. The clearance was affected by body weight, age and sex. No influences of interacting drugs were found. DNA samples were used for pharmacogenetic explorations. The MDR1 polymorphism 3435C>T appears to affect the disposition of imatinib. Large inter-individual variability remained unexplained by the demographic covariates considered, both on clearance (40%) and distribution volume (71%). Together with intra-patient variability (34%), this translates into an 8-fold width of the 90%-prediction interval of plasma concentrations expected under a fixed dosing regimen ! This is a strong argument to further investigate the possible usefulness of a therapeutic drug monitoring programme for imatinib. It may help to individualise the dosing regimen before overt disease progression or observation of treatment toxicity, thus improving both the long-term therapeutic effectiveness and tolerability of this drug.

Tunable reporter signal production in feedback-uncoupled arsenic bioreporters.

Escherichia coli-based bioreporters for arsenic detection are typically based on the natural feedback loop that controls ars operon transcription. Feedback loops are known to show a wide range linear response to the detriment of the overall amplification of the incoming signal. While being a favourable feature in controlling arsenic detoxification for the cell, a feedback loop is not necessarily the most optimal for obtaining highest sensitivity and response in a designed cellular reporter for arsenic detection. Here we systematically explore the effects of uncoupling the topology of arsenic sensing circuitry on the developed reporter signal as a function of arsenite concentration input. A model was developed to describe relative ArsR and GFP levels in feedback and uncoupled circuitry, which was used to explore new ArsR-based synthetic circuits. The expression of arsR was then placed under the control of a series of constitutive promoters, which differed in promoter strength, and which could be further modulated by TetR repression. Expression of the reporter gene was maintained under the ArsR-controlled Pars promoter. ArsR expression in the systems was measured by using ArsR-mCherry fusion proteins. We find that stronger constitutive ArsR production decreases arsenite-dependent EGFP output from Pars and vice versa. This leads to a tunable series of arsenite-dependent EGFP outputs in a variety of systematically characterized circuitries. The higher expression levels and sensitivities of the response curves in the uncoupled circuits may be useful for improving field-test assays using arsenic bioreporters.

Wavelet analysis residual kriging vs. neural network residual kriging

This paper deals with the problem of spatial data mapping. A new method based on wavelet interpolation and geostatistical prediction (kriging) is proposed. The method - wavelet analysis residual kriging (WARK) - is developed in order to assess the problems rising for highly variable data in presence of spatial trends. In these cases stationary prediction models have very limited application. Wavelet analysis is used to model large-scale structures and kriging of the remaining residuals focuses on small-scale peculiarities. WARK is able to model spatial pattern which features multiscale structure. In the present work WARK is applied to the rainfall data and the results of validation are compared with the ones obtained from neural network residual kriging (NNRK). NNRK is also a residual-based method, which uses artificial neural network to model large-scale non-linear trends. The comparison of the results demonstrates the high quality performance of WARK in predicting hot spots, reproducing global statistical characteristics of the distribution and spatial correlation structure.

BME-based uncertainty assessment of the Chernobyl fallout

The vast territories that have been radioactively contaminated during the 1986 Chernobyl accident provide a substantial data set of radioactive monitoring data, which can be used for the verification and testing of the different spatial estimation (prediction) methods involved in risk assessment studies. Using the Chernobyl data set for such a purpose is motivated by its heterogeneous spatial structure (the data are characterized by large-scale correlations, short-scale variability, spotty features, etc.). The present work is concerned with the application of the Bayesian Maximum Entropy (BME) method to estimate the extent and the magnitude of the radioactive soil contamination by 137Cs due to the Chernobyl fallout. The powerful BME method allows rigorous incorporation of a wide variety of knowledge bases into the spatial estimation procedure leading to informative contamination maps. Exact measurements (?hard? data) are combined with secondary information on local uncertainties (treated as ?soft? data) to generate science-based uncertainty assessment of soil contamination estimates at unsampled locations. BME describes uncertainty in terms of the posterior probability distributions generated across space, whereas no assumption about the underlying distribution is made and non-linear estimators are automatically incorporated. Traditional estimation variances based on the assumption of an underlying Gaussian distribution (analogous, e.g., to the kriging variance) can be derived as a special case of the BME uncertainty analysis. The BME estimates obtained using hard and soft data are compared with the BME estimates obtained using only hard data. The comparison involves both the accuracy of the estimation maps using the exact data and the assessment of the associated uncertainty using repeated measurements. Furthermore, a comparison of the spatial estimation accuracy obtained by the two methods was carried out using a validation data set of hard data. Finally, a separate uncertainty analysis was conducted that evaluated the ability of the posterior probabilities to reproduce the distribution of the raw repeated measurements available in certain populated sites. The analysis provides an illustration of the improvement in mapping accuracy obtained by adding soft data to the existing hard data and, in general, demonstrates that the BME method performs well both in terms of estimation accuracy as well as in terms estimation error assessment, which are both useful features for the Chernobyl fallout study.

Statistical learning theory for geospatial data. Case study: Aral sea

In recent years there has been an explosive growth in the development of adaptive and data driven methods. One of the efficient and data-driven approaches is based on statistical learning theory (Vapnik 1998). The theory is based on Structural Risk Minimisation (SRM) principle and has a solid statistical background. When applying SRM we are trying not only to reduce training error ? to fit the available data with a model, but also to reduce the complexity of the model and to reduce generalisation error. Many nonlinear learning procedures recently developed in neural networks and statistics can be understood and interpreted in terms of the structural risk minimisation inductive principle. A recent methodology based on SRM is called Support Vector Machines (SVM). At present SLT is still under intensive development and SVM find new areas of application (www.kernel-machines.org). SVM develop robust and non linear data models with excellent generalisation abilities that is very important both for monitoring and forecasting. SVM are extremely good when input space is high dimensional and training data set i not big enough to develop corresponding nonlinear model. Moreover, SVM use only support vectors to derive decision boundaries. It opens a way to sampling optimization, estimation of noise in data, quantification of data redundancy etc. Presentation of SVM for spatially distributed data is given in (Kanevski and Maignan 2004).

The Smiling Project: development of a new program for gait and balance training

Introduction: The SMILING project, a multicentric project fundedby the European Union, aims to develop a new gait and balance trainingprogram to prevent falls in older persons. The program includes the"SMILING shoe", an innovative device that generates mechanical perturbationwhile walking by changing the soles' inclination. Induced perturbationschallenge subjects' balance and force them to react to avoidfalls. By training specifically the complex motor reactions used to maintainbalance when walking on irregular ground, the program will improvesubjects' ability to react in situation of unsteadiness and reduce theirrisk of falling. Methods: The program will be evaluated in a multicentric,cross-over randomized controlled trial. Overall, 112 subjects (aged≥65 years, ≥1 falls, POMA score 22-26/28) will be enrolled. Subjectswill be randomised in 2 groups: group A begin the training with active"SMILING shoes", group B with inactive dummy shoes. After 4 weeksof training, group A and B will exchange the shoes. Supervised trainingsessions (30 minutes twice a week for 8 weeks) include walkingtasks of progressive difficulties.To avoid a learning effect, "SMILINGshoes" perturbations will be generated in a non-linear and chaotic way.Gait performance, fear of falling, and acceptability of the program willbe assessed. Conclusion: The SMILING program is an innovative interventionfor falls prevention in older persons based on gait and balancetraining using chaotic perturbations. Because of the easy use of the"SMILING shoes", this program could be used in various settings, suchas geriatric clinics or at home.

Dynamic changes in brain functional connectivity during concurrent dual-task performance.

This study investigated the spatial, spectral, temporal and functional proprieties of functional brain connections involved in the concurrent execution of unrelated visual perception and working memory tasks. Electroencephalography data was analysed using a novel data-driven approach assessing source coherence at the whole-brain level. Three connections in the beta-band (18-24 Hz) and one in the gamma-band (30-40 Hz) were modulated by dual-task performance. Beta-coherence increased within two dorsofrontal-occipital connections in dual-task conditions compared to the single-task condition, with the highest coherence seen during low working memory load trials. In contrast, beta-coherence in a prefrontal-occipital functional connection and gamma-coherence in an inferior frontal-occipitoparietal connection was not affected by the addition of the second task and only showed elevated coherence under high working memory load. Analysis of coherence as a function of time suggested that the dorsofrontal-occipital beta-connections were relevant to working memory maintenance, while the prefrontal-occipital beta-connection and the inferior frontal-occipitoparietal gamma-connection were involved in top-down control of concurrent visual processing. The fact that increased coherence in the gamma-connection, from low to high working memory load, was negatively correlated with faster reaction time on the perception task supports this interpretation. Together, these results demonstrate that dual-task demands trigger non-linear changes in functional interactions between frontal-executive and occipitoparietal-perceptual cortices.

Breakpoints in ventilation, cerebral and muscle oxygenation, and muscle activity during an incremental cycling exercise.

The aim of this study was to locate the breakpoints of cerebral and muscle oxygenation and muscle electrical activity during a ramp exercise in reference to the first and second ventilatory thresholds. Twenty-five cyclists completed a maximal ramp test on an electromagnetically braked cycle-ergometer with a rate of increment of 25 W/min. Expired gazes (breath-by-breath), prefrontal cortex and vastus lateralis (VL) oxygenation [Near-infrared spectroscopy (NIRS)] together with electromyographic (EMG) Root Mean Square (RMS) activity for the VL, rectus femoris (RF), and biceps femoris (BF) muscles were continuously assessed. There was a non-linear increase in both cerebral deoxyhemoglobin (at 56 ± 13% of the exercise) and oxyhemoglobin (56 ± 8% of exercise) concomitantly to the first ventilatory threshold (57 ± 6% of exercise, p > 0.86, Cohen's d < 0.1). Cerebral deoxyhemoglobin further increased (87 ± 10% of exercise) while oxyhemoglobin reached a plateau/decreased (86 ± 8% of exercise) after the second ventilatory threshold (81 ± 6% of exercise, p < 0.05, d > 0.8). We identified one threshold only for muscle parameters with a non-linear decrease in muscle oxyhemoglobin (78 ± 9% of exercise), attenuation in muscle deoxyhemoglobin (80 ± 8% of exercise), and increase in EMG activity of VL (89 ± 5% of exercise), RF (82 ± 14% of exercise), and BF (85 ± 9% of exercise). The thresholds in BF and VL EMG activity occurred after the second ventilatory threshold (p < 0.05, d > 0.6). Our results suggest that the metabolic and ventilatory events characterizing this latter cardiopulmonary threshold may affect both cerebral and muscle oxygenation levels, and in turn, muscle recruitment responses.

Electrochemical As(III) whole-cell based biochip sensor.

The development of a whole-cell based sensor for arsenite detection coupling biological engineering and electrochemical techniques is presented. This strategy takes advantage of the natural Escherichia coli resistance mechanism against toxic arsenic species, such as arsenite, which consists of the selective intracellular recognition of arsenite and its pumping out from the cell. A whole-cell based biosensor can be produced by coupling the intracellular recognition of arsenite to the generation of an electrochemical signal. Hereto, E. coli was equipped with a genetic circuit in which synthesis of beta-galactosidase is under control of the arsenite-derepressable arsR-promoter. The E. coli reporter strain was filled in a microchip containing 16 independent electrochemical cells (i.e. two-electrode cell), which was then employed for analysis of tap and groundwater samples. The developed arsenic-sensitive electrochemical biochip is easy to use and outperforms state-of-the-art bacterial bioreporters assays specifically in its simplicity and response time, while keeping a very good limit of detection in tap water, i.e. 0.8ppb. Additionally, a very good linear response in the ranges of concentration tested (0.94ppb to 3.75ppb, R(2)=0.9975 and 3.75 ppb to 30ppb, R(2)=0.9991) was obtained, complying perfectly with the acceptable arsenic concentration limits defined by the World Health Organization for drinking water samples (i.e. 10ppb). Therefore, the proposed assay provides a very good alternative for the portable quantification of As (III) in water as corroborated by the analysis of natural groundwater samples from Swiss mountains, which showed a very good agreement with the results obtained by atomic absorption spectroscopy.

Accelerated Microstructure Imaging via Convex Optimization (AMICO) from diffusion MRI data.

Microstructure imaging from diffusion magnetic resonance (MR) data represents an invaluable tool to study non-invasively the morphology of tissues and to provide a biological insight into their microstructural organization. In recent years, a variety of biophysical models have been proposed to associate particular patterns observed in the measured signal with specific microstructural properties of the neuronal tissue, such as axon diameter and fiber density. Despite very appealing results showing that the estimated microstructure indices agree very well with histological examinations, existing techniques require computationally very expensive non-linear procedures to fit the models to the data which, in practice, demand the use of powerful computer clusters for large-scale applications. In this work, we present a general framework for Accelerated Microstructure Imaging via Convex Optimization (AMICO) and show how to re-formulate this class of techniques as convenient linear systems which, then, can be efficiently solved using very fast algorithms. We demonstrate this linearization of the fitting problem for two specific models, i.e. ActiveAx and NODDI, providing a very attractive alternative for parameter estimation in those techniques; however, the AMICO framework is general and flexible enough to work also for the wider space of microstructure imaging methods. Results demonstrate that AMICO represents an effective means to accelerate the fit of existing techniques drastically (up to four orders of magnitude faster) while preserving accuracy and precision in the estimated model parameters (correlation above 0.9). We believe that the availability of such ultrafast algorithms will help to accelerate the spread of microstructure imaging to larger cohorts of patients and to study a wider spectrum of neurological disorders.