Biomedical subjects

Ilpo Vattulainen

Publications and source records attributed to Ilpo Vattulainen.

At least 19 recordsLinked to original sources

PubMed2007

Ion leakage through transient water pores in protein-free lipid membranes driven by transmembrane ionic charge imbalance.

We have employed atomic-scale molecular dynamics simulations to address ion leakage through transient water pores in protein-free phospholipid membranes. Our results for phospholipid membranes in aqueous solution with NaCl and KCl salts show that the formation of transient water pores and the consequent ion leakage can be induced and be driven by a transmembrane ionic charge imbalance, an inherent feature in living cells. These processes take place if the gradient is large enough to develop a sufficiently significant potential difference across the membrane. The transport of cations and anions through the water pores is then seen; it discharges the transmembrane potential, considerably reduces the size of a water pore, and makes the water pore metastable, leading eventually to its sealing. The ion transport is found to be sensitive to the type of ions. It turns out that Na(+) and Cl(-) ions leak through a membrane at approximately the same ratio despite the fact that Na(+) ions are expected to experience a lower potential barrier for the permeation through the pore. This is because of strong interactions of sodium ions with the carbonyl region of a phospholipid membrane as well as with lipid headgroups forming pore "walls," considerably slowing down the permeation of sodium ions. In contrast, we observed a pronounced selectivity of a phospholipid membrane to the permeation of potassium ions as compared to chloride ions: Potassium ions, being larger than sodium ions, interact only weakly with phospholipid headgroups, so that these interactions are not able to compensate for a large difference in free-energy barriers for permeation of K(+) and Cl(-) ions. These findings are found to be robust to a choice of force-field parameters for ions (tested by Gromacs and Charmm force-fields for ions). What is more, a potassium ion is found to be able to permeate a membrane along an alternate, "water-defect-mediated" pathway without actual formation of a pore. The "water-defect-mediated" leakage involves formation of a single water defect only and is found to be at least one order of magnitude faster than the pore-mediated ion leakage.

Cell Membrane↗

PubMed2006

Tilt: major factor in sterols' ordering capability in membranes.

Using extensive atomistic simulations, we show that there is a single experimentally accessible parameter--the sterol tilt--that can be used to determine a sterol's capability to induce order, and thus to promote, e.g., formation of lipid rafts. The observations also facilitate designing new effective sterols.

Membranes, Artificial↗

PubMed2006

Atomic-scale structure and electrostatics of anionic palmitoyloleoylphosphatidylglycerol lipid bilayers with Na+ counterions.

Anionic palmitoyloleoylphosphatidylglycerol (POPG) is one of the most abundant lipids in nature, yet its atomic-scale properties have not received significant attention. Here we report extensive 150-ns molecular dynamics simulations of a pure POPG lipid membrane with sodium counterions. It turns out that the average area per lipid of the POPG bilayer under physiological conditions is approximately 19% smaller than that of a bilayer built from its zwitterionic phosphatidylcholine analog, palmitoyloleoylphosphatidylcholine. This suggests that there are strong attractive interactions between anionic POPG lipids, which overcome the electrostatic repulsion between negative charges of PG headgroups. We demonstrate that interlipid counterion bridges and strong intra- and intermolecular hydrogen bonding play a key role in this seemingly counterintuitive behavior. In particular, the substantial strength and stability of ion-mediated binding between anionic lipid headgroups leads to complexation of PG molecules and ions and formation of large PG-ion clusters that act in a concerted manner. The ion-mediated binding seems to provide a possible molecular-level explanation for the low permeability of PG-containing bacterial membranes to organic solvents: highly polar interactions at the water/membrane interface are able to create a high free energy barrier for hydrophobic molecules such as benzene.

Anions↗

PubMed2006

Insight into the putative specific interactions between cholesterol, sphingomyelin, and palmitoyl-oleoyl phosphatidylcholine.

The effects of cholesterol (Chol) on phospholipid bilayers include ordering of the fatty acyl chains, condensing of the lipids in the bilayer plane, and promotion of the liquid-ordered phase. These effects depend on the type of phospholipids in the bilayer and are determined by the nature of the underlying molecular interactions. As for Chol, it has been shown to interact more favorably with sphingomyelin than with most phosphatidylcholines, which in given circumstances leads to formation of lateral domains. However, the exact origin and nature of Chol-phospholipid interactions have recently been subjects of speculation. We examine interactions between Chol, palmitoylsphingomyelin (PSM) and palmitoyl-oleoyl-phosphatidylcholine (POPC) in hydrated lipid bilayers by extensive atom-scale molecular dynamics simulations. We employ a tailored lipid configuration: Individual PSM and Chol monomers, as well as PSM-Chol dimers, are embedded in a POPC lipid bilayer in the liquid crystalline phase. Such a setup allows direct comparison of dimeric and monomeric PSMs and Chol, which ultimately shows how the small differences in PSM and POPC structure can lead to profoundly different interactions with Chol. Our analysis shows that direct hydrogen bonding between PSM and Chol does not provide an adequate explanation for their putative specific interaction. Rather, a combination of charge-pairing, hydrophobic, and van der Waals interactions leads to a lower tilt in PSM neighboring Chol than in Chol with only POPC neighbors. This implies improved Chol-induced ordering of PSM's chains over POPC's chains. These findings are discussed in the context of the hydrophobic mismatch concept suggested recently.

Cholesterol↗

PubMed2006

Influence of pyrene-labeling on fluid lipid membranes.

We elucidate the influence of pyrene-labeled phospholipids on the structural properties of a fluid dipalmitoylphosphatidylcholine lipid membrane. To this end, we employ extensive atomic-scale molecular dynamics simulations with varying concentrations of pyrene-linked lipids. We find pyrene labeling to perturb the membrane structure significantly in the vicinity of the probe, the correlation length in the bilayer plane being about 1.0-1.5 nm. The local perturbations lead to enhanced ordering and packing of lipid acyl chains located in the vicinity of the probe. Surprisingly, this holds true not only for lipids that reside in the same leaflet as the pyrene-labeled probe but also for lipids in the opposite monolayer. The latter is due to substantial interdigitation of the pyrene moiety into the opposite leaflet, suggesting that occasional excimer formation may take place for probes in different leaflets. As a related issue, we also discuss the location and conformational orientation of the pyrene moieties. In particular, the orientational distribution of pyrene turns out to be more broad and diverse than the distribution of the corresponding acyl tails of nonlabeled lipids.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2006

Interaction of fusidic acid with lipid membranes: Implications to the mechanism of antibiotic activity.

We have studied the effects of cholesterol and steroid-based antibiotic fusidic acid (FA) on the behavior of lipid bilayers using a variety of experimental techniques together with atomic-scale molecular dynamics simulations. Capillary electrophoretic measurements showed that FA was incorporated into fluid 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine membranes. Differential scanning calorimetry in turn showed that FA only slightly altered the thermodynamic properties of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers, whereas cholesterol abolished all endotherms when the mole fraction of cholesterol (X(chol)) was >0.20. Fluorescence spectroscopy was then used to further characterize the influence of these two steroids on DPPC large unilamellar vesicles. In the case of FA, our result strongly suggested that FA was organized into lateral microdomains with increased water penetration into the membrane. For cholesterol/DPPC mixtures, fluorescence spectroscopy results were compatible with the formation of the liquid-ordered phase. A comparison of FA and cholesterol-induced effects on DPPC bilayers through atomistic molecular dynamics simulations showed that both FA and cholesterol tend to order neighboring lipid chains. However, the ordering effect of FA was slightly weaker than that of cholesterol, and especially for deprotonated FA the difference was significant. Summarizing, our results show that FA is readily incorporated into the lipid bilayer where it is likely to be enriched into lateral microdomains. These domains could facilitate the association of elongation factor-G into lipid rafts in living bacteria, enhancing markedly the antibiotic efficacy of FA.

Anti-Bacterial Agents↗

PubMed2006

Molecular dynamics study of charged dendrimers in salt-free solution: effect of counterions.

Polyamidoamine dendrimers, being protonated under physiological conditions, represent a promising class of nonviral, nanosized vectors for drug and gene delivery. We performed extensive molecular dynamics simulations of a generic model dendrimer in a salt-free solution with dendrimer's terminal beads positively charged. Solvent molecules as well as counterions were explicitly included as interacting beads. We find that the size of the charged dendrimer depends nonmonotonically on the strength of electrostatic interactions demonstrating a maximum when the Bjerrum length equals the diameter of a bead. Many other structural and dynamic characteristics of charged dendrimers are also found to follow this pattern. We address such a behavior to the interplay between repulsive interactions of the charged terminal beads and their attractive interactions with oppositely charged counterions. The former favors swelling at small Bjerrum lengths and the latter promotes counterion condensation. Thus, counterions can have a dramatic effect on the structure and dynamics of charged dendrimers and, under certain conditions, cannot be treated implicitly.

Computer Simulation↗

PubMed2006

Atomistic simulation studies of cholesteryl oleates: model for the core of lipoprotein particles.

We have conducted molecular dynamics simulations to gain insight into the atomic-scale properties of an isotropic system of cholesteryl oleate (CO) molecules. Cholesteryl esters are major constituents of low density lipoprotein particles, the key players in the formation of atherosclerosis, as well as the storage form of cholesterol. Here the aim is to clarify structural and dynamical properties of CO molecules under conditions, which are suggestive of those in the core of low density lipoprotein particles. The simulations in the fluid phase indicate that the system of CO molecules is characterized by an absence of translational order, as expected, while the orientational order between distinct CO molecules is significant at short distances, persisting over a molecular size. As for intramolecular properties, the bonds along the oleate chain are observed to be weakly ordered with respect to the sterol structure, unlike the bonds along the short hydrocarbon chain of cholesterol where the ordering is significant. The orientational distribution of the oleate chain as a whole with respect to the sterol moiety is of broad nature, having a major amount of extended and a less considerable proportion of bended structures. Distinct transient peaks at specific angles also appear. The diffusion of CO molecules is found to be a slow process and characterized by a diffusion coefficient of the order of 2x10(-9) cm2/s. This is considerably slower than diffusion, e.g., in ordered domains of lipid membranes rich in sphingomyelin and cholesterol. Analysis of the rotational diffusion rates and trans-to-gauche transition rates yield results consistent with experiments.

Atherosclerosis↗

PubMed2005

Pore formation coupled to ion transport through lipid membranes as induced by transmembrane ionic charge imbalance: atomistic molecular dynamics study.

We have employed atomic-scale molecular dynamics simulations to address ion transport through transient water pores in phospholipid membranes. The formation of a water pore is induced by a transmembrane ionic charge imbalance, which gives rise to a significant potential difference across the membrane. The subsequent transport of ions through the pore discharges the transmembrane potential and makes the water pore metastable, leading eventually to its sealing. The findings highlight the importance of ionic charge fluctuations in spontaneous pore formation and their role in ion leakage through protein-free lipid membranes.

Computer Simulation↗

PubMed2005

Under the influence of alcohol: the effect of ethanol and methanol on lipid bilayers.

Extensive microscopic molecular dynamics simulations have been performed to study the effects of short-chain alcohols, methanol and ethanol, on two different fully hydrated lipid bilayer systems (POPC and DPPC) in the fluid phase at 323 K. It is found that ethanol has a stronger effect on the structural properties of the membranes. In particular, the bilayers become more fluid and permeable: ethanol molecules are able to penetrate through the membrane in typical timescales of approximately 200 ns, whereas for methanol that timescale is considerably longer, at least of the order of microseconds. A closer examination exposes a number of effects due to ethanol. Hydrogen-bonding analysis reveals that a large fraction of ethanols is involved in hydrogen bonds with lipids. This in turn is intimately coupled to the ordering of hydrocarbon chains: we find that binding to an ethanol decreases the order of the chains. We have also determined the dependence of lipid-chain ordering on ethanol concentration and found that to be nonmonotonous. Overall, we find good agreement with NMR and micropipette studies.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Influence of chain length and unsaturation on sphingomyelin bilayers.

Sphingomyelins (SMs) are among the most common phospholipid components of plasma membranes, usually constituting a mixture of several molecular species with various fatty acyl chain moieties. In this work, we utilize atomistic molecular dynamics simulations to study the differences in structural and dynamical properties of bilayers comprised of the most common natural SM species. Keeping the sphingosine moiety unchanged, we vary the amide bonded acyl chain from 16 to 24 carbons in length and examine the effect of unsaturation by comparing lipids with saturated and monounsaturated chains. As for structural properties, we find a slight decrease in average area per lipid and a clear linear increase in bilayer thickness with increasing acyl chain length both in saturated and unsaturated systems. Increasing the acyl chain length is found to further the interdigitation across the bilayer center. This is related to the dynamics of SM molecules, as the lateral diffusion rates decrease slightly for an increasing acyl chain length. Interdigitation also plays a role in interleaflet friction, which is stronger for unsaturated chains. The effect of the cis double bond is most significant on the local order parameters and rotation rates of the chains, though unsaturation shows global effects on overall lipid packing and dynamics as well. Regarding hydrogen bonding or properties related to the lipid/water interface region, no significant effects were observed due to varying chain length or unsaturation. The significance of the findings presented is discussed.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Effect of monovalent salt on cationic lipid membranes as revealed by molecular dynamics simulations.

An atomic-scale understanding of cationic lipid membranes is required for development of gene delivery agents based on cationic liposomes. To address this problem, we recently performed molecular dynamics (MD) simulations of mixed lipid membranes comprised of cationic dimyristoyltrimethylammonium propane (DMTAP) and zwitterionic dimyristoylphosphatidylcholine (DMPC) (Biophys. J. 2004, 86, 3461-3472). Given that salt ions are always present under physiological conditions, here we focus on the effects of monovalent salt (NaCl) on cationic (DMPC/DMTAP) membranes. Using atomistic MD simulations, we found that salt-induced changes in membranes depend strongly on their composition. When the DMTAP mole fraction is small (around 6%), the addition of monovalent salt leads to a considerable compression of the membrane and to a concurrent enhancement of the ordering of lipid acyl chains. That is accompanied by reorientation of phosphatidylcholine headgroups in the outward normal direction and slight changes in electrostatic properties. We attribute these changes to complexation of DMPC lipids with Na(+) ions which penetrate deep into the membrane and bind to the carbonyl region of the DMPC lipids. In contrast, at medium and high molar fractions of cationic DMTAP (50 and 75%) a substantial positive surface charge density of the membranes prevents the binding of Na(+) ions, making such membranes almost insensitive to monovalent salt. Finally, we compare our results to the Poisson-Boltzmann theory. With the exception of the immediate vicinity of the bilayer plane, we found excellent agreement with the theory. This is as expected since unlike in the theoretical description the surface is now structured due to its atomic scale nature.

Computer Simulation↗

PubMed2005

Significance of sterol structural specificity. Desmosterol cannot replace cholesterol in lipid rafts.

Desmosterol is an immediate precursor of cholesterol in the Bloch pathway of sterol synthesis and an abundant membrane lipid in specific cell types. The significance of the difference between the two sterols, an additional double bond at position C24 in the tail of desmosterol, is not known. Here, we provide evidence that the biophysical and functional characteristics of the two sterols differ and that this is because the double bond at C24 significantly weakens the sterol ordering potential. In model membranes, desmosterol was significantly weaker than cholesterol in promoting the formation or stability of ordered domains, and in mammalian cell membranes, desmosterol associated less avidly than cholesterol with detergent-resistant membranes. Atomic scale molecular dynamics simulations showed that the double bond gives rise to additional stress in the tail, creating a rigid structure between C24 and C27 and favoring tilting of desmosterol distinct from cholesterol. Functional effects of desmosterol in cell membranes were assessed upon acutely exchanging approximately 70% of cholesterol to desmosterol. This led to impaired raft-dependent signaling via the insulin receptor, whereas non-raft-dependent protein secretion was not affected. We suggest that the choice of cholesterol synthesis route may provide a physiological mechanism to modulate raft-dependent functions in cells.

Albumins↗

PubMed2005

Response to comment by Almeida et al.: free area theories for lipid bilayers--predictive or not?

Free area theories for lateral diffusion in lipid bilayers are reviewed and discussed. It has been suggested by Almeida et al. that free area theories yield quantitative predictions for lateral diffusion coefficients of lipids. We investigate the plausibility of this suggestion by first sketching what is to be expected of a quantitative theory with predictive power, and subsequently examining whether existing free area theories comply with these expectations. Our conclusion is that current free area theories for lipid bilayers are not quantitative theories with predictive power. They involve a number of adjustable parameters, all of which are not estimated independently, but derived from fitting the theory to the very data whose behavior the theory is supposed to predict. Further, the interpretation and behavior of some of the parameters are ambiguous. The best example is the so-called activation barrier, whose qualitative behavior with the cholesterol concentration in a DMPC bilayer varies depending on the experimental method used to generate the input data and the exact assumptions made to formulate the theory. Independent determination of the activation barrier from numerical simulations or experiments appears to be very difficult.

Animals↗

PubMed2005

DNA overstretching transition: ionic strength effects.

As double-stranded DNA is stretched to its B-form contour length, models of polymer elasticity can describe the dramatic increase in measured force. When the molecule is stretched beyond the contour length, it further shows a highly cooperative overstretching transition. We have developed a theoretical description for this transition by coupling the two-state model and the elasticity theory proposed earlier by others. Furthermore, we have extended this model to account for monovalent salt effects on elastic moduli during the transition. We find that this theoretical description is in very good agreement with recent measurements for the salt dependence of the overstretching transition, allowing us to gain insight into the mechanisms that govern the transition. In double-stranded DNA, the effective length per unit charge varies with salt in agreement with the Manning and Poisson-Boltzmann models for thin polyelectrolyte rods, whereas the other model parameters describing structural features have barely any salt dependence. The results thus suggest that the electrostatic component of force-induced overstretching is mediated mesoscopically via elasticity.

Computer Simulation↗

PubMed2005

Influence of DPH on the structure and dynamics of a DPPC bilayer.

We have conducted extensive molecular dynamics (MD) simulations together with differential scanning calorimetry (DSC) and nuclear magnetic resonance (NMR) experiments to quantify the influence of free 1,6-diphenyl-1,3,5-hexatriene (DPH) fluorescent probes on the structure and dynamics of a dipalmitoylphosphatidylcholine bilayer. Atomistic MD simulations show that in the membrane-water interface the influence of DPH is minor, whereas in the acyl-chain region DPH gives rise to major perturbations. In the latter case, DPH is found to influence a wide range of membrane properties, such as the packing and ordering of hydrocarbon tails and the lateral diffusion of lipid molecules. The effects are prominent but of local nature, i.e., the changes observed in the properties of lipid molecules are significant in the vicinity of DPH, but reduce rapidly as the distance from the probe increases. Long-range perturbations due to DPH are hence not expected. Detailed DSC and (2)H NMR measurements support this view. DSC shows only subtle perturbation to the cooperative behavior of the membrane system in the presence of DPH, and (2)H NMR shows that DPH gives rise to a slight increase in the lipid chain order, in agreement with MD simulations. Potential effects of other probes such as pyrene are briefly discussed.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Exploring the effect of xenon on biomembranes.

We report the initial findings of 100 ns molecular dynamics simulations of the role of cellular membranes in general anaesthesia. The effect of xenon on hydrated dipalmitoylphosphatidylcholine bilayers is described. The xenon atoms were found to prefer the interfacial and central regions of the bilayer. The presence of xenon was observed to lead to a small increase in the surface area, membrane thickness, and order of the acyl chains.

Anesthetics, Inhalation↗

PubMed2005

Modeling glycolipids: take one.

Molecular dynamics simulations of glycolipid bilayers consisting of 1,2-di-O-palmitoyl-3-O-beta-D-glucosyl-sn-glycerol were performed using five different force field parameterizations. Comparing the results with experimental data revealed that only the all-atom model correctly reproduces both the phase behavior and the surface area per lipid. Only one of the united atom models studied reproduces the correct phase behavior.

Computer Simulation↗