Protein/Membrane system size
This page collects notes on how to decide how to initially setup a membrane protein simulation. Please do not follow these notes blindly — they are more guidelines than actual rules.
Once you've decided on box dimensions and number of lipids you can use any of the Transmembrane protein insertion methods methods to set up the system.
- Decide on the bilayer composition. Have you got force field parameters for the lipids?
- System size:
- rule of thumb 1: have at least 2-3 lipid layers between protein and periodic box boundary in the x-y plane to ensure(?) membrane-like behaviour; more may be necessary: you will need to monitor the hydrophobic mismatch and make the system big enough so that the bilayer at the box edges has its natural thickness.
- rule of thumb 2: have at least 1-1.5nm of water between protein and z-boundary (so that Coulomb forces are effectively screened)
- Especially for coarse-grained self-assembly simulations you also need to take into account the lipid-to-water ratio and its affect on membrane structure.
- Add ions at e.g. 100 mM (+counter ions) for additional Coulomb screening
- Remember that run time scales like N log N (or even N2) and N scales with the volume, i.e. L3! Thus: keep your system small and pretty (i.e. the smallest system that still behaves like a big system... invariably you are trading size-artifacts versus speed and thus sampling)
Number of lipids
Estimate the cross sectional area of the protein, e.g. using the radius of gyration
Aprot = π RG
and with the approximate area per lipid Al = 64 Å the number of lipids in two leaflets is for a box with dimensions Lx and Ly
N = 2(LxLy - Aprot) / Al
If you are simulating a mixed bilayer, e.g. POPE:POPG 4:1 then you will need 4N/5 POPE lipids and N/5 POPG ones.
Number of ions
Once the system is solvated you can calculate the numbers of ions <math>N</math> needed for a given concentration c (in addition to any counter ions):
c = n/V = N/(NAV)
(where NA = 6.02214179 × 1023mol-1 is the Avogadro constant).
The problem is to get the water volume of the inhomogenous system, consisting of membrane and protein with the water. An approximation is to use the volume of a water molecule at standard conditions, vw = 30 Å3 and calculate V = Nw vw from the number of waters in the system, Nw. The equations become even simpler if we use the standard concentration of water, cw = nw/(Nw vw) = 55.5 mol/l.
Using all this we get
N+ = c+/V = c+/cw Nw
for the (monovalent) cations, and correspondingly for monovalent anions
N- = N+
For a concentration of 100 mM = 0.1 M and 8000 water molecules we will need to add
N- = N+ = 0.1 M/55.5 M × 8000 = 14
ions of each kind (e.g. 14 sodium and 14 chloride and whatever counterions are required to make the system charge neutral).