Interface water-induced hydrophobic carbon chain unfolding in water

Zhang Xie,Zheng Li,Gang Lou,Qing Liang,Jiang-Xing Chen,Jianlong Kou,*and Gui-Na Wei

1 Institute of Condensed Matter Physics,Zhejiang Normal University,Jinhua 321004,China

2 School of Petroleum Engineering,China University of Petroleum (East China),Qingdao 266580,China

3 Department of Physics,Hanghzou Dianzi University,Hangzhou 310018,China

4 Division of Nephrology,The Second Affiliated Hospital of Zhejiang University,School of Medicine,Hangzhou 310009,China

Abstract The folding and unfolding of the carbon chain,which is the basic constitutional unit of polymers,are important to the performance of the material.However,it is difficult to regulate conformational transition of the carbon chain,especially in an aqueous environment.In this paper,we propose a strategy to regulate the conformational transition of the carbon chain in water based on the all-atom molecular dynamics simulations.It is shown that the unfolded carbon chain will spontaneously collapse into the folded state,while the folded carbon chain will unfold with an external electric field.The regulation ability of the electric field is attributed to the electric field-induced redistribution of interface water molecules near the carbon chain.The demonstrated method of regulating conformational transition of the carbon chain in water in this study provides an insight into regulating hydrophobic molecules in water,and has great potential in drug molecule design and new polymer material development.

Keywords: carbon chain unfolding,electric field,interface water,molecular dynamics simulations

1.Introduction

The regulation of conformational transition of the carbon chain in aqueous solutions is of fundamental importance in various fields such as biology,chemistry and physical processes,[1-5] because whether the carbon chain is folded or unfolded can greatly influence the properties of complex carbon products [6-9].Understanding the mechanism behind this regulation behavior is a prerequisite for its application and popularization [10-15].

Figure 1.Two typical initial conformations without hydrogen atoms at carbon chain end,(a)the folded state and(b)the unfolded state.The carbon atoms are depicted by spheres in cyan.The water molecules are shown as cloud-like.Gyration radius(Rg)of the hydrophobic carbon chain as a function of simulation time without (black) and with (red) electric field for using (c) initial conformation (a),and (d) initial conformation (b),respectively.The dashed blue lines are the critical Rgof the unfolded and folded carbon chain.

As the simplest hydrophobic molecules and the constitutional unit of many complex soft matters,like the polymer chain,[16] the carbon chain can largely affect the molecule performance by folding or unfolding [17,18].Generally,the preparation of macromolecules containing carbon chains is conducted in aqueous solutions [19,20].In aqueous solutions,the carbon chain usually has low solubility because carbon atoms are usually regarded as hydrophobic [21,22].The unfolded carbon chain will spontaneously collapse into the folded state [21,22].The driving force provided by the hydrophobic interaction has been considered to be one of the main reasons for the self-assembly of the carbon chain [23].Similarly,some proteins and DNAs will also be folded driven by the interaction between its internal hydrophobic groups[24-26].To obtain an unfolded carbon chain,different methods were proposed.[27-34]P.Leuenberger et al[27]systematically analyzed the sensitivity of the protein,indicating that higher temperature can make the protein unfolded.S.N.Innes-Gold et al [28] found that a charged polymer would change shape continuously by increasing ion concentrations in the solution.R.Zangi et al [29] showed by molecular dynamics (MD) simulations that a urea solution can change the folded carbon chain into an unfolded one.R.Godawat et al [30] found that adding guanidine chloride to the aqueous solution could overcome the hydrophobic interaction of the carbon chain,and the conformation would change.However,the introduced molecules are difficult to recycle after the system reaches stability,so they will affect the properties of the soft matter.

Electric field is an effective and controllable method to control the properties of materials[35-43].Water confined in a nanochannel can freeze with a weak electric field,[35] and it will melt with a strong electric field[36].A gradient electric field could drive nonpolar nanoparticles in water motion in one direction [39] The orientation of graphene and carbon nanotubes in water can be controlled by an electric field[40,41].Based on these insights,here we study conformational transition of the carbon chain in water using molecular dynamics simulations.We demonstrate that the conformational transition of the carbon chain in water can be regulated by applying an external electric field,i.e.,the carbon chain in water can spontaneously collapse into the folded conformation,while the folded carbon chain will unfold with an external electric field.

2.Methodology and simulation methods

It has been reported that the carbon chain is the basic unit of polymer materials [16].Here,various carbon chains are prepared,i.e.,40 carbon atoms with/without hydrogen atoms,and 50 carbon atoms with/without hydrogen atoms.For the carbon chains without hydrogen atoms,carbon atoms are like beads that string them together in a chain.For the carbon chains with hydrogen atoms,the chain is like a saturated hydrocarbon molecule with hydrogen atoms on each carbon atom.Figure 1 shows the simulation framework,which consists of one carbon chain of 40 carbon atoms without hydrogen atoms and 7121 water molecules.Two typical initial conformations,i.e.,folding topology and unfolding topology,were prepared,as shown in figures 1(a)and(b).To investigate the conformational transition of the carbon chain,we chose an electric field with strength below 2.5V nm-1,which is based on the previous reports that water molecules could be instantaneously dissociated when the electric field strength was greater than 3.5 V nm-1[44].In addition,the selected electric field strength is comparable to that of water molecules near the surfaces of biopolymers or in the cracks of amino acid crystals [45,46].Similar electric field strengths have been used successfully for studying water transportation,biological gating and DNA detection [35-41,43].

In our simulations,the OPLS-AA force field was used.For the carbon chain without hydrogen atoms,the carbon atoms in the chain are thought to be uncharged.The uncharged carbon atoms were modeled as Lennard-Jones(L-J)particles with εCC=0.276 kJ mol-1and σCC=0.35 nm.The two-carbon bond was 0.153 nm and the three-carbon angle was 112.7° before relaxation,which were maintained by harmonic potentials with spring constants of 224 262.4 kJ mol-1nm-2and 488.273 kJ mol-1deg-2,respectively.For the carbon chain with hydrogen atoms,the carbon atoms in the chain are thought to be charged and the whole chain is electrically neutral.The detailed force-field parameters can be found in our paper [47].The van der Waals interactions were described with L-J potential with a cut-off distance of 1.4 nm,and particle mesh Ewald method was used to calculate the electrostatic interaction,in which the mesh spacing of fast fourier transform was 0.12 nm.All simulations were carried out by using the Gromacs-5.1.5 molecular dynamics program in the NPT ensemble [48].The TIP4P water model was used,[49]which has been used successfully to study phase transition of water with an electric field,[36,50,51]and the TIP3P water model was also tested [52].The pressure was controlled by Parrinello-Rahman scheme to 1.0 bar with a compressibility coefficient of 4.5×10-5bar-1[53].The V-rescale thermostat was used to keep the temperature of the system at 300 K [54].The periodic boundary conditions were employed in all directions.For each simulation,a total of 60 ns was performed,and the data of the last 40 ns were collected for analysis.For each condition,three independent simulations were performed.

The free energy landscape can be determined by calculating the probability distribution of carbon chain conformation [55] P(x)=exp[-RT G(x)]/Z,where x is the conformation of the carbon chain and Z is the partition function.The relative free energy can be expressed as ΔG=G(x2) -G(x1) = -Here,we used the gyration radius (Rg) and root mean square deviation(RMSD) to describe the conformation of the carbon chain.

3.Results and discussion

We performed molecular dynamics to study conformational transition of the carbon chain.In order to reduce the effect of charge of the carbon chain on the conformational transition of the carbon chain with an electric field,we first selected the carbon chain model of 40 carbon atoms without hydrogen atoms.For the convenience of discussion,the gyration radius(Rg) of the carbon chain is defined as the average distance between the coordinates of all atoms and the center of mass.When the Rgis smaller than 0.6 nm,the conformation is considered to be in a folded state,as shown in figure 1(a).The conformation is regarded as an unfolded state when the Rgis greater than 1.1 nm,as shown in figure 1(b).The dashed blue lines in figures 1(c) and (d) point out the critical Rgof the unfolded and folded carbon chain.

Figures 1(c)and(d)show the Rgof the carbon chain as a function of time in water without and with an electric field for different initial conformations.When a folded carbon chain is used as initial conformation(see figure 1(a)),the Rgis almost unchanged without electric field,as shown in figure 1(c)(black line)or movie 1 in the Supporting Information(stacks.iop.org/CTP/73/055602/mmedia).This indicates that the folded carbon chain is still in a folded state,which is consistent with previous findings [29].This is a result of the competing balance of the carbon-carbon interaction,the rigidity of the carbon chain and the carbon-water interaction[29].When an electric field is applied to the system,the Rgof the carbon chain increases rapidly,exceeding the critical value of 1.1nm,as shown in figure 1(c)(red line)or movie 2 in the Supporting Information (stacks.iop.org/CTP/73/055602/mmedia).This means that the folded carbon chain is unfolded.A longer simulation time of 200 ns was also performed,as shown in figure S1 of the Supporting Information.The trend of folding and unfolding will not change with increasing simulation time.This confirms that the folded carbon chain is still in the folded state without electric field,and the folded carbon chain can be unfolded with electric field.To check the reliability of the initial conformation of the carbon chain,an initially unfolded carbon chain (figure 1(b))is used in the molecular dynamics simulations.In the absence of the electric field,the unfolded carbon chain turns into the folded state spontaneously.When an electric field is applied to the system,the unfolded carbon chain will remain unfolded,as shown in figure 1(d).Moreover,we performed the molecular dynamics simulations with the TIP3P water model[52].The carbon chain exhibits similar behaviors of folding and unfolding,as shown in figure S2.This suggests that the electric field can unfold the carbon chain in water.

We calculated the radial distribution,g(r),of the oxygen atom in water molecules around the proximal carbon atoms of the carbon chain,as shown in figure S3.Water molecules exhibit a layer-like distribution along the perpendicular direction to the carbon chain.At 0.5 nm from the carbon chain,there is a highest-density water layer.The region from the carbon chain to 0.5 nm is named the interface-region.Other regions are named as bulk-region.We calculated the number of water molecules in the interface-region with varying electric field strengths,as shown in figure 2(a).With increasing electric strength,more water molecules appeared in the interface-region.This indicates that water molecules in the interface-region redistribute with an electric field and the surface of the carbon chain becomes hydrophilic.

The average energies of carbon-water interactions,carbon-carbon interactions and the rigidity of the carbon chain are important to conformation transition of the carbon chain[24,25].We fixed the configuration of the carbon chain and calculated the carbon-water interaction energy for the folded carbon chain under different electric field strengths,and the results are illustrated in figure 2(b).In the absence of an electric field,the carbon-carbon interaction is greater than the sum of the rigidity of the carbon chain and the carbon-water interaction,leading to the folding of the carbon chain.When an electric field is applied to the system,the carbon-water interaction increases [41].This means that the surface of carbon-water becomes more hydrophilic,and the folded chain will be opened to ensure maximum contact area between the water molecules and the surface of the hydrophilic chain.In particular,the sum of rigidity of the carbon chain and the carbon-water interaction is greater than the carbon-carbon interaction,and the carbon chain will break away from selfrestraint to open.

Figure 2.Electric field dependence of(a)the average number(N)of water molecules in the interface-region,and(b)the average energy(Pcw)of the carbon-water interaction.

Figure 3.(a) Electric field dependence of the average number of hydrogen bonds of water molecules in interface-region.(b) Schematic diagrams of(I)a bending carbon chain with an angle α,(II)spatial orientation of a water molecule and the electric field,a water molecule on the surface of carbon chain (III) without electric field and (IV) with an electric field.

We calculated the average number of hydrogen bonds of a water molecule(NHbond)in the interface-region,as shown in figure 3(a).A hydrogen bond is formed when the distance of O-O is less than 0.35 nm and the angle O-H-O is less than 30°between two water molecules [58].It can be see that the NHbondincreases with increasing electric field strength.In the absence of an electric field,the application of a hydrophobic carbon chain destroys the hydrogen bond network between water molecules.The hydrophobic carbon chain will collapse to reduce the interface area between the water molecules and carbon chain and the destruction of the hydrogen bond[56,57].When an electric field is applied,the arrangement of water molecules in the interface-region will change,resulting in the change of the hydrogen bond network in the interfaceregion.The increase of hydrogen bonds from figure 3(a)indicates that the hydrophobic interaction is weakened,which means that the surface of the carbon chain is more hydrophilic.Therefore,the carbon chain will be open to obtain a larger carbon-water interface area,which is consistent with the above discussion.We calculated the solvent accessible surface area (SASA) of the folded carbon chain without and with the electric field,as shown in figure S4 of the Supporting Information.We can see that the SASA keeps a small value without the electric field,while the SASA keeps a large value with an electric field.

To find out why the carbon chain breaks away from selfrestraint to an open chain with direction and what happens to the interaction between water molecules and the carbon chain with the electric field,we select one part of the carbon chain with specific bending angle θ for analysis,as shown in figure 3(b)(I).We define the dipole(θ)as the angle formed by the dipole moment of the water molecule and the reverse direction of the electric field,as shown in figure 3(b)(II).The direction of the carbon chain is defined as the direction of the connection between the first and the last carbon atoms.For the electric field-free case at the bend location,the oxygen atom of the water molecule is toward the carbon atom of the chain,while the hydrogen atoms of the water molecule are faraway from the carbon atom,as shown in figure 3(b)(III).When an electric field is applied,the water molecule will obtain a torque as,M=E×P,[40,62]where P is the dipole moment of water molecule and E is the applied electric field.The torque can be rewritten as,M=-P Esinθ.The dipole direction of the water molecule will turn to the direction of the electric field [40,41,59,62].In addition,the electrostatic energy of the dipole moment is calculated as W=-P·E.The electrostatic energy of the dipole moment can be rewritten as,W=P E cosθ.When the angle θ is 180°,the electrostatic energy is the minimum,as shown in figure 3(b)(IV).In other words,when the direction of the dipole moment is along the direction of the electric field,the water molecule is the most stable.Applying an electric field will cause the rotation of the water molecule dipole,which induces a torque,resulting in the extension of the carbon chain.

Figure 4.Dynamics of water molecules under different electric field strengths.(a) The snapshots of distribution of water molecules of unfolded carbon chain and folded carbon chain in the interface region.Inset arrow indicates the direction of the electric field.The average probability distribution of the dipole of water molecules in (b) the bulk-region,and (c) the interface-region.

To understand unfolding carbon chain with an electric field in water,we calculated the dipole of water molecules.A folded carbon chain and an unfolded carbon chain are selected and kept in initial state,as shown in figure 4(a).Figures 4(b)and (c) show the dipole distribution of water molecules both the bulk-region and the interface-region in different electric fields.It is well known that water molecules flip along the direction of the electric field[40,41,59-62].We can see from the half height width of figures 4(b) and (c) that the dipole direction of water molecule is closer to the direction of electric field with a larger electric field strength (or see figure 4(a)),which is consistent with the discussion in the previous paragraph.For the bulk-region,whether the carbon chain is folded or not,the water molecules keep the same dipole distribution,as shown in figure 4(b).This indicates that carbon chain can not affect the dipole distribution of water molecules in the bulk-region,only the dipole direction of water molecule is induced by the electric field.For the interface-region,the dipole distribution of water molecules of the folded carbon chain is different from the unfolded carbon chain,as shown in figure 4(c).When the electric field strength is equal to zero,the dipole distribution of water molecules in folded carbon chain is higher than that of unfolded carbon chain.Therefore,there is a tendency that the dipole distribution of water molecules in the unfolded carbon chain system changes to that of the folded carbon chain system.As a result,the unfolded carbon chain is folded spontaneously.When the electric field is applied,the dipole distribution of water molecules of the unfolded carbon chain is higher than that of the folded carbon chain,as shown in figure 4(c).Comparing figures 4(b) and (c),we also see that the dipole distribution of the unfolded carbon chain in the interfaceregion is closer to that of the bulk-region.Based on the previous knowledge,the electric field can unfold the carbon chain in water,i.e.,the unfolded carbon chain is more stable.In other words,there is a tendency that the dipole distribution of water molecules in the folded carbon chain changes to that of the unfolded carbon chain.In this way,rotating the dipole direction of the water molecules will apply a force on the folded carbon chain to make it turn,so that the conformational transition of the carbon chain will take place from the folded state to the unfolded state.

Figure 5.The probability distribution of (a) the average angle between the direction of the electric field and the direction of carbon chain,(b) the gyration radius of the carbon chain with different electric field strength.

Furthermore,we investigate the relationship between the electric field direction and the carbon chain direction.The probability distribution of the angle between the carbon chain and the electric field is calculated,as shown in figure 5(a).The angle between the carbon chain and the electric field is defined as the angle between the direction of the connection between the first and the last carbon atoms and the electric field.We find that the angle between the carbon chain and the electric field is not completely the same under an electric field.There is a small angle between the carbon chain and the electric field.This means that the extended carbon chain will be folded at some time,which is consistent with figure 1 in that there are some fluctuations in the gyration radius spectrum,indicating that the extended carbon chain folds occasionally,or the folded carbon chain extends occasionally.This shows that the carbon chain in the water is not stable.When a larger electric field is applied,the angle will decrease,but is still not zero.To clarify the stability,we calculated the probability distribution of the gyration radii of the carbon chain with varying electric field strength,as shown in figure 5(b).It can be seen that when the electric field strength is lower than 0.5 V nm-1,the effect of the electric field on the carbon chain is almost negligible and the carbon chain will remain folded as for that without an electric field.When the electric field strength increases,the Rgwill be enlarged.In particular,when the electric field strength is higher than 1.5 V nm-1,the Rgis almost greater than 1.1 nm corresponding to the unfolded carbon chain.It can be understood that a larger electric field can more easily control the dipole orientation of water molecules and provide a higher torque.When the electric field strength is between 0.5 V nm-1and 1.5 V nm-1,the unstable carbon chain exhibits two maximum probabilities at the gyration radius of 0.50 nm and 1.25 nm,respectively.In such an electric field,the folded state and the unfolded state will transform to each other continuously.

To probe the physical forces governing the conformational transition of the carbon chain,we computed the free energy in the process of the transition.Although we have defined the gyration radius which can quantitatively describe the folding degree of the carbon chain,there are different conformations with the same gyration radius (Rg).We calculated root mean square deviation (RMSD) at a given gyration radius,as shown in figure S5 of the Supporting Information.Initial states of the unfolded conformation and folded conformation are taken as the references for figures S5(a) and (b),respectively.Figure 6 shows the free energy landscape of the carbon chain in different gyration radii corresponding to the RMSD without and with applied electric fields.To observe the spontaneous transformation process of the carbon chain,an unfolded carbon chain is used without electric field (inset in the lower-right corner of figure 6(a)),while a folded carbon chain is used with electric field(inset in the lower-left corner of figure 6(b)),and the initial state is taken as the reference in calculating RMSD and free energies.One gyration radius can correspond to different conformations with different free energies.It can be seen that the free energies display the opposite trend with variation of the gyration radius,as shown in figure 6.This is consistent with the discussion above that the extended carbon chain will be folded,while the folded carbon chain will be unfolded under the effect of the electric field.It is noted that without electric field(see figure 6(a)),the folded carbon chain has a lower free energy than that of the unfolded carbon chain.This means that the unfolded carbon chain prefers to be folded.Therefore,the unfolded carbon chain will change to the folded carbon chain spontaneously,which is consistent with the previous findings [29].On the contrary,when an electric field is applied to the folded carbon chain with low free energy(without electric field)(see figure 6(b)),the free energy of the folded carbon chain increases,while the free energy of the unfolded carbon chain decreases.As a result,the free energy of the folded carbon chain is higher than that of the unfolded carbon chain with an electric field.The relatively lower free energy of the unfolded carbon chain will allow it to remain unfolded.This is the reason why the folded carbon chain can be converted into the unfolded carbon chain by the electric field in water.

Figure 6.Free energy landscape of the carbon chain as a function of Rg and RMSD (a) without electric field and (b) with an electric field.Insets in the lower-right corner of figure 6(a) and inset in the lower-left corner of figure 6(b) are the initial conformations.

Figure 7.Snapshots of five carbon chains suspended in water in various field strengths.(a)Unfolded carbon chain conformation with electric field.(b) Folded carbon chain conformation without electric field.

Generally,a polymer material is usually composed of many carbon chains.Here,we also constructed five carbon chains,each of which is consistent with the above carbon chain.Initially,the five carbon chains are arranged side by side in the water.When the electric field is applied to the system,the five carbon chains will be arranged closely and keep the same direction as the electric field,as shown in figure 7(a).With the increase of the electric field strength,the carbon chains will unfold and gradually extend completely.Note that in the absence of an electric field,five carbon chains will fold spontaneously and become a ball,as shown in figure 7(b).If an electric field is applied to the system with folded conformation,all carbon chains will be unfolded again.This expected fast change of carbon chain conformation under electric field may be exploited to develop new rapid response materials.

Finally,it is significant to explore the universality of the electric field-induced carbon chain unfolding,which may bridge the connection with various applications.We have chosen the various conformations,e.g.the longer carbon chain with 50 carbon atoms,and the carbon chain with hydrogen atoms,as shown in figure S6.With increasing length of the carbon chain,the carbon chain is more stable with folding or unfolding; meanwhile,the carbon chain unfolding with an electric field and the carbon chain folding spontaneously can still be observed.The existence of hydrogen atoms on the carbon chain will increase its rigidity,resulting in a difference in the spontaneous carbon chain folding.These studies indicate that the electric field-induced carbon chain unfolding in water is rather universal.

4.Conclusions

In summary,we have investigated the conformational transition of the carbon chain in water by molecular dynamics simulations.It is found that the unfolded carbon chain can spontaneously collapse into the folded state,while the folded carbon chain can be unfolded by an external electric field.Different analyses are given to understand the mechanism of folding and unfolding of the carbon chain.The electric field can enhance the interaction between carbon and water,and reduce the free energy at the unfolded state of the carbon chain.This work presents a new method to regulate the carbon chain conformation,which has applications in the design of ordered materials,and can provide new insights into understanding the conformational transition of the hydrophobic molecule chain.The method with the electric field is safe for the conformational transition of the carbon chain.It does not introduce chemicals and avoids the permanent destruction of the polymer structure.In addition,if the applied electric field is removed,the unfolded polymer will recover to the folded state.

Acknowledgments

This work was partially supported by the Zhejiang Provincial Natural Science Foundation of China (No.LR21A020001)and the National Natural Science Foundation of China (No.11 774 313).J.K.was supported by the Start-up fund of Zhejiang Normal University.

Supporting Information

Long simulation time (figure S1),water model test (figure S2),radial distribution function(figure S3),solvent accessible surface area (figure S4),root mean square deviation (figure S5),and the reliability of the carbon chain length and hydrogen on carbon chain (figure S6).Movies of conformational transition of the carbon chain without the electric field(movie 1) and with the electric field (movie 2).

ORCID iDs