Next Article in Journal
Micron-Sized Pored Membranes Based on Polyvinylidene Difluoride Hexafluoropropylene Prepared by Phase Inversion Techniques
Next Article in Special Issue
Synthesis and Application of Aurophilic Poly(Cysteine) and Poly(Cysteine)-Containing Copolymers
Previous Article in Journal
Incorporation of Calcium Containing Mesoporous (MCM-41-Type) Particles in Electrospun PCL Fibers by Using Benign Solvents
Previous Article in Special Issue
Vesicles in Multiple Shapes: Fine-Tuning Polymersomes’ Shape and Stability by Setting Membrane Hydrophobicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cluster Formation of Polyphilic Molecules Solvated in a DPPC Bilayer

by
Xiang-Yang Guo
,
Christopher Peschel
,
Tobias Watermann
,
Guido Falk von Rudorff
and
Daniel Sebastiani
*
Institute of Chemistry, MLU Halle-Wittenberg, von-Danckelmann-Platz 4, 06120 Halle, Germany
*
Author to whom correspondence should be addressed.
Polymers 2017, 9(10), 488; https://doi.org/10.3390/polym9100488
Submission received: 15 September 2017 / Revised: 29 September 2017 / Accepted: 29 September 2017 / Published: 6 October 2017
(This article belongs to the Special Issue From Amphiphilic to Polyphilic Polymers)

Abstract

:
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.

1. Introduction

Molecules which have the ability of self-assembling are of huge interest for biochemical (lipid bilayers), and nanosized materials [1]. To understand such behavior is crucial for the rational design of model and advanced systems. Many molecule types have been inserted into lipid bilayers in experiment and in simulations to understand their behavior. This ranges from early attempts by inserting alkanes into a lipid bilayer system to see where they are located inside the system [2], up to the latest investigations of fluorinated alkanes and alcohols [3]. In extend to those publications we investigated a very complex compound combining many philicities in one molecule. These so called polyphilic molecules, as the name suggests, are compounds consisting of fragments with different philicities. Special attention is paid here to polyphiles that contain a rigid rod-shaped aromatic core with opposing end groups and lateral groups with different philicities. The two end groups are typically highly polar, allowing the formation of multiple hydrogen bonds, which is only possible in the headgroup region of the bilayer and in the aqueous phase, whilst the lateral groups are alkyl, partially fluorinated or perfluorinated chains [4]. This novel class of molecules have received significant attention in recent years [5,6,7,8,9,10]. Polyphilic molecules can be used to modify the phase transitions temperature of a lipid bilayer [11] and also serve as a drug delivery agent [12,13]. Furthermore, they cause effects like compression or stretching of bilayer systems which was recently shown [11,14,15,16,17,18]. Fluorocarbon compounds are also studied for influencing the metabolism of rats [19,20] and for in vitro synthesis of lipid bilayer proteins [21].
To fully understand why this type of molecules influences the bilayer properties, the polyphilic molecules themselves have to be studied on a molecular scale. The molecule (B16/10) being investigated in this work possesses three philicities, namely fluorophilic, hydrophilic and lipophilic parts [6,7,8,9,10,22]. The rigid aromatic phenylene–ethynylene-backbone forms its frame structure, and two hydrophilic groups are terminating this backbone chain. In the middle of the backbone, two side chains are attached, one of which is a perfluoro-n-alkane, the other is a regular n-alkane (see Figure 1). The length of the side chains can be derived from the name (16 carbon atoms for the alkyl chain and 10 carbon atoms for the perfluorinated alkyl groups). This molecular structure ensures a trans-bilayer orientation of the backbone, yielding an anchor point for the alkyl/perfluoroalkyl chains at the center of the bilayer. Mobile aliphatic side chains and polar end groups offer many possibilities of functionalization to tailor the interactions with the bilayers. As the fundamental building blocks of cellular membranes, phospholipid bilayers play a decisive role in many of their biological functions. The modification of lipid bilayer functions via interactions with biomolecules such as proteins or peptides has been widely investigated [23,24,25,26]. Due to experimental investigations of the influence of purely synthetic molecules on DPPC bilayers [14,15,27], it serves perfectly as a model system for thorough investigation of the influence of our polyphilic molecules on the bilayer properties. The concept of polyphilicity was even recently used for directly modifying the lipid bilayer itself by perfluorinating the end of the lipid tails [28].
The available experimental results show that, when bolapolyphile molecules (BP) are incorporated into gel phase lipid (DPPC) bilayers, the formation of large BP domains within the bilayer and a separation into different lamellar species can be observed [16,17]. The thermal behavior of the lipid bilayers was drastically altered upon BP incorporation and several endothermic transitions above Tm of pure DPPC bilayer occurred [11]. In the liquid phase, the BPs were homogeneously distributed in the lipid bilayer plane [11,16].
Lipid bilayer simulations have reached an exciting point, where the time and length scales of simulations are approaching experimental resolutions and can be used to interpret experiments on increasingly complex model bilayers. Within molecular dynamics (MD) of hybrid-bilayer systems, one is able to get an insight into the dynamical behavior on a molecular scale. Therefore, these simulations provide complementary information to experiments [25,29,30,31,32,33,34]. Furthermore, they can yield molecular-level insight into the structure and dynamics of these systems with a spatial resolution and time-scale that may not be feasible experimentally. Summarized, they serve as rich sources of quantitative data on molecular flexibility, lipid diffusion, ordering and atomic interactions. A detailed understanding of lipid bilayer properties is necessary to fully understand their important biophysical characteristics.
Previously, we have reported a MD study of one single B16/10 molecule inserted as a trans-bilayer agent into a dipalmitoylphosphatidylcholine (DPPC) bilayer [35]. The results showed that B16/10 is commensurate with the bilayer, and at the same time, a certain intramolecular bending and an inclination with respect to the bilayer plane is observed. While the lipophilic groups remain bilayer-centered, the fluorophilic parts tend to orient towards the phosphate headgroups, which is due to a slight size mismatch between the bilayer and the lipophilic backbone of the molecule.
In the present work, MD simulations with full atomic detail were carried out for a small trans-bilayer cluster containing six B16/10 molecules embedded into a DPPC bilayer on a large time scale. The difference of the dynamical and structural properties of the DPPC bilayer between pure bilayer system and the mixture of the lipid bilayer with B16/10 molecules is illustrated. The configurations and dynamics of incorporated B16/10 are studied. The results are compared with available experimental and literature data.

2. System Setup and Computational Details

A small cluster of six B16/10 molecules was embedded into a DPPC bilayer consisting 288 lipid molecules (144 per leaflet). The mixed molecular system was hydrated using 8756 TIP3 water molecules. The resulting periodic box has a dimension of 98 × 98 × 68 Å 3 . A snapshot of the system is shown in Figure 2.
The initial configurations of B16/10 molecules was selected based on the experimentally observed phase formation of similar structured molecules. According to the experimental results, similar structured B12 molecules spontaneously self-organize in lipid bilayers (DPPC), forming ordered snowflake like structures with 6-fold symmetry in giant unilamellar vesicles [16,17]. Furthermore, the initial structure completely vanished after an equilibration time of 20 ns. The first 20 ns have been omitted in all analyses. The structure after the equilibration is shown in Figure 3 as a top view onto the system.
All MD simulations in this work were performed by software package namd 2.9 using the CHARMM force field [32,36,37]. The detailed parameters for B16/10 are presented in a previous work of our group [38]. Three dimensional periodic boundary conditions were used. The system is kept at a constant pressure of 1 bar and a constant specified temperature (isobaric-isothermal NpT ensemble) using a modified Nose–Hoover method in which Langevin dynamics is used to control fluctuations in the barostat. The semi-isotropic pressure coupling was applied separately for the bilayer plane and bilayer normal with a coupling constant of 1 bar. Experimentally measured phase transition temperature of a DPPC bilayer is between 313 and 315 K [16,17]. It is experimentally observed that the presence of polyphile molecules increases this phase transition temperature. In order to avoid a simulation in a gel phase, we set the simulation temperature to 335 K in this work.
The system was simulated for 1 μ s with a time step of 2 fs. The bond lengths were constrained using the SHAKE algorithm. The vdW cutoff follows the force field specifications for lipid bilayers. A particle mesh Ewald summation was used to calculate the electrostatic interactions. Cutoff radius for van der Waals interactions was set to 1.0 nm. Particle mesh Ewald (PME) summations were applied for long-range electrostatic interactions with a grid spacing of 0.12 nm and a cutoff radius of 1.0 nm was employed for real space summation. For the purpose of comparison, simulations of a pure DPPC bilayer was carried out under same conditions. The TIP3P water model was used to solvate the system. Data analysis of the trajectories were done by using VMD plugins [39], the python module MDAnalysis [40], the freeware program package TRAVIS [41] and our own codes.

3. Results and Discussion

Statistical analysis was carried out to characterize the dynamical and structural properties of B16/10 trans-bilayer molecules inside a DPPC bilayer. In particular, we (i) calculated the lateral diffusion coefficients of DPPC and B16/10 molecules and (ii) investigated the axial location and orientation of B16/10 molecules and internal structure of B16/10 cluster. The results are compared with a pure DPPC bilayer system, available experimental data and former results of our group.

3.1. Lateral Diffusion

Investigation of DPPC and B16/10 molecules lateral mobility in a planar lipid bilayer are carried out and the results are compared with the pure bilayer system. The averaged self-diffusion coefficients of B16/10 and DPPC molecules are calculated from the mean squared displacement (MSD) r 2 ( t ) using the Einstein relation:
M S D ( t ) = r 2 ( t ) = [ r ( t ) r ( 0 ) ] 2
D = 1 2 d l i m t d d t [ r ( t ) r ( 0 ) ] 2 ,
where d is the number of dimension which is 2 for our calculation. The MSD is calculated with respect to the center of mass (COM) of all lipid molecules to avoid artifacts from water layer movement. The averaged MSD for DPPC and B16/10 molecules are shown in Figure 4.
The calculated self-diffusion coefficient of pure DPPC at 335K as presented in Table 1 is in very good agreement with the experimentally determined value of 14.2 ± 1.2 × 10 12 m 2 / s [42]. The diffusion coefficient of B16/10 molecules is slightly lower than the diffusion coefficient of the lipids in a pure bilayer. For the mixed system, we find high diffusivity values for the lipid molecules, which is unexpected.

3.2. B16/10 Radial Distribution Functions

We calculated the pair correlation function, g ( r ) , as the probability of finding a pair of B16/10 molecules at distance r apart, relative to the probability expected for a completely random distribution at the same density [43]. The g ( r ) of the center of mass for each central phenylene ring of a B16/10 molecule (COR-COR) and of its terminal groups (CH3 and CF3) is shown in Figure 5.
We observe a clustering effect documented by the two peak nature in the radial distribution function. The first peak at 7 Å distance shows that two B16/10 molecules are directly adjacent to each other, as there is no space for interlaced molecules. Nevertheless the lipid molecules force the B16/10 molecules to tilt a little away from a coplanar orientation and causing this distance of 7 Å. Due to the tilt and high distance, pi–pi stacking does not frequently occur. Nevertheless, the radial distribution function has also non-zero values below 7 Å, which at least can be interpreted as possible pi–pi stacking. As there are only six B16/10 molecules in the system, effects of a big ensemble like in experiment [11,17] cannot be captured and therefore a strong pi–pi stacking behavior cannot be excluded even though it is unfavorable in our simulation. The second peak at 12 Å shows that a third B16/10 molecule is forming a little cluster with the other two. As the second peak is not double of the first one it can be seen that they are not completely aligned but forming a triangle. The clustering can be seen in Figure 6. Comparing Figure 3 and Figure 6, one can see that the molecules inside the cluster are exchanging within the simulation.
Having this in mind it is even more surprising that the diffusion coefficient discussed in the previous section of 11 × 10 12 m 2 s for B16/10 (see Table 1) is rather high. This clustering behavior of x-shaped molecules like the non-perfluorinated B12 molecule was also observed experimentally [11,17]. Therefore it can be assumed that by perfluorinating one side chain, the clustering behavior is maintained. For the terminal groups CH3 and CF3 there is only one well defined peak at around 3 Å for CH3 and at 5 Å for CF3. Beyond this there is a smooth decay to unity apart from residual structure but with no recognizable pattern. All in all this shows that there is a correlation of the backbones of the B16/10 molecules forming a little cluster whereas there is no correlation of the side chains aside from pairwise orientation.

3.3. B16/10 Backbone Angles Distribution

We use three angels α 1 , α 2 and β to characterize the orientation of B16/10 backbone within the bilayer as shown in Figure 7. The distributions of the three angels are displayed in Figure 8.
The calculated angular distributions of the B16/10 backbone are in very good agreement with our former results of a single B16/10 molecule inside a DPPC bilayer [35]. The B16/10 molecules are well incorporated into the lipid bilayer, cross the whole bilayer with their well-matched hydrophobic core length and adapting into the bilayer by a slight tilt.

3.4. B16/10 Terminal Group Integration and Side Chain Orientations

In Figure 9 the averaged position distribution of CF3 and CH3 terminal groups of B16/10 molecules relative to the center of mass of the system is shown. For avoiding artifacts, the center of mass movement of the whole system has been subtracted. The center of the lipid bilayer is denoted by z = 0 and z is the axis along the normal of the bilayer plane. The bilayer itself has a thickness of approximately 40 ± 2.5 Å. As one can see in Figure 9, both side chains approach the head group region of the bilayer. Nevertheless the CF side chain tends to stay more in the headgroup region than the CH chain. The combination of higher solubility inside the leaflets (which was recently shown by Brehm et al. [3]) and sterical issues force the CF3 terminal group to stay in the headgroup region as in direct comparison the CF side chain is nearly inflexible whereas the CH side is more flexible and thinner which makes it easier for the CH chain to find its way through the lipid bilayer. Not surprising is, that the CH3 terminal group has a high occurrence in the middle of the bilayer, as this the most non polar region of the bilayer which has been observed in experiment before for alkanes [2]. Of course sterical issues are more dominant because not every polyphile can arrange like this when incorporated into a cluster and therefore it also stays frequently inside the leaflets. On first sight the occurrence of the CF3 terminal group in the middle of the bilayer seems astonishing. It is easy explained by attempts to flip the sides of the bilayer. As it is only ninth of the occurrence of the highest peak this rarely happens and is never a stationary state. By closer inspection we found that the flipping of the CF side chain happens on a nanosecond time scale.
The distribution of the side chain angles can be used to confirm these statements as seen in Figure 10. The angles are defined as the angle of the vector from the terminal groups to the center of the phenylene ring in the middle of the backbone and the bilayer normal. Therefore 90° denotes directly aligned between the leaflets in the bilayer plane. The prominent angles around 30° and 150° for the CF3 terminal group confirms that it mostly stays in head group region of the bilayer. The very slight occurrence around 90° also shows the behavior to flip rather than staying between the leaflets. More or less Figure 9 and Figure 10 resemble each other in appearance and therefore also confirm the statements for CH terminal head group mentioned above.

4. Conclusions

In this paper, we report a study of six polyphilic molecules embedded in a lipid bilayer/water system. Our results give insights not only into the the conformational preference of the B16/10 molecules inside the DPPC bilayer, but also on the dynamics of the B16/10 molecules in the DPPC bilayer environment. Within 1 ms of simulation the stability of the bilayer maintained upon insertion of the additive B16/10 molecules. The diffusion of the lipid molecules is increased compared to the pure lipid bilayer and B16/10 molecules move only slightly slower than the lipids in the pure bilayer. For the intermolecular structure of the B16/10 molecules a clustering effect can be observed in the RDFs. The difference in the side chain orientation and configuration is also in good agreement with the simulation results we obtained before [35]. In general, the conclusions drawn from simulations are consistent with experimentally observed effects [16,17,42]. We anticipate that these findings will be important for understanding the role of polyphilic molecules in modulating and modifying bilayer properties.

Acknowledgments

This work has been supported by the German Research Foundation (DFG) within the Forschergruppe FOR1145 (grant number Se 1008/9-1). The authors thank the China Scholarship Council (CSC) for providing a scholarship.

Author Contributions

Xiang-Yang Guo, Tobias Watermann and Guido Falk von Rudorff planned and set up the simulation. Xiang-Yang Guo and Christopher Peschel analyzed the data and wrote the manuscript. Daniel Sebastiani supervised this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ariga, K.; Hill, J.P.; Lee, M.V.; Vinu, A.; Charvet, R.; Acharya, S. Challenges and breakthroughs in recent research on self-assembly. Sci. Technol. Adv. Mater. 2008, 9, 014109. [Google Scholar] [CrossRef] [PubMed]
  2. White, S.H.; King, G.I.; Cain, J.E. Location of hexane in lipid bilayers determined by neutron diffraction. Nature 1981, 290, 161–163. [Google Scholar] [CrossRef]
  3. Brehm, M.; Saddiq, G.; Watermann, T.; Sebastiani, D. Influence of Small Fluorophilic and Lipophilic Organic Molecules on Dipalmitoylphosphatidylcholine Bilayers. J. Phys. Chem. B 2017, 121, 8311–8321. [Google Scholar] [CrossRef] [PubMed]
  4. Bates, M.; Walker, M. Dissipative particle dynamics simulation of T- and X-shaped polyphilic molecules exhibiting honeycomb columnar phases. Soft Matter 2009, 5, 346–353. [Google Scholar] [CrossRef]
  5. Crane, A.J.; Martínez-Veracoechea, F.J.; Escobedo, F.A.; Müller, E.A. Molecular dynamics simulation of the mesophase behaviour of a model bolaamphiphilic liquid crystal with a lateral flexible chain. Soft Matter 2008, 4, 1820–1829. [Google Scholar] [CrossRef]
  6. Hill, E.; Stratton, K.; Whitten, D.G.; Evans, D.G. Molecular Dynamics Simulation Study of the Interaction of Cationic Biocides with Lipid Bilayers: Aggregation Effects and Bilayer Damage. Langmuir 2012, 28, 14849–14854. [Google Scholar] [CrossRef] [PubMed]
  7. Schulz, M.; Olubummo, A.; Binder, W.H. Beyond the lipid-bilayer: Interaction of polymers and nanoparticles with membranes. Soft Matter 2012, 8, 4849–4864. [Google Scholar] [CrossRef]
  8. Hinks, J.; Wang, Y.; Poh, W.H.; Donose, B.C.; Thomas, A.W.; Wuertz, S.; Loo, S.C.J.; Bazan, G.C.; Kjelleberg, S.; Mu, Y.; et al. Modeling Cell Membrane Perturbation by Molecules Designed for Transmembrane Electron Transfer. Langmuir 2014, 30, 2429–2440. [Google Scholar] [CrossRef] [PubMed]
  9. James, P.V.; Sudeep, P.K.; Suresh, C.H.; Thomas, K.G. Photophysical and Theoretical Investigations of Oligo(p-phenyleneethynylene)s: Effect of Alkoxy Substitution and Alkyne–Aryl Bond Rotations. J. Phys. Chem. A 2006, 110, 4329–4337. [Google Scholar] [CrossRef] [PubMed]
  10. Hill, E.H.; Whitten, D.G.; Evans, D.G. Computational Study of Bacterial Membrane Disruption by Cationic Biocides: Structural Basis for Water Pore Formation. J. Phys. Chem. B 2014, 118, 9722–9732. [Google Scholar] [CrossRef] [PubMed]
  11. Lechner, B.-D.; Ebert, H.; Prehm, M.; Werner, S.; Meister, A.; Hause, G.; Beerlink, A.; Saalwächter, K.; Bacia, K.; Tschierske, C.; et al. Temperature-Dependent In-Plane Structure Formation of an X-Shaped Bolapolyphile within Lipid Bilayers. Langmuir 2015, 31, 2839–2850. [Google Scholar] [CrossRef] [PubMed]
  12. Mukhopadhyay, S.; Maitra, U. Chemistry and biology of bile acids. Curr. Sci. 2004, 87, 1666–1683. [Google Scholar]
  13. Krafft, M.P. Fluorocarbons and fluorinated amphiphiles in drug delivery and biomedical research. Adv. Drug Deliv. Rev. 2001, 47, 209–228. [Google Scholar] [CrossRef]
  14. Schwieger, C.; Achilles, A.; Scholz, S.; Rüger, J.; Bacia, K.; Saalwächter, K.; Kressler, J.; Blume, A. Binding of amphiphilic and triphilic block copolymers to lipid model membranes: The role of perfluorinated moieties. Soft Matter 2014, 10, 6147–6160. [Google Scholar] [CrossRef] [PubMed]
  15. Scholtysek, P.; Achilles, A.; Hoffmann, C.-V.; Lechner, B.-D.; Meister, A.; Tschierske, C.; Saalwächter, K.; Edwards, K.; Blume, A. A T-Shaped Amphiphilic Molecule Forms Closed Vesicles in Water and Bicelles in Mixtures with a Membrane Lipid. J. Phys. Chem. B 2012, 116, 4871–4878. [Google Scholar] [CrossRef] [PubMed]
  16. Achilles, A.; Bärenwald, R.; Lechner, B.; Werner, S.; Ebert, H.; Tschierske, C.; Blume, A.; Bacia, K.; Saalwächter, K. Self-Assembly of X-Shaped Bolapolyphiles in Lipid Membranes: Solid-State NMR Investigations. Langmuir 2016, 32, 673–682. [Google Scholar] [CrossRef] [PubMed]
  17. Werner, S.; Ebert, H.; Lechner, B.-D.; Lange, F.; Achilles, A.; Bärenwald, R.; Poppe, S.; Blume, A.; Saalwächter, K.; Tschierske, C.; et al. Dendritic Domains with Hexagonal Symmetry Formed by X-Shaped Bolapolyphiles in Lipid Membranes. Chem. Eur. J. 2015, 21, 8840–8850. [Google Scholar] [CrossRef] [PubMed]
  18. Bärenwald, R.; Achilles, A.; Lange, F.; Mendes, T.F.; Saalwächter, K. Applications of Solid-State NMR Spectroscopy for the Study of Lipid Membranes with Polyphilic Guest (Macro)Molecules. Polymers 2016, 8, 439. [Google Scholar] [CrossRef]
  19. Permadi, H.; Lundgren, B.; Andersson, K.; DePierre, J.W. Effects of perfluoro fatty acids on xenobiotic-metabolizing enzymes, enzymes which detoxify reactive forms of oxygen and lipid peroxidation in mouse liver. Biochem. Pharmacol. 1992, 44, 1183–1191. [Google Scholar] [CrossRef]
  20. Heuvel, J.P.V.; Kuslikis, B.I.; Van Rafelghem, M.J.; Peterson, R.E. Tissue distribution, metabolism, and elimination of perfluorooctanoic acid in male and female rats. J. Biochem. Toxicol. 1991, 6, 83–92. [Google Scholar] [CrossRef]
  21. Park, K.-H.; Berrier, C.; Lebaupain, F.; Pucci, B.; Popot, J.-L.; Ghazi, A.; Zito, F. Fluorinated and hemifluorinated surfactants as alternatives to detergents for membrane protein cell-free synthesis. Biochem. J. 2007, 403, 183–187. [Google Scholar] [CrossRef] [PubMed]
  22. Brehm, M.; Weber, H.; Thomas, M.; Hollóczki, O.; Kirchner, B. Domain Analysis in Nanostructured Liquids: A Post-Molecular Dynamics Study at the Example of Ionic Liquids. ChemPhysChem 2015, 16, 3271–3277. [Google Scholar] [CrossRef] [PubMed]
  23. Nagle, J.F.; Tristram-Nagle, S. Structure of lipid bilayers. Biochim. Biophys. Acta 2000, 1469, 159–195. [Google Scholar] [CrossRef]
  24. Liu, B.; Hoopes, M.I.; Karttunen, M. Molecular Dynamics Simulations of DPPC/CTAB Monolayers at the Air/Water Interface. J. Phys. Chem. B 2014, 118, 11723–11737. [Google Scholar] [CrossRef] [PubMed]
  25. Repáková, J.; Holopainen, J.M.; Morrow, M.R.; McDonald, M.C.; Capková, P.; Vattulainen, I. Influence of DPH on the Structure and Dynamics of a DPPC Bilayer. Biophys. J. 2005, 88, 3398–3410. [Google Scholar] [CrossRef] [PubMed]
  26. Hughes, Z.E.; Mark, A.E.; Mancera, R.L. Molecular Dynamics Simulations of the Interactions of DMSO with DPPC and DOPC Phospholipid Membranes. J. Phys. Chem. B 2012, 116, 11911–11923. [Google Scholar] [CrossRef] [PubMed]
  27. Zeng, X.; Kieffer, R.; Glettner, B.; Nürnberger, C.; Liu, F.; Pelz, K.; Prehm, M.; Baumeister, U.; Hahn, H.; Lang, H.; et al. Complex Multicolor Tilings and Critical Phenomena in Tetraphilic Liquid Crystals. Science 2011, 331, 1302–1306. [Google Scholar] [CrossRef] [PubMed]
  28. Yoshino, M.; Kikukawa, T.; Takahashi, H.; Takagi, T.; Yokoyama, Y.; Amii, H.; Baba, T.; Kanamori, T.; Sonoyama, M. Physicochemical Studies of Bacteriorhodopsin Reconstituted in Partially Fluorinated Phosphatidylcholine Bilayers. J. Phys. Chem. B 2013, 117, 5422–5429. [Google Scholar] [CrossRef] [PubMed]
  29. Bennett, W.F.D.; Tieleman, D.P. Computer simulations of lipid membrane domains. Biochim. Biophys. Acta 2013, 1828, 1765–1776. [Google Scholar] [CrossRef] [PubMed]
  30. Lehmler, H.-J.; Bummer, P.M. Mixing of perfluorinated carboxylic acids with dipalmitoylphosphatidylcholine. Biochim. Biophys. Acta 2004, 1664, 141–149. [Google Scholar] [CrossRef] [PubMed]
  31. Korlach, J.; Schwille, P.; Webb, W.W.; Feigenson, G.W. Characterization of lipid bilayer phases by confocal microscopy and fluorescence correlation spectroscopy. Proc. Natl. Acad. Sci. USA 1999, 96, 8461–8466. [Google Scholar] [CrossRef] [PubMed]
  32. Klauda, J.B.; Venable, R.M.; Freites, J.A.; O’Connor, J.W.; Tobias, D.J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A.D.; Pastor, R.W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114, 7830–7843. [Google Scholar] [CrossRef] [PubMed]
  33. Krafft, M.P. Controlling phospholipid self-assembly and film properties using highly fluorinated components—Fluorinated monolayers, vesicles, emulsions and microbubbles. Biochimie 2012, 10, 11–25. [Google Scholar] [CrossRef] [PubMed]
  34. Feller, S.E.; Venable, R.M.; Pastor, R.W. Computer Simulation of a DPPC Phospholipid Bilayer: Structural Changes as a Function of Molecular Surface Area. Langmuir 1997, 13, 6555–6561. [Google Scholar] [CrossRef]
  35. Von Rudorff, G.F.; Watermann, T.; Guo, X.Y.; Sebastiani, D. Conformational Space of a Polyphilic Molecule with a Fluorophilic Side Chain Integrated in a DPPC Bilayer. J. Comput. Chem. 2017, 38, 576–583. [Google Scholar] [CrossRef] [PubMed]
  36. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kalé, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [PubMed]
  37. Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; et al. CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields. J. Comput. Chem. 2010, 31, 671–690. [Google Scholar] [CrossRef] [PubMed]
  38. Von Rudorff, G.F.; Watermann, T.; Sebastiani, D. Perfluoroalkane Force Field for Lipid Membrane Environments. J. Phys. Chem. B. 2014, 118, 12531–12540. [Google Scholar] [CrossRef] [PubMed]
  39. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  40. Michaud-Agrawal, N.; Denning, E.J.; Woolf, T.B.; Beckstein, O. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations. J. Comput. Chem. 2011, 32, 2319–2327. [Google Scholar] [CrossRef] [PubMed]
  41. Brehm, M.; Kirchner, B. TRAVIS—A Free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics Trajectories. J. Chem. Inf. Model. 2011, 51, 2007–2023. [Google Scholar] [CrossRef] [PubMed]
  42. Junglas, M.; Danner, B.; Bayerl, T.M. Molecular Order Parameter Profiles and Diffusion Coefficients of Cationic Lipid Bilayers on a Solid Support. Langmuir 2003, 19, 1914–1917. [Google Scholar] [CrossRef]
  43. McQuarrie, D.A. Statistical Mechanics; University Science Books: Davis, CA, USA, 1976; p. 641. [Google Scholar]
Figure 1. Molecular structure of the polyphilic molecule B16/10 and dipalmitoylphosphatidylcholine (DPPC) molecule studied in this report.
Figure 1. Molecular structure of the polyphilic molecule B16/10 and dipalmitoylphosphatidylcholine (DPPC) molecule studied in this report.
Polymers 09 00488 g001
Figure 2. A snapshot of the simulated system containing six B16/10 molecules, 288 DPPC molecules and 8756 water molecules. Periodic boundary conditions were used in all directions. Atoms of B16/10 molecules and Phosphorus atoms of DPPC headgroups are represented by solid spheres. Phosphorus atoms are green, Fluorine atom of B16/10 side chains are pink, carbon atoms are cyan. All the lipid tails of DPPC are represented by cyan lines. One can see, that the backbone of the polyphile (cyan spheres) is bent and not straight.
Figure 2. A snapshot of the simulated system containing six B16/10 molecules, 288 DPPC molecules and 8756 water molecules. Periodic boundary conditions were used in all directions. Atoms of B16/10 molecules and Phosphorus atoms of DPPC headgroups are represented by solid spheres. Phosphorus atoms are green, Fluorine atom of B16/10 side chains are pink, carbon atoms are cyan. All the lipid tails of DPPC are represented by cyan lines. One can see, that the backbone of the polyphile (cyan spheres) is bent and not straight.
Polymers 09 00488 g002
Figure 3. A snapshot of the simulated system as a top view after an equilibration time of 20 ns. The polyphiles are colored and a clustering can be observed for some of the molecules. The simulation cell dimensions are shown as blue rectangle.
Figure 3. A snapshot of the simulated system as a top view after an equilibration time of 20 ns. The polyphiles are colored and a clustering can be observed for some of the molecules. The simulation cell dimensions are shown as blue rectangle.
Polymers 09 00488 g003
Figure 4. Mean square displacement of lipids in a pure DPPC bilayer (red), the lipids (blue) and B16/10 molecules (green) in the mixed system (lipids/B16/10).
Figure 4. Mean square displacement of lipids in a pure DPPC bilayer (red), the lipids (blue) and B16/10 molecules (green) in the mixed system (lipids/B16/10).
Polymers 09 00488 g004
Figure 5. Radial distribution functions of the center of the central phenylene ring (COR), CF3 and CH3 terminal groups of B16/10 molecules within the DPPC bilayer.
Figure 5. Radial distribution functions of the center of the central phenylene ring (COR), CF3 and CH3 terminal groups of B16/10 molecules within the DPPC bilayer.
Polymers 09 00488 g005
Figure 6. Snapshot of the top view onto the system showing the clustering of the B16/10 molecules. Namely the red, yellow and black colored molecules form a cluster, whereas the green and orange molecules are surrounding the cluster. Compared to Figure 3 (same coloring of the molecules), one can see that the cluster exchanges the molecules within the simulation. The simulation cell dimensions are shown as blue rectangle.
Figure 6. Snapshot of the top view onto the system showing the clustering of the B16/10 molecules. Namely the red, yellow and black colored molecules form a cluster, whereas the green and orange molecules are surrounding the cluster. Compared to Figure 3 (same coloring of the molecules), one can see that the cluster exchanges the molecules within the simulation. The simulation cell dimensions are shown as blue rectangle.
Polymers 09 00488 g006
Figure 7. Angels between the bilayer normal and B16/10 backbone ( α 1 and α 2 ) and backbone bending angle β .
Figure 7. Angels between the bilayer normal and B16/10 backbone ( α 1 and α 2 ) and backbone bending angle β .
Polymers 09 00488 g007
Figure 8. Probability distribution of backbone angles α 1 , α 2 and β .
Figure 8. Probability distribution of backbone angles α 1 , α 2 and β .
Polymers 09 00488 g008
Figure 9. Position of the CF3 and CH3 terminal groups along the z-axis relative to the center of mass of the whole system.
Figure 9. Position of the CF3 and CH3 terminal groups along the z-axis relative to the center of mass of the whole system.
Polymers 09 00488 g009
Figure 10. Probability distribution of side chains angles.
Figure 10. Probability distribution of side chains angles.
Polymers 09 00488 g010
Table 1. Lateral self-diffusion coefficients (D) at 335 K obtained from MD simulation (MD) and experimental data from literature (exp) [42]. Uncertainties are given as three times the standard deviation.
Table 1. Lateral self-diffusion coefficients (D) at 335 K obtained from MD simulation (MD) and experimental data from literature (exp) [42]. Uncertainties are given as three times the standard deviation.
D/( 10 12 m 2 / s )Lipids (Pure DPPC)Lipids (mix)B16/10 (mix)
MD13.9 ± 0.324 ± 0.610.7 ± 0.3
exp14.2 ± 1.2

Share and Cite

MDPI and ACS Style

Guo, X.-Y.; Peschel, C.; Watermann, T.; Rudorff, G.F.v.; Sebastiani, D. Cluster Formation of Polyphilic Molecules Solvated in a DPPC Bilayer. Polymers 2017, 9, 488. https://doi.org/10.3390/polym9100488

AMA Style

Guo X-Y, Peschel C, Watermann T, Rudorff GFv, Sebastiani D. Cluster Formation of Polyphilic Molecules Solvated in a DPPC Bilayer. Polymers. 2017; 9(10):488. https://doi.org/10.3390/polym9100488

Chicago/Turabian Style

Guo, Xiang-Yang, Christopher Peschel, Tobias Watermann, Guido Falk von Rudorff, and Daniel Sebastiani. 2017. "Cluster Formation of Polyphilic Molecules Solvated in a DPPC Bilayer" Polymers 9, no. 10: 488. https://doi.org/10.3390/polym9100488

APA Style

Guo, X. -Y., Peschel, C., Watermann, T., Rudorff, G. F. v., & Sebastiani, D. (2017). Cluster Formation of Polyphilic Molecules Solvated in a DPPC Bilayer. Polymers, 9(10), 488. https://doi.org/10.3390/polym9100488

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop