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

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.European Commission - Seventh Framework Programme (FP7)Irish Research CouncilIntel Irelan

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

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

Explicit-pH Coarse-Grained Molecular Dynamics Simulations Enable Insights into Restructuring of Intestinal Colloidal Aggregates with Permeation Enhancers

Permeation enhancers (PEs) can increase the bioavailability of drugs. The mechanisms of action of these PEs are complex, but, typically, when used for oral administration, they can transiently induce the alteration of trans- and paracellular pathways, including increased solubilization and membrane fluidity, or the opening of the tight junctions. To elucidate these mechanistic details, it is important to understand the aggregation behavior of not only the PEs themselves but also other molecules already present in the intestine. Aggregation processes depend critically on, among other factors, the charge state of ionizable chemical groups, which is affected by the pH of the system. In this study, we used explicit-pH coarse-grained molecular dynamics simulations to investigate the aggregation behavior and pH dependence of two commonly used PEs—caprate and SNAC—together with other components of fasted- and fed-state simulated intestinal fluids. We also present and validate a coarse-grained molecular topology for the bile salt taurocholate suitable for the Martini3 force-field. Our results indicate an increase in the number of free molecules as a function of the system pH and for each combination of FaSSIF/FeSSIF and PEs. In addition, there are differences between caprate and SNAC, which are rationalized based on their different molecular structures and critical micelle concentrations

Explicit-pH Coarse-Grained Molecular Dynamics Simulations Enable Insights into Restructuring of Intestinal Colloidal Aggregates with Permeation Enhancers

Permeation enhancers (PEs) can increase the bioavailability of drugs. The mechanisms of action of these PEs are complex, but, typically, when used for oral administration, they can transiently induce the alteration of trans- and paracellular pathways, including increased solubilization and membrane fluidity, or the opening of the tight junctions. To elucidate these mechanistic details, it is important to understand the aggregation behavior of not only the PEs themselves but also other molecules already present in the intestine. Aggregation processes depend critically on, among other factors, the charge state of ionizable chemical groups, which is affected by the pH of the system. In this study, we used explicit-pH coarse-grained molecular dynamics simulations to investigate the aggregation behavior and pH dependence of two commonly used PEs—caprate and SNAC—together with other components of fasted- and fed-state simulated intestinal fluids. We also present and validate a coarse-grained molecular topology for the bile salt taurocholate suitable for the Martini3 force-field. Our results indicate an increase in the number of free molecules as a function of the system pH and for each combination of FaSSIF/FeSSIF and PEs. In addition, there are differences between caprate and SNAC, which are rationalized based on their different molecular structures and critical micelle concentrations

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

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)P-2]-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