Molecular dynamics simulations of glycocholate-oleic acid mixed micelle assembly

Fundamental Aspects of Lipid-Based Excipients in Lipid-Based Product Development

Poor aqueous solubility of drugs is still a foremost challenge in pharmaceutical product development. The use of lipids in designing formulations provides an opportunity to enhance the aqueous solubility and consequently bioavailability of drugs. Pre-dissolution of drugs in lipids, surfactants, or mixtures of lipid excipients and surfactants eliminate the dissolution/dissolving step, which is likely to be the rate-limiting factor for oral absorption of poorly water-soluble drugs. In this review, we exhaustively summarize the lipids excipients in relation to their classification, absorption mechanisms, and lipid-based product development. Methodologies utilized for the preparation of solid and semi-solid lipid formulations, applications, phase behaviour, and regulatory perspective of lipid excipients are discussed.

Characterization of the self-assembly and size dependent structural properties of dietary mixed micelles by molecular dynamics simulations

The bile salts and phospholipids are secreted by the gallbladder to form dietary mixed micelles in which the solvation of poorly absorbed lipophilic drugs and nutraceuticals take place. A comprehensive understanding of the micellization and structure of the mixed micelles are crucial to design effective delivery systems for such substances. In this study, the evolution of the dietary mixed micelle formation under physiologically relevant concentrations and the dependence of structural properties on micelle size were investigated through coarse-grained molecular dynamics simulations. The MARTINI force field was used to model cholate and POPC as the representative bile salt and phospholipid, respectively. The micellization behavior was similar under both fasted and fed state concentrations. Total lipids concentration and the micelle size did not affect the internal structure of the micelles. All the micelles were slightly ellipsoidal in shape independent of their size. The extent of deviation from spherical geometry was found to depend on the micellar POPC/cholate ratio. We also found that the surface and core packing density of the micelles increased with micelle size. The former resulted in more perpendicular alignments of cholates with respect to the surface, while the latter resulted in an improved alignment of POPC tails with the radial direction and more uniform core density.

Computational Models of the Gastrointestinal Environment. 1. The Effect of Digestion on the Phase Behavior of Intestinal Fluids

Improved models of the gastrointestinal environment have great potential to assist the complex process of drug formulation. Molecular dynamics (MD) is a powerful method for investigating phase behavior at a molecular level. In this study we use multiple MD simulations to calculate phase diagrams for bile before and after digestion. In these computational models, undigested bile is represented by mixtures of palmitoyl-oleoylphosphatidylcholine (POPC), sodium glycodeoxycholate (GDX), and water. Digested bile is modeled using a 1:1 mixture of oleic acid and palmitoylphosphatidylcholine (lysophosphatidylcholine, LPC), GDX, and water. The computational phase diagrams of undigested and digested bile are compared, and we describe the typical intermolecular interactions that occur between phospholipids and bile salts. The diffusion coefficients measured from MD simulation are compared to experimental diffusion data measured by DOSY-NMR, where we observe good qualitative agreement. In an additional set of simulations, the effect of different ionization states of oleic acid on micelle formation is investigated.

Computational prediction of formulation strategies for beyond-rule-of-5 compounds

The physicochemical properties of some contemporary drug candidates are moving towards higher molecular weight, and coincidentally also higher lipophilicity in the quest for biological selectivity and specificity. These physicochemical properties move the compounds towards beyond rule-of-5 (B-r-o-5) chemical space and often result in lower water solubility. For such B-r-o-5 compounds non-traditional delivery strategies (i.e. those other than conventional tablet and capsule formulations) typically are required to achieve adequate exposure after oral administration. In this review, we present the current status of computational tools for prediction of intestinal drug absorption, models for prediction of the most suitable formulation strategies for B-r-o-5 compounds and models to obtain an enhanced understanding of the interplay between drug, formulation and physiological environment. In silico models are able to identify the likely molecular basis for low solubility in physiologically relevant fluids such as gastric and intestinal fluids. With this baseline information, a formulation scientist can, at an early stage, evaluate different orally administered, enabling formulation strategies. Recent computational models have emerged that predict glass-forming ability and crystallisation tendency and therefore the potential utility of amorphous solid dispersion formulations. Further, computational models of loading capacity in lipids, and therefore the potential for formulation as a lipid-based formulation, are now available. Whilst such tools are useful for rapid identification of suitable formulation strategies, they do not reveal drug localisation and molecular interaction patterns between drug and excipients. For the latter, Molecular Dynamics simulations provide an insight into the interplay between drug, formulation and intestinal fluid. These different computational approaches are reviewed. Additionally, we analyse the molecular requirements of different targets, since these can provide an early signal that enabling formulation strategies will be required. Based on the analysis we conclude that computational biopharmaceutical profiling can be used to identify where non-conventional gateways, such as prediction of 'formulate-ability' during lead optimisation and early development stages, are important and may ultimately increase the number of orally tractable contemporary targets.

Mechanism of enhanced antiosteoporosis effect of circinal-icaritin by self-assembled nanomicelles in vivo with suet oil and sodium deoxycholate

BACKGROUND

Circinal-icaritin (CIT), one new active aglycone of Epimedium, can exert a beneficial effect on osteoporotic bone. However, its low bioavailability limits its clinical efficacy for the treatment of osteoporosis.

MATERIALS AND METHODS

In this paper, suet oil (SO) was used to improve the oral bioavailability of CIT and enhance its antiosteoporosis effect and absorption. After oral administration of CIT together with SO, the CIT and SO self-assembled into nanomicelles under the action of sodium deoxycholate (DOC) by bile secretion. The antiosteoporosis effects of the CIT-SO-DOC nanomicelles were evaluated in osteoporotic rats by bone mineral density, serum biochemical markers, bone microarchitecture, bone biomechanical properties, and related protein and gene expressions. We examined the bioavailability of CIT and its nanomicelles in vivo, and subsequently the nanomicelles were verified using transmission electron microscopy. Finally, we evaluated absorption across a rat intestinal perfusion model.

RESULTS

Compared with CIT, in the CIT-SO groups, protein and messenger ribonucleic acid expressions of osteoprotegerin were increased, while expressions of receptor activator of nuclear factor-κB ligand in bone tissue were decreased; bone-turnover markers in serum of hydroxyproline, alkaline phosphatase, tartrate-resistant acid phosphatase 5b, and receptor activator of nuclear factor-κB ligand levels were decreased, while osteoprotegerin and osteocalcin levels were increased; and trabecular bone mass, microarchitecture, and bone biomechanical strength were enhanced. The relative bioavailabilities of CIT-SO high dosage, CIT-SO medium dosage, and CIT-SO low dosage (area under concentration-time curve [AUC]0-∞) compared with that of raw CIT high dosage, CIT medium dosage, and CIT low dosage (AUC0-∞) were 127%, 121%, and 134%, respectively. The average particle size of CIT-DOC was significantly decreased after adding SO (P<0.01), and the intestinal permeability coefficients of CIT-SO-DOC nanomicelles in the duodenum, jejunum, ileum, and colon were all significantly improved (P<0.01).

CONCLUSION

The increased antiosteoporosis effects and bioavailability of CIT-SO-DOC self-assembled nanomicelles were due to an increase in absorption of CIT by reducing the particle sizes of CIT. SO may be a practical oral carrier for antiosteoporosis drugs with low bioavailability.

Constant pH simulations with the coarse-grained MARTINI model — Application to oleic acid aggregates

Long chain fatty acids are biologically important molecules with complex and pH sensitive aggregation behavior. The carboxylic head group of oleic acid is ionizable, with the pKa shifting to larger values, even above a value of 7, in certain aggregate states. While experiments have determined the macroscopic phase behavior, we have yet to understand the molecular level details for this complex behavior. This level of detail is likely required to fully appreciate the role of fatty acids in biology and for nanoscale biotechnological and industrial applications. Here, we introduce the use of constant pH molecular dynamics (MD) simulations with the coarse-grained MARTINI model and apply the method to oleic acid aggregates and a model lipid bilayer. By running simulations at different constant pH values, we determined titration curves and the resulting pKa for oleic acid in different environments. The coarse-grained model predicts positive pKa shifts, with a shift from 4.8 in water to 6.5 in a small micelle, and 6.6 in a dioleoylphosphatidylcholine (DOPC) bilayer, similar to experimental estimates. The size of the micelles increased as the pH increased, and correlated with the fraction of deprotonated oleic acid. We show this combination of constant pH MD and the coarse-grained MARTINI model can be used to model pH-dependent surfactant phase behavior. This suggests a large number of potential new applications of large-scale MARTINI simulations in other biological systems with ionizable molecules.

Understanding pharmacokinetic food effects using molecular dynamics simulation coupled with physiologically based pharmacokinetic modeling

In this study, a molecular dynamics (MD) method is assessed as a new front-end tool for deriving relevant drug-micelle partitioning rates for use in conjunction with a compartmental-style gastrointestinal absorption model. A refined mechanistic approach for handling micelle-associated vs unbound drug is presented and examined in terms of its utility for projecting human oral pharmacokinetic food/formulation effects. Similar to predecessor oral absorption models, the intestinal drug absorption rate is formulated as a function of the combined permeability through the unstirred water layer and the epithelial membrane, however, an additional diffusion coefficient adjustment is applied to account for the viscosity changes of the postprandial small intestine. Bulk passage of drug particles through the GI tract is simulated by compartmental transit through a network of nine compartments comprising the stomach, small intestine and colon. The tandem MD simulation/compartmental absorption algorithm is applied to identify factors influencing the fed vs fasted absorption ratios of a structurally diverse set of orally administered drugs. The data illustrate the interplay and apparent compromise between micelle solubilization and permeability.

Aggregation behavior of ibuprofen, cholic acid and dodecylphosphocholine micelles

Non-steroidal anti-inflammatory drugs (NSAIDs) are frequently used to treat chronic pain and inflammation. However, prolonged use of NSAIDs has been known to result in Gastrointestinal (GI) ulceration/bleeding, with a bile-mediated mechanism underlying their toxicity to the lower gut. Bile acids (BAs) and phosphatidylcholines (PCs), the major components of bile, form mixed micelles to reduce the membrane disruptive actions of monomeric BAs and simple BA micelles. NSAIDs are suspected to alter the BA/PC balance in the bile, but the molecular interactions of NSAID-BA or NSAID-BA-PC remain undetermined. In this work, we used a series of all-atom molecular dynamics simulations of cholic acid (CA), ibuprofen (IBU) and dodecylphosphocholine (DPC) mixtures to study the spontaneous aggregation of CA and IBU as well as their adsorption on a DPC micelle. We found that the size of CA-IBU mixed micelles varies with their molar ratio in a non-linear manner, and that micelles of different sizes adopt similar shapes but differ in composition and internal interactions. These observations are supported by NMR chemical shift changes, NMR ROESY crosspeaks between IBU and CA, and dynamic light scattering experiments. Smaller CA-IBU aggregates were formed in the presence of a DPC micelle due to the segregation of CA and IBU away from each other by the DPC micelle. While the larger CA-IBU aggregates arising from higher IBU concentrations might be responsible for NSAID-induced intestinal toxicity, the absence of larger CA-IBU aggregates in the presence of DPC micelles may explain the observed attenuation of NSAID toxicity by PCs.

Using molecular dynamics to study liquid phase behavior: simulations of the ternary sodium laurate/sodium oleate/water system

The prediction of surfactant phase behavior has applications in a wide range of areas. An accurate modeling of liquid phase behavior can aid our understanding of colloidal process or be used to design phases that respond in a defined way to their environment. In this work, we use molecular dynamics to model the phase behavior of the ternary sodium laurate/sodium oleate/water system and compare the simulation results to experimental data. Simulations were performed with the GROMOS 53A6 united-atom force field and cover the entire ternary phase diagram, producing micellar, hexagonal, and lamellar phases. The aggregate simulation time for the 33 simulations performed during this study is 4.4 μs. We find that the simulations were able to model the experimentally observed liquid phase behavior accurately, showing that the carboxylate and lipid parameters of the 53A6 force field give very good quality results for the in silico prediction of liquid system phase behavior.

Practical considerations for building GROMOS-compatible small-molecule topologies

Molecular dynamics simulations are being applied to increasingly complex systems, including those involving small endogenous compounds and drug molecules. In order to obtain meaningful and accurate data from these simulations, high-quality topologies for small molecules must be generated in a manner that is consistent with the derivation of the force field applied to the system. Often, force fields are designed for use with macromolecules such as proteins, making their transferability to other species challenging. Investigators are increasingly attracted to automated topology generation programs, although the quality of the resulting topologies remains unknown. Here we assess the applicability of the popular PRODRG server that generates small-molecule topologies for use with the GROMOS family of force fields. We find that PRODRG does not reproduce topologies for even the most well-characterized species in the force field due to inconsistent charges and charge groups. We assessed the effects of PRODRG-derived charges on several systems: pure liquids, amino acids at a hydrophobic-hydrophilic interface, and an enzyme-cofactor complex. We found that partial atomic charges generated by PRODRG are largely incompatible with GROMOS force fields, and the behavior of these systems deviates substantially from that of simulations using GROMOS parameters. We conclude by proposing several points as "best practices" for parametrization of small molecules under the GROMOS force fields.