Perfluoroalkane force field for lipid membrane environments

Binding of Per- and Polyfluoroalkyl Substances to β-Lactoglobulin from Bovine Milk

Per- and polyfluoroalkyl substances (PFAS) are known for their high environmental persistence and potential toxicity. The presence of PFAS has been reported in many dairy products. However, the mechanisms underlying the accumulation of PFAS in these products remain unclear. Here, we used native mass spectrometry and molecular dynamics simulations to probe the interactions between 19 PFAS of environmental concern and two isoforms of the major bovine whey protein β-lactoglobulin (β-LG). We observed that six of these PFAS bound to both protein isoforms with low- to mid-micromolar dissociation constants. Based on quantitative, competitive binding experiments with endogenous ligands, PFAS can bind orthosterically and preferentially to β-LG's hydrophobic ligand-binding calyx. β-Cyclodextrin can also suppress binding of PFAS to β-LG owing to the ability of β-cyclodextrin to directly sequester PFAS from solution. This research sheds light on PFAS-β-LG binding, suggesting that such interactions could impact lipid-fatty acid transport in bovine mammary glands at high PFAS concentrations. Furthermore, our results highlight the potential use of β-cyclodextrin in mitigating PFAS binding, providing insights toward the development of strategies to reduce PFAS accumulation in dairy products and other biological systems.

Controlling the self-assembly of perfluorinated surfactants in aqueous environments

Surface active per- and polyfluoroalkyl substances (PFAS) released in the environment generate great concern in the US and worldwide. The sequestration of PFAS amphiphiles from aqueous media can be limited by their strong tendency to form micelles that plug the pores in the adsorbent material, rendering most of the active surface inaccessible. A joint experimental and simulation approach has been used to investigate the structure of perfluorooctanoate ammonium (PFOA) micelles in aqueous solutions, focusing on the understanding of ethanol addition on PFOA micelle formation and structure. Structurally compact and slightly ellipsoidal in shape, PFOA micelles in pure water become more diffuse with increasing ethanol content, and break into smaller PFOA clusters in aqueous solutions with high ethanol concentration. A transition from a co-surfactant to a co-solvent behavior with the increase of ethanol concentration has been observed by both experiments and simulations, while the latter also provide insight on how to achieve co-solvent conditions with other additives. An improved understanding of how to modulate PFAS surfactant self-assembly in water can inform the fate and transport of PFAS in the environment and the PFAS sequestration from aqueous media.

Molecular Dynamics Simulations of Water Condensation on Surfaces with Tunable Wettability

Water condensation plays a major role in a wide range of industrial applications. Over the past few years, many studies have shown interest in designing surfaces with enhanced water condensation and removal properties. It is well known that heterogeneous nucleation outperforms homogeneous nucleation in the condensation process. Because heterogeneous nucleation initiates on a surface at a small scale, it is highly desirable to characterize water-surface interactions at the molecular level. Molecular dynamics (MD) simulations can provide direct insight into heterogeneous nucleation and advance surface designs. Existing MD simulations of water condensation on surfaces were conducted by tuning the solid-water van der Waals interaction energy as a substitute for modeling surfaces with different wettabilities. However, this approach cannot reflect the real intermolecular interactions between the surface and water molecules. Here, we report MD simulations of water condensation on realistic surfaces of alkanethiol self-assembled monolayers with different head group chemistries. We show that decreasing surface hydrophobicity significantly increases the electrostatic forces between water molecules and the surface, thus increasing the water condensation rate. We observe a strong correlation between our rate of condensation results and the results from other surface characterization metrics, such as the interfacial thermal conductance, contact angle, and the molecular-scale wettability metric of Garde and co-workers. This work provides insight into the water condensation process at the molecular scale on surfaces with tunable wettability.

Conformational Space of a Polyphilic Molecule with a Fluorophilic Side Chain Integrated in a DPPC Bilayer

We investigate the conformational space of a polyphilic molecule with hydrophilic, lipophilic and fluorophilic parts inserted as a transmembrane agent into a dipalmitoylphosphatidylcholine bilayer by means of all-atom molecular dynamics simulations. Special focus is put on the competing structural driving forces arising from the hydrophilic, lipophilic and fluorophilic side chains and the aromatic backbone of the polyphile. We observe a significant difference between the lipophilic and the fluorophilic side chains regarding their intramembrane distribution. While the lipophilic groups remain membrane-centered, the fluorophilic parts tend to orient toward the phosphate headgroups. This trend is important for understanding the influence of polyphile agents on the properties of phospholipid membranes. From a fundamental point of view, our computed distribution functions of the side chains are related to the interplay of sterical, enthalpic and entropic driving forces. Our findings illustrate the potential of rationally designed membrane additives which can be exploited to tune the properties of phospholipid membranes. © 2017 Wiley Periodicals, Inc.

Effect of Perfluoroalkyl Endgroups on the Interactions of Tri-Block Copolymers with Monofluorinated F-DPPC Monolayers

We studied the interaction of amphiphilic and triphilic polymers with monolayers prepared from F-DPPC (1-palmitoyl-2-(16-fluoropalmitoyl)-sn-glycero-3-phosphocholine), a phospholipid with a single fluorine atom at the terminus of the sn-2 chain, an analogue of dipalmitoyl-phosphatidylcholine (DPPC). The amphiphilic block copolymers contained a hydrophobic poly(propylene oxide) block flanked by hydrophilic poly(glycerol monomethacrylate) blocks (GP). F-GP was derived from GP by capping both termini with perfluoro-n-nonyl segments. We first studied the adsorption of GP and F-GP to lipid monolayers of F-DPPC. F-GP was inserted into the monolayer up to a surface pressure Π of 42.4 mN m⁻¹, much higher than GP (32.5 mN m⁻¹). We then studied isotherms of lipid-polymer mixtures co-spread at the air-water interface. With increasing polymer content in the mixture a continuous shift of the onset of the liquid-expanded (LE) to liquid-condensed (LC) transition towards higher molecular and higher area per lipid molecule was observed. F-GP had a larger effect than GP indicating that it needed more space. At a Π-value of 32 mN m⁻¹, GP was excluded from the mixed monolayer, whereas F-GP stayed in F-DPPC monolayers up to 42 mN m⁻¹. F-GP is thus more stably anchored in the monolayer up to higher surface pressures. Images of mixed monolayers were acquired using different fluorescent probes and showed the presence of perfluorinated segments of F-GP at LE-LC domain boundaries.

Cluster Formation of Polyphilic Molecules Solvated in a DPPC Bilayer

We analyse the initial stages of cluster formation of polyphilic additive molecules which are solvated in a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer. Our polyphilic molecules comprise an aromatic (trans-bilayer) core domain with (out-of-bilayer) glycerol terminations, complemented with a fluorophilic and an alkyl side chain, both of which are confined within the aliphatic segment of the bilayer. Large-scale molecular dynamics simulations (1 μ s total duration) of a set of six of such polyphilic additives reveal the initial steps towards supramolecular aggregation induced by the specific philicity properties of the molecules. For our intermediate system size of six polyphiles, the transient but recurrent formation of a trimer is observed on a characteristic timescale of about 100 ns. The alkane/perfluoroalkane side chains show a very distinct conformational distribution inside the bilayer thanks to their different philicity, despite their identical anchoring in the trans-bilayer segment of the polyphile. The diffusive mobility of the polyphilic additives is about the same as that of the surrounding lipids, although it crosses both bilayer leaflets and tends to self-associate.

Polyphilic Interactions as Structural Driving Force Investigated by Molecular Dynamics Simulation (Project 7)

We investigated the effect of fluorinated molecules on dipalmitoylphosphatidylcholine (DPPC) bilayers by force-field molecular dynamics simulations. In the first step, we developed all-atom force-field parameters for additive molecules in membranes to enable an accurate description of those systems. On the basis of this force field, we performed extensive simulations of various bilayer systems containing different additives. The additive molecules were chosen to be of different size and shape, and they included small molecules such as perfluorinated alcohols, but also more complex molecules. From these simulations, we investigated the structural and dynamic effects of the additives on the membrane properties, as well as the behavior of the additive molecules themselves. Our results are in good agreement with other theoretical and experimental studies, and they contribute to a microscopic understanding of interactions, which might be used to specifically tune membrane properties by additives in the future.