Data from molecular dynamics simulations for COSMOmic

Introduction

COSMOmic is an extension of COSMO-RS for anisotropic systems such as membranes and micelles, it was introduced by Klamt et al. and detailed information can be found in their original article [DOI: 10.1021/jp801736k ]. The COSMOmic method is implemented in the COSMOtherm program of COSMOlogic.

COSMOmic requires as input the structural composition of the anisotropic system. More precisely, the atomic distribution for each atom type of the system is needed, usually calculated from a snapshot of a molecular dynamics (MD) simulation. COSMOtherm takes as input the coordinates of the atoms as pdb file (single snapshot). However, we have shown that using averaged (over a MD trajectory) atomic distributions are more reliable and scientifically more reasonable [DOI: 10.1002/jcc.23262]. Whereas some groups provide pdb files of lipid membranes from MD simulations, atomic distributions suitable for COSMOtherm or trajectories from which the distributions could be calculated can usually not be found. Therefore, here we provide depth-dependent membrane and micelle compositions averaged over MD simulations carried out in our group. In the following, we assume that you know our article which discusses the connection of MD simulations and COSMOmic:

S. Jakobtorweihen, T. Ingram, and I. Smirnova,
Combination of COSMOmic and molecular dynamics simulations for the calculation of membrane-water partition coefficients,
J. Comput. Chem. 2013, 34, 1332-1340.

Provided data

Below you will find a table that contains data and information on different lipid bilayer systems. A snapshot is shown in the first column. Information about the system and the MD simulation can be found in the second column, the simulation time refers to the time over which the atomic distribution was averaged. If a salt concentration is given, it refers to the concentration in the water phase (outer layers), and not the concentration in the whole system. Averaging was started only after an appropriate equilibration time. Three different files are provided for each system: (1) the averaged atomic distribution in COSMOtherm mic format, (2) a single system snapshot in pdb format, and (3) a lipid/surfactant conformer taken from the MD simulation. The latter is provided, as COSMO-RS needs at least one conformer for each species in the system and it has been shown that conformers of amphiphilic molecules should be taken out of the correct environment (see DOI: 10.1002/jcc.23262 and references therein). Therefore, the lipid conformer was chosen from the lipid bilayer MD trajectory as described in our article [DOI: 10.1002/jcc.23262]. Analog, the surfactant conformers were taken from a micelle simulation. The lipid/surfactant conformer can be used to calculate a cosmo file to be used in the COSMO-RS calculations.

Usage in COSMOtherm

Below you find mic files (COSMOtherm 14 and later) each containing an averaged atomic distribution. Before you can use these files, you have to adjust the list of cosmo files. If you are using COSMOmic via COSMOthermX, you can load the mic files with Read Micelle. If you are using the command line, the mic file name has to be given in the input file (RMIC=).

The data in the mic files is organized as follows: The first line defines variables needed to characterize the bilayer or micelle, see the COSMOmic manual for further details. Note that the temperature given in this line is the temperature the distribution was created at. The temperature of the COSMOmic calculations is specified elsewhere (inp file). On the second line the list of cosmo files starts. For each compound making up the system a cosmo file is needed. Here you have to add the correct name and location of your cosmo files. It is important that you do not change the order of the file list, as it has to correspond to the atomic distributions. Nr lines follow this list of cosmo files, where nr is the number of layers. That is, each line defines one layer, starting from the membrane/micelle center to the water phase. Each line has as many columns as atom types in the bilayer or micelle. For example for a DMPC/water system you will find 121 columns, as DMPC has 118 atoms and water has 3 atoms. In all provided files the first columns are the lipid or surfactant atoms and the last three columns are the water atoms (hence, the third last column represents the oxygen atoms of water). The order of the lipid/surfactant atoms can be seen in the conformer file. Each column sums up to one, as the values are the fractions of the atom types that are found in a specific layer. The lipid/surfactant conformer used in the calculation has to fit to the distributions; the order of the atoms should be the same. If you use your own conformer with the provided data, you have to take care that the atom order is correct. The provided conformer has the correct order. However, you have to generate a cosmo file (e.g., single-point Turbomole calculation) to use it for a COSMO-RS calculation. Please note, the pdb files contain the atom names as used in the MD force field; in addition the element symbol is included in the conformer files. Sometimes older versions of TmoleX have problems with these atom names (ERROR: No radius for … found). If you encounter this problem, you might wane use Open Babel to convert the pdb file to a xyz file and use that as input.

References

If you use averaged atomic distributions or material from this web page, please cite:

Data for lipid bilayers


System Properties Files Comments

DMPC/water

  • water/lipid = 30.6
  • Nlipids=128
  • T=303 K
  • simulation time 80 ns
  • force field: CHARMM36

POPC/water

  • water/lipid = 31.6
  • Nlipids=128
  • T=303 K
  • simulation time 30 ns
  • force field: CHARMM36

DOPC/water

  • water/lipid = 41.9
  • Nlipids=128
  • T=298 K
  • simulation time 30 ns
  • force field: CHARMM36

SOPC/water

  • water/lipid = 30.4
  • Nlipids=128
  • T=303 K
  • simulation time 40 ns
  • force field: CHARMM36

DMPC/water/KCl

  • water/lipid = 79.2
  • Nlipids=128
  • CKCl(zn) = 0.1 M
  • T=303 K
  • simulation time 30 ns
  • force field: CHARMM36
  • please cite DOI: 10.1002/jcc.23262
  • note that CKCl(zn) is the salt conc. in the water layer



Micelles

For micelles, instead of the simulation time, the number of configurations Nconf over which the atomic distribution was averaged is shown. The reason is that micelles sizes can fluctuate over the simulation and the averaging is done over all micelles of the same size found in a simulation [DOI: 10.1016/j.fluid.2016.03.006]. The specific micelle size provided here is chosen according to the criteria discussed in our article [DOI: 10.1016/j.fluid.2016.03.006], where a spherical shape is one criterion. Hence, it is not necessarily the experimental aggregation number. The averaged eccentricity ε as a measure for the spherical shape of the micelle is also provided in the table. Furthermore, the table shows the aggregation number of the micelles Nagg for which the atomic distribution was calculated. In addition, the number of surfactants in the simulation Nsurf is shown. In some cases the number of surfactants in the system is higher than the aggregation number (Nagg ≠ Nsurf), see also [DOI: 10.1016/j.fluid.2016.03.006]. Hence, the system snapshot contains Nsurf surfactants and the atomic distribution is calculated for micelles with Nagg.

If you use averaged atomic distributions for micelles or material from the table below, please consider this article:

Force field parameters for Triton were introduced in this article:

Data for micelles


System Properties Files Comments

Triton X-114/water

  • water/surfactant = 251.8
  • T = 298
  • Nsurf = 40
  • Nagg = 40
  • Nconf = 576
  • ε = 0.12 +/- 0.05
  • force field: CHARMM36

Triton X-100/water

  • water/surfactant = 247.5
  • T = 298
  • Nsurf = 40
  • Nagg = 40
  • Nconf = 346
  • ε = 0.14 +/- 0.04
  • force field: CHARMM36

Brij35/water

  • water/surfactant = 227.3
  • T = 298 K
  • Nsurf = 40
  • Nagg = 39
  • Nconf = 198
  • ε = 0.11 +/- 0.04
  • force field: CHARMM36

C12E10/water

  • water/surfactant = 53.3
  • T = 298 K
  • Nsurf = 60
  • Nagg = 59
  • Nconf = 151
  • ε = 0.11 +/- 0.05
  • force field: CHARMM36

SDS/water

  • water/surfactant = 79.0
  • T = 298
  • Nsurf = 60
  • Nagg = 60
  • Nconf = 359
  • ε = 0.13 +/- 0.04
  • force field: CHARMM36

CTAB/water

  • water/surfactant = 73.7
  • T = 298
  • Nsurf = 90
  • Nagg = 90
  • Nconf = 591
  • ε = 0.15 +/- 0.06
  • force field: CHARMM36