POxload: Machine Learning Estimates Drug Loadings of Polymeric Micelles.

Block copolymers, composed of poly(2-oxazoline)s and poly(2-oxazine)s, can serve as drug delivery systems; they form micelles that carry poorly water-soluble drugs. Many recent studies have investigated the effects of structural changes of the polymer and the hydrophobic cargo on drug loading. In this work, we combine these data to establish an extended formulation database. Different molecular properties and fingerprints are tested for their applicability to serve as formulation-specific mixture descriptors. A variety of classification and regression models are built for different descriptor subsets and thresholds of loading efficiency and loading capacity, with the best models achieving overall good statistics for both cross- and external validation (balanced accuracies of 0.8). Subsequently, important features are dissected for interpretation, and the DrugBank is screened for potential therapeutic use cases where these polymers could be used to develop novel formulations of hydrophobic drugs. The most promising models are provided as an open-source software tool for other researchers to test the applicability of these delivery systems for potential new drug candidates.

Unravel the Tangle: Atomistic Insight into Ultrahigh Curcumin-Loaded Polymer Micelles.

Amphiphilic ABA-triblock copolymers, comprised of poly(2-oxazoline) and poly(2-oxazine), can solubilize poorly water-soluble molecules in a structure-dependent manner forming micelles with exceptionally high drug loading. All-atom molecular dynamics simulations are conducted on previously experimentally characterized, curcumin-loaded micelles to dissect the structure-property relationships. Polymer-drug interactions for different levels of drug loading and variation in polymer structures of both the inner hydrophobic core and outer hydrophilic shell are investigated. In silico, the system with the highest experimental loading capacity shows the highest number of drug molecules encapsulated by the core. Furthermore, in systems with lower loading capacity outer A blocks show a greater extent of entanglement with the inner B blocks. Hydrogen bond analyses corroborate previous hypotheses: poly(2-butyl-2-oxazoline) B blocks, found experimentally to have reduced loading capacity for curcumin compared to poly(2-propyl-2-oxazine), establish fewer but longer-lasting hydrogen bonds. This possibly results from different sidechain conformations around the hydrophobic cargo, which is investigated by unsupervised machine learning to cluster monomers in smaller model systems mimicking different micelle compartments. Exchanging poly(2-methyl-2-oxazoline) with poly(2-ethyl-2-oxazoline) leads to increased drug interactions and reduced corona hydration; this suggests an impairment of micelle solubility or colloidal stability. These observations can help driving forward a more rational a priori nanoformulation design.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations.

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE-CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

In Vitro and in Vivo Behavior of Liposomes Decorated with PEGs with Different Chemical Features.

The colloidal stability, in vitro toxicity, cell association, and in vivo pharmacokinetic behavior of liposomes decorated with monomethoxy-poly(ethylene glycol)-lipids (mPEG-lipids) with different chemical features were comparatively investigated. Structural differences of the mPEG-lipids used in the study included: (a) surface-anchoring moiety [1,2-distearoyl-sn-glycero-3-phosphoethanolamine (DSPE), cholesterol (Chol), and cholane (Chln)]; (b) mPEG molecular weight (2 kDa mPEG45 and 5 kDa mPEG114); and (c) mPEG shape (linear and branched PEG). In vitro results demonstrated that branched (mPEG114)2-DSPE confers the highest stealth properties to liposomes (∼31-fold lower cell association than naked liposomes) with respect to all PEGylating agents tested. However, the pharmacokinetic studies showed that the use of cholesterol as anchoring group yields PEGylated liposomes with longer permeance in the circulation and higher systemic bioavailability among the tested formulations. Liposomes decorated with mPEG114-Chol had 3.2- and ∼2.1-fold higher area under curve (AUC) than naked liposomes and branched (mPEG114)2-DSPE-coated liposomes, respectively, which reflects the high stability of this coating agent. By comparing the PEGylating agents with same size, namely, linear 5 kDa PEG derivatives, linear mPEG114-DSPE yielded coated liposomes with the best in vitro stealth performance. Nevertheless, the in vivo AUC of liposomes decorated with linear mPEG114-DSPE was lower than that obtained with liposomes decorated with linear mPEG114-Chol. Computational molecular dynamics modeling provided additional insights that complement the experimental results.

Rigorous Computational Study Reveals What Docking Overlooks: Double Trouble from Membrane Association in Protein Kinase C Modulators.

Increasing protein kinase C (PKC) activity is of potential therapeutic value. Its activation involves an interaction between the C1 domain and diacylglycerol (DAG) at intracellular membrane surfaces; DAG mimetics hold promise as new drugs. We previously developed the isophthalate derivative HMI-1a3, an effective but highly lipophilic (clogP = 6.46) DAG mimetic. Although a less lipophilic pyrimidine analog, PYR-1gP (clogP = 3.30), gave positive results in computational docking, it unexpectedly presented greatly diminished binding to PKC in vitro. Through more rigorous computational molecular modeling, we reveal that, unlike HMI-1a3, PYR-1gP forms an intramolecular hydrogen bond, which both obstructs binding and reorients PYR-1gP in the membrane in a fashion that prevents it from correctly accessing the PKC C1 domain. Our results highlight the great value of molecular dynamics simulations as a key component for the drug design process of ligands targeting weakly membrane-associated proteins, where simulation in the relevant membrane environment is crucial for obtaining biologically applicable results.

Biophysical Characterization of Supported Lipid Bilayers Using Parallel Dual-Wavelength Surface Plasmon Resonance and Quartz Crystal Microbalance Measurements.

Supported lipid bilayers (SLBs) have been used extensively as an effective model of biological membranes, in the context of in vitro biophysics research, and the membranes of liposomes, in the context of the development of nanoscale drug delivery devices. Despite numerous surface-sensitive techniques having been applied to their study, the comprehensive optical characterization of SLBs using surface plasmon resonance (SPR) has not been conducted. In this study, Fresnel multilayer analysis is utilized to effectively calculate layer parameters (thickness and refractive indices) with the aid of dual-wavelength and dispersion coefficient analysis, in which the linear change in the refractive index as a function of wavelength is assumed. Using complementary information from impedance-based quartz crystal microbalance experiments, biophysical properties, for example, area-per-lipid-molecule and the quantity of lipid-associated water molecules, are calculated for different lipid types and mixtures, one of which is representative of a raft-forming lipid mixture. It is proposed that the hydration layer beneath the bilayer is, in fact, an integral part of the measured optical signal. Also, the traditional Jung model analysis and the ratio of SPR responses are investigated in terms of assessing the structure of the lipid layer that is formed.

Control of Peptide Aggregation and Fibrillation by Physical PEGylation.

Peptide therapeutics have the potential to self-associate, leading to aggregation and fibrillation. Noncovalent PEGylation offers a strategy to improve their physical stability; an understanding of the behavior of the resulting polymer/peptide complexes is, however, required. In this study, we have performed a set of experiments with additional mechanistic insight provided by in silico simulations to characterize the molecular organization of these complexes. We used palmitoylated vasoactive intestinal peptide (VIP-palm) stabilized by methoxy-poly(ethylene glycol)5kDa-cholane (PEG-cholane) as our model system. Homogeneous supramolecular assemblies were found only when complexes of PEG-cholane/VIP-palm exceeded a molar ratio of 2:1; at and above this ratio, the simulations showed minimal exposure of VIP-palm to the solvent. Supramolecular assemblies formed, composed of, on average, 9-11 PEG-cholane/VIP-palm complexes with 2:1 stoichiometry. Our in silico results showed the structural content of the helical conformation in VIP-palm increases when it is complexed with the PEG-cholane molecule; this behavior becomes yet more pronounced when these complexes assemble into larger supramolecular assemblies. Our experimental results support this: the extent to which VIP-palm loses helical structure as a result of thermal denaturation was inversely related to the PEG-cholane:VIP-palm molar ratio. The addition of divalent buffer species and increasing the ionic strength of the solution both accelerate the formation of VIP-palm fibrils, which was partially and fully suppressed by 2 and >4 mol equivalents of PEG-cholane, respectively. We conclude that the relative freedom of the VIP-palm backbone to adopt nonhelical conformations is a key step in the aggregation pathway.

Lipid Architectonics for Superior Oral Bioavailability of Nelfinavir Mesylate: Comparative in vitro and in vivo Assessment.

Nelfinavir mesylate (NFV), a human immunodeficiency virus (HIV) protease inhibitor, is an integral component of highly active anti retro viral therapy (HAART) for management of AIDS. NFV possesses pH-dependent solubility and has low and variable bioavailability hampering its use in therapeutics. Lipid-based particulates have shown to improve solubility of poorly water soluble drugs and oral absorption, thereby aiding in improved bioavailability. The current study compares potential of vesicular and solid lipid nanocarriers of NFV with drug nanocrystallites and microvesicular systems like cochleates in improving bioavailability of NFV. The paper outlines investigation of systems using in vitro models like in vitro lipolysis, in vitro release, and permeation through cell lines to predict the in vivo potential of nanocarriers. Finally, in vivo pharmacokinetic study is reported which provided proof of concept in sync with results from in vitro studies. Graphical Abstract ᅟ.

Rational design of liposomal drug delivery systems, a review: Combined experimental and computational studies of lipid membranes, liposomes and their PEGylation.

Combined experimental and computational studies of lipid membranes and liposomes, with the aim to attain mechanistic understanding, result in a synergy that makes possible the rational design of liposomal drug delivery system (LDS) based therapies. The LDS is the leading form of nanoscale drug delivery platform, an avenue in drug research, known as "nanomedicine", that holds the promise to transcend the current paradigm of drug development that has led to diminishing returns. Unfortunately this field of research has, so far, been far more successful in generating publications than new drug therapies. This partly results from the trial and error based methodologies used. We discuss experimental techniques capable of obtaining mechanistic insight into LDS structure and behavior. Insight obtained purely experimentally is, however, limited; computational modeling using molecular dynamics simulation can provide insight not otherwise available. We review computational research, that makes use of the multiscale modeling paradigm, simulating the phospholipid membrane with all atom resolution and the entire liposome with coarse grained models. We discuss in greater detail the computational modeling of liposome PEGylation. Overall, we wish to convey the power that lies in the combined use of experimental and computational methodologies; we hope to provide a roadmap for the rational design of LDS based therapies. Computational modeling is able to provide mechanistic insight that explains the context of experimental results and can also take the lead and inspire new directions for experimental research into LDS development. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Effects of Membrane PEGylation on Entry and Location of Antifungal Drug Itraconazole and Their Pharmacological Implications.

Itraconazole (ITZ) is an antifungal agent used clinically to treat mycotic infections. However, its therapeutic effects are limited by low solubility in aqueous media. Liposome-based delivery systems (LDS) have been proposed as a delivery mechanism for ITZ to alleviate this problem. Furthermore, PEGylation, the inclusion in the formulation of a protective "stealth sheath" of poly(ethylene glycol) around carrier particles, is widely used to increase circulation time in the bloodstream and hence efficacy. Together, these themes highlight the importance of mechanistic and structural understanding of ITZ incorporation into liposomes both with and without PEGylation because it can provide a potential foundation for the rational design of LDS-based systems for delivery of ITZ, using alternate protective polymers or formulations. Here we have combined atomistic simulations, cryo-TEM, Langmuir film balance, and fluorescence quenching experiments to explore how ITZ interacts with both pristine and PEGylated liposomes. We found that the drug can be incorporated into conventional and PEGylated liposomes for drug concentrations up to 15 mol % without phase separation. We observed that, in addition to its protective properties, PEGylation significantly increases the stability of liposomes that host ITZ. In a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer without PEGylation, ITZ was found to reside inside the lipid bilayer between the glycerol and the double-bond regions of POPC, adopting a largely parallel orientation along the membrane surface. In a PEGylated liposome, ITZ partitions mainly to the PEG layer. The results provide a solid basis for further development of liposome-based delivery systems.

Interactions between Chloramphenicol, Carrier Polymers, and Bacteria-Implications for Designing Electrospun Drug Delivery Systems Countering Wound Infection.

Antibacterial drug-loaded electrospun nano- and microfibrous dressings are of major interest as novel topical drug delivery systems in wound care. In this study, chloramphenicol (CAM)-loaded polycaprolactone (PCL) and PCL/poly(ethylene oxide) (PEO) fiber mats were electrospun and characterized in terms of morphology, drug distribution, physicochemical properties, drug release, swelling, cytotoxicity, and antibacterial activity. Computational modeling together with physicochemical analysis helped to elucidate possible interactions between the drug and carrier polymers. Strong interactions between PCL and CAM together with hydrophobicity of the system resulted in much slower drug release compared to the hydrophilic ternary system of PCL/PEO/CAM. Cytotoxicity studies confirmed safety of the fiber mats to murine NIH 3T3 cells. Disc diffusion assay demonstrated that both fast and slow release fiber mats reached effective concentrations and had similar antibacterial activity. A biofilm formation assay revealed that both blank matrices are good substrates for the bacterial attachment and formation of biofilm. Importantly, prolonged release of CAM from drug-loaded fibers helps to avoid biofilm formation onto the dressing and hence avoids the treatment failure.

Indocyanine Green-Loaded Liposomes for Light-Triggered Drug Release.

Light-triggered drug delivery systems enable site-specific and time-controlled drug release. In previous work, we have achieved this with liposomes containing gold nanoparticles in the aqueous core. Gold nanoparticles absorb near-infrared light and release the energy as heat that increases the permeability of the liposomal bilayer, thus releasing the contents of the liposome. In this work, we replaced the gold nanoparticles with the clinically approved imaging agent indocyanine green (ICG). The ICG liposomes were stable at storage conditions (4-22 °C) and at body temperature, and fast near-infrared (IR) light-triggered drug release was achieved with optimized phospholipid composition and a 1:50 ICG-to-lipid molar ratio. Encapsulated small molecular calcein and FITC-dextran (up to 20 kDa) were completely released from the liposomes after light exposure for 15 s. Location of ICG in the PEG layer of the liposomes was simulated with molecular dynamics. ICG has important benefits as a light-triggering agent in liposomes: fast content release, improved stability, improved possibility of liposomal size control, regulatory approval to use in humans, and the possibility of imaging the in vivo location of the liposomes based on the fluorescence of ICG. Near-infrared light used as a triggering mechanism has good tissue penetration and safety. Thus, ICG liposomes are an attractive option for light-controlled and efficient delivery of small and large drug molecules.

Hydrophobin film structure for HFBI and HFBII and mechanism for accelerated film formation.

Hydrophobins represent an important group of proteins from both a biological and nanotechnological standpoint. They are the means through which filamentous fungi affect their environment to promote growth, and their properties at interfaces have resulted in numerous applications. In our study we have combined protein docking, molecular dynamics simulation, and electron cryo-microscopy to gain atomistic level insight into the surface structure of films composed of two class II hydrophobins: HFBI and HFBII produced by Trichoderma reesei. Together our results suggest a unit cell composed of six proteins; however, our computational results suggest P6 symmetry, while our experimental results show P3 symmetry with a unit cell size of 56 Å. Our computational results indicate the possibility of an alternate ordering with a three protein unit cell with P3 symmetry and a smaller unit cell size, and we have used a Monte Carlo simulation of a spin model representing the hydrophobin film to show how this alternate metastable structure may play a role in increasing the rate of surface coverage by hydrophobin films, possibly indicating a mechanism of more general significance to both biology and nanotechnology.

Activity and potential mechanisms of action of persimmon tannins according to their structures: A review.

One distinguishing feature of the persimmon, that differentiates it from other fruits, is its high proanthocyanidins content, known as persimmon tannin (PT). Despite the poor absorption of PT in the small intestine, results from animal studies demonstrate that PT has many health benefits. Our goal in this review is to summarize the literature that elucidates the relationship between PT structure and activity. In addition, we also summarize the potential mechanisms underlying the health benefits that result from PT consumption; this includes the hypolipidemic, hypoglycemic, antioxidant, anti-inflammatory, antiradiation, antibacterial and antiviral, detoxification effects on snake venom, and the absorption of heavy metals and dyes. Studies show that PT is a structurally distinct proanthocyanidins that exhibits a high degree of polymerization. It is galloylation-rich and possesses unique A-type interflavan linkages in addition to the more common B-type interflavan bonds. Thus, PT is converted into oligomeric proanthocyanidins by depolymerization strategies, including the nucleophilic substitution reaction, acid hydrolysis, and hydrogenolysis. In addition, multiple health benefits exerted by PT mainly involve the inactivation of lipogenic and intracellular inflammatory signaling pathways, activation of the fatty acid oxidation signaling pathway, regulation of gut microbiota, and highly absorptive properties.

Expanding the Paradigm of Structure-Based Drug Design: Molecular Dynamics Simulations Support the Development of New Pyridine-Based Protein Kinase C-Targeted Agonists.

Protein kinase C (PKC) modulators hold therapeutic potential for various diseases, including cancer, heart failure, and Alzheimer's disease. Targeting the C1 domain of PKC represents a promising strategy; the available protein structures warrant the design of PKC-targeted ligands via a structure-based approach. However, the PKC C1 domain penetrates the lipid membrane during binding, complicating the design of drug candidates. The standard docking-scoring approach for PKC lacks information regarding the dynamics and the membrane environment. Molecular dynamics (MD) simulations with PKC, ligands, and membranes have been used to address these shortcomings. Previously, we observed that less computationally intensive simulations of just ligand-membrane interactions may help elucidate C1 domain-binding prospects. Here, we present the design, synthesis, and biological evaluation of new pyridine-based PKC agonists implementing an enhanced workflow with ligand-membrane MD simulations. This workflow holds promise to expand the approach in drug design for ligands targeted to weakly membrane-associated proteins.

Strong preferences of dopamine and l-dopa towards lipid head group: importance of lipid composition and implication for neurotransmitter metabolism.

The interactions of the neurotransmitter dopamine, and its precursor l-dopa, with membrane lipids were investigated through a set of molecular dynamic simulations with all atom resolution. The results obtained indicate that both dopamine and l-dopa have a pronounced association with the lipid head groups, predominantly mediated through H-bonds. As a result the molecules are anchored to the interfacial region of the membrane. The strength of this interaction is dependent on lipid composition - the presence of phosphatidylserine leads to an increase in the strength of this interaction, resulting in an H-bond network with a lifetime much longer than the timescale of our simulations. Also, bilayers that include sphingomieline and cholesterol interact strongly with dopamine and l-dopa. We postulate that the high membrane association that we have observed for both dopamine and l-dopa could have the following effects: 1) when on the plasma membrane exterior, favour the availability of these compounds for cell membrane uptake processes and, 2) when on an internal membrane surface, accentuate the importance of membrane-bound metabolizing enzymes over their soluble counterparts.

Insights into the behavior of unsaturated diacylglycerols in mixed lipid bilayers in relation to protein kinase C activation-A molecular dynamics simulation study.

The lipid second messenger diacylglycerol (DAG) is known for its involvement in many types of cellular signaling, especially as an endogenous agonist for protein kinase C (PKC). Evidence has emerged that the degree of saturation of the DAG molecules can affect PKC activation. DAG molecules with different acyl chain saturation have not only been observed to induce varying extents of PKC activation, but also to express selectivity towards different PKC isozymes. Both qualities are important for precise therapeutic activation of PKC; understanding DAG behavior at the molecular level in different environments has much potential in the development of drugs to target PKC. We used molecular dynamics simulations to study the behavior of two different unsaturated DAG species in lipid environments with varying degrees of unsaturation. We focus on phosphatidylethanolamine (PE) instead of phosphatidylcholine (PC) to more accurately model the relevant biomembranes. The effect of cholesterol (CHOL) on these systems was also explored. We found that both high level of unsaturation in the acyl chains of the DAG species and presence of CHOL in the surrounding membrane increase DAG molecule availability at the lipid-water interface. This can partially explain the previously observed differences in PKC activation strength and specificity, the complete mechanism is, however, likely to be more complex. Our simulations coupled with the current understanding of lipids highlight the need for more simulations of biologically accurate lipid environments in order to determine the correct correlations between molecular mechanisms and biological behavior when studying PKC activation.

Study of PEGylated lipid layers as a model for PEGylated liposome surfaces: molecular dynamics simulation and Langmuir monolayer studies.

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design.

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.