1 |
– |
%\documentclass[aps,pre,twocolumn,amssymb,showpacs,floatfix]{revtex4} |
2 |
– |
%\documentclass[aps,pre,preprint,amssymb]{revtex4} |
1 |
|
\documentclass[12pt]{article} |
2 |
+ |
\usepackage{graphicx} |
3 |
|
\usepackage{times} |
4 |
|
\usepackage{mathptm} |
5 |
|
\usepackage{tabularx} |
6 |
|
\usepackage{setspace} |
7 |
|
\usepackage{amsmath} |
8 |
|
\usepackage{amssymb} |
10 |
– |
\usepackage{graphicx} |
9 |
|
\usepackage[ref]{overcite} |
10 |
|
\pagestyle{plain} |
11 |
|
\pagenumbering{arabic} |
16 |
|
\renewcommand{\baselinestretch}{1.2} |
17 |
|
\renewcommand\citemid{\ } % no comma in optional reference note |
18 |
|
|
19 |
+ |
|
20 |
|
\begin{document} |
21 |
|
%\renewcommand{\thefootnote}{\fnsymbol{footnote}} |
22 |
|
%\renewcommand{\theequation}{\arabic{section}.\arabic{equation}} |
317 |
|
\centering |
318 |
|
\includegraphics[width=4in]{2lipidModel} |
319 |
|
\caption{The parameters defining the behavior of the lipid |
320 |
< |
models. $l / d$ is the ratio of the head group to body diameter. |
320 |
> |
models. $\sigma_h / d$ is the ratio of the head group to body diameter. |
321 |
|
Molecular bodies had a fixed aspect ratio of 3.0. The solvent model |
322 |
|
was a simplified 4-water bead ($\sigma_w \approx d$) that has been |
323 |
< |
used in other coarse-grained (DPD) simulations. The dipolar strength |
323 |
> |
used in other coarse-grained simulations. The dipolar strength |
324 |
|
(and the temperature and pressure) were the only other parameters that |
325 |
|
were varied systematically.\label{fig:lipidModel}} |
326 |
|
\end{figure} |
364 |
|
R_{ij} = \sqrt {r_{ij}^2 + \frac{{d_i^2 }}{4} + \frac{{d_j^2 |
365 |
|
}}{4}}. |
366 |
|
\end{equation} |
367 |
< |
Here, $d_i$ and $d_j$ are effect charge separation distances |
368 |
< |
associated with each of the two dipolar sites. This approximation to |
369 |
< |
the multipole expansion maintains the fast fall-off of the multipole |
370 |
< |
potentials but lacks the normal divergences when two polar groups get |
371 |
< |
close to one another. |
367 |
> |
Here, $d_i$ and $d_j$ are charge separation distances associated with |
368 |
> |
each of the two dipolar sites. This approximation to the multipole |
369 |
> |
expansion maintains the fast fall-off of the multipole potentials but |
370 |
> |
lacks the normal divergences when two polar groups get close to one |
371 |
> |
another. |
372 |
|
|
373 |
|
For the interaction between nonequivalent uniaxial ellipsoids (in this |
374 |
|
case, between spheres and ellipsoids), the spheres are treated as |
387 |
|
replicate the dielectric properties of water. Note that although we |
388 |
|
are using larger cutoff and switching radii than Marrink {\it et al.}, |
389 |
|
our solvent density at 300 K remains at 0.944 g cm$^{-3}$, and the |
390 |
< |
solvent diffuses at 0.43 $\AA^2 ps^{-1}$ (roughly twice as fast as |
391 |
< |
liquid water). |
390 |
> |
solvent diffuses at 0.43 $\AA^2 ps^{-1}$ (only twice as fast as liquid |
391 |
> |
water). |
392 |
|
|
393 |
|
\begin{table*} |
394 |
|
\begin{minipage}{\linewidth} |
421 |
|
The parameters that were systematically varied in this study were the |
422 |
|
size of the head group ($\sigma_h$), the strength of the dipole moment |
423 |
|
($\mu$), and the temperature of the system. Values for $\sigma_h$ |
424 |
< |
ranged from 5.5 \AA\ to 6.5 \AA\ . If the width of the tails is taken |
424 |
> |
ranged from 5.5 \AA\ to 6.5 \AA. If the width of the tails is taken |
425 |
|
to be the unit of length, these head groups correspond to a range from |
426 |
|
$1.2 d$ to $1.41 d$. Since the solvent beads are nearly identical in |
427 |
|
diameter to the tail ellipsoids, all distances that follow will be |
537 |
|
entirely on the size of the head bead relative to the molecular body. |
538 |
|
These values are tabulated in table \ref{tab:property}. Kucera {\it |
539 |
|
et al.} have measured values for the head group spacings for a number |
540 |
< |
of PC lipid bilayers that range from 30.8 \AA (DLPC) to 37.8 (DPPC). |
540 |
> |
of PC lipid bilayers that range from 30.8 \AA\ (DLPC) to 37.8 \AA\ (DPPC). |
541 |
|
They have also measured values for the area per lipid that range from |
542 |
|
60.6 |
543 |
|
\AA$^2$ (DMPC) to 64.2 \AA$^2$ |
787 |
|
We have computed translational diffusion constants for lipid molecules |
788 |
|
from the mean-square displacement, |
789 |
|
\begin{equation} |
790 |
< |
D = \lim_{t\rightarrow \infty} \frac{1}{6 t} \langle {|\left({\bf r}_{i}(t) - {\bf r}_{i}(0) \right)|}^2 \rangle, |
790 |
> |
D = \lim_{t \rightarrow \infty} \frac{1}{6 t} \langle {|\left({\bf r}_{i}(t) - {\bf r}_{i}(0) \right)|}^2 \rangle, |
791 |
|
\end{equation} |
792 |
|
of the lipid bodies. Translational diffusion constants for the |
793 |
|
different head-to-tail size ratios (all at 300 K) are shown in table |
796 |
|
polynomial correlation function, |
797 |
|
\begin{equation} |
798 |
|
C_{\ell}(t) = \langle P_{\ell}\left({\bf \mu}_{i}(t) \cdot {\bf |
799 |
< |
\mu}_{i}(0) \right) |
799 |
> |
\mu}_{i}(0) \right) \rangle |
800 |
|
\end{equation} |
801 |
|
of the head group dipoles. The orientational correlation functions |
802 |
|
appear to have multiple components in their decay: a fast ($12 \pm 2$ |
833 |
|
\caption{Fit values for the rotational correlation times for the head |
834 |
|
groups ($\tau^h$) and molecular bodies ($\tau^b$) as well as the |
835 |
|
translational diffusion constants for the molecule as a function of |
836 |
< |
the head-to-body width ratio (all at 300 K). In all of the phases, |
837 |
< |
the head group correlation functions decay with an fast librational |
838 |
< |
contribution ($12 \pm 1$ ps). There are additional moderate |
839 |
< |
($\tau^h_{\rm mid}$) and slow $\tau^h_{\rm slow}$ contributions to |
840 |
< |
orientational decay that depend strongly on the phase exhibited by the |
841 |
< |
lipids. The symmetric ripple phase ($\sigma_h / d = 1.35$) appears to |
842 |
< |
exhibit the slowest molecular reorientation.} |
836 |
> |
the head-to-body width ratio. All correlation functions and transport |
837 |
> |
coefficients were computed from microcanonical simulations with an |
838 |
> |
average temperture of 300 K. In all of the phases, the head group |
839 |
> |
correlation functions decay with an fast librational contribution ($12 |
840 |
> |
\pm 1$ ps). There are additional moderate ($\tau^h_{\rm mid}$) and |
841 |
> |
slow $\tau^h_{\rm slow}$ contributions to orientational decay that |
842 |
> |
depend strongly on the phase exhibited by the lipids. The symmetric |
843 |
> |
ripple phase ($\sigma_h / d = 1.35$) appears to exhibit the slowest |
844 |
> |
molecular reorientation.} |
845 |
|
\begin{tabular}{lcccc} |
846 |
|
\hline |
847 |
|
$\sigma_h / d$ & $\tau^h_{\rm mid} (ns)$ & $\tau^h_{\rm |
848 |
< |
slow} (\mu s)$ & $\tau_b (\mu s)$ & $D (\times 10^{-11} m^2 s^{-1})$ \\ |
848 |
> |
slow} (\mu s)$ & $\tau^b (\mu s)$ & $D (\times 10^{-11} m^2 s^{-1})$ \\ |
849 |
|
\hline |
850 |
|
1.20 & $0.4$ & $9.6$ & $9.5$ & $0.43(1)$ \\ |
851 |
|
1.28 & $2.0$ & $13.5$ & $3.0$ & $5.91(3)$ \\ |
920 |
|
computable for the all-atom and coarse-grained simulations that have |
921 |
|
been published in the literature.\cite{deVries05} |
922 |
|
|
923 |
+ |
Experimental verification of our predictions of dipolar orientation |
924 |
+ |
correlating with the ripple direction would require knowing both the |
925 |
+ |
local orientation of a rippled region of the membrane (available via |
926 |
+ |
AFM studies of supported bilayers) as well as the local ordering of |
927 |
+ |
the membrane dipoles. Obtaining information about the local |
928 |
+ |
orientations of the membrane dipoles may be available from |
929 |
+ |
fluorescence detected linear dichroism (LD). Benninger {\it et al.} |
930 |
+ |
have recently used axially-specific chromophores |
931 |
+ |
2-(4,4-difluoro-5,7-dimethyl-4-bora-3a,4a-diaza-s-indacene-3-pentanoyl)-1-hexadecanoyl-sn-glycero-3-phospocholine |
932 |
+ |
($\beta$-BODIPY FL C5-HPC or BODIPY-PC) and 3,3' |
933 |
+ |
dioctadecyloxacarbocyanine perchlorate (DiO) in their |
934 |
+ |
fluorescence-detected linear dichroism (LD) studies of plasma |
935 |
+ |
membranes of living cells.\cite{Benninger:2005qy} The DiO dye aligns |
936 |
+ |
its transition moment perpendicular to the membrane normal, while the |
937 |
+ |
BODIPY-PC transition dipole is parallel with the membrane normal. |
938 |
+ |
Without a doubt, using fluorescence detection of linear dichroism in |
939 |
+ |
concert with AFM surface scanning would be difficult experiments to |
940 |
+ |
carry out. However, there is some hope of performing experiments to |
941 |
+ |
either verify or falsify the predictions of our simulations. |
942 |
+ |
|
943 |
|
Although our model is simple, it exhibits some rich and unexpected |
944 |
|
behaviors. It would clearly be a closer approximation to reality if |
945 |
|
we allowed bending motions between the dipoles and the molecular |