Computer Simulations of Phosphatidylcholine Monolayers at Air/Water and CCl4/Water Interfaces

Modeling Gas-Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air-Water Interface

Adsorption of surfactants at gas-liquid interfaces that causes reduction in the surface tension is a classical problem in colloid and interface science with multiple practical applications in oil and gas recovery, separations, cosmetics, personal care, and biomedicine. Here, we develop an original coarse-grained model of the liquid-gas interface within the conventional dissipative particle dynamics (DPD) framework with the goal of quantitatively predicting the surface tension in the presence of surfactants. As a practical case-study example, we explore the adsorption of the cationic surfactant cetyl trimethyl ammonium bromide (CTAB) on the air-water interface. The gas phase is modeled as a DPD fluid composed of fictitious hard-core "gas" beads with exponentially decaying repulsive potentials to prevent penetration of the liquid phase components. A rigorous parametrization scheme is proposed based on matching the bulk and interfacial properties of water and octane taken as the reference compounds. Quantitative agreement between the simulated and experimental surface tension of CTAB solutions is found for a wide range of bulk surfactant concentrations (∼10^-3 to ∼1 mmol/L) with the reduction of the surface tension from ∼72 mN/m (pure water) to the limiting value of ∼37.5 mN/m at the critical micelle concentration. The gas phase DPD model with the proposed parametrization scheme can be extended and applied to modeling various gas-liquid interfaces with surfactant and lipid monolayers, such as bubble suspensions, foams, froths, etc.

Size-Dependent Interactions of Lipid-Coated Gold Nanoparticles: Developing a Better Mechanistic Understanding Through Model Cell Membranes and in vivo Toxicity

Introduction

Humans are intentionally exposed to gold nanoparticles (AuNPs) where they are used in variety of biomedical applications as imaging and drug delivery agents as well as diagnostic and therapeutic agents currently in clinic and in a variety of upcoming clinical trials. Consequently, it is critical that we gain a better understanding of how physiochemical properties such as size, shape, and surface chemistry drive cellular uptake and AuNP toxicity in vivo. Understanding and being able to manipulate these physiochemical properties will allow for the production of safer and more efficacious use of AuNPs in biomedical applications.

Methods and Materials

Here, AuNPs of three sizes, 5 nm, 10 nm, and 20 nm, were coated with a lipid bilayer composed of sodium oleate, hydrogenated phosphatidylcholine, and hexanethiol. To understand how the physical features of AuNPs influence uptake through cellular membranes, sum frequency generation (SFG) was utilized to assess the interactions of the AuNPs with a biomimetic lipid monolayer composed of a deuterated phospholipid 1.2-dipalmitoyl-d62-sn-glycero-3-phosphocholine (dDPPC).

Results and Discussion

SFG measurements showed that 5 nm and 10 nm AuNPs are able to phase into the lipid monolayer with very little energetic cost, whereas, the 20 nm AuNPs warped the membrane conforming it to the curvature of hybrid lipid-coated AuNPs. Toxicity of the AuNPs were assessed in vivo to determine how AuNP curvature and uptake influence cell health. In contrast, in vivo toxicity tested in embryonic zebrafish showed rapid toxicity of the 5 nm AuNPs, with significant 24 hpf mortality occurring at concentrations ≥20 mg/L, whereas the 10 nm and 20 nm AuNPs showed no significant mortality throughout the five-day experiment.

Conclusion

By combining information from membrane models using SFG spectroscopy with in vivo toxicity studies, a better mechanistic understanding of how nanoparticles (NPs) interact with membranes is developed to understand how the physiochemical features of AuNPs drive nanoparticle-membrane interactions, cellular uptake, and toxicity.

Molecular Ordering in Lipid Monolayers: An Atomistic Simulation

We report on atomistic simulations of DPPC lipid monolayers using the CHARMM36 lipid force field (and also the Slipid force field as a control case), combined with a four-point OPC water model. The entire two-phase region where domains of the "liquid-condensed" (LC) phase coexist with domains of the "liquid-expanded" (LE) phase has been explored. The simulations are long enough that the complete phase-transition stage, with two domains coexisting in the monolayer, is reached in all cases. Also, system sizes used are larger than those in previous works. As expected, domains of the minority phase are elongated, emphasizing the importance of anisotropic van der Waals and/or electrostatic dipolar interactions in the monolayer plane. The molecular structure is quantified in terms of distribution functions for the hydrocarbon chains and the PN dipoles. In contrast to previous work, where average distributions are calculated, distributions are here extracted for each of the coexisting phases by first identifying lipid molecules that belong to either LC or LE regions. In the case of the CHARMM36 force field, the three-dimensional distributions show that the average tilt angle of the chains with respect to the normal outward direction is (39.0 ± 0.1)° in the LC phase and (48.1 ± 0.5)° in the LC phase. In the case of the PN dipoles, the distributions indicate a tilt angle of (110.8 ± 0.5)° in the LC phase and (112.5 ± 0.5)° in the LE phase. These results are quantitatively different from those in previous works, which indicated a smaller normal component of the PN dipole. Also, the distributions of the monolayer-projected chains and PN dipoles have been calculated. Chain distributions peak along a particular direction in the LC domains, while they are uniform in the LE phase. Long-range ordering associated with the projected PN dipoles is absent in both phases. These results strongly suggest that LC domains do not exhibit dipolar ordering in the plane of the monolayer, the effect of these components being averaged out at short distances. Therefore, the only relevant component of the molecular dipoles, with regard to both intra- and long-range interdomain interactions, is normal to the monolayer. Also, the local orientation of chain projections is almost constant in LC domains and points in the direction along which domains are elongated, suggesting that the line tension driving the phase transition might be anisotropic with respect to the interfacial domain boundary.

Anisotropic line tension of domains in lipid monolayers

We formulate a simple effective model to describe molecular interactions in a lipid monolayer and calculate the line tension between coexisting domains. The model represents lipid molecules in terms of two-dimensional anisotropic particles on the plane of the monolayer. These particles interact through forces that are believed to be relevant for the understanding of fundamental properties of the monolayer: van der Waals interactions originating from lipid chains and dipolar forces between dipole groups in the molecular heads. The model stresses the liquid-crystalline nature of the ordered phase in lipid monolayers and explains coexistence properties between ordered and disordered phases in terms of molecular parameters. Thermodynamic and interfacial properties of the model are analyzed using density-functional theory. In particular, the line tension at the interface between ordered and disordered phases turns out to be highly anisotropic with respect to the angle between the nematic director and the interface separating the coexisting phases. This important feature mainly results from the tilt angle of lipid chains and, to a lesser extent, from dipolar interactions perpendicular to the monolayer. The role of the two dipolar components, parallel and perpendicular to the monolayer, is assessed by comparing with computer simulation results for lipid monolayers.

Characterisation of N-(Octadecyl)-1,8-naphthalimide Monolayer Compression Using Molecular Dynamics and Experimental Approaches

The development of luminescent surfaces is an active area of supramolecular chemistry, particularly for the development of new sensing platforms. One particularly useful surface deposition method is the Langmuir-Blodgett technique where organic amphiphilic fluorophores (e.g. 1,8-naphthalimides) can form ordered monolayers at an air-water interface before being deposited onto solid supports. The ability to simulate monolayer formation and consequently develop predictability over film formation would allow for significant advances in the luminescent materials field where synthesis might be directed by simulation data. Here, we compare pressure-area isotherms of N-(octadecyl)-1,8-naphthalimide determined experimentally, using the Langmuir-Blodgett technique, and computationally, using three different simulation techniques. We find that all three simulation techniques are able to describe the liquid-condensed/liquid-expanded region of the isotherm, and that the isotherms are highly similar in this region, although the NγT ensemble performs best. Experimental isotherms showed film formation properties that align with the simulation data, suggesting that simulations are a viable means to direct synthesis. Investigation of the underlying structural details disclosed by the simulations reveals the compression-induced ordering at atomic-level detail, which will allow prediction of how functionalisation of the naphthalimides will alter the monolayer compression and mounting process.

Zwitterionic and Charged Lipids Form Remarkably Different Structures on Nanoscale Oil Droplets in Aqueous Solution

The molecular structure of zwitterionic and charged monolayers on small oil droplets in aqueous solutions is determined using a combined second harmonic and sum frequency study. From the interfacial vibrational signature of the acyl chains and phosphate headgroups as well as the response of the hydrating water, we find that zwitterionic and charged lipids with identical acyl chains form remarkably different monolayers. Zwitterionic phospholipids form a closely packed monolayer with highly ordered acyl tails. In contrast, the charged phospholipids form a monolayer with a low number density and disordered acyl tails. The charged headgroups are oriented perpendicular to the monolayer rather than parallel, as is the case for zwitterionic lipids. These significant differences between the two types of phospholipids indicate important roles of phospholipid headgroups in the determination of properties of cellular membranes and lipid droplets. The observed behavior of charged phospholipids is different from expectations based on studies performed on extended planar interfaces, at which condensed monolayers are readily formed. The difference can be explained by nanoscale related changes in charge condensation behavior that has its origin in a different balance of interfacial intermolecular interactions.

Computer simulations of lung surfactant

Lung surfactant lines the gas-exchange interface in the lungs and reduces the surface tension, which is necessary for breathing. Lung surfactant consists mainly of lipids with a small amount of proteins and forms a monolayer at the air-water interface connected to bilayer reservoirs. Lung surfactant function involves transfer of material between the monolayer and bilayers during the breathing cycle. Lipids and proteins are organized laterally in the monolayer; selected species are possibly preferentially transferred to bilayers. The complex 3D structure of lung surfactant and the exact roles of lipid organization and proteins remain important goals for research. We review recent simulation studies on the properties of lipid monolayers, monolayers with phase coexistence, monolayer-bilayer transformations, lipid-protein interactions, and effects of nanoparticles on lung surfactant. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Three Dimensional Nano "Langmuir Trough" for Lipid Studies

A three-dimensional-phospholipid monolayer with tunable molecular structure was created on the surface of oil nanodroplets from a mixture of phospholipids, oil, and water. This simple nanoemulsion preparation technique generates an in situ prepared membrane model system with controllable molecular surface properties that resembles a lipid droplet. The molecular interfacial structure of such a nanoscopic system composed of hexadecane, 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), and water was determined using vibrational sum frequency scattering and second harmonic scattering techniques. The droplet surface structure of DPPC can be tuned from a tightly packed liquid condensed phase like monolayer to a more dilute one that resembles the liquid condensed/liquid expanded coexistence phase by varying the DPPC/oil/water ratio. The tunability of the chemical structure, the high surface-to-volume ratio, and the small sample volume make this system an ideal model membrane for biochemical research.

Coarse-graining poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO) block copolymers using the MARTINI force field

The MARTINI coarse-grain (CG) force field is extended for a class of triblock block copolymers known as Pluronics. Existing MARTINI bead types are used to model the non-bonded part of the potential while single chain properties for both homopolymers, poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO), are used to develop the bonded interactions. The new set of force field parameters reproduces structural and dynamical properties of high molecular weight homo- and copolymers. The CG model is moderately transferable in solvents of different polarity and concentration; however, the PEO homopolymer model presents a reduced thermodynamic transferability especially in water probably due to the lack of hydrogen bonds with the solvent. Our simulations of a monolayer of Pluronic L44 show polymer-brush-like characteristics for the PEO segments which protrude into the aqueous phase. Other membrane properties not easily accessible using experimental techniques such as its membrane thickness are also calculated.

Simulations of lipid monolayers

A lipid monolayer lining a boundary between two immiscible phases forms a complex interface with inhomogeneous distribution of forces. Unlike lipid bilayers, monolayers are formed in asymmetric environment and their properties depend strongly on lipid surface density. The monolayer properties are also affected significantly by the representation of the pure interface. Here we give a brief theoretical introduction and describe methods to simulate lipid monolayers starting from force-fields and system setup to reproducing state points on the surface tension (pressure)-area isotherms and transformations between them.

Molecular dynamics simulation of hydrated DPPC monolayers using charge equilibration force fields

We present results of molecular dynamics simulations of a model DPPC-water monolayer using charge equilibration (CHEQ) force fields, which explicitly account for electronic polarization in a classical treatment of intermolecular interactions. The surface pressure, determined as the difference between the monolayer and pure water surface tensions at 323 K, is predicted to be 22.92 ±1.29 dyne/cm, just slightly below the broad range of experimental values reported for this system. The surface tension for the DPPC-water monolayer is predicted to be 42.35 ±1.16 dyne/cm, in close agreement with the experimentally determined value of 40.9 dyne/cm. This surface tension is also consistent with the value obtained from DPPC monolayer simulations using state-of-the-art nonpolarizable force fields. The current results of simulations predict a monolayer-water potential difference relative to the pure water-air interface of 0.64 ±0.02 Volts, an improved prediction compared to the fixed-charge CHARMM27 force field, yet still overestimating the experimental range of 0.3 to 0.45 Volts. As the charge equilibration model is a purely charge-based model for polarization, the current results suggest that explicitly modeled polarization effects can offer improvements in describing interfacial electrostatics in such systems.

Variations in the condensing effect of cholesterol on saturated versus unsaturated phosphatidylcholines at low and high sterol concentration

In this work, we have investigated the condensing and ordering effect induced by cholesterol on phosphatidylcholines (PCs). To perform the studies systematically, for the experiments we have selected phospholipids differing only in the number of cis monounsaturated chains (1,2-distearoyl-sn-glycero-3-phosphocholine--DSPC, 1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine--SOPC, 1,2-dioleoyl-sn-glycero-3-phosphocholine--DOPC) or in the length (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine--POPC vs SOPC) of sn-1 acyl chain. Because the cholesterol concentration in mammalian membranes can be as high as 70 mol % of total lipids, the investigations were performed in a wide range of the sterol content. The results of the Langmuir monolayer experiments evidence that the relation between the structure of hydrophobic part of PC and the magnitude of the effects induced by cholesterol found at lower sterol content is different from that observed at higher sterol concentration. At a lower concentration of sterol (up to 30%), the condensing effect of cholesterol is stronger on saturated DSPC than on PCs containing monounsaturated chain(s), which is consistent with the conclusions drawn by other authors. However, at higher sterol content (≥50%), saturated DSPC is less susceptible to the influence of sterol than the investigated unsaturated PCs. To explain these irregularities, we have considered the strength of van der Waals interactions as well as the influence of sterol on the tilt of polar heads of PCs. It was also found that in the whole range of sterol concentration the ordering effect is stronger on saturated DSPC as compared to unsaturated phospholipids. However, at lower sterol content (up to 30%) the ordering effect induced on unsaturated PCs is rather weak, and the ordering does not change drastically in comparison with pure PCs film.

Structural inhomogeneity of interfacial water at lipid monolayers revealed by surface-specific vibrational pump-probe spectroscopy

We report vibrational lifetime measurements of the OH stretch vibration of interfacial water in contact with lipid monolayers, using time-resolved vibrational sum frequency (VSF) spectroscopy. The dynamics of water in contact with four different lipids are reported and are characterized by vibrational relaxation rates measured at 3200, 3300, 3400, and 3500 cm(-1). We observe that the water molecules with an OH frequency ranging from 3300 to 3500 cm(-1) all show vibrational relaxation with a time constant of T(1) = 180 ± 35 fs, similar to what is found for bulk water. Water molecules with OH groups near 3200 cm(-1) show distinctly faster relaxation dynamics, with T(1) < 80 fs. We successfully model the data by describing the interfacial water containing two distinct subensembles in which spectral diffusion is, respectively, rapid (3300-3500 cm(-1)) and absent (3200 cm(-1)). We discuss the potential biological implications of the presence of the strongly hydrogen-bonded, rapidly relaxing water molecules at 3200 cm(-1) that are decoupled from the bulk water system.

Interfacial water structure associated with phospholipid membranes studied by phase-sensitive vibrational sum frequency generation spectroscopy

Phase-sensitive vibrational sum frequency generation is employed to investigate the water structure at phospholipid/water interfaces. Interfacial water molecules are oriented preferentially by the electrostatic potential imposed by the phospholipids and have, on average, their dipole pointing toward the phospholipid tails for all phospholipids studied, dipalmitoyl phosphocholine (DPPC), dipalmitoyl phosphoethanolamine (DPPE), dipalmitoyl phosphate (DPPA), dipalmitoyl phosphoglycerol (DPPG), and dipalmitoyl phospho-l-serine (DPPS). Zwitterionic DPPC and DPPE reveal weaker water orienting capability relative to net negative DPPA, DPPG, and DPPS. Binding of calcium cations to the lipid phosphate group reduces ordering of the water molecules.

Observation of buried water molecules in phospholipid membranes by surface sum-frequency generation spectroscopy

We investigate the structure and orientation of water molecules at the water-lipid interface, using vibrational sum-frequency generation in conjunction with a maximum entropy phase retrieval method. We find that interfacial water molecules have an orientation opposite to that predicted by electrostatics and thus are likely localized between the lipid headgroup and its apolar alkyl chain. This type of water molecule is observed for phospholipids but not for structurally simpler surfactants.

Interactions of dimethylsulfoxide with a dipalmitoylphosphatidylcholine monolayer studied by vibrational sum frequency generation

The interactions between phospholipid monolayers and dimethylsulfoxide (DMSO) molecules were investigated by vibrational sum frequency generation (VSFG) spectroscopy in a Langmuir trough system. Both the head and the tail groups of dipalmitoylphosphatidylcholine (DPPC) as well as DMSO were probed to provide a comprehensive understanding of the interactions between DPPC and DMSO molecules. A condensing effect is observed for the DPPC monolayer on a concentrated DMSO subphase (>20 mol %). This effect results in a well-ordered conformation for the DPPC alkyl chains at very large mean molecular areas. Interactions between DMSO and DPPC headgroups were also studied. DMSO-induced dehydration of the DPPC phosphate group is revealed at DMSO concentration above 10 mol %. The average orientation of DMSO with DPPC versus dipalmitoylphosphate sodium salt (DPPA) monolayers was compared. The comparison revealed that DMSO molecules are perturbed and reorient because of the interfacial electric field created by the charged lipid headgroups. The orientation of the DPPC alkyl chains remains nearly unchanged in the liquid condensed phase with the addition of DMSO. This suggests that DMSO molecules are expelled from the condensed monolayer. In addition, implications for the DMSO-induced permeability enhancement of biological membranes from this work are discussed.

Molecular dynamics simulation study of a pulmonary surfactant film interacting with a carbonaceous nanoparticle

This article reports an all-atom molecular dynamics simulation to study a model pulmonary surfactant film interacting with a carbonaceous nanoparticle. The pulmonary surfactant is modeled as a dipalmitoylphosphatidylcholine monolayer with a peptide consisting of the first 25 residues from surfactant protein B. The nanoparticle model with a chemical formula C188H53 was generated using a computational code for combustion conditions. The nanoparticle has a carbon cage structure reminiscent of the buckyballs with open ends. A series of molecular-scale structural and dynamical properties of the surfactant film in the absence and presence of nanoparticle are analyzed, including radial distribution functions, mean-square displacements of lipids and nanoparticle, chain tilt angle, and the surfactant protein B peptide helix tilt angle. The results show that the nanoparticle affects the structure and packing of the lipids and peptide in the film, and it appears that the nanoparticle and peptide repel each other. The ability of the nanoparticle to translocate the surfactant film is one of the most important predictions of this study. The potential of mean force for dragging the particle through the film provides such information. The reported potential of mean force suggests that the nanoparticle can easily penetrate the monolayer but further translocation to the water phase is energetically prohibitive. The implication is that nanoparticles can interact with the lung surfactant, as supported by recent experimental data by Bakshi et al.

Comparing experimental and simulated pressure-area isotherms for DPPC

Although pressure-area isotherms are commonly measured for lipid monolayers, it is not always appreciated how much they can vary depending on experimental factors. Here, we compare experimental and simulated pressure-area isotherms for dipalmitoylphosphatidylcholine (DPPC) at temperatures ranging between 293.15 K and 323.15 K, and explore possible factors influencing the shape and position of the isotherms. Molecular dynamics simulations of DPPC monolayers using both coarse-grained (CG) and atomistic models yield results that are in rough agreement with some of the experimental isotherms, but with a steeper slope in the liquid-condensed region than seen experimentally and shifted to larger areas. The CG lipid model gives predictions that are very close to those of atomistic simulations, while greatly improving computational efficiency. There is much more variation among experimental isotherms than between isotherms obtained from CG simulations and from the most refined simulation available. Both atomistic and CG simulations yield liquid-condensed and liquid-expanded phase area compressibility moduli that are significantly larger than those typically measured experimentally, but compare well with some experimental values obtained under rapid compression.

Pressure-area isotherm of a lipid monolayer from molecular dynamics simulations

We calculated the pressure-area isotherm of a dipalmitoyl-phosphatidylcholine (DPPC) lipid monolayer from molecular dynamics simulations using a coarse-grained molecular model. We characterized the monolayer structure, geometry, and phases directly from the simulations and compared the calculated isotherm to experiments. The calculated isotherm shows liquid-expanded and liquid-condensed phases and their coexistence plateau. At high pressure, the monolayer surface is rippled; upon further compression, the monolayer undergoes a collapse. We studied the effect of temperature and system size on the isotherm slope and phase coexistence region. Thermodynamic and dynamic properties of the monolayer phases were also investigated.

Protein-directed assembly of binary monolayers at the interface and surface patterns of protein on the monolayers

Ferritin-directed assembly of binary monolayers of zwitterionic dipalmitoylphosphatidylcholine and cationic dioctadecyldimethylammonium bromide (DOMA) at the interface and surface patterns of ferritin on the monolayers have been investigated using a combination of infrared reflection absorption spectroscopy, surface plasmon resonance, and atomic force microscopy. Ferritin binding to the binary monolayers at the air-water interface at the surface pressure 30 mN/m, primarily driven by the electrostatic interaction, gives rise to a change in tilt angle of hydrocarbon chains from 15 degrees +/- 1 degrees to 10 degrees +/- 1 degrees with respect to the normal of the monolayer at the mole fraction of DOMA (XDOMA) of 0.1. The chains at XDOMA = 0.3 are oriented vertical to the water surface before and after protein binding. A new mechanism for protein binding to the binary monolayers is proposed. The secondary structures of the adsorbed ferritin are prevented from changing to some extent due to the existence of the monolayers. The amounts of the bound protein on the monolayers at the air-water interface are increased in comparison with those on the pre-immobilized monolayers at low XDOMA. The increased amounts and different patterns of the adsorbed protein at the monolayers are mostly attributed to the formation of multiple binding sites available for ferritin, which is due to the lateral reorganization of the lipid components in the monolayers induced by the protein in the subphase. The created multiple binding sites on the monolayer surfaces through the protein-directed assembly can be preserved for subsequent protein binding.

Adsorption of Octyl Cyanide at the Free Water Surface as Studied by Monte Carlo Simulation

Monte Carlo simulations of the adsorption layer of octyl cyanide have been performed on the canonical (N, V, T) ensemble at 300 K. The systems simulated cover the range of octyl cyanide surface densities from 0.27 to 7.83 mumol/m2. The surface density value at which the saturation of the adsorption layer occurs is estimated to be 1.7 mumol/m2. At low surface densities, the main driving force of the adsorption is found to be the formation of hydrogen bonds between the water and octyl cyanide molecules, whereas at higher surface concentrations, the dipole-dipole attraction between the neighboring adsorbed octyl cyanide molecules becomes more important. At low surface concentrations, the water-octyl cyanide hydrogen bonds prefer tilted alignments relative to the interface; however, in the case of the saturated adsorption layer, the number of such hydrogen bonds is maximized, leading to the preference of these bonds for the orientation perpendicular to the interface. Contrary to nonionic surfactants of multiple hydrogen bonding abilities (e.g., 1-octanol, C8E3), the increasing surface concentration of octyl cyanide was not found to lead to considerable competition of the molecules for positions of optimal arrangement. As a consequence, the energy and geometry of the water-octyl cyanide hydrogen bonds are found to be insensitive to the octyl cyanide surface concentration.

Counterion and Surface Density Dependence of the Adsorption Layer of Ionic Surfactants at the Vapor−Aqueous Solution Interface: A Computer Simulation Study

To test the validity of currently used adsorption theories and understand the origin of the lack of their ability of adequately describing existing surface tension measurement data, we have performed a series of molecular dynamics simulations of the adsorption layer of alkali decyl sulfate at the vapor/aqueous solution interface. The simulations have been performed with five different cations (i.e., Li+, Na+, K+, Rb+, and Cs+) at two different surface concentrations (i.e., 2 micromol/m2 and 4 micromol/m2). The obtained results clearly show that the thickness of the outer Helmholtz plate, a key quantity of the various adsorption theories, depends on two parameters, that is, the size of the cations and the surface density of the anionic surfactant. Namely, with increasing surface concentration, the electrostatic attraction between the two, oppositely charged, layers becomes stronger, leading to a considerable shrinking of the outer Helmholtz plate. Furthermore, this layer is found to be thicker in the presence of larger cations. The former effect could be important in understanding the anomalous shape of the adsorption isotherms of alkali alkyl sulfate surfactants, while the second effect seems to be essential in explaining the cation specificity of these isotherms.

Equation of state for a coarse-grained DPPC monolayer at the air/water interface

Pulmonary surfactant, a complex mixture of phospholipids and proteins, secreted by the type II epithelial cells in the lungs, is crucial to reducing the effort required for breathing. A lack of adequate amounts of pulmonary surfactant in underdeveloped lungs of premature infants results in infant respiratory distress syndrome (RDS). Surfactant replacement therapy (SRT) is the approved method of mitigating the effects of RDS. The development of new SRT regimens requires a fundamental understanding of the links between surfactant biochemistry and functionality. We use a coarse-grained (CG) model to predict the surface pressure-surface concentration relationship (equation of state) for pure DPPC, which is a major component of endogenous and synthetic pulmonary surfactant mixtures. We show that the model can be efficiently used to predict the equation of state in excellent agreement with experimental results. We also study the structure of the monolayer as a function of the surface tension of the system. We show that a decrease in the applied surface tension results in an increase in order in the head group region and a decrease in order in the tail region of DPPC. This technique may be useful in the prediction of equations of state for surfactant replacements.

Computer Studies on the Effects of Long Chain Alcohols on Sodium Dodecyl Sulfate (SDS) Molecules in SDS/Dodecanol and SDS/Hexadecanol Monolayers at the Air/Water Interface

Molecular dynamics simulations of sodium dodecyl sulfate (SDS)/dodecanol and SDS/hexadecanol monolayers at the air/water interface were investigated where the monolayer mixtures were prepared by two different configurations. In the first configuration, all of the dodecanol (or hexadecanol) molecules were placed together and also the SDS molecules were placed together in the surface area. In the second configuration, the dodecanol (or hexadecanol) molecules were uniformly distributed with the SDS molecules, forming a homogeneous mixture. The results showed that the alcohol tails are more ordered and thicker than the SDS tails in monolayers where the alcohol molecules are close to each other and separated from the SDS. However, the reverse trend is observed in monolayers where the SDS and alcohol molecules are well mixed; that is, the alcohol tails seem to have less order. Studies of how the SDS tails are affected by the presence of long chain alcohols are also discussed. Basically, by increasing the alcohol chain length, the order and the thickness of the SDS tails increased when those molecules were placed all together in a region of the surface area. When both surfactants were well mixed, the order and thickness of the SDS chains decreased as the alcohol chain length increased. Comparisons of the present results with actual experiments of similar systems were performed, and they showed similar tendencies.

DPPC Langmuir monolayer at the air-water interface: probing the tail and head groups by vibrational sum frequency generation spectroscopy

Dipalmitoylphosphatidylcholine (DPPC) is the predominant lipid component in lung surfactant. In this study, the Langmuir monolayer of deuterated dipalmitoylphosphatidylcholine (DPPC-d62) in the liquid-expanded (LE) phase and the liquid-condensed (LC) phase has been investigated at the air-water interface with broad bandwidth sum frequency generation (BBSFG) spectroscopy combined with a Langmuir film balance. Four moieties of the DPPC molecule are probed by BBSFG: the terminal methyl (CD3) groups of the tails, the methylene (CD2) groups of the tails, the choline methyls (CH3) in the headgroup, and the phosphate in the headgroup. BBSFG spectra of the four DPPC moieties provide information about chain conformation, chain orientation, headgroup orientation, and headgroup hydration. These results provide a comprehensive picture of the DPPC phase behavior at the air-water interface. In the LE phase, the DPPC hydrocarbon chains are conformationally disordered with a significant number of gauche configurations. In the LC phase, the hydrocarbon chains are in an all-trans conformation and are tilted from the surface normal by 25 degrees. In addition, the orientations of the tail terminal methyl groups are found to remain nearly unchanged with the variation of surface area. Qualitative analysis of the BBSFG spectra of the choline methyl groups suggests that these methyl groups are tilted but lie somewhat parallel to the surface plane in both the LE and LC phases. The dehydration of the phosphate headgroup due to the LE-LC phase transition is observed through the frequency blue shift of the phosphate symmetric stretch in the fingerprint region. In addition, implications for lung surfactant function from this work are discussed.

Structure of the nonionic surfactant triethoxy monooctylether C8E3 adsorbed at the free water surface, as seen from surface tension measurements and Monte Carlo simulations

The adsorption layer of the nonionic surfactant triethoxy monooctylether C8E3 has been investigated at the free water surface by means of both experimental and computer simulation methods. The surface tension of the aqueous solution of C8E3 has been measured by pendant drop shape analysis in the entire concentration range in which C8E3 is soluble in water. The data obtained from these measurements are used to derive the adsorption isotherm. The critical micellar concentration and the surface excess concentration of the saturated adsorption layer are found to be 7.48 mM and 4.03 micromol/m2, respectively, the latter value corresponding to the average area per molecule of 41 A2. In order to analyze the molecular level structure of the unsaturated adsorption layer, Monte Carlo simulations have been performed at four different surface concentration values, i.e., 0.68, 1.36, 2.04, and 2.72 micromol/m2, respectively. It has been found that the water surface is already almost fully covered at the lowest surface density value investigated, and the adsorbed molecules show a strong preference for lying parallel with the interface in elongated conformations. No sign of the penetration of the hydrophilic triethoxy headgroups into the aqueous phase to any extent has been observed. With increasing surface densities the preferential orientation of the apolar octyl tails gradually turns from lying parallel with the interface to pointing toward the vapor phase by their CH3 end, whereas the conformation of the adsorbed molecules becomes gradually less elongated. Both of these changes lead to the increase of the number of C8E3 molecules being in a direct contact (i.e., forming hydrogen bonds) with water. However, the increasing number of the C8E3 molecules hydrogen bonded to water is found to be accompanied by the weakening of this binding, i.e., the decrease of both the number of hydrogen bonds a bound C8E3 molecule forms with water and the magnitude of the average binding energy of the adsorbed C8E3 molecules.

Adsorption of 1-octanol at the free water surface as studied by Monte Carlo simulation

The adsorption of 1-octanol at the free water surface has been investigated by Monte Carlo computer simulation. Six different systems, built up by an aqueous and a vapor phase, the latter also containing various number of octanol molecules, have been simulated. The number of the octanol molecules has been chosen in such a way that the octanol surface density varies in a broad range, between 0.27 and 7.83 micromol/m(2) in the six systems simulated. For reference, the interfacial system containing bulk liquid octanol in the apolar phase has also been simulated. The results have shown that the formation of hydrogen bonds between the interfacial water and adsorbed octanol molecules is of key importance in determining the properties of the adsorbed layer. At low octanol surface concentration values all the octanol molecules are strongly (i.e., by hydrogen bonds) bound to the aqueous phase, whereas their interaction with each other is negligibly small. Hence, they are preferentially oriented in such a way that their own binding energy (and thus their own free energy) is minimized. In this preferred orientation the O-H bond of the octanol molecule points flatly toward the aqueous phase, declining by about 30 degrees from the interfacial plane, irrespectively from whether the octanol molecule is the H-donor or the H-acceptor partner in the hydrogen bond. Hence, in its preferred orientation the octanol molecule can form at least two low energy hydrogen bonds with water: one as a H-donor and another one as a H-acceptor. Moreover, the preferred orientation of the hydrogen bonded water partners is close to one of the two preferred interfacial water alignments, in which the plane of the water molecule is parallel with the interface. When increasing the octanol surface density, the water surface gets saturated with hydrogen bonded octanols, and hence any further octanol molecule can just simply condense to the layer of the adsorbed octanols. The surface density value at which this saturation occurs is estimated to be about 1.7 micromol/m(2). Above this surface density value the hydrogen bonded octanols and their water partners are oriented in such a way that the number of the water-octanol hydrogen bonds is maximized. Hence, the preferred alignment of the O...O axes of these hydrogen bonds is perpendicular to the interface. This orientation is far from the optimal alignment of the individual octanol molecules, which is also reflected in the observed fact that, unlike in the case of many other adsorbents, the average molecular binding energy of the adsorbed octanol molecules increases (i.e., becomes less negative) with increasing octanol surface density.