Quantitative imaging of doxorubicin diffusion and cellular uptake in biomimetic gels with human liver tumor cells.

Novel tumor-on-a-chip approaches are increasingly used to investigate tumor progression and potential treatment options. To improve the effect of any cancer treatment it is important to have an in depth understanding of drug diffusion, penetration through the tumor extracellular matrix and cellular uptake. In this study, we have developed a miniaturized chip where drug diffusion and cellular uptake in different hydrogel environments can be quantified at high resolution using live imaging. Diffusion of doxorubicin was reduced in a biomimetic hydrogel mimicking tissue properties of cirrhotic liver and early stage hepatocellular carcinoma (373 ± 108 µm2/s) as compared to an agarose gel (501 ± 77 µm2/s, p = 0.019). The diffusion was further lowered to 256 ± 30 µm2/s (p = 0.028) by preparing the biomimetic gel in cell media instead of phosphate buffered saline. The addition of liver tumor cells (Huh7 or HepG2) to the gel, at two different densities, did not significantly influence drug diffusion. Clinically relevant and quantifiable doxorubicin concentration gradients (1-20 µM) were established in the chip within one hour. Intracellular increases in doxorubicin fluorescence correlated with decreasing fluorescence of the DNA-binding stain Hoechst 33342 and based on the quantified intracellular uptake of doxorubicin an apparent cell permeability (9.00 ± 0.74 × 10-4 µm/s for HepG2) was determined. Finally, the data derived from the in vitro model were applied to a spatio-temporal tissue concentration model to evaluate the potential clinical impact of a cirrhotic extracellular matrix on doxorubicin diffusion and tumor cell uptake.

Effect of fluid velocity and particle size on the hydrodynamic diffusion layer thickness.

The aim of this study was to determine the thickness of the hydrodynamic diffusion layer (hHDL) of three poor water-soluble compounds under laminar fluid flow using a single particle dissolution technique. The single particle dissolution experiments were performed in a flowing aqueous medium using four different fluid velocities (v), ranging from 46 to 103 mm/s. The particles used had an initial radius (r) of 18.8 to 52.3 µm. The determined hHDL values were calculated from both dissolution experiments and computational fluid dynamics (CFD) simulation. In this study, single particle dissolution experiments gave, with one exception, hHDL values in the range of 2.09 to 8.85 µm and corresponding simulations gave hHDL values in the range of 2.53 to 4.38 µm. Hence, we found a semi-quantitative concordance between experimental and simulated determined hHDL values. Also, a theoretical relation between the dependence of hHDL on particle radius and flow velocity of the medium was established by a series of CFD simulations in a fluid velocity range of 10-100 mm/s and particle size (radius) range of 5-40 µm. The outcome suggests a power law relation of the form hHDL∝r3/5v-2/5. In addition, the hHDL seems to be independent of the solubility, while it has a diffusion coefficient dependence. In conclusion, the hHDL values were determined under well-defined conditions; hence, this approach can be used to estimate the hHDL under different conditions to increase the understanding of the mass transfer mechanisms during the dissolution process.

Drug diffusion in biomimetic hydrogels: importance for drug transport and delivery in non-vascular tumor tissue.

Hydrogels of varying complexity are routinely used as scaffolds and 3D structures for in vitro tumor models to increase physiological relevance within pre-clinical cancer research. Relatively simple hydrogels such as agarose are well characterised, meanwhile biomimetic gels containing collagen and fibrin(ogen) have been studied to a much lesser extent. In this study, hydrogels mimicking the biophysical characteristics of liver cancer progression were investigated in terms of their UV-properties and influence on diffusion coefficients of different substances. UV-imaging technology was used to both visualize and quantify the diffusion process in a simple and rapid way. In general, agarose gel diffusion agreed well with predictions using the Stokes-Einstein equation meanwhile the biomimetic gels reduced diffusion coefficients by up to 70%. For doxorubicin, spatio-temporal tissue concentration modelling was used to translate in vitro diffusion to the more clinical context of tumor penetration in a solid liver tumor supplied by arterial blood.

In Vitro Cell Toxicity and Intracellular Uptake of Doxorubicin Exposed as a Solution or Liposomes: Implications for Treatment of Hepatocellular Carcinoma.

Cytostatic effects of doxorubicin in clinically applied doses are often inadequate and limited by systemic toxicity. The main objective of this in vitro study was to determine the anti-tumoral effect (IC50) and intracellular accumulation of free and liposomal doxorubicin (DOX) in four human cancer cell lines (HepG2, Huh7, SNU449 and MCF7). The results of this study showed a correlation between longer DOX exposure time and lower IC50 values, which can be attributed to an increased cellular uptake and intracellular exposure of DOX, ultimately leading to cell death. We found that the total intracellular concentrations of DOX were a median value of 230 times higher than the exposure concentrations after exposure to free DOX. The intracellular uptake of DOX from solution was at least 10 times higher than from liposomal formulation. A physiologically based pharmacokinetic model was developed to translate these novel quantitative findings to a clinical context and to simulate clinically relevant drug concentration-time curves. This showed that a liver tumor resembling the liver cancer cell line SNU449, the most resistant cell line in this study, would not reach therapeutic exposure at a standard clinical parenteral dose of doxorubicin (50 mg/m2), which is serious limitation for this drug. This study emphasizes the importance of in-vitro to in-vivo translations in the assessment of clinical consequence of experimental findings.

Hydrogen bond effects in multimode nuclear dynamics of acetic acid observed via resonant x-ray scattering.

A theoretical and experimental study of the gas phase and liquid acetic acid based on resonant inelastic x-ray scattering (RIXS) spectroscopy is presented. We combine and compare different levels of theory for an isolated molecule for a comprehensive analysis, including electronic and vibrational degrees of freedom. The excitation energy scan over the oxygen K-edge absorption reveals nuclear dynamic effects in the core-excited and final electronic states. The theoretical simulations for the monomer and two different forms of the dimer are compared against high-resolution experimental data for pure liquid acetic acid. We show that the theoretical model based on a dimer describes the hydrogen bond formation in the liquid phase well and that this bond formation sufficiently alters the RIXS spectra, allowing us to trace these effects directly from the experiment. Multimode vibrational dynamics is accounted for in our simulations by using a hybrid time-dependent stationary approach for the quantum nuclear wave packet simulations, showing the important role it plays in RIXS.

Vibrational resonant inelastic X-ray scattering in liquid acetic acid: a ruler for molecular chain lengths.

Quenching of vibrational excitations in resonant inelastic X-ray scattering (RIXS) spectra of liquid acetic acid is observed. At the oxygen core resonance associated with localized excitations at the O-H bond, the spectra lack the typical progression of vibrational excitations observed in RIXS spectra of comparable systems. We interpret this phenomenon as due to strong rehybridization of the unoccupied molecular orbitals as a result of hydrogen bonding, which however cannot be observed in x-ray absorption but only by means of RIXS. This allows us to address the molecular structure of the liquid, and to determine a lower limit for the average molecular chain length.

Long-Term Physical (In)Stability of Spray-Dried Amorphous Drugs: Relationship with Glass-Forming Ability and Physicochemical Properties.

This study shows the importance of the chosen method for assessing the glass-forming ability (GFA) and glass stability (GS) of a drug compound. Traditionally, GFA and GS are established using in situ melt-quenching in a differential scanning calorimeter. In this study, we included 26 structurally diverse glass-forming drugs (i) to compare the GFA class when the model drugs were produced by spray-drying with that when melt-quenching was used, (ii) to investigate the long-term physical stability of the resulting amorphous solids, and (iii) to investigate the relationship between physicochemical properties and the GFA of spray-dried solids and their long-term physical stability. The spray-dried solids were exposed to dry (<5% RH) and humid (75% RH) conditions for six months at 25 °C. The crystallization of the spray-dried solids under these conditions was monitored using a combination of solid-state characterization techniques including differential scanning calorimetry, Raman spectroscopy, and powder X-ray diffraction. The GFA/GS class assignment for 85% of the model compounds was method-dependent, with significant differences between spray-drying and melt-quenching methods. The long-term physical stability under dry condition of the compounds was predictable from GFA/GS classification and glass transition and crystallization temperatures. However, the stability upon storage at 75% RH could not be predicted from the same data. There was no strong correlation between the physicochemical properties explored and the GFA class or long-term physical stability. However, there was a slight tendency for compounds with a relatively larger molecular weight, higher glass transition temperature, higher crystallization temperature, higher melting point and higher reduced glass transition temperature to have better GFA and better physical stability. In contrast, a high heat of fusion and entropy of fusion seemed to have a negative impact on the GFA and physical stability of our dataset.

Supersaturation Potential of Amorphous Active Pharmaceutical Ingredients after Long-Term Storage.

This study explores the effect of physical aging and/or crystallization on the supersaturation potential and crystallization kinetics of amorphous active pharmaceutical ingredients (APIs). Spray-dried, fully amorphous indapamide, metolazone, glibenclamide, hydrocortisone, hydrochlorothiazide, ketoconazole, and sulfathiazole were used as model APIs. The parameters used to assess the supersaturation potential and crystallization kinetics were the maximum supersaturation concentration (Cmax,app), the area under the curve (AUC), and the crystallization rate constant (k). These were compared for freshly spray-dried and aged/crystallized samples. Aged samples were stored at 75% relative humidity for 168 days (6 months) or until they were completely crystallized, whichever came first. The solid-state changes were monitored with differential scanning calorimetry, Raman spectroscopy, and powder X-ray diffraction. Supersaturation potential and crystallization kinetics were investigated using a tenfold supersaturation ratio compared to the thermodynamic solubility using the µDISS Profiler. The physically aged indapamide and metolazone and the minimally crystallized glibenclamide and hydrocortisone did not show significant differences in their Cmax,app and AUC when compared to the freshly spray-dried samples. Ketoconazole, with a crystalline content of 23%, reduced its Cmax,app and AUC by 50%, with Cmax,app being the same as the crystalline solubility. The AUC of aged metolazone, one of the two compounds that remained completely amorphous after storage, significantly improved as the crystallization kinetics significantly decreased. Glibenclamide improved the most in its supersaturation potential from amorphization. The study also revealed that, besides solid-state crystallization during storage, crystallization during dissolution and its corresponding pathway may significantly compromise the supersaturation potential of fully amorphous APIs.

Comminution-amorphisation relationships during ball milling of lactose at different milling conditions.

The purpose of the study was to investigate the relationship between comminution and amorphisation of α-lactose monohydrate particles during ball milling under different milling conditions, including ball-to-powder mass ratio, milling time and ball diameter. The results revealed that at a constant ball filling ratio, ball-to-powder mass ratio of 25:1 resulted in the lowest minimum particle diameter of ∼5µm and the highest degree of apparent amorphous content of 82%. The rate of comminution was high during early stage of milling whereas the degree of apparent amorphous content increased gradually at a slow rate. An increased ball-to-powder mass ratio during milling increased both the rate of comminution and the rate of amorphisation. Using a given ball-to-powder mass ratio, the ball diameter affected the degree of apparent amorphous content of the particles while the particle diameter remained unchanged. The relationship between comminution and amorphisation could be described as consisting of two stages, i.e. comminution dominated and amorphisation dominated stage. It was proposed that the rate constant of comminution and amorphisation are controlled by stress energy distribution in the milling jar and the stress energy distribution is regulated by the ball motion pattern that can be affected by the process parameter used.

Considerations on the quantitative analysis of apparent amorphicity of milled lactose by Raman spectroscopy.

The main purpose of the study was to evaluate various pre-processing and quantification approaches of Raman spectrum to quantify low level of amorphous content in milled lactose powder. To improve the quantification analysis, several spectral pre-processing methods were used to adjust background effects. The effects of spectral noise on the variation of determined amorphous content were also investigated theoretically by propagation of error analysis and were compared to the experimentally obtained values. Additionally, the applicability of calibration method with crystalline or amorphous domains in the estimation of amorphous content in milled lactose powder was discussed. Two straight baseline pre-processing methods gave the best and almost equal performance. By the succeeding quantification methods, PCA performed best, although the classical least square analysis (CLS) gave comparable results, while peak parameter analysis displayed to be inferior. The standard deviations of experimental determined percentage amorphous content were 0.94% and 0.25% for pure crystalline and pure amorphous samples respectively, which was very close to the standard deviation values from propagated spectral noise. The reasonable conformity between the milled samples spectra and synthesized spectra indicated representativeness of physical mixtures with crystalline or amorphous domains in the estimation of apparent amorphous content in milled lactose.

Local electronic structure of functional groups in glycine as anion, zwitterion, and cation in aqueous solution.

Nitrogen and oxygen K emission spectra of glycine in the form of anions, zwitterions, and cations in aqueous solution are presented. It is shown that protonation has a dramatic influence on the local electronic structure and that the functional groups give a distinct spectral fingerprint.

Electronic structure of water molecules confined in a micelle lattice.

Oxygen K absorption and emission spectra of water molecules confined in dodecyltrimethyl ammonium chloride micelle structures are presented. The local electronic structure of the water molecules is found to be dramatically different from the electronic structure of water molecules in the gas-phase as well as in liquid water. Hybridization with states of the ions in the surrounding ions is directly observed, and evidence for stabilization of the water molecules relative to molecules in bulk water is found.

Implications of regular solution theory on the release mechanism of catanionic mixtures from gels.

The aim of this study was to apply the regular solution theory of mixed micelles to gain new insights on the drug release mechanism, when using catanionic mixtures as a method of obtaining prolonged release from gels. Synergistic effects were investigated at equilibrium and quantified in terms of regular solution theory interaction parameters. The drug release from catanionic aggregates was studied both in a polymer free environment, using dialysis membranes, and in gels, using a modified USP paddle method. The drug release kinetics was modelled theoretically by combining the regular solution theory with Fick's diffusion laws assuming a contribution to the transport only from monomeric species (stationary aggregates). The theoretical predictions were found to be in reasonably good agreement with experiments. An analysis of the calculated distribution of species between aggregated and monomeric states was shown to provide further insights into the release mechanism.

The multiple depletion curves method provides accurate estimates of intrinsic clearance (CLint), maximum velocity of the metabolic reaction (Vmax), and Michaelis constant (Km): accuracy and robustness evaluated through experimental data and Monte Carlo simulations.

The use of multiple depletion curves for the estimation of maximum velocity of the metabolic reaction (V(max)), the Michaelis constant (K(m)), and intrinsic clearance (CL(int)) was thoroughly evaluated by means of experimental data and through a series of Monte Carlo simulations. The enzyme kinetics of seven compounds were determined using the multiple depletion curves method (MDCM), the traditional initial formation rate of metabolite method (IFRMM), and the "in vitro t(1/2)" method, and the parameter estimates that were derived from the three methods were compared. The impact of a change in enzyme activity during the incubation period on the parameter estimates and the possibility to correct for this were also investigated. The MDCM was in good overall agreement with the IFRMM. Correction for a change in enzyme activity was possible and resulted in a better concordance in CL(int) estimates. The robustness of the method in coping with different rates of substrate turnover and variable starting concentrations were also demonstrated through Monte Carlo simulations. Furthermore, the limitations imposed by assumptions inherent in the in vitro t(1/2) method were demonstrated both experimentally and by simulations. This study demonstrates that the MDCM is a robust and efficient method for estimating enzyme kinetic variables with high accuracy and precision. The method may potentially be used in a wide range of applications, from pure enzyme kinetics to in vitro-based predictions of the pharmacokinetics of compounds with multiple and/or unknown metabolic pathways.

Photoinduced formation of N2 molecules in ammonium compounds.

Via fluorescence yield (FY) and resonant inelastic scattering spectroscopy in the soft X-ray range we find that soft X-rays induce formation of N2 molecules in solid NH4Cl and in related compounds. The nitrogen molecules form weak bonds in NH4Cl, so that a substantial fraction of the molecules remains in the sample. From measurements of the FY as a function of exposure and temperature, the rates for the photochemical processes are estimated. At elevated temperatures (363 K), several nitrogen atoms are removed from the sample per incoming photon. At lower temperatures (233 K), the rate is reduced to around 0.02 nitrogen atoms for each incoming photon. Virtually all these atoms form N2 molecules which are bound in the sample. The generality and implications of these results are briefly discussed.

Conductivity percolation in loosely compacted microcrystalline cellulose: An in situ study by dielectric spectroscopy during densification.

The present study aims at contributing to a complete understanding of the water-induced ionic charge transport in cellulose. The behavior of this transport in loosely compacted microcrystalline cellulose (MCC) powder was investigated as a function of density utilizing a new type of measurement setup, allowing for dielectric spectroscopy measurement in situ during compaction. The ionic conductivity in MCC was found to increase with increasing density until a leveling-out was observed for densities above approximately 0.7 g/cm3. Further, it was shown that the ionic conductivity vs density followed a percolation type behavior signifying the percolation of conductive paths in a 3D conducting network. The density percolation threshold was found to be between approximately 0.2 and 0.4 g/cm3, depending strongly on the cellulose moisture content. The observed percolation behavior was attributed to the forming of interparticulate bonds in the MCC and the percolation threshold dependence on moisture was linked to the moisture dependence of particle rearrangement and plastic deformation in MCC during compaction. The obtained results add to the understanding of the density-dependent water-induced ionic transport in cellulose showing that, at given moisture content, the two major parameters determining the magnitude of the conductivity are the connectedness of the interparticluate bonds and the connectedness of pores with a diameter in the 5-20 nm size range. At densities between approximately 0.7 and 1.2 g/cm3 both the bond and the pore networks have percolated, facilitating charge transport through the MCC compact.

Regional levels of drug transporters along the human intestinal tract: co-expression of ABC and SLC transporters and comparison with Caco-2 cells.

A vast number of drugs are subjected to active or facilitated transport and multiple transport mechanism may contribute to the net flux during drug absorption. The main objective of this study was to quantify the regional mRNA expression and determine the co-expression of drug transporters from the ABC (Pgp, BCRP, MRP2, MRP3) and SLC (PEPT1, MCT1, OATPB, OCTN2, OCT1) families along the human intestine (duodenum, jejunum, ileum, and colon). A second objective was to compare the transporter expression between the different intestinal regions and Caco-2 cells. Eight out of nine of the investigated transporters exhibited significant regional differences in expression. OATPB was the only transporter that did not show a region-dependency in the expression along the human intestinal canal. The expression of Pgp, BCRP, OCTN2 and MCT1 differed along the small intestine, but the expression differences were greater than five-fold only for Pgp. The rank order of transcript prevalence was identical in the ileum and the jejunum. Between the ileum and colon, seven transcripts were differentially expressed, and MCT1, OCTN2 and MRP3 were expressed at higher levels in the colon than in the small intestine. The expression of transporters in Caco-2 was closest to the expression pattern in the small intestine, although the expression of OATPB, BCRP and MRP2 differed more than five-fold between the Caco-2 cells and ileum. In conclusion, this study provides quantitative data on the expression of transporters from the ABC and SLC families along the human intestine, which can be useful in the interpretation of clinical studies where more than one intestinal transporter contribute to the net transport and in the computer modelling of drug absorption.

Direct nose-to-brain transfer of morphine after nasal administration to rats.

PURPOSE: The aim of this study was to quantify the olfactory transfer of morphine to the brain hemispheres by comparing brain tissue and plasma morphine levels after nasal administration with those after intravenous administration. METHODS: Morphine (1.0 mg/kg body weight) was administered via the right nostril or intravenously as a 15-min constant-rate infusion to male rats. The content of morphine and its metabolite morphine-3-glucuronide in samples of the olfactory bulbs, brain hemispheres, and plasma was assessed using high-performance liquid chromatography, and the areas under the concentration-time curves (AUC) were calculated. RESULTS: At both 5 and 15 min after administration, brain hemisphere morphine concentrations after nasal administration were similar to those after i.v. administration of the same dose, despite lower plasma concentrations after nasal administration. The brain hemispheres/plasma morphine AUC ratios for the 0-5 min period were thus approximately 3 and 0.1 after nasal and i.v. administration, respectively, demonstrating a statistically significant early distribution advantage of morphine to the brain hemispheres via the nasal route. CONCLUSION: Morphine is transferred via olfactory pathways to the brain hemispheres, and drug transfer via this route significantly contributes to the early high brain concentrations after nasal administration to rats.

Caco-2 permeability of weakly basic drugs predicted with the double-sink PAMPA pKa(flux) method.

The aim of this study was to analyze pH-dependent permeability of cationic drugs in Caco-2 cell monolayers using the pK(a)(flux) method and to correlate the results with those obtained in PAMPA (parallel artificial membrane permeability assay). The pH-dependent permeability of verapamil and propranolol was studied in Caco-2 cell monolayers. The data were subsequently processed using software developed for the PAMPA pK(a)(flux) method. Literature values for an additional nine cationic drugs were also analyzed. Double-Sink PAMPA data were also obtained for the same cationic drugs, to compare with the Caco-2 data. The Algorithm Builder program was then used to develop a predictive model of Caco-2 permeability based on PAMPA permeability and calculated Abraham molecular descriptors. From the relationship between permeability and pH it was shown that in PAMPA only the uncharged form of the drugs permeated across the membrane barrier, while charged and ionized forms of the drugs were significantly permeable in Caco-2. The charged-form permeability, P(i), was therefore determined and subsequently subtracted from all permeability coefficients in Caco-2 prior to the comparison with PAMPA. The resulting intrinsic permeability coefficients (P(o)) obtained in Caco-2 were successfully related to those derived from the PAMPA model. In this study we have shown that permeability coefficients obtained in PAMPA can predict the passive transcellular permeability in Caco-2.

Formulating insulin for oral administration: preparation of hyaluronan-insulin complex.

PURPOSE: To investigate the behaviour of peptides and hyaluronan in strong acid solutions containing electrolytes in the preparation of a new formulation of insulin, hyaluronan-insulin complex, and to evaluate the in vivo oral activity of the formulation. METHODS: Individual processing parameters in the preparation of the insulin complex were first refined, and two formulations were subsequently investigated. The chemical structure, particle size and hydrophilic/hydrophobic properties of the insulin complex in these formulations were studied using light scattering techniques, amino acid analysis, atomic force microscopy and cryo-transmission electron microscopy. The in vivo activity of oral hyaluronan-insulin complex was then evaluated by measuring the decrease in blood glucose concentrations in streptozotocin diabetic rats. RESULTS: Five of seven batches of the two insulin complex formulations fit the baseline criteria for approval of the new formulation. The formulation consists of a transparent aqua sol containing a solid hydrophobic phase as precipitate. Glucose-lowering activity was demonstrated after oral administration of the insulin complex'to diabetic rats. CONCLUSION: A new insulin formulation, a hyaluronan-insulin complex, has been developed and oral activity has been demonstrated.