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.

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.

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.

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.

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.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery.

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

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.

Cholesterol Reduces Partitioning of Antifungal Drug Itraconazole into Lipid Bilayers.

Cholesterol plays a crucial role in modulating the physicochemical properties of biomembranes, both increasing mechanical strength and decreasing permeability. Cholesterol is also a common component of vesicle-based delivery systems, including liposome-based drug delivery systems (LDSs). However, its effect on the partitioning of drug molecules to lipid membranes is very poorly recognized. Herein, we performed a combined experimental/computational study of the potential for the use of the LDS formulation for the delivery of the antifungal drug itraconazole (ITZ). We consider the addition of cholesterol to the lipid membrane. Since ITZ is only weakly soluble in water, its bioavailability is limited. Use of an LDS has thus been proposed. We studied lipid membranes composed of cholesterol, 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphocholine (POPC), and ITZ using a combination of computational molecular dynamics (MD) simulations of lipid bilayers and Brewster angle microscopy (BAM) experiments of monolayers. Both experimental and computational results show separation of cholesterol and ITZ. Cholesterol has a strong preference to orient parallel to the bilayer normal. However, ITZ, a long and relatively rigid molecule with weakly hydrophilic groups along the backbone, predominantly locates below the interface between the hydrocarbon chain region and the polar region of the membrane, with its backbone oriented parallel to the membrane surface; the orthogonal orientation in the membrane could be the cause of the observed separation. In addition, fluorescence measurements demonstrated that the affinity of ITZ for the lipid membrane is decreased by the presence of cholesterol, which is thus probably not a suitable formulation component of an LDS designed for ITZ delivery.

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.

Insight into the antimicrobial mechanism of action of ß2,2-amino acid derivatives from molecular dynamics simulation: Dancing the can-can at the membrane surface.

The development of antimicrobial agents that target and selectively disrupt biofilms is a pressing issue since, so far, no antibiotics have been developed that achieve this effectively. Previous experimental work has found a promising set of antibacterial peptides: ß2,2-amino acid derivatives, relatively small molecules with common structural elements composed of a polar head group and two non-polar hydrocarbon arms. In order to develop insight into possible mechanisms of action of these novel antibacterial agents, we have performed an in silico investigation of four leading ß2,2-amino acid derivatives, interacting with models of both bacterial (target) and eukaryotic (host) membranes, using molecular dynamics simulation with a model with all-atom resolution. We found an unexpected result that could shed light on the mechanism of action of these antimicrobial agents: the molecules assume a conformation where one of the hydrophobic arms is directed downward into the membrane core while the other is directed upwards, out of the membrane and exposed above the position of the membrane headgroups; we dubbed this conformation the "can-can pose". Intriguingly, the can-can pose was most closely linked to the choice of headgroup. Also, the compound previously found to be most effective against biofilms displayed the strongest extent of this behavior and, additionally, this behavior was more pronounced for this compound in the bacterial than in the eukaryotic membrane. We hypothesize that adopting the can-can pose could possibly disrupt the protective peptidoglycan macronet found on the exterior of the bacterial membrane.

Determinants of Orexin Receptor Binding and Activation-A Molecular Dynamics Study.

We assess the stability of two previously suggested binding modes for the neuropeptide orexin-A in the OX2 receptor through extensive molecular dynamics simulations. As the activation determinants of the receptor remain unknown, we simulated an unliganded receptor and two small-molecular ligands, the antagonist suvorexant and the agonist Nag26 for comparison. Each system was simulated in pure POPC membrane as well as in the 25% cholesterol-POPC membrane. In total, we carried out 36 µs of simulations. Through this set of simulations, we report a stable binding mode for the C-terminus of orexin-A. In addition, we suggest interactions that would promote orexin receptor activation, as well as others that would stabilize the inactive state.

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 ᅟ.

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.

The effect of light sensitizer localization on the stability of indocyanine green liposomes.

Light triggered drug delivery systems offer attractive possibilities for sophisticated therapy, providing both temporal and spatial control of drug release. We have developed light triggered liposomes with clinically approved indocyanine green (ICG) as the light sensitizing compound. Amphiphilic ICG can be localized in different compartments of the liposomes, but the effect of its presence, on both triggered release and long term stability, has not been studied. In this work, we report that ICG localization has a significant effect on the properties of the liposomes. Polyethylene glycol (PEG) coating of the liposomes leads to binding and stabilization of the ICG molecules on the surface of the lipid bilayer. This formulation showed both good storage stability in buffer solution (at +4-37 °C) and adequate stability in serum and vitreous (at +37 °C). The combination of ICG within the lipid bilayer and PEG coating lead to poor stability at elevated temperatures of +22 °C and + 37 °C. The mechanisms of the increased instability due to ICG insertion in the lipid bilayer was elucidated with molecular dynamics simulations. Significant PEG insertion into the bilayer was induced in the presence of ICG in the lipid bilayer. Finally, feasibility of freeze-drying as a long term storage method for the ICG liposomes was demonstrated. Overall, this is the first detailed study on the interactions of lipid bilayer, light sensitizer (ICG) and PEG coating on the liposome stability. The localization of the light triggering agent significantly alters the structure of the liposomes and it is important to consider these aspects in triggered drug delivery system design.

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.

Comparison of time-gated surface-enhanced raman spectroscopy (TG-SERS) and classical SERS based monitoring of Escherichia coli cultivation samples.

The application of Raman spectroscopy as a monitoring technique for bioprocesses is severely limited by a large background signal originating from fluorescing compounds in the culture media. Here, we compare time-gated Raman (TG-Raman)-, continuous wave NIR-process Raman (NIR-Raman), and continuous wave micro-Raman (micro-Raman) approaches in combination with surface enhanced Raman spectroscopy (SERS) for their potential to overcome this limit. For that purpose, we monitored metabolite concentrations of Escherichia coli bioreactor cultivations in cell-free supernatant samples. We investigated concentration transients of glucose, acetate, AMP, and cAMP at alternating substrate availability, from deficiency to excess. Raman and SERS signals were compared to off-line metabolite analysis of carbohydrates, carboxylic acids, and nucleotides. Results demonstrate that SERS, in almost all cases, led to a higher number of identifiable signals and better resolved spectra. Spectra derived from the TG-Raman were comparable to those of micro-Raman resulting in well-discernable Raman peaks, which allowed for the identification of a higher number of compounds. In contrast, NIR-Raman provided a superior performance for the quantitative evaluation of analytes, both with and without SERS nanoparticles when using multivariate data analysis. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1533-1542, 2018.

Membrane bound COMT isoform is an interfacial enzyme: general mechanism and new drug design paradigm.

The enzyme catechol-O-methyltransferase (COMT) has water soluble (S-COMT) and membrane associated (MB-COMT), bitopic, isoforms. Of these MB-COMT is a drug target in relation to the treatment of Parkinson's disease. Using a combination of computational and experimental protocols, we have determined the substrate selection mechanism specific to MB-COMT. We show: (1) substrates with preferred affinity for MB-COMT over S-COMT orient in the membrane in a fashion conducive to catalysis from the membrane surface and (2) binding of COMT to its cofactor ADOMET induces conformational change that drives the catalytic surface of the protein to the membrane surface, where the substrates and Mg2+ ions, required for catalysis, are found. Bioinformatics analysis reveals evidence of this mechanism in other proteins, including several existing drug targets. The development of new COMT inhibitors with preferential affinity for MB-COMT over S-COMT is now possible and insight of broader relevance, into the function of bitopic enzymes, is provided.