Cholesterol is not all bad. It is important in Membranes. We are studying Cholesterol Microstructures in Membranes.


Biophysics is one of the fastest growing areas of physics.  That’s part of the reason.  Another part is becasuse the most interesting questions in modern physics are in biology.

Biophysics is that branch of knowledge that applies the principles of physics and chemistry and the methods of mathematical analysis and computer modeling to understand the mechanisms of biological systems in terms of fundamental principles. It is one of the most diverse and interdisciplinary areas of study.

Biophysics is primarily a molecular science. It seeks to explain biological function in terms of the molecular structures and properties of specific molecules. The size of these molecules varies dramatically. The problems studied also vary from foundational studies of proteins to studies in the neurosciences to pharmacology.

Careers in Biophysics

Biophysics offers career opportunities in teaching and research at colleges, universities, medical centers, research institutes, and government laboratories, in research and development at pharmaceutical, biotechnology, nanotechnology, and other biologically-oriented high-tech companies, and in regulation and public policy at agencies of state and federal governments and at private organizations.

Biophysical Society Website:

Biophysics Research at Goshen College

Experiments:   We do experiments on model biological membranes and ion channels.

Theoretical: We do mathematical modeling of cholesterol microstructures and of ion channel dynamics using Monte Carlo simulations.


Introduction: Cell membranes are double layers (bilayers) of lipid molecules closely packed, as shown here. [Membranes shown here were adapted from part of Figure 5: Weber et al. (1998) Cell 92:759-72.]

Fig. 1: Phospholipid bilayer. Single lipid molecule (a), bilayer of lipid molecules with polar heads on outside of bilayer and acyl chains together (b).

At the junctions between nerve cells (neurons) the presynaptic membrane of the axon is separated from the postsynaptic membrane of the dendrite of the next neuron by a narrow synaotic cleft as shown here.

Fig. 2: Membranes at the Synapse between neurons.

The axon of a neuron contains sodium and potassium channels. These channels are opened sequentially down the axon of the neuron after excitation of the cell. This causes depolarization of the cell and constitutes the electrical pulse, which propagates along the cell.

At the synapse this electrical pulse results in release of calcium. The calcium causes a pore to form between vesicles containing the neurotransmitter docked at the presynaptic membrane and the synaptic cleft.  The sequence leading to release of neurotransmitter is shown in Fig. 3.

Fig. 3: Pore formation and endocytosis. Docked vesicle with neurotransmitter (a), Pore formed after Ca++ addition (b), endocytosis and total release of neurotransmtter (c).

A detailed picture of the transmission at the synapse is shown in the drawing in Fig. 4.

Fig. 4: Synapse. The axon has vesicles containing the neurotransmitter. Vesicles release neurotransmitter into synaptic cleft when excitation pulse reaches synapse. Excitation of dendrite results from neurotransmitter binding to receptors.

The receptors are ligand gated ion channels or metabotropic receptors.

Vesicle fusions can be studied using a model system in the laboratory. A phospholipid bilayer (membrane) is formed over a small (150-300μm) hole separating two salt solutions. The apparatus is shown in Fig. 5.

Fig. 5: Apparatus for vesicle fusion studies. The black plastic base is called the “block.” The chamber and the cup hold the two salt solutions. The vesicles are introduced into the chamber and stirred by the flea bar visible at the bottom of the Chamber.


The chamber for vesicles contains salt solution 1 with vesicles and the cup contains salt solution 2 (see Fig. 8). The hole in the cup, indicated by the white dot, is where the lipid bilayer is formed. The window is for viewing while “brushing” on the bilayer.

The bolck and cup are located inside a Faraday cage that shields against outside electrical interference from lights and other apparatus. Fig. 6 shows the apparatus with the Faraday cage open.

Fig. 6: Apparatus with block and cup in Faraday cage. The stir bar motor rotates a magnet that drives the flea stirbar in the forward (cis) chamber. The Vibration free table is air cushoned to prevent any mechanical disturbances.

The entire apparatus including amplifier and computer is shown in Fig. 7.


Fig. 7: Entire apparatus. The Faraday cage is closed. Visible are the pump for the vibration free table, the amplifier and filter, and the computer for running the experiment and recording the data.

In an experiment vesicles are introduced into the salt solution on one side of the bilayer. This is shown in Fig. 8.

Fig. 8: Phospholipid bilayer separating two salt solutions. One solution contains vesicles. These vesicles fuse with the bilayer.

The flea stir bar keeps the vesicles moving in the chamber and they will collide with the bilayer. Fusion will occur

If we put ion channels in the vesicles each vesicle fusion will transport those ion channels to the bilayer (see Fig. 3(c)). If we also impose an electrical potential across the bilayer we can measurethe electrical current resulting from the ion channels. We can then see each vesicle fusion as jump in the electrical current through the bilayer.

This is shown schematically in Fig. 9.

Fig. 9: Fusion of a vesicle containing channels with phospholipid bilayer. The graph below is a plot of the current through the bilayer as a function of time corresponding to the physical fusion above.

In our experiments we use the plant equivalent of cholesterol, which is ergosterol, because it is more effective.  The nystatin/ergosterol channels close as the ergosterol diffuses into the membrane.  This results in decay in the bilayer current in Fig. 9.

In Fig. 10 we show actual data for current spikes.

Fig. 10: Bilayer current spikes from fusion of vesicles containing nystatin/ergosterol channels.

The Big Surprise

We were first only interested in the number of fusions because we were interested in modeling the transmission of neuro transmitter by vesicles. The form of the current spikes was of initially of no particular interest. We expected the decay of the bilayer current to be exponential.

But we noticed that the decay of the bilayer current was linear. Data are shown in Fig. 11.

Fig. 11: Bilayer current spikes when mol fraction of ergosterol in vesicles is χerg = 0.4. The bilayer potential was held at 60mV. Red lines are linear fits. At χerg = 0.4 spikes of short and long duration were observed. Short duration (a), long duration (b).

The linear decay scheme was the real surprise. But we soon found that there was also a bimodal distribution in the spike durations. In retropspect this bimodal duration may be seen in the spikes in Fig. 10. The extreme difference in the time scales is more evident in the spikes in panels (a) and (b) of Fig. 11.

In Fig. 12 we have plotted the duration of the spikes τD in seconds against the mol fraction of ergosterol in the vesicles.


Fig. 12: Semilogarithmic plot of bilayer current spike duration as a function of mol fraction ergosterol in the vesicles. Error bars have been excluded for the sake of visualization.

Between χerg = 0.25 and χerg = 0.50 we see spikes of constant short duration and spikes with a long duration that increases approximately exponentially up to χerg = 0.4 and then drops precipitously at χerg = 0.5.

Neither the observation of a linear decay scheme nor the observation of a bimodal distribution of duration times is consistent with a simple picture of ergosterol diffusing away from fused vesicle remnants on thebilayer. We have proposed a model to account for these observations.

Our Model

The spike decay would be exponential if the nystatin/ergosterol channels were independent of one another. The linear decay indicates that they must be linked. We proposed a linkage of the nystatin to domains containing high concentration of ergosterol. Each domain then links a group of nystatin/ergosterol channels together and the channels dissolve as the linkage is lost. We proposed that the nystatin/ergosterol channels formed around the perimeter of the ersgosterol domains.

We also proposed that the domains contained erosterol in regular arrays known as superlattices. These were first proposed by Parkson Chong’s group at Temple University and had been studied by Juyang Huang at Texas Tech. In Fig. 13 we show the basic form of these superlattices


Fig. 13: Ergosterol (cholesterol) superlattices. The sterol occupies an acyl chain site in the lipid bilayer. Two open circles form each lipid. A single dark circle is a sterol

The model we have proposed is illustrated in 14.

Fig. 14: Our model of the arrangement of nystatin/ergosterol channels linked to an ergosterol rich domain. The open circles are acyl chains of the lipids forming the membrane, the dark filled circles are ergosterol molecules, and the large circles represent nystatin/ergosterol channels.

In Fig. 14 we have included a number of sterols diffusing out of the ergosterol-rich domain. In our model the ergosterol-rich domain dissolves through the diffusion of sterols into the open bilayer.

The actual form of the nystatin/ergosterol channels is unknown at this time. Manuel Prieto and his coworkers at the CQFM, Instituto Superior Te´ cnico, Lisbon, Portugal , have shown that nystatin/ergosterol channels contain between 4 and 12 nystatins. The nystatin monomers, lying on top of the membrane, rotate and penetrate the membrane to form channels. This is shown diagrammatically in Fig. 15.

Fig. 15: Diagram of nystatin/ergosterol channel formation from nystatin monomers.

In Fig. 15 we have made no attempt to indicate the presence of the necessary ergosterols for support of the channel.

In Fig. 16 we show a molecular modeling result for a nystatin channel in which there are six nystatin monomers.



Fig. 16: Nystatin channel formed by six nystatin monomers in a lipid bilayer. Nystatins are colored red, blue, green, yellow, orange and grey for visibility. Lipid heads are dull green. Insert gives side view through the bilayer. Pore is visible in the center of the channel [Dixon Woodbury, Brigham Young University

Because of the uncertainty in the actual form of the nystatin/ergosterol channel, particularly in the interaction with the sterol, we have made no attempt to include the details of channel structure in our model. Our model includes

·          Nystatin/ergosterol channels form around the perimeter of ergosterol-rich domains in lipid bilayers.

·          Ergosterol-rich domains are superlattice structures.

·          The ergosterol-rich domains dissolve by diffusion of the ergosterol into the empty bilayer.

The model has been tested against our own data and those of Dixon Woodbury. This has been pubished [Helrich, C. S., J. A. Schmucker, and D. J. Woodbury. 2006. Evidence that Nystatin Channels Form at the Boundaries, Not the Interiors of Lipid Domains. Biophys. J. 91:1116-1127].

We are presently pusuing tests of our model against further data. This is the subject of the next section.

Theoretical Work

We initially realized that our model was basically correct when we conducted a statistical mechanical analysis. We then decided to carry out numerical simulations of this using the Monte Carlo (MC) approach.

MC simulations are based on what is called a Markovian system. In a Markovian system each molecule moves from its location at a particular time to another location based on probabilities. But the molecule system retains no memory of the past locations of the molecules. So the beginning of each MC step requires ony information about where the molecules are. In a MC step each molecule is considered separately. A possible move is chosen and then “permission” is considered. The decisions of the possible move and the permission are based on random numbers. That is, the MC simulation is like gambling with a weighted “permission.” Hence the name “Monte Carlo”.

In our experiments vesicles are prepared containing a known percentage of ergosterol and sufficient nystatin to form channels. So the sesicle lipid membrane contains ergosterol-rich domains (Ergosterol Superlattics Domains, ESLDs), which are presumably surrounded by nystatin channels. When such a vesicle fuses with a planar bilayer, the vesicle membrane is transported to the bilayer. The bilayer is then an open lipid bilayer with a small remnant (patch) of ergosterol-rich domains surrounded by channels.

To model this situation we first form the remnat from a MC simulation at the sterol concentration desired. Then we place this remnant in the center of a large open bilayer and perform a second MC simulation to study the rate at which the ESLDs dissolve due to diffusion. We equate the number of channels to the perimeter, based on the assumption that there is a distinct linear density of channels around the perimeter of each ESLD.


In Fig. 17 we have an example of a vesicle bilayer resulting from a MC simulation.



Fig. 17: Remnant from vesicle fusion for χerg = 0.43. The remnant containins a mosaic of superlattices corresponding to χerg = 0.4 and 0.5.

When we place this remnant in a larger lattice of acyl chain sites (an open bilayer) and conduct a MC simulation, we obtain the dynamics of the breakup of the mosaic as shown in the snapshots in Fig. 18.


Fig. 18: Breakup of a fusion remnant placed in an open bilayer. The initial remnant, a mosaic of ESLDs, and time development is shown. The diffusion of sterols from the ESLDs can be seen.

Bilayer Current

A plot of the perimeter of the ESLDs as a function of MC step is shown in Fig. 19, right panel, and a corresponding experimental result is shown in the left panel. In both panels χerg = 0.43.

Fig. 19: Experimental data for bilayer current, left panel, and ESLD perimeter (number of channels) from MC simulation, right panel. In the right panel the two perimeters for the ESLDs corresponding to χerg = 0.4 and 0.5 are shown separately as well as the total perimeter.

We see that the MC simulation matches the experimental data in Fig. 19. We also note that separation of experimental spikes into a fast initial peak followed by a slower portion, which is observed in many experiments, can be understood from the distinct rates of decay of the two separate ESLDs.


Parkson Chong has conducted experiments on the fluorescence of dehydroergosterol (DHE) in vesicle bilayers. The presence of regular sterol regions, decreases fluorescence because the DHE molecules are then widely separated. Chong found fluorecence minima at particular values of χerg. Many of these minima are sharp and correspond to the critical values of (Cr) of χerg for which there are superlattices.

Our MC simulations of the mosaic appearing in the vesicle membrane (see Fig. 17) provide predictions of the number of sterols that are in superalttices. In Fig. 20 we show the result of our calculation compared to Chong’s experimental data. We have plotted our results for the number of sterols in the superlattice subtracted from an arbitrary number, because Chong’s data show minima when most of the sterols are, presumably, in superlattices.



Fig. 20: Comaprison of MC simulations of vesicle membranes from our model with Chong’s data. Both plots reflect disordered clustering of sterols in the membrane. Downward peaks show order in the domains.

Our results show the distinct ordering at χerg≈ 0.25. This corresponds to a superlattice (see Fig. 13) predicted geometrically. Snapshots of our MC results show this distinctly. There is also a geometrically predicted superlattice at χerg≈ 0.33. Both plots show a dip at this point, but neither shows a sharp peak. Snapshots of our results show a mosaic of the form shown in Fig. 17, but composed of χerg≈ 0.4 superlattices.

If we plot the percent sterols in each of theESLD types (superlattices) that we observed between χerg≈ 0.22 – 0.42 we obtain the plot shown in Fig. 21.


Fig. 21: Percent sterols in ESLDs.

According to Chong there should be a maximum in the percent sterols appearing on the boundaries of the ESLDs (of the superlattices) at Cr. n Fig. 22 we have plotted the number of sterols on ESLD boundaries as a function of χerg.


Fig. 22: Percent sterols on ESLD boundaries.

The requirement for a sterol to be inside an ESLD is more stringent than that for a boundary. This accounts for the fact that there are sterols on ESLD boundaries in regions where there are no sterols within ESLDs.


We have shown that our model produces the linear decay scheme we have consistently observed. We are also able to understand the fast and slow processes observed in terms of diffusion of sterols from different types of ESLDs.

Our model also predicts the results of the fluorescent experiments using DHE conducted by Chong. We also verify Chong’s theory regarding the percent of the sterols on the ESLD perimeters.

Our model explains why the fluorescent data at χerg = 0.33 do not form a sharp peak. The predicted superlattice at χerg = 0.33 is not observed. This peak results from an increase in the density of the χerg = 0.4 ESLDs.

In 1999 Dixon Woodbury noted that the average fusion spike current is a linear function of vesicle diameter not of vesicle surface area [Woodbury, D. J. 1999. Nystatin/Ergosterol method for reconstituting ion channels into planar lipid bilayers. In Methods in Enzymology. Ion Channels Part C. P. M. Conn, editor, 294(17):319-339] . This is shown in Fig. 23.


Fig. 23: Peak current as a function of vesicle diameter.

If the channels were evenly distributed over the vesicles the peak current, which is proportional to the number of channels, would be a linear function of the square of the diameter (area of the vesicle). However, if the channels are distributed around the perimeter of the ESLDs the number of channels would vary linearly with vesicle diameter. Our model then explains these data.


At this stage in the research our model has ben successful in explaining

  • The linear decay of bilayer current.
  • The fast and slow processes we have observed.
  • The fluorescent data of Chong.
  • The data for peak current of Woodbury.

And our model has verified the Chong theory regarding the percent sterols on the ESLD boundaries.

Juyang Huang has noted that we do not see a peak in the density of sterols in ESLDs at χerg = 0.4. He has noted a distinctive rise in the value of the chemical potential at χerg = 0.4. We presently have results from our model for chemical potential using the method of Kirkwood. These do show an increase in chemical potential at χerg = 0.4. We will report these results at the annual meeting of the Biophysical Society in 2009.

We then have rather strong evidence that the arrangement of sterols (cholesterol and ergosterol) in membranes is such that superlattices are formed.

Why is this important? It is possible, and maybe even plausible, that sterol rich regions are necessary for the incorporation of proteins in cell membranes, and for other processes.