Biomedical subjects
Publications and source records attributed to Mikko Karttunen.
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.
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.
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.
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.
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.
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.
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.
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.
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.
We generalize the Bergman-Milton spectral representation, originally derived for a two-component composite, to extract the spectral density function for the effective dielectric constant of a graded composite. This work has been motivated by a recent study of the optical absorption spectrum of a graded metallic film [Huang and Yu, Appl. Phys. Lett. 85, 94 (2004)] in which a broad surface-plasmon absorption band was shown to be responsible for enhanced nonlinear optical response and an attractive figure of merit. It turns out that, unlike in the case of homogeneous constituent components, the characteristic function of a graded composite is a continuous function because of the continuous variation of the dielectric function within the constituent components. Analytical generalization to three-dimensional graded composites is discussed, and numerical calculations for multilayer composites are given as a simple application.
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.
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.
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.
Free volume pockets or voids are important to many biological processes in cell membranes. Free volume fluctuations are a prerequisite for diffusion of lipids and other macromolecules in lipid bilayers. Permeation of small solutes across a membrane, as well as diffusion of solutes in the membrane interior are further examples of phenomena where voids and their properties play a central role. Cholesterol has been suggested to change the structure and function of membranes by altering their free volume properties. We study the effect of cholesterol on the properties of voids in dipalmitoylphosphatidylcholine (DPPC) bilayers by means of atomistic molecular dynamics simulations. We find that an increasing cholesterol concentration reduces the total amount of free volume in a bilayer. The effect of cholesterol on individual voids is most prominent in the region where the steroid ring structures of cholesterol molecules are located. Here a growing cholesterol content reduces the number of voids, completely removing voids of the size of a cholesterol molecule. The voids also become more elongated. The broad orientational distribution of voids observed in pure DPPC is, with a 30% molar concentration of cholesterol, replaced by a distribution where orientation along the bilayer normal is favored. Our results suggest that instead of being uniformly distributed to the whole bilayer, these effects are localized to the close vicinity of cholesterol molecules.
A first-principles numerical model for crumpling of a stiff tethered membrane is introduced. This model displays wrinkles, ridge formation, ridge collapse, and initiation of stiffness divergence. The amplitude and wavelength of the wrinkles and the scaling exponent of the stiffness divergence are consistent with both theory and experiment. Close to the stiffness divergence further buckling is hindered by the nonzero thickness of the membrane, and its elastic behavior becomes similar to that of dry granular media. No change in the distribution of contact forces can be observed at the crossover, implying that the network of ridges is then simultaneously a granular force-chain network.
It is well known that in two dimensions Turing systems produce spots, stripes and labyrinthine patterns, and in three dimensions lamellar and spherical structures, or their combinations, are observed. In this paper we study transitions between these states in both two and three dimensions. First, we derive the regions of stability for different patterns using nonlinear bifurcation analysis. Then, we apply large scale computer simulations to analyze the pattern selection in a bistable system by studying the effect of parameter selection on morphological clustering and the appearance of topological defects. The method elaborated in this paper presents a probabilistic approach for studying pattern selection in a bistable reaction-diffusion system.
We construct a coarse-grained (CG) model for dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayers and apply it to large-scale simulation studies of lipid membranes. Our CG model is a two-dimensional representation of the membrane, where the individual lipid and sterol molecules are described by pointlike particles. The effective intermolecular interactions used in the model are systematically derived from detailed atomic-scale molecular dynamics simulations using the Inverse Monte Carlo technique, which guarantees that the radial distribution properties of the CG model are consistent with those given by the corresponding atomistic system. We find that the coarse-grained model for the DPPC/cholesterol bilayer is substantially more efficient than atomistic models, providing a speedup of approximately eight orders of magnitude. The results are in favor of formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations, in agreement with the experimental phase diagram of the system. We also explore the limits of the coarse-grained model, and discuss the general validity and applicability of the present approach.
The frequency-dependent nonlinear dielectric increment of dipolar fluids in nonpolar fluids is often measured by using a stationary relaxation method in which two electric fields are used: The static direct current (dc) field of high strength causing the dielectric nonlinearity, and the probing alternating current (ac) field of low strength and high frequency. When a nonlinear composite is subjected to a sinusoidal electric field, the electric response in the composite will, in general, consist of ac fields at frequencies of higher-order harmonics. Based on the Fröhlich model, we present a theory to investigate the nonlinear ac responses of dipolar fluids containing both polarizable monomers and dimers. In the case of monomers only, our theory reproduces the known results. We obtain the fundamental, second-, and third-order harmonics of the Fröhlich field by performing a perturbation expansion. The even-order harmonics are induced by the coupling between the ac and dc fields, although the system under consideration has a cubic nonlinearity only. The harmonics of the Fröhlich field can be affected by the field frequency, temperature, dispersion strength, and the characteristic frequency of the dipolar fluid, as well as the dielectric constant of the nonpolar fluid. The results are found to be in agreement with recent experimental observations.