Molecular Dynamics Simulations of Self-Assembling Colloids in Fed-State Human Intestinal Fluids and Their Solubilization of Lipophilic Drugs

Understanding drug solubilization in intestinal mixed micelles through molecular dynamics simulations

HYPOTHESIS

Solubilization is a fundamental process that underpins various technologies in the pharmaceutical and chemical industry. However, knowledge of the location, orientation and interactions of solubilized molecules in the micelles is still limited. We expect all-atom molecular dynamics simulations to improve the molecular-level understanding of solubilization and to enable its in silico prediction.

METHODS

The solubilization of six drugs in intestinal mixed micelles composed of taurocholate and dioleoyl phosphatidylcholine was simulated by molecular dynamics in explicit water and measured experimentally by liquid chromatography. The location and orientation of the solubilized drugs were visualized by cumulative radial distribution functions and interactions were characterized by radial distribution function ratios and hydrogen bonding.

FINDINGS

A new simulation-derived parameter was defined, which accounts for drug-micelle and drug-water interactions and correlates (R2 = 0.83) with the experimentally measured solubilization. Lipophilicity was found to govern the location of all drugs in the micelle (hydrophobic core, palisade layer or on the surface), while hydrogen bonding was crucial for orientation and solubilization of two of the molecules. The study demonstrates that explicit, hydrogen bond-forming water molecules are vital for accurate prediction of solubilization and provides a comprehensive framework for quantitative studies of drug location and orientation within the micelles.

Roles of Supersaturation and Liquid-Liquid Phase Separation for Enhanced Oral Absorption of Poorly Soluble Drugs from Amorphous Solid Dispersions

Amorphous solid dispersion (ASD) is one of the most important enabling formulation technologies for the development of poorly soluble drugs. Because of its thermodynamically unstable nature in both solid and wet states, the evaluation and optimization of the formulation performance involves some difficulties. The dissolution process is sensitively influenced by various factors, including the applied dose, medium composition, and pH. Supersaturated solutions can cause liquid-liquid phase separation (LLPS) and/or crystallization, which complicates the comprehension of the dissolution process. However, LLPS should be evaluated carefully because it is closely related to oral absorption. As LLPS concentration is analogous to amorphous solubility, it can be a key factor in predicting oral absorption from ASDs, if absorption is limited by solubility. Moreover, LLPS droplets are expected to increase transmembrane flux by increasing the drug concentration near the epithelial cell membrane. In this review, recently updated knowledge on the dissolution, membrane permeation, and oral absorption behaviors of ASDs is discussed with an emphasis on LLPS behavior.

Bridging the gap between in vitro and in vivo solubility-permeability interplay

Use of solubilization carriers for poorly soluble drugs may disturb transmembrane absorption by lowering the activity of drug molecules, which is known as the solubility-permeability interplay. However, although many in vitro studies have indicated the negative impacts of use of solubilization carriers for oral absorption, in vivo studies that showed the interplay effect are limited. This study provides systematic in vitro, in situ, and in vivo investigation of the interplay effect of cyclodextrin using dexamethasone as a model drug. The evaluation methods included permeation through polymeric, artificial lipid, cell, and intestinal closed-loop membranes. Then, the results were compared with oral administration studies in mice and dogs. Although the interplay effect was clearly observed in the in vitro studies, no obvious interplay was found in the in vivo studies, suggesting that the interplay effect is more prominent in the in vitro permeation studies. Absence of in vivo interplay was attributed to the dilution effect in the gastrointestinal tract, interaction of the drug with living components, and clearance of the drug after membrane permeation. Overall, this investigation clearly demonstrated the applicability and limitations of in vitro permeation studies for predicting the interplay effects of solubilizers after the oral administration.

Coarse-Grained Molecular Dynamics Simulation of Heterogeneous Polysorbate 80 Surfactants and their Interactions with Small Molecules and Proteins

Polysorbate 80 (PS80) is widely used in pharmaceutical formulations, and its commercial grades exhibit certain levels of structural heterogeneity. The objective of this study was to apply coarse-grained molecular dynamics simulations to better understand the effect of PS80 heterogeneity on micelle self-assembly, the loading of hydrophobic small molecules into the micelle core, and the interactions between PS80 and a protein, bovine serum albumin (BSA). Four representative PS80 variants with different head and tail structures were studied. Our simulations found that PS80 structural heterogeneity could affect blank micelle properties such as solvent-accessible surface area, aggregation number, and micelle aspect ratio. It was also found that hydrophobic small molecules such as ethinyl estradiol preferentially partitioned into the PS80 micelle core and PS80 dioleates formed a more hydrophobic core compared to PS80 monooleates. Furthermore, multiple PS80 molecules could bind to BSA, and PS80 heterogeneity profoundly changed the binding ratio as well as the surfactant-protein contact area.

Development of nano-delivery systems for loaded bioactive compounds: using molecular dynamics simulations

Over the past decade, a remarkable surge in the development of functional nano-delivery systems loaded with bioactive compounds for healthcare has been witnessed. Notably, the demanding requirements of high solubility, prolonged circulation, high tissue penetration capability, and strong targeting ability of nanocarriers have posed interdisciplinary research challenges to the community. While extensive experimental studies have been conducted to understand the construction of nano-delivery systems and their metabolic behavior in vivo, less is known about these molecular mechanisms and kinetic pathways during their metabolic process in vivo, and lacking effective means for high-throughput screening. Molecular dynamics (MD) simulation techniques provide a reliable tool for investigating the design of nano-delivery carriers encapsulating these functional ingredients, elucidating the synthesis, translocation, and delivery of nanocarriers. This review introduces the basic MD principles, discusses how to apply MD simulation to design nanocarriers, evaluates the ability of nanocarriers to adhere to or cross gastrointestinal mucosa, and regulates plasma proteins in vivo. Moreover, we presented the critical role of MD simulation in developing delivery systems for precise nutrition and prospects for the future. This review aims to provide insights into the implications of MD simulation techniques for designing and optimizing nano-delivery systems in the healthcare food industry.

Nanoscale Prediction of the Thermal, Mechanical, and Transport Properties of Hydrated Clay on 106- and 1015-Fold Larger Length and Time Scales

Coupled thermal, hydraulic, mechanical, and chemical (THMC) processes, such as desiccation-driven cracking or chemically driven fluid flow, significantly impact the performance of composite materials formed by fluid-mediated nanoparticle assembly, including energy storage materials, ordinary Portland cement, bioinorganic nanocomposites, liquid crystals, and engineered clay barriers used in the isolation of hazardous wastes. These couplings are particularly important in the isolation of high-level radioactive waste (HLRW), where heat generated by radioactive decay can drive the temperature up to at least 373 K in the engineered barrier. Here, we use large-scale all-atom molecular dynamics simulations of hydrated smectite clay nanoparticle assemblages to predict the fundamental THMC properties of hydrated compacted clay over a wide range of temperatures (up to 373 K) and dry densities relevant to HLRW management. Equilibrium simulations of clay-water mixtures at different hydration levels are analyzed to quantify material properties, including thermal conductivity, heat capacity, thermal expansion, suction, water and ion self-diffusivity, and hydraulic conductivity. Predictions are validated against experimental results for the properties of compacted bentonite clay. Our results demonstrate the feasibility of using atomistic-level simulations of assemblages of clay nanoparticles on scales of tens of nanometers and nanoseconds to infer the properties of compacted bentonite on scales of centimeters and days, a direct upscaling over 6 orders of magnitude in space and 15 orders of magnitude in time.