Molecular dynamics simulations of structure-property relationships of Tween 80 surfactants in water and at interfaces

Coarse-Grained Molecular Dynamics Simulation of Heterogeneous Polysorbate 80 Surfactants and their Interactions with Small Molecules and Proteins

Polysorbate 80 (PS80) is widely used in pharmaceutical formulations, and its commercial grades exhibit certain levels of structural heterogeneity. The objective of this study was to apply coarse-grained molecular dynamics simulations to better understand the effect of PS80 heterogeneity on micelle self-assembly, the loading of hydrophobic small molecules into the micelle core, and the interactions between PS80 and a protein, bovine serum albumin (BSA). Four representative PS80 variants with different head and tail structures were studied. Our simulations found that PS80 structural heterogeneity could affect blank micelle properties such as solvent-accessible surface area, aggregation number, and micelle aspect ratio. It was also found that hydrophobic small molecules such as ethinyl estradiol preferentially partitioned into the PS80 micelle core and PS80 dioleates formed a more hydrophobic core compared to PS80 monooleates. Furthermore, multiple PS80 molecules could bind to BSA, and PS80 heterogeneity profoundly changed the binding ratio as well as the surfactant-protein contact area.

Development and Evaluation of Docetaxel-Loaded Nanostructured Lipid Carriers for Skin Cancer Therapy

This study focuses on the design, characterization, and optimization of nanostructured lipid carriers (NLCs) loaded with docetaxel for the treatment of skin cancer. Employing a systematic formulation development process guided by Design of Experiments (DoE) principles, key parameters such as particle size, polydispersity index (PDI), zeta potential, and entrapment efficiency were optimized to ensure the stability and drug-loading efficacy of the NLCs. Combined XRD and cryo-TEM analysis were employed for NLC nanostructure evaluation, confirming the formation of well-defined nanostructures. In vitro kinetics studies demonstrated controlled and sustained docetaxel release over 48 h, emphasizing the potential for prolonged therapeutic effects. Cytotoxicity assays on human umbilical vein endothelial cells (HUVEC) and SK-MEL-24 melanoma cell line revealed enhanced efficacy against cancer cells, with significant selective cytotoxicity and minimal impact on normal cells. This multidimensional approach, encompassing formulation optimization and comprehensive characterization, positions the docetaxel-loaded NLCs as promising candidates for advanced skin cancer therapy. The findings underscore the potential translational impact of these nanocarriers, paving the way for future preclinical investigations and clinical applications in skin cancer treatment.

Catalytic degradation of carbamazepine by surface-modified zerovalent copper via the activation of peroxymonosulphate: mechanism, degradation pathways and ecotoxicity

In this research work, surface-modified nano zerovalent copper (nZVC) was prepared using a simple borohydride reduction method. The spectroscopic and crystallographic results revealed the successful synthesis of surface-modified nano zerovalent copper (nZVC) using solvents such as ethanol (ETOH), ethylene glycol (EG) and tween80 (T80). The as-synthesized material was fully characterized for morphological surface and crystal structural properties. The results indicated that EG provides an excellent synthesis environment to nZVC compared to ETOH and T80 in terms of good dispersion, high surface area and excellent catalytic properties. The catalytic efficiency of nZVC/EG was investigated alone and with peroxymonosulphate (PMS) in the absence of light. The degradation results demonstrated that the involvement of PMS synergistically boosted the catalytic efficiency of synthesized nZVC/EG material. Furthermore, the degradation products (DPs) of CBZ were determined by GC-MS and subsequently, the degradation pathways were proposed. The ecotoxicity analysis of the DPs was also explored. The proposed (nZVC/EG/PMS) system is economical and efficient and thus could be applied for the degradation of CBZ from an aquatic system after altering the degradation pathways in such a way that results in harmless products.

Linseed Oil-Based Oleogel Vehicles for Hydrophobic Drug Delivery-Physicochemical and Applicative Properties

In this study, a methodology for synthesizing oleogels based on linseed oil and emulsifiers, such as beeswax and Tween 20 and Tween 80, was developed. Linseed oil served as the main oil phase, while beeswax acted as a gelling and emulsifying agent. Tween compounds are non-ionic surfactants composed of hydrophobic and hydrophilic parts, allowing for the formation of a stable system with promising properties. Surface wetting analysis of the obtained oleogels, FT-IR spectroscopy, and determination of relative and absolute humidity over time, as well as optical microscope analysis and rheological analysis of the obtained oleogels, were conducted as part of the research. The results indicate that increasing the amount of Tween 20 decreases the hydrophilicity of the oleogel, while Tween 80 exhibits the opposite effect. Surface energy analysis suggests that a higher content of Tween 20 may lead to a reduction in the surface energy of the oleogels, which may indicate greater material stability. Changes in relative humidity and FT-IR spectral analysis confirm the influence of emulsifiers on the presence of characteristic functional groups in the structure of the oleogels. Additionally, microscopic analysis suggests that an emulsifier with a longer hydrophobic tail leads to a denser material structure.

Computational predictions of interfacial tension, surface tension, and surfactant adsorption isotherms

All-atom (AA) molecular dynamics (MD) simulations are employed to predict interfacial tensions (IFT) and surface tensions (ST) of both ionic and non-ionic surfactants. The general AMBER force field (GAFF) and variants are examined in terms of their performance in predicting accurate IFT/ST, γ, values for chosen water models, together with the hydration free energy, ΔGhyd, and density, ρ, predictions for organic bulk phases. A strong correlation is observed between the quality of ρ and γ predictions. Based on the results, the GAFF-LIPID force field, which provides improved ρ predictions is selected for simulating surfactant tail groups. Good γ predictions are obtained with GAFF/GAFF-LIPID parameters and the TIP3P water model for IFT simulations at a water-triolein interface, and for GAFF/GAFF-LIPID parameters together with the OPC4 water model for ST simulations at a water-vacuum interface. Using a combined molecular dynamics-molecular thermodynamics theory (MD-MTT) framework, a mole fraction of C12E6 molecule of 1.477 × 10-6 (from the experimental critical micelle concentration, CMC) gives a simulated surface excess concentration, ΓMAX, of 76 C12E6 molecules at a 36 nm2 water-vacuum surface (3.5 × 10-10 mol cm-2), which corresponds to a simulated ST of 35 mN m-1. The results compare favourably with an experimental ΓMAX of C12E6 of 3.7 × 10-10 mol cm-2 (80 surfactants for a 36 nm2 surface) and experimental ST of C12E6 of 32 mN m-1 at the CMC.

Kinetic description of water transport during spontaneous emulsification induced by Span 80

Spontaneous emulsification is a phenomenon that forms nanometer-sized droplets (nanodroplets) without the application of any external force, and the mechanism has been actively studied for application to various technologies. In this study, we analyzed the kinetics of spontaneous emulsification induced by Span 80. The measurement of water concentration in Span 80 hexadecane solution indicated that the chemical potential of water in the nanodroplets decreased as the amount of water in the nanodroplets decreased. Based on this result, water transport between the aqueous phase and nanodroplets in which the chemical potential of water was controlled was quantitatively investigated by using a microfluidic device. The results demonstrate that the kinetics of water transport during spontaneous emulsification induced by Span 80 was described by a model of osmotic transport through an organic liquid film between the aqueous phase and nanodroplets.

Ecological impact of surfactant Tween-80 on plankton: High-scale analyses reveal deeper hazards

The ecological risks of surfactants have been largely neglected because of their low toxicity. Multiscale studies have indicated that even if a pollutant causes no acute toxicity in a test species, it may alter interspecific interactions and community characteristics through sublethal impacts on test organisms. Therefore, we investigated the lethal and sublethal responses of the plankton species Scenedesmus quadricauda, Chlorella vulgaris, and Daphnia magna, to surfactant Tween-80. Then, high-scale responses in grazer life-history traits and stability of the D. magna-larval damselfly system were further explored. The results showed that discernible adverse effects on the growth or survival of the three plankton species were evident only at exceptionally high concentrations (≥100 mg L-1). However, 10 mg L-1 of Tween-80 notably affected the MDA concentration in grazer species, simultaneously displaying a tendency to diminish grazer's heartbeat and swimming frequency. Furthermore, Tween-80 reduced the grazer reproductive capacity and increased its predation risk by larval damselflies, which ultimately jeopardized the stability of the D. magna-larval damselfly system at much lower concentrations (10-100 fold lower) than the individual-scale responses. This study provides evidence that high-scale traits are far more sensitive to Tween-80, compared with individual-scale traits for plankton organisms, suggesting that the ecological risks of Tween-80 demand careful reassessment. SYNOPSIS: The concentration of Tween-80 needed to induce changes in community characteristics is markedly lower than that needed to produce individual-scale consequences. Thus, high-scale analyses have broad implications for understanding the hazardous effects of surfactants compared with an individual-scale analysis.

HRMAS-NMR and simulation study of the self-assembly of surfactants on carbon nanotubes

Polyethoxylated surfactants, such as those of the Tween and Pluronic series, are commonly used to disperse carbon nanotubes (CNTs) and other nanoparticles. However, the current understanding of the nature of interactions between these surfactants and CNTs is limited. The nature of the interactions between surfactants (Tween-80 [T80] and Pluronic F68 [PF68]) and CNTs was investigated using high-resolution magic angle spinning nuclear magnetic resonance (HRMAS-NMR) and coarse-grained molecular dynamics (MD) simulations. HRMAS-NMR revealed that T80 molecules interact with single-walled CNTs (SWCNTs) and multi-walled CNTs (MWCNTs) via the oleyl chain, whereas PF68 molecules interact with the surface of SWCNTs and MWCNTs via the polypropylene oxide residues. The polyethylene oxide chains were oriented towards the external aqueous environment. The HRMAS-NMR results were supported by MD simulations, and the latter provided further insights into the nature of the interactions.

Tween-80 on Water/Oil Interface: Structure and Interfacial Tension by Molecular Dynamics Simulations

Surfactants are used in many fields of the chemical industry in a wide range of applications. Generally, surfactants on water/oil interfaces reduce interfacial tension, enabling the formation of emulsions or providing greater stability to the emulsion formed. Although molecular dynamics has been extensively used and has achieved remarkable success in describing thermodynamic and molecular properties of systems with surface-active compounds, the traditional molecular simulation force fields considerably constrain the system size and the time scale of simulations. Here, we propose a coarse-grained model of polysorbate 80, a nonionic surfactant commercially known as Tween-80. Based on the proposed coarse-grained model, we evaluate the influence of the more internal ethylene oxide chain on the properties of the water/Tween-80/decane system. We verify with the simulation results that the model can reproduce the expected decrease in interfacial tension as the surfactant quantity increases in the simulations. Furthermore, we observe changes in the surfactant orientation as their quantity in the interface increases, indicating a preferential orientation for these molecules in the adsorption layer. We also assessed the partition coefficient of the surfactant between the two bulk phases by performing free energy calculations, which showed a higher affinity of Tween-80 surfactants with the water phase.

Impact of Surfactants on the Performance of Clopidogrel-Copovidone Amorphous Solid Dispersions: Increased Drug Loading and Stabilization of Nanodroplets

PURPOSE

Surfactants are increasingly being added to amorphous solid dispersion (ASDs) formulations to enhance processability and release performance. The goal of the current work was to investigate the impact of cationic, anionic and non-ionic surfactants on the rate and extent of clopidogrel (CPD) release from copovidone-based ASDs.

METHODS

CPD release was evaluated for ASDs with different drug loadings using a surface normalized intrinsic dissolution apparatus. Studies were also carried out using dynamic light scattering, zeta potential measurements, and nuclear magnetic resonance spectroscopy to probe the impact of surfactants on drug-rich nanodroplet physical stability and clopidogrel-surfactant interactions.

RESULTS

CPD ASDs showed good release for drug loadings as high as 40%, before the release fell off a cliff at higher drug loadings. Only sodium dodecyl sulfate, added at a 5% level, was able to improve the release at 50% drug loading, with other surfactants proving to be ineffective. However, some of the surfactants evaluated did show some benefits in improving nanodroplet stability against size enlargement. Ionic and non-ionic surfactants were observed to interact differently with CPD-rich nanodroplets, and variations in the kinetics and morphology of water-induced phase separation were noted in the presence and absence of surfactants in ASD films.

CONCLUSIONS

In summary, addition of surfactants to ASD formulations may lead to some improvements in formulation performance, but predictive capabilities and mechanisms of surfactant effect still require further studies.

Computational and Experimental Models of Type III Lipid-Based Formulations of Loratadine Containing Complex Nonionic Surfactants

Type III lipid-based formulations (LBFs) combine poorly water-soluble drugs with oils, surfactants, and cosolvents to deliver the drugs into the systemic circulation. However, the solubility of the drug can be influenced by the colloidal phases formed in the gastrointestinal tract as the formulation is dispersed and makes contact with bile and other materials present within the GI tract. Thus, an understanding of the phase behavior of LBFs in the gut is critical for designing efficient LBFs. Molecular dynamics (MD) simulation is a powerful tool for the study of colloidal systems. In this study, we modeled the internal structures of five type III LBFs of loratadine containing poly(ethylene oxide) nonionic surfactants polysorbate 80 and polyoxyl hydrogenated castor oil (Kolliphor RH40) using long-timescale MD simulations (0.4-1.7 μs). We also conducted experimental investigations (dilution of formulations with water) including commercial Claritin liquid softgel capsules. The simulations show that LBFs form continuous phase, water-swollen reverse micelles, and bicontinuous and phase-separated systems at different dilutions, which correlate with the experimental observations. This study supports the use of MD simulation as a predictive tool to determine the fate of LBFs composed of medium-chain lipids, polyethylene oxide surfactants, and polymers.

Characterization of the Partition Rate of 2-Naphthol (NAP) and Ritonavir (RTV) Across the Water-Octanol Interface and the Influence of Common Pharmaceutical Excipients

The partition rate of 2 model compounds, 2-naphthol (NAP) and ritonavir (RTV), across the water-octanol (W/O) interface was determined from aqueous solutions with and without the presence of a variety of common pharmaceutical excipients. The NAP molecule was present as either a unionized species in the phosphate buffer solution or as an anion in 0.01 M NaOH solution while RTV was present as dications in 0.1 M HCl. Excipients examined include different type of polymers (e.g., PVP-VA, HPMC, HPMCAS, PVP, PEG 8000), small molecules (e.g., glucose, lactose, maltoheptaose), and surfactants (e.g., Tween 80 and SDS). A noticeable to significant reduction of the partition rate of both NAP and RTV across the W/O interface was observed with aqueous media containing 0.01-0.1 μg/mL polymers including PVP-VA, HPMC, and HPMCAS. Reduction of NAP and RTV partition rate associated with such extremely low concentration of polymers in the aqueous media affirms the presence of a surface excess of adsorbed excipients at the W/O interface. Pi values of NAP and RTV across the W/O interface are found to be sensitive to (1) the molecule surface area, molecular weight and the degree of the polymerization of these pharmaceutical excipients, and (2) the ionization state of the model compound.

Improvement of Foaming Ability of Surfactant Solutions by Water-Soluble Polymers: Experiment and Molecular Dynamics Simulation

Aqueous foam is widely used in fire extinguishing and dust suppression technologies. Improving the foaming ability is the key to reducing the added concentration of foaming agents as well as the economic cost. In this work, the effect of a water-soluble polymer (polyvinyl alcohol, PVA) on the foaming ability of anionic surfactant (sodium dodecyl ether sulfate, SDES) was studied by an experiment and molecular dynamics simulation. The experimental results showed that PVA greatly improves the foaming ability of SDES solutions when the surfactant concentration is less than 0.1%, which is attributed to the fact that the polymer can enhance the stability of bubble films and reduce the bubble rupture rate during the foam generation process. The simulation results indicate that PVA can enhance the hydration of surfactant head groups and contribute to the formation of a three-dimensional hydrogen bond network between surfactants, polymers, and water molecules, thus greatly improving the stability of bubble liquid films. The above results suggest that water-soluble polymers can be used to improve the foaming ability of surfactant solutions by enhancing the bubble film stability, which is beneficial as it reduces the added concentration of foaming agents in aqueous foam applications.

Deep Learning to Speed up the Development of Structure-Property Relations For Hexagonal Boron Nitride and Graphene

Structure-property maps play a key role in accelerated materials discovery. The current norm for developing these maps includes computationally expensive physics-based simulations. Here, the capabilities of deep learning agents are explored such as convolutional neural networks (CNNs) and multilayer perceptrons (MLPs) to predict structure-property relations and reduce dependence on simulations. This study contains simulated hexagonal boron nitride (h-BN) microstructures damaged by various levels of radiation and temperature, with the objective to predict their residual strengths from the final atomic positions. By developing low dimensional physical descriptors to statistically describe the defects, these results show that purpose-specific microstructure representation can help in achieving a good prediction accuracy at low computational cost. Furthermore, the adaptability of the trained deep learning agents is explored to predict structure-property maps of other 2D materials using transfer learning. It is shown that in order to achieve good predictions accuracy (≈95% R² ), an agent that is training for the first time ("learning from scratch") requires 23-45% of simulated data, whereas an agent adapting to a different material ("transfer learning") requires only about 10% or less. This suggests that transfer learning is a potential game changer in material discovery and characterization approaches.

Oil-in-Water fL Droplets by Interfacial Spontaneous Fragmentation and Their Electrical Characterization

Inkjet printing is here employed for the first time as a method to produce femtoliter-scale oil droplets dispersed in water. In particular, picoliter-scale fluorinated oil (FC40) droplets are printed in the presence of perfluoro-1-octanol surfactant at a velocity higher than 5 m/s. Femtoliter-scale oil droplets in water are spontaneously formed through a fragmentation process at the water/air interface using minute amounts of nonionic surfactant (down to 0.003% v/v of Tween 80). This fragmentation occurs by a Plateau-Rayleigh mechanism at a moderately high Weber number (10¹). A microfluidic chip with integrated microelectrodes allows droplets characterization in terms of number and diameter distribution (peaked at about 3 μm) by means of electrical impedance measurements. These results show an unprecedented possibility to scale oil droplets down to the femtoliter scale, which opens up several perspectives for a tailored oil-in-water emulsion fabrication for drug encapsulation, pharmaceutic preparations, and cellular biology.

Adsorption and Aggregation Properties of Some Polysorbates at Different Temperatures

Measurements of the surface tension of aqueous solutions of polysorbates (Tween 20, Tween 60 and Tween 80) at 293, 303 and 313 K were made. On the basis of the obtained results the Gibbs surface excess concentration of the Tweens at the water–air interface and critical micelle concentrations were determined. Knowing the Gibbs surface excess concentration and taking into account the difference between the limiting area occupied by water and Tween molecules at the water–air interface, the fraction occupied by Tween molecules was established. The limiting area occupied by the Tween molecule was calculated by applying the Joos equation. The area determined in such a way was confirmed by the calculations of cross section of Tween molecules based on the bond lengths and the angles between them as well as the average distance between the molecules, taking into account their different conformations. This area was used for calculation of the standard Gibbs energy of adsorption using the Langmuir equation. The standard Gibbs energy of Tweens adsorption at the water–air interface was also calculated from the hydrophobic part of Tween molecule–water interface tension and that of hydrophobic part. Using the determined values of standard Gibbs energy of adsorption at different temperatures, the standard enthalpy and entropy values were deduced. The standard thermodynamic functions of micellization were also determined and compared to the Gibbs energy of Tween molecules interactions through the water phase.

Molecular-Scale Description of SPAN80 Desorption from a Squalane-Water Interface

Extensive all-atom molecular dynamics calculations on the water-squalane interface for nine different loadings with sorbitan monooleate (SPAN80), at T = 300 K, are analyzed for the surface tension equation of state, desorption free-energy profiles as they depend on loading, and to evaluate escape times for adsorbed SPAN80 into the bulk phases. These results suggest that loading only weakly affects accommodation of a SPAN80 molecule by this squalane-water interface. Specifically, the surface tension equation of state is simple through the range of high tension to high loading studied, and the desorption free-energy profiles are weakly dependent on loading here. The perpendicular motion of the centroid of the SPAN80 headgroup ring is well-described by a diffusional model near the minimum of the desorption free-energy profile. Lateral diffusional motion is weakly dependent on loading. Escape times evaluated on the basis of a diffusional model and the desorption free energies are 7 × 10^-2 s (into the squalane) and 3 × 10² h (into the water). The latter value is consistent with desorption times of related lab-scale experimental work.

Reverse micellar extraction of lactoferrin from its synthetic solution using CTAB/n-heptanol system

The partitioning of Lactoferrin (LF) into the reverse micellar phase formed by a cationic surfactant, cetyltrimethylammonium bromide (CTAB) in n-heptanol from the synthetic solution of LF was studied. The solubilization behaviour of LF into the reverse micellar phase and back extraction using a fresh stripping phase were improved by studying the effect of processing parameters, including surfactant concentration, solution pH, electrolyte salt concentration and addition of alcohol as co-solvent. Forward extraction of 100% was achieved at CTAB concentration of 50 mM in n-heptanol solvent, pH of 10 and 1 M NaCl. The electrostatic force and hydrophobic interaction have major influence on LF extraction during forward and back extraction respectively. The size of the reverse micelles and their corresponding water content were measured at different operating conditions to assess their role on the LF extraction. The present reverse micellar system has potential to solubilise almost all the LF into the reverse micelles during forward extraction and could able to allow back extraction from the reverse micellar phase with addition of small amount of co-solvent.

Combined Molecular Dynamics Simulation-Molecular-Thermodynamic Theory Framework for Predicting Surface Tensions

A molecular modeling approach is presented with a focus on quantitative predictions of the surface tension of aqueous surfactant solutions. The approach combines classical Molecular Dynamics (MD) simulations with a molecular-thermodynamic theory (MTT) [ Y. J. Nikas, S. Puvvada, D. Blankschtein, Langmuir 1992 , 8 , 2680 ]. The MD component is used to calculate thermodynamic and molecular parameters that are needed in the MTT model to determine the surface tension isotherm. The MD/MTT approach provides the important link between the surfactant bulk concentration, the experimental control parameter, and the surfactant surface concentration, the MD control parameter. We demonstrate the capability of the MD/MTT modeling approach on nonionic alkyl polyethylene glycol surfactants at the air-water interface and observe reasonable agreement of the predicted surface tensions and the experimental surface tension data over a wide range of surfactant concentrations below the critical micelle concentration. Our modeling approach can be extended to ionic surfactants and their mixtures with both ionic and nonionic surfactants at liquid-liquid interfaces.

Micelles, Rods, Liposomes, and Other Supramolecular Surfactant Aggregates: Computational Approaches

Surfactants are an interesting class of compounds characterized by the segregation of polar and apolar domains in the same molecule. This peculiarity makes possible a whole series of microscopic and macroscopic effects. Among their features, their ability to segregate particles (fluids or entire domains) and to reduce the surface/interfacial tension is the utmost important. The interest in the chemistry of surfactants never weakened; instead, waves of increasing interest have occurred every time a new field of application of these molecules has been discovered. All these special characteristics depend largely on the ability of surfactants to self-assemble and constitute supramolecular structures where their chemical properties are amplified. The possibility to obtain structural and energy information and, above all, the possibility of forecast the self-organizing mechanisms of surfactants have had a significant boost via computational chemistry. The molecular dynamics models, initially coarse-grained and subsequently (with the increasing computer power) using more accurate models, allowed, over the years, to better understand different aspects of the processes of dispersion, self-assembly, segregation of surfactant. Moreover, several other aspects have been investigated as the effect of the counterions of many ionic surfactants in defining the final supramolecular structures, the mobility of side chains, and the capacity of some surfactant to envelope entire proteins. This review constitutes a perspective/prospective view of these results. On the other hand, some comparison of in silico results with experimental information recently acquired through innovative analytical techniques such as ion mobility mass spectrometry which have been introduced.

Insoluble layer deposition and dilatational rheology at a microscale spherical cap interface

The dilatational properties of insoluble monolayers are important for understanding the mechanics of biological systems and consumer products, but isolating the dilatational response of an interface is challenging due to the difficulties in separating dilatation from shear and other deformation modes. Oscillations of a microscale bubble radius are useful for generating purely dilatational flows, but the current deposition methods for insoluble layers onto fluid interfaces are not easily scaled down. In this paper, we describe a miscible solvent exchange procedure for generating insoluble layers at an air-water interface pinned at the tip of a capillary tens of micrometers in diameter. We show that the amount of surfactant adsorbed at the interface can be controlled by the initial concentration dissolved in isopropanol (the starting solvent) and the volumetric flow rate of solvent exchange. Surface pressure-area isotherms and dilatational moduli are measured concurrently for three insoluble surfactants: palmitic acid (PA), dimyristoylphosphatidylcholine (DMPC) and dipalmitoylphosphatidylcholine (DPPC). The isotherms measured on the microscale interface compare well with previous experiments performed on a Langmuir trough. However, the magnitudes of the dilatational moduli differ from those measured on either Langmuir trough or pendant drop apparatuses. Several possible reasons for the observed differences are discussed. A comparison of the dilatational modulus with the Gibbs elasticity is used to determine the presence of dilatational extra stresses at the interface. The isotherm and dilatational modulus of the insoluble component of the industrial surfactant Tween 80 are measured using this approach. The methods developed here also open the possibility for future study of the important role of finite size effects on microstructure formation and the resulting interfacial mechanics.

Interfacial tension of oil/water emulsions with mixed non-ionic surfactants: comparison between experiments and molecular simulations

Smaller Span molecules occupy the free spaces between bulkier Tween molecules thus lowering interfacial tension as compared to those obtained for single surfactant systems.

Solubilization Behavior of Polyene Antibiotics in Nanomicellar System: Insights from Molecular Dynamics Simulation of the Amphotericin B and Nystatin Interactions with Polysorbate 80

Amphotericin B (AmB) and Nystatin (Nys) are the drugs of choice for treatment of systemic and superficial mycotic infections, respectively, with their full clinical potential unrealized due to the lack of high therapeutic index formulations for their solubilized delivery. In the present study, using a coarse-grained (CG) molecular dynamics (MD) simulation approach, we investigated the interaction of AmB and Nys with Polysorbate 80 (P80) to gain insight into the behavior of these polyene antibiotics (PAs) in nanomicellar solution and derive potential implications for their formulation development. While the encapsulation process was predominantly governed by hydrophobic forces, the dynamics, hydration, localization, orientation, and solvation of PAs in the micelle were largely controlled by hydrophilic interactions. Simulation results rationalized the experimentally observed capability of P80 in solubilizing PAs by indicating (i) the dominant kinetics of drugs encapsulation over self-association; (ii) significantly lower hydration of the drugs at encapsulated state compared with aggregated state; (iii) monomeric solubilization of the drugs; (iv) contribution of drug-micelle interactions to the solubilization; (v) suppressed diffusivity of the encapsulated drugs; (vi) high loading capacity of the micelle; and (vii) the structural robustness of the micelle against drug loading. Supported from the experimental data, our simulations determined the preferred location of PAs to be the core-shell interface at the relatively shallow depth of 75% of micelle radius. Deeper penetration of PAs was impeded by the synergistic effects of (i) limited diffusion of water; and (ii) perpendicular orientation of these drug molecules with respect to the micelle radius. PAs were solvated almost exclusively in the aqueous poly-oxyethylene (POE) medium due to the distance-related lack of interaction with the core, explaining the documented insensitivity of Nys solubilization to drug-core compatibility in detergent micelles. Based on the obtained results, the dearth of water at interior sites of micelle and the large lateral occupation space of PAs lead to shallow insertion, broad radial distribution, and lack of core interactions of the amphiphilic drugs. Hence, controlled promotion of micelle permeability and optimization of chain crowding in palisade layer may help to achieve more efficient solubilization of the PAs.

Multiscale Molecular Dynamics Simulations of Model Hydrophobically Modified Ethylene Oxide Urethane Micelles

The flower-like micelles of various aggregation numbers of a model hydrophobically modified ethylene oxide urethane (HEUR) molecule, C16E45C16, and their corresponding starlike micelles, containing the surfactants C16E22 and C16E23, were studied by atomistic and coarse-grained molecular dynamic (MD) simulations. We used free energies from umbrella sampling to calculate the size distribution of micelle sizes and the average time for escape of a hydrophobic group from the micelle. Using the coarse-grained MARTINI force field, the most probable size of the model HEUR molecule was thereby determined to be about 80 hydrophobes per micelle and the average hydrophobe escape time to be about 0.1 s, both of which are consistent with previous experimental studies. Atomistic simulations reveal that hydrogen bond formation and the mean lifetime of hydration waters of the poly(ethylene oxide) (or PEO) groups are location-dependent in the HEUR micelle, with PEO groups immediately adjacent to the C16 groups forming the fewest hydrogen bonds with water and having hydration waters with longer lifetimes than those of the PEO groups located further away from the C16 groups.

Reversible and Irreversible Adsorption Energetics of Poly(ethylene glycol) and Sorbitan Poly(ethoxylate) at a Water/Alkane Interface

We simulate poly(ethylene glycol) (PEG) oligomers and model Tween 80 (polyoxyethylene sorbitan monooleate) molecules at water/alkane interfaces. Using the weighted histogram analysis method (WHAM), including an extension of WHAM to two reaction coordinates to remove hysteresis, we calculate interfacial potentials of mean force (PMFs) for PEG and Tween 80 using three force fields: the atomistic GROMOS 53a6OXY+D and two coarse-grained (CG) MARTINI force fields. Because the force fields have not yet been validated for PEO adsorption to hydrophobic interfaces, we calculate PMFs for alcohol ethoxylates C12E2 and C12E8 and find that they agree with semiempirical results of Mulqueen and Blankschtein [Langmuir 2002, 18 (2), 365-376] for the GROMOS 53a6OXY+D force field, whereas for both MARTINI force fields, PEO adsorbs too weakly to a clean hydrophobic interface. One MARTINI force field incorrectly shows depletion rather than adsorption to a clean hydrophobic interface. We find that the adsorption free energy for PEG oligomers at a clean, planar water/alkane interface is around 1.3 kBT per monomer for the atomistic force field but is less than half of this for the two CG force fields. With the newly validated GROMOS 53a6OXY+D force field, we bracket the dilute adsorption free energy for a model Tween 80 molecule at the clean water/squalane interface. We also calculate the pressure-area isotherm. We exploit these data with the Nikas-Mulqueen-Blankschtein (NMB) theory and a simple transport model to demonstrate a transition from irreversible to reversible adsorption with increasing surface coverage, consistent with experimental results of Reichert and Walker [Langmuir 2013, 29 (6), 1857-1867].

Sequential adsorption of an irreversibly adsorbed nonionic surfactant and an anionic surfactant at an oil/aqueous interface

Aerosol-OT (AOT) and Tween 80 are two of the main surfactants in commercial dispersants used in response to oil spills. Understanding how multicomponent surfactant systems interact at oil/aqueous interfaces is crucial for improving both dispersant design and application efficacy. This is true of many multicomponent formulations; a lack of understanding of competition for the oil/water interface hinders formulation optimization. In this study, we have characterized the sequential adsorption behavior of AOT on squalane/aqueous interfaces that have been precoated with Tween 80. A microtensiometer is used to measure the dynamic interfacial tension of the system. Tween 80 either partially or completely irreversibly adsorbs to squalane/aqueous interfaces when rinsed with deionized water. These Tween 80 coated interfaces are then exposed to AOT. AOT adsorption increases with AOT concentration for all Tween 80 coverages, and the resulting steady-state interfacial tension values are interpreted using a Langmuir isotherm model. In the presence of 0.5 M NaCl, AOT adsorption significantly increases due to counterion charge screening of the negatively charged head groups. The presence of Tween 80 on the interface inhibits AOT adsorption, reducing the maximum surface coverage as compared to a clean interface. Tween 80 persists on the interface even after exposure to high concentrations of AOT.

Molecular dynamics investigation on adsorption layer of alcohols at the air/brine interface

Alcohols are a significant group of surfactants which have been employed extensively in industry to improve the interfacial effects. Recently, the change in surface potential (ΔV) of two isomeric hexanols, methyl isobutyl carbinol (MIBC) and 1-hexanol, was investigated by using an ionizing (241)Am electrode. It clearly showed the opposite effects between MIBC and 1-hexanol in the interfacial zone: one enhanced the presence of cations, whereas the other enhanced the presence of anions. This study employs molecular dynamics simulation to provide new insights into the interactions between alcohol molecules and ions as well as water at the molecular level. The results qualitatively agreed with the experimental data and verified the significance of MIBC branching structure on the molecular arrangement within the interfacial zone. The results also highlighted the role of the second water layer on the interfacial properties.