1 |
mmeineke |
226 |
\documentclass[11pt]{article} |
2 |
|
|
|
3 |
|
|
\usepackage{graphicx} |
4 |
|
|
\usepackage{color} |
5 |
|
|
\usepackage{floatflt} |
6 |
|
|
\usepackage{amsmath} |
7 |
|
|
\usepackage{amssymb} |
8 |
|
|
\usepackage{subfigure} |
9 |
|
|
\usepackage{palatino} |
10 |
|
|
\usepackage[ref]{overcite} |
11 |
|
|
|
12 |
|
|
|
13 |
|
|
|
14 |
|
|
\pagestyle{plain} |
15 |
|
|
\pagenumbering{arabic} |
16 |
|
|
\oddsidemargin 0.0cm \evensidemargin 0.0cm |
17 |
|
|
\topmargin -21pt \headsep 10pt |
18 |
|
|
\textheight 9.0in \textwidth 6.5in |
19 |
|
|
\brokenpenalty=10000 |
20 |
|
|
\renewcommand{\baselinestretch}{1.2} |
21 |
|
|
\renewcommand\citemid{\ } % no comma in optional reference note |
22 |
|
|
|
23 |
|
|
|
24 |
|
|
\begin{document} |
25 |
|
|
|
26 |
|
|
|
27 |
|
|
\title{A Mesoscale Model for Phospholipid Simulations} |
28 |
|
|
|
29 |
|
|
\author{Matthew A. Meineke, Charles F. Vardeman II, and J. Daniel Gezelter\\ |
30 |
|
|
Department of Chemistry and Biochemistry\\ |
31 |
|
|
University of Notre Dame\\ |
32 |
|
|
Notre Dame, Indiana 46556} |
33 |
|
|
|
34 |
|
|
\date{\today} |
35 |
|
|
\maketitle |
36 |
|
|
|
37 |
|
|
\section{Model and Methodology} |
38 |
|
|
|
39 |
|
|
\subsection{The Phospholipid Model} |
40 |
|
|
\label{sec:lipidModel} |
41 |
|
|
|
42 |
|
|
\begin{figure} |
43 |
|
|
\centering |
44 |
|
|
\includegraphics[angle=-90,width=80mm]{lipidModel.epsi} |
45 |
|
|
\caption{A representation of the lipid model. $\phi$ is the torsion angle, $\theta$ is the bend angle, $\mu$ is the dipole moment of the head group, and n is the chain length.} |
46 |
|
|
\label{fig:lipidModel} |
47 |
|
|
\end{figure} |
48 |
|
|
|
49 |
|
|
The lipid molecules in our simulations are unified atom models. Figure |
50 |
|
|
\ref{fig:lipidModel} shows a schematic for one of our |
51 |
|
|
lipids. The head group of the phospholipid is replaced by a single |
52 |
|
|
Lennard-Jones sphere with a freely oriented dipole placed at it's |
53 |
|
|
center. The magnitude of the dipole moment is 20.6 D, chosen to match |
54 |
|
|
that of DPPC\cite{Cevc87}. The tail atoms are unified $\text{CH}_2$ |
55 |
|
|
and $\text{CH}_3$ atoms and are also modeled as Lennard-Jones |
56 |
|
|
spheres. The total potential for the lipid is represented by Equation |
57 |
|
|
\ref{eq:lipidModelPot}. |
58 |
|
|
|
59 |
|
|
\begin{equation} |
60 |
|
|
V_{\text{lipid}} = |
61 |
|
|
\sum_{i}V_{i}^{\text{internal}} |
62 |
|
|
+ \sum_i \sum_{j>i} \sum_{\alpha_i} |
63 |
|
|
\sum_{\beta_j} |
64 |
|
|
V_{\text{LJ}}(r_{\alpha_{i}\beta_{j}}) |
65 |
|
|
+\sum_i\sum_{j>i}V_{\text{dp}}(r_{1_i,1_j},\Omega_{1_i},\Omega_{1_j}) |
66 |
|
|
\label{eq:lipidModelPot} |
67 |
|
|
\end{equation} |
68 |
|
|
where, |
69 |
|
|
\begin{equation} |
70 |
|
|
V_{i}^{\text{internal}} = |
71 |
|
|
\sum_{\text{bends}}V_{\text{bend}}(\theta_{\alpha\beta\gamma}) |
72 |
|
|
+ \sum_{\text{torsions}}V_{\text{tors.}}(\phi_{\alpha\beta\gamma\zeta}) |
73 |
|
|
+ \sum_{\alpha_i} \sum_{\beta_i > (\alpha_i + 4)}V_{\text{LJ}} |
74 |
|
|
(r_{\alpha_i \beta_i}) |
75 |
|
|
\label{eq:lipidModelPotInternal} |
76 |
|
|
\end{equation} |
77 |
|
|
|
78 |
|
|
The non-bonded interactions, $V_{\text{LJ}}$ and $V_{\text{dp}}$, are |
79 |
|
|
the Lennard-Jones and dipole-dipole interactions respectively. For the |
80 |
|
|
bonded potentials, only the bend and the torsional potentials are |
81 |
|
|
calculated. The bond potential is not calculated, and the bond lengths |
82 |
|
|
are constrained via RATTLE.\cite{leach01:mm} The bend potential is of |
83 |
|
|
the form: |
84 |
|
|
\begin{equation} |
85 |
|
|
V_{\text{bend}}(\theta_{\alpha\beta\gamma}) |
86 |
|
|
= k_{\theta}\frac{(\theta_{\alpha\beta\gamma} - \theta_0)^2}{2} |
87 |
|
|
\label{eq:bendPot} |
88 |
|
|
\end{equation} |
89 |
|
|
Where $k_{\theta}$ sets the stiffness of the bend potential, and $\theta_0$ |
90 |
|
|
sets the equilibrium bend angle. The torsion potential was given by: |
91 |
|
|
\begin{equation} |
92 |
|
|
V_{\text{tors.}}(\phi_{\alpha\beta\gamma\zeta}) |
93 |
|
|
= c_1 [1+\cos\phi_{\alpha\beta\gamma\zeta}] |
94 |
|
|
+ c_2 [1 - \cos(2\phi_{\alpha\beta\gamma\zeta})] |
95 |
|
|
+ c_3 [1 + \cos(3\phi_{\alpha\beta\gamma\zeta})] |
96 |
|
|
\label{eq:torsPot} |
97 |
|
|
\end{equation} |
98 |
|
|
All parameters for bonded and non-bonded potentials in the tail atoms |
99 |
|
|
were taken from TraPPE.\cite{Siepmann1998} The bonded interactions for |
100 |
|
|
the head atom were also taken from TraPPE, however it's dipole moment |
101 |
|
|
and mass were based on the properties of the phosphatidylcholine head |
102 |
|
|
group. The Lennard-Jones parameter for the head group was chosen such |
103 |
|
|
that it was roughly twice the size of a $\text{CH}_3$ atom, and it's |
104 |
|
|
well depth was set to be approximately equal to that of $\text{CH}_3$. |
105 |
|
|
|
106 |
|
|
|
107 |
|
|
\end{document} |
108 |
|
|
|