Analysis of stabilization mechanisms in β-lactoglobulin-based amorphous solid dispersions by experimental and computational approaches

Our previous work shows that β-lactoglobulin-stabilized amorphous solid dispersion (ASD) loaded with 70 % indomethacin remains stable for more than 12 months. The stability is probably due to hydrogen bond networks spread throughout the ASD, facilitated by the indomethacin which has both hydrogen donors and acceptors. To investigate the stabilization mechanisms further, here we tested five other drug molecules, including two without any hydrogen bond donors. A combination of experimental techniques (differential scanning calorimetry, X-ray power diffraction) and molecular dynamics simulations was used to find the maximum drug loadings for ASDs with furosemide, griseofulvin, ibuprofen, ketoconazole and rifaximin. This approach revealed the underlying stabilization factors and the capacity of computer simulations to predict ASD stability. We searched the ASD models for crystalline patterns, and analyzed diffusivity of the drug molecules and hydrogen bond formation. ASDs loaded with rifaximin and ketoconazole remained stable for at least 12 months, even at 90 % drug loading, whereas stable drug loadings for furosemide, griseofulvin and ibuprofen were at a maximum of 70, 50 and 40 %, respectively. Steric confinement and hydrogen bonding to the proteins were the most important stabilization mechanisms at low drug loadings (≤ 40 %). Inter-drug hydrogen bond networks (including those with induced donors), ionic interactions, and a high Tg of the drug molecule were additional factors stabilizing the ASDs at drug loading greater than 40 %.

Molecular dynamics study on micelle-small molecule interactions: developing a strategy for an extensive comparison

Theoretical predictions of the solubilizing capacity of micelles and vesicles present in intestinal fluid are important for the development of new delivery techniques and bioavailability improvement. A balance between accuracy and computational cost is a key factor for an extensive study of numerous compounds in diverse environments. In this study, we aimed to determine an optimal molecular dynamics (MD) protocol to evaluate small-molecule interactions with micelles composed of bile salts and phospholipids. MD simulations were used to produce free energy profiles for three drug molecules (danazol, probucol, and prednisolone) and one surfactant molecule (sodium caprate) as a function of the distance from the colloid center of mass. To address the challenges associated with such tasks, we compared different simulation setups, including freely assembled colloids versus pre-organized spherical micelles, full free energy profiles versus only a few points of interest, and a coarse-grained model versus an all-atom model. Our findings demonstrate that combining these techniques is advantageous for achieving optimal performance and accuracy when evaluating the solubilization capacity of micelles. All-atom (AA) and coarse-grained (CG) umbrella sampling (US) simulations and point-wise free energy (FE) calculations were compared to their efficiency to computationally analyze the solubilization of active pharmaceutical ingredients in intestinal fluid colloids.

Mechanisms of Drug Solubility Enhancement Induced by β-Lactoglobulin-Based Amorphous Solid Dispersions

Protein-based amorphous solid dispersions (ASDs) have emerged as a promising approach for enhancing solubility in comparison to crystalline drugs. The dissolution behavior of protein-based amorphous solid dispersions (ASDs) was investigated in various pH media. ASDs of four poorly soluble model drugs with acidic (furosemide and indomethacin), basic (carvedilol), and neutral (celecoxib) properties were prepared by spray drying at 30 wt % drug loading with the protein β-lactoglobulin (BLG). The effect of spray-dried BLG (SD-BLG) solubility and protein binding ability with dissolved drugs in solution were investigated to retrieve the mechanisms governing the improvement of drug solubility from the BLG-based ASDs. Powder dissolution results showed that all ASDs obtained a higher maximum concentration (Cmax) compared to the respective pure crystalline drugs. It was found that the solubility increase of the drugs from the ASDs was to a large extent dependent on the solubility of the pure SD-BLG at the investigated pH values (low solubility at pH near the isoelectric point (pI) of BLG). Furthermore, drug-protein interactions in a solution were observed, in particular at pH values where the drugs were neutral. These drug-protein interactions also resulted, to some extent, in the stabilization of the drug in supersaturation.

Stabilizing Mechanisms of β-Lactoglobulin in Amorphous Solid Dispersions of Indomethacin

Proteins, and in particular whey proteins, have recently been introduced as a promising excipient class for stabilizing amorphous solid dispersions. However, despite the efficacy of the approach, the molecular mechanisms behind the stabilization of the drug in the amorphous form are not yet understood. To investigate these, we used experimental and computational techniques to study the impact of drug loading on the stability of protein-stabilized amorphous formulations. β-Lactoglobulin, a major component of whey, was chosen as a model protein and indomethacin as a model drug. Samples, prepared by either ball milling or spray drying, formed single-phase amorphous solid dispersions with one glass transition temperature at drug loadings lower than 40-50%; however, a second glass transition temperature appeared at drug loadings higher than 40-50%. Using molecular dynamics simulations, we found that a drug-rich phase occurred at a loading of 40-50% and higher, in agreement with the experimental data. The simulations revealed that the mechanisms of the indomethacin stabilization by β-lactoglobulin were a combination of (a) reduced mobility of the drug molecules in the first drug shell and (b) hydrogen-bond networks. These networks, formed mostly by glutamic and aspartic acids, are situated at the β-lactoglobulin surface, and dependent on the drug loading (>40%), propagated into the second and subsequent drug layers. The simulations indicate that the reduced mobility dominates at low (<40%) drug loadings, whereas hydrogen-bond networks dominate at loadings up to 75%. The computer simulation results agreed with the experimental physical stability data, which showed a significant stabilization effect up to a drug fraction of 70% under dry storage. However, under humid conditions, stabilization was only sufficient for drug loadings up to 50%, confirming the detrimental effect of humidity on the stability of protein-stabilized amorphous formulations.

Endothelial glycocalyx permeability for nanoscale solutes

Glycocalyx has a great impact on the accessibility of the endothelial cell membranes. Although the specific interactions play a crucial role in cross-membrane solute transport, nonspecific interactions cannot be neglected. In this work, we used computational modeling to quantify the nonspecific interactions that control the distribution of nanosized solutes across the endothelial glycocalyx. We evaluated the probabilities of various nanoparticles' passage through the luminal layer to the membrane. The calculations demonstrate that excluded volume and electrostatic interactions are decisive for the solute transport as compared with van der Waals and hydrodynamic interactions. Damaged glycocalyx models showed a relatively weak efficiency in sieving plasma solutes. We estimated the energy barriers and corresponding mean first passage times for nanoscale solute transport through the model glycocalyx.

Pharmaceutical profiling and molecular dynamics simulations reveal crystallization effects in amorphous formulations

Robust and reliable in vivo performance of medicines based on amorphous solid dispersions (ASDs) depend on maintenance of physical stability and efficient supersaturation. However, molecular drivers of these two kinetic processes are poorly understood. Here we used molecular dynamics (MD) simulations coupled with experimental assessments to explore supersaturation, nucleation, and crystal growth. The effect of drug loading on physical stability and supersaturation potential was highly drug specific. Storage under humid conditions influenced crystallization, but also resulted in morphological changes and particle fusion. This led to increased particle size, which significantly reduced dissolution rate. MD simulations identified the importance of nano-compartmentalization in the crystallization rate of the ASDs. Nucleation during storage did not inherently compromise the ASD. Rather, the poorer performance resulted from a combination of properties of the compound, nanostructures formed in the formulation, and crystallization.

Molecular Dynamics Simulations Reveal Membrane Interactions for Poorly Water-Soluble Drugs: Impact of Bile Solubilization and Drug Aggregation

Molecular transport mechanisms of poorly soluble hydrophobic drug compounds to lipid membranes were investigated using molecular dynamics (MD) simulations. The model compound danazol was used to investigate the mechanism(s) by which bile micelles delivered it to the membrane. The interactions between lipid membrane and pure drug aggregates-in the form of amorphous aggregates and nanocrystals-were also studied. Our simulations indicate that bile micelles formed in the intestinal fluid may facilitate danazol incorporation into cellular membranes through two different mechanisms. The micelle may be acting as: i) a shuttle that presents the danazol directly to the membrane or ii) an elevator that moves the solubilized danazol with it as the colloidal structure itself becomes incorporated and solubilized within the membrane. The elevator hypothesis was supported by complementary lipid monolayer adsorption experiments. In these experiments, colloidal structures formed with simulated intestinal fluid were observed to rapidly incorporate into the monolayer. Simulations of membrane interaction with drug aggregates showed that both the amorphous aggregates and crystalline nanostructures incorporated into the membrane. However, the amorphous aggregates solubilized more quickly than the nanocrystals into the membrane, thereby improving the danazol absorption.

Molecular simulation as a computational pharmaceutics tool to predict drug solubility, solubilization processes and partitioning

In this review we will discuss how computational methods, and in particular classical molecular dynamics simulations, can be used to calculate solubility of pharmaceutically relevant molecules and systems. To the extent possible, we focus on the non-technical details of these calculations, and try to show also the added value of a more thorough and detailed understanding of the solubilization process obtained by using computational simulations. Although the main focus is on classical molecular dynamics simulations, we also provide the reader with some insights into other computational techniques, such as the COSMO-method, and also discuss Flory-Huggins theory and solubility parameters. We hope that this review will serve as a valuable starting point for any pharmaceutical researcher, who has not yet fully explored the possibilities offered by computational approaches to solubility calculations.

Hydrodynamic lift forces on solutes in a tilted nanopillar array: A computer simulation study

We study solute transport in a microfluidic channel, where the walls hold an array of tilted rigid nanopillars. By solving numerically the flow equations in the channel, we show that a combination of hydrodynamic effects with excluded volume interactions between the solute particles and the pillars leads to a hydrodynamic lift effect, which varies with the particle size, and depends in a strongly nonlinear fashion on the flow rate. We show that the lift force can be sufficiently strong to drive the solute accumulation or removal from the pillar region and can be switched to the opposite direction by variation of the shear rate or driving pressure. We also demonstrate that the nanopillar array can be used to selectively attract particles of certain size and enhance solute trapping at the surface.