Potential of mean force of polyethylenimine-mediated DNA attraction.

The aggregation of DNA molecules induced by cationic polymers is of importance to applications in gene delivery. In this work, we performed a series of umbrella sampling molecular dynamics simulations to calculate the potential of mean force (PMF) between two DNA molecules in the presence of polyethylenimine (PEI) molecules using the weighted histogram analysis method. The distance between the centers of mass of the two DNAs was chosen as the reaction coordinate, and the location and depth of the global minimum in the PMF curve were used to gauge the compactness and stability of the formed aggregate. The effects of PEI to DNA charge ratio (N/P charge ratio) and protonation state of the PEI were investigated. The DNA aggregation was found to be more favorable at higher PEI/DNA charge ratios and higher PEI protonation ratios. As compared with small multivalent ions, PEIs give rise to stronger DNA attraction even with an N/P charge ratio at which the DNAs are overneutralized. The effect of changing the protonation ratio on the compactness of the aggregate is more prominent than changing the N/P charge ratio, while the depth of the PMF well is more strongly influenced by the N/P charge ratio.

[1]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[2]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[3]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[4]  S. Broyde,et al.  Base flipping free energy profiles for damaged and undamaged DNA. , 2010, Chemical research in toxicology.

[5]  G. S. Manning.  The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides , 1978, Quarterly Reviews of Biophysics.

[6]  R. W. Wilson,et al.  Counterion-induced condesation of deoxyribonucleic acid. a light-scattering study. , 1979, Biochemistry.

[7]  Serdar Kuyucak,et al.  Potential of mean force calculations of ligand binding to ion channels from Jarzynski's equality and umbrella sampling. , 2008, The Journal of chemical physics.

[8]  J. Behr,et al.  Optimized galenics improve in vitro gene transfer with cationic molecules up to 1000-fold. , 1996, Gene therapy.

[9]  L. Monaco,et al.  Nanoscopic structure of DNA condensed for gene delivery. , 1997, Nucleic acids research.

[10]  Wenjin Guo,et al.  Efficient Gene Delivery Using Anionic Liposome-Complexed Polyplexes (LPDII) , 2000, Bioscience reports.

[11]  R. Golestanian,et al.  Condensation of DNA-actin polyelectrolyte mixtures driven by ions of different valences. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

[13]  F. Oosawa.  Interaction between parallel rodlike macroions , 1968 .

[14]  H. Uludaǧ,et al.  Molecular dynamics simulations of PEI mediated DNA aggregation. , 2011, Biomacromolecules.

[15]  H. Mao,et al.  Charge density and molecular weight of polyphosphoramidate gene carrier are key parameters influencing its DNA compaction ability and transfection efficiency. , 2010, Biomacromolecules.

[16]  A. Kornyshev,et al.  Electrostatic Zipper Motif for DNA Aggregation , 1999 .

[17]  Martin Egli,et al.  The Dickerson-Drew B-DNA Dodecamer Revisited-At Atomic Resolution , 1998 .

[18]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[19]  F. Livolant,et al.  DNA condensed by protamine: a "short" or "long" polycation behavior. , 2009, Biomacromolecules.

[20]  B. Shklovskii.  WIGNER CRYSTAL MODEL OF COUNTERION INDUCED BUNDLE FORMATION OF RODLIKE POLYELECTROLYTES , 1998, cond-mat/9809429.

[21]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[22]  S. Khalid,et al.  DNA and lipid bilayers: self-assembly and insertion , 2008, Journal of The Royal Society Interface.

[23]  D. Scherman,et al.  A versatile vector for gene and oligonucleotide transfer into cells in culture and in vivo: polyethylenimine. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Yasuhiko Tabata,et al.  N/P ratio significantly influences the transfection efficiency and cytotoxicity of a polyethylenimine/chitosan/DNA complex. , 2009, Biological & pharmaceutical bulletin.

[25]  D. Fischer,et al.  A Novel Non-Viral Vector for DNA Delivery Based on Low Molecular Weight, Branched Polyethylenimine: Effect of Molecular Weight on Transfection Efficiency and Cytotoxicity , 1999, Pharmaceutical Research.

[26]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[27]  W. L. Jorgensen.  Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water , 2002 .

[28]  Hsin‐Lung Chen,et al.  Effect of divalent cations on DNA condensed on the surface of rigid cationic membrane , 2011 .

[29]  Serdar Kuyucak,et al.  Accurate determination of the binding free energy for KcsA-charybdotoxin complex from the potential of mean force calculations with restraints. , 2011, Biophysical journal.

[30]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[31]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[32]  T. Andresen,et al.  Elucidating the interplay between DNA-condensing and free polycations in gene transfection through a mechanistic study of linear and branched PEI. , 2011, Biomaterials.

[33]  A. Mikos,et al.  Size matters: molecular weight affects the efficiency of poly(ethylenimine) as a gene delivery vehicle. , 1999, Journal of biomedical materials research.

[34]  V. Bloomfield.  DNA condensation by multivalent cations. , 1997, Biopolymers.

[35]  T. Hellweg,et al.  Interaction of short DNA fragments with the cationic polyelectrolyte poly(ethylene imine): a dynamic light scattering study , 2000 .

[36]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

[37]  Jesse D. Ziebarth,et al.  Molecular dynamics simulations of DNA-polycation complex formation. , 2009, Biophysical journal.

[38]  V. Labhasetwar,et al.  Size-dependency of nanoparticle-mediated gene transfection: studies with fractionated nanoparticles. , 2002, International journal of pharmaceutics.

[39]  K. Yoshikawa,et al.  Direct observation of the biphasic conformational change of DNA induced by cationic polymers , 1991, FEBS letters.

[40]  J. DeRouchey,et al.  Cation charge dependence of the forces driving DNA assembly. , 2010, Biophysical journal.

[41]  M. Cotten,et al.  Transferrin-polycation-DNA complexes: the effect of polycations on the structure of the complex and DNA delivery to cells. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[42]  J. Dubochet,et al.  Brownian dynamics simulation of DNA condensation. , 1999, Biophysical journal.

[43]  L. Nordenskiöld,et al.  A Monte Carlo simulation study of electrostatic forces between hexagonally packed DNA double helices , 1986 .

[44]  Richard Lavery,et al.  Base pair opening within B-DNA: free energy pathways for GC and AT pairs from umbrella sampling simulations. , 2003, Nucleic acids research.

[45]  H. Uludaǧ,et al.  Molecular dynamics simulations of DNA/PEI complexes: effect of PEI branching and protonation state. , 2011, Biophysical journal.

[46]  Feng Yan,et al.  The use of polyethylenimine-grafted graphene nanoribbon for cellular delivery of locked nucleic acid modified molecular beacon for recognition of microRNA. , 2011, Biomaterials.

[47]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[48]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[49]  J. Schellman,et al.  Compact form of DNA induced by spermidine , 1976, Nature.

[50]  Y. Oh,et al.  Polyethylenimine-mediated cellular uptake, nucleus trafficking and expression of cytokine plasmid DNA , 2002, Gene Therapy.

[51]  N. Kasyanenko,et al.  DNA Self-Assembling Nanostructures Induced by Trivalent Ions and Polycations , 2008 .

[52]  Molecular dynamics simulation of multivalent-ion mediated attraction between DNA molecules. , 2008, Physical review letters.

[53]  Gianni De Fabritiis,et al.  Optimized Potential of Mean Force Calculations for Standard Binding Free Energies. , 2011, Journal of chemical theory and computation.

[54]  V A Parsegian,et al.  Direct measurement of the intermolecular forces between counterion-condensed DNA double helices. Evidence for long range attractive hydration forces. , 1992, Biophysical journal.