Understanding the Protonation Behavior of Linear Polyethylenimine in Solutions through Monte Carlo Simulations
- 2 December 2009
- journal article
- research article
- Published by American Chemical Society (ACS) in Biomacromolecules
- Vol. 11 (1), 29-38
- https://doi.org/10.1021/bm900842d
Abstract
The success of polyethyleneimine (PEI) as a nonviral-based gene delivery vector has been attributed to its proton buffering capacity. Despite the great interest in PEI for its use in nonviral-based gene delivery, the protonation behavior of PEI in solution is not well understood. Earlier experimental studies have reported inconsistent values of the protonation state of PEI. In this work, we report our investigation of the protonation behavior of a realistic linear PEI (lPEI) with computational approaches. Reported experimental pKa values of several diamine compounds are first examined. A screened Coulombic interaction with a distance dependence dielectric is shown to reproduce the shifted pKa values of the model diamine compounds. Then atomistic molecular dynamic simulations of lPEI chain with 20 repeating units are performed and the results are used to provide parameters for a coarse-grained polyamine model. The screened Coulombic interaction is then incorporated in the coarse-grained lPEI chain and computational titrations are performed. The obtained computational titration curves of lPEI in solutions were found to be in best agreement with experimental results by Smits et al., but the computational titration curves have too strong of a dependence on salt concentration compared to the experimental results by Smits et al. Disregarding the discrepancy in the salt dependence, our computational titrations reveal that approximately 55% of the lPEI amine groups are protonated under physiological conditions in solution with a nearly alternating arrangement of protonated and nonprotonated amines. Titrations of lPEI in the presence of a polyanion are also performed to determine how the charge state of lPEI could be affected by complexation with DNA in gene therapy preparations. While the presence of the polyanion increases the degree of protonation of the PEI, many of PEI amines remain unprotonated under physiological conditions, providing evidence that PEI complexed with DNA could still have proton buffering capacity. Potential sources of error that have resulted in the inconsistency of previously reported protonation states of PEI were also discussed.Keywords
This publication has 53 references indexed in Scilit:
- Molecular Dynamics Simulations of DNA-Polycation Complex FormationBiophysical Journal, 2009
- Titration of hydrophobic polyelectrolytes using Monte Carlo simulationsThe Journal of Chemical Physics, 2005
- Adsorption of a weakly charged polymer on an oppositely charged colloidal particle: Monte Carlo simulations investigationPolymer, 2005
- Determination of the effective dielectric constant from the accurate solution of the Poisson equationJournal of Computational Chemistry, 2002
- Proton Binding Characteristics of Branched PolyelectrolytesMacromolecules, 1997
- A Monte Carlo Study of Titrating Polyelectrolytes in the Presence of SaltMacromolecules, 1996
- Ising Models of Polyprotic Acids and BasesThe Journal of Physical Chemistry, 1994
- A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP modelThe Journal of Physical Chemistry, 1993
- Monte Carlo study of titration of linear polyelectrolytesThe Journal of Chemical Physics, 1992
- Distance‐dependent dielectric constants and their application to double‐helical DNAPeptide Science, 1991