Lipid transfer between vesicles: effect of high vesicle concentration.
(73/16154)
The problem of the desorption of a lipid molecule from a lipid vesicle (donor) and its incorporation into another vesicle (acceptor) at high acceptor concentrations, which has been investigated experimentally (Jones, J. D. and Thompson, T. E., 1990. Biochemistry, 29:1593-1600), is analyzed here from a theoretical point of view, formulated in terms of the diffusion equation with appropriate boundary conditions. The goal is to determine whether or not the observed acceleration of the off-rate from a donor is caused by interaction with an acceptor vesicle at short range, or is simply the result of statistical effects due the proximity of the acceptor and its influence on the probability of the test lipid returning to the donor. We establish a correspondence between the theoretical parameters and the experimental, thermodynamic and dynamic variables entering the problem. The solution shows that, because of the extremely high Gibbs activation energy for desorption of a phospholipid, the process would always be first-order, even at very high vesicle concentrations. This means that acceleration of the off-rate must be due to donor-acceptor interactions at short distances, as proposed in the experimental work. (+info)
Simulation study of a gramicidin/lipid bilayer system in excess water and lipid. II. Rates and mechanisms of water transport.
(74/16154)
A gramicidin channel in a fluid phase DMPC bilayer with excess lipid and water has been simulated. By use of the formal correspondence between diffusion and random walk, a permeability for water through the channel was calculated, and was found to agree closely with the experimental results of Rosenberg and Finkelstein (Rosenberg, P.A., and A. Finkelstein. 1978. J. Gen. Physiol. 72:327-340; 341-350) for permeation of water through gramicidin in a phospholipid membrane. By using fluctuation analysis, components of resistance to permeation were computed for movement through the channel interior, for the transition step at the channel mouth where the water molecule solvation environment changes, and for the process of diffusion up to the channel mouth. The majority of the resistance to permeation appears to occur in the transition step at the channel mouth. A significant amount is also due to structurally based free energy barriers within the channel. Only small amounts are due to local friction within the channel or to diffusive resistance for approaching the channel mouth. (+info)
Vesicle deformation by an axial load: from elongated shapes to tethered vesicles.
(75/16154)
A sufficiently large force acting on a single point of the fluid membrane of a flaccid phospholipid vesicle is known to cause the formation of a narrow bilayer tube (tether). We analyze this phenomenon by means of general mathematical methods allowing us to determine the shapes of strongly deformed vesicles including their stability. Starting from a free vesicle with an axisymmetric, prolate equilibrium shape, we consider an axial load that pulls (or pushes) the poles of the vesicle apart. Arranging the resulting shapes of strained vesicles in dependence of the axial deformation and of the area difference of monolayers, phase diagrams of stable shapes are presented comprising prolate shapes with or without equatorial mirror symmetry. For realistic values of membrane parameters, we study the force-extension relation of strained vesicles, and we demonstrate in detail how the initially elongated shape of an axially stretched vesicle transforms into a shape involving a membrane tether. This tethering transition may be continuous or discontinuous. If the free vesicle is mirror symmetric, the mirror symmetry is broken as the tether forms. The stability analysis of tethered shapes reveals that, for the considered vesicles, the stable shape is always asymmetric (polar), i.e., it involves only a single tether on one side of the main vesicle body. Although a bilayer tube formed from a closed vesicle is not an ideal cylinder, we show that, for most practical purposes, it is safe to assume a cylindrical geometry of tethers. This analysis is supplemented by the documentation of a prototype experiment supporting our theoretical predictions. It shows that the currently accepted model for the description of lipid-bilayer elasticity (generalized bilayer couple model) properly accounts for the tethering phenomenon. (+info)
Dynamical properties of phospholipid bilayers from computer simulation.
(76/16154)
We present the results of a 10-ns molecular dynamics simulation of a dipalmitoylphosphatidylcholine/water system. The main emphasis of the present study is on the investigation of the stability over a long time and the dynamic properties of the water/membrane system. The motion of the lipid molecules is characterized by the center of mass movement and the displacement of individual atom groups. Because of the slow movement of the headgroup atoms, their contributions to the dipole potential vary slowly and with a large amplitude. Nevertheless, the water molecules compensate the strong fluctuations and maintain an almost constant total dipole potential. From the lateral displacement of the center of masses, we calculate the lateral diffusion coefficient to be Dlat = (3 +/- 0.6) x 10(-7) cm2/s, in agreement with neutron scattering results. The rotational motion is also investigated in our simulations. The calculated value for the rotational diffusion coefficient parallel to the molecular long axis, D = (1.6 +/- 0.1) x 10(8) s-1, is in good agreement with the experiment. (+info)
Monte Carlo simulation of two-component bilayers: DMPC/DSPC mixtures.
(77/16154)
In this paper, we describe a relatively simple lattice model of a two-component, two-state phospholipid bilayer. Application of Monte Carlo methods to this model permits simulation of the observed excess heat capacity versus temperature curves of dimyristoylphosphatidylcholine (DMPC)/distearoylphosphatidylcholine (DSPC) mixtures as well as the lateral distributions of the components and properties related to these distributions. The analysis of the bilayer energy distribution functions reveals that the gel-fluid transition is a continuous transition for DMPC, DSPC, and all DMPC/DSPC mixtures. A comparison of the thermodynamic properties of DMPC/DSPC mixtures with the configurational properties shows that the temperatures characteristics of the configurational properties correlate well with the maxima in the excess heat capacity curves rather than with the onset and completion temperatures of the gel-fluid transition. In the gel-fluid coexistence region, we also found excellent agreement between the threshold temperatures at different system compositions detected in fluorescence recovery after photobleaching experiments and the temperatures at which the percolation probability of the gel clusters is 0.36. At every composition, the calculated mole fraction of gel state molecules at the fluorescence recovery after photobleaching threshold is 0.34 and, at the percolation threshold of gel clusters, it is 0.24. The percolation threshold mole fraction of gel or fluid lipid depends on the packing geometry of the molecules and the interchain interactions. However, it is independent of temperature, system composition, and state of the percolating cluster. (+info)
A microscopic interaction model of maximum solubility of cholesterol in lipid bilayers.
(78/16154)
We recently reported the equilibrium maximum solubility of cholesterol in a lipid bilayer, chi*chol, to be 0.66 in four different phosphatidylcholines, and 0.51 in a phosphatidylethanolamine (Huang, J.,J.T. Buboltz, and G. W. Feigenson. 1999. Biochim. Biophys. Acta. in press). Here we present a model of cholesterol-phospholipid mixing that explains these observed values of chi*chol. Monte Carlo simulations show that pairwise-additivity of nearest-neighbor interactions is inadequate to describe all the chi*chol values. Instead, if cholesterol multibody interactions are assigned highly unfavorable energy, then jumps occur in cholesterol chemical potential that lead to its precipitation from the bilayer. Cholesterol precipitation is most likely to occur near three discrete values of cholesterol mole fraction, 0.50, 0.57, and 0.67, which correspond to cholesterol/phospholipid mole ratios of 1/1, 4/3, and 2/1, respectively. At these solubility limits, where cholesterol chemical potential jumps, the cholesterol-phospholipid bilayer mixture forms highly regular lipid distributions in order to minimize cholesterol-cholesterol contacts. This treatment shows that dramatic structural and thermodynamic changes can occur at particular cholesterol mole fractions without any stoichiometric complex formation. The physical origin of the unfavorable cholesterol multibody interaction is explained by an "umbrella model": in a bilayer, nonpolar cholesterol relies on polar phospholipid headgroup coverage to avoid the unfavorable free energy of cholesterol contact with water. Thus, at high cholesterol mole fraction, this unfavorable free energy, not any favorable cholesterol-phospholipid interaction, dominates the mixing behavior. This physical origin also explains the "cholesterol condensing effect" and the increase in acyl chain order parameter in cholesterol-phospholipid mixtures. (+info)
Time-resolved absorption and photothermal measurements with sensory rhodopsin I from Halobacterium salinarum.
(79/16154)
An expansion accompanying the formation of the first intermediate in the photocycle of transducer-free sensory rhodopsin I (SRI) was determined by means of time-resolved laser-induced optoacoustic spectroscopy. For the native protein (SRI-WT), the absolute value of the expansion is approximately 5.5 mL and for the mutant SRI-D76N, approximately 1.5 mL per mol of phototransformed species (in 0.5 M NaCl), calculated by using the formation quantum yield for the first intermediate (S610) of Phi610 = 0.4 +/- 0.05 for SRI-WT and 0.5 +/- 0.05 for SRI-D76N, measured by laser-induced optoacoustic spectroscopy and by laser flash photolysis. The similarity in Phi610 and in the determined value of the energy level of S610, E610 = (142 +/- 12) kJ/mol for SRI-WT and SRI-D76N indicates that Asp76 is not directly involved in the first step of the phototransformation. The increase with pH of the magnitude of the structural volume change for the formation of S610 in SRI-WT and in SRI-D76N upon excitation with 580 nm indicates also that amino acids other than Asp76, and other than those related to the Schiff base, are involved in the process. The difference in structural volume changes as well as differences in the activation parameters for the S610 decay should be attributed to differences in the rigidity of the cavity surrounding the chromophore. Except for the decay of the first intermediate, which is faster than in the SRI-transducer complex, the rate constants of the photocycle for transducer-free SRI in detergent suspension are strongly retarded with respect to wild-type membranes (this comparison should be done with great care because the preparation of both samples is very different). (+info)
Reversibility and hierarchy of thermal transition of hen egg-white lysozyme studied by small-angle x-ray scattering.
(80/16154)
To clarify mechanisms of folding and unfolding of proteins, many studies of thermal denaturation of proteins have been carried out at low protein concentrations because in many cases thermal denaturation accompanies a great tendency of aggregation. As small-angle x-ray scattering (SAXS) measurements are liable to use low-concentration solutions of proteins to avoid aggregation, SAXS has been regarded as very difficult to observe detailed features of thermal structural transitions such as intramolecular structural changes. By using synchrotron radiation SAXS, we have found that the presence of repulsive interparticle interaction between proteins can maintain solute particles separately to prevent further aggregation in thermal denaturation processes and that under such conditions the thermal structural transition of hen egg-white lysozyme (HEWL) holds high reversibility even at 5% w/v HEWL below pH approximately 5. Because of the use of the high concentration of the solutions, the scattering data has enough high-statistical accuracy to discuss the thermal structural transition depending on the structural hierarchy. Thus, the tertiary structural change of HEWL starts from mostly the onset temperature determined by the differential scanning calorimetry measurement, which accompanies a large heat absorption, whereas the intramolecular structural change, corresponding to the interdomain correlation and polypeptide chain arrangement, starts much prior to the above main transition. The present finding of the reversible thermal structural transitions at the high protein concentration is expected to enable us to analyze multiplicity of folding and unfolding processes of proteins in thermal structural transitions. (+info)