Revealing the interaction between peptide drugs and permeation enhancers in the presence of intestinal bile salts.

Permeability enhancer-based formulations offer a promising approach to enhance the oral bioavailability of peptides. We used all-atom molecular dynamics simulations to investigate the interaction between two permeability enhancers (sodium caprate, and SNAC), and four different peptides (octreotide, hexarelin, degarelix, and insulin), in the presence of taurocholate, an intestinal bile salt. The permeability enhancers exhibited distinct effects on peptide release based on their properties, promoting hydrophobic peptide release while inhibiting water-soluble peptide release. Lowering peptide concentrations in the simulations reduced peptide-peptide interactions but increased their interactions with the enhancers and taurocholates. Introducing peptides randomly with enhancer and taurocholate molecules yielded dynamic molecular aggregation, and reduced peptide-peptide interactions and hydrogen bond formation compared to peptide-only systems. The simulations provided insights into molecular-level interactions, highlighting the specific contacts between peptide residues responsible for aggregation, and the interactions between peptide residues and permeability enhancers/taurocholates that are crucial within the mixed colloids. Therefore, our results can provide insights into how modifications of these critical contacts can be made to alter drug release profiles from peptide-only or mixed peptide-PE-taurocholate aggregates. To further probe the molecular nature of permeability enhancers and peptide interactions, we also analyzed insulin secondary structures using Fourier transform infrared spectroscopy. The presence of SNAC led to an increase in ß-sheet formation in insulin. In contrast, both in the absence and presence of caprate, α-helices, and random structures dominated. These molecular-level insights can guide the design of improved permeability enhancer-based dosage forms, allowing for precise control of peptide release profiles near the intended absorption site.

Numerical simulation of peristalsis to study co-localization and intestinal distribution of a macromolecular drug and permeation enhancer.

In this work, simulations of intestinal peristalsis are performed to investigate the intraluminal transport of macromolecules (MMs) and permeation enhancers (PEs). Properties of insulin and sodium caprate (C10) are used to represent the general class of MM and PE molecules. Nuclear magnetic resonance spectroscopy was used to obtain the diffusivity of C10, and coarse-grain molecular dynamics simulations were carried out to estimate the concentration-dependent diffusivity of C10. A segment of the small intestine with the length of 29.75 cm was modeled. Peristaltic speed, pocket size, release location, and occlusion ratio of the peristaltic wave were varied to study the effect on drug transport. It was observed that the maximum concentration at the epithelial surface for the PE and the MM increased by 397 % and 380 %, respectively, when the peristaltic wave speed was decreased from 1.5 to 0.5 cm s-1. At this wave speed, physiologically relevant concentrations of PE were found at the epithelial surface. However, when the occlusion ratio is increased from 0.3 to 0.7, the concentration approaches zero. These results suggest that a slower-moving and more contracted peristaltic wave leads to higher efficiency in transporting mass to the epithelial wall during the peristalsis phases of the migrating motor complex.

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.

Aggregation Behavior of Structurally Similar Therapeutic Peptides Investigated by 1H NMR and All-Atom Molecular Dynamics Simulations.

Understanding of peptide aggregation propensity is an important aspect in pharmaceutical development of peptide drugs. In this work, methodologies based on all-atom molecular dynamics (AA-MD) simulations and 1H NMR (in neat H2O) were evaluated as tools for identification and investigation of peptide aggregation. A series of structurally similar, pharmaceutically relevant peptides with known differences in aggregation behavior (D-Phe6-GnRH, ozarelix, cetrorelix, and degarelix) were investigated. The 1H NMR methodology was used to systematically investigate variations in aggregation with peptide concentration and time. Results show that 1H NMR can be used to detect the presence of coexisting classes of aggregates and the inclusion or exclusion of counterions in peptide aggregates. Interestingly, results suggest that the acetate counterions are included in aggregates of ozarelix and cetrorelix but not in aggregates of degarelix. The peptides investigated in AA-MD simulations (D-Phe6-GnRH, ozarelix, and cetrorelix) showed the same rank order of aggregation propensity as in the NMR experiments. The AA-MD simulations also provided molecular-level insights into aggregation dynamics, aggregation pathways, and the influence of different structural elements on peptide aggregation propensity and intermolecular interactions within the aggregates. Taken together, the findings from this study illustrate that 1H NMR and AA-MD simulations can be useful, complementary tools in early evaluation of aggregation propensity and formulation development for peptide drugs.

In Silico-Based Experiments on Mechanistic Interactions between Several Intestinal Permeation Enhancers with a Lipid Bilayer Model.

Oral administration of drugs is generally considered convenient and patient-friendly. However, oral administration of biological drugs exhibits low oral bioavailability (BA) due to enzymatic degradation and low intestinal absorption. A possible approach to circumvent the low BA of oral peptide drugs is to coformulate the drugs with permeation enhancers (PEs). PEs have been studied since the 1960s and are molecules that enhance the absorption of hydrophilic molecules with low permeability over the gastrointestinal epithelium. In this study, we investigated the impact of six PEs on the structural properties of a model membrane using molecular dynamics (MD) simulations. The PEs included were the sodium salts of the medium chain fatty acids laurate, caprate, and caprylate and the caprylate derivative SNAC─all with a negative charge─and neutral caprate and neutral sucrose monolaurate. Our results indicated that the PEs, once incorporated into the membrane, could induce membrane leakiness in a concentration-dependent manner. Our simulations suggest that a PE concentration of at least 70-100 mM is needed to strongly affect transcellular permeability. The increased aggregation propensity seen for neutral PEs might provide a molecular-level mechanism for the membrane disruptions seen at higher concentrations in vivo. The ability for neutral PEs to flip-flop across the lipid bilayer is also suggestive of possible intracellular modes of action other than increasing membrane fluidity. Taken together, our results indicate that MD simulations are useful for gaining insights relevant to the design of oral dosage forms based around permeability enhancer molecules.

TIRF Microscopy-Based Monitoring of Drug Permeation Across a Lipid Membrane Supported on Mesoporous Silica.

There is an urgent demand for analytic approaches that enable precise and representative quantification of the transport of biologically active compounds across cellular membranes. In this study, we established a new means to monitor membrane permeation kinetics, using total internal reflection fluorescence microscopy confined to a ≈500 nm thick mesoporous silica substrate, positioned underneath a planar supported cell membrane mimic. This way, we demonstrate spatiotemporally resolved membrane permeation kinetics of a small-molecule model drug, felodipine, while simultaneously controlling the integrity of, and monitoring the drug binding to, the cell membrane mimic. By contrasting the permeation behaviour of pure felodipine with felodipine coupled to the permeability enhancer caprylate (C8), we provide evidence for C8-facilitated transport across lipid membranes, thus validating the potential for this approach to successfully quantify carrier system-induced changes to cellular membrane permeation.

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.

Influence of Bile Composition on Membrane Incorporation of Transient Permeability Enhancers.

Transient permeability enhancers (PEs), such as caprylate, caprate, and salcaprozate sodium (SNAC), improve the bioavailability of poorly permeable macromolecular drugs. However, the effects are variable across individuals and classes of macromolecular drugs and biologics. Here, we examined the influence of bile compositions on the ability of membrane incorporation of three transient PEs-caprylate, caprate, and SNAC-using coarse-grained molecular dynamics (CG-MD). The availability of free PE monomers, which are important near the absorption site, to become incorporated into the membrane was higher in fasted-state fluids than that in fed-state fluids. The simulations also showed that transmembrane perturbation, i.e., insertion of PEs into the membrane, is a key mechanism by which caprylate and caprate increase permeability. In contrast, SNAC was mainly adsorbed onto the membrane surface, indicating a different mode of action. Membrane incorporation of caprylate and caprate was also influenced by bile composition, with more incorporation into fasted- than fed-state fluids. The simulations of transient PE interaction with membranes were further evaluated using two experimental techniques: the quartz crystal microbalance with dissipation technique and total internal reflection fluorescence microscopy. The experimental results were in good agreement with the computational simulations. Finally, the kinetics of membrane insertion was studied with CG-MD. Variation in micelle composition affected the insertion rates of caprate monomer insertion and expulsion from the micelle surface. In conclusion, this study suggests that the bile composition and the luminal composition of the intestinal fluid are important factors contributing to the interindividual variability in the absorption of macromolecular drugs administered with transient PEs.

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.