Structural and dynamical properties of Polyethylenimine in explicit water at different protonation states: A Molecular Dynamics Study

Similar documents
SUPPLEMENTAL MATERIAL

All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010

Supplemntary Infomation: The nanostructure of. a lithium glyme solvate ionic liquid at electrified. interfaces

k θ (θ θ 0 ) 2 angles r i j r i j

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

The Molecular Dynamics Method

Structural Bioinformatics (C3210) Molecular Mechanics

Electronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani

An introduction to Molecular Dynamics. EMBO, June 2016

Coarse-Grained Models!

Molecular mechanics. classical description of molecules. Marcus Elstner and Tomáš Kubař. April 29, 2016

Lecture 1. Conformational Analysis in Acyclic Systems

Force Fields for MD simulations

The change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:

Supporting Information

q lm1 q lm2 q lm3 (1) m 1,m 2,m 3,m 1 +m 2 +m 3 =0 m 1 m 2 m 3 l l l

Molecular dynamics simulation of Aquaporin-1. 4 nm

Ionic Behavior in Highly Concentrated Aqueous Solutions Nanoconfined between Discretely Charged Silicon Surfaces

Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?

Mechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations

Citation for published version (APA): Sok, R. M. (1994). Permeation of small molecules across a polymer membrane: a computer simulation study s.n.

Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information

Supporting Information. How Different are the Characteristics of Aqueous Solutions of tert-

Lecture 11: Potential Energy Functions

Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure

Bioengineering 215. An Introduction to Molecular Dynamics for Biomolecules

Crystal Polymorphism in Hydrophobically Nanoconfined Water. O. Vilanova, G. Franzese Universitat de Barcelona

Appendix C - Persistence length 183. Consider an ideal chain with N segments each of length a, such that the contour length L c is

Morphology of Supported Polymer Electrolyte. Ultra-thin Films: a Numerical Study. Supporting Information

ACMAC s PrePrint Repository

Chapter 6 Cyclic urea - a new central unit in bent-core compounds

Supporting Information

arxiv: v1 [cond-mat.soft] 22 Oct 2007

Supporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes

MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE

Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together.

Subject of the Lecture:

The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity

A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane

This semester. Books

arxiv:cond-mat/ v1 2 Feb 94

The Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii

Autocorrelation Study for a Coarse-Grained Polymer Model

Molecular Dynamics Simulation of a Nanoconfined Water Film

Chapter 4. Glutamic Acid in Solution - Correlations

Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate

Adsorption of gases on solids (focus on physisorption)

Why Is Molecular Interaction Important in Our Life

Dihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769

Supporting Information

PH 548 Atomistic Simulation Techniques

Statistical Thermodynamics Exercise 11 HS Exercise 11

Molecular Modelling. part of Bioinformatik von RNA- und Proteinstrukturen. Sonja Prohaska. Leipzig, SS Computational EvoDevo University Leipzig

Pancake Bonding in 1,2,3,5-dithiadiazolyl and 1,2,3,5-diselenadiazolyl Radical Dimers and their Derivatives

Part III : M6 Polymeric Materials

arxiv: v1 [cond-mat.stat-mech] 6 Jan 2014

Supporting Information

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

Lattice protein models

a) DIATOMIC ELEMENTS e.g. . MEMORIZE THEM!!

Nanotube AFM Probe Resolution

Intermolecular distance and density scaling of dynamics in molecular liquids

Physics 9 Summer 2010 Midterm

Biomolecular modeling I

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea

CECAM Tutorial Stuttgart, Oct. 2012

Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement

The micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction

Journal of Physics: Conference Series PAPER. To cite this article: Jiramate Kitjanon et al 2017 J. Phys.: Conf. Ser

PHASE TRANSITIONS IN SOFT MATTER SYSTEMS

Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations

Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole

Statistical and hydrodynamic properties of topological polymers for various graphs showing enhanced short-range correlation

LAMMPS Performance Benchmark on VSC-1 and VSC-2

Free energy calculations of small molecules in dense amorphous polymers. Effect of the initial guess configuration in molecular dynamics studies

Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials

Pharmaceutical co-crystals of diflunisal and. diclofenac with theophylline

Physics 211B : Problem Set #0

Journal of Chemical and Pharmaceutical Research, 2012, 4(3): Research Article

Structure of the First and Second Neighbor Shells of Water: Quantitative Relation with Translational and Orientational Order.

Lecture 10. Born-Oppenheimer approximation LCAO-MO application to H + The potential energy surface MOs for diatomic molecules. NC State University

Applications of Molecular Dynamics

Homework Problem Set 1 Solutions

Atomic Structure & Interatomic Bonding

Water structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials

Chapter 2 Experimental sources of intermolecular potentials

Electron Configuration & Structure

A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE

Supporting Information. First principles kinetic study on the effect of zeolite framework on 1-butanol dehydration

U N I T T E S T P R A C T I C E

Polymer-specific effects of bulk relaxation and stringlike correlated motion in the dynamics of a supercooled polymer melt

L718Q mutant EGFR escapes covalent inhibition by stabilizing. a non-reactive conformation of the lung cancer drug. osimertinib

MOLECULAR INTERACTIONS NOTES

49 th INTERNATIONAL CHEMISTRY OLYMPIAD 2017 UK Round One MARK SCHEME

PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS

Practice Exam 2 Key Fall 2014

Advanced Subsidiary Unit 1: The Core Principles of Chemistry

EGN 3365 Review on Bonding & Crystal Structures by Zhe Cheng

Potential Energy (hyper)surface

Transcription:

This journal is The Royal Society of Chemistry Structural and dynamical properties of Polyethylenimine in explicit water at different protonation states: A Molecular Dynamics Study Chandan Kumar Choudhury a and Sudip Roy a Received Xth XXXXXXXXXX XX, Accepted Xth XXXXXXXXX XX First published on the web Xth XXXXXXXXXX X DOI:.39/bx a Physical Chemistry Division, ational Chemical Laboratory, Pune, India. Fax: +9 () 5965; Tel: +9 () 59 735; E-mail: s.roy@ncl.res.in

This journal is The Royal Society of Chemistry Table S Coulombic parameters (q i, see Equation ) for PEI chain. See Figure S for atom numbers. Atom o. Atom name Atom type Charge, q i amber99 38 () -.74, 3, amber99 7 ().5 4 CA amber99 (CT) -.4 5, 6, 8, 9 A, A, B, B amber99 9 ().5 7 CB amber99 (CT) -.8 amber99 8 (T) -.593, amber99 7 ().69, 5 CA, CB amber99 (CT) -.86 3, 4, 6, 7 A, A, B, B amber99 9 ().4 8 amber99 8 (T) -.5 9, amber99 7 ().8 CA amber99 (CT) -.8, A, A amber99 9 ().5 3 CB amber99 (CT) -.3 4, 5, 6 B, B, B3 amber99 8 (C).7

This journal is The Royal Society of Chemistry Table S Coulombic parameters (q i, see Equation ) for alternate protonated PEI chain. See Figure S for atom numbers. Atom o. Atom name Atom type Charge, q i amber99 8 (T) -.47, 3, 4,, 3 amber99 7 ().36 5 CA amber99 (CT).56 6, 7 A, A amber99 9 ().87 8 CB amber99 (CT) -.47 9, B, B amber99 9 ().7 amber99 8 (T) -.593 amber99 7 ().69 3, 6 CA, CB amber99 (CT) -.86 4, 5, 7, 8 A, A, B, B amber99 9 ().4 9 amber99 8 (T) -.5,, amber99 7 ().34, 5 CA, CB amber99 (CT) -.3 3, 4, 6, 7 A, A, B, B amber99 9 ().9 8 amber99 8 (T) -.5 9 amber99 7 ().8 3 CA amber99 (CT) -.8 3, 3 A, A amber99 9 ().5 33 CB amber99 (CT) -.3 34, 35, 36 B, B, B3 amber99 8 (C).7 3

This journal is The Royal Society of Chemistry Table S3 Coulombic parameters (q i, see Equation ) for completely protonated PEI chain. See Figure S(c) for atom numbers. Atom o. Atom name Atom type Charge, q i amber99 8 (T) -.47, 3, 4,, 3 amber99 7 ().36 5 CA amber99 (CT).56 6, 7 A, A amber99 9 ().87 8 CB amber99 (CT) -.47 9, B, B amber99 9 ().7 amber99 8 (T) -.5, 3, amber99 7 ().34 4, 7 CA, CB amber99 (CT) -.3 5, 6, 8, 9 A, A, B, B amber99 9 ().9 amber99 8 (T) -.35,, amber99 7 ().33 3 CA amber99 (CT) -.67 4, 5 A, A amber99 9 (). 6 CB amber99 (CT) -.34 7, 8, 9 B, B, B3 amber99 8 (C).6 4

This journal is The Royal Society of Chemistry Table S4 Lennard Jones parameters (σ and ε, see Equation ) for the atom types mentioned in Table S, Table S, and Table S3. Atom type σ, nm ε, kj mol amber99 (CT).339967.45773 amber99 7 ().698.656888 amber99 8 (C).64953.656888 amber99 9 ().4735.656888 amber99 38 ().35.78 amber99 8 (T).35.78 Table S5 Bond parameters (b and ki bonds j, see Equation ) for the bond between the atom types (i j) in Table S, Table S and Table S3. Bond (Atom types, i j) b, nm ki bonds j, kj mol nm -. 3637. - CT.463 8.6 CT -.9 845. CT - CT.56 5948. CT - T.47 375.6 T -. 3637. CT - C.9 845. 5

This journal is The Royal Society of Chemistry Table S6 Bonded angle parameters (θ and k angle i jk ) for angles formed by the atom types (i jk) mentioned in Table S, Table S and Table S3. Angles (Atom types, i jk) θ, k angle i jk, kj mol rad - -. 9.88 - - CT 8.4 48.4 - CT - 9.5 48.4 - CT - CT. 669.4 CT - CT - 9.5 48.4 - CT - 9.5 9.8 CT - CT - T. 669.4 CT - T - 9.5 48.4 T - CT - 9.5 48.4 CT - T - CT 9.5 48.4 CT - CT - C 9.5 48.4 C - CT - C 9.5 9.8 - T - 9.5 9.88 Table S7 Dihedral parameters (φ s and ki dihedral jkl, see Equation ) for armonic functions and the values of C n for Ryckaert-Bellemans functions. X denotes any undefined atom types in the harmonic functions in Table S, Table S and Table S3 Dihedral (Atom types), i jkl φ s ki dihedral jkl (kj mol ) n - - CT - CT. 3.7656 3 - T - CT - CT. 3.7656 3 - T - CT -. 3.7656 3 CT - CT - T - CT..676 3 CT - CT - T - CT 8..4 Dihedral (Atom types), i jkl, kj mol C C C C 3 C 4 C 5 X - - CT - X...... X - CT - CT - X.6584.9553. -.6338.. Table S8 Possible nearest neighbor (a, b) and next nearest neighbour (c, d and e) conformations for the PEI chain (c) (d) (e) 6

This journal is The Royal Society of Chemistry 5 6 Repeating Unit 6 7 4 4 7 5 n 8 3 5 3 8 9 3 4 9 6 -Terminal C-Terminal 4 3 6 7 5 8 9 4 Repeating Unit 7 8 3 3 6 9 n 5 4 5 n 9 33 8 3 6 7 3 3 34 35 36 -Terminal C-Terminal 4 3 (c) 6 7 Repeating Unit 3 8 9 4 5 8 4 7 n 6 9 5 3 5 6 7 9 8 -Terminal C-Terminal Fig. S Structure of PEI with atom number. Unprotonated chain, Alternate protonated, and (c) chain. For and chain n = 8 for -mer and 48 for 5-mer, while for alternated protonated chain, n = 9 and 4 for and 5-mer, respectively. The atoms are numbered for Table S,Table S and Table S3. 7

This journal is The Royal Society of Chemistry 4 36 -mer 5-mer p(r) Density (kg m -3 ) 3 8 4 7 7 7 73 74 75 76 77 78 79 8 Time (ns) Fig. S The variation of density of and 5-mer PEI with time. 8 6 End-to-end distance (nm) 8 7 6 5 4 3 5 5 3 35 4 45 5 Time (ns) 4 p(r)..5.3.35.4.45.5 Distance between, r (nm) 4 8 6 4.4.3.. 3 4 5 6 7 End-to-end distance (nm) Fig. S4 End-to-end distance and its distribution of -mer PEI for different protonation states..3.4.5.6.7.8 Distance between, r (nm) Fig. S3 The nearest neighbour, and The next nearest neighbour, distributions of the nitrogen atoms for the -mer PEI. 8

This journal is The Royal Society of Chemistry Radius of gyration (R g ) (nm).5.5.5.4.3.. unprototnated...3.4.5.6.7.8.9...3 Ratio (syn/trans).5 5 5 3 35 4 45 5 Time (ns).4.3.5.4...3. 4 6 8 4 6 8 Ratio (syn/trans)..6.8..4.6.8 Radius of gyration (R g ) (nm).5.4 (c) Fig. S5 Radius of Gyration and its distribution of -mer PEI for different protonation states..3.. Probability.5.4.3.. un [ -C-C- ] un [ C-C--C ] alternate [ -C-C- ] alternate [ C-C--C ] completely [ -C-C- ] completely [ C-C--C ]..5..5.3.35.4.45.5.55.6 Ratio (syn/trans).5.4.3. (d) - -5 - -5 5 5 Dihedral angle (Degrees) Fig. S6 Dihedral distribution for -mer PEI for different protonation states.. 3 4 5 6 7 8 9 3 4 alternate syn-anti pairs Fig. S7 Distribution of ratio (syn/anti) when all, -C-C- and (c) C-C--C dihedrals were considered. (d) Alternate syn-anti occurrences of all the dihedrals for -mer PEI. 9

This journal is The Royal Society of Chemistry g(r) 3.5.5.5 Residence time probability.4.35.3.5..5..5 st solvation st solvation protonated nd solvation nd solvation protonated..3.4.5.6.7.8.9.. r (nm) 5 5 5 3 35 4 45 5 Time (ps) g(r).5.5 Residence time probability.5.45.4.35.3.5..5. st solvation st solvation protonated nd solvation nd solvation protonated.5..3.4.5.6.7.8.9.. r (nm) 5 5 5 3 35 4 45 5 Time (ps) g(r).4..8.6.4. (c)..3.4.5.6.7.8.9.. r (nm) Fig. S8 The pair distribution plot of of PEI and O of Water, C of PEI and O of Water and (c) and C of PEI and O of Water for -mer PEI. Residence time probability.5.45.4.35.3.5..5..5 (c) st solvation st solvation protonated nd solvation nd solvation protonated 5 5 5 3 35 4 45 5 Time (ps) Fig. S9 Residence time probability of the water molecules in different solvation shells of the polymer chain (5-mer) at 3 K 3 K and (c) 33 K.

This journal is The Royal Society of Chemistry MSD (nm ).5..5. 3 K 3 K 3 K 33 K 3 K 3 K 3 K 33 K.5 4 6 8 4 6 8 Time (ps) MSD (nm ).5.4.3. 3 K 3 K 3 K 33 K 3 K 3 K 3 K 33 K. 4 6 8 4 6 8 Time (ps) Fig. S Mean square displacement of water molecules around the st solvation shell, and nd solvation shell at different temperatures for (solid lines) and protonated (dotted - dashed lines) states of 5-mer PEI.