Next Article in Journal
Systematic Studies on Anti-Cancer Evaluation of Stilbene and Dibenzo[b,f]oxepine Derivatives
Next Article in Special Issue
A MD Simulation Prediction for Regulation of N-Terminal Modification on Binding of CD47 to CD172a in a Force-Dependent Manner
Previous Article in Journal
Bioactive Compounds and Antioxidant Composition of Nut Bars with Addition of Various Edible Insect Flours
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 8
10.3390/molecules28083557
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Study on the Aggregation Behavior of Triton X Micelles with Different PEO Chain Lengths in Aqueous Solution

by
** Peng
,
**aoju Song
,
**n Li
,
Yongkang Jiang
,
Guokui Liu
,
Yaoyao Wei
* and
Qiying **a
*
School of Chemistry and Chemical Engineering, Linyi University, Linyi 276000, China
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(8), 3557; https://doi.org/10.3390/molecules28083557
Submission received: 11 March 2023 / Revised: 2 April 2023 / Accepted: 16 April 2023 / Published: 18 April 2023
(This article belongs to the Special Issue Advances in Molecular Modeling in Chemistry)
Browse Figures
Versions Notes

Abstract

:
The aggregation structure of Triton X (TX) amphiphilic molecules in aqueous solution plays an important role in determining the various properties and applications of surfactant solutions. In this paper, the properties of micelles formed by TX-5, TX-114, and TX-100 molecules with different poly(ethylene oxide) (PEO) chain lengths in TX series of nonionic surfactants were studied via molecular dynamics (MD) simulation. The structural characteristics of three micelles were analyzed at the molecular level, including the shape and size of micelles, the solvent accessible surface area, the radial distribution function, the micelle configuration, and the hydration numbers. With the increase of PEO chain length, the micelle size and solvent accessible surface area also increase. The distribution probability of the polar head oxygen atoms on the surface of the TX-100 micelle is higher than that in the TX-5 or TX-114 micelle. In particular, the tail quaternary carbon atoms in the hydrophobic region are mainly located at the micelle exterior. For TX-5, TX-114, and TX-100 micelles, the interactions between micelles and water molecules are also quite different. These structures and comparisons at the molecular level contribute to the further understanding of the aggregation and applications of TX series surfactants.

1. Introduction

Micelles aggregated from surfactant molecules have been widely used owing to their ability to dissolve hydrophobic compounds effectively. People’s interest in micelle solutions comes from their application potential as functionally molecular assemblies [1]. Different cationic, anionic, zwitterionic, and non-ionic surfactants have been widely studied. Compared with ionic surfactants, non-ionic surfactants are much less toxic and have more efficient surface active properties [2]. Among them, the series of Triton X (TX) non-ionic surfactants have been studied in depth and characterized [3]. Because of their unique molecular structure and amphiphilic nature, the micelles formed by TX surfactants have very flexible surface-active properties. The critical micelle concentration (CMC) of TX surfactant is very low, especially in aqueous solution. Different TX surfactants have been widely used in the fields of microbiology and biomedicine [4].
The TX-100 as one typical nonionic surfactant [5,6,7], has been extensively studied [8,9,10]. This surfactant consists of one hydrophilic chain of 9–10 ethylene oxide units linked to a benzene ring with an octyl chain [6]. Many researchers have explored the structural characteristics of TX-100 micelles in aqueous solution and their interaction with water by experimental techniques such as 2D NOESY NMR [11,12,13], pulsed field gradient NMR [8], solvent paramagnetic relaxation enhancement [7], fluorescence spectra [14,15], surface tension [15], light scattering [16], static and dynamic light scattering [15,17], turbidimetric method [18], small-angle X-ray scattering [13], quasi-elastic light scattering spectroscopy [19]. It is very important to study the hydration and the size of TX-100 micelles [20], because the effectiveness of micellar applications depends on their size and their effects on the properties of the solution [21]. Robson and Dennis [16], Paradies [13], Streletzky and Phillies [19], Phillies et al. [22] used different experimental methods to determine the hydration degree of TX-100 micelles. Other relevant researchers studied the size of the TX-100 micelle by measuring its radius. The aggregation characteristics of TX-100 micelles are also the key aspects that affect their potential applications [7]. Yuan et al. [12] revealed the dependence of the conformation of TX-100 micelles on the environment, and emphasized that the spatial arrangement of PEO chains in micelles is to a certain extent very compact. In addition, the self aggregation and supramolecular micelle structure of TX-100 surfactant molecules in aqueous solution determined by Denkova et al. [11] shows that TX-100 may aggregate in a double-layer or multilayer spherical conformation interlaced between PEO chains and octyl phenyl parts. In order to further explore the structural characteristics, Zhang et al. [7] proposed a more detailed aggregation mode of TX-100 micelles, pointing out that this micelle is more likely to be a multilayer staggered spherical micelle and octyl phenyl is partially dispersed in different layers.
With the development of the computer and algorithm, molecular dynamics (MD) simulation has been widely applied in scientific studies from a micro-perspective. MD simulation is a method of calculating and simulating the time evolution of a group of interacting atoms using Newtonian equations of motion. By solving classical Newton’s equations of motion, the motion trajectory of the system can be obtained. This trajectory includes a large number of samplings of molecular configurations, as well as information of the position and velocity of the particles. Further analysis of the obtained configuration ensemble can provide macroscopic observable properties (density, surface tension, solubility, viscosity, and thermodynamic information) and microscopic properties (intuitive structure, conformational distribution, degree of structural fluctuation, position and intensity of non-bonded interactions, particle spatial distribution) of the system. The description (force field) of the interactions between molecules and atoms is the key for this method to obtain reliable results. This calculation method shows a high degree of relevance in the detailed characterization of self-assembly related systems, including complementarity with experimental data, optimization of experimental design, and prediction of related properties of chemical systems. This may be expensive or difficult for experimentation [23]. Moreover, MD simulations have been effectively applied to investigate the self-assembly and interfacial adsorption of surfactant systems.
In addition to many experimental studies, some literature studies [24,25,26,27,28,29] further discussed the aggregation characteristics of TX-100 micelles at the molecular level by MD simulation. Among these literature studies, the simulation results of several of these literature studies [24,25,28,29] show that the shape and behavior of micelles mainly depend on the aggregation number of aggregates. An important application of TX-100 micelles is that the hydrophobic part can dissolve insoluble organic compounds [30]. A detailed understanding of the pure micelle system is a prerequisite for the study of organic matter and micelle composite systems [31]. In order to better understand this solubilization mechanism, researchers [26,27] explored the basic structure of micelles under a specific aggregation number. Farafonov et al. [26] first observed that the surface of the micelles was highly irregular when they were finally balanced in aqueous solution. Generally, when the temperature was 298 K, the aggregation number of the TX-100 micelle measured by experiments [8,19,22] was about 100. Therefore, Pacheco-Blas and Vicente [27] studied the morphology of TX-100 micelles with an aggregation number of 100 with MD simulation. It was observed that the micelles appeared as a quasi-sphere. The hydrophobic parts of the surfactants preferentially gather in the core of the micelles; however, there is some probability that these hydrophobic groups may be located at the outer region of the micelle.
Micelle structure in surfactant solutions is an important factor affecting the properties and application fields of self-assembly [11]. In this paper, we used all-atom MD simulations to study the self-aggregation properties of TX-5, TX-114, and TX-100 surfactant systems. By analyzing micelle shape and size, the solvent accessible surface area (SASA) of the micelle, the micelle morphology, as well as the interactions between micelles and water molecules, the differences and connections of these three micelles could be compared.

2. Results and Discussion

Figure 1 shows the time evolutions of radius of gyration (Rg), the solvent accessible surface area (SASA), and the energy of the system during a simulation run to determine the equilibriums of the systems. As can be seen from Figure 1a,b, the Rg and SASA of three micelles tend to be stable after 150 ns. The total energies of the three systems hardly change during the whole process. These changes of Rg, SASA, and total energies indicated that our simulated systems all reach equilibrium after 150 ns, and the last 30 ns trajectory of each system was used for result analyses in this paper.

2.1. Micelle Shape and Size

2.1.1. Micelle Shape

The ratio of the moments of inertia Imax/Imin is usually used to characterize the shape of the micelle, where Imax and Imin are the largest and the smallest moment of inertia along the x, y, or z axis, respectively. When Imax/Imin is equal to 1, it is a perfect spherical micelle [32]. According to the results obtained in Table 1, the ellipsoidal degree of TX-100 micelle is significantly greater than that of TX-5 and TX-114 micelles. The TX-114 micelle is slightly more spherical than the TX-5 micelle. In addition, eccentricity (e) is also an important indicator to judge the shape of the micelle, which is defined as e = 1−Imin/Iavg [32]. The Iavg is the average moment of inertia of micelle. When e is equal to 0, the micelle shows a perfect spherical shape. When micelle shape is stable, the eccentricity will converge near a certain value and small fluctuations will appear [33]. In order to further verify the accuracy of the calculation results, we calculated the e values of the three micelles. The e of the TX-100 micelle is about 0.12, which is close to the value (0.08) of the TX-100 micelle obtained by Pacheco-Blas et al. [27]. The results show that the ellipsoidal degree of the TX-100 micelle is the largest, and the spheroidal degree of the TX-114 micelle is slightly larger than the TX-5 micelle.

2.1.2. Micelle Size

The size of the micelle is one of the important characteristics of micelle structure, in which Rg is a standard to characterize the size of the micelle [34]. It can be seen from Table 1 that Rg increases with the increase of the PEO chain length. In addition, the average micelle radius (Rs) is also used to judge the size of the micelle, and its definition [35] is as follows:
R s = 5 3   R g
According to the results obtained in Table 1, the change trend of Rs is the same as that of Rg. The radius value of the TX-100 micelle (≈2.83 nm) obtained in this paper is consistent with the radius value of the TX-100 micelle (≈2.87 nm) simulated by Pacheco-Blas et al. [27] and the experimental result (3.1 nm) obtained by Brown et al. [8]. Our calculated results are therefore considered to be reliable.

2.1.3. Solvent Accessible Surface Area

Solvent accessible surface area (SASA) is the surface area of the TX micelle in contact with solvent, which can well reflect the properties of the micelle surface. To calculate the SASA of the micelle, the double cubic lattice method [36] was applied. With this method, all water molecules are first removed from the system [35], and then one spherical probe molecule with a radius of 1.4 Å is used to simulate the rolling of water molecules on the surface of the micelle [37]. Table 1 lists the total SASA and the hydrophilic and hydrophobic SASA values of the three micelles. We find that all SASA values increase with the increase of PEO chain length, which is consistent with the increasing trend of micelle size. The hydrophobic surface area of each micelle is much larger than the hydrophilic surface area, and the proportion of the hydrophobic area increases with the rise of the polar head chain. This is on account of the polar head of the research system being polymerized by PEO, and each PEO contains a hydrophobic C2H4 group and a hydrophilic O group. Therefore, the hydrophobic part accounts for a large proportion. In addition, some hydrophobic tails are distributed outside of the micelle. These hydrophobic groups will also increase the proportion of the hydrophobic surface area of the micelle.

2.2. Micelle Structure

The spatial distribution of atoms in a micelle can characterize the structural properties of the micelle [38]. The main difference between the non-ionic surfactant molecules of the TX family lies in the number of ethylene oxide units or the chain length [29]. In order to analyze the structures of the three micelles TX-5, TX-114, and TX-100 studied in this paper, we calculated the probability distributions of O atoms at different positions in each micelle system and the carbon atoms at the end of the hydrophobic tail chain relative to the micelle center of mass (COM). All the results are shown in Figure 2.
In general, the nonpolar alkyl chain of amphiphilic surfactant molecules usually aggregates in the hydrophobic region of the micelle, while the polar chain is exposed in the hydrophilic region. Sodium dodecyl sulfate (SDS) is a commonly used ionic surfactant and has been widely studied. Shelley et al. explored the structure of the SDS micelle in aqueous solution by MD simulation. It was found that most of the terminal C atoms of the SDS hydrophobic chain were distributed in the interior of the micelle, and only a small part were distributed on the surface of the micelle [39]. Gao et al. [40], Palazzesi et al. [38], MacKerell et al. [41], and Liu et al. [35] also found a similar phenomenon. It should be noted that these studies are ionic surfactants with small tail chains and polar heads. For TX micelle, there is a relatively large benzene ring near the end of the hydrophobic tail chain and long polar PEO chain head, which is very different from the SDS ionic surfactant. Pacheco-Blas et al. [27] used MD simulation to study the behavior of the non-ionic surfactant Triton X-100 in extracting metal ion Cd2+ in an aqueous environment. In the paper, it was found that the hydrophobic chain had a relatively high distribution in the outside of the micelle. In addition, Zhang et al. [7] considered that the TX-100 micelle was more likely to be a multilayer spherical micelle. Molecules from different layers were staggered, and octyl phenyl was dispersed in different layers. In this study, we found that the probability distribution of the tail C of the TX-100 micelle with respect to micelle COM is relatively wide, ranging from the micelle core to the outside of the micelle. In particular, the probability distribution of the terminal C in the outer region of the micelle is higher. The TX-5 and TX-114 micelles also show a similar phenomenon. This may be due to the large benzene rings in the hydrophobic region of the micelles, which are spatially exclusive and widely distributed. Although the tail C has a high probability in the outer region of the micelle, its distribution is also located on the left side of the distribution of the outermost polar O atoms. This means that most tail C atoms are included by the external polar O atoms, to avoid contact with environmental water as much as possible. The distributions of both hydrophobic terminal C and polar head O atoms are multimodal, indicating the existence of the micelle multilayer structure.
In order to show the differences among the three micelle structures intuitively, we highlighted the C atom at the end of each micelle and the outermost O atom on the polar head to observe their internal structures more clearly. It can be seen from Figure 3 that the micelle shape is consistent with that obtained by the ratio of the moment of inertia and the eccentricity. In agreement with the multimodal probability distributions of terminal C atoms from the micelle center to the micelle surface, some terminal C atoms of the hydrophobic chains are distributed widely in the whole micelle.

2.3. Interactions between Micelle and Water

2.3.1. Hydration Numbers

The interactions between the micelle and water are characterized by calculating the hydration numbers of C and O atoms at different positions in TX-5, TX-114, and TX-100 surfactant molecule. For convenient analysis, the selected C and O atoms are grouped. Among them, a C atom with a similar chemical environment was selected and averaged for calculation, as shown in Figure 4. The hydration numbers were obtained by calculating the radial distribution function (RDF) integral of water molecules around the selected C atom or O atom in the range of 0.35 nm [40,42]. It can be seen from Figure 4 that the three systems have a similar distribution. From the end of the hydrophobic chain to the position of the polar head, the hydration numbers of the selected atoms show increasing trends. The hydration numbers of C atoms 1–4 on the hydrophobic chain are very low. The hydration numbers of C atoms in the benzene ring increase slightly. For C atoms and O atoms in the PEO chain, the hydration numbers present a zigzag trend. The hydration number of the O atom in each ethylene oxide (EO) unit is greater than that of the C atom. Moreover, the hydration numbers of the O atom closer to the position of the polar head is larger. The hydration numbers of the outermost O atoms of the polar head are significantly larger than those of other atoms (>3) with the decreased sequence of TX-100 > TX-114 ≈ TX-5. The hydrations numbers of H atoms on the OH groups also show this sequence. In addition, the hydration numbers of the C atoms at the same position in the three micelles are similar.

2.3.2. Number of Hydrogen Bonds

In order to further characterize the interactions between the three micelles and water, we calculated the average numbers of hydrogen bonds between the O atoms at different positions and the water molecules. The geometric criterion for the existence of hydrogen bonds is that the distance between donor and acceptor pairs is within 3.5 Å and the angle between hydroxyl and hydrogen atoms is less than 120° [43]. As shown in Figure 5, the numbers of hydrogen bonds of the three systems show a significant upward trend from inner O atom to micelle surface O atom. Among them, the numbers of hydrogen bonds formed between the O atom close to the hydrophobic moiety and the water molecules are the lowest. In previous discussions, the hydration numbers of the outmost polar head O atoms of the three systems were the largest, ca. 4. However, the number of hydrogen bonds formed between these O atoms and water molecules is ca. 0.3, which means that not all the hydrated water molecules can form hydrogen bonds.

2.3.3. Time Correlation Function of the Hydrogen Bond

One method to analyze the dynamic properties of the hydrogen bond is to calculate the intermittent hydrogen bond time correlation function (CHB(t)), which is defined as follows:
C HB ( t ) = 〈 h ( t ) h ( 0 ) 〉 ( h )
This function is considered to be the criterion for judging the formation or fracture of hydrogen bonds, and CHB(t) allows the hydrogen bonds to break and reform within the time interval t [44]. The faster CHB(t) decays, the more unstable is the hydrogen bond [45].
In this paper, the CHB(t) changes of the hydrogen bonds between O atoms at different positions and water molecules and between the outermost hydroxyl group and O atoms of water molecules in the three systems were calculated by this function. As shown in Figure 6, TX-5, TX-114, and TX-100 systems all show three fast decay curves for O atoms connected to hydrophobic chain, the outermost O and H atoms of the polar head. Among these three curves, the outermost H atoms decay fastest. The hydration numbers of the outermost H atoms are the largest and they have a spatial advantage, so that the water molecules that form hydrogen bonds with them exchange faster. This may result in the faster decay curve of the hydrogen bond. The outmost O atoms also behave the same. For internal O atoms, they may behave due to the environmental impact of being connected to the benzene ring of the hydrophobic chain and the spatial limits.

3. Materials and Methods

In this paper, we chose the aggregation number of TX-5, TX-114, and TX-100 micelles as 100 according to the experimental values [46,47]. The molecular structures of TX-5, TX-114, and TX-100 were optimized by Gaussian 09 software [48] at the 6-31g(d) level. Then, the force field parameters of all studied molecules were obtained through the automated force field topology builder (ATB) [49,50]. With Packmol software [51], 100 pre-aggregated surfactant molecules for the three micelles were constructed. Then, the TX-5, TX-114, and TX-100 systems were treated with Gromacs 2019 software [52]. Pre-constructed micelle was placed in the center of a 10 × 10 × 10 nm3 cube box for each system by inbuilt tools. With the gmx solvate tool, water molecules were randomly added to fill the box. This inbuilt tool can ensure that water molecules are added in the box at locations larger than the volume of water molecules. Subsequently, the tension of the system can be released through pre-equilibrium processes. The simulated systems can achieve equilibrium and obtain a stable conformation. The single point charge (SPC) model was selected to describe water molecules. Specific compositions of all the studied systems are shown in Table 2. After minimizing the energy of the system, a reasonable initial structure was obtained for each system. With the V-rescale method [53], the 2 ns NVT equilibrium was performed to reach the required 298 K. The constant τT of 0.1 ps was used. Then, the NPT equilibrium with 2 ns was carried out to stabilize the system at 1 bar with a Berendsen pressure controller [54]. In this step, a constant τp of 2.0 ps was applied. After balancing the temperature and pressure of the system, the 200 ns NVT simulations were carried out to obtain the MD production trajectories. The same temperature controlling method as the NVT equilibrium step was used in the MD production step. In the process of simulation, periodic boundary conditions were used in all directions of the simulation system. The LINCS algorithm [55] was used to limit hydrogen-included bonds, and the particle-mesh Ewald (PME) [56] summation method was used for electrostatic interactions. The 1.4 nm cutoff value was used for van der Waals interactions and short-range electrostatic interactions. The integration time step was 2 fs. For every 10 ps interval, the MD sampling configurations were stored.

4. Conclusions

Through MD simulations, we studied the structural characteristics of TX-5, TX-114, and TX-100 micelles at the same aggregation number and for the related properties of each micelle in aqueous solution. With the increase of the PEO chain length, the degree of micelle sphericity was reduced while the micelle radius increased. It is worth noting that the tail C atoms of each micelle were distributed widely with multimodal distributions extending from the micelle COM to the micelle surface. Large parts of the hydrophobic tail of the C atoms were located in the outer region of the micelle. Some polar O atoms close to the hydrophobic chains also showed multimodal distributions. These distributions indicate the existence of a multilayer structure of the TX micelles. The distributions of the tail chains and the polar heads of these three micelles are significantly different from those of traditional cationic and anionic surfactant micelles. For TX-5, TX-114, and TX-100 micelles, the outermost O and H atoms have the largest hydration numbers and hydrogen bond numbers, but the decay of the hydrogen bonds is the fastest. The outermost OH groups of TX-100 show stronger interactions with the surrounding water molecules than those of TX-5 and TX-114. More importantly, the representative structures of TX-5, TX-114, and TX-100 micelles can serve as the structural basis for understanding the properties of the relevant micelles. This study provides a reference for further expanding the exploration and application of TX-5, TX-114, and TX-100 micelles in related fields and it is of great significance for the exploration of other non-ionic surfactant micelle systems.

Author Contributions

Conceptualization, J.P., G.L., Y.W. and Q.X.; methodology, J.P. and X.S.; software, J.P. and X.L.; validation, G.L., Y.W. and Q.X.; formal analysis, J.P., X.S., X.L. and Y.J.; investigation, J.P., X.S., X.L., Y.J., G.L., Y.W. and Q.X.; data curation, G.L., Y.W. and Q.X.; writing—original draft preparation, J.P. and X.S.; writing—review and editing, G.L., Y.W. and Q.X.; supervision, G.L., Y.W. and Q.X.; funding acquisition, G.L. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Shandong Province, grant number ZR2021QB153 and ZR2022QB043.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented will be made available on request by the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Barzykin, A.V.; Tachiya, M. Interpretation of passive permeability measurements on lipid-bilayer vesicles. Effect of fluctuations. BBA Biomembr. 1997, 1330, 121–126. [Google Scholar] [CrossRef] [PubMed]
  2. Azum, N.; Rub, M.A.; Azim, Y.; Asiri, A.M. Micellar and spectroscopic studies of amphiphilic drug with nonionic surfactant in the presence of ionic liquids. J. Mol. Liq. 2020, 315, 113732. [Google Scholar] [CrossRef]
  3. Jaiswal, S.; Mondal, R.; Paul, D.; Mukherjee, S. Investigating the micellization of the triton-X surfactants: A non-invasive fluorometric and calorimetric approach. Chem. Phys. Lett. 2016, 646, 18–24. [Google Scholar] [CrossRef]
  4. Abu-Ghunmi, L.; Badawi, M.; Fayyad, M. Fate of Triton X-100 Applications on Water and Soil Environments: A Review. J. Surfactants. Deterg. 2014, 17, 833–838. [Google Scholar] [CrossRef]
  5. Molina-Bolívar, J.A.; Aguiar, J.; Ruiz, C.C. Growth and Hydration Of Triton X-100 Micelles In Monovalent Alkali Salts: A Light Scattering Study. J. Phys. Chem. B 2002, 106, 870–877. [Google Scholar] [CrossRef]
  6. Parra, J.G.; Iza, P.; Dominguez, H.; Schott, E.; Zarate, X. Effect of Triton X-100 surfactant on the interfacial activity of ionic surfactants SDS, CTAB and SDBS at the air/water interface: A study using molecular dynamic simulations. Colloid. Surface. A 2020, 603, 125284. [Google Scholar] [CrossRef]
  7. Zhang, L.; Chai, X.; Sun, P.; Yuan, B.; Jiang, B.; Zhang, X.; Liu, M. The Study of the Aggregated Pattern of TX100 Micelle by Using Solvent Paramagnetic Relaxation Enhancements. Molecules 2019, 24, 1649. [Google Scholar] [CrossRef]
  8. Brown, W.; Rymden, R.; Van Stam, J.; Almgren, M.; Svensk, G. Static and dynamic properties of nonionic amphiphile micelles: Triton X-100 in aqueous solution. J. Phys. Chem. 1989, 93, 2512–2519. [Google Scholar] [CrossRef]
  9. Li, M.; Rharbi, Y.; Huang, X.; Winnik, M.A. Small variations in the composition and properties of Triton X-100. J. Colloid Interface Sci. 2000, 230, 135–139. [Google Scholar] [CrossRef]
  10. Thakkar, K.; Patel, V.; Ray, D.; Pal, H.; Aswal, V.K.; Bahadur, P. Interaction of imidazolium based ionic liquids with Triton X-100 micelles: Investigating the role of the counter ion and chain length. Rsc. Adv. 2016, 6, 36314–36326. [Google Scholar] [CrossRef]
  11. Denkova, P.S.; Lokeren, L.V.; Verbruggen, I.; Willem, R. Self-Aggregation and Supramolecular Structure Investigations of Triton X-100 and SDP2S by NOESY and Diffusion Ordered NMR Spectroscopy. J. Phys. Chem. B 2008, 112, 10935–10941. [Google Scholar] [CrossRef] [PubMed]
  12. Yuan, H.-Z.; Cheng, G.-Z.; Zhao, S.; Miao, X.-J.; Yu, J.-Y.; Shen, L.-F.; Du, Y.-R. Conformational Dependence of Triton X-100 on Environment Studied by 2D NOESY and 1H NMR Relaxation. Langmuir 2000, 16, 3030–3035. [Google Scholar] [CrossRef]
  13. Paradies, H.H. Shape and size of a nonionic surfactant micelle. Triton X-100 in aqueous solution. J. Phys. Chem. 1980, 84, 599–607. [Google Scholar] [CrossRef]
  14. Kumbhakar, M.; Goel, T.; Mukherjee, T.; Pal, H. Role of Micellar Size and Hydration on Solvation Dynamics: A Temperature Dependent Study in Triton-X-100 and Brij-35 Micelles. J. Phys. Chem. B 2004, 108, 19246–19254. [Google Scholar] [CrossRef]
  15. Ruiz, C.C.; Molina-Bolívar, J.A.; Aguiar, J.; MacIsaac, G.; Moroze, S.; Palepu, R. Thermodynamic and Structural Studies of Triton X-100 Micelles in Ethylene Glycol-Water Mixed Solvents. Langmuir 2001, 17, 6831–6840. [Google Scholar] [CrossRef]
  16. Robson, R.J.; Dennis, E.A. The size, shape, and hydration of nonionic surfactant micelles. Triton X-100. J. Phys. Chem. 1977, 81, 1075–1078. [Google Scholar] [CrossRef]
  17. Kumbhakar, M.; Nath, S.; Mukherjee, T.; Pal, H. Solvation dynamics in triton-X-100 and triton-X-165 micelles: Effect of micellar size and hydration. J. Chem. Phys. 2004, 121, 6026–6033. [Google Scholar] [CrossRef]
  18. Fedyaeva, O.A.; Poshelyuzhnaya, E.G. Dimensions and Orientation of Triton X-100 Micelles in Aqueous Solutions, According to Turbidimetric Data. Russ. J. Phys. Chem. A 2019, 93, 2559–2561. [Google Scholar] [CrossRef]
  19. Streletzky, K.; Phillies, G.D. Temperature dependence of Triton X-100 micelle size and hydration. Langmuir 1995, 11, 42–47. [Google Scholar] [CrossRef]
  20. Mandal, A.; Ray, S.; Biswas, A.; Moulik, S. Physicochemical studies on the characterization of Triton X 100 micelles in an aqueous environment and in the presence of additives. J. Phys. Chem. 1980, 84, 856–859. [Google Scholar] [CrossRef]
  21. Phillies, G.D.; Yambert, J.E. Solvent and solute effects on hydration and aggregation numbers of Triton X-100 micelles. Langmuir 1996, 12, 3431–3436. [Google Scholar] [CrossRef]
  22. Phillies, G.D.; Stott, J.; Ren, S. Probe diffusion in the presence of nonionic amphiphiles: Triton X 100. J. Phys. Chem. 1993, 97, 11563–11568. [Google Scholar] [CrossRef]
  23. Leach, A.R. Molecular Modelling: Principles and Applications; Pearson Education: London, UK, 2001. [Google Scholar]
  24. De Nicola, A.; Kawakatsu, T.; Rosano, C.; Celino, M.; Rocco, M.; Milano, G. Self-Assembly of Triton X-100 in Water Solutions: A Multiscale Simulation Study Linking Mesoscale to Atomistic Models. J. Chem. Theory. Comput. 2015, 11, 4959–4971. [Google Scholar] [CrossRef] [PubMed]
  25. Murakami, W.; De Nicola, A.; Oya, Y.; Takimoto, J.-I.; Celino, M.; Kawakatsu, T.; Milano, G. Theoretical and Computational Study of the Sphere-to-Rod Transition of Triton X-100 Micellar Nanoscale Aggregates in Aqueous Solution: Implications for Membrane Protein Purification and Membrane Solubilization. ACS Appl. Nano Mater. 2021, 4, 4552–4561. [Google Scholar] [CrossRef]
  26. Farafonov, V.; Lebed, A.; Khimenko, N.; Mchedlov-Petrossyan, N. Molecular Dynamics Study of an Acid-Base Indicator Dye in Triton X-100 Non-Ionic Micelles. Vopr. Khimii I Khimicheskoi Tekhnologii 2020, 1, 97–103. [Google Scholar] [CrossRef]
  27. Pacheco-Blas, M.D.A.; Vicente, L. Molecular dynamics study of the behaviour of surfactant Triton X-100 in the extraction process of Cd2+. Chem. Phys. Lett. 2020, 739, 136920. [Google Scholar] [CrossRef]
  28. Ritter, E.; Yordanova, D.; Gerlach, T.; Smirnova, I.; Jakobtorweihen, S. Molecular dynamics simulations of various micelles to predict micelle water partition equilibria with COSMOmic: Influence of micelle size and structure. Fluid Phase Equilib. 2016, 422, 43–55. [Google Scholar] [CrossRef]
  29. Yordanova, D.; Smirnova, I.; Jakobtorweihen, S. Molecular Modeling of Triton X Micelles: Force Field Parameters, Self-Assembly, and Partition Equilibria. J. Chem. Theory Comput. 2015, 11, 2329–2340. [Google Scholar] [CrossRef]
  30. Rauf, M.A.; Hisaindee, S.; Graham, J.P.; Al-Zamly, A. Effect of various solvents on the absorption spectra of dithizone and DFT calculations. J. Mol. Liq. 2015, 211, 332–337. [Google Scholar] [CrossRef]
  31. Faramarzi, S.; Bonnett, B.; Scaggs, C.A.; Hoffmaster, A.; Grodi, D.; Harvey, E.; Mertz, B. Molecular Dynamics Simulations as a Tool for Accurate Determination of Surfactant Micelle Properties. Langmuir 2017, 33, 9934–9943. [Google Scholar] [CrossRef]
  32. Wei, Y.; Wang, H.; Liu, G.; Wang, Z.; Yuan, S. A molecular dynamics study on two promising green surfactant micelles of choline dodecyl sulfate and laurate. RSC Adv. 2016, 6, 84090–84097. [Google Scholar] [CrossRef]
  33. Allen, D.T.; Saaka, Y.; Lawrence, M.J.; Lorenz, C.D. Atomistic description of the solubilisation of testosterone propionate in a sodium dodecyl sulfate micelle. J. Phys. Chem. B 2014, 118, 13192–13201. [Google Scholar] [CrossRef] [PubMed]
  34. Wei, Y.; Liu, G.; Wang, Z.; Yuan, S. Molecular dynamics study on the aggregation behaviour of different positional isomers of sodium dodecyl benzenesulphonate. RSC Adv. 2016, 6, 49708–49716. [Google Scholar] [CrossRef]
  35. Liu, G.; Zhang, H.; Liu, G.; Yuan, S.; Liu, C. Tetraalkylammonium interactions with dodecyl sulfate micelles: A molecular dynamics study. Phys. Chem. Chem. Phys. 2016, 18, 878–885. [Google Scholar] [CrossRef] [PubMed]
  36. Eisenhaber, F.; Lijnzaad, P.; Argos, P.; Sander, C.; Scharf, M. The double cubic lattice method: Efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies. J. Comput. Chem. 1995, 16, 273–284. [Google Scholar] [CrossRef]
  37. Tang, X.; Koenig, P.H.; Larson, R.G. Molecular dynamics simulations of sodium dodecyl sulfate micelles in water-the effect of the force field. J. Phys. Chem. B 2014, 118, 3864–3880. [Google Scholar] [CrossRef]
  38. Palazzesi, F.; Calvaresi, M.; Zerbetto, F. A molecular dynamics investigation of structure and dynamics of SDS and SDBS micelles. Soft Matter. 2011, 7, 9148–9156. [Google Scholar] [CrossRef]
  39. Shelley, J.; Watanabe, K.; Klein, M.L. Simulation of a sodium dodecylsulfate micelle in aqueous solution. Int. J. Quantum Chem. 1990, 38, 103–117. [Google Scholar] [CrossRef]
  40. Gao, J.; Ge, W.; Hu, G.; Li, J. From Homogeneous Dispersion to Micelles A Molecular Dynamics Simulation on the Compromise of the Hydrophilic and Hydrophobic Effects of Sodium Dodecyl Sulfate in Aqueous Solution. Langmuir 2005, 21, 5223–5229. [Google Scholar] [CrossRef]
  41. MacKerell, A.D. Molecular Dynamics Simulation Analysis of a Sodium Dodecyl Sulfate Micelle in Aqueous Solution: Decreased Fluidity of the Micelle Hydrocarbon Interior. J. Phys. Chem. 1995, 99, 1846–1855. [Google Scholar] [CrossRef]
  42. Chun, B.J.; Choi, J.I.; Jang, S.S. Molecular dynamics simulation study of sodium dodecyl sulfate micelle: Water penetration and sodium dodecyl sulfate dissociation. Colloid. Surface. A 2015, 474, 36–43. [Google Scholar] [CrossRef]
  43. Liu, G.; Li, R.; Wei, Y.; Gao, F.; Wang, H.; Yuan, S.; Liu, C. Molecular dynamics simulations on tetraalkylammonium interactions with dodecyl sulfate micelles at the air/water interface. J. Mol. Liq. 2016, 222, 1085–1090. [Google Scholar] [CrossRef]
  44. Pal, S.; Bagchi, B.; Balasubramanian, S. Hydration Layer of a Cationic Micelle, C10TAB:  Structure, Rigidity, Slow Reorientation, Hydrogen Bond Lifetime, and Solvation Dynamics. J. Phys. Chem. B 2005, 109, 12879–12890. [Google Scholar] [CrossRef]
  45. Song, X.; Zhang, X.; Peng, J.; Li, Y.; Leng, X.; Liu, G.; **a, Q.; Wei, Y. Molecular dynamics study on the aggregation behaviours of Platonic micelle in different NaCl solutions. J. Mol. Liq. 2022, 353, 118828. [Google Scholar] [CrossRef]
  46. Pal, A.; Pillania, A. Modulations in surface and aggregation properties of non-ionic surfactant Triton X-45 on addition of ionic liquids in aqueous media. J. Mol. Liq. 2017, 233, 243–250. [Google Scholar] [CrossRef]
  47. Qiao, L.; Easteal, A.J. Mass transport in Triton X series nonionic surfactant solutions: A new approach to solute-solvent interactions. Colloid Polym. Sci. 1996, 274, 974–980. [Google Scholar] [CrossRef]
  48. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision, B.01; Gaussian, Inc.: Wallingford, CT, USA, 2010. [Google Scholar]
  49. Koziara, K.B.; Stroet, M.; Malde, A.K.; Mark, A.E. Testing and validation of the Automated Topology Builder (ATB) version 2.0: Prediction of hydration free enthalpies. J. Comput. Aided Mol. Des. 2014, 28, 221–233. [Google Scholar] [CrossRef] [PubMed]
  50. Malde, A.K.; Zuo, L.; Breeze, M.; Stroet, M.; Poger, D.; Nair, P.C.; Oostenbrink, C.; Mark, A.E. An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0. J. Chem. Theory Comput. 2011, 7, 4026–4037. [Google Scholar] [CrossRef]
  51. Martínez, L.; Andrade, R.; Birgin, E.G.; Martínez, J.M. PACKMOL: A package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157–2164. [Google Scholar] [CrossRef]
  52. Lindahl, E.; Abraham, M.; Hess, B.; van der Spoel, D. GROMACS 2019.3 Source Code (2019.3). Zenodo. 2019. Available online: https://doi.org/10.5281/zenodo.3243833 (accessed on 10 March 2023).
  53. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef]
  54. Berendsen, H.J.C.; Postma, J.P.M.V.; Van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef]
  55. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  56. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef]
Molecules 28 03557 g001 550
Figure 1. (a) The radius of gyration of TX-5, TX-114, and TX-100 micelles with time evolution, (b) the solvent accessible surface area of TX-5, TX-114, and TX-100 micelles with time evolution, and (c) the total energy of TX-5, TX-114, and TX-100 micelles with time evolution.
Figure 1. (a) The radius of gyration of TX-5, TX-114, and TX-100 micelles with time evolution, (b) the solvent accessible surface area of TX-5, TX-114, and TX-100 micelles with time evolution, and (c) the total energy of TX-5, TX-114, and TX-100 micelles with time evolution.
Molecules 28 03557 g001
Molecules 28 03557 g002 550
Figure 2. Probability distribution of O atoms at different position and the carbon at the end of the hydrophobic tail chain in TX-5, TX-114, and TX-100 micelle systems with respect to micelle COM.
Figure 2. Probability distribution of O atoms at different position and the carbon at the end of the hydrophobic tail chain in TX-5, TX-114, and TX-100 micelle systems with respect to micelle COM.
Molecules 28 03557 g002
Molecules 28 03557 g003 550
Figure 3. Representative structures of TX-5, TX-114 and TX-100 micelles. Note: the dark green ball in the figure symbolizes C atom, the red ball symbolizes O atom, and the blue auxiliary line symbolizes the outline.
Figure 3. Representative structures of TX-5, TX-114 and TX-100 micelles. Note: the dark green ball in the figure symbolizes C atom, the red ball symbolizes O atom, and the blue auxiliary line symbolizes the outline.
Molecules 28 03557 g003
Molecules 28 03557 g004 550
Figure 4. Hydration numbers of C and O atoms at different positions in TX-5, TX-114, and TX-100 surfactant molecules. Note: the Triton X molecular diagram on the right describes the atoms corresponding to each number. Dark green sphere represents C atoms, red sphere represents O atoms, and gray-white sphere represents H atoms.
Figure 4. Hydration numbers of C and O atoms at different positions in TX-5, TX-114, and TX-100 surfactant molecules. Note: the Triton X molecular diagram on the right describes the atoms corresponding to each number. Dark green sphere represents C atoms, red sphere represents O atoms, and gray-white sphere represents H atoms.
Molecules 28 03557 g004
Molecules 28 03557 g005 550
Figure 5. The average numbers of hydrogen bonds formed between O atoms at different positions and water molecules in a surfactant molecule in TX-5, TX-114, and TX-100 micelles. The O atom label is same as the schematic diagram of Figure 2.
Figure 5. The average numbers of hydrogen bonds formed between O atoms at different positions and water molecules in a surfactant molecule in TX-5, TX-114, and TX-100 micelles. The O atom label is same as the schematic diagram of Figure 2.
Molecules 28 03557 g005
Molecules 28 03557 g006 550
Figure 6. Time correlation function between O atoms at different positions and water molecules, and between the outermost hydroxyl group and O atoms in water molecules in TX-5, TX-114, and TX-100 micelle systems. The O atom label is the same as the schematic diagram of Figure 2. Note: OW represents the O atom in a water molecule.
Figure 6. Time correlation function between O atoms at different positions and water molecules, and between the outermost hydroxyl group and O atoms in water molecules in TX-5, TX-114, and TX-100 micelle systems. The O atom label is the same as the schematic diagram of Figure 2. Note: OW represents the O atom in a water molecule.
Molecules 28 03557 g006
Table
Table 1. Structural characteristics of micelles Triton X-5, Triton X-114, Triton X-100.
Table 1. Structural characteristics of micelles Triton X-5, Triton X-114, Triton X-100.
TX-5TX-114TX-100
Imax/Imin1.15 ± 0.031.12 ± 0.041.26 ± 0.03
e0.07 ± 0.020.06 ± 0.020.12 ± 0.02
Rg (nm)1.96 ± 0.012.07 ± 0.012.19 ± 0.01
Rs (nm)2.53 ± 0.012.67 ± 0.012.83 ± 0.01
SASA(total) (nm2)132.24 ± 2.77147.52 ± 3.25159.54 ± 2.97
SASA(hydrophilic) (nm2)30.62 ± 2.3332.72 ± 2.2735.81 ± 1.96
SASA(hydrophobic) (nm2)101.66 ± 2.19114.77 ± 2.29123.71 ± 2.54
Table
Table 2. Specific composition of micelle system Triton X-5, Triton X-114, Triton X-100.
Table 2. Specific composition of micelle system Triton X-5, Triton X-114, Triton X-100.
MoleculeTX-5 (n = 5)TX-114 (n = 7)TX-100 (n = 9)
Nagg100100100
H2O30,83030,26429,751
n represents oxyethylene groups.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peng, J.; Song, X.; Li, X.; Jiang, Y.; Liu, G.; Wei, Y.; **a, Q. Molecular Dynamics Study on the Aggregation Behavior of Triton X Micelles with Different PEO Chain Lengths in Aqueous Solution. Molecules 2023, 28, 3557. https://doi.org/10.3390/molecules28083557

AMA Style

Peng J, Song X, Li X, Jiang Y, Liu G, Wei Y, **a Q. Molecular Dynamics Study on the Aggregation Behavior of Triton X Micelles with Different PEO Chain Lengths in Aqueous Solution. Molecules. 2023; 28(8):3557. https://doi.org/10.3390/molecules28083557

Chicago/Turabian Style

Peng, **, **aoju Song, **n Li, Yongkang Jiang, Guokui Liu, Yaoyao Wei, and Qiying **a. 2023. "Molecular Dynamics Study on the Aggregation Behavior of Triton X Micelles with Different PEO Chain Lengths in Aqueous Solution" Molecules 28, no. 8: 3557. https://doi.org/10.3390/molecules28083557

Article Metrics

Back to TopTop