Mechanical equilibrium of thick, hollow, liquid membrane cylinders.

The mechanical equilibrium of bilayer membrane cylinders is analyzed. The analysis is motivated by the observation that mechanically formed membrane strands (tethers) can support significant axial loads and that the tether radius varies inversely with the axial force. Previously, thin shell theory has been used to analyze the tether formation process, but this approach is inadequate for describing and predicting the equilibrium state of the tether itself. In the present work the membrane is modeled as two adjacent, thick, anisotropic liquid shells. The analysis predicts an inverse relationship between axial force and tether radius, which is consistent with experimental observation. The area expansivity modulus and bending stiffness of the tether membrane are calculated using previously measured values of tether radii. These calculated values are consistent with values of membrane properties measured previously. Application of the analysis to precise measurements of the relationship between tether radius and axial force will provide a novel method for determining the mechanical properties of biomembrane.     

Spectrin, human erythrocyte shapes, and mechanochemical properties.

Physical studies of human erythrocyte spectrin indicate that isolated spectrin dimers and tetramers in solution are worm-like coils with a persistence length of approximately 20 nm. This finding, the known polyelectrolytic nature of spectrin, and other structural information about spectrin and the membrane skeleton molecular organization have lead us to the hypothesis that the human erythrocyte membrane skeleton constitutes a two-dimensional ionic gel (swollen ionic elastomer). This concept is incorporated in what we refer to as the protein gel-lipid bilayer membrane model. The model accounts quantitatively for red elastic shear modulus and the maximum elastic extension ratio reported for the human erythrocytes membrane. Gel theory further predicts that depending on the environmental conditions, the membrane skeleton modulus of area compression may be small or large relative to the membrane elastic shear modulus. Our analyses show that the ratio between these two parameters affects both the geometry and the stability of the favored cell shapes and that the higher the membrane skeleton compressibility the smaller the values of the gel tension needed to induce cell shape transformations. The main virtue of the protein gel-lipid bilayer membrane model is that it offers a novel theoretical and molecular basis for the various mechanical properties of the membrane skeleton such as the membrane skeleton modulus of area compression and osmotic tension, and the effects of these properties on local membrane skeleton density, cell shape, and shape transformations.     

An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins

Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available.     

Membrane adsorption,folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

History of the membrane (pump) theory of the living cell from its beginning in mid-19th century to its disproof 45 years ago--though still taught worldwide today as established truth

The concept that the basic unit of all life, the cell, is a membrane-enclosed soup of (free) water, (free) K+ (and native) proteins is called the membrane theory. A careful examination of past records shows that this theory has no author in the true sense of the word. Rather, it grew mostly out of some mistaken ideas made by Theodor Schwann in his Cell Theory. (This is not to deny that there is a membrane theory with an authentic author but this authored membrane theory came later and is much more narrowly focussed and accordingly can at best be regarded as an offshoot of the broader and older membrane theory without an author.) However, there is no ambiguity on the demise of the membrane theory, which occurred more than 60 years ago, when a flood of converging evidence showed that the asymmetrical distribution of K+ and Na+ observed in virtually all living cells is not the result of the presence of a membrane barrier that permits some solutes like water and K+ to move in and out of the cell, while barring--absolutely and permanently--the passage of other solutes like Na+. To keep the membrane theory afloat, submicroscopic pumps were installed across the cell membrane to maintain, for example, the level of Na+ in the cell low and the level of K+ high by the ceaseless pumping activities at the expense of metabolic energy. Forty-five year ago this version of the membrane theory was also experimentally disproved. In spite of all these overwhelming evidence against the membrane-pump theory, it still is being taught as verified truth in all high-school and biology textbooks known to us today. Meanwhile, almost unnoticed, a new unifying theory of the living cell, called the association-induction hypothesis came into being some 40 years ago. Also little noticed was the fact that it has received extensive confirmation worldwide and has shown an ability to provide self-consistent interpretations of most if not all known experimental observations that are contradicting the membrane-pump theory as well as other observations that seem to support the membrane pump theory.     

A method to relate steady-state ionic currents, conductances, and membrane potential in ion exchange membranes with unknown thermodynamic properties

A method is presented by which the steady-state properties of an homogeneous, permselective membrane at uniform temperature can be predicted without knowledge of its thermodynamic properties other than assuming that they are functions only of local mole fractions in the membrane. By making this assumption, it is shown how the ionic conductances can be calculated at any point in the membrane from two sets of measurements, (a) R(symm), the steady-state resistance of the membrane measured between identical solutions and (b) V(0), the potential difference between nonidentical solutions for zero current. These two parameters are measured at different external solution compositions (e.g. a varying sodium-potassium ratio ranging from zero to infinity). From these measurements it is shown how the flux equations may be integrated without a knowledge of mobilities, activity coefficients, and other interior membrane parameters. The application of the method to fixed site membranes with variable mobilities is described and the theory for this particular case has also been verified experimentally in glass membranes.1 A possible application to biological membranes is discussed and a comparison is made between the present treatment and previous treatments used to calculate the steady-state properties of cell membranes, notably the theory of Teorell, Meyer, and Sievers and the constant field theory.     

The Energetics of Membrane Fusion from Binding,through Hemifusion,Pore Formation,and Pore Enlargement

The main steps of viral membrane fusion are local membrane approach, hemifusion, pore formation, and pore enlargement. Experiments and theoretical analyses have helped determine the relative energies required for each step. Key protein structures and conformational changes of the fusion process have been identified. The physical deformations of monolayer bending and lipid tilt have been applied to the steps of membrane fusion. Experiment and theory converge to strongly indicate that, contrary to former conceptions, the fusion process is progressively more energetically difficult: hemifusion has a relatively low energy barrier, pore formation is more energy-consuming, and pore enlargement is the most difficult to achieve.This revised version was published online in August 2005 with a corrected cover date.     

A theoretical formalism for aggregation of peroxidized lipids and plasma membrane stability during photolysis.

The objective of this investigation was to examine, from a theoretical perspective, the mechanism underlying the lysis of plasma membranes by photoinduced, chemically mediated damage such as is found in photolysis. Toward this end, a model is presented which relates the membrane lifetime to the thermodynamic parameters of the membrane components based upon the kinetic theory of aggregate formation. The formalism includes a standard birth/death process for the formation of damaged membrane components (i.e., peroxidized lipids) as well as a terminating condensation process for the formation of aggregates of peroxidized plasma membrane lipids. Our theory predicts that 1) the membrane lifetime is inversely correlated with predicted rate of membrane damage; 2) an upper limit on the duration of membrane damage exists, above which the mean and variance of the membrane lifetime is independent of further membrane damage; and 3) both the mean and variance of the time of membrane lifetime distribution are correlated with the number of sites that may be damaged to form a single membrane defect. The model provides a framework to optimize the lysis of cell membranes by photodynamic therapy.     

Basal metabolic rate: history, composition, regulation, and usefulness

The concept of basal metabolic rate (BMR) was developed to compare the metabolic rate of animals and initially was important in a clinical context as a means of determining thyroid status of humans. It was also important in defining the allometric relationship between body mass and metabolic rate of mammals. The BMR of mammals varies with body mass, with the same allometric exponent as field metabolic rate and with many physiological and biochemical rates. The membrane pacemaker theory proposes that the fatty acid composition of membrane bilayers is an important determinant of a species BMR. In both mammals and birds, membrane polyunsaturation decreases and monounsaturation increases with increasing body mass and a decrease in mass-specific BMR. The secretion and production of thyroid hormones in mammals are related to body mass, with the allometric exponent similar to BMR; yet there is no body size-related variation in either total or free concentrations of thyroid hormones in plasma of mammals. It is suggested that in different-sized mammals, the secretion/production of thyroid hormones is a result of BMR differences rather than their cause. BMR is a useful concept in some situations but not in others.     

Acidic lipids, H-ATPases, and mechanism of oxidative phosphorylation. Physico-chemical ideas 30 years after P. Mitchell's Nobel Prize award

Peter D. Mitchell, who was awarded the Nobel Prize in Chemistry 30 years ago, in 1978, formulated the chemiosmotic theory of oxidative phosphorylation. This review initially analyzes the major aspects of this theory, its unresolved problems, and its modifications. A new physico-chemical mechanism of energy transformation and coupling of oxidation and phosphorylation is then suggested based on recent concepts regarding proteins, including ATPases that work as molecular motors, and acidic lipids that act as hydrogen ion (H⁺) carriers. According to this proposed mechanism, the chemical energy of a redox substrate is transformed into nonequilibrium states of electron-transporting chain (ETC) coupling proteins. This leads to nonequilibrium pumping of H⁺ into the membrane. An acidic lipid, cardiolipin, binds with this H⁺ and carries it to the ATP-synthase along the membrane surface. This transport generates gradients of surface tension or electric field along the membrane surface. Hydrodynamic effects on a nanolevel lead to rotation of ATP-synthase and finally to the release of ATP into aqueous solution. This model also explains the generation of a transmembrane protonmotive force that is used for regulation of transmembrane transport, but is not necessary for the coupling of electron transport and ATP synthesis.     

Effects of divalent cations,temperature, osmotic pressure gradient,and vesicle curvature on phosphatidylserine vesicle fusion

Summary Fusion of phosphatidylserine vesicles induced by divalent cations, temperature and osmotic pressure gradients across the membrane was studied with respect to variations in vesicle size. Vesicle fusion was followed by two different methods: 1) the Tb/DPA fusion assay, whereby the fluorescent intensity upon mixing of the internal aqueous contents of fused lipid vesicles was monitored, and 2) measurement of the changes in turbidity of the vesicle suspension due to vesicle fusion. It was found that the threshold concentration of divalent cations necessary to induce vesicle fusion depended on the size of vesicles; as the diameter of the vesicle increased, the threshold value increased and the extent of fusion became less. For the osmotic pressure-induced vesicle fusion, the larger the diameter of vesicles, the smaller was the osmotic pressure gradient required to induce membrane fusion. Divalent cations, temperature increase and vesicle membrane expansion by osmotic pressure gradient all resulted in increase in surface energy (tension) of the membrane. The degree of membrane fusion correlated with the corresponding surface energy changes of vesicle membranes due to the above fusion-inducing agents. The increase in surface energy of 9.5 dyn/cm from the reference state corresponded to the threshold point of phosphatidylserine membrane fusion. An attempt was made to explain the factors influencing fusion phenomena on the basis of a single unifying theory.     

Filtration,diffusion, and molecular sieving through porous cellulose membranes

1. A study has been made of the diffusion and filtration of a graded series of molecules (including tritium-labelled water, urea, glucose, antipyrine, sucrose, raffinose, and hemoglobin) in aqueous solution through porous cellulose membranes of three degrees of porosity. 2. Experimental results were in close agreement with predictions based on the membrane pore theory of Pappenheimer et al. (1,2). Restriction to molecular diffusion is a function of pore radius and molecular radius described by equation (11) in the text. Molecular sieving during ultrafiltration is a function of total pore area per unit path length, pore radius, molecular radius, and filtration rate given by equations (16) and (19). 3. Estimates of average pore radius made by means of this theory were considerably larger than estimates made by the method of Elford and Ferry (3) (Table II). Sources of error in the latter method are discussed and a new method of membrane calibration is proposed in which the total cross-sectional area of the pores is measured by direct diffusion of isotope-labelled water. 4. Steady-state osmotic pressures of solutions of sucrose and raffinose measured during molecular sieving through cellulose membranes were found to be close to the "ideal" osmotic pressures calculated by van't Hoff's law. Thus the present experimental data support the methods used by Pappenheimer et al. in their studies on living capillary walls as well as their theory of membrane pore permeability.     

Trapping by Clusters of Channels,Receptors, and Transporters: Quantitative Description

Various membrane functional units such as receptors, transporters, and channels, whose action necessarily involves capturing diffusing molecules, are often organized into multimeric complexes forming clusters on the cell and organelle membranes. These functional units themselves are usually oligomers of several integral proteins, which have their own symmetry. Depending on the symmetry, they form clusters on different packing lattices. Moreover, local membrane inhomogeneities, e.g., the so-called membrane domains, rafts, stalks, etc., lead to different patterns even within the structures on the same packing lattice. Units in the cluster compete for diffusing molecules and screen each other. Here we propose a general approach that allows one to quantify the screening effects. The approach is used to derive simple approximate formulas giving the trapping rates of diffusing molecules by clusters of absorbers on lattices of different packing symmetries. The obtained results describe smooth variation of the trapping rate from the sum of the rates of individual absorbers forming the cluster to the effective collective rate. The latter shows how the trapping efficiency of an individual absorber decreases as the number of absorbers in the cluster increases and/or the inter-absorber distance decreases. Numerical tests demonstrate good agreement between the rates predicted by the theory and obtained from Brownian dynamics simulations for clusters of different shapes and sizes.     

A Multiscale Red Blood Cell Model with Accurate Mechanics, Rheology, and Dynamics

Red blood cells (RBCs) have highly deformable viscoelastic membranes exhibiting complex rheological response and rich hydrodynamic behavior governed by special elastic and bending properties and by the external/internal fluid and membrane viscosities. We present a multiscale RBC model that is able to predict RBC mechanics, rheology, and dynamics in agreement with experiments. Based on an analytic theory, the modeled membrane properties can be uniquely related to the experimentally established RBC macroscopic properties without any adjustment of parameters. The RBC linear and nonlinear elastic deformations match those obtained in optical-tweezers experiments. The rheological properties of the membrane are compared with those obtained in optical magnetic twisting cytometry, membrane thermal fluctuations, and creep followed by cell recovery. The dynamics of RBCs in shear and Poiseuille flows is tested against experiments and theoretical predictions, and the applicability of the latter is discussed. Our findings clearly indicate that a purely elastic model for the membrane cannot accurately represent the RBC's rheological properties and its dynamics, and therefore accurate modeling of a viscoelastic membrane is necessary.     

Generation of membrane potential beyond the conceptual range of Donnan theory and Goldman-Hodgkin-Katz equation

Donnan theory and Goldman-Hodgkin-Katz equation (GHK eq.) state that the nonzero membrane potential is generated by the asymmetric ion distribution between two solutions separated by a semipermeable membrane and/or by the continuous ion transport across the semipermeable membrane. However, there have been a number of reports of the membrane potential generation behaviors in conflict with those theories. The authors of this paper performed the experimental and theoretical investigation of membrane potential and found that (1) Donnan theory is valid only when the macroscopic electroneutrality is sufficed and (2) Potential behavior across a certain type of membrane appears to be inexplicable on the concept of GHK eq. Consequently, the authors derived a conclusion that the existing theories have some limitations for predicting the membrane potential behavior and we need to find a theory to overcome those limitations. The authors suggest that the ion adsorption theory named Ling’s adsorption theory, which attributes the membrane potential generation to the mobile ion adsorption onto the adsorption sites, could overcome those problems.     

Coarse-grained molecular dynamics simulations of phase transitions in mixed lipid systems containing LPA, DOPA, and DOPE lipids

The mechanisms that mediate biomembrane shape transformations are of considerable interest in cell biology. Recent in vitro experiments show that the chemical transformation of minor membrane lipids can induce dramatic shape changes in biomembranes. Specifically, it was observed that the addition of DOPA to DOPE has no effect on the stability of the bilayer structure of the membrane. In contrast, the addition of LPA to DOPE stabilizes the bilayer phase of DOPE, increasing the temperature of a phase transition from the bilayer to the inverted hexagonal phase. This result suggests that the chemical conversion of DOPA to LPA is sufficient for triggering a dramatic change in the shape of biomembranes. The LPA/DOPA/DOPE mixture of lipids provides a simple model system for understanding the molecular events driving the shape change. In this work, we used coarse-grained molecular dynamics simulations to study the phase transitions of this lipid mixture. We show that despite the simplicity of the coarse-grained model, it reproduces the experimentally observed phase changes of: 1), pure LPA and DOPA with respect to changes in the concentration of cations; and 2), LPA/DOPE and DOPA/DOPE mixtures with respect to temperature. The good agreement between the model and experiments suggests that the computationally inexpensive coarse-grained approach can be used to infer macroscopic membrane properties. Furthermore, analysis of the shape of the lipid molecules demonstrates that the phase behavior of single-lipid systems is consistent with molecular packing theory. However, the phase stability of mixed lipid systems exhibits significant deviations from this theory, which suggests that the elastic energy of the lipids, neglected in the packing theory, plays an important role.     

Crowding-induced membrane remodeling: Interplay of membrane tension,polymer density,architecture

The plasma membrane hosts a wide range of biomolecules, mainly proteins and carbohydrates, that mediate cellular interactions with its environment. The crowding of such biomolecules regulates cellular morphologies and cellular trafficking. Recent discoveries have shown that the structure and density of cell surface polymers and hence the signaling machinery change with the state of the cell, especially in cancer progression. The alterations in membrane-attached glycocalyx and glycosylation of proteins and lipids are common features of cancer cells. The overexpression of glycocalyx polymers, such as mucin and hyaluronan, strongly correlates with cancer metastasis. Here, we present a mesoscale biophysics-based model that accounts for the shape regulation of membranes by crowding of membrane-attached biopolymer-glycocalyx and actin networks. Our computational model is based on the dynamically triangulated Monte Carlo model for membranes and coarse-grained representations of polymer chains. The model allows us to investigate the crowding-induced shape transformations in cell membranes in a tension- and graft polymer density-dependent manner. Our results show that the number of membrane protrusions and their shape depend on membrane tension, with higher membrane tension inducing more tubular protrusions than the vesicular shapes formed at low tension at high surface coverage of polymers. The shape transformations occur above the threshold density predicted by the polymer brush theory, but this threshold also depends on the membrane tension. Increasing the size of the polymer, either by changing the length or by adding side chains, is shown to increase the crowding-induced curvature. The effect of crowding is more prominent for flexible polymers than for semiflexible rigid polymers. We also present an extension of the model that incorporates properties of the actin-like filament networks and demonstrate how tubular structures can be generated by biopolymer crowding on the cytosolic side of cell membranes.     

Donnan membrane equilibrium,sedimentation equilibrium,and coil expansion of DNA in salt solutions

This is a review of applications of the McMillan-Mayer-Hill virial theory and the ionic double-layer theory to dilute colloidal solutions, in particular, solutions of DNA. Interactions of highly charged colloidal rods are developed in terms of the second virial coefficients between two rods, and between one rod and one small co-ion. The relevant cluster integrals are evaluated with interaction potentials based on the Poisson-Boltzmann equation. The treatment is extended to the intrachain repulsion responsible for the statistical swelling of coiled DNA (excluded volume effect). The theory is compared with three sets of experimental data: The salt distribution in Donnan membrane equilibria of DNA-salt solutions, sedimentation equilibria of short DNA fragments at different ionic strengths, and the intrinsic viscosity of T7 DNA in NaCl solutions. In all cases the theory agrees well with the experiments. The agreement is not convincing for the sedimentation equilibrium at low ionic strength, because here the experimental DNA concentration is too high for the truncated dilute solution expansion of the DNA-salt repulsion.