Next Article in Journal
Naturally-Derived Amphiphilic Polystyrenes Prepared by Aqueous Controlled/Living Cationic Polymerization and Copolymerization of Vinylguaiacol with R–OH/BF3·OEt2
Previous Article in Journal
Inkjet Printable and Self-Curable Disperse Dyes/P(St-BA-MAA) Nanosphere Inks for Both Hydrophilic and Hydrophobic Fabrics
Previous Article in Special Issue
Multiscale Simulation of Branched Nanofillers on Young’s Modulus of Polymer Nanocomposites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantifying Mg2+ Binding to ssDNA Oligomers: A Self-Consistent Field Theory Study at Varying Ionic Strengths and Grafting Densities

1
Department of Chemical Engineering, University of South Carolina, Columbia, SC 29208, USA
2
Biomedical Engineering Program, University of South Carolina, Columbia, SC 29208, USA
*
Author to whom correspondence should be addressed.
Polymers 2018, 10(12), 1403; https://doi.org/10.3390/polym10121403
Submission received: 15 November 2018 / Revised: 12 December 2018 / Accepted: 14 December 2018 / Published: 18 December 2018
(This article belongs to the Special Issue Multiscale Modeling of Polymers)

Abstract

:
The performance of aptamer-based biosensors is crucially impacted by their interactions with physiological metal ions, which can alter their structures and chemical properties. Therefore, elucidating the nature of these interactions carries the utmost importance in the robust design of highly efficient biosensors. We investigated Mg 2 + binding to varying sequences of polymers to capture the effects of ionic strength and grafting density on ion binding and molecular reorganization of the polymer layer. The polymers are modeled as ssDNA aptamers using a self-consistent field theory, which accounts for non-covalent ion binding by integrating experimentally-derived binding constants. Our model captures the typical polyelectrolyte behavior of chain collapse with increased ionic strength for the ssDNA chains at low grafting density and exhibits the well-known re-entrant phenomena of stretched chains with increased ionic strength at high grafting density. The binding results suggest that electrostatic attraction between the monomers and Mg 2 + plays the dominant role in defining the ion cloud around the ssDNA chains and generates a nearly-uniform ion distribution along the chains containing varying monomer sequences. These findings are in qualitative agreement with recent experimental results for Mg 2 + binding to surface-bound ssDNA.

Graphical Abstract

1. Introduction

Aptamers are an important class of biomolecules consisting of single-stranded DNA (ssDNA), RNA, or peptides that can fold into unique secondary and tertiary structures for shape-specific target recognition [1]. Due to the highly specific and selective nature of their target binding, aptamers are widely studied for a range of applications from biosensing [2,3] to drug design [4,5,6]. A recent work reported a major breakthrough in biosensor research by using aptamers with field-effect transistors to overcome the Debye length limitations [7]. Aptamers are polyelectrolytic in nature with their monomer units (nucleobases) participating in acid-base equilibrium and counterion binding reactions with the surrounding solution environments. Charge regulation and counterion binding in aptamers, or polyelectrolytes in general, are modulated by the metal ions present in the system, which can non-trivially alter their chemical and structural properties [8,9,10,11,12,13]. The presence of metal ions affects the performance of the aptamers as biosensing probes or therapeutics [14,15,16] due to the electrostatic screening of the charges on their surface, which changes their structure and chemistry. These interactions of aptamers with metal ions are complicated in nature owing to the fact that multiple binding sites on the nucleobases are capable of such interactions, following different binding pathways and thus having varying energy landscapes [17,18]. In this work, we have particularly addressed magnesium ion (Mg 2 + ) binding because of its relevance to almost all nucleic acid-related biological processes in the intracellular environment [19,20,21].
A myriad of computational studies has been conducted with Molecular Dynamics (MD) [22,23,24,25,26,27,28,29,30] and Monte Carlo (MC) [31] simulations to elucidate the nature of metal ion binding to nucleic acids. Most of such theoretical studies are based on double-stranded DNA (dsDNA)–monovalent cation (such as Na + , K + ) interactions [25,27,29,30]. Among the few that include multivalent cations, Hayes et al. [26] employed a hybrid structure-based MD model to count explicitly the number of excess Mg 2 + ions bound to RNA sequences in the presence of background potassium chloride with a Manning condensation estimated by a non-linear Poisson–Boltzmann equation. Li et al. [28] used implicit Mg 2 + binding to dsDNA sequences with classical MD simulation to study the effect of counterion condensation on DNA structure and conformational dynamics. While atomistic MD simulations give a full distribution of the ion atmosphere around the nucleic acids, they suffer from drawbacks due to the enormous computational cost and the choice of force fields that might lead to over- or under-estimation of the same ion cloud [27,32]. These studies also rely heavily on parameterization to match experimental studies, which imposes unrestrained bias toward their agreement with experimental results [32]. On the other hand, almost all theoretical studies consider nucleic acid chains in bulk conditions; therefore, characteristics of nucleic acid strands end-tethered to a surface remain elusive.
Along this line, we address metal ion binding to surface-anchored nucleic acid oligomers with a Self-Consistent Field Theory (SCFT) approach to construct a comprehensive and statistically-robust model for quantifying the number of Mg 2 + ions bound to each chain, while capturing the ion-binding effect on their structure and properties. The molecular model analyzes Mg 2 + binding to nucleic acid oligomers containing adenine (A) and guanine (G) nucleobases, while capturing the details of experimental studies for similar systems as precisely as possible. Metal ion binding to the monomers is explicitly included with equilibrium binding reactions by using experimentally-derived [33] binding free energies for relevant binding modes [34]. The molecular model characterizes the spatial variation of the structure and properties of the oligonucleotide chains along the distance from the grafting surface, at varying ionic strength and grafting densities, and quantifies the number of bound ions at thermodynamic equilibrium with the oligonucleotides. The model explicitly accounts for the thermodynamic, structural, and electrostatic properties of all the species involved in the system, while remaining free of adjustable parameters. Despite the many benefits of using an SCFT approach to studying metallic interactions with oligonucleotides, there are several drawbacks. The first is that a form of the mean-field approximation must be used to model the intermolecular interactions. While much of the relevant physics related to correlations is present in our model, not all of the correlations are completely captured. As discussed in much greater detail later, this calculation also uses a model for building the configurational space of the oligonucleotides that could possibly be better represented by using simulation results, and we are performing a one-dimensional calculation that only treats inhomogeneities in the direction normal to the surface. Lateral heterogeneity would have to be accounted for in situations where there are appreciable defects in the grafting of the oligonucleotides to the surface, or if there is a need to capture microphase separation, where the system would be able to adopt a periodic two-dimensional structure. This field theoretic model helps to set up the foundation for future studies involving secondary and tertiary structures of aptamers interacting with multivalent metal ions.

2. Theoretical Methodology

The theoretical model in this work is developed to represent surface-grafted ssDNA oligomers as a co-polymer with two monomer units, adenine (A) and guanine (G), at a coarse-grained level. We explicitly consider the physical and chemical properties of these nucleobases in a solution environment to capture their behavior as accurately as possible. We study three chain sequences: the diblock co-polymer of A and G with the A-end grafted to the surface (A 6 G 6 , a block of six adenine monomers followed by a block of six guanine monomers), the diblock of A and G with the G-end grafted to the surface (G 6 A 6 , a block of six guanine monomers followed by a block of six adenine monomers), and an alternating sequence of A and G ((AG) 6 , one adenine monomer followed by a guanine monomer in an alternating manner) (Figure 1). Each of these chains contains 12 monomers at varying grafting densities and MgCl 2 concentrations, while keeping the solution temperature and background NaCl concentration fixed at 298 K and 10 mM, respectively. NaCl is added to the system to comply with the relevant experimental study by Holland et al. [34], which hypothesizes that Mg 2 + might replace Na + in a probable binding pathway.
The theoretical model is constructed using a Self-Consistent Field Theory (SCFT) approach for a single polyelectrolyte chain in a field of interacting species [8,35,36,37]. The polyelectrolyte chains are end-tethered to a surface and submerged in a salt and water bath, containing both NaCl and MgCl 2 salts. The motivation behind our study is to discover the extent of Mg 2 + ion binding to the polyelectrolytes and how it changes the structure and properties of the polyelectrolyte chains. In this molecular model, N p polyelectrolyte chains are end-grafted to a surface with cross-sectional area A. We assume the system to be homogeneous in the x and y directions and heterogeneous in the z direction, along the distance perpendicular to the grafting surface. Within the field theory framework, the heterogeneity normal to the surface is accounted for by discretizing the system space into a number of layers. The concentrations of the salts are converted to a density field to determine their individual contribution to the field, as a function of distance from the grafting surface. Cation binding to the polyelectrolytes is considered within the scope of reaction equilibrium calculations, rather than condensation near the charged monomers [38], with binding reactions relevant to the experimental study of a similar system [34]. We start constructing the model by calculating the total Helmholtz free energy of the system, which is given by,
F = T S c o n f T S m i x + F c h e m + F e l e c t + E r e p
Here, S c o n f is the conformational entropy of the grafted polymer chains, and S m i x is the mixing or translational entropy of all the free species: water (w), protons (H + ), hydroxyl ions (OH ), cations (Na + , Mg 2 + ), and anions (Cl , OH ). F c h e m is the free energy associated with the equilibrium reactions that the monomers undergo in this system. We have explicitly considered three equilibrium reactions for each monomer: (1) protonation and deprotonation via acid-base equilibrium reactions, (2) Mg 2 + binding, and (3) Na + binding. This yields a total of six equilibrium reactions accounted for in our model:
A + H + A H
G + H + G H
A + M g 2 + A M g +
G + M g 2 + G M g +
A + N a + A N a
G + N a + G N a
F e l e c t is the total electrostatic energy due to the charged species and E r e p is the repulsive interaction experienced between all species due to steric hindrance. T is the temperature of the system, which is held constant at 298 K. Expansion of all the energy and entropy terms gives the total Helmholtz free energy of the system,
β F A = σ p α P ( α ) ln P ( α ) + ρ w ( z ) ( ln ρ w ( z ) v w 1 ) d z + ρ H + ( z ) ( ln ρ H + ( z ) v w 1 + β μ H + 0 ) d z + ρ O H ( z ) ( ln ρ O H ( z ) v w 1 + β μ O H 0 ) d z + ρ N a + ( z ) ( ln ρ N a + ( z ) v w 1 + β μ N a + ) d z + ρ M g 2 + ( z ) ( ln ρ M g 2 + ( z ) v w 1 + β μ M g 2 + ) d z + ρ C l ( z ) ( ln ρ C l ( z ) v w 1 + β μ C l ) d z + β [ ρ q ( z ) ψ ( z ) + 1 2 ϵ w ( d ψ ( z ) d z ) 2 ] d z + ρ A ( z ) [ f A ( z ) ln f A ( z ) + f A H ( z ) ln f A H ( z ) + f A N a ( z ) ln f A N a ( z ) + f A M g + ( z ) ln f A M g + ( z ) + f A ( z ) β μ A 0 + f A H ( z ) β μ A H 0 + f A N a ( z ) β μ A N a 0 + f A M g + ( z ) β μ A M g + 0 ] d z + ρ G ( z ) [ f G ( z ) ln f G ( z ) + f G H ( z ) ln f G H ( z ) + f G N a ( z ) ln f G N a ( z ) + f G M g + ( z ) ln f G M g + ( z ) + f G ( z ) β μ G 0 + f G H ( z ) β μ G H 0 + f G N a ( z ) β μ G N a 0 + f G M g + ( z ) β μ G M g + 0 ] d z
The first term in Equation (8) stands for the structural or conformational energy ( T S c o n f ) of the polyelectrolyte chains, where P ( α ) is the probability of the polymer to be in a conformational state, α , and σ p = N p A is the grafting density on the surface. The second through seventh term represent the total mixing energy ( T S m i x ) of all the mobile species in the system. ρ i ( z ) is the position-dependent density of species i. μ 0 i and v i are the standard chemical potential and the volume of species i, respectively. The second term is the translational energy of the undissociated water molecules. The third and fourth terms represent the translational entropy of H + and OH ions, respectively, from the dissociation of water. The fifth, sixth, and seventh terms describe the translational entropy of the anions and cations (Mg 2 + , Na + , and Cl ) produced by the dissociation of the salts (MgCl 2 and NaCl). We assume that the salts in this system are completely dissociated. The eighth term represents the contribution of the electrostatic energy, F e l e c t . ψ ( z ) is the electrostatic potential, and ϵ w is the dielectric constant of water. ρ q ( z ) is the ensemble average charge density of the system, given by,
ρ q ( z ) = i = N a + , M g 2 + , C l , H + , O H e z i ρ i ( z ) e ρ A ( z ) ( f A ( z ) f A M g + ( z ) ) e ρ G ( z ) ( f G ( z ) f G M g + ( z ) )
Here, f i ( z ) is the fraction of a monomer in different states, such as deprotonated or charged (A , G ), protonated ( A H , G H ) or cation-bound (AMg + , GMg + ), etc., along the distance from the grafting surface.
The ninth and tenth terms are derived from the reaction equilibrium of Equations (2)–(5), respectively. The reaction equilibrium energies include the entropy of mixing between charged and uncharged groups, as well as the standard chemical potentials of these groups. Any volume change of the polyelectrolyte segments due to protonation and ion-binding is neglected. The system is in contact with a bath of anions, cations, protons, and hydroxyl ions and, therefore, best described by a semi-grand canonical ensemble. The proper thermodynamic potential (W) is given by,
W = F γ = w , N a + , M g 2 + , C l , H + , O H μ γ N γ
Here, N γ is the total number of molecules of species γ and μ γ is the exchange chemical potential of those species.

Extremization of the Free Energy

The free energy equation is extremized subjected to two constraints. One is the incompressibility constraint, originating from the repulsive interactions between different species,
ϕ p ( z ) + ϕ w ( z ) + ϕ N a + ( z ) + ϕ M g 2 + ( z ) + ϕ H + ( z ) + ϕ O H ( z ) + ϕ C l ( z ) = 1
Here, ϕ p ( z ) is the ensemble average polymer volume fraction, which is given by,
ϕ p ( z ) = ρ A ( z ) ( f A ( z ) v A + f A H ( z ) v A H + f A N a ( z ) v A N a + f A M g + ( z ) v A M g + ) + ρ G ( z ) ( f G ( z ) v G + f G H ( z ) v G H + f G N a ( z ) v G N a + f G M g + ( z ) v G M g + )
ρ i ( z ) is the average number density of monomer i at position z, expressed as,
ρ i = A , G ( z ) = σ p α P ( α ) n i = A , G ( α ; z ) v i = A , G
ϕ i ( z ) is the volume fraction of free species i in the system as a function of the distance from the grafting surface and defined as ϕ i ( z ) = ρ i v i .
The second constraint arises from the fact that the total number of each species of monomer is fixed. Hence, monomer fractions of different derivative species sum to unity,
f A ( z ) + f A H ( z ) + f A N a ( z ) + f A M g + ( z ) = 1
f G ( z ) + f G H ( z ) + f G N a ( z ) + f G M g + ( z ) = 1
Under these two constraints, the total thermodynamic potential per unit area is given by,
β W A = β F A k = w , C l , O H β μ k ρ k ( z ) d z β μ H + [ ρ H + ( z ) + f A H ( z ) ρ A ( z ) + f G H ( z ) ρ G ( z ) ] d z β μ N a + [ ρ N a + ( z ) + f A N a ( z ) ρ A ( z ) + f G N a ( z ) ρ G ( z ) ] d z β μ M g 2 + [ ρ M g 2 + ( z ) + f A M g + ( z ) ρ A ( z ) + f G M g + ( z ) ρ G ( z ) ] d z + β λ 1 ( z ) ρ A ( z ) ( f A ( z ) + f A H ( z ) + f A N a ( z ) + f A M g + ( z ) 1 ) d z + β λ 2 ( z ) ρ G ( z ) ( f G ( z ) + f G H ( z ) + f G N a ( z ) + f G M g + ( z ) 1 ) d z + β π ( z ) ( ϕ p ( z ) + ϕ w ( z ) + ϕ N a + ( z ) + ϕ M g 2 + ( z ) + ϕ H + ( z ) + ϕ O H ( z ) + ϕ C l ( z ) 1 ) d z
Here, π ( z ) , λ 1 ( z ) , and λ 2 ( z ) are Lagrange multipliers to incorporate the constraints into the free energy equation.
Extremization of the thermodynamic potential Equation (16), with respect to the relevant variables, gives the expressions for the equilibrium values of ρ i ( z ) , f i ( z ) , ψ ( z ) , and P ( α ) . From functional extremization, the density profiles of the free species are given by the following equations,
ρ w ( z ) w w = exp [ β π ( z ) v w ]
ρ H + ( z ) v w = exp [ β μ H + β μ H + 0 β π ( z ) v H + β ψ ( z ) e ]
ρ O H ( z ) v w = exp [ β μ O H β μ O H 0 β π ( z ) v O H + β ψ ( z ) e ]
ρ N a + ( z ) v w = exp [ β μ N a + β μ N a + 0 β π ( z ) v N a + β ψ ( z ) e ]
ρ M g 2 + ( z ) v w = exp [ β μ M g 2 + β μ M g 2 + 0 β π ( z ) v M g 2 + 2 β ψ ( z ) e ]
ρ C l ( z ) v w = exp [ β μ C l β μ C l 0 β π ( z ) v C l + β ψ ( z ) e ]
Functional extremization with respect to the monomer fractions, f i ( z ) , yields the governing equations for the chemical equilibrium reactions of both the monomers, A and G:
f A ( z ) f A H ( z ) = K A H 0 exp ( β π ( z ) Δ v A H ) ρ H + ( z ) v w
f A ( z ) f A N a ( z ) = K A N a 0 exp ( β π ( z ) Δ v A N a ) ρ N a + ( z ) v w
f A ( z ) f A M g + ( z ) = K A M g + 0 exp ( β π ( z ) Δ v A M g + ) ρ M g 2 + ( z ) v w
f G ( z ) f G H ( z ) = K G H 0 exp ( β π ( z ) Δ v G H ) ρ H + ( z ) v w
f G ( z ) f G N a ( z ) = K G N a 0 exp ( β π ( z ) Δ v G N a ) ρ N a + ( z ) v w
f G ( z ) f G M g + ( z ) = K G M g + 0 exp ( β π ( z ) Δ v G M g + ) ρ M g 2 + ( z ) v w
The quantity, K i 0 = exp ( β Δ G i 0 ) , corresponds to the chemical equilibrium constant that is derived from the standard chemical free energy, Δ G i 0 , of the respective formation reactions for A H , G H , A N a , G N a , A M g + , or G M g + . Δ v i denotes the volume change due to the reactions. The change in the standard free energy for the reaction, A + M g 2 + A M g + , is given by: Δ G A M g + 0 = μ A M g + 0 μ A 0 μ M g 2 + 0 , and the volume change of reaction is Δ v A M g + = v A M g + v A v M g 2 + . The reaction constants and the change in the volumes for the other reactions can be derived in a similar manner.
Extremization of the free energy with respect to the electrostatic potential yields the Poisson equation,
ϵ w d 2 ψ ( z ) d z 2 = ρ q ( z )
ϵ w d ψ ( z ) d z | z = 0 = 0 , lim r ψ ( z ) = 0
The probability distribution function (pdf) is derived from the functional minimization with P ( α ) ,
P ( α ) = 1 e exp [ n A ( α ; z ) v A ( ln f A ( z ) + β μ A 0 + β π ( z ) v A β e ψ ( z ) ) d z n G ( α ; z ) v G ( ln f G ( z ) + β μ G 0 + β π ( z ) v G β e ψ ( z ) ) ] d z
Equations (17)–(31) are solved simultaneously following the procedure described in the previous publications using this general approach [8,35,38]. These integro-differential equations are solved numerically by discretizing the space for a discretization length of 0.3 nm for 100 discrete layers. The solution of these sets of non-linear coupled equations yields the unknowns of the model, which are the Lagrange multiplier, π ( z ) , and the electrostatic potential, ψ ( z ) . The inputs necessary to solve the system of equations are the bulk concentrations of the salts, bulk pH, grafting density, volumes of different species, a set of polymer conformations, and the equilibrium reaction constants. The equilibrium constants for the binding reactions between the monomers (adenine and guanine) and the Mg 2 + cations are obtained from the binding free energies of the Second Harmonic Generation (SHG) and Atomic Force Microscopy (AFM) studies of Holland et al. [34], by using the previously-mentioned equation for equilibrium constant, K i 0 . The binding free energies are 32.1 KJ/mol and 35.6 KJ/mol for adenine and guanine, respectively. The pKa values for G and A are 1.6 and 3.5 , respectively [39]. The volume of Mg a q 2 + is 0.18 nm 3 , and the volumes for Na a q + and Cl a q are 0.05 nm 3 . We did not include a change in volume upon binding in this analysis. The set of polymer conformations is derived using a Rotational Isomeric State (RIS) model [40]. The RIS model enumerates the configurations by changing the dihedral angles between connected monomers along the chain. There are three possible configurations: trans, or gauche-plus, or gauche-minus. We then use solid body rotations of each finished configuration to provide us with several million physical configurations to use in the calculation.

3. Results and Discussions

To elucidate the effects of Mg 2 + binding on the structure and properties of a short chain surface-grafted A-G oligomer, we studied three different chain sequences: A-G diblock with the A-end grafted to the surface (A 6 G 6 ), A-G diblock with the G-end grafted to the surface (G 6 A 6 ), and A-G co-polymer with alternating A and G along the chain ((AG) 6 ) (referring to Figure 1). We studied the variation of the MgCl 2 concentration and the polymer surface coverage for the same chain types, as well as quantifying the number of bound Mg 2 + ions to each chain in different solution conditions, while keeping the background NaCl concentration and pH constant, at 10 mM and 7.0 , respectively. We present the quantitative results for all the chain sequences, and the qualitative results for the A-grafted chains only, as the results for the other chain systems do not deviate quantitatively from the A-grafted system. Our results provide valuable insight into the molecular details of this polyelectrolyte system in various biologically-relevant environments.

3.1. Effect of Sequence Heterogeneity on Mg 2 + Binding

The number of bound Mg 2 + per chain is calculated from the volume fraction of the polymer and the fraction of the Mg 2 + -bound monomers by using the following equation:
N M g 2 + / c h a i n = ϕ p ( z ) f P M g + ( z ) d z σ p v p
Figure 2a,b shows the number of bound Mg 2 + with A 6 G 6 , G 6 A 6 and (AG) 6 sequences at different grafting densities for 3 mM and 180 mM MgCl 2 content. The data in both plots demonstrate that, despite the sequence heterogeneity and variation on the grafting ends, all three sequences bind similar amounts of Mg 2 + ions at different grafting densities and salt concentrations, except for the 0.5 chains/nm 2 and 3 mM condition. The relatively low cation binding at this condition can be attributed to the unavailability of enough Mg 2 + ions at the low concentration of the MgCl 2 salt and high grafting density of the polymer chains. These results are in agreement with the experimental findings that Mg 2 + ions do not aggregate around the strongest binder (guanine), and are uniformly distributed throughout the length of the chains [33]. This agrees with the notion that non-specific interactions due to the electrostatics of the sugar-phosphate backbone dominating the ion binding to the ssDNA chains, as the non-specific part of the free energy for adenine and guanine ( 21 KJ/mol, from Holland et al. [34]) is much higher than the specific part of the free energies ( 11.1 KJ/mol and 14.6 KJ/mol, respectively, from Holland et al. [34]), and the specific free energies are all within a few KJ of each other, hence making the ion binding less distinctive.

3.2. Effect of Ionic Strength and Grafting Density

Figure 3a–c shows the volume fraction profiles of the A-grafted chains (A 6 G 6 ), along the distance from the grafting surface for low ( 0.005 chains/nm 2 ), medium ( 0.05 chains/nm 2 ), and high ( 0.5 chains/nm 2 ) grafting densities, for a range of MgCl 2 concentrations. At low grafting density (Figure 3a), as we increase the MgCl 2 salt concentration from 3 mM–180 mM, the chain structures contract towards the grafting surface, and we get a distinct peak of the highest volume fractions at a distance close to the surface. This change in the chain structure is due to the reduction in negative charge of the polyelectrolytes by binding to the Mg 2 + ions that lowers the repulsive interaction between the monomers.
Although at first glance, the collapse does not seem to be very prominent, it is significant for the length scales of the polyelectrolytes in our system (chain length of only 12 monomers). However, for higher grafting density (Figure 3b), the peak is less distinctive at 3 mM and 50 mM MgCl 2 concentrations than 180 mM and nearly plateaus down from 0.8 nm–2 nm. Herein, the high grafting density creates steric repulsion, and the system faces a competition between charge screening and steric repulsion to stabilize the polyelectrolyte structures. At the low MgCl 2 salt concentrations (3 mM and 50 mM), the system energetically favors chain stretching to accommodate both charge repulsion and steric hindrance. However, when we further increase the MgCl 2 salt concentration to 180 mM, charge screening by the Mg 2 + ion dominates over the steric repulsion, and we get a distinct chain collapse near the grafting surface. Figure 3c corresponds to the volume fraction profiles of the polyelectrolyte chains at high grafting density ( 0.5 chains/nm 2 ) and shows a clearly opposite picture than Figure 3a,b, with the chains stretching, while we increase the salt concentration. In this case, at 3 mM MgCl 2 , the number of the bound Mg 2 + ions is sufficient for mitigating the charge repulsion inside the brush. A further increase in the salt concentration induces additional steric hindrance due to the volume exclusion, and the chains pay in conformational entropy to stretch the chains and accommodate more Mg 2 + ions inside the brush. This phenomenon of polyelectrolyte chain stretching with increasing counterion concentration is well known as the “re-entrant phenomena” and has been experimentally verified by several studies [41,42].
Figure 4 represents the profiles of the deprotonated (negatively-charged) polymer fraction ( f P ( z ) ). This quantity is calculated from the individual monomer-species fractions by,
f P ( z ) = f A ( z ) ϕ A ( z ) + f G ( z ) ϕ G ( z ) ϕ p ( z )
In Figure 4a, at the 0.005 chains/nm 2 grafting density, almost 100 % of the polymers are bound to Mg 2 + for 50 mM and 180 mM salts, and even for 3 mM salt, a very small fraction is deprotonated ( f P ( z ) 0.0001 ), with the rest of the polymers being Mg 2 + -bound. At 0.05 chains/nm 2 (Figure 4b), the polymers are still (nearly) entirely bound to Mg 2 + at 50 mM and 180 mM salt concentrations. However, for 3 mM salt, the fraction of the negatively-charged polymers slightly increases.
As we further increase the grafting density up to 0.5 chains/nm 2 in Figure 4c, the ionic strength shows a more prominent effect on the chain structure and chemistry. Now, at 3 mM MgCl 2 , about 50 % of the polymers are negatively charged. The increase of the ionic strength to 50 mM and 180 mM decreases the negative charge by binding to more Mg 2 + to reduce the negative charge inside the brush.
The “re-entrant behavior” shown in Figure 3c can be further explained using Figure 4c and Figure 5. At 3 mM MgCl 2 , about half of the polymers are negatively charged, and another half is positively charged due to the formation of A M g + and G M g + complexes with the monomers. Hence, the brush environment is nearly neutralized with low or no residual positive or negative charges, and the brush becomes collapsed for that grafting density. Therefore, increasing the salt concentration does not contribute to the global charge neutralization inside the brush; rather, the availability of more Mg 2 + ions, that have very high affinity towards binding to the negatively-charged monomers, creates high steric hindrance. As a result, we see stretching of the chains, and the system pays in conformational entropy to reduce steric repulsion inside the brush.
Figure 5 represents the volume fractions of free Mg 2 + ions inside the brush and the surrounding medium. At low grafting density ( 0.005 chains/nm 2 , Figure 5a), there are significant amounts of free Mg 2 + ions inside the brush at all salt concentrations, which is depicted by the small deviation in the free Mg 2 + volume fractions from the bulk value.
Here, along with the ion binding, localization of the free Mg 2 + ions screens the negative charges inside the brush, which results in a collapse of the polyelectrolyte chains (referring to Figure 3a). However, as we increase the grafting density to 0.05 chains/nm 2 (Figure 5b), the amount of the free Mg 2 + ions inside the brush drops noticeably from the bulk value, compared to the low grafting density brush. However, the system still entropically favors localization of free Mg 2 + ions while staying at a collapsed state (Figure 3b). A further increase in the grafting density up to 0.5 chains/nm 2 gives rise to the above-mentioned re-entrant phenomena (Figure 3c). At this point, the brush region acts like a barrier to the free Mg 2 + ions, and the brush is completely devoid of the free Mg 2 + ions at all the salt concentrations (Figure 5c). For a grafting density this high, the brush is so dense that the binding to the Mg 2 + ions increases the volume of the monomers, which also increases the steric hindrance, even at a low MgCl 2 salt concentration. Increasing the MgCl 2 salt concentration promotes binding to more Mg 2 + ions, which is energetically more favorable than bringing the free ions from the bulk and localizing them inside the brush. Hence, the system pays in conformational entropy to stretch the chains to avoid high steric repulsion in the brush, and the re-entrant phenomena arises.
Figure 6a,b presents the change in local pH inside the brush along the distance from the grafting surface for low ( 0.005 chains/nm 2 ) and high ( 0.5 chains/nm 2 ) grafting densities, respectively, for varying ionic strengths. The pH profiles, when no MgCl 2 was added in the system, are included in the figures to compare the system’s response with the added MgCl 2 salt. At no added MgCl 2 (0 mM), the system only had 10 mM NaCl as a background electrolyte to modulate charge regulation inside the brush. At the low grafting density, the pH inside the brush (Figure 6a) only deviates by ∼0.5 from the bulk. However, for the high grafting density, the available Na + ions are not enough to reduce the negative charges on the polyelectrolytes, and the local pH is much lower with the lowest value of ∼4.5 (Figure 6b). When we add MgCl 2 to the system, even as low as 3 mM, the negatively-charged monomers readily bind to the Mg 2 + to form A M g + and G M g + complexes that are positively charged. The production of the positively-charged complexes creates a significant change in the pH inside the brush for all the grafting densities. For the low grafting density (Figure 6a), at 3 mM MgCl 2 , the pH change is only for one unit. As the MgCl 2 concentration is increased, the pH surge subsides, and at 180 mM, the charge of the system is almost neutral. However, for a high grafting density (Figure 6b), the addition of only 3 mM MgCl 2 salt creates a pH increase of more than four units (from pH = 4.5 to pH = 9.0 ). As the MgCl 2 content is increased, unlike the lower grafting density case, the pH further increases up to ∼9.5. This dramatic change in the local pH with the increase in MgCl 2 is consistent with the change in the free chloride volume fraction profiles in Figure 7 and Figure 8.
In Figure 7, volume fractions of free chloride ions along the distance from the grafting surface are reported at a low grafting density of the brush for varying MgCl 2 concentrations. When there is no MgCl 2 in the system (Figure 7a), the free chloride ions from the dissociation of NaCl salt are excluded from the system and reside in the bulk, due to the repulsion of highly negatively-charged monomers inside the brush. However, when the MgCl 2 is added, the positively-charged monomer-cation complex takes the place of the deprotonated monomers. At this stage, the free chloride ions contribute as counterions to minimize the repulsion due to the positively-charged species inside the brush. At 50 mM and 180 mM MgCl 2 (Figure 7c,d, respectively), the free chloride ion concentration inside the brush is significantly higher than the bulk, and the local pH is close to neutral. However, for the MgCl 2 salt concentration as low as 3 mM, the available chloride ions are not enough to neutralize the positive charge in the brush, and the pH is highest for the specific grafting density.
Figure 8 presents the free chloride volume fractions inside the chain as a function of distance from the grafting surface at high grafting density ( 0.5 chains/nm 2 ). In the absence of MgCl 2 , almost all the chloride ions are excluded from the system due to the repulsion by the highly negatively-charged chains, and the counterions (Na + ) are excluded due to the high steric hindrance of the crowded brush. As we add MgCl 2 salt in the system, the formation of the positively-charged complex now gives rise to electrostatic repulsion inside the brush. It would be thermodynamically favorable for the system to employ negatively-charged chloride counterions to mitigate the electrostatic repulsion. However, the system is already highly dense with the added Mg 2 + ions bound to the polymers, and further accommodation of any more species creates high steric hindrance. Hence, the system prefers to stay in a positively-charged condition inside the brush, with a maximum pH of 10.0 , and the free chloride ions mostly reside in the outer periphery of the brush, with high volume fractions near the bulk (Figure 8b–d).
Figure 9 presents the water volume fraction profiles around the oligomer chains at different grafting densities and varying MgCl 2 salt concentrations. These plots depict the water displacement around the oligomer chains due to the Mg 2 + ion binding with increasing MgCl 2 salt concentration. In Figure 9a, at low grafting density ( 0.005 chains/nm 2 ) and at all the MgCl 2 concentrations, the water volume fractions are nearly uniform throughout the polymer layer and the bulk, with a small amount of water displacement near the grafting surface. The water displacement is higher at the high salt concentrations (50 mM and 180 mM MgCl 2 ) compared to when no MgCl 2 salt is added. At a medium grafting density ( 0.05 chains/nm 2 ), the water volume fractions inside the polymer chains varies significantly compared to the bulk at all of the concentrations of the MgCl 2 salt. Herein, to accommodate more Mg 2 + ions inside the brush, more water is being displaced. Similar to Figure 9a,b, Figure 9c shows that more water is displaced with increasing MgCl 2 salt concentration at a high grafting density ( 0.5 chains/nm 2 ), and the difference in the water volume fraction is very high between the bulk and inside the chain. The high grafting density of the oligomer chains creates high steric repulsion inside the brush that lowers the water content around the oligomers.
Figure 10 depicts the volume fraction profiles of free Na + ions as a function of distance from the grafting surface at varying grafting densities with the presence of varying concentrations of MgCl 2 salt. When there is no MgCl 2 salt present (the blue lines in the plots in Figure 10a–c), the NaCl salt in the system solely contributes to the charge neutralization inside the oligomer brushes, and hence, the volume fractions of the Na + ions in the brush layer are very high compared to the bulk. However, as we add the MgCl 2 salt in our system, the divalent Mg 2 + ions displace the Na + ions, and we see a significant reduction in the volume fractions of the Na + ions inside the brush at all the grafting densities and all the concentrations of the MgCl 2 salt. At the low grafting density (Figure 10a), we still see a considerable amount of the Na + ions inside the chain region, due to low steric repulsion inside the brush. Nevertheless, for higher grafting densities of the chains (Figure 10b,c), almost all the Na + ions are replaced by the divalent Mg 2 + cations at all the MgCl 2 salt concentrations. These results theoretically support the hypothesis from the relevant experimental study by Holland et al. [34], which states that Mg 2 + displaces Na + through a possible binding pathway with surface-grafted ssDNA oligomers.
In summary, the results in this section portray the complex and coupled interplay between the grafting density of the chains and the ionic strength to govern the structure, ion binding, and local environment inside the brush. The system generates the thermodynamic equilibrium state by adjusting between the conformational entropy, electrostatic potential, and repulsive interactions.
It is important to note that this work considers the Mg 2 + ion explicitly, without taking into account the solvation effect. Mg 2 + has a strong hydration shell compared to the bulk and binds to the nucleotides via a solvent-mediated-ion pairing [43]. However, to capture the thermodynamics of Mg 2 + -nucleotide interactions accurately, it is necessary to treat Mg 2 + explicitly [26]. This work provides a basis for further theoretical study involving the Mg 2 + -solvation effect on the tethered ssDNA-Mg 2 + binding. Furthermore, this work does not include the possibility of secondary or tertiary structure formation for these particular sequences of nucleic acid oligomers, as no structure formation is evidenced by a relevant experimental study [33]. However, the possibility of secondary or tertiary structure formation will be included in our future studies of these polyelectrolyte systems.

4. Conclusions

Divalent metal ion binding to surface-grafted nucleic acid oligomers is investigated by studying the effects of the ionic strength and grafting density on the oligomer structure and chemistry, with a field theoretic molecular model. The cation binding is explicitly included in the model by utilizing the experimentally-derived binding free energies of the relevant reactions. Quantitative assessment of the ion cloud around the oligomers shows a uniform distribution of the ions around different sequences and reinforces the dominance of non-specific electrostatic attraction between the nucleobases and the cations as the driving force for the cation binding [33,44]. The analysis of the system with the variation in the ionic strength and polymer grafting density shows a complex coupling between the chain conformation and the ion cloud to maintain the stability of the system by achieving the minimum energy state. At low grafting density, when the polymers are sparsely grafted, cation binding and ion condensation around the charged oligomers leads to charge neutralization inside the brush, which is accompanied by a chain collapse. At high grafting density, however, cation binding results in the reversal of the oligomer charge, which can no longer be neutralized by the anions due to anion exclusion from the brush to avoid steric repulsion, and hence, we get a highly-stretched polymer brush. Our results also show that the ionic strength has a more prominent effect on the structure and properties of the oligomer brushes when they are densely grafted, compared to their sparsely-grafted counterpart. In its current state, this model can serve as a foundation for field theoretic studies of more complex systems to dissect the ion binding scenario around aptamers and single-stranded nucleic acids. Although, originally constructed to mimic surface-anchored nucleic acid aptamers for the robust design of aptamer-based biosensors and therapeutics, this molecular model can also be employed to understand the molecular level interactions of other natural or synthetic polyelectrolytes with metal ions in a solution environment, for applications ranging from colloid chemistry to drug design for controlled release.

Author Contributions

Conceptualization, M.J.U.; Methodology, M.J. and M.J.U.; Software, M.J.U.; Validation, M.J. and M.J.U; Formal Analysis, M.J. and M.J.U.; Investigation, M.J.; Resources, M.J.U.; Data Curation, M.J.; Writing-Original Draft Preparation, M.J.; Writing-Review & Editing, M.J. and M.J.U.; Visualization, M.J.; Supervision, M.J.U.; Project Administration, M.J.U.; Funding Acquisition, M.J.U.

Funding

This research was funded by National Science Foundation Grant NSF-DMR-1609939.

Acknowledgments

The authors acknowledge Rikkert Nap of Northwestern University for useful discussions and Celeste Kennard for a critical reading of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keefe, A.D.; Pai, S.; Ellington, A. Aptamers as therapeutics. Nat. Rev. Drug Discov. 2010, 9, 537. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, R.; Li, W.; Cai, T.; Deng, Y.; Ding, Z.; Liu, Y.; Zhu, X.; Wang, X.; Liu, J.; Liang, B. TiO2 Nanolayer-Enhanced Fluorescence for Simultaneous Multiplex Mycotoxin Detection by Aptamer Microarrays on a Porous Silicon Surface. ACS Appl. Mater. Interfaces 2018, 10, 14447–14453. [Google Scholar] [CrossRef] [PubMed]
  3. Cho, E.J.; Lee, J.-W.; Ellington, A.D. Applications of aptamers as sensors. Annu. Rev. Anal. Chem. 2009, 2, 241–264. [Google Scholar] [CrossRef] [PubMed]
  4. Dua, P.; Kang, S.; Shin, H.-S.; Kim, S.; Lee, D.-K. Cell-SELEX-Based Identification of a Human and Mouse Cross-Reactive Endothelial Cell-Internalizing Aptamer. Nucleic Acid Ther. 2018, 28. [Google Scholar] [CrossRef] [PubMed]
  5. Zhou, G.; Latchoumanin, O.; Hebbard, L.; Duan, W.; Liddle, C.; George, J.; Qiao, L. Aptamers as targeting ligands and therapeutic molecules for overcoming drug resistance in cancers. Adv. Drug Deliv. Rev. 2018, 134, 107–121. [Google Scholar] [CrossRef] [PubMed]
  6. Foster, A.; DeRosa, M.C. Development of a biocompatible layer-by-layer film system using aptamer technology for smart material applications. Polymers 2014, 6, 1631–1654. [Google Scholar] [CrossRef]
  7. Nakatsuka, N.; Yang, K.-A.; Abendroth, J.M.; Cheung, K.M.; Xu, X.; Yang, H.; Zhao, C.; Zhu, B.; Rim, Y.S.; Yang, Y. Aptamer–field-effect transistors overcome Debye length limitations for small-molecule sensing. Science 2018, 362, 319–324. [Google Scholar] [CrossRef] [PubMed]
  8. Uline, M.J.; Rabin, Y.; Szleifer, I. Effects of the salt concentration on charge regulation in tethered polyacid monolayers. Langmuir 2011, 27, 4679–4689. [Google Scholar] [CrossRef] [PubMed]
  9. Nap, R.J.; Szleifer, I. Effect of calcium ions on the interactions between surfaces end-grafted with weak polyelectrolytes. J. Chem. Phys. 2018, 149, 163309. [Google Scholar] [CrossRef] [PubMed]
  10. Nap, R.J.; Solveyra, E.G.; Szleifer, I. The interplay of nanointerface curvature and calcium binding in weak polyelectrolyte-coated nanoparticles. Biomater. Sci. 2018, 6, 1048–1058. [Google Scholar] [CrossRef] [PubMed]
  11. Zwanikken, J.W.; Guo, P.; Mirkin, C.A.; Olvera de la Cruz, M. Local ionic environment around polyvalent nucleic acid-functionalized nanoparticles. J. Phys. Chem. C 2011, 115, 16368–16373. [Google Scholar] [CrossRef]
  12. Kumar, R.; Sumpter, B.G.; Kilbey, S.M. Charge regulation and local dielectric function in planar polyelectrolyte brushes. J. Chem. Phys. 2012, 136, 234901. [Google Scholar] [CrossRef] [PubMed]
  13. Lewis, T.; Pandav, G.; Omar, A.; Ganesan, V. Complexation between weakly basic dendrimers and linear polyelectrolytes: effects of grafts, chain stiffness, and pOH. Soft Matter 2013, 9, 6955–6969. [Google Scholar] [CrossRef]
  14. Liu, J.; Cao, Z.; Lu, Y. Functional nucleic acid sensors. Chem. Rev. 2009, 109, 1948–1998. [Google Scholar] [CrossRef] [PubMed]
  15. Zhou, W.; Huang, P.-J.J.; Ding, J.; Liu, J. Aptamer-based biosensors for biomedical diagnostics. Analyst 2014, 139, 2627–2640. [Google Scholar] [CrossRef] [PubMed]
  16. Lee, J.-S.; Han, M.S.; Mirkin, C.A. Colorimetric detection of mercuric ion (Hg2+) in aqueous media using DNA-functionalized gold nanoparticles. Angew. Chem. Int. Ed. 2007, 46, 4093–4096. [Google Scholar] [CrossRef] [PubMed]
  17. Saenger, W. Principles of Nucleic Acid Structure; Springer: Berlin, Germany, 1984; pp. 201–219. [Google Scholar]
  18. Reshetnikov, R.V.; Sponer, J.; Rassokhina, O.I.; Kopylov, A.M.; Tsvetkov, P.O.; Makarov, A.A.; Golovin, A.V. Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process. Nucleic Acids Res. 2011, 39, 9789–9802. [Google Scholar] [CrossRef] [PubMed]
  19. Ono, A.; Torigoe, H.; Tanaka, Y.; Okamoto, I. Binding of metal ions by pyrimidine base pairs in DNA duplexes. Chem. Soc. Rev. 2011, 40, 5855–5866. [Google Scholar] [CrossRef] [PubMed]
  20. Anastassopoulou, J.; Theophanides, T. Magnesium–DNA interactions and the possible relation of magnesium to carcinogenesis. Irradiation and free radicals. Critical Rev. Oncol./Hematol. 2002, 42, 79–91. [Google Scholar] [CrossRef]
  21. Pascal, A.; Grover, N.; Westhof, E. Metal ion binding to RNA. Met. Ions Life Sci. 2011, 9, 1–35. [Google Scholar]
  22. Mocci, F.; Laaksonen, A. Insight into nucleic acid counterion interactions from inside molecular dynamics simulations is ‘worth its salt’. Soft Matter 2012, 8, 9268–9284. [Google Scholar] [CrossRef]
  23. Murthy, C.S.; Bacquet, R.J.; Rossky, P.J. Ionic distributions near polyelectrolytes. A comparison of theoretical approaches. J. Phys. Chem. 1985, 89, 701–710. [Google Scholar] [CrossRef]
  24. Patra, C.N.; Yethiraj, A. Density functional theory for the distribution of small ions around polyions. J. Phys. Chem. B 1999, 103, 6080–6087. [Google Scholar] [CrossRef]
  25. Howard, J.J.; Lynch, G.C.; Pettitt, B.M. Ion and solvent density distributions around canonical B-DNA from integral equations. J. Phys. Chem. B 2010, 115, 547–556. [Google Scholar] [CrossRef] [PubMed]
  26. Hayes, R.L.; Noel, J.K.; Whitford, P.C.; Mohanty, U.; Sanbonmatsu, K.Y.; Onuchic, J.N. Reduced model captures Mg2+-RNA interaction free energy of riboswitches. Biophys. J. 2014, 106, 1508–1519. [Google Scholar] [CrossRef] [PubMed]
  27. Savelyev, A.; MacKerell, A.D., Jr. Competition among Li+, Na+, K+ and Rb+ monovalent ions for DNA in molecular dynamics simulations using the additive CHARMM36 and Drude polarizable force fields. J. Phys. Chem. B 2015, 119, 4428–4440. [Google Scholar] [CrossRef] [PubMed]
  28. Li, W.; Nordenskiöld, L.; Mu, Y. Sequence-Specific Mg2+ –DNA Interactions: A Molecular Dynamics Simulation Study. J. Phys. Chem. B 2011, 115, 14713–14720. [Google Scholar] [CrossRef] [PubMed]
  29. Gebala, M.; Bonilla, S.; Bisaria, N.; Herschlag, D. Does cation size affect occupancy and electrostatic screening of the nucleic acid ion atmosphere? J. Am. Chem. Soc. 2016, 138, 10925–10934. [Google Scholar] [CrossRef] [PubMed]
  30. Gebala, M.; Giambasu, G.M.; Lipfert, J.; Bisaria, N.; Bonilla, S.; Li, G.; York, D.M.; Herschlag, D. Cation–anion interactions within the nucleic acid ion atmosphere revealed by ion counting. J. Am. Chem. Soc. 2015, 137, 14705–14715. [Google Scholar] [CrossRef] [PubMed]
  31. Mills, P.; Anderson, C.F.; Record, M.T., Jr. Monte Carlo studies of counterion-DNA interactions. Comparison of the radial distribution of counterions with predictions of other polyelectrolyte theories. J. Phys. Chem. 1985, 89, 3984–3994. [Google Scholar] [CrossRef]
  32. Jacobson, D.R.; Saleh, O.A. Counting the ions surrounding nucleic acids. Nucleic Acids Res. 2016, 45, 1596–1605. [Google Scholar] [CrossRef] [PubMed]
  33. Holland, J.G.; Malin, J.N.; Jordan, D.S.; Geiger, F.M. Specific and Nonspecific Metal Ion- Nucleotide Interactions at Aqueous/Solid Interfaces Functionalized with Adenine, Thymine, Guanine, and Cytosine Oligomers. J. Am. Chem. Soc. 2011, 133, 2567–2570. [Google Scholar] [CrossRef] [PubMed]
  34. Holland, J.G.; Jordan, D.S.; Geiger, F.M. Divalent metal cation speciation and binding to surface-bound oligonucleotide single strands studied by second harmonic generation. J. Phys. Chem. B 2011, 115, 8338–8345. [Google Scholar] [CrossRef] [PubMed]
  35. Szleifer, I.; Carignano, M. Tethered polymer layers. Adv. Chem. Phys. Polym. Syst. 1996, 94, 165–260. [Google Scholar]
  36. Uline, M.J.; Meng, S.; Szleifer, I. Surfactant driven surface anchoring transitions in liquid crystal thin films. Soft Matter 2010, 6, 5482–5490. [Google Scholar] [CrossRef]
  37. Van der Munnik, N.P.; Sajib, M.S.J.; Moss, M.A.; Wei, T.; Uline, M.J. Determining the Potential of Mean Force for Amyloid-β Dimerization: Combining Self-Consistent Field Theory with Molecular Dynamics Simulation. J. Chem. Theory Comput. 2018, 14, 2696–2704. [Google Scholar] [CrossRef] [PubMed]
  38. Nap, R.J.; Park, S.H.; Szleifer, I. Competitive calcium ion binding to end-tethered weak polyelectrolytes. Soft Matter 2018, 14, 2365–2378. [Google Scholar] [CrossRef] [PubMed]
  39. Bloomfield, V.; Crothers, D.M. Nucleic Acids: Structures, Properties and Functions; No. 574.192 B52; Universidad Autónoma Chapingo: Texcoco, Mexico, 2000. [Google Scholar]
  40. Flory, P.J.; Volkenstein, M. Statistical Mechanics of Chain Molecules; Interscience: New York, NY, USA, 1969. [Google Scholar]
  41. Yu, J.; Mao, J.; Yuan, G.; Satija, S.; Jiang, Z.; Chen, W.; Tirrell, M. Structure of polyelectrolyte brushes in the presence of multivalent counterions. Macromolecules 2016, 49, 5609–5617. [Google Scholar] [CrossRef]
  42. Wu, T.; Gong, P.; Szleifer, I.; Vlcěk, P.; Šubr, V.; Genzer, J. Behavior of surface-anchored poly (acrylic acid) brushes with grafting density gradients on solid substrates: 1. Experiment. Macromolecules 2007, 40, 8756–8764. [Google Scholar] [CrossRef]
  43. Lee, Y.; Thirumalai, D.; Hyeon, C. Ultrasensitivity of water exchange kinetics to the size of metal ion. J. Am. Chem. Soc. 2017, 139, 12334–12337. [Google Scholar] [CrossRef] [PubMed]
  44. Bai, Y.; Greenfeld, M.; Travers, K.J.; Chu, V.B.; Lipfert, J.; Doniach, S.; Herschlag, D. Quantitative and comprehensive decomposition of the ion atmosphere around nucleic acids. J. Am. Chem. Soc. 2007, 129, 14981–14988. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic representation of an end-grafted polymer in a salt solution environment (left) and chain sequences used for molecular modeling (right).
Figure 1. Schematic representation of an end-grafted polymer in a salt solution environment (left) and chain sequences used for molecular modeling (right).
Polymers 10 01403 g001
Figure 2. No. of bound Mg 2 + to different ssDNA sequences at varying grafting densities for (a) 3 mM MgCl 2 and (b) 180 mM MgCl 2 . The color bars correspond to the A-grafted chain (blue), G-grafted chain (yellow), and A-G alternate chain (red).
Figure 2. No. of bound Mg 2 + to different ssDNA sequences at varying grafting densities for (a) 3 mM MgCl 2 and (b) 180 mM MgCl 2 . The color bars correspond to the A-grafted chain (blue), G-grafted chain (yellow), and A-G alternate chain (red).
Polymers 10 01403 g002
Figure 3. Total polymer volume fraction profiles as a function of distance from the grafting surface at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 . Blue lines correspond to 3 mM MgCl 2 , yellow lines to 50 mM MgCl 2 , and red lines to 180 mM MgCl 2 .
Figure 3. Total polymer volume fraction profiles as a function of distance from the grafting surface at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 . Blue lines correspond to 3 mM MgCl 2 , yellow lines to 50 mM MgCl 2 , and red lines to 180 mM MgCl 2 .
Polymers 10 01403 g003
Figure 4. Deprotonated polymer fraction profiles at 3 mM (blue lines), 50 mM (yellow lines), and 180 mM (red lines) MgCl 2 concentrations. (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 .
Figure 4. Deprotonated polymer fraction profiles at 3 mM (blue lines), 50 mM (yellow lines), and 180 mM (red lines) MgCl 2 concentrations. (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 .
Polymers 10 01403 g004
Figure 5. Free Mg 2 + volume fraction profiles as a function of distance from the grafting surface at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 grafting densities for 3 mM (blue lines), 50 mM (Yellow lines), and 180 mM (green lines) MgCl 2 concentrations.
Figure 5. Free Mg 2 + volume fraction profiles as a function of distance from the grafting surface at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 grafting densities for 3 mM (blue lines), 50 mM (Yellow lines), and 180 mM (green lines) MgCl 2 concentrations.
Polymers 10 01403 g005
Figure 6. pH profiles along the distance from the grafting surface at (a) 0.005 chains/nm 2 and (b) 0.5 chains/nm 2 for 0 mM (blue lines), 3 mM (yellow lines), 50 mM (red lines), and 180 mM (green lines) MgCl 2 concentrations.
Figure 6. pH profiles along the distance from the grafting surface at (a) 0.005 chains/nm 2 and (b) 0.5 chains/nm 2 for 0 mM (blue lines), 3 mM (yellow lines), 50 mM (red lines), and 180 mM (green lines) MgCl 2 concentrations.
Polymers 10 01403 g006
Figure 7. Chloride volume fractions at 0.005 chains/nm 2 for 0 mM (a), 3 mM (b), 50 mM (c), and 180 mM (d) MgCl 2 concentrations.
Figure 7. Chloride volume fractions at 0.005 chains/nm 2 for 0 mM (a), 3 mM (b), 50 mM (c), and 180 mM (d) MgCl 2 concentrations.
Polymers 10 01403 g007
Figure 8. Chloride volume fractions at 0.5 chains/nm 2 for 0 mM (a), 3 mM (b), 50 mM (c), and 180 mM (d) MgCl 2 concentrations.
Figure 8. Chloride volume fractions at 0.5 chains/nm 2 for 0 mM (a), 3 mM (b), 50 mM (c), and 180 mM (d) MgCl 2 concentrations.
Polymers 10 01403 g008
Figure 9. Water volume fraction profiles as a function of distance from the grafting surface at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 grafting densities for 0 mM (blue lines), 50 mM (green lines), and 180 mM (yellow lines) MgCl 2 salt concentrations. Note the different scales of the plots (a,b,c).
Figure 9. Water volume fraction profiles as a function of distance from the grafting surface at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 grafting densities for 0 mM (blue lines), 50 mM (green lines), and 180 mM (yellow lines) MgCl 2 salt concentrations. Note the different scales of the plots (a,b,c).
Polymers 10 01403 g009
Figure 10. Volume fraction profiles of free Na + ions as a function of distance from the grafting surface in the presence of 10 mM NaCl salt and at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 grafting densities for 0 mM (blue lines), 3 mM (red lines) 50 mM (green lines), and 180 mM (yellow lines) MgCl 2 salt concentrations. Note the different scales of the plots (a,b,c).
Figure 10. Volume fraction profiles of free Na + ions as a function of distance from the grafting surface in the presence of 10 mM NaCl salt and at (a) 0.005 chains/nm 2 , (b) 0.05 chains/nm 2 , and (c) 0.5 chains/nm 2 grafting densities for 0 mM (blue lines), 3 mM (red lines) 50 mM (green lines), and 180 mM (yellow lines) MgCl 2 salt concentrations. Note the different scales of the plots (a,b,c).
Polymers 10 01403 g010

Share and Cite

MDPI and ACS Style

Jahan, M.; Uline, M.J. Quantifying Mg2+ Binding to ssDNA Oligomers: A Self-Consistent Field Theory Study at Varying Ionic Strengths and Grafting Densities. Polymers 2018, 10, 1403. https://doi.org/10.3390/polym10121403

AMA Style

Jahan M, Uline MJ. Quantifying Mg2+ Binding to ssDNA Oligomers: A Self-Consistent Field Theory Study at Varying Ionic Strengths and Grafting Densities. Polymers. 2018; 10(12):1403. https://doi.org/10.3390/polym10121403

Chicago/Turabian Style

Jahan, Merina, and Mark J. Uline. 2018. "Quantifying Mg2+ Binding to ssDNA Oligomers: A Self-Consistent Field Theory Study at Varying Ionic Strengths and Grafting Densities" Polymers 10, no. 12: 1403. https://doi.org/10.3390/polym10121403

APA Style

Jahan, M., & Uline, M. J. (2018). Quantifying Mg2+ Binding to ssDNA Oligomers: A Self-Consistent Field Theory Study at Varying Ionic Strengths and Grafting Densities. Polymers, 10(12), 1403. https://doi.org/10.3390/polym10121403

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop