Computer simulations on the pH-sensitive tri-block copolymer containing zwitterionic sulfobetaine as a novel anti-cancer drug carrier

pH-induced morphological transition of aggregates formed by miktoarm star polymers in dilute solution: a mesoscopic simulation study

The self-assembly of miktoarm star polymers μ-A i (B(D)) j C k in a neutral solution and the pH-responsive behaviors of vesicles and spherical micelles in an acidic solution have been investigated by DPD simulation. The results show that the self-assembled morphologies can be regulated by the lengths of pH-responsive arm B and hydrophilic arm C, leading to the formation of vesicles, discoidal micelles, and spherical micelles in a neutral solution. The dynamic evolution pathways of vesicles and spherical micelles are categorized into three stages: nucleation, coalescence, and growth. Subsequently, the pH-responsive behaviors of vesicles and spherical micelles have been explored by tuning the protonation degree of pH-responsive arm B. The vesicles evolves from nanodisks to nanosheets, then to nanoribbons, as the protonation degree increases, corresponding to a decrease in pH value, while the spherical micelles undergoes a transition into worm-like micelles, nanosheets, and nanoribbons. Notably, the electrostatic interaction leads the counterions to form a regular hexagonal pattern in nanosheets, while an alternative distribution of charged beads has been observed in nanoribbons. Furthermore, the role of the electrostatic interaction in the morphological transition has been elucidated through the analysis of the distribution of positive and negative charges, as well as the electrostatic potential for associates.

Investigation of morphology and structure of drug-loaded PLA-b-PEO-b-PLA polymeric micelle: A dissipative particle dynamics simulations study

The dissipative particle dynamics (DPD) simulation was used to study the morphologies and structures of the paclitaxel-loaded PLA-b-PEO-b-PLA polymeric micelle. We focused on the influences of PLA block length, PLA-b-PEO-b-PLA copolymer concentration, paclitaxel drug content on morphologies and structures of the micelle. Our simulations show that: (i) with the PLA block length increase, the self-assemble structure of PLA-b-PEO-b-PLA copolymers with paclitaxel vary between onion-like structure (core-middle layer-shell) to spherical core-shell structure. The PEO shell thins and the size of the PLA core increases. The onionlike structures are comprised of the PEO hydrophilic core, the PLA hydrophobic middle layer, and the PEO hydrophilic shell, the distribution of the paclitaxel drug predominantly occurs within the hydrophobic intermediate layer; (ii) The system forms a spherical core-shell structure when a small amount of the drug is added, and within a certain range, the size of the spherical structure increases as the drug amount increases. When the drug contents (volume fraction) cdrug = 10%, it can be observed that the PLA4-b-PEO19-b-PLA4 spherical structures connect to form rod-shaped structures. With the length of PLA block NPLA = 8, as the paclitaxel drug concentrations cdrug = 4%, PEO has been insufficient to completely encapsulate the PLA and paclitaxel drug beads. To enhance drug loading capacity while maintaining stability of the system in aqueous solution, the optimal composition for loading paclitaxel is PLA4-b-PEO19-b-PLA4; the drug content is not higher than 4%; (iii) The paclitaxel-loaded PLA4-b-PEO19-b-PLA4 micelle undergo the transition from onionlike (core-middle layer-shell) to spherical (core-shell) to rod-shaped and lamellar structure as the PLA4-b-PEO19-b-PLA4 copolymer concentration increases from ccp = 10% to 40%.

A dissipative particle dynamics simulation of controlled loading and responsive release of theranostic agents from reversible crosslinked triblock copolymer vesicles

To control the transport stability and release efficiency of loaded theranostic drugs in triblock copolymer carriers, the reversible crosslinking ability is of great significance. A molecular level exploration of such a function is needed to extend existing stabilizing and responsive dissociation mechanisms of carriers. Here, dissipative particle dynamics simulations were used to first demonstrate the formation of triblock copolymer vesicular carriers. Chemical crosslinking was used to strengthen the structural stability of the vesicle shell to avoid drug leakage. Reversible decrosslinking along with dissociation of the vesicle and release of loaded drugs were then explored. The structural, energetic and dynamical properties of the system were discussed at the molecular level. The regulation mechanism of drug release patterns was revealed by systematically exploring the effect of intra and intermolecular repulsive interactions. The results indicate that the chemical crosslinking of copolymers enhanced the compactness of the vesicle shell with a strengthened microstructure, increased binding energy, and limited chain migration, thus achieving more stable delivery of drugs. In terms of drug release, we clarified how the pairwise interactions of beads in the solution system affect the responsive dissociation of the vesicle and associated release patterns (speed and amount) of drugs. More efficient delivery and smart release of theranostic drugs are achieved using such reversible crosslinked triblock copolymer vesicles.

Modeling Polyzwitterion-Based Drug Delivery Platforms: A Perspective of the Current State-of-the-Art and Beyond

Drug delivery platforms are anticipated to have biocompatible and bioinert surfaces. PEGylation of drug carriers is the most approved method since it improves water solubility and colloid stability and decreases the drug vehicles' interactions with blood components. Although this approach extends their biocompatibility, biorecognition mechanisms prevent them from biodistribution and thus efficient drug transfer. Recent studies have shown (poly)zwitterions to be alternatives for PEG with superior biocompatibility. (Poly)zwitterions are super hydrophilic, mainly stimuli-responsive, easy to functionalize and they display an extremely low protein adsorption and long biodistribution time. These unique characteristics make them already promising candidates as drug delivery carriers. Furthermore, since they have highly dense charged groups with opposite signs, (poly)zwitterions are intensely hydrated under physiological conditions. This exceptional hydration potential makes them ideal for the design of therapeutic vehicles with antifouling capability, i.e., preventing undesired sorption of biologics from the human body in the drug delivery vehicle. Therefore, (poly)zwitterionic materials have been broadly applied in stimuli-responsive "intelligent" drug delivery systems as well as tumor-targeting carriers because of their excellent biocompatibility, low cytotoxicity, insignificant immunogenicity, high stability, and long circulation time. To tailor (poly)zwitterionic drug vehicles, an interpretation of the structural and stimuli-responsive behavior of this type of polymer is essential. To this end, a direct study of molecular-level interactions, orientations, configurations, and physicochemical properties of (poly)zwitterions is required, which can be achieved via molecular modeling, which has become an influential tool for discovering new materials and understanding diverse material phenomena. As the essential bridge between science and engineering, molecular simulations enable the fundamental understanding of the encapsulation and release behavior of intelligent drug-loaded (poly)zwitterion nanoparticles and can help us to systematically design their next generations. When combined with experiments, modeling can make quantitative predictions. This perspective article aims to illustrate key recent developments in (poly)zwitterion-based drug delivery systems. We summarize how to use predictive multiscale molecular modeling techniques to successfully boost the development of intelligent multifunctional (poly)zwitterions-based systems.

DPD Modelling of the Self- and Co-Assembly of Polymers and Polyelectrolytes in Aqueous Media: Impact on Polymer Science

This review article is addressed to a broad community of polymer scientists. We outline and analyse the fundamentals of the dissipative particle dynamics (DPD) simulation method from the point of view of polymer physics and review the articles on polymer systems published in approximately the last two decades, focusing on their impact on macromolecular science. Special attention is devoted to polymer and polyelectrolyte self- and co-assembly and self-organisation and to the problems connected with the implementation of explicit electrostatics in DPD numerical machinery. Critical analysis of the results of a number of successful DPD studies of complex polymer systems published recently documents the importance and suitability of this coarse-grained method for studying polymer systems.

Temperature-responsive morphology formation of a PS-b-PI copolymer: a dissipative particle dynamics simulation study

Self-assembly responsiveness to stimuli of polystyrene-block-polyisoprene (PS-b-PI) diblock copolymer materials is explored by means of classical molecular dynamics (MD) and dissipative particle dynamics (DPD) simulations. A concerted relationship between the parameters achieved from atomistic and DPD simulations is obtained for this molecular recognition as clearly pronounced in a phase transition. Effects of temperature, model size and composition on the morphological formation were systematically investigated for the diblock copolymeric system. Structural changes resulting in the evolution of rheology as well as an equilibrium ordered structure were analyzed in terms of order parameters and radial distribution functions. From our models, various morphologies were observed including discrete clusters (sphere-liked morphology), connected clusters (gyroid-liked morphology), hexagonally packed cylinders (HEX), connected cylinders, irregular cylinders, perfect lamellae, perforated lamellae and defected lamellae. Based on this finding, a bottom-up multi-scale simulation of the PS-b-PI diblock copolymer provides a link between equilibrium copolymeric morphologies and the crucial parameters.

A novel anticaries agent, honokiol-loaded poly(amido amine) dendrimer, for simultaneous long-term antibacterial treatment and remineralization of demineralized enamel

OBJECTIVE

Existing agents to induce enamel self-repair and inhibit the progression of dental caries in the early stage have been proven to be inadequate and far from satisfactory. In this study, a honokiol-loaded poly(amido amine) (PAMAM) dendrimer (PAMH) was constructed to combat early caries lesions in enamel.

METHODS

PAMH was prepared via a codissolution method. Computational simulation analysis was used to explore the mechanism of honokiol release. The cytotoxicity of PAMH was tested. The antibacterial effects of PAMH were tested by planktonic growth assays and biofilm formation inhibition assays. The remineralization effect of PAMH was examined via transverse microradiography and scanning electron microscopy after a pH cycling model. The in vivo anti-caries effect of PAMH was carried out in a rat model.

RESULTS

Honokiol released from PAMH was slower but more durable in a cariogenic pH environment than in a neutral pH environment, which could be explained through the computational simulation analysis results. Under electrostatic action, P3 beads with the same charge repelled each other and extended outwards, resulting in the rapid expansion of the PAMAM dendrimer and accelerating the release of the drug. At a low pH of 5.5, the protonated P3 beads were not charged and the protonated P1 beads were positively charged. However, the electrostatic repulsive interaction between protonated P1 beads was restricted by the P3 beads in the outermost layer of the PAMAM dendrimer, so the swelling rate was relatively slow, resulting in the slow release of drug molecules in the acidic environment. The cytotoxicity demonstration and the biocompatibility experiment in animal study showed that PAMH is biologically safe. PAMH showed excellent enamel remineralizing ability after pH cycling and showed a long-term antibacterial effect in vitro. Meanwhile, PAMH showed long-term anticaries efficacy in vivo.

SIGNIFICANCE

Our findings indicated that PAMH had great potential to combat early caries lesions in enamel for future clinical application.

Computer Simulations on a pH-Responsive Anticancer Drug Delivery System Using Zwitterion-Grafted Polyamidoamine Dendrimer Unimolecular Micelles

Unimolecular micelles have attracted wide attention in the field of drug delivery because of their thermodynamic stability and uniform size distribution. However, their drug loading/release mechanisms at the molecular level have been poorly understood. In this work, the stability and drug loading/release behaviors of unimolecular micelles formed using generation-5 polyamidoamine-graft-poly(carboxybetaine methacrylate) (PAMAM(G5)-PCBMA) were studied by dissipative particle dynamics simulations. In addition, the unimolecular micelles formed using generation-5 polyamidoamine-graft-poly(ethyleneglycol methacrylate) (PAMAM(G5)-PEGMA) were used as a comparison. The simulation results showed that PAMAM(G5)-PCBMA can spontaneously form core-shell unimolecular micelles. The PAMAM(G5) dendrimer constitutes a hydrophobic core to load the doxorubicin (DOX), while the zwitterionic PCBMA serves as a protective shell to improve the stability of the unimolecular micelle. The DOX can be encapsulated into the cavity of PAMAM(G5) at the physiological pH 7.4. The drug loading efficiency and drug loading content showed some regularities with the increase in the drug concentration. At the acidic pH 5.0, the loaded DOX can be released gradually from the hydrophobic core. The comparison of DOX-loaded morphologies between the PAMAM(G5)-PCBMA system and PAMAM(G5)-PEGMA system showed that the former has better monodisperse stability. This work could offer theoretical guidance for the design and development of promising unimolecular micelles for drug delivery.

Computer simulation of zwitterionic polymer brush grafted silica nanoparticles to modify polyvinylidene fluoride membrane

Dissipative particle dynamics (DPD) simulations was adopted to investigate the modification of polyvinylidene fluoride (PVDF) membrane by adding zwitterionic polymer brush poly(sulfobetaine methacrylate)- tetraethyl orthosilicate (PSBMA-TEOS) grafted silicon nanoparticles (SNPs) to the casting solution. The effects of polymer concentration and grafting architecture (PSBMA length and SNPs grafting ratio) on membrane morphology are discussed. When the polymer concentration reaches 40%, part of the SNPs is embedded in the membrane; the optimal polymer concentration is around 25-30%. In the SNPs system with the grafting ratio of 1, some SNPs are eluted into solution during phase separation. Compared with different grafting architectures, M_8-5, M_10-5 and M_12-5 system (M_x-y, where x represents the length of the zwitterionic polymer brush and y represents the grafting ratio of the silica nanoparticles) exhibited stable membrane morphologies. This work can provide guidance for the design and modification of organic-inorganic composite membrane and help understand the distribution of modified materials on the membrane surface.

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.

Biomimetic Design of Peptide Neutralizer of Ebola Virus with Molecular Simulation

Ebola virus (EBOV) belongs to the Filoviridae family, which can cause severe hemorrhagic fever in humans and nonprimates. The neutralization of EBOV by monoclonal antibody (mAb) ADI-15946 was reported recently. In the present study, the molecular interactions between the receptor GPcl of EBOV and ADI-15946 were studied by molecular dynamics (MD) simulation and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis. Hydrophobic interaction was identified as the main driving force for the binding of ADI-15946 on EBOV. Moreover, the contribution of each amino acid residue for the binding was evaluated. Then, an affinity binding model (ABM) was constructed using the residues favorable for the binding, including Y107, F108, D109, W110, and R113. The biomimetic design of neutralizer against EBOV according to the ABM of ADI-15946 was then performed, followed by screening using docking, structural similarity. Two neutralizers YFDWHMR and YFDWRYR were obtained, which were proven to be capable of strong binding on GPcl and then neutralizing GPcl. These results would be helpful for the development of neutralizers for Ebola virus.

pH-Responsive Zwitterionic Copolymer DHA-PBLG-PCB for Targeted Drug Delivery: A Computer Simulation Study

In this work, the self-assembled behaviors of zwitterionic copolymer docosahexaenoic acid- b-poly(γ-benzyl-l-glutamate)- b-poly(carboxybetaine methacrylate) (DHA-PBLG-PCB) and the loading and release mechanism of the anticancer drug doxorubicin (DOX) was investigated via computer simulations. The effects of polymer concentration, drug content, and pH on polymeric micelles were explored by dissipative particle dynamics (DPD) simulations. Simulation results show that DHA-PBLG₁₅-PCB₁₀ can self-assemble into core-shell micelles; in addition, the drug-loaded micelles have a pH-responsive feature. DOX can be encapsulated into the core-shell micelle under normal physiological pH conditions, whereas it can be released under acidic pH conditions. The self-assembled behaviors of copolymer DHA-PBLG-PEG were also studied to have a comparison with those of DHA-PBLG-PCB. The DHA-PBLG₁₅-PCB₁₀ system has a stable structure and it has a great potential to serve as drug delivery vehicles for targeted drug delivery.

4/6-Herto-arm and 4/6-mikto-arm star-shaped block polymeric drug-loaded micelles and their pH-responsive controlled release properties: a dissipative particle dynamics simulation

Star-shaped polymers have received significant attention and have been widely developed for prospective applications in drug delivery owing to their topological structure and unique physiochemical characteristics.

Computational modelling of drug delivery to solid tumour: Understanding the interplay between chemotherapeutics and biological system for optimised delivery systems

Drug delivery to solid tumour involves multiple physiological, biochemical and biophysical processes taking place across a wide range of length and time scales. The therapeutic efficacy of anticancer drugs is influenced by the complex interplays among the intrinsic properties of tumours, biophysical aspects of drug transport and cellular uptake. Mathematical and computational modelling allows for a well-controlled study on the individual and combined effects of a wide range of parameters on drug transport and therapeutic efficacy, which would not be possible or economically viable through experimental means. A wide spectrum of mathematical models has been developed for the simulation of drug transport and delivery in solid tumours, including PK/PD-based compartmental models, microscopic and macroscopic transport models, and molecular dynamics drug loading and release models. These models have been used as a tool to identify the limiting factors and for optimal design of efficient drug delivery systems. This article gives an overview of the currently available computational models for drug transport in solid tumours, together with their applications to novel drug delivery systems, such as nanoparticle-mediated drug delivery and convection-enhanced delivery.

Exploring Poly(ethylene glycol)-Polyzwitterion Diblock Copolymers as Biocompatible Smart Macrosurfactants Featuring UCST-Phase Behavior in Normal Saline Solution

Nonionic-zwitterionic diblock copolymers are designed to feature a coil-to-globule collapse transition with an upper critical solution temperature (UCST) in aqueous media, including physiological saline solution. The block copolymers that combine presumably highly biocompatible blocks are synthesized by chain extension of a poly(ethylene glycol) (PEG) macroinitiator via atom transfer radical polymerization (ATRP) of sulfobetaine and sulfabetaine methacrylates. Their thermoresponsive behavior is studied by variable temperature turbidimetry and ¹H NMR spectroscopy. While the polymers with polysulfobetaine blocks exhibit phase transitions in the physiologically interesting window of 30⁻50 °C only in pure aqueous solution, the polymers bearing polysulfabetaine blocks enabled phase transitions only in physiological saline solution. By copolymerizing a pair of structurally closely related sulfo- and sulfabetaine monomers, thermoresponsive behavior can be implemented in aqueous solutions of both low and high salinity. Surprisingly, the presence of the PEG blocks can affect the UCST-transitions of the polyzwitterions notably. In specific cases, this results in "schizophrenic" thermoresponsive behavior displaying simultaneously an UCST and an LCST (lower critical solution temperature) transition. Exploratory experiments on the UCST-transition triggered the encapsulation and release of various solvatochromic fluorescent dyes as model "cargos" failed, apparently due to the poor affinity even of charged organic compounds to the collapsed state of the polyzwitterions.

Hierarchically porous metal–organic frameworks with single-crystal structures and their enhanced catalytic properties

Stable hierarchically porous metal–organic frameworks (HP-MOFs) have been successfully synthesized under hydrothermal conditions using a template strategy.

Understanding the Cellular Uptake of pH-Responsive Zwitterionic Gold Nanoparticles: A Computer Simulation Study

Surface functionalization of nanoparticles (NPs) with stealth polymers (e.g., hydrophilic and zwitterionic polymers) has become a common strategy to resist nonspecific protein adsorption recently. Understanding the role of surface decoration on NP-biomembrane interactions is of great significance to promote the application of NPs in biomedical fields. Herein, using coarse-grained molecular dynamics (CGMD) simulations, we investigate the interactions between stealth polymer-coated gold nanoparticles (AuNPs) and lipid membranes. The results show that AuNPs grafted with zwitterionic polymers can more easily approach the membrane surface than those coated with hydrophilic poly(ethylene glycol) (PEG), which can be explained by the weak dipole-dipole interaction between them. For zwitterionic AuNPs which can undergo pH-dependent charge conversion, different interaction modes which depend on the polymer protonation degree are found. When the protonation degree is low, the particles just adsorb on the membrane surface; at moderate protonation degrees, the particles can directly translocate across the lipid membrane through a transient hydrophilic pore formed on the membrane surface; the particles are fully wrapped by the curved lipid membrane at high protonation degrees, which may lead to endocytosis. Finally, the effect of polymer chain length on the cellular uptake of zwitterionic polymer-coated AuNPs is considered. The results demonstrate that longer polymer chain length will block the translocation of AuNPs across the lipid membrane when the protonation degree is not high; however, it can improve the transmembrane efficiency of AuNPs at high protonation degrees. We expect that these findings are of immediate interest to the design and synthesis of pH-responsive nanomaterials based on zwitterionic polymers and can prompt their further applications in the field of biomedicine.

Mesoscopic Structures of Poly(carboxybetaine) Block Copolymer and Poly(ethylene glycol) Block Copolymer in Solutions

The antifouling property of exogenous materials is vital for their in vivo applications. In this work, dissipative particle dynamics simulations are performed to study the self-assembled morphologies of two copolymer systems containing poly(ethylene glycol) (PEG) and poly(carboxybetaine) (PCB) in aqueous solutions. Effects of polymer composition and polymer concentration on the self-assembled structures of the two copolymers (PLA-PEG and PLA-PCB) are investigated, respectively [PLA represents poly(lactic acid)]. Results show that whatever the copolymer composition is, PLA-PEG systems will self-assemble into core-shell structures, whereas onion-like and vesicle structures are also found for the PLA-PCB systems. Different morphologies are obtained at different polymer concentrations in both copolymer systems. Simulation results demonstrate that PCB is more stable than PEG in maintaining self-assembled spherical structures of copolymer systems because PLA-PEG forms dumbbell-like structures whereas PLA-PCB is spherical under the same polymer concentration. Although both copolymer systems can self-assemble into core-shell nanoparticles when the block ratio of PLA:PEG or PLA:PCB is 80:20, the core-shell structures of the nanoparticles are quite different. The shell layers formed by PEG in PLA-PEG nanoparticles are inhomogeneous in size because of the amphiphilicity of PEG, whereas the shell layers in PLA-PCB nanoparticles are homogenous because of the strong hydrophilicity of the zwitterionic PCB polymer block.