Molecular Modeling of PEGylated Peptides, Dendrimers, and Single-Walled Carbon Nanotubes for Biomedical Applications

Molecular Dynamics Simulations for Rationalizing Polymer Bioconjugation Strategies: Challenges, Recent Developments, and Future Opportunities

The covalent modification of proteins with polymers is a well-established method for improving the pharmacokinetic properties of therapeutically valuable biologics. The conjugated polymer chains of the resulting hybrid represent highly flexible macromolecular structures. As the dynamics of such systems remain rather elusive for established experimental techniques from the field of protein structure elucidation, molecular dynamics simulations have proven as a valuable tool for studying such conjugates at an atomistic level, thereby complementing experimental studies. With a focus on new developments, this review aims to provide researchers from the polymer bioconjugation field with a concise and up to date overview of such approaches. After introducing basic principles of molecular dynamics simulations, as well as methods for and potential pitfalls in modeling bioconjugates, the review illustrates how these computational techniques have contributed to the understanding of bioconjugates and bioconjugation strategies in the recent past and how they may lead to a more rational design of novel bioconjugates in the future.

Charge shielding effects of PEG bound to NH2-terminated PAMAM dendrimers - an experimental approach

Cationic poly(amido amine) (PAMAM) dendrimers exhibit great potential for use in drug delivery, but their high charge density leads to an inherent cytotoxicity. To increase biocompatibility, many studies have attached poly(ethylene glycol) (PEG) chains to the dendrimer surface. It is unclear how these tethered PEG chains influence the physicochemical properties of the dendrimer. Here, we develop a fluorescence-based assay utilizing anionic biological tissue to quantify the electrostatic binding affinity of a library of PEG-PAMAM conjugates with various PEG chain lengths and grafting densities. We find that covalently bound PEG chains reduce the electrostatic binding affinity more significantly than what can be achieved through covalent bonds only. Contrary to previous thought, this reduction is not explained by the steric hindrance effects of PEG chains, suggesting that other, non-covalent interactions between PEG and PAMAM are present. Using acetylated PAMAM conjugates, we convert electrostatic binding affinity to the number of charged amines accessible to the physiological environment. These data, coupled with ¹H-NMR, allows us to study more closely the non-covalent interactions between PEG and PAMAM. We find that increasing PEG chain length increases the number of non-covalent interactions. Additionally, at low grafting densities, increasing the number of PEG chains on the PAMAM surface also increases the non-covalent interactions. At higher grafting densities, however, PEG chains sterically repel one another, forcing chains to elongate away from the surface and reducing the number of interactions between PAMAM and individual PEG chains. The data presented here provides a framework for a more precise mechanistic understanding of how the length and density of tethered PEG chains on PAMAM dendrimers influence drug delivery properties.

Coarse-Grained Molecular Dynamics Simulations of Paclitaxel-Loaded Polymeric Micelles

A continuous manufacturing technology based on coaxial turbulent jet in coflow was previously developed to produce paclitaxel-loaded polymeric micelles. Herein, coarse-grained molecular dynamics (CG-MD) simulations were implemented to better understand the effect of the material attributes (i.e., the drug-polymer ratio and the ethanol concentration) and process parameters (i.e., temperature) on the self-assembly process of polymeric micelles as well as to provide molecular details on micelle instability. An all-atom (AA) poly (ethylene glycol)-poly (lactic acid) (PEG-PLA) polymer model was developed as the reference for parameterizing a coarse-grained (CG) model, and the AA polymer model was further validated with experimental glass transition temperature (T_g). The model transferability was verified by comparing structural properties between the AA and CG models. The CG model was further validated with experimental data, including micelle particle size measurements and drug encapsulation efficiency. Furthermore, the encapsulation of paclitaxel into the polymeric micelles was included in the simulations, taking into consideration the interactions between the paclitaxel and the polymers. The results from various points of view demonstrated a strong dependence of the shape of the micelles on the drug encapsulation, with micelles transitioning from spherical to ellipsoidal structures with an increasing paclitaxel amount. Simulation data were also used to identify the critical aggregation number (i.e., the number of polymer and drug molecules required for transition from one shape to another). Improved micellar structural stability was found with a larger micellar size and less solvent accessibility. Lastly, an evaluation was performed on the micellar dissociation free energy using a steered molecular dynamics simulation over a range of temperatures and ethanol concentrations. These simulations revealed that at higher ethanol and temperature conditions, micelles become destabilized, resulting in greater paclitaxel release. The increased drug release was determined to originate from the solvation of the hydrophobic core, which promoted micellar swelling and an associated reduction in hydrophobic interactions, leading to a loosely packed micellar structure.

Structure and dynamics of an active polymer adsorbed on the surface of a cylinder

The structure and dynamics of an active polymer on a smooth cylindrical surface are studied by Brownian dynamics simulations. The effect of an active force on the polymer adsorption behavior and the combined effect of chain mobility, length N, rigidity κ, and cylinder radius, R, on the phase diagrams are systemically investigated. We find that complete adsorption is replaced by the irregular alternative adsorption/desorption process at a large driving force. Three typical (spiral, helix-like, and rod-like) conformations of the active polymer are observed, dependent on N, κ, and R. Dynamically, the polymer shows rotational motion in the spiral state, snake-like motion in the intermediate state, and straight translational motion without turning back in the rod-like state. In the spiral state, we find that the rotation velocity ω and the chain length follow a power-law relation ω ∼ N^-0.42, consistent with the torque-balance theory of general Archimedean spirals. And the polymer shows super-diffusive behavior along the cylinder for a long time in the helix-like and rod-like states. Our results highlight that the mobility, rigidity, and curvature of surface can be used to regulate the polymer behavior.

PEGylation enhances the antibacterial and therapeutic potential of amphibian host defence peptides

Aurein 2.1, aurein 2.6 and aurein 3.1 are amphibian host defence peptides that kill bacteria via the use of lytic amphiphilic α-helical structures. The C-terminal PEGylation of these peptides led to decreased antibacterial activity (Minimum Lethal Concentration (MLCs) ↓ circa one and a half to threefold), reduced levels of amphiphilic α-helical structure in solvents (α-helicity ↓ circa 15.0%) and lower surface activity (Δπ ↓ > 1.5 mN m^-1). This PEGylation of aureins also led to decreased levels of amphiphilic α-helical structure in the presence of anionic membranes and zwitterionic membranes (α-helicity↓ > 10.0%) as well as reduced levels of penetration (Δπ ↓ > 3.0 mN m^-1) and lysis (lysis ↓ > 10.0%) of these membranes. Based on these data, it was proposed that the antibacterial action of PEGylated aureins involved the adoption of α-helical structures that promote the lysis of bacterial membranes, but with lower efficacy than their native counterparts. However, PEGylation also reduced the haemolytic activity of native aureins to negligible levels (haemolysis ↓ from circa 10% to 3% or less) and improved their relative therapeutic indices (RTIs ↑ circa three to sixfold). Based on these data, it is proposed that PEGylated aureins possess the potential for therapeutic development; for example, to combat infections due to multi-drug resistant strains of S. aureus, designated as high priority by the World Health Organization.

Analysis of Atomistic Potentials for Poly(ethylene glycol) Ethers

Two different potentials, the modified TraPPE-UA model of Fischer, J. [ J. Phys. Chem. B 2008, 112, 2388-2398] and the modified general AMBER force-field (GAFF) model of Barbosa, N. S. V. [ J. Mol. Model. 2017, 23, 194]are tested for a variety of temperature-dependent properties for neat and aqueous solutions for small poly(oxyethylene) (PEO) oligomers and larger polymers. A set of charges for PEO of arbitrary size is adapted for the modified GAFF model. Both models accurately reproduce experimental properties, but we find that the modified TraPPE-UA model is more accurate.

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.

Molecular Simulations of PEGylated Biomolecules, Liposomes, and Nanoparticles for Drug Delivery Applications

Since the first polyethylene glycol (PEG)ylated protein was approved by the FDA in 1990, PEGylation has been successfully applied to develop drug delivery systems through experiments, but these experimental results are not always easy to interpret at the atomic level because of the limited resolution of experimental techniques. To determine the optimal size, structure, and density of PEG for drug delivery, the structure and dynamics of PEGylated drug carriers need to be understood close to the atomic scale, as can be done using molecular dynamics simulations, assuming that these simulations can be validated by successful comparisons to experiments. Starting with the development of all-atom and coarse-grained PEG models in 1990s, PEGylated drug carriers have been widely simulated. In particular, recent advances in computer performance and simulation methodologies have allowed for molecular simulations of large complexes of PEGylated drug carriers interacting with other molecules such as anticancer drugs, plasma proteins, membranes, and receptors, which makes it possible to interpret experimental observations at a nearly atomistic resolution, as well as help in the rational design of drug delivery systems for applications in nanomedicine. Here, simulation studies on the following PEGylated drug topics will be reviewed: proteins and peptides, liposomes, and nanoparticles such as dendrimers and carbon nanotubes.

Structure-function relationships of nonviral gene vectors: Lessons from antimicrobial polymers

In recent years, substantial advances have been achieved in the design and synthesis of nonviral gene vectors. However, lack of effective and biocompatible vectors still remains a major challenge that hinders their application in clinical settings. In the past decade, there has been a rapid expansion of cationic antimicrobial polymers, due to their potent, rapid, and broad-spectrum biocidal activity against resistant microbes, and biocompatible features. Given that antimicrobial polymers share common features with nonviral gene vectors in various aspects, such as membrane affinity, functional groups, physicochemical characteristics, and unique macromolecular architectures, these polymers may provide us with inspirations to overcome challenges in the design of novel vectors toward more safe and efficient gene delivery in clinic. Building off these observations, we provide here an overview of the structure-function relationships of polymers for both antimicrobial applications and gene delivery by elaborating some key structural parameters, including functional groups, charge density, hydrophobic/hydrophilic balance, MW, and macromolecular architectures. By borrowing a leaf from antimicrobial agents, great advancement in the development of newer nonviral gene vectors with high transfection efficiency and biocompatibility will be more promising. STATEMENT OF SIGNIFICANCE: The development of gene delivery is still in the preclinical stage for the lack of effective and biocompatible vectors. Given that antimicrobial polymers share common features with gene vectors in various aspects, such as membrane affinity, functional groups, physicochemical characteristics, and unique macromolecular architectures, these polymers may provide us with inspirations to overcome challenges in the design of novel vectors toward more safe and efficient gene delivery in clinic. In this review, we systematically summarized the structure-function relationships of antimicrobial polymers and gene vectors, with which the design of more advanced nonviral gene vectors is anticipated to be further boosted in the future.

Molecular Self-Assembly Strategy for Encapsulation of an Amphipathic α-Helical Antimicrobial Peptide into the Different Polymeric and Copolymeric Nanoparticles

Encapsulation of peptide and protein-based drugs in polymeric nanoparticles is one of the fundamental fields in controlled-release drug delivery systems. The molecular mechanisms of absorption of peptides to the polymeric nanoparticles are still unknown, and there is no precise molecular data on the encapsulation process of peptide and protein-based drugs. Herein, the self-assembly of different polymers and block copolymers with combinations of the various molecular weight of blocks and the effects of resultant polymer and copolymer nanomicelles on the stability of magainin2, an α-helical antimicrobial peptide, were investigated by means of all-atom molecular dynamics (MD) simulation. The micelle forming, morphology of micellar aggregations and changes in the first hydration shell of the micelles during micelles formation were explored as well. The results showed that the peptide binds to the polymer and copolymer micelles and never detaches during the MD simulation time. In general, all polymers and copolymers simultaneously encapsulated the peptide during micelles formation and had the ability to maintain the helical structure of the peptide, whereas the first hydration shell of the peptide remained unchanged. Among the micelles, the polyethylene glycol (PEG) micelles completely encapsulated magainin2 and, surprisingly, the NMR structure of the peptide was perfectly kept during the encapsulation process. The MD results also indicated that the aromatic and basic residues of the peptide strongly interact with polymers/copolymers and play important roles in the encapsulation mechanism. This research will provide a good opportunity in the design of polymer surfaces for drug delivery applications such as controlled-release peptide delivery systems.

Noncovalent Pyrene-Polyethylene Glycol Coatings of Carbon Nanotubes Achieve in Vitro Biocompatibility

Single-walled carbon nanotubes (SWNTs) have become increasingly exploited in biological applications, such as imaging and drug delivery. The application of SWNTs in biological settings requires the surface chemistry to remain through the low solubility in aqueous media. In this research, a facile approach for the preparation of a polyethylene glycol (PEG)-coated SWNT-based nanocarrier was reported. We focused on the effect of PEG chain length and SWNT size on the cytotoxicity of PEG-coated SWNTs as a superior drug delivery nanovector. First, all-atom molecular dynamics (MD) simulations were employed to explore the stability and behavior of SWNT/pyrene-PEG (SWNT/Pyr-PEG) structures at a molecular level that is not attainable with experiments. The MD studies revealed that (i) π-π stacking interactions between the pyrene bearing PEG molecules and SWNTs are maintained in bulky situations, regardless of PEG molecular weight or SWNT size; (ii) pyrene molecules diffuse over the SWNT surface without detaching; and (iii) both short and long dynamic Pyr-PEG chains have the capability of effectively coating the SWNT surface. In light of the simulations, noncovalent (π-π stacking) assemblies of SWNT/Pyr-PEG with different molecular weights of PEG ( M_w = 2000, 5000, and 12000) were successfully fabricated and characterized. For longer PEG chains, more effective coating of SWNTs was obtained, resulting in more biocompatible SWNT/Pyr-PEG nanomaterials. The number of SWNTs coated by Pyr-PEG was highly dependent on the length of pyrene bearing PEG polymers. Moreover, the short SWNTs showed a higher amount of PEG coating with respect to the long SWNTs. Cell viability results demonstrated a dose-dependent cytotoxicity of coated SWNTs. Short SWNTs coated with longer PEG chains have low cytotoxicity to be used in in vivo studies.

Coarse-Grained Simulations of Aqueous Thermoresponsive Polyethers

Thermoresponsive polymers can change structure or solubility as a function of temperature. Block co-polymers of polyethers have a response that depends on polymer molecular weight and co-polymer composition. A coarse-grained model for aqueous polyethers is developed and applied to polyethylene oxide and polyethylene oxide-polypropylene oxide-polyethylene oxide triblock co-polymers. In this model, no interaction sites on hydrogen atoms are included, no Coulombic interactions are present, and all interactions are short-ranged, treated with a combination of two- and three-body terms. Our simulations find that The triblock co-polymers tend to associate at temperatures above 350 K. The aggregation is stabilized by contact between The hydrophobic methyl groups on The propylene oxide monomers and involves a large, favorable change in entropy.

Solid Phase Stepwise Synthesis of Polyethylene Glycols

Polyethylene glycol (PEG) and derivatives with eight and twelve ethylene glycol units were synthesized by stepwise addition of tetraethylene glycol monomers on a polystyrene solid support. The monomer contains a tosyl group at one end and a dimethoxytrityl group at the other. The Wang resin, which contains the 4-benzyloxy benzyl alcohol function, was used as the support. The synthetic cycle consists of deprotonation, Williamson ether formation (coupling), and detritylation. Cleavage of PEGs from solid support was achieved with trifluoroacetic acid. The synthesis including monomer synthesis was entirely chromatography-free. PEG products including those with different functionalities at the two termini were obtained in high yields. The products were analyzed with ESI and MALDI-TOF MS and were found close to monodispersity.

In vitro blood cell viability profiling of polymers used in molecular assembly

Biocompatible polymers have been extensively applied to molecular assembly techniques on a micro- and nanoscale to miniaturize functional devices for biomedical uses. However, cytotoxic assessments of developed devices are prone to partially focus on non-specific cells or cells associated with the specific applications. Thereby, since toxicity is dependent on the type of cells and protocols, we do not fully understand the relative toxicities of polymers. Additionally, we need to ensure the blood cell biocompatibility of developed devices prior to that of targeted cells because most of the devices contact the blood before reaching the targeted regions. Motivated by this issue, we focused on screening cytotoxicity of polymers widely used for the layer-by-layer assembly technique using human blood cells. Cytotoxicity at the early stage was investigated on twenty types of polymers (positively charged, negatively charged, or neutral) and ten combination forms via hemolysis, cell viability, and AnnexinV-FITC/PI staining assays. We determined their effects on the cell membrane depending on their surface chemistry by molecular dynamics simulations. Furthermore, the toxicity of LbL-assembled nanofilms was assessed by measuring cell viability. Based on this report, researchers can produce nanofilms that are better suited for drug delivery and biomedical applications by reducing the possible cytotoxicity.

Novel Biodegradable Polymer with Redox-Triggered Backbone Cleavage Through Sequential 1,6-Elimination and 1,5-Cyclization Reactions

In the past decade, the self-immolative biodegradable polymer arose as a novel paradigm for its efficient degradation mechanism and vast potential for advanced biomedical applications. This study reports successful synthesis of a novel biodegradable polymer capable of self-immolative backbone cleavage. The monomer is designed by covalent conjugations of both pendant redox-trigger (p-nitrobenzyl alcohol) and self-immolative linker (p-hydroxybenzyl alcohol) to the cyclization spacer (n-2-(hydroxyethyl)ethylene diamine), which serves as the structural backbone. The polymerization of the monomer with hexamethylene diisocyanate yields a linear redox-sensitive polymer that can systemically degrade via sequential 1,6-elimination and 1,5-cyclization reactions within an effective timeframe. Ultimately, the polymer's potential for biomedical application is simulated through in vitro redox-triggered release of paclitaxel from polymeric nanoparticles.

Effect of PEG Architecture on the Hybridization Thermodynamics and Protein Accessibility of PEGylated Oligonucleotides

PEGylation is an attractive approach to modifying oligonucleotides intended for therapeutic purposes. PEG conjugation reduces protein interactions with the oligonucleotide, and helps to overcome their intrinsic biopharmaceutical shortcomings, such as poor enzymatic stability, rapid body clearance, and unwanted immunostimulation. However, the effect of PEG architecture and the manner in which the PEG component interferes with the hybridization of the oligonucleotide remain poorly understood. In this study, we systematically compare the hybridization thermodynamics and protein accessibility of several DNA conjugates involving linear, Y-shaped, and brush-type PEG. It is found that PEGylated DNA experiences two opposing effects: local excluded volume effect and chemical interactions, the strengths of which are architecture-dependent. Notably, the brush architecture is able to offer significantly greater protein shielding capacity than its linear or Y-shaped counterparts, while maintaining nearly identical free energy for DNA hybridization compared with free DNA.

Dynamical Interactions of 5-Fluorouracil Drug with Dendritic Peptide Vectors: The Impact of Dendrimer Generation, Charge, Counterions, and Structured Water

Molecular dynamics simulations are utilized to investigate the interactions between the skin cancer drug 5-fluorouracil (5FU) and peptide-based dendritic carrier systems. We find that these drug-carrier interactions do not conform to the traditional picture of long-time retention of the drug within a hydrophobic core of the dendrimer carrier. Rather, 5FU, which is moderately soluble in its own right, experiences weak, transient chattering interactions all over the dendrimer, mediated through multiple short-lived hydrogen bonding and close contact events. We find that charge on the periphery of the dendrimer actually has a negative effect on the frequency of drug-carrier interactions due to a counterion screening effect that has not previously been observed. However, charge is nevertheless an important feature since neutral dendrimers are shown to have a significant mutual attraction that can lead to clustering or agglomeration. This clustering is prevented due to charge repulsion for the titrated dendrimers, such that they remain independent in solution.

Computational and experimental approaches for investigating nanoparticle-based drug delivery systems

Most therapeutic agents suffer from poor solubility, rapid clearance from the blood stream, a lack of targeting, and often poor translocation ability across cell membranes. Drug/gene delivery systems (DDSs) are capable of overcoming some of these barriers to enhance delivery of drugs to their right place of action, e.g. inside cancer cells. In this review, we focus on nanoparticles as DDSs. Complementary experimental and computational studies have enhanced our understanding of the mechanism of action of nanocarriers and their underlying interactions with drugs, biomembranes and other biological molecules. We review key biophysical aspects of DDSs and discuss how computer modeling can assist in rational design of DDSs with improved and optimized properties. We summarize commonly used experimental techniques for the study of DDSs. Then we review computational studies for several major categories of nanocarriers, including dendrimers and dendrons, polymer-, peptide-, nucleic acid-, lipid-, and carbon-based DDSs, and gold nanoparticles. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

In vivo drug delivery of gemcitabine with PEGylated single-walled carbon nanotubes

Gemcitabine (GEM) is an anticancer agent widely used in non-small cell lung and pancreatic cancers. The clinical use of GEM has been limited by its rapid metabolism and short plasma half-life. These restrictions lead to frequent administration of high drug doses which can cause severe side effects. Therefore, new delivery strategies are needed aiming toward improved therapeutic effects. Single-walled carbon nanotubes (SWCNTs) are emerging as promising carriers for drug delivery due to their unique properties including high drug loading capacities, notable cell membrane penetrability and prolonged circulation times. In this work, pristine SWCNTs were functionalized through carboxylation, acylation, amination, PEGylation and finally GEM conjugation. The prepared SWCNT-GEM and SWCNT-PEG-GEM conjugates were characterized by FTIR, NMR, DSC and TEM to confirm the successful functionalization. The amount of GEM bound to the conjugates was 43.14% (w/w) for the SWCNT-GEM and 37.32% for the SWCNT-PEG-GEM, indicating high loading capacity. MTT assay on the human lung carcinoma cell line (A549) and the human pancreatic carcinoma cell line (MIA PaCa-2) demonstrated that the SWCNT-GEM was more cytotoxic than SWCNT-PEG-GEM and GEM. The SWCNT-PEG-GEM conjugates afford higher efficacy in suppressing tumor growth than SWCNT-GEM and GEM in B6 nude mice. The results demonstrate that the new formulation of GEM is useful strategy for improving the antitumor efficacy of GEM.