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.

Intrinsic lipolysis rate for systematic design of lipid-based formulations.

Lipid-based formulations (LBFs) are used by the pharmaceutical industry in oral delivery systems for both poorly water-soluble drugs and biologics. Digestibility is key for the performance of LBFs and in vitro lipolysis is commonly used to compare the digestibility of LBFs. Results from in vitro lipolysis experiments depend highly on the experimental conditions and formulation characteristics, such as droplet size (which defines the surface area available for digestion) and interfacial structure. This study introduced the intrinsic lipolysis rate (ILR) as a surface area-independent approach to compare lipid digestibility. Pure acylglycerol nanoemulsions, stabilized with polysorbate 80 at low concentration, were formulated and digested according to a standardized pH-stat lipolysis protocol. A methodology originally developed to calculate the intrinsic dissolution rate of poorly water-soluble drugs was adapted for the rapid calculation of ILR from lipolysis data. The impact of surfactant concentration on the apparent lipolysis rate and lipid structure on ILR was systematically investigated. The surfactant polysorbate 80 inhibited lipolysis of tricaprylin nanoemulsions in a concentration-dependent manner. Coarse-grained molecular dynamics simulations supported these experimental observations. In the absence of bile and phospholipids, tricaprylin was shielded from lipase at 0.25% polysorbate 80. In contrast, the inclusion of bile salt and phospholipid increased the surfactant-free area and improved the colloidal presentation of the lipids to the enzyme, especially at 0.125% polysorbate 80. At a constant and low surfactant content, acylglycerol digestibility increased with decreasing acyl chain length, decreased esterification, and increasing unsaturation. The calculated ILR of pure acylglycerols was successfully used to accurately predict the IRL of binary lipid mixtures. The ILR measurements hold great promise as an efficient method supporting pharmaceutical formulation scientists in the design of LBFs with specific digestion profiles.

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.

Endothelial glycocalyx plays multiple roles in the vascular system: it regulates vascular permeability, modulates the interactions between blood and endothelial cells and controls the shear stress produced by the blood flow. The defense it provides against nanoscale blood solutes such as viruses and nanoparticles is based on nonspecific interactions. Being a natural sieve, it also influences the delivery of medicines transported by the blood flow. In this work, using a computational study of the permeability of endothelial glycocalyx, we demonstrate that the mesh formed by glycan fibers forms an effective mechanical barrier against penetration by nanoscale solutes. We evaluate the characteristic energies and times for the penetration of solutes of various sizes through the glycocalyx and assess the role of different contributions to nanoparticle­glycocalyx interactions. The obtained results may help clarify glycocalyx dysfunction and design drug nanocarriers and other nanoparticle-based medical applications.

Membrane insertion mechanism of the caveola coat protein Cavin1.

Caveolae are small plasma membrane invaginations, important for control of membrane tension, signaling cascades, and lipid sorting. The caveola coat protein Cavin1 is essential for shaping such high curvature membrane structures. Yet, a mechanistic understanding of how Cavin1 assembles at the membrane interface is lacking. Here, we used model membranes combined with biophysical dissection and computational modeling to show that Cavin1 inserts into membranes. We establish that initial phosphatidylinositol (4, 5) bisphosphate [PI(4,5)P2]-dependent membrane adsorption of the trimeric helical region 1 (HR1) of Cavin1 mediates the subsequent partial separation and membrane insertion of the individual helices. Insertion kinetics of HR1 is further enhanced by the presence of flanking negatively charged disordered regions, which was found important for the coassembly of Cavin1 with Caveolin1 in living cells. We propose that this intricate mechanism potentiates membrane curvature generation and facilitates dynamic rounds of assembly and disassembly of Cavin1 at the membrane.

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.

Structure and elasticity of bush and brush-like models of the endothelial glycocalyx.

The endothelial glycocalyx (EG), a sugar-rich layer that lines the luminal surface of blood vessels, is an important constituent of the vascular system. Although the chemical composition of the EG is fairly well known, there is no consensus regarding its ultrastructure. While previous experiments probed the properties of the layer at the continuum level, they did not provide sufficient insight into its molecular organisation. In this work, we investigate the EG mechanics using two simple brush and bush-like simulation models, and use these models to describe its molecular structure and elastic response to indentation. We analyse the relationship between the mechanical properties of the EG layer and several molecular parameters, including the filament bending rigidity, grafting density, and the type of ultrastructure . We show that variations in the glycan density determine the elasticity of the EG for small deformations, and that the normal stress may be effectively dampened by the EG layer, preventing the stress from being transferred to the cell membrane. Furthermore, our bush-like model allows us to evaluate the forces and energies required to overcome the mechanical resistance of the EG.

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.