DNA Binding to the Silica Surface - The Journal of Physical Chemistry

May 12, 2015 - Eric R. Castro , Mark D. Tarn , Pavlína Ginterová , Hanliang Zhu , Ying Xu ... Bobo Shi , Yun Kyung Shin , Ali A. Hassanali , Sherwin...
0 downloads 0 Views 5MB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

DNA Binding to the Silica Surface Bobo Shi, Yun Kyung Shin, Ali A. Hassanali, and Sherwin J Singer J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.5b01983 • Publication Date (Web): 12 May 2015 Downloaded from http://pubs.acs.org on May 14, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

DNA Binding to the Silica Surface Bobo Shi‡ , Yun Kyung Shin†,a , Ali A. Hassanali‡,b , Sherwin J. Singer†,‡,∗ †Department of Chemistry and Biochemistry and ‡Biophysics Program Ohio State University, Columbus, OH 43210 May 10, 2015

Abstract We investigate the DNA-silica binding mechanism using molecular dynamics simulations. This system is of technological importance, and also of interest to explore how negatively charged DNA can bind to a silica surface, which is also negatively charged at pH values above its isoelectric point near pH 3. We find that the two major binding mechanisms are attractive interactions between DNA phosphate and surface silanol groups and hydrophobic bonding between DNA base and silica hydrophobic region. Umbrella sampling and the Weighted Histogram Analysis Method (WHAM) are used to calculate the free energy surface for detachment of DNA from a binding configuration to a location far from the silica surface. Several factors explain why single-stranded DNA (ssDNA) has been observed to be more strongly attracted to silica than double-stranded (dsDNA): 1) ssDNA is more flexible and therefore able to maximize the number of binding interactions. 2) ssDNA has free unpaired bases to form hydrophobic attachment to silica while dsDNA has to break hydrogen bonds with base partners to get free bases. 3) The linear charge density of dsDNA is twice that of ssDNA. We devise a procedure to approximate the atomic forces between biomolecules and amorphous silica to enable large-scale biomolecule-silica simulations as reported here. Keywords: silica, DNA adsorption, umbrella sampling, simulation

* Corresponding author: Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH 43210 614-292-8909 [email protected] Current address: Department of Mechanical and Nuclear Engineering, Pennsylvania State University, University Park PA 16802 b Current address: International Center for Theoretical Physics, Condensed Matter and Statistical Physics 34151, Trieste, Italy a

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1 Introduction Currently there is great interest in the development of micro- and nano-scale devices for analysis and delivery of bio-molecules.1–6 Moving from the macro- to the micro- and nano-scale provides several advantages: less sample required, greater sensitivity because of increased contact area between analyte solution and device surface, and the possibility of developing portable, or even implantable devices through miniaturization. Silica – amorphous silicon dioxide – is one of the common materials used to fabricate micro- and nano-scale biomedical devices. Some of the benefits of silica include low toxicity of the amorphous oxide,7, 8 cost, and an accumulated expertise in fabrication methods. The latter includes expertise in silicon device fabrication because fluids come in contact with a surface layer of amorphous oxide in silicon-based devices.9–11 Silicon dioxide is the principle component of various materials referred to as “glass,” such as soda lime or borosilicate glass. Hence the behavior of bio-molecules near the amorphous silica surface is of both scientific and technological importance. In this work, we examine the behavior of a DNA oligomer near an amorphous silica surface with a charge density typical of neutral pH. It is known that DNA binds to the silica surface, even though both the silica surface and DNA are negatively charged. DNA binding to silica in the presence of chaotropic salts, such as guanidinium thiocyanate and sodium perchlorate, is a standard purification method.12–18 There are also many reports of DNA binding to silica in the absence of chaotropic slates, and at pH values where the silica surface is negatively charged. In the absence of chaotropic salts, Fujiwara et al. found weak binding of double-stranded DNA to silica under neutral or basic pH values. However the binding increased as the pH was lowered toward the isoelectric point of silica.19 Also, the binding of DNA to silica nanoparticles in a pH 7, salt-free solution was detected by Marcie and Savoie20 via a small change in the intensity and frequency of a Raman spectral band at 1092 cm−1 , which was assigned to the DNA phosphate group. To access the type of interactions likely to occur when DNA interacts with silica, Kwon et al.21 performed atomic force microscope measurements in which the AFM probe was a silica particle and the surface was coated with phosphorylated dextran. At neutral pH, the pull-off force for phosphorylated dextran was 80% that of normal dextran, but still exhibited a measurable attraction. Vandeventer et al.22 monitored bulk found adsorption of DNA to silica particles in the presence and absence of chaotropic salts. Additional experimental evidence for DNA binding to silica

2

ACS Paragon Plus Environment

Page 2 of 39

Page 3 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

at neutral pH was recently provided by Vandeventer et al.22 and Smerkova et al.23 Single-stranded DNA is found to bind more strongly than double-stranded DNA to silica,24–26 with appreciable binding at pH values above the isoelectric point of silica up to neutral pH.24, 25 This has been explained by increased opportunity for hydrophobic interactions between the unpaired bases of single-stranded DNA and the silica surface.24–26 In this paper we confirm this mechanism, but suggest additional mechanisms that lead to stronger binding of single-stranded DNA. There are many experiments in which DNA is introduced into silica nanopores, often with the goals of chromatographic separation of different oligomers, furnishing a medium for delivery of DNA, or DNA sequencing. The effect of DNA on the ion current is often the experimental method to track the presence, and ultimately the sequence, of DNA. Viovy,27 Baldessari and Santiago28 Branton et al.29 and Zwolak and Di Ventra30 provide general reviews of progress in electrokinetic separation of biomolecules in nano-channels. Since biomolecules are never far from the surface in nano-pores, a key property is the adsorption/desorption behavior of biomolecules at channel walls, which we explore in this work for DNA/silica systems. The numerous reported observations of DNA in micro- and nano-pores include electrokinetic separation of DNA of different lengths in glass nano-channels,31–33 capillary electrophoresis in silica micro-channels to separate DNA fragments of varying size,34 tracking the distribution of DNA in silica micro-channels under various conditions, including the presence of Poiseuille flow.35 Both increase and decrease of ion current as DNA passes through silica nano-pores has been observed, depending on salt concentration.36 Other experiments probe the diffusivity of DNA in silica slit nano-pores with gap heights between 75 and 545nm,37 and the electrophoretic force on a single DNA molecule passing through a SiO2 /SiN/SiO2 pore.38 Zhang and coworkers39 observed translocation of DNA through cylindrical pores in silica. They found translocation in the direction of EOF for the smallest 9nm diameter pores, while translocation was in the direction of DNA electrophoresis for 75nm diameter pores.40 Pure silica nano-pores have not met with the same degree of success as other inorganic or biological pore structure in constructing potential DNA sequencing devices, but silica has been employed in combination with other materials.40–42 Given the extent of experimental interest in systems involving intimate contact between DNA and silica, theoretical investigations are warranted to provide molecular insight into various outstanding issues, such

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

as the mechanism by which negatively charged DNA binds to negatively charged silica, and why singlestranded DNA binds preferentially over double-stranded DNA. To date, there have been some theoretical efforts directed toward DNA in nanochannels, most often using continuum theory43–48 There have been several molecularly detailed simulations of DNA near surfaces, but none with a realistic description of DNA in contact with amorphous silica, as developed in this work. Aksimentiev and co-workers when they studied translocation of DNA in Si3 N4 nanopores.49–52 Luan et al. simulated motion of DNA in a channel with regions of metallic and amorphous SiO2 walls.53–55 However, the SiO2 in contact with water was not hydroxylated and charged, and the DNA restrained from contact with the walls. DNA confined near graphene56 and carbon nanotubes57 have also been explored in molecular simulations. Herein, we describe our efforts to elucidate the DNA-amorphous silica binding mechanism for both double-stranded DNA (dsDNA) and single-stranded DNA (ssDNA). Since there are many important applications where biomolecules interact with amorphous silica, our goal in section 2 was to devise a general procedure to estimate interactions between those biomolecules and a silica surface, which we intend to apply . Then in section 3 we specify the methods used in our molecular dynamics and free energy calculations. Results are presented in section 4, followed by section 5 where discussion and conclusions are given.

2 Interaction model development Our strategy to develop a model for biomolecules at the silica/water interface is to interpolate ab initio quantum chemical results between various probe molecules and silica fragments. We interpolate in two senses. First, like any force field, we assume that parameters for atoms in similar bonding situations are transferable. Secondly, we use standard combining rules58 to infer interactions where we have not generated quantum chemical data. These approximations are in standard use for current force fields,59–61 and can capture qualitative trends for systems that are currently beyond the reach of more rigorous methods. Interactions of probe molecules (Fig. 1) with four fragments from the silica surface (Fig. 2) were investigated. The representative probe molecules contain functional groups often found in biomolecules – methane (CH4 ), methanol (CH3 OH), ammonium (NH+4 ), acetate (CH3 COO− ) and benzene (C6 H6 ) This set of molecules includes non-polar, polar, charged, and aromatic species. Some silica fragments contained no 4

ACS Paragon Plus Environment

Page 4 of 39

Page 5 of 39

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 39

a dissociated silanol group is exposed to a probe molecule (Fig. 2c), and the largest fragment containing 5 silicon atoms used for benzene as a probe in which two silanol groups interact with the probe (Fig. 2d). All cases contained siloxane bonds in addition to possible silanols or dissociated silanols. It should be noted that in the last fragment, the benzene probe had significant exposure to siloxane groups as well as two silanol groups. All results of ab initio calculations are provided in the Supplementary Information. Table 1: Lennard-Jones potential parameters of siloxane (OX ), silanol (OH ), and dissociated silanol (OM ) oxygens that were optimized to best reproduce ab initio data from fragment calculations. The charge parameters are derived from the BKS potential.69 silica atom types charge (e) σ(Å) ǫ(eV) OX (siloxane) −1.2 3.02 4.691 × 10−2 OH (silanol) −1.2 3.06 9.086 × 10−3 OM (dissociated silanol) −1.6 2.95 2.352 × 10−2 Si 2.4 – 0 H (silanol) 0.6 0 0 Using a set of Lennard-Jones parameters for the interaction potential between like atoms, we employed the Lorentz-Berthelot combination rules58 to generate Lennard-Jones parameters for unlike atom pairs.

σi j =

 1 √ σii + σ j j , ǫi j = ǫii ǫ j j 2

(1)

In Eq. (1), σii , σ j j , ǫii , ǫ j j are like-atom Lennard-Jones parameters which are used to estimate σi j , ǫi j , the unlike-atom parameters. Our strategy is to parametrize the interaction model with like-atom Lennard-Jones parameters for siloxane (OX ), silanol (OH ), and dissociated silanol (OM ) atom types, and take analogous likeatom Lennard-Jones parameters for biomolecule atom types from available force fields. The full interaction potential is then generated with the combining rules. We optimized Lennard-Jones parameters for the three atom types of the silica surface (OX ,OH ,OM ) other than silicon to achieve a best fit to the ab initio data when using Lorentz-Berthelot combining rules.58 Details are given in the Supplementary Information, and the resulting parameters are summarized in Table 1.

6

ACS Paragon Plus Environment

Page 7 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3 Simulation procedures The DL POLY simulation package70–72 was used to equilibrate the silica slab. All other simulations were performed with the GROMACS program version 4.5.5.73 Interactions within DNA and with water were described with the AMBER99sb force field,74, 75 originally based on the ffamber port to GROMACS76, 77 and now a standard part of GROMACS. A double-stranded DNA with the sequence d(CGCGAATTCGCG)·d(CGCGAATTCGCG), the Drew-Dickerson dodecamer,78 and a single-stranded DNA with the sequence (CGCGAATTCGCG) were used in the simulations. The amorphous silica was made from cycle I-IV of Huff and co-workers.79 To generate the silica slab, the periodic simulation cell was lengthened in one direction, cleaving the silica and introducing an empty region. The exposed silica surface was subsequently annealed, hydroxylated by replacing 2M rings to vicinal silanols and nonbrigding oxygen atoms (NBOs) to geminal silanol pairs, and re-annealed as described in ref. 67. Simulations were conducted using SPC/E water model, and equilibrated after introduction of water. Each simulation is performed in a 7.96nm × 8.61nm × 9.74nm box, the long dimension perpendicular to the silica surface. There are 1611416333 water molecules in the box and 130 (119) sodium ions which neutralize the ds-DNA (ss-DNA) and silica. The number of water molecules was adjusted to achieve a mass density of 1.0g cm−3 in the channel region. The thickness of the silica wall is about 2.7nm and the height of the channel is about 7.0nm. The Lennard-Jones potential parameters for sodium ions are from AMBER99sb force field.74, 75

ssDNA dsDNA

binding events 7 8

phosphate 19 11

binding sites base (hydrophobic) deoxyribose sugar 7 4 1 1

Table 2: Binding statistics for 8 independent simulations for a ssDNA oligomer, 8 independent simulations for a dsDNA oligomer. All of the dsDNA simulations resulted in binding to an undissociated silica surface, while 7 of the 8 ssDNA runs led to binding.

Neutral silica surfaces with no dissociated silanols were employed to construct initial bound configurations. For both dsDNA and ssDNA, 8 simulations, each of length 30ns, were performed in which DNA was introduced 7Å from one of the neutral silica surfaces with axis parallel to the surface. After evolving in time for 30ns, all of the 8 dsDNA simulations, and 7 out of 8 ssDNA simulations were closely associated with the surface. The prevalence of binding modes, which are discussed in detail in section 4, are summarized in 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 2. The number of binding sites in each event ranged from one to seven. In experimental situations, the silica surface charge density depends on solution pH and electrolyte concentration.80, 81 Since most experiments are conducted near neutral pH, we studied a surface with a charge density of −0.126C/m2 , typical of the neutral range.82 Our silica slab was dissociated by deleting the proton of randomly chosen silanol groups. Then the same number of sodium ions were added to make the system neutral. Several binding configurations were chosen for detailed analysis. They were further equilibrated on the dissociated silica surface for 10ns. In all cases, the DNA oligomer remained bound after charging the surface, although in some cases the bound state was metastable. In the discussion section, we argue this likely underestimates the degree of DNA binding to negatively charged silica. The binding free energy was calculated by umbrella sampling83 combined with the weighted histogram analysis method (WHAM).84–87 A series of simulation windows centered about steadily increasing distances between the binding group, such as phosphate groups or bases, and surface were performed. The umbrella potentials were based on a reaction coordinate which was the distance from the center of mass of the detaching group, either a phosphate or a base, and a silicon atom directly under the group in the bound configuration. In some cases, histograms involving two reaction coordinates were accumulated. When the oligomer is close to its relaxed binding position, the interaction with the silica surface is strong. Hence, a high spring force constant was used to influence the DNA in the presence of strong interactions near the surface. A sharp position distribution is generated by a high spring constant, so the distance between window centers was small near the silica surface. Alternatively, when the oligomer is far from its binding position, the interaction between DNA and silica surface is weaker. A weaker spring force constant is enough to influence the DNA, which gives wider position distribution. The distance between window centers was allowed to grow with increasing distance from the surface. Distance between windows and force constants are included in the supporting information. The configuration obtained at the end of window i was used as the starting point for the next window i + 1. Each window was equilibrated for 1ns and then data was collected for another 1ns. During the sampling of each window, the values of reaction coordinate r, generally the distance between a point on the DNA and another point on the surface, were obtained in order to generate (b) (b) a biased probability histogram, P(b) i (r) in window i. We verified that for all windows i, Pi (r) and Pi+1 (r),

histograms for adjacent windows, had significant overlap. 8

ACS Paragon Plus Environment

Page 8 of 39

Page 9 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Following the WHAM procedure,84–87 the biased probability histograms were combined to obtain the unbiased probability.

P(r) = C

X i

P

j nj

ni P(b) i (r) −β(w (r)− f ) j j e

(2)

In the above equation, wi is the biasing potential, ni is the simulation length in window i, and fi is the free energy difference between ith biased system and unbiased system. The fi are obtained by iteratively solving the following self-consistency condition.

e−β fk = C

Z

dr e−βwk (r)

X i

ni P(b) i (r) P

j n je

−β[w j (r)− f j ]

.

(3)

In this equation the integration over r is replaced by a sum over histogram elements. The reaction cooridnate r usually stands for a single variable, but in two cases we studied hydrophobic base-silica interaction using 2D umbrella sampling.

4 Results

Figure 3: A phosphate of a dsDNA is bound to two silanol groups of the silica surface. This binding configuration corresponds to the curve labeled dsDNA1 in Fig. 4. The inset in the lower left is an expanded view of the phosphatesilanol interaction site, where the two silanols interacting with a phosphate oxygen are indicated with a dashed line.

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4.1

Phosphate-silanol interaction

For both dsDNA and ssDNA, we frequently observed phosphate groups of the DNA binding to one or two silanols on the silica surface, as shown in Fig. 3. Umbrella sampling and WHAM provided the the free energy profiles in Fig. 4 for five cases of phosphate-silanol binding. The reaction coordinate was the distance from the phosphate center of mass to a silicon atom directly beneath it in the bound configuration. In all but the run labeled ssDNA2, the phosphate-silanol bond was located toward one end of the oligomer, while the other end of the DNA was sufficiently far from the surface so that it was not interacting. As the biasing potential was changed, both ends of the oligomer moved out of interaction range with the surface. Therefore, pulling apart the phosphate-silanol bond provided an estimate of the overall free energy surface for binding via a single phosphate-silanol bond. In two cases, labeled dsDNA1 and ssDNA1 in Fig. 4, a single phosphate oxygen was in close proximity to two silanol groups. To confirm that the entire range of interaction was included in Fig. 4, we extended the umbrella sampling range, this time adding a biasing potential to two sites to keep the oligomer parallel to the surface. Over an extended range which was 0.5nm more than the range shown in Fig. 4, we observed no change in the free energy. The two dsDNA cases shown in Fig. 4, dsDNA1 and dsDNA2, where binding occurred at just a single phosphate-silanol site, were converted into ssDNA1 and ssDNA3, respectively, by removal of the chain that did not participate in binding with the surface. The free energy profile for binding was re-determined for the remaining chain using umbrella sampling and WHAM. Thus the pairs, dsDNA1-ssDNA1 and dsDNA2ssDNA3, furnish a direct comparison of binding of double- and single-stranded DNA with other factors held nearly the same. The run labeled ssDNA2 in Fig. 4 was somewhat different than the rest. In this run, the ssDNA was attached to the surface by two phosphate silanol interactions, one at the 4th nucleotide from the 3′ end, and another phosphate at the 5′ end. The ssDNA2 curve in Fig. 4 is the free energy to break the phosphate silanol interaction at the 4th nucleotide while the chain remained bound to the surface at the 5′ end. In the bound state, the phosphate-surface distance, as described in the preceding paragraph, was 0.42 to 0.46nm from the surface. Out until 0.54nm from the surface, specific attractive interactions between the phosphate groups and surface contribute strongly to the free energy surface. In the range between 0.55nm

10

ACS Paragon Plus Environment

Page 10 of 39

Page 11 of 39

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

from the center of mass of the phosphate group to a silicon lying directly below it, not the closest approach between atoms of the phosphate and the silica surface.) While the repulsive portion of the free energy surfaces in Fig. 4 are quite similar among the ssDNA and dsDNA groups, the overall binding free energy varies for different binding locations. Each phosphate has a binding free energy in the range of +20kJ mol−1 (metastable binding) to −41kJ mol−1 . Here the terms stable and metastable are used to distinguish when the bound configuration lies lower or higher in free energy, as compared with the free DNA oligomer. Of course, a dynamic equilibrium between bound and unbound DNA will exist, as directly observed, for example, by Kang et al.24 More precisely, stability or metastability predicts the relative amount of bound and unbound DNA. Three factors cause the large range of overall binding free energy: First, because phosphate location is constrained by its position in the DNA polymer, an individual phosphate oxygen may not be able to achieve its most favorable configuration relative to its silanol bonding partner. Second, the local charge distribution of the silica surface near the binding sites may affect the binding free energy. If there are more negatively charged sites around the binding site, the binding is weaker. For example, the phosphate-silanol binding site of ssDNA3 in Fig. 4 lies close to four dissociated silanol groups that provide negative charge around the binding site. As a consequence, this is the one instance of ssDNA for which binding is metastable. Finally, the inherent differences between ssDNA and dsDNA produce a range of binding energies. ssDNA tends to give stronger binding than dsDNA because ssDNA is more flexible and able to optimize its interactions with the surface, and because the more negatively charged dsDNA produces greater electrostatic repulsion to silica surface than ssDNA. Using B3LYP quantum chemical methods, two groups have estimated the binding energy of a phosphate group to an uncharged silica fragment. Murashov and Leszczynski calculated the gas phase binding energies of an orthosilicic acid unit and either a phosphate or dimethyl phosphate ion.89 They found an average binding energy of 59kJ mol−1 per phosphate-silanol bond. Kwon et al. calculated the binding enthalpy of methyl glucosidic phosphate to an octa-hydroxy silsesquioxane fragment. They obtained a binding enthapy of −129kJ mol−1 for a fragment with two phosphate-silanol bonds (along with other interactions which are much weaker).21 Again, the gas phase binding energy came out to roughly −60kJ mol−1 . These sets of

12

ACS Paragon Plus Environment

Page 12 of 39

Page 13 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

results compare nicely with the largest drop in free energy from the maximum near 0.55nm to the bottom of the well observed in Fig. 4. The quantum chemical calculations21, 89 discussed above employed a neutral silica fragment. Therefore, the presence of nearby negative charges on the silica and on the DNA can lessen the magnitude of the binding energy.

4.2

Hydrophobic bonding

Bases of both single stranded DNA (ssDNA) and double stranded DNA (dsDNA) exhibit a tendency to lay adjacent to regions of the silica surface which, because of natural statistical variations, are depleted of silanol groups. We have previously shown that such silanol-free regions of the silica surface are hydrophobic, i.e. that the density of water near these regions is less than the bulk density, and considerably less than the water density near silanol groups.67 The hydrophobic behavior of regions lacking silanol groups was verified by ab initio molecular dynamics simulations.68

(a)

(b)

Figure 5: Guanine and a cytosine of ssDNA oligomer bound to hydrophobic regions of silica surface. (a) and (b) are the same picture but DNA is invisible in (b). Larger red surface atoms are dissociated silanols and white dots are hydrogen atoms. In both snapshots, water and counter-ions are not shown.

Fig. 5 shows an example of ssDNA binding to silica via hydrophobic interactions. Two adjacent bases, guanine and cytosine, were in contact with the silica surface, facing regions depleted of silanol groups. In the figures, a dashed line indicates the binding region. Figs. 5(a) and 5(b) show the same physical configuration. 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Water and ions are left out of the picture in Fig. 5(a), and additionally the DNA is not depicted in Fig. 5(b), allowing one to see that the binding regions are depleted of silanol groups. Fig. 6 shows an example of hydrophobic binding of dsDNA. In the double-stranded case, the hydrophobic interaction requires a base to break its hydrogen bonds with its complement, incurring a free energy cost that is not present in the single-stranded case. In our studies, we only observed the end base pair of dsDNA breaking apart to enable hydrophobic interaction with the surface.

(a)

(b)

Figure 6: A cytosine of dsDNA is bound to hydrophobic region of silica surface. (a) and (b) are the same picture but DNA is invisible in (b). Large red spheres are dissociated silanols and white dots are hydrogen atoms.

(a)

(b)

Figure 7: Hydrogen bonding between bases and silica surface concurrent with hydrophobic interactions. (a) Hydrogen bonding between oxygen atom on cytosine and hydrogen atom on silanol group. (b) Hydrogen bonding between NH2 on cytosine and oxygen on silanol group.

14

ACS Paragon Plus Environment

Page 14 of 39

Page 15 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In most cases, concurrent with hydrophobic binding, some hydrogen bonds are formed between bases and the silica surface at the edge of the hydrophobic region. Fig. 7(a) shows a hydrogen bond between the oxygen atom of a hydrophobically bound cytosine from ssDNA and the hydrogen of a silanol group. Simultaneous hydrophobic attraction and hydrogen bonding at the periphery is possible via a ketone oxygen for G, C and T, and via an -NH or NH2 group for all bases, provided the hydrophobic site contains the required hydrogen bonding sites on its periphery. Fig. 7(b) shows a hydrogen bonding between the NH2 group on a cytosine and oxygen on a silanol group. Adenine and guanine also contain a NH2 group, and may engage in simultaneous hydrophobic/hydrophilic binding. We also found one case, discussed below, where the hydroxyl group of a deoxyribose forms a hydrogen bond close to the site of hydrophobic binding. We now examine the free energy surfaces, as obtained from umbrella sampling and WHAM, which govern the approach of the DNA oligomer to the silica surface. For the ssDNA case depicted in Fig. 5, two bases, guanine and cytosine at the 5′ end, were constrained by a harmonic potential in each sampling window. Because ssDNA is so flexible, the z-component of the centers of mass of bases 3 and 7 counting from the 5′ end were constrained to prevent new attachment forming during this process. The x- and y-components of bases 3 and 7 were free to adjust during the detachment of the guanine and cytosine at the end of the chain to prevent strain developing during the detachment. The change in reaction coordinates and free energy with successive umbrella sampling windows is shown in Fig. 8. Detachment from the silica surface was accomplished in three steps. First, the distance between the surface and guanine center of mass was steadily increased by 0.913nm in sampling windows 1-39 while distance between surface and cytosine was fixed, as traced by the green and red curves in Fig. 8. This breaks most of the guanine-surface interaction without significant distortion of the DNA. There are no associated hydrogen bonds broken during this initial stage of the detachment. At the end of this series of windows, the free energy levels off, as can be seen in the plateau of the purple curve in windows 33-39. After guanine was detached from the surface, the cytosine center of mass was steadily moved by 1.01nm from the surface in the series of sampling windows 40-83 while guanine was fixed. Again, this breaks most of the cytosine-surface interaction without distortion of the DNA. During the latter stages of cytosine detachment from the surface a hydrogen bond involving a deoxyribose hydroxyl group is broken (blue curve in Fig. 8). The rise of the sugar OH-surface distance in Fig. 8 be15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 16 of 39

Page 17 of 39

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 18 of 39

Page 19 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

group of the silica surface. A binding free energy of 8.25kJ mol−1 for a sugar OH group was established in section 4.2 (see Fig. 10).

Figure 11: Hydrogen bonding for deoxyribose sugar of ssDNA.

5 Discussion and Conclusion One of the main conclusions of our work is that phosphate-silanol and hydrophobic interactions are the principal mechanisms for binding of DNA to silica. These short-range attractions can be sufficiently strong to overcome the electrostatic repulsion between negatively charged DNA and a negatively charged silica surface. While DNA and silica above its isoelectric point (∼ pH 3) are both negatively charged, the repulsion between silica and DNA is weak because it is effectively screened by the counter-ions near both surfaces. Both silica and DNA are composed of functional groups that may be either positive or negative, hydrophobic or hydrophilic. Binding is a matter of overcoming charge repulsion by aligning functional groups with complementary properties. Hydrophobic attachments, where a nucleobase lies against a region of the silica surface free from hydrophilic silanols, each contributed from −73 to −190kJ mol−1 . Including cases that were not discussed in section 4.1, each silanol-phosphate interaction, where the oxygen of a phosphate group closely approaches the hydrogen of a surface silanol, can contribute over −65kJ mol−1 to the overall binding energy. Cases where the DNA oligomer formed only one phosphate-silanol interaction were presented in section 4.1. We also encountered

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

cases where the oligomer formed more than one phosphate-silanol attachment. Per attachment, we generally found hydrophobic interactions to be the strongest interactions. However, phosphate-silanol binding interactions strongly outnumbered hydrophobic interactions (Table 2), especially for dsDNA. Hydrophobic interactions were so infrequent in our dsDNA simulations because only a terminal base pair unfolded to allow one member of the pair to lie against the surface. Our observation that the hydrogen bonds between terminal base pairs of dsDNA were most labile is in accord with previous simulations. In prior work on the same Drew-Dickerson dsDNA dodecamer, the terminus of the dsDNA was observed to occasionally fray or open,92, 93 or for the gap between base pairs to widen near the terminals.94, 95 Hence phosphate-silanol and hydrophobic interactions were of roughly equal importance for ssDNA, while phosphate silanol interactions are more important for dsDNA. In comparison, auxiliary polar interactions at the periphery of a hydrophobic site, and surface-sugar hydroxyl interactions played a minor role. In other simulation work, hydrophobic interactions between DNA bases and ordered surfaces have been reported. Aksimetiev, Schulten and co-workers reported hydrophobic attachment of DNA bases to crystalline Si3 N4 .49–51 Graphene56 and carbon nanotubes57 were also found to participate in hydrophobic interactions with DNA bases. Another main result is the identification of several factors which explain why ssDNA is attracted to silica more than dsDNA: 1) ssDNA is more flexible, and therefore able to maximize the number of complementary DNA-silica interactions. In a series of runs in which DNA was placed near a neutral DNA surface, ssDNA formed 4 binding sites per run, while dsDNA formed less than 2. The number of each type of interaction – phosphates-silanol, hydrophobic, and sugar hydryoxyl – was less for dsNDA compared to ssDNA. 2) ssDNA has an increased ability for hydrophobic attachment to silica because the bases of ssDNA do not have to unpair with their hydrogen bond partners, as do the bases of dsDNA. 3) The linear charge density of dsDNA is twice that of ssDNA. Our calculated binding free energies are likely underestimates. To explore DNA binding to silica, lengthy molecular dynamics simulations are required, including extensive sampling for free energy calculations. Out of necessity, this work is based on empirical force fields. In addition, after DNA was allowed to bind to silica, we randomly dissociated silanol groups by removing the hydrogens of randomly selected silanol

20

ACS Paragon Plus Environment

Page 20 of 39

Page 21 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

groups and changing the oxygen species from hydroxyl oxygen to O− . This procedure sometimes left surface O− underneath the DNA chain. We probably underestimated the binding free energy because in actuality the dissociation of silanols would not be random. Instead, silanol dissociation would react to the DNA placement so as to minimize the overall free energy. Experimental evidence for DNA binding to silica under conditions where the silica surface carries a negative charge has come in various forms, as summarized in the introduction. Typically, some of the factors listed above have been proposed to explain binding in any one work. Yeung and co-workers have proposed that hydrophobic interactions and hydrogen bonding of unpaired bases are the main driving force of adsorption.26 Mao and co-workers studied DNA binding to silica by Fourier-transform infrared spectroscopy.96 The results were explained as the formation of phosphate-silanol hydrogen bonds. Later, Mercier and Savoie studied the interaction of DNA with various types of silica particles by infrared and Raman spectroscopy.20 They disputed the spectroscopic assignments of Mao et al., although they also concluded that the spectrum indicated phosphate-silanol interactions. Kwon and co-workers21 also found that phosphate groups formed strong H-bonds with neutral silanols ( 129kJmol−1 ). Finally, Kwon et al. pulled DNA from a phosphorylated dextran surface in atomic force microscopy experiments. They interpreted their data in terms of phosphatesilanol attachments. We conclude with perhaps a provocative statement. Much has been written about the uniqueness of the water-biomolecule interface.97–100 Silica is a material that is not considered to be an analog to biological molecules. Yet, the silica surface is a disordered patchwork of hydrophilic and hydrophobic regions, much like the surface of a protein. Indeed, the interaction of HIV nucleocapsid protein and DNA has been discussed as a combination of hydrogen bonding and hydrophobic interactions, in terms very similar to our treatment of the DNA-silca interaction.101 Therefore, comparison of water near biological molecules with water near silica may be a useful way to ascertain which properties of water near biological molecules are truly unique.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Acknowledgements This material is based upon work supported by the National Science Foundation under Grant No. EEC0914790. The calculations reported here were made possible by a grant of resources from the Ohio Supercomputer Center. This information is available free of charge via the Internet at http://pubs.acs.org

22

ACS Paragon Plus Environment

Page 22 of 39

Page 23 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

References [1] Desai, T. A.; Hansford, D. J.; Kulinsky, L.; Nashat, A. H.; Rasi, G.; Tu, J.; Wang, Y.; Zhang, M.; Ferrari, M. Nanopore Technology for Biomedical Applications. Biomedical Microdevices, 1999. 2, 11–40 [2] Malsch, N. H. Biomedical Nanotechnology. CRC Press, Boca Raton, 2005 [3] Desai, T. A.; et al. Nanoporous Implants for Controlled Drug Delivery. Desai, T. A.; Bhatia, S., editors, Therapeutic Micro/Nanotechnology, volume III of BioMEMS and Biomedical Nanotechnology. Springer, 2006. 263–286 [4] Meng, E. Biomedical Microsystems. CRC Press, Boca Raton, 2010 [5] Bagchi, D. Bio-Nanotechnology: A Revolution in Food, Biomedical and Health Sciences. WileyBlackwell, Chichester, West Sussex, 2013 [6] Berezin, B. D. Nanotechnology for Biomedical Imaging and Diagnostics: From Nanoparticle Design to Clinical Applications. Wiley, Hoboken, New Jersey, 2014 [7] Fubini, B. Health effects of silica. Legrand, A. P., editor, The Surface Properties of Silicas. Wiley, New York, 1998. 415–464 [8] Saffiotti, U. Silicosis and Lung Cancer: a Fifty-Year Perspective. Acta Bio-Med. Ateneo Parmense, 2005. 76, 30–37 [9] Sigmon, T. W.; Chu, W. K.; Lugujjo, E.; Mayer, J. W. Stoichiometry of Thin Silicon Oxide Layers on Silicon. Appl. Phys. Lett., 1974. 24, 105–107 [10] Dabrowski, J.; M¨ussig, H.-J. Silicon Surfaces and Formation of Interfaces: Basic Science in the Industrial World. World Scientific, River Edge, NJ, 2000 [11] Yao, S.; Myers, A. M.; Posner, J. D.; Rose, K. A.; Santiago, J. G. Electroosmotic Pumps Fabricated From Porous Silicon Membranes. J. Microelectromech. Syst., 2006. 15, 717–728

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[12] Melzak, K. A.; Sherwood, C. S.; Turner, R. F. B.; Haynes, C. A. Driving Forces for DNA Adsorption to Silica in Perchlorate Solutions. J. Colloid Interface Sci., 1996. 181, 635–644 [13] Boom, R.; Sol, C. J.; Salimans, M. M.; Jansen, C. L.; Wertheim-van Dillen, P. M.; van der Noordaa, J. Rapid and Simple Method for Purification of Nucleic Acids. J. Clin. Microbiol., 1990. 28, 495–503 [14] Towner, P. Purification of DNA. Brown, T. A., editor, Essential Molecular Biology, volume 1, chapter 3. Oxford, London, 2nd edition, 2000. 55–67 [15] Cady, N. C.; Stelick, S.; Batt, C. A. Nucleic Acid Purification Using Microfabricated Silicon Structures. Biosens. Bioelectron., 2003. 19, 59–66 [16] Nanassy, O. Z.; Haydock, P. V.; Reed, M. W. Capture of Genomic DNA on Glass Microscope Slides. Anal. Biochem., 2007. 365, 240–245 [17] Li, X.; Zhang, J.; Gu, H. Adsorption and Desorption Behaviors of DNA with Magnetic Mesoporous Silica Nanoparticles. Langmuir, 2011. 27, 6099–6106 [18] Li, X.; Zhang, J.; Gu, H. Study on the Adsorption Mechanism of DNA with Mesoporous Silica Nanoparticles in Aqueous Solution. Langmuir, 2012. 28, 2827–2834 [19] Fujiwara, M.; Yamamoto, F.; Okamoto, K.; Shiokawa, K.; Nomura, R. Adsorption of Duplex DNA on Mesoporous Silicas: Possibility of Inclusion of DNA into Their Mesopores. Anal. Chem., 2005. 77, 8138–8145 [20] Mercier, P.; Savoie, R. Interaction of DNA with Silica Particles: A Vibrational Spectroscopic Study. Biospectroscopy, 1997. 3, 299–306 [21] Kwon, K. D.; Vadillo-Rodriguez, V.; Logan, B. E.; Kubicki, J. D. Interactions of Biopolymers with Silica Surfaces: Force Measurements and Electronic Structure Calculation Studies. Geochimica et Cosmochimica Acta, 2006. 70, 3803–3819 [22] Vandeventer, P. E.; Lin, J. S.; Zwang, T. J.; Nadim, A.; Johal, M. S.; Niemz, A. Multiphasic DNA Adsorption to Silica Surfaces Under Varying Buffer, pH, and Ionic Strength Conditions. J. Phys. Chem. B, 2012. 116, 5661–5670 24

ACS Paragon Plus Environment

Page 24 of 39

Page 25 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[23] Smerkova, K.; Dostalova, S.; Vaculovicova, M.; Kynicky, J.; Trnkova, L.; Kralik, M.; Adam, V.; Hubalek, J.; Provaznik, I.; Kizek, R. Investigation of Interaction Between Magnetic Silica Particles and Lambda Phage DNA Fragment. J. Pharm. Biomed. Anal., 2013. 86, 65–72 [24] Kang, S. H.; Shortreed, M. R.; Yeung, E. S. Real-Time Dynamics of Single-DNA Molecules Undergoing Adsorption and Desorption At Liquid-Solid Interfaces. Anal. Chem., 2001. 73, 1091–1099 [25] Li, H.-W.; Park, H.-Y.; Porter, M. D.; Yeung, E. S. Single DNA Molecules as Probes of Chromatographic Surfaces. Anal. Chem., 2005. 77, 3256–3260 [26] Isailovic, S.; Li, H.-W.; Yeung, E. S. Adsorption of Single DNA Molecules At the Water/fused-Silica Interface. J. Chromatogr., 2007. A1150, 259–266 [27] Viovy, J.-L. Electrophoresis of DNA and Other Polyelectrolytes: Physical Mechanisms. Rev. Mod. Phys., 2000. 72, 813–872 [28] Baldessari, F.; Santiago, J. Electrophoresis in Nanochannels: Brief Review and Speculation. J. Nanobiotechnol., 2006. 4, 12 [29] Branton, D.; et al. The Potential and Challenges of Nanopore Sequencing. Nat. Biotechnol., 2008. 26, 1146–1153 [30] Zwolak, M.; Di Ventra, M. Colloquium: Physical Approaches to DNA Sequencing and Detection. Rev. Mod. Phys., 2008. 80, 141–165 [31] Culbertson, C. T.; Ramsey, R. S.; Ramsey, J. M. Electroosmotically Induced Hydraulic Pumping on Microchips: Differential Ion Transport. Analytical Chemistry, 2000. 72, 2285–2291 [32] Ramsey, J. M.; Alarie, J. P.; Jacobson, S. C.; Peterson, N. J. Molecular Transport Through Nanometer Confined Channels. Baba, Y.; Berg, A.; Shoji, S., editors, Micro Total Analysis Systems, volume 1. Springer Netherlands, Dordrecht, 314–316. (Proceedings of a symposium in Nara, Japan, 3-7 November 2002) [33] Petersen, N. J.; Alarie, J. P.; Jacobson, S. C.; Ramsey, J. M. Polyelectrolyte Transport in Nanoconfined Channels. Northrup, M. A.; Jensen, K. F.; Harrison, D. J., editors, Proceedings of Micro Total 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analysis Systems, volume 2. Transducers Research Foundation, Squaw Valley, California, USA, 701– 703 [34] Zabzdyr, J. L.; Lillard, S. J. Novel Elution Strategy for Monitoring DNA Counter-Migration in the Presence of Electroosmotic Flow. J. Chromatogr. A, 2004. 1040, 283–290 [35] Zheng, J.; Yeung, E. S. Anomalous Radial Migration of Single DNA Molecules in Capillary Electrophoresis. Anal. Chem., 2002. 74, 4536–4547 [36] Smeets, R. M. M.; Keyser, U. F.; Krapf, D.; Wu, M.; Dekker, N. H.; Dekker, C. Salt Dependence of Ion Transport and DNA Translocation Through Solid-State Nanopores. Nano Lett., 2006. 6, 89–95 [37] Balducci, A.; Mao, P.; Han, J.; Doyle, P. S. Double-Stranded DNA Diffusion in Slitlike Nanochannels. Macromolecules, 2006. 39, 6273–6281 [38] Storm, A. J.; Chen, J. H.; Ling, X. S.; Zandbergen, H. W.; Dekker, C. Fabrication of Solid-State Nanopores with Single-Nanometre Precision. Nat. Mater., 2003. 2, 537–540 [39] Zhang, B.; Wood, M.; Lee, H. A Silica Nanochannel and Its Applications in Sensing and Molecular Transport. Anal. Chem., 2009. 81, 5541–5548 [40] Liang, X.; Morton, K. J.; Austin, R. H.; Chou, S. Y. Single Sub-20 Nm Wide, Centimeter-Long Nanofluidic Channel Fabricated by Novel Nanoimprint Mold Fabrication and Direct Imprinting. Nano Lett., 2007. 7, 3774–3780 [41] Iqbal, S. M.; Akin, D.; Bashir, R. Solid-State Nanopore Channels with DNA Selectivity. Nat. Nanotech., 2007. 2, 243–248 [42] Liang, X.; Chou, S. Y. Nanogap Detector Inside Nanofluidic Channel for Fast Real-Time Label-Free DNA Analysis. Nano Lett., 2008. 8, 1472–1476 [43] Ghosal, S. Electrokinetic-flow-induced viscous drag on a tethered DNA inside a nanopore. Phys. Rev. E, 2007. 76, 061916 [44] Wong, C. T. A.; Muthukumar, M. Polymer Capture by Electro-Osmotic Flow of Oppositely Charged Nanopores. J. Chem. Phys., 2007. 126, 164903 26

ACS Paragon Plus Environment

Page 26 of 39

Page 27 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[45] Das, S.; Chakraborty, S. Transverse Electrodes for Improved DNA Hybridization in Microchannels. AIChE J., 2007. 53, 1086–1099 [46] Duong-Hong, D.; Han, J.; Wang, J.; Hadjiconstantinou, N.; Chen, Y.; Liu, G. Realistic Simulations of Combined DNA Electrophoretic Flow and EOF in Nano-Fluidic Devices. Electrophoresis, 2008. 29, 4880–4886 [47] Bonthuis, D. J.; Meyer, C.; Stein, D.; Dekker, C. Conformation and Dynamics of DNA Confined in Slitlike Nanofluidic Channels. Phys. Rev. Lett., 2008. 101, 108303 [48] Kowalczyk, S. W.; Grosberg, A. Y.; Rabin, Y.; Dekker, C. Modeling the Conductance and DNA Blockade of Solid-State Nanopores. Nanotechnology, 2011. 22, 315101 [49] Aksimentiev, A.; Heng, J. B.; Timp, G.; Schulten, K. Microscopic Kinetics of DNA Translocation Through Synthetic Nanopores. Biophys. J., 2004. 87, 2086–2097 [50] Heng, J. B.; Ho, C.; Kim, T.; Timp, R.; Aksimentiev, A.; Grinkova, Y. V.; Sligar, S.; Schulten, K.; Timp, G. Sizing DNA Using a Nanometer-Diameter Pore. Biophys. J., 2004. 87, 2905–2911 [51] Heng, J. B.; Aksimentiev, A.; Ho, C.; Marks, P.; Grinkova, Y. V.; Sligar, S.; Schulten, K.; Timp, G. Stretching DNA using the electric field in a synthetic nanopore. Nano Lett., 2005. 5, 1883–1888 [52] Heng, J. B.; Aksimentiev, A.; Ho, C.; Marks, P.; Grinkova, Y. V.; Sligar, S.; Schulten, K.; Timp, G. The Electromechanics of DNA in a Synthetic Nanopore. Biophys. J., 2006. 90, 1098–1106 [53] Luan, B.; Peng, H.; Polonsky, S.; Rossnagel, S.; Stolovitzky, G.; Martyna, G. Base-By-Base Ratcheting of Single Stranded DNA Through a Solid-State Nanopore. Phys. Rev. Lett., 2010. 104, 238103. WOS:000278608000009 [54] Luan, B.; Martyna, G.; Stolovitzky, G. Characterizing and Controlling the Motion of SsDNA in a Solid-State Nanopore. Biophys. J., 2011. 101, 2214–2222 [55] Luan, B.; Stolovitzky, G.; Martyna, G. Slowing and controlling the translocation of DNA in a solidstate nanopore. Nanoscale, 2012. 4, 1068–1077

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[56] Wells, D. B.; Belkin, M.; Comer, J.; Aksimentiev, A. Assessing Graphene Nanopores for Sequencing DNA. Nano Lett., 2012. 12, 4117–4123 [57] Yeh, I.-C.; Hummer, G. Nucleic Acid Transport Through Carbon Nanotube Membranes. Proceedings of the National Academy of Sciences of the United States of America, 2004. 101, 12177–12182 [58] Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids. Clarendon Press, New York, 1987 [59] Case, D. A.; Cheatham III, T. E.; Darden, T.; Gohlke, H.; Luo, R.; Merz Jr., K. M.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R. J. The AMBER Biomolecular Simulation Programs. J. Comput. Chem., 2005. 26, 1668 [60] Daura, X.; Mark, A. E.; van Gunsteren, W. F. Parametrization of Aliphatic CHn United Atoms of GROMOS96 force field. J. Comput. Chem., 1998. 19, 535–547 [61] Daura, X.; Suter, R.; van Gunsteren, W. F. Validation of Molecular Simulation by Comparison with Experiment: Rotational Reorientation of Tryptophan in Water. J. Chem. Phys., 1999. 110, 3049 [62] Møller, C.; Plesset, M. Note on an approximation treatment for many-electron systems. Phys. Rev., 1934. 46, 618–622 [63] Bartlett, R. J.; Silver, D. M. Correlation energy in Lithium Hydride, Boron Monohydride, and Hydrogen Fluoride with Many-Body Perturbation Theory Ssing Slater-Type Atomic Orbitals. Int. J. Quantum Chem., 1974. 8, 271–276 [64] Bartlett, R. J.; Silver, D. M. Diagrammatic Perturbation Theory. Int. J. Quantum Chem., 1975. 9, 183–198 [65] Binkley, J. S.; Pople, J. A. Moeller-Plesset Theory for Atomic Ground State Energies. Int. J. Quantum Chem., 1975. 9, 229–236 [66] Bartlett, R. J.; Silver, D. M. Correlation Energy in LiH, BH, and HF with Many-Body Perturbation Theory Using Slater-Type Atomic Orbitals. Int. J. Quantum Chem., 1976. 10, 185–186 [67] Hassanali, A. A.; Singer, S. J. A Model for the Water/amorphous Silica Interface: The Undissociated Surface. J. Phys. Chem. B, 2007. 111, 11181–11193 28

ACS Paragon Plus Environment

Page 28 of 39

Page 29 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[68] Hassanali, A. A.; Zhang, H.; Knight, C.; Shin, Y. K.; Singer, S. J. The Dissociated Amorphous Silica Surface: Model Development and Evaluation. J. Chem. Theory Comput., 2010. 6, 3456–3471 [69] van Beest, B. W. H.; Kramer, G. J.; van Santen, R. A. Force Fields for Silicas and Aluminophosphates Based on Ab Initio Calculations. Phys. Rev. Lett., 1990. 64, 1955–1958 [70] Todorov, I. T.; Smith, W.; Trachenko, K.; Dove, M. T. DL poly 3: New Dimensions in Molecular Dynamics Simulations Via Massive Parallelism. J. Mater. Chem., 2006. 16, 1911–1918 [71] Smith, W.; Todorov, I. T. A Short Description of DL POLY. Mol. Simul., 2006. 32, 935–943 [72] Todorov, I.; Smith, W. The DL POLY-4 User Manual, 2014 [73] van der Spoel, D.; Lindahl, E.; Hess, B. Gromacs user manual version 4.6. 2013 [74] Wang, J.; Cieplak, P.; Kollman, P. A. How Well Does a Restrained Electrostatic Potential (RESP) Model Perform in Calculating Conformational Energies of Organic and Biological Molecules? J. Comput. Chem., 2000. 21, 1049–1074 [75] Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of Multiple Amber Force Fields and Development of Improved Protein Backbone Parameters. Proteins: Struct., Funct., Bioinf., 2006. 65, 712–725 [76] Sorin, E. J.; Pande, V. S. Exploring the Helix-Coil Transition Via All-Atom Equilibrium Ensemble Simulations. Biophys. J., 2005. 88, 2472–2493 [77] DePaul, A. J.; Thompson, E. J.; Patel, S. S.; Haldeman, K.; Sorin, E. J. Equilibrium Conformational Dynamics in An RNA Tetraloop From Massively Parallel Molecular Dynamics. Nucleic Acids Research, 2010. 38, 4856–4867. PMID: 20223768 [78] Wing, R.; Drew, H.; Takano, T.; Broka, C.; Tanaka, S.; Itakura, K.; Dickerson, R. E. Crystal Structure Analysis of a Complete Turn of B-DNA. Nature, 1980. 287, 755–758 [79] Huff, N. T.; Demiralp, E.; Cagin, T.; Goddard III, W. A. Factors Affecting Molecular Dynamics Simulated Vitreous Silica Structures. J . Non-Cryst. Solids, 1999. 253, 133–142 29

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[80] Lyklema, J. Fundamentals of Interface and Colloid Science, volume II. Solid-Liquid Interfaces. Academic, New York, 1995 [81] Hoogerheide, D. P.; Garaj, S.; Golovchenko, J. A. Probing Surface Charge Fluctuations with SolidState Nanopores. Phys. Rev. Lett., 2009. 102, 256804 [82] Tadros, T.; Lyklema, J. Adsorption of Potential-Determining Ions At the Silica-Aqueous Electrolyte Interface and the Role of Some Cations. J. Electroanal. Chem., 1968. 17, 267–275 [83] Torrie, G. M.; Valleau, J. P. Nonphysical Sampling Distributions in Monte Carlo Free-Energy Estimation: Umbrella Sampling. J. Comput. Phys., 1977. 23, 187–199 [84] Kumar, S.; Rosenberg, J. M.; Bouzida, D.; Swendsen, R. H.; Kollman, P. A. The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules. I. The Method. J. Comput. Chem., 1992. 13, 1011–1021 [85] Boczko, E. M.; Brooks, C. L. Constant-Temperature Free Energy Surfaces for Physical and Chemical Processes. J. Phys. Chem., 1993. 97, 4509–4513 [86] Boczko, E. M.; Brooks, C. L. First-Principles Calculation of the Folding Free Energy of a Three-Helix Bundle Protein. Science, 1995. 269, 393–396. PMID: 7618103 [87] Souaille, M.; Roux, B. Extension to the Weighted Histogram Analysis Method: Combining Umbrella Sampling with Free Energy Calculations. Comput. Phys. Commun., 2001. 135, 40–57 [88] Flyvbjerg, H.; Petersen, H. G. Error Estimates on Averages of Correlated Data. J. Chem. Phys., 1989. 91, 461–466 [89] Murashov, V. V.; Leszczynski, J. Adsorption of the Phosphate Groups on Silica Hydroxyls: An Ab Initio Study. J. Phys. Chem. A, 1999. 103, 1228–1238 ˇ [90] Sponer, J.; Jureˇcka, P.; Hobza, P. Accurate Interaction Energies of Hydrogen-Bonded Nucleic Acid Base Pairs. J. Am. Chem. Soc., 2004. 126, 10142–10151

30

ACS Paragon Plus Environment

Page 30 of 39

Page 31 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[91] Cornell, W. D.; Cieplak, P.; Bayly, C. I.; Gould, I. R.; Merz, K. M.; Ferguson, D. M.; Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. J. Am. Chem. Soc., 1995. 117, 5179–5197 [92] Cieplak, P.; Cheatham, T. E.; Kollman, P. A. Molecular Dynamics Simulations Find That 3’ Phosphoramidate Modified DNA Duplexes Undergo a BtoA Transition and Normal DNA Duplexes An AtoB Transition. J. Am. Chem. Soc., 1997. 119, 6722–6730 [93] Cheatham, T. E.; Young, M. A. Molecular Dynamics Simulation of Nucleic Acids: Successes, Limitations, and Promise. Biopolymers, 2000. 56, 232–256 [94] Young, M. A.; Ravishanker, G.; Beveridge, D. L. A 5-Nanosecond Molecular Dynamics Trajectory for B-DNA: Analysis of Structure, Motions, and Solvation. Biophys. J., 1997. 73, 2313–2336 [95] Reddy, S. Y.; Leclerc, F.; Karplus, M. DNA Polymorphism: A Comparison of Force Fields for Nucleic Acids. Biophys. J., 2003. 84, 1421–1449 [96] Mao, Y.; Daniel, L. N.; Whittaker, N.; Saffiotti, U. DNA Binding to Crystalline Silica Characterized by Fourier-Transform Infrared Spectroscopy. Environ. Health Perspect., 1994. 102, 165–171 [97] Nandi, N.; Bagchi, B. Dielectric Relaxation of Biological Water. J. Phys. Chem. B, 1997. 101, 10954–10961 [98] S. K. Pal and J. Peon and A. H. Zewail. Biological Water At the Protein Surface: Dynamical Solvation Probed Directly with Femtosecond Resolution. Proc. Nat. Acad. Sci. U.S.A., 2002. 99, 1763–1768 [99] S. K. Pal and J. Peon and B. Bagchi and A. H. Zewail. Biological Water: Femtosecond Dynamics of Macromolecular Hydration. J. Phys. Chem. B, 2002. 106, 12376–12395 [100] Pal, S.; Zewail, A. Dynamics of Water in Biological Recognition. Chem. Rev., 2004. 104, 2099 [101] Bazzi, A.; Zargarian, L.; Chaminade, F.; Boudier, C.; Rocquigny, H. D.; Ren´e, B.; M´ely, Y.; Foss´e, P.; Mauffret, O. Structural Insights into the CTAR DNA Recognition by the HIV-1 Nucleocapsid Protein: Role of Sugar Deoxyriboses in the Binding Polarity of NC. Nucleic Acids Research, 2011. 39, 3903–3916 31

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 1: Lennard-Jones potential parameters of siloxane (OX ), silanol (OH ), and dissociated silanol (OM ) oxygens that were optimized to best reproduce ab initio data from fragment calculations. The charge parameters are derived from the BKS potential.69 silica atom types charge (e) σ(Å) ǫ(eV) OX (siloxane) −1.2 3.02 4.691 × 10−2 OH (silanol) −1.2 3.06 9.086 × 10−3 OM (dissociated silanol) −1.6 2.95 2.352 × 10−2 Si 2.4 – 0 H (silanol) 0.6 0 0

Table 2: Binding statistics for 8 independent simulations for a ssDNA oligomer, 8 independent simulations for a dsDNA oligomer. All of the dsDNA simulations resulted in binding to an undissociated silica surface, while 7 of the 8 ssDNA runs led to binding.

ssDNA dsDNA

binding events 7 8

phosphate 19 11

binding sites base (hydrophobic) deoxyribose sugar 7 4 1 1

32

ACS Paragon Plus Environment

Page 32 of 39

Page 33 of 39

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 34 of 39

Page 35 of 39

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(a)

(b)

Figure 5: Guanine and a cytosine of ssDNA oligomer bound to hydrophobic regions of silica surface. (a) and (b) are the same picture but DNA is invisible in (b). Larger red surface atoms are dissociated silanols and white dots are hydrogen atoms. In both snapshots, water and counter-ions are not shown.

(a)

(b)

Figure 6: A cytosine of dsDNA is bound to hydrophobic region of silica surface. (a) and (b) are the same picture but DNA is invisible in (b). Large red spheres are dissociated silanols and white dots are hydrogen atoms.

35

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 36 of 39

Page 37 of 39

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 11: Hydrogen bonding for deoxyribose sugar of ssDNA.

38

ACS Paragon Plus Environment

Page 38 of 39

Page 39 of 39

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment