In Vivo Mechanisms of Intestinal Drug Absorption from Aprepitant Nanoformulations.

Over recent decades there has been an increase in the proportion of BCS class II and IV drug candidates in industrial drug development. To overcome the biopharmaceutical challenges associated with the less favorable properties of solubility and/or intestinal permeation of these substances, the development of formulations containing nanosuspensions of the drugs has been suggested. The intestinal absorption of aprepitant from two nanosuspensions (20 µM and 200 µM total concentrations) in phosphate buffer, one nanosuspension (200 µM) in fasted-state simulated intestinal fluid (FaSSIF), and one solution (20 µM) in FaSSIF was investigated in the rat single-pass intestinal perfusion model. The disappearance flux from the lumen (Jdisapp) was faster for formulations containing a total concentration of aprepitant of 200 µM than for those containing 20 µM, but was unaffected by the presence of vesicles. The flux into the systemic circulation (Japp) and, subsequently, the effective diffusion constant (Deff) were calculated using the plasma concentrations. Japp was, like Jdisapp, faster for the formulations containing higher total concentrations of aprepitant, but was also faster for those containing vesicles (ratios of 2 and 1.5). This suggests that aprepitant is retained in the lumen when presented as nanoparticles in the absence of vesicles. In conclusion, increased numbers of nanoparticles and the presence of vesicles increased the rate of transport and availability of aprepitant in plasma. This effect can be attributed to an increased rate of mass transport through the aqueous boundary layer (ABL) adjacent to the gut wall.

Controlling the structure of spin-coated multilayer ethylcellulose/hydroxypropylcellulose films for drug release.

Porous phase-separated ethylcellulose/hydroxypropylcellulose (EC/HPC) films are used to control drug transport out of pharmaceutical pellets. Water-soluble HPC leaches out and forms a porous structure that controls the drug transport. Industrially, the pellets are coated using a fluidized bed spraying device, and a layered film exhibiting varying porosity and structure after leaching is obtained. A detailed understanding of the formation of the multilayered, phase-separated structure during production is lacking. Here, we have investigated multilayered EC/HPC films produced by sequential spin-coating, which was used to mimic the industrial process. The effects of EC/HPC ratio and spin speed on the multilayer film formation and structure were investigated using advanced microscopy techniques and image analysis. Cahn-Hilliard simulations were performed to analyze the mixing behavior. A gradient with larger structures close to the substrate surface and smaller structures close to the air surface was formed due to coarsening of the layers already coated during successive deposition cycles. The porosity of the multilayer film was found to vary with both EC/HPC ratio and spin speed. Simulation of the mixing behavior and in situ characterization of the structure evolution showed that the origin of the discontinuities and multilayer structure can be explained by the non-mixing of the layers.

Frequency-dependent signaling in cardiac myocytes.

Background: Recent experimental data support the view that signaling activity at the membrane depends on its geometric parameters such as surface area and curvature. However, a mathematical, biophysical concept linking shape to receptor signaling is missing. The membranes of cardiomyocytes are constantly reshaped due to cycles of contraction and relaxation. According to constant-volume behavior of cardiomyocyte contraction, the length shortening is compensated by Z-disc myofilament lattice expansion and dynamic deformation of membrane between two adjacent Z-discs. Both morphological changes are strongly dependent on the frequency of contraction. Here, we developed the hypothesis that dynamic geometry of cardiomyocytes could be important for their plasticity and signaling. This effect may depend on the frequency of the beating heart and may represent a novel concept to explain how changes in frequency affect cardiac signaling. Methods: This hypothesis is almost impossible to answer with experiments, as the in-vitro cardiomyocytes are almost two-dimensional and flattened rather than being in their real in-vivo shape. Therefore, we designed a COMSOL multiphysics program to mathematically model the dynamic geometry of a human cardiomyocyte and explore whether the beating frequency can modulate membrane signal transduction. Src kinase is an important component of cardiac mechanotransduction. We first presented that Src mainly localizes at costameres. Then, the frequency-dependent signaling effect was studied mathematically by numerical simulation of Src-mediated PDGFR signaling pathway. The reaction-convection-diffusion partial differential equation was formulated to simulate PDGFR pathway in a contracting sarcomeric disc for a range of frequencies from 1 to 4 Hz. Results: Simulations exhibits higher concentration of phospho-Src when a cardiomyocyte beats with higher rates. The calculated phospho-Src concentration at 4, 2, and 1 Hz beat rates, comparing to 0 Hz, was 21.5%, 9.4%, and 4.7% higher, respectively. Conclusion: Here we provide mathematical evidence for a novel concept in biology. Cell shape directly translates into signaling, an effect of importance particularly for the myocardium, where cells continuously reshape their membranes. The concept of locality of surface-to-volume ratios is demonstrated to lead to changes in membrane-mediated signaling and may help to explain the remarkable plasticity of the myocardium in response to biomechanical stress.

In silico prediction of drug solubility: 4. Will simple potentials suffice?

In view of the extreme importance of reliable computational prediction of aqueous drug solubility, we have established a Monte Carlo simulation procedure which appears, in principle, to yield reliable solubilities even for complex drug molecules. A theory based on judicious application of linear response and mean field approximations has been found to reproduce the computationally demanding free energy determinations by simulation while at the same time offering mechanistic insight. The focus here is on the suitability of the model of both drug and solvent, i.e., the force fields. The optimized potentials for liquid simulations all atom (OPLS-AA) force field, either intact or combined with partial charges determined either by semiempirical AM1/CM1A calculations or taken from the condensed-phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field has been used. The results illustrate the crucial role of the force field in determining drug solubilities. The errors in interaction energies obtained by the simple force fields tested here are still found to be too large for our purpose but if a component of this error is systematic and readily removed by empirical adjustment the results are significantly improved. In fact, consistent use of the OPLS-AA Lennard-Jones force field parameters with partial charges from the COMPASS force field will in this way produce good predictions of amorphous drug solubility within 1 day on a standard desktop PC. This is shown here by the results of extensive new simulations for a total of 47 drug molecules which were also improved by increasing the water box in the hydration simulations from 500 to 2000 water molecules.

In silico prediction of drug solubility: 1. Free energy of hydration.

As a first step in the computational prediction of drug solubility the free energy of hydration, DeltaG*(vw) in TIP4P water has been computed for a data set of 48 drug molecules using the free energy of perturbation method and the optimized potential for liquid simulations all-atom force field. The simulations were performed in two steps, where first the Coulomb and then the Lennard-Jones interactions between the solute and the water molecules were scaled down from full to zero strength to provide physical understanding and simpler predictive models. The results have been interpreted using a theory assuming DeltaG*(vw) = A(MS)gamma + E(LJ) + E(C)/2 where A(MS) is the molecular surface area, gamma is the water-vapor surface tension, and E(LJ) and E(C) are the solute-water Lennard-Jones and Coulomb interaction energies, respectively. It was found that by a proper definition of the molecular surface area our results as well as several results from the literature were found to be in quantitative agreement using the macroscopic surface tension of TIP4P water. This is in contrast to the surface tension for water around a spherical cavity that previously has been shown to be dependent on the size of the cavity up to a radius of approximately 1 nm. The step of scaling down the electrostatic interaction can be represented by linear response theory.

In silico prediction of drug solubility: 2. Free energy of solvation in pure melts.

The solubility of drugs in water is investigated in a series of papers and in the current work. The free energy of solvation, DeltaG*(vl), of a drug molecule in its pure drug melt at 673.15 K (400 degrees C) has been obtained for 46 drug molecules using the free energy perturbation method. The simulations were performed in two steps where first the Coulomb and then the Lennard-Jones interactions were scaled down from full to no interaction. The results have been interpreted using a theory assuming that DeltaG*(vl) = DeltaG(cav) + E(LJ) + E(C)/2 where the free energy of cavity formation, DeltaG(cav), in these pure drug systems was obtained using hard body theories, and E(LJ) and E(C) are the Lennard-Jones and Coulomb interaction energies, respectively, of one molecule with the other ones. Since the main parameter in hard body theories is the volume fraction, an equation of state approach was used to estimate the molecular volume. Promising results were obtained using a theory for hard oblates, in which the oblate axial ratio was calculated from the molecular surface area and volume obtained from simulations. The Coulomb term, E(C)/2, is half of the Coulomb energy in accord with linear response, which showed good agreement with our simulation results. In comparison with our previous results on free energy of hydration, the Coulomb interactions in pure drug systems are weaker, and the van der Waals interactions play a more important role.

In silico prediction of drug solubility. 3. Free energy of solvation in pure amorphous matter.

The solubility of drugs in water is investigated in a series of papers. In this work, we address the process of bringing a drug molecule from the vapor into a pure drug amorphous phase. This step enables us to actually calculate the solubility of amorphous drugs in water. In our general approach, we, on one hand, perform rigorous free energy simulations using a combination of the free energy perturbation and thermodynamic integration methods. On the other hand, we develop an approximate theory containing parameters that are easily accessible from conventional Monte Carlo simulations, thereby reducing the computation time significantly. In the theory for solvation, we assume that DeltaG* = DeltaGcav + ELJ + EC/2, where the free energy of cavity formation, DeltaGcav, in pure drug systems is obtained using a theory for hard-oblate spheroids, and ELJ and EC are the Lennard-Jones and Coulomb interaction energies between the chosen molecule and the others in the fluid. The theoretical predictions for the free energy of solvation in pure amorphous matter are in good agreement with free energy simulation data for 46 different drug molecules. These results together with our previous studies support our theoretical approach. By using our previous data for the free energy of hydration, we compute the total free energy change of bringing a molecule from the amorphous phase into water. We obtain good agreement between the theory and simulations. It should be noted that to obtain accurate results for the total process, high precision data are needed for the individual subprocesses. Finally, for eight different substances, we compare the experimental amorphous and crystalline solubility in water with the results obtained by the proposed theory with reasonable success.

In silico predictions of gastrointestinal drug absorption in pharmaceutical product development: application of the mechanistic absorption model GI-Sim.

Oral drug delivery is the predominant administration route for a major part of the pharmaceutical products used worldwide. Further understanding and improvement of gastrointestinal drug absorption predictions is currently a highly prioritized area of research within the pharmaceutical industry. The fraction absorbed (fabs) of an oral dose after administration of a solid dosage form is a key parameter in the estimation of the in vivo performance of an orally administrated drug formulation. This study discloses an evaluation of the predictive performance of the mechanistic physiologically based absorption model GI-Sim. GI-Sim deploys a compartmental gastrointestinal absorption and transit model as well as algorithms describing permeability, dissolution rate, salt effects, partitioning into micelles, particle and micelle drifting in the aqueous boundary layer, particle growth and amorphous or crystalline precipitation. Twelve APIs with reported or expected absorption limitations in humans, due to permeability, dissolution and/or solubility, were investigated. Predictions of the intestinal absorption for different doses and formulations were performed based on physicochemical and biopharmaceutical properties, such as solubility in buffer and simulated intestinal fluid, molecular weight, pK(a), diffusivity and molecule density, measured or estimated human effective permeability and particle size distribution. The performance of GI-Sim was evaluated by comparing predicted plasma concentration-time profiles along with oral pharmacokinetic parameters originating from clinical studies in healthy individuals. The capability of GI-Sim to correctly predict impact of dose and particle size as well as the in vivo performance of nanoformulations was also investigated. The overall predictive performance of GI-Sim was good as >95% of the predicted pharmacokinetic parameters (C(max) and AUC) were within a 2-fold deviation from the clinical observations and the predicted plasma AUC was within one standard deviation of the observed mean plasma AUC in 74% of the simulations. GI-Sim was also able to correctly capture the trends in dose- and particle size dependent absorption for the study drugs with solubility and dissolution limited absorption, respectively. In addition, GI-Sim was also shown to be able to predict the increase in absorption and plasma exposure achieved with nanoformulations. Based on the results, the performance of GI-Sim was shown to be suitable for early risk assessment as well as to guide decision making in pharmaceutical formulation development.

Nucleation and crystal growth in supersaturated solutions of a model drug.

The crystallization process in aqueous solutions of the drug bicalutamide and the effect of the polymer polyvinylpyrrolidone (PVP) have been studied. Results showed that PVP decreased the crystallization rate significantly in a system with PVP concentrations as low as 0.01% (w/w), without changing the polymorph formed. The crystal habit was altered already at PVP concentrations as low as 0.001% (w/w). Measurements made with self-diffusion NMR indicated that the decrease in crystallization rate was not because of a reduced supersaturation due to bicalutamide binding to PVP in solution. Furthermore, in experiments designed to specifically study crystal nucleation, the same nucleation rate was found in the absence and presence of PVP. Instead, PVP adsorbs to the crystals formed in solution and by doing so, the crystal growth rate is reduced. This was confirmed in separate experiments using bicalutamide nanocrystals. By combining theories describing classical nucleation and crystal growth, with some modifications, a consistent description of several independent experiments performed in polymer-free systems was obtained. From these experiments a crystal-water interfacial tension of 22.1 mN/m was extracted. We also analyze the interfacial tension of other crystalline organic solids and find that it varies approximately as the logarithm of the solubility. This finding is discussed within the framework of the Bragg-Williams regular solution theory where we also compare with the tension of liquid alkanes.

Amorphous drug nanosuspensions. 3. Particle dissolution and crystal growth.

In the present paper, we have studied particle dissolution and crystal growth of the poorly water soluble drug felodipine, using fluorescence as a probe for the amount of crystalline material. Dissolution kinetics is essentially diffusion-controlled, while the rate of crystal growth is significantly slower compared to the diffusion-controlled limit. The deviation from diffusion control was characterized by the effective length, lambda, related to the kinetics of a surface integration process. Amorphous nanoparticles may be highly unstable in the presence of small amounts of crystalline particles. This is due to the fact that the molecular solubility from the amorphous nanoparticles often is at least an order of magnitude higher than the corresponding crystalline solubility. In a mixed system where crystalline nanoparticles have been added to an amorphous nanosuspension, the bulk will have a monomer concentration intermediate between the amorphous and crystalline solubilities, and is thus supersaturated with respect to the crystalline particles while being undersaturated with respect to the amorphous particles. As a consequence, the amorphous particles spontaneously dissolve, while crystalline particles grow, in a combined process which is similar to Ostwald ripening. By knowing the parameters describing dissolution and crystal growth, respectively, it was possible to simulate the outcome of controlled seeding experiments, where a small amount of crystalline nanoparticles was added to a dispersion of amorphous nanoparticles. A good agreement between model calculations and experiments was obtained including how the crystal growth rate varied with the amounts of added crystalline seeds.