Exploring the Folding Pattern of a Polymer Chain in a Single Crystal

Mar 5, 2013 - Understanding the folding pattern of a single polymer chain within its single crystal .... Chinese Journal of Polymer Science 2014 32, 1...
0 downloads 0 Views 2MB Size
Letter pubs.acs.org/Langmuir

Exploring the Folding Pattern of a Polymer Chain in a Single Crystal by Combining Single-Molecule Force Spectroscopy and Steered Molecular Dynamics Simulations Yu Song,† Wei Feng,†,§ Kai Liu,‡ Peng Yang,† Wenke Zhang,*,† and Xi Zhang‡ †

State Key Laboratory of Supramolecular Structure and Materials, College of Chemistry, and §Institute of Atomic and Molecular Physics, Jilin University, Changchun 130012, P. R. China ‡ Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, P. R. China S Supporting Information *

ABSTRACT: Understanding the folding pattern of a single polymer chain within its single crystal will shed light on the mechanism of crystallization. Here, we use the combined techniques of atomic force microscopy (AFM)-based singlemolecule force spectroscopy (SMFS) and steered molecular dynamics (SMD) simulations to study the folding pattern of a polyethylene oxide (PEO) chain in its single crystal. Our results show that the folding pattern of a PEO chain in the crystal formed in dilute solution follows the adjacent re-entry folding model. While in the crystal obtained from the melt, the nonadjacent folding with large and irregular loops contributes to big force fluctuations in the force−extension curves. The method established here can offer a novel strategy to directly unravel the chain-folding pattern of polymer single crystals at single-molecule level.



unfolding of multiple polyproteins,9,11 the sequence-dependent unbinding of RNA−protein complexes,12 or detachment of multiple polymer loops from the substrate;8 long plateaucontaining force curves may reflect the unwinding of the helical structures13,14 or detachment of train-like polymer fragments from the substrate.8,15 This information has deepened our understanding of the relationship between molecular structure and mechanical properties of individual polymer chains. Apart from the great progress in single-molecule detection techniques, great advancement has also been made in the field of steered molecular dynamic (SMD) simulations, in which the force-induced conformation transition can be revealed at the atomic level, representing a powerful tool complementing and guiding single-molecule experiments in vitro.16−18 Here, we used AFM-based single-molecule force spectroscopy and steered molecular dynamics (SMD) simulations to explore the folding pattern of a polyethylene oxide (PEO) chain in its single crystal prepared by solution or melt crystallization. SMD simulations have been performed on adjacent re-entry19 and nonadjacent folds20,21 of PEO crystal

INTRODUCTION Analysis and characterization of the folding pattern of polymer chains in their single crystals is an exceedingly challenging task, mainly due to the long-chain nature of the polymer molecules leading to highly complex and various growth forms. Despite techniques such as neutron scattering1 and IR spectroscopy2 having been used to indirectly probe the chain folding pattern of polymer single crystals, single-molecule-level experimental verification of the chain folding pattern has only been gained so far in silico3 due to the complexity of the sample system and the limitation of detecting methods. With the appearance of those elegant single-molecule manipulation techniques, such as atomic force microscopy (AFM),4,5 optical tweezers,6 magnetic beads,7 and so on, the study of molecular interactions at the individual molecule level becomes possible. Among those single-molecule manipulation techniques, AFM can work in both aqueous solution and organic solvent, making it a powerful technique in the investigation of molecular interactions in biological and synthetic polymer systems.8−10 Fruitful information can be derived from the force spectrum of polymers. Due to the fact that some force spectra can be correlated well with the force-induced structure/conformation change of a polymer chain, the force spectrum can also be termed a force fingerprint of the specific polymer. For example, the sawtooth pattern in the force spectrum may represent the © 2013 American Chemical Society

Received: December 29, 2012 Revised: February 1, 2013 Published: March 5, 2013 3853

dx.doi.org/10.1021/la305157p | Langmuir 2013, 29, 3853−3857

Langmuir

Letter

PEO fragments are connected by small amorphous loops containing only a single repeat unit (for SMD simulations results on adjacent re-entry folds with bigger loops can be found in Figure S2 in Supporting Information) of PEO (−OCH2CH2−) (shown in Figure 1A), while in nonadjacent

systems. SMFS results on PEO single crystals grown from dilute solution and the melt phase were compared with those SMD simulations, and different folding patterns were identified for those two types of polymer crystals.



EXPERIMENTAL SECTION

Materials and Sample Preparation. Thiol-terminated polyethylene oxide (HS-PEO-OCH3, Mn = 48.5k g/mol, PDI = 1.05) was purchased from Polymer Source Inc. Amyl acetate and dimethylsulfoxide were purchased from Aldrich. N-(α-Maleimidoacetoxy)succinimide ester (AMAS) was purchased from Thermo Scientific (Pierce). 3-Aminopropyldimethylmethoxysilane was purchased from Fluorochem, UK. All materials were used as received. PEO single crystals were prepared from solution and melt by the self-seeding method, which was adapted from the previous work by Li et al.22 and Wang et al.23 The process of immobilization of the PEO single crystals has been described in detail in the Supporting Information and our previous work.24 SMD Simulations. The steered molecular dynamics simulations were performed in the NPT ensemble (1 atm press and 298 K) using the programs VMD1.8.7,25 NAMD2.7,26 and the CHARMM32 ether force field,27,28 with an integration time step of 1 fs. Van der Waals (vdW) energies were calculated using a smooth (12 Å) cutoff. We set up the crystal structure model by MS Modeling (the direction of crystal lattice c corresponding to z axis); 5 times length of one crystal lattice along z axis29 was used so the length was nearly equal to the thickness of the crystal as estimated in the AFM experiment (∼10 nm), and constrained atoms in the outside chains (along x, −x, y, and −y axis) to mimic the expanded crystal. After energy minimization and a 2 ns equilibration of the simulation system, we randomly fixed four terminal carbon atoms (as the fixed atom) at one side of the crystal, and pulled one chain of the opposite surface by attaching a harmonic spring (with a force constant of 7 kcal/mol/Å2) in between the terminal carbon (as the SMD atom) atom and the stretching point at a constant speed of 20 m/s along z axis. Single-Molecule Force Spectroscopy. All AFM experiments were carried out on a NanoWizardII BioAFM (JPK instrument AG, Berlin, Germany) using gold-coated Si3N4 AFM tips (Biolever, Olympus Japan) and Si3N4 tips (MSCT, Bruker Corp.). Cantilevers were calibrated in solution using the thermal noise method30 to obtain the spring constant (typically around 0.06 N/m for MSCT and 0.03 N/m for Biolever). All experiments were performed in amyl acetate at room temperature with a pulling velocity of 2 μm/s. The experimental details have been described in our previous work.24 For a random picking experiment, the same AFM tip can be used up to three times (i.e., extracting three individual molecules out of the single crystal) without significantly affecting the experimental result, while for the gold nanoparticle positioning experiment, the AFM tip can normally be used only once. All data were analyzed using custom software written in Igor Pro (Wavemetrics).

Figure 1. SMD simulations of adjacent re-entry and nonadjacent reentry. (A) Adjacent re-entry folding mode. The crystallized region (PEO fragments) is colored in red. The amorphous region (amorphous loops) colored in yellow contains only one single repeat unit of (−OCH2CH2−). (B) Nonadjacent re-entry folding mode. The crystallized region (PEO fragments) is colored in blue. The amorphous region (amorphous loops) contains 1 (yellow), 3 (black), and 5 (green) repeat units of (−OCH2CH2−). (C) Stretching of three folds of PEO fragments out of its single crystal. For the sake of clarity, the surrounding PEO fragments were removed. (D) Comparison of force−extension curves obtained in adjacent re-entry (colored in red) with nonadjacent re-entry (colored in blue) systems by constant velocity SMD simulations, in which 17 PEO stems were pulled out of the crystals. The arrows (a, b, c) show the unfolding of top loops that contain 5, 3, and 1 repeat units, respectively.



RESULTS AND DISCUSSION SMD Simulations on the Pulling of Adjacent Re-Entry and Nonadjacent Re-Entry Folding Chains Produce Different Force Patterns. To describe the folding patterns of polymer chains in their crystals, different models have been suggested, such as the adjacent re-entry folding,19 the nonadjacent folding (including the switchboard),20,21 and so on. There are many arguments on the correctness of those models. To test those suggested models, SMD simulations were performed by assuming ideal adjacent re-entry and nonadjacent re-entry folding for the PEO chains within the single crystal, respectively. It needs to be pointed out that although the SMD simulation was performed at high pulling speed (due to the limitation of computing power) this will not affect the accuracy of the pattern of force changes along the lateral pulling direction.16−18 For the adjacent re-entry folding, the adjacent

re-entry folding mode the nonadjacent PEO fragments are connected by the amorphous loops containing various number of repeat units 1, 3, and 5, respectively(shown in Figure 1B). Figure 1C shows the pulling of six stems of PEO fragments (three folds) out of its crystal, which produced three main wave-like peaks. The peak regions, which are relatively flat, correspond to the traveling of the single PEO stem and bottom loop within the crystal phase, as well as the unfolding of a fold (top left panel in Figure 1C; see also the Supporting 3854

dx.doi.org/10.1021/la305157p | Langmuir 2013, 29, 3853−3857

Langmuir

Letter

Information Movie). The valley regions, which are still well above zero, correspond to the unraveling of the loop at the amorphous region that locates on the top crystal surface. The patterns on the unfolding curves are quite regular for the adjacent re-entry folding mode, though more complicated nonadjacent re-entry folding mode is observed (see Figure 1D). On the valley regions, the unbinding force decreases (c→a shown in Figure 1D) with the increase of loop sizes (from 1 to 5) on the top surface of the crystal, while on the peak regions, the force is clearly enhanced by the bigger bottom loops. When these facts are put together, we can understand the bigger force fluctuations observed in nonadjacent re-entry folds, that is, the length of the loops located on both crystal surfaces significantly affects the regularity of patterns on the stretching curves. Chain Folding in Solution-Growth PEO Single Crystal Is Different from That in Melt-Crystallized One. Apart from SMD simulations, we have also used single-molecule force spectroscopy to experimentally investigate the folding patterns of PEO chain in single crystals prepared from the dilute solution and melt, respectively. To pull a single PEO chain out of its crystal, we have developed a method for the immobilization of the single crystal and for the attachment of the polymer chain end to the AFM tip as described in our previous study.24 In brief, thiol-labeled PEO single crystals, which were prepared from its dilute solution, were immobilized by one side on silicon substrate via a short cross-linker (N-(αmaleimidoacetoxy)succinimide ester, AMAS), leaving the other side of the crystal surface exposed for picking up by the AFM tip. The SMFS pulling experiments were performed on the crystal surface with a gold-coated AFM tip under the guidance of AFM imaging (Figure S1 show the AFM images of PEO single crystal prepared from dilute solution and the melt). In other words, by AFM imaging, a specific single crystal can be found, then the AFM tip is brought into contact with the specific area on the crystal surface then separated from it. During contact, the thiol end of the polymer chain can get attached to the AFM tip via gold−thiol chemistry. The PEO molecule can then be pulled out of the crystal during separation. Force−extension curves with a sawtooth-like plateau were obtained during the pulling experiments, as shown in Figure 2A. Such curves take about 90% in total. Since it is widely accepted that single crystals formed in dilute solution should fold more regularly, we intend to compare our SMFS result with those SMD simulations on adjacent re-entry folding systems. Due to the huge difference in stretching velocity between the SMD and AFM-SMFS experiments, it is impossible to quantitatively compare the forces (i.e., the absolute forces are not comparable). However, the patterns on the stretching curves are quite similar and the force plateaus are relatively flat (shown in Figure 2B), which indicate that the sizes of the loops on the crystal surfaces should be relatively small and regular, supporting the adjacent re-entry folding model. It needs to be pointed out, however, that the periodicity of the sawtooth pattern on stretching curves obtained in SMFS experiments is not as regular/obvious as that in SMD simulations. This may due to the following facts: (1) for the current SMD simulation on the pulling of adjacent re-entry folds (Figure 1), the folding mode is idealized, and the size of loops in the amorphous region is very small and uniform; (2) in real solution-crystallized single crystals, the polymer chains can adopt the adjacent re-entry folding mode, but the sizes of loops that connect the adjacent stems may show some variation;19 (3) the force fluctuation in the sawtooth pattern of the

Figure 2. Single-molecule force spectroscopy experimental results obtained on the PEO single crystal prepared from solution. (A) Typical stretching curves. (B) Comparison of force−extension curves obtained by SMFS on the PEO single crystal prepared from solution (colored in black) with that obtained by SMD simulations on the unfolding of adjacent re-entry (colored in red) PEO chain. The inset shows the adjacent re-entry folding model. (C) A typical stretching curve with deviation at short and long extensions. The inset shows irregular folding at the ends of the polymer chain.

stretching curve obtained by SMFS is of similar magnitude to that of thermal noise; as a result, some patterns are not recognizable. SMD simulation results of the pulling of adjacent re-entry folding but with bigger loop sizes (containing six repeat units of PEO) show that the loop-size change (from one to six) will affect the gap between adjacent peaks (i.e., the periodicity of the force pattern); please see Figure S2 in Supporting Information for more detail. This is because the originally bigger loop will get shortened to the minimum length before it enters the crystal phase during stretching as shown in the Supporting Information Movie. In addition, there are some force−extension curves (14% in 1662 curves) which show deviation from the force plateau at short and long extensions (shown in Figure 2A and C). This may be ascribed to irregular folding at the ends of the polymer chain in agreement with earlier assumptions.19,20 For the single crystal obtained from melt, PEO chains were individually pulled out of its single crystal, which was crystallized on the silicon substrate. Typical force curves show a plateau with very large fluctuations (see blue curves in Figure 3A,B) similar to the SMD simulations on irregularly folded systems (shown in Figure 1). Therefore, it is not difficult to imagine that part of the loops on such crystal surfaces should 3855

dx.doi.org/10.1021/la305157p | Langmuir 2013, 29, 3853−3857

Langmuir

Letter

be relatively large and irregular, supporting the nonadjacent reentry folding model. To find out the magnitude of fluctuation on the force plateau more accurately, more than 200 stretching curves were used to construct histograms of force distribution on the retraction baseline (the right region of the stretching curve after the force plateau, as shown in Figure 3A) and the force plateaus, respectively,10 as shown in Figure 3C. From the figure, we can see that the noise level on solution- and melt-crystallized samples is almost the same since the force distribution of the retraction baseline is indistinguishable (see left histograms marked with dotted line in Figure 3C), while the force distribution of the force plateau from the melt-crystallized sample is broader than that of the solution-crystallized one (see right solid histograms in Figure 3C) indicative of bigger fluctuations. By superimposing the center of the histogram of the baseline with the histogram of the force plateau, we can see more clearly the difference of the fluctuation magnitude. The force fluctuations on solution-crystallized sample (colored in red) are close to that of the baseline (i.e., the noise level), as shown in the Figure 3D top traces, while the fluctuation from the melt-crystallized sample is much bigger since the force distribution is much broader (see bottom traces in blue in Figure 3D). To put all these results together, we can estimate that the adjacent re-entry folding mode plays a predominant role in the crystal grown from solution and nonadjacent re-entry with large and irregular loops occurs in the crystal prepared from the melt.



CONCLUSIONS In summary, the folding pattern of a PEO chain in its single crystal has been explored for the first time at single molecule level by a good combination of single-molecule force spectroscopy and SMD simulations. The results demonstrate that adjacent re-entry folding is the predominant mode in the crystal prepared from dilute solution, whereas nonadjacent reentry with large and irregular loops on the crystal surfaces occurs in the crystal from the melt. We believe that this study has provided a new physical insight into the polymer crystallization and the combined techniques represent a novel strategy to directly investigate the folding pattern/mechanism in other polymer single crystals as well.



ASSOCIATED CONTENT

S Supporting Information *

Description of sample preparation, AFM images on single crystals prepared form solution and melt, together with SMD simulations on adjacent re-entry folds with different loops sizes. This material is available free of charge via the Internet at http://pubs.acs.org.

Figure 3. Comparison of single-molecule force spectroscopy experimental results obtained on PEO single crystal prepared from the solution and melt, respectively. (A) Typical stretching curves obtained on PEO single crystal prepared from the melt. The green arrows show the valley regions corresponding to unraveling of the big top loops. The red arrows show the peak regions corresponding to unfolding of large bottom loops. The inset shows the nonadjacent reentry folding mode with large loops on the top surface (colored in green) and on bottom surface (colored in red). (B) Direct comparison of typical stretching curves obtained on crystals prepared from the solution (colored in red) and the melt (colored in blue). (C) Histograms of force distributions of the baseline (the dotted curve on the left) and the force plateau region (the solid curve on the right). The results from solution-crystallized and melt-crystallized samples are marked in red and blue, respectively. (D) The center positions of the histograms (in (C) above) of baseline (noise level) are shifted rightward to get superimposed with that of the histogram from the force plateau. For clarity purposes, histograms for solution- and meltcrystallized samples were offset from each other.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was funded by the NSFC (20974039, 21221063, 91127031) and NSFC-DFG joint project (TRR61), and the National Basic Research Program (2013CB834503). The work was also supported by SRF for ROCS (SEM), the Program for 3856

dx.doi.org/10.1021/la305157p | Langmuir 2013, 29, 3853−3857

Langmuir

Letter

nonequilibrium to near equilibrium. Macromolecules 2006, 39, 324− 329. (24) Liu, K.; Song, Y.; Feng, W.; Liu, N. N.; Zhang, W. K.; Zhang, X. Extracting a single polyethylene oxide chain from a single crystal by a combination of atomic force microscopy imaging and single-molecule force spectroscopy: toward the investigation of molecular interactions in their condensed states. J. Am. Chem. Soc. 2011, 133, 3226−3229. (25) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graphics 1996, 14, 33−38. (26) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781− 1802. (27) Vorobyov, I.; Anisimov, V. M.; Greene, S.; Venable, R. M.; Moser, A.; Pastor, R. W.; MacKerell, A. D. J. Additive and classical drude polarizable force fields for linear and cyclic ethers. J. Chem. Theory Comput. 2007, 3, 1120−1133. (28) Lee, H.; Venable, R. M.; MacKerell, A. D. J.; Pastor, R. W. Molecular dynamics studies of polyethylene oxide and polyethylene glycol: hydrodynamic radius and shape anisotropy. Biophys. J. 2008, 95, 1590−1599. (29) Takahashi, Y.; Tadokoro, H. Structural studies of polyethers, (-(CH2)m-O-)n. X. crystal structure of poly(ethylene oxide). Macromolecules 1973, 6, 672−675. (30) Florin, E.-L.; Rief, M.; Lehmann, H.; Ludwig, M.; Dornmair, C.; Moy, V. T.; Gaub, H. E. Sensing specific molecular interactions with the atomic force microscope. Biosens. Bioelectron. 1995, 10, 895−901.

New Century Excellent Talents in University (NCET), the State Key Laboratory, and the High Performance Computing Center of Jilin University. We thank Prof. Zhigang Wang (JLU) for helpful discussions and suggestions on SMD simulations.



REFERENCES

(1) Sadler, D. Y.; Harris, R. A neutron scattering analysis: how does the conformation in polyethylene crystals reflect that in the melt? J. Polym. Sci. Polym. Phys. Ed. 1982, 20, 561−578. (2) Spells, S. J.; Keller, A.; Sadler, D. M. I.R. study of solution-grown crystals of polyethylene: correlation with the model from neutron scattering. Polymer 1984, 25, 749−758. (3) Mansfield, M. L. Monte carlo study of chain folding in meltcrystallized polymers. Macromolecules 1983, 16, 914−920. (4) Binig, G.; Quate, C. F.; Gerber, C. Atomic force microscope. Phys. Rev. Lett. 1986, 56, 930−933. (5) Rief, M.; Gautel, M.; Oesterhelt, F.; Fernandez, J. M.; Gaub, H. E. Reversible unfolding of individual titin immunoglobulin domains by AFM. Science 1997, 276, 1109−1112. (6) Ashkin, A.; Schütze, K.; Dziedzic, J. M.; Euteneuer, U.; Schliwa, M. Force generation of organelle transprot meared in vivo by an infrared laser trap. Nature 1990, 348, 346−348. (7) Smith, S. B.; Finzi, L.; Bustamante, C. Direct mechanical measurements of the elasticity of single DNA molecules by using magenetic beads. Science 1992, 258, 1122−1126. (8) Zhang, W. K.; Zhang, X. Single molecule mechanochemistry of macromolecules. Prog. Polym. Sci. 2003, 28, 1271−1295. (9) Li, H. B.; Cao, Y. Protein mechanics: from single molecules to functional biomaterials. Acc. Chem. Res. 2010, 43, 1331−1341. (10) Li, I. T. S.; Walker, G. C. Interfacial free energy governs single polystyrene chain collapse in water and aqueous solutions. J. Am. Chem. Soc. 2010, 132, 6530−6540. (11) Oesterhelt, F.; Oesterhelt, D.; Pfeiffer, M.; Engel, A.; Gaub, H. E.; Müller, D. J. Unfolding pathways of individual bacteriorhodopsins. Science 2000, 288, 143−147. (12) Liu, N. N.; Peng, B.; Lin, Y.; Su, Z. H.; Niu, Z. W.; Wang, Q.; Zhang, W. K.; Li, H. B.; Shen, J. C. Pulling genetic RNA out of tobacco mosaic virus using single-molecule force spectroscopy. J. Am. Chem. Soc. 2010, 132, 11036−11038. (13) Zhang, Q. M.; Lu, Z. Y.; Hu, H.; Yang, W. T.; Marszalek, P. E. Direct detection of the formation of V-amylose helix by single molecule force spectroscopy. J. Am. Chem. Soc. 2006, 128, 9387−9393. (14) Rief, M.; Clausen-Schaumann, H.; Gaub, H. E. Sequencedependent mechanics of single DNA molecules. Nat. Struct. Biol. 1999, 6, 346−349. (15) Geisler, M.; Netz, R. R.; Hugel, T. Pulling a single polymer molecule off a substrate reveals the binding thermodynamics of cosolutes. Angew. Chem., Int. Ed. 2010, 49, 4730−4733. (16) Sotomayor, M.; Schulten, K. Single-molecule experiments in vitro and in silico. Science 2007, 316, 1144−1148. (17) Carrion-Vazquez, M.; Li, H. B.; Lu, H.; Marszalek, P. E.; Oberhauser, A. F.; Fernandez, J. M. The mechanical stability of ubiquitin is linkage dependent. Nat. Struct. Biol. 2003, 10, 739−743. (18) Cui, S. X.; Yu, J.; Kühner, F.; Schulten, K.; Gaub, H. E. Doublestranded DNA dissociates into single strands when dragged into a poor solvent. J. Am. Chem. Soc. 2007, 129, 14710−14716. (19) Keller, A. Polymer crystals. Rep. Prog. Phys. 1968, 31, 623−704 and references therein. (20) Frank, F. C. General introduction. Discuss. Faraday Soc. 1979, 68, 7−13. (21) Yoon, D. Y.; Flory, P. J. Small-angle neutron scattering by semicrystalline polyethylene. Polymer 1977, 18, 509−513. (22) Li, B.; Li, C. Y. Immobilizing Au nanoparticles with polymer single crystals, patterning and asymmetric functionalization. J. Am. Chem. Soc. 2007, 129, 12−13. (23) Zhai, X. M.; Wang, W.; Zhang, G. L.; He, B. L. Crystal pattern formation and transitions of PEO monolayers on solid substrates from 3857

dx.doi.org/10.1021/la305157p | Langmuir 2013, 29, 3853−3857