28 |
|
hydrophobic interior of the membrane, and for the hydrophobic tails |
29 |
|
to be exposed to the aqueous environment\cite{Sasaki2004}. A number |
30 |
|
of studies indicate that the flipping of phospholipids occurs |
31 |
< |
rapidly in the eukaryotic ER and the bacterial cytoplasmic membrane |
32 |
< |
via a bi-directional, facilitated diffusion process requiring no |
33 |
< |
metabolic energy input. Another system of interest would be the |
34 |
< |
distribution of sites occupied by inhaled anesthetics in membrane. |
35 |
< |
Although the physiological effects of anesthetics have been |
36 |
< |
extensively studied, the controversy over their effects on lipid |
37 |
< |
bilayers still continues. Recent deuterium NMR measurements on |
38 |
< |
halothane in POPC lipid bilayers suggest the anesthetics are |
39 |
< |
primarily located at the hydrocarbon chain region\cite{Baber1995}. |
40 |
< |
Infrared spectroscopy experiments suggest that halothane in DMPC |
41 |
< |
lipid bilayers lives near the membrane/water |
42 |
< |
interface\cite{Lieb1982}. |
31 |
> |
rapidly in the eukaryotic endoplasmic reticulum and the bacterial |
32 |
> |
cytoplasmic membrane via a bi-directional, facilitated diffusion |
33 |
> |
process requiring no metabolic energy input. Another system of |
34 |
> |
interest is the distribution of sites occupied by inhaled |
35 |
> |
anesthetics in membrane. Although the physiological effects of |
36 |
> |
anesthetics have been extensively studied, the controversy over |
37 |
> |
their effects on lipid bilayers still continues. Recent deuterium |
38 |
> |
NMR measurements on halothane on POPC lipid bilayers suggest the |
39 |
> |
anesthetics are primarily located at the hydrocarbon chain |
40 |
> |
region\cite{Baber1995}. However, infrared spectroscopy experiments |
41 |
> |
suggest that halothane in DMPC lipid bilayers lives near the |
42 |
> |
membrane/water interface\cite{Lieb1982}. |
43 |
|
|
44 |
|
Molecular dynamics simulations have proven to be a powerful tool for |
45 |
|
studying the functions of biological systems, providing structural, |
95 |
|
w'(r_{ij} ,\Omega _i ,\Omega _j ) = (\cos \theta _{ij} - 0.6)^2 (\cos \theta _{ij} + 0.8)^2 - w_0 \\ |
96 |
|
\end{array} |
97 |
|
\] |
98 |
< |
Although dipole-dipole and sticky interactions are more |
98 |
> |
Although the dipole-dipole and sticky interactions are more |
99 |
|
mathematically complicated than Coulomb interactions, the number of |
100 |
|
pair interactions is reduced dramatically both because the model |
101 |
|
only contains a single-point as well as "short range" nature of the |
102 |
< |
higher order interaction. |
102 |
> |
more expensive interaction. |
103 |
|
|
104 |
|
\subsection{\label{lipidSection:coarseGrained}The Coarse-Grained Phospholipid Model} |
105 |
|
|
117 |
|
\caption[A representation of coarse-grained phospholipid model]{A |
118 |
|
representation of coarse-grained phospholipid model. The lipid |
119 |
|
headgroup consists of $\text{{\sc PO}}_4$ group (yellow), |
120 |
< |
$\text{{\sc NC}}_4$ group (blue) and a united $\text{{\sc C}} atom |
120 |
> |
$\text{{\sc NC}}_4$ group (blue) and a united $\text{{\sc C}}$ atom |
121 |
|
(gray) $ with a dipole, while the glycerol backbone includes dipolar |
122 |
< |
$\text{{\sc CE}}$ (read) and $\text{{\sc CK}}$ (pink) groups. Alkyl |
122 |
> |
$\text{{\sc CE}}$ (red) and $\text{{\sc CK}}$ (pink) groups. Alkyl |
123 |
|
groups in hydrocarbon chains are simply represented by gray united |
124 |
|
atoms.} \label{lipidFigure:coarseGrained} |
125 |
|
\end{figure} |
135 |
|
\] |
136 |
|
where $N_A$ and $N_B$ are the number of point charges in the two |
137 |
|
molecular species. Originally developed to study ionic crystals, the |
138 |
< |
Ewald summation method mathematically transforms this |
139 |
< |
straightforward but conditionally convergent summation into two more |
140 |
< |
complicated but rapidly convergent sums. One summation is carried |
141 |
< |
out in reciprocal space while the other is carried out in real |
142 |
< |
space. An alternative approach is a multipole expansion, which is |
143 |
< |
based on electrostatic moments, such as charge (monopole), dipole, |
144 |
< |
quadruple etc. |
138 |
> |
Ewald sum method mathematically transforms this straightforward but |
139 |
> |
conditionally convergent summation into two more complicated but |
140 |
> |
rapidly convergent sums. One summation is carried out in reciprocal |
141 |
> |
space while the other is carried out in real space. An alternative |
142 |
> |
approach is the multipole expansion, which is based on electrostatic |
143 |
> |
moments, such as charge (monopole), dipole, quadrupole etc. |
144 |
|
|
145 |
|
Here we consider a linear molecule which consists of two point |
146 |
|
charges $\pm q$ located at $ ( \pm \frac{d}{2},0)$. The |
218 |
|
|
219 |
|
\subsection{One Lipid in Sea of Water Molecules} |
220 |
|
|
221 |
< |
To exclude the inter-headgroup interaction, atomistic models of one |
222 |
< |
lipid (DMPC or DLPE) in sea of water molecules (TIP3P) were built |
223 |
< |
and studied using atomistic molecular dynamics. The simulation was |
224 |
< |
analyzed using a set of radial distribution functions, which give |
225 |
< |
the probability of finding a pair of molecular species a distance |
226 |
< |
apart, relative to the probability expected for a completely random |
227 |
< |
distribution function at the same density. |
221 |
> |
To tune our parameters without the inter-headgroup interactions, |
222 |
> |
atomistic models of one lipid (DMPC or DLPE) in sea of water |
223 |
> |
molecules (TIP3P) were built and studied using atomistic molecular |
224 |
> |
dynamics. The simulation was analyzed using a set of radial |
225 |
> |
distribution functions, which give the probability of finding a pair |
226 |
> |
of molecular species a distance apart, relative to the probability |
227 |
> |
expected for a completely random distribution function at the same |
228 |
> |
density. |
229 |
|
|
230 |
|
\begin{equation} |
231 |
|
g_{AB} (r) = \frac{1}{{\rho _B }}\frac{1}{{N_A }} < \sum\limits_{i |
247 |
|
in DMPC is around the $\text{{\sc PO}}_4$ Group. When the water |
248 |
|
molecules are far from the headgroup, the distribution of the two |
249 |
|
angles should be uniform. The correlation close to center of the |
250 |
< |
headgroup dipole in both plots tell us that in the closely-bound |
250 |
> |
headgroup dipole in both plots tells us that in the closely-bound |
251 |
|
region, the dipoles of the water molecules are preferentially |
252 |
|
anti-aligned with the dipole of headgroup. When the water molecules |
253 |
|
are far from the headgroup, the distribution of the two angles |
326 |
|
period of 2 ns. During the equilibration period, the system was |
327 |
|
initially simulated at constant volume for 1 ns. Once the system was |
328 |
|
equilibrated at constant volume, the cell dimensions of the system |
329 |
< |
was relaxed by performing under NPT conditions using Nos¨¦-Hoover |
329 |
> |
was relaxed by performing under NPT conditions using Nos\'{e}-Hoover |
330 |
|
extended system isothermal-isobaric dynamics. After equilibration, |
331 |
|
different properties were evaluated over a production run of 5 ns. |
332 |
|
|
424 |
|
|
425 |
|
\section{\label{lipidSection:Conclusion}Conclusion} |
426 |
|
|
427 |
< |
Atomistic simulations are used in this study to determine the |
427 |
> |
Atomistic simulations have been used in this study to determine the |
428 |
|
preferred orientation and location of water molecules relative to |
429 |
|
the location and orientation of the PC and PE lipid headgroups. |
430 |
< |
Based on the result from all-atom simulations, we developed a simple |
431 |
< |
coarse-grained model capturing essential features of the headgroup |
432 |
< |
solvation which is crucial to transport process in membrane system. |
433 |
< |
In addition, the model has been explored in a bilayer system which |
434 |
< |
is shown to have reasonable electron density profile, |
435 |
< |
$\text{S}_{\text{{\sc cd}}}$ order parameter and other structural |
436 |
< |
properties. The accuracy of this model is achieved by matching |
437 |
< |
atomistic result. It is also easy to represent different |
438 |
< |
phosphorlipids by changing a few parameters of the model. Another |
430 |
> |
Based on the results from our all-atom simulations, we developed a |
431 |
> |
simple coarse-grained model which captures the essential features of |
432 |
> |
the headgroup solvation which is crucial to transport process in |
433 |
> |
membrane system. In addition, the model has been explored in a |
434 |
> |
bilayer system was shown to have reasonable electron density |
435 |
> |
profiles, $\text{S}_{\text{{\sc cd}}}$ order parameter and other |
436 |
> |
structural properties. The accuracy of this model is achieved by |
437 |
> |
matching atomistic result. It is also easy to represent different |
438 |
> |
phospholipids by changing a few parameters of the model. Another |
439 |
|
important characteristic of this model distinguishing itself from |
440 |
|
other models\cite{Goetz1998,Marrink2004} is the computational speed |
441 |
< |
gaining by introducing short range electrostatic approximation. |
441 |
> |
gained by introducing a short range electrostatic approximation. |
442 |
|
Further studies of this system using z-constraint method could |
443 |
|
extend the time length of the simulations to study transport |
444 |
|
phenomena in large-scale membrane systems. |