Molecular Dynamics Simulations of PEI Mediated DNA Aggregation

Reversible stabilization of DNA/PEI complexes by reducible click-linkage between DNA and polymer. A new polyplex concept for lowering polymer quantity

Nonviral transfection of mammalian cells can be performed with electrostatic complexes (polyplexes) between a plasmid DNA (pDNA) encoding a foreign gene and a cationic polymer. However, an excess of the cationic polymer is required for pDNA condensation and polyplexes formation, which generate in vivo toxicity. Here, we present a new concept of polyplexes preparation aiming to reduce the polymer quantity. pDNA was functionalized with 3,6,9-trioxaundecan-1- {4 - [(2-chloroethyl) ethylamino)] - benzylamino}, 11-azide, and polyethyleneimine (lPEI) with reducible dibenzocyclooctyl (SS-DBCO) groups allowing azide-alkyne cycloaddition between pDNA and lPEI after condensation. The size of polyplexes with DBCO-SS-lPEI was smaller than with lPEI due to a stronger DNA condensation thanks to linkages between polymer and pDNA preventing dissociation until disulfide bridges reduction. In vitro transfection showed that the amount of DBCO-SS-lPEI leading to the most efficient polyplexes was three times lower than lPEI. As expected, toxicity in mice was significantly reduced upon intravenous injection of DBCO-SS-lPEI polyplexes at doses where the lPEI polyplexes killed mice. This is probably due to the high stability of the DBCO-SS-lPEI polyplexes which prevented their aggregation in the pulmonary capillaries. Overall, this new concept of polyplexes with DBCO-SS-lPEI offering the possibility of administering higher doses of polyplexes than lPEI and their ability to pass the pulmonary barrier could be favorably exploited for transfection of distant organs or tissues, such as tumors.

In silico study of PEI-PEG-squalene-dsDNA polyplex formation: the delicate role of the PEG length in the binding of PEI to DNA

Biocompatible hydrophilic polyethylene glycol (PEG) is widely used in biomedical applications, such as drug or gene delivery, tissue engineering or as an antifouling component in biomedical devices. Experimental studies have shown that the size of PEG can weaken polycation-polyanion interactions, like those between branched polyethyleneimine (b-PEI) and DNA in gene carriers, but details of its cause and underlying interactions on the atomic scale are still not clear. To better understand the interaction mechanisms in the formation of polyplexes between b-PEI-PEG based carriers and DNA, we have used a combination of in silico tools and experiments on three multicomponent systems differing in PEG MW. Using the PEI-PEG-squalene-dsDNA systems of the same size, both in the all-atom MD simulations and in experimental in-gel electrophoresis measurements, we found that the binding between DNA and the vectors is highly influenced by the size of PEG, with the binding efficiency increasing with a shorter PEG length. The mechanism of how PEG interferes with the binding between PEI and DNA is explained using a two-step MD simulation protocol that showed that the DNA-vector interactions are influenced by the PEG length due to the hydrogen bond formation between PEI and PEG. Although computationally demanding we find it important to study molecular systems of the same size both in silico and in a laboratory and to simulate the behaviour of the carrier prior to the addition of bioactive molecules to understand the molecular mechanisms involved in the formation of the polyplex.

Controlled DNA Delivery Using Poly(lactide) Nanoparticles and Understanding the Binding Interactions

Cationic polymer-based gene delivery vectors suffer from several limitations such as low DNA-loading capacity, poor transfection, toxicity, environmental degradations, etc. Again, very limited works are available addressing the binding interactions in detail at the atomic level explaining the loading capacity, protection ability against harsh environments, and controlled release behavior of the DNA-encapsulated vehicles. Here, a poly(l-lactide) (PLA) nanoparticle-based controlled DNA release system is proposed. The developed vehicle possesses a high DNA-loading capacity and can release the loaded DNA in a controlled manner. Spectroscopic, physicochemical, and molecular simulation techniques (AM1 and atomistic molecular dynamics) have been employed to understand the binding interactions between PLA and DNA molecules enabling high DNA loading, protection against external harsh environments, and controlled DNA release behavior. Methyl thiazolyl tetrazolium (MTT) assay experiments confirm the biocompatible nature of the vehicle. Cellular uptake efficiency and endo-lysosomal escape capabilities have been investigated against HeLA cells. This study, therefore, demonstrates the development of a promising nonviral DNA delivery vector and includes a detailed investigation of the atomic-level interaction behavior between PLA and DNA molecules.

Polymeric Delivery of Therapeutic Nucleic Acids

The advent of genome editing has transformed the therapeutic landscape for several debilitating diseases, and the clinical outlook for gene therapeutics has never been more promising. The therapeutic potential of nucleic acids has been limited by a reliance on engineered viral vectors for delivery. Chemically defined polymers can remediate technological, regulatory, and clinical challenges associated with viral modes of gene delivery. Because of their scalability, versatility, and exquisite tunability, polymers are ideal biomaterial platforms for delivering nucleic acid payloads efficiently while minimizing immune response and cellular toxicity. While polymeric gene delivery has progressed significantly in the past four decades, clinical translation of polymeric vehicles faces several formidable challenges. The aim of our Account is to illustrate diverse concepts in designing polymeric vectors towards meeting therapeutic goals of in vivo and ex vivo gene therapy. Here, we highlight several classes of polymers employed in gene delivery and summarize the recent work on understanding the contributions of chemical and architectural design parameters. We touch upon characterization methods used to visualize and understand events transpiring at the interfaces between polymer, nucleic acids, and the physiological environment. We conclude that interdisciplinary approaches and methodologies motivated by fundamental questions are key to designing high-performing polymeric vehicles for gene therapy.

Structure and Dynamics of Thermosensitive pDNA Polyplexes Studied by Time-Resolved Fluorescence Spectroscopy

Combining multiple stimuli-responsive functionalities into the polymer design is an attractive approach to improve nucleic acid delivery. However, more in-depth fundamental understanding how the multiple functionalities in the polymer structures are influencing polyplex formation and stability is essential for the rational development of such delivery systems. Therefore, in this study the structure and dynamics of thermosensitive polyplexes were investigated by tracking the behavior of labeled plasmid DNA (pDNA) and polymer with time-resolved fluorescence spectroscopy using fluorescence resonance energy transfer (FRET). The successful synthesis of a heterofunctional poly(ethylene glycol) (PEG) macroinitiator containing both an atom transfer radical polymerization (ATRP) and reversible addition-fragmentation chain-transfer (RAFT) initiator is reported. The use of this novel PEG macroinitiator allows for the controlled polymerization of cationic and thermosensitive linear triblock copolymers and labeling of the chain-end with a fluorescent dye by maleimide-thiol chemistry. The polymers consisted of a thermosensitive poly(N-isopropylacrylamide) (PNIPAM, N), hydrophilic PEG (P), and cationic poly(2-(dimethylamino)ethyl methacrylate) (PDMAEMA, D) block, further referred to as NPD. Polymer block D chain-ends were labeled with Cy3, while pDNA was labeled with FITC. The thermosensitive NPD polymers were used to prepare pDNA polyplexes, and the effect of the N/P charge ratio, temperature, and composition of the triblock copolymer on the polyplex properties were investigated, taking nonthermosensitive PD polymers as the control. FRET was observed both at 4 and 37 °C, indicating that the introduction of the thermosensitive PNIPAM block did not compromise the polyplex structure even above the polymer's cloud point. Furthermore, FRET results showed that the NPD- and PD-based polyplexes have a less dense core compared to polyplexes based on cationic homopolymers (such as PEI) as reported before. The polyplexes showed to have a dynamic character meaning that the polymer chains can exchange between the polyplex core and shell. Mobility of the polymers allow their uniform redistribution within the polyplex and this feature has been reported to be favorable in the context of pDNA release and subsequent improved transfection efficiency, compared to nondynamic formulations.

Collapse of DNA in packaging and cellular transport

The dawn of molecular biology and recombinant DNA technology arose from our ability to manipulate DNA, including the process of collapse of long extended DNA molecules into nanoparticles of approximately 100 nm diameter. This condensation process is important for the packaging of DNA in the cell and for transporting DNA through the cell membrane for gene therapy. Multivalent cations, such as natural polyamines (spermidine and spermine), were initially recognized for their ability to provoke DNA condensation. Current research is targeted on molecules such as linear and branched polymers, oligopeptides, polypeptides and dendrimers that promote collapse of DNA to nanometric particles for gene therapy and on the energetics of DNA packaging.

Interfacial analysis of siRNA complexes with poly-ethylenimine (PEI) or PAMAM dendrimers in gene delivery

Solution and interfacial analysis has been employed to gain insight into the complexation of siRNA using either G4 PAMAM dendrimers or 25kDa branched poly-ethylenimine (bPEI). The size, charge and shape/structure of the complexing agents were probed using atomic force microscopy (AFM), circular dichroism spectrometry (CD), dynamic light scattering (DLS), and gel electrophoresis (GE). The binding capability of these cationic polymers to the siRNA molecule, subsequently controls the surface/adsorption behaviour of the complexes to a negatively charged surface. G4 PAMAM dendrimers bind to the major groove of the siRNA structure, while bPEI binds to both major and minor groove. PAMAM-siRNA complexes form a thin uniform surface film with adsorption of monomeric particles, whilst bPEI-siRNA complexes adsorb as particles in random orientations at low bPEI concentration and form network structures across the surface at high charge ratio. This is due to their ability to bind to both regions within siRNA. This new understanding of the interfacial behaviour of siRNA complexes correlates with observations of cellular transfection and can be used in the design of optimal transfection agents.

Biodegradable Poly(Amino Ester) with Aromatic Backbone as Efficient Nonviral Gene Delivery Vectors

The development of gene delivery vectors with high efficiency and biocompatibility is one of the critical points of gene therapy. Two biodegradable poly(amino ester)s were synthesized via ring-opening polymerization between low molecular weight (LMW) PEI and diepoxide. The molecular weights of poly(amino ester)s were measured by GPC. Agarose gel retardation assays showed that these materials have good DNA-binding ability and can retard the electrophoretic mobility of plasmid DNA (pDNA) at a weight ratio of 1. The formed polyplexes have proper sizes of around 200 nm and zeta-potential values of about 30-40 mV for cellular uptake. In vitro experiments revealed that polymer P2 gave higher transfection efficiency than PEI 25KDa and Lipofectamine 2000 with less toxicity, especially in 293 cells. Results demonstrate that such a type of degradable poly(amino ester) may serve as a promising non-viral gene vector.

Polymer-Nucleic Acid Interactions

Gene therapy is an important therapeutic strategy in the treatment of a wide range of genetic disorders. Polymers forming stable complexes with nucleic acids (NAs) are non-viral gene carriers. The self-assembly of polymers and nucleic acids is typically a complex process that involves many types of interaction at different scales. Electrostatic interaction, hydrophobic interaction, and hydrogen bonds are three important and prevalent interactions in the polymer/nucleic acid system. Electrostatic interactions and hydrogen bonds are the main driving forces for the condensation of nucleic acids, while hydrophobic interactions play a significant role in the cellular uptake and endosomal escape of polymer-nucleic acid complexes. To design high-efficiency polymer candidates for the DNA and siRNA delivery, it is necessary to have a detailed understanding of the interactions between them in solution. In this chapter, we survey the roles of the three important interactions between polymers and nucleic acids during the formation of polyplexes and summarize recent understandings of the linear polyelectrolyte-NA interactions and dendrimer-NA interactions. We also review recent progress optimizing the gene delivery system by tuning these interactions.

Structural Comparisons of PEI/DNA and PEI/siRNA Complexes Revealed with Molecular Dynamics Simulations

Polyplexes composed of polyethyleneimine (PEI) and DNA or siRNA have attracted great attention for their use in gene therapy. Although many physicochemical characteristics of these polyplexes remain unknown, PEI/DNA complexes have been repeatedly shown to be more stable than their PEI/siRNA counterparts. Here, we examine potential causes for this difference using atomistic molecular dynamics simulations of complexation between linear PEI and DNA or siRNA duplexes containing the same number of bases. The two types of polyplexes are stabilized by similar interactions, as PEI amines primarily interact with nucleic acid phosphate groups but also occasionally interact with groove atoms of both nucleic acids. However, the number of interactions in PEI/DNA complexes is greater than in comparable PEI/siRNA complexes, with interactions between protonated PEI amines and DNA being particularly enhanced. These results indicate that structural differences between DNA and siRNA may play a role in the increased stability of PEI/DNA complexes. In addition, we investigate the binding of PEI chains to polyplexes that have a net positive charge. The binding of PEI to these overcharged complexes involves interactions between PEI and areas on the nucleic acid surface that have maintained a negative electrostatic potential and is facilitated by the release of water from the nucleic acid.

Quantifying the Interactions between PEI and Double-Stranded DNA: Toward the Understanding of the Role of PEI in Gene Delivery

Poly(ethylene imine) (PEI) is one of the most efficient nonviral vectors, and its binding mode/strength with double-stranded DNA (dsDNA), which is still not clear, is a core area of transfection studies. In this work we used the atomic force microscopy (AFM)-based single molecule force spectroscopy (SMFS) to detect the interaction between branched PEI and dsDNA quantitatively by using a long chain DNA as a probe. Our results indicate that PEI binds to phosphoric acid skeletons of dsDNA mainly via electrostatic interactions, no obvious groove-binding or intercalation has happened. The interaction strength is about 24-25 pN, and it remains unchanged at pH 5.0 and 7.4, which correspond to the pH values in lysosomes and in the cytoplasmic matrix, respectively. However, the interaction is found to be sensitive to the ionic strength of the environment. In addition, the unbinding force shows no obvious loading rate dependence indicative of equilibrium binding/unbinding.

Evidence for plasmid DNA exchange after polyplex mixing

The self-assembly of a plasmid DNA (pDNA) with cationic polymers or cationic liposomes forms nanosized supramolecular structures called lipoplexes, polyplexes and lipopolyplexes. Here, we report that when two polyplex preparations made using the same polymer and the same pDNA but labelled with two different fluorophores are mixed together, pDNA molecules are exchanged. Indeed, when Flu-pDNA complexed with histidinylated lPEI (Flu-pDNA/His-lPEI) polyplexes are mixed with Cy5-pDNA complexed with histidinylated lPEI (Cy5-pDNA/His-lPEI) polyplexes, a high quantity of polyplexes emitting dual fluorescence is observed and FRET indicates that one single polyplex contains two kinds of fluorescent pDNA molecules. This phenomenon depends on the polymer-type and the strength of the pDNA/polymer interaction. No exchange is observed with polylysine polyplexes, caged His-lPEI polyplexes, lipoplexes, lipopolyplexes or when His-lPEI polyplexes are mixed with lipoplexes. Our results suggest that aggregation or collapse of polyplexes occurs after their interaction leading to their unpackaging followed by the formation of new polyplexes with the exchange of pDNA.

Live Cell Characterization of DNA Aggregation Delivered through Lipofection

DNA trafficking phenomena, such as information on where and to what extent DNA aggregation occurs, have yet to be fully characterised in the live cell. Here we characterise the aggregation of DNA when delivered through lipofection by applying the Number and Brightness (N&B) approach. The N&B analysis demonstrates extensive aggregation throughout the live cell with DNA clusters in the extremity of the cell and peri-nuclear areas. Once within the nucleus aggregation had decreased 3-fold. In addition, we show that increasing serum concentration of cell media results in greater cytoplasmic aggregation. Further, the effects of the DNA fragment size on aggregation was explored, where larger DNA constructs exhibited less aggregation. This study demonstrates the first quantification of DNA aggregation when delivered through lipofection in live cells. In addition, this study has presents a model for alternative uses of this imaging approach, which was originally developed to study protein oligomerization and aggregation.

Branched-linear polyion complexes investigated by Monte Carlo simulations

Complexes formed by one charged and branched copolymer with an oppositely charged and linear polyion have been investigated by Monte Carlo simulations. A coarse-grained description has been used, in which the main chain of the branched polyion and the linear polyion possess the same absolute charge and charge density. The spatial extension and other structural properties, such as bond-angle orientational correlation function, asphericity, and scaling analysis of formed complexes, at varying branching density and side-chain length of the branched polyion, have been explored. In particular, the balance between cohesive Coulomb attraction and side-chain repulsions resulted in two main structures of a polyion complex. These structures are (i) a globular polyion core surrounded by side chains appearing at low branching density and (ii) an extended polyion core with side chains still being expelled at high branching density. The globule-to-extended transition occurred at a crossover branching density being practically independent of the side chain length.

Lipid-modified polyethylenimine-mediated DNA attraction evaluated by molecular dynamics simulations

The effect of lipid modification on polyethylenimine (PEI)-mediated DNA attraction was studied by performing umbrella sampling molecular dynamics simulations that involved PEIs modified with three different types of lipids: oleic acid (OA), linoleic acid (LA), and caprylic acid (CA). The potential of mean force between two DNA molecules in the presence of these lipid-modified PEIs was calculated using the weighted histogram analysis method, and it predicted the stability and size of the DNA aggregate. When compared to native PEI, lipid modification was found to enhance the stability of DNA aggregation in the case of long lipids (LA and OA) but reduce the stability in the case of a short lipid (CA). In addition, LA-substituted PEI was shown to form stronger DNA aggregate than OA-substituted PEI, which correlates positively with previous experimental observations.

Molecular modeling of polynucleotide complexes

Delivery of polynucleotides into patient cells is a promising strategy for treatment of genetic disorders. Gene therapy aims to either synthesize desired proteins (DNA delivery) or suppress expression of endogenous genes (siRNA delivery). Carriers constitute an important part of gene therapeutics due to limitations arising from the pharmacokinetics of polynucleotides. Non-viral carriers such as polymers and lipids protect polynucleotides from intra and extracellular threats and facilitate formation of cell-permeable nanoparticles through shielding and/or bridging multiple polynucleotide molecules. Formation of nanoparticulate systems with optimal features, their cellular uptake and intracellular trafficking are crucial steps for an effective gene therapy. Despite the great amount of experimental work pursued, critical features of the nanoparticles as well as their processing mechanisms are still under debate due to the lack of instrumentation at atomic resolution. Molecular modeling based computational approaches can shed light onto the atomic level details of gene delivery systems, thus provide valuable input that cannot be readily obtained with experimental techniques. Here, we review the molecular modeling research pursued on critical gene therapy steps, highlight the knowledge gaps in the field and providing future perspectives. Existing modeling studies revealed several important aspects of gene delivery, such as nanoparticle formation dynamics with various carriers, effect of carrier properties on complexation, carrier conformations in endosomal stages, and release of polynucleotides from carriers. Rate-limiting steps related to cellular events (i.e. internalization, endosomal escape, and nuclear uptake) are now beginning to be addressed by computational approaches. Limitations arising from current computational power and accuracy of modeling have been hindering the development of more realistic models. With the help of rapidly-growing computational power, the critical aspects of gene therapy are expected to be better investigated and direct comparison between more realistic molecular modeling and experiments may open the path for design of next generation gene therapeutics.

A molecular dynamics simulation study on the effect of lipid substitution on polyethylenimine mediated siRNA complexation

Polycations have been explored as non-viral carriers for effective delivery of small interfering RNA (siRNA). Modifying polycations such as polyethylenimine (PEI) with lipid substitution was found to improve the siRNA delivery efficiency of polycationic carriers. However, the role of such lipid modification is not clear and remains to be probed at the atomistic level. In this work, we elucidate the role of lipid modification through a series of all-atom molecular dynamics simulations on siRNA complexation mediated by a native PEI and four analogous obtained by different lipid modifications. The lipid modification does not affect PEI's capability of neutralizing the siRNA charge, neither does it affect the polyion bridging which plays an important role in siRNA complexation. Significant linkages among the lipid modified PEIs via association of lipid side-groups are observed and this results in more stable and compact PEI/siRNA polyplexes. The lipid associations between short lipids form and break frequently while the lipid associations between long lipids are more stable. For PEIs modified with short lipids, increasing the lipid substitution level results in more compact and stable siRNA structure. For PEIs modified with long lipids, increasing the lipid substitution does not change the amount of PEI linkage via lipid association, and has a reverse effect on compacting siRNA structure due to increased steric hindrance brought by the lipid association on individual PEIs. The simulation results generally correlate well with experimental data and suggest a framework of designing and systematic evaluation of polycation-based siRNA carriers using molecular dynamics simulations.

Multiscale Simulation as a Framework for the Enhanced Design of Nanodiamond-Polyethylenimine-based Gene Delivery

Nanodiamonds (NDs) are emerging carbon platforms with promise as gene/drug delivery vectors for cancer therapy. Specifically, NDs functionalized with the polymer polyethylenimine (PEI) can transfect small interfering RNAs (siRNA) in vitro with high efficiency and low cytotoxicity. Here we present a modeling framework to accurately guide the design of ND-PEI gene platforms and elucidate binding mechanisms between ND, PEI, and siRNA. This is among the first ND simulations to comprehensively account for ND size, charge distribution, surface functionalization, and graphitization. The simulation results are compared with our experimental results both for PEI loading onto NDs and for siRNA (C-myc) loading onto ND-PEI for various mixing ratios. Remarkably, the model is able to predict loading trends and saturation limits for PEI and siRNA, while confirming the essential role of ND surface functionalization in mediating ND-PEI interactions. These results demonstrate that this robust framework can be a powerful tool in ND platform development, with the capacity to realistically treat other nanoparticle systems.

Molecular dynamics simulations for complexation of DNA with 2 kDa PEI reveal profound effect of PEI architecture on complexation

A series of all-atom molecular dynamics (MD) simulations of the complexation between DNA and 2 kDa branched and linear polyethylenimines (PEIs) are reported in this study. The simulations revealed distinct binding modes of branched and linear PEIs to DNA, with branched PEIs adhering to the DNA surface like beads and linear PEIs adhering to the DNA surface like cords. The dynamics of each PEI's binding state to the DNA during the simulations and how the PEIs neutralize the DNA were quantified. For both branched and linear PEIs, the addition of salt ions similar to physiological conditions were found to have only a small effect on DNA/PEI complexation compared to salt-free conditions. The simulation results reported here will be helpful to understand the mechanism of action for the PEI-based gene carriers.