Cell Shape Model

Background

One of our research themes focuses on the ability of L-form bacteria to fit the mould and/or fill spaces by growing into them. We would like to explore the biophysical properties of the cell and the limits of the cell size the L-form can grow to before dividing by constructing a 2-dimentional cell model.

The shape of the bacterial cell is usually determined by the structures of the cell wall and bacterial cytoskeleton. The cell wall of B.subtilis is a thick layer of peptidoglycan which contains cell-wall proteins, techoic and lipotechoic acids. Peptidoglycan is a rigid structure which helps to protect the cell from osmotic pressure and mechanical damage. It also helps the cell to maintain a constant shape, while not specifically defining it. This role is played by the cytoskeleton. The bacterial cytoskeleton is largely similar to that of eukaryotes. It is comprised of three types of protein structures: FtsZ, MreB and crescentin (only present in certain species) which are homologous to eukaryotic tubulin, actin and intermediate filaments respectively. These proteins also perform other functions in the cell; usually they are involved in cell growth and division and therefore are often essential for the cell's survival (Shih and Rothfield, 2006). The communication between the two regulatory components is assured by a series of membrane-associated enzymes (penicillin binding proteins)(Cabeen, Jacobs-Wagner, 2005).

Because the L-forms lack the peptidoglycan cell wall, cells naturally adopt the most energetically favourable shape (i.e. sphere) regardless of the cytoskeleton structure. With protoplasts (cells with chemically digested cell wall), the growth, variation in shape and division is restricted due to the limited amount of the cell membrane which simply ruptures if the pressure inside the cell increases. The L-forms created with our bio-brick will have a mutation in the MurE gene which causes spontaneous increase in production of the IspA-like gene product which results in increased synthesis of the cell membrane. This stabilises the cell and allows for growth and division. Also, according to Svatina and Zeks(2002), different ratios between the surface area of the membrane and the cell volume yield different cell shapes due to the energy balance involved.

Approach

For the purposes of the study the complex model of the growing cell inside of the confined space can be broken down to simpler models of the system at three phases. The first phase would be a constantly growing cell, followed by a model of the cell when it touches the boundaries and therefore changes the trajectory of the membrane growth and starts adopting the shape of the boundaries.

Unlimited Cell Growth Model

Physiological basis

Since the L-Forms do not have the cell wall, the shape of the cell is determined by the laws of Newtonian physics and therefore will be as close to the shape of a sphere as possible. For the purposes of this model we will assume that the cell is spherical to start with and therefore will expand evenly in all directions. L-Forms of Bacillus subtilis made with our BioBrick will have an induced mutation in the ispA gene, which allows the cell to constitutively produce cell membrane phospholipids and therefore the surface area of the membrane will increase over time as the cell grows. Provided that the cell has a sufficient nutrient supply, phospholipids will be produced at a constant rate p, and transported to the surface of the membrane at rate T. Once at the surface, they are inserted into the membrane at rate i, which is determined by the maximum velocity of an enzyme x. These processes will result in an increase of the sphere’s radius at a rate dR(t)/dt=f(R,T,i,p). In the paper written by D. Fanelli and A.J. McKane on "Thermodynamics of vesicle growth and instability" it is explained that the growth of a vesicle is quasistatic in nature, so an equation for growth of the radius is given:

, where

However, this model seems physiologically unrealistic. As to allow for the exponential growth of the radius, the cell would have to produce lipids at an exponentially increasing rate, which is larger than the rate of increase of the radius by a factor of 2π. This seems highly unrealistic, as even at infinite concentration of source lipids inside the cell, the rate at which they are converted to phospholipids and transported to the membrane must be limited by the rate of diffusion, concentration of enzymes, and their turnover rate. while rate of diffusion are measurable constants, enzyme concentration is expected to be more or less constant to support homeostasis.

Therefore we propose a simplistic model of growth which is based on following assumptions:

• the cell is a perfect sphere
• lipids from which the membrane is made are uniformly distributed throughout the cell and the membrane has a constant thickness
• the cell is in a growth phase throughout the experiment
• rate of conversion into membrane phospholipids is constant
• the insertion of the synthesised membrane lipids is constant and spontaneous
• pressure inside and outside of the cell is quasistatic.
• there's sufficients nutrient supply to sustain growth

Limitation of Growth by a Set Boundary

As the maximum volume of the l-form cell before it divides is unknown, for the purposes of this modelling exercise we assume it's larger than the volume of the chamber and is negligible. Here we work under assumption that the cell grows uniformly in all directions as there's no evidence in the papers that would suggest otherwise. The cell membrane would not adhere to the sides of the chamber, therefore the new membrane would be redistributed towards the spare space, gradually filling the whole chamber. We used a square shape in our model.

Model

The software of choice for this particular model is MathWorks MATLAB as it allows a quick visualisation, graphing and mathematical operations.

The model was first constructed on paper, then adapted for MATLAB.

Derivation

Given that a cell is a sphere with a volume V which grows over time , it has a surface area which is dependant on time . Membrane precursor lipids are being constantly synthesised, and their total mass is equal to , where is the initial mass of membrane lipids in the cell and is a rate constant. Therefore the area of the cell will grow proportionally to he concentration of precursor lipids in the cell, and a rate constant λ.

Hence, , so , or , where

This can be solved for r(t):

, where a A is an integration constant, determined by setting , so

, so

MATLAB code

Below is the code of the model shown in MATLAB notation. The meaning of the lines is explained in green comments. Some of the parameters used in this model were taken from the ncbi database, whereas others were estimated.

%plots a time dependent circle with a growing radius which stops growing

	cellv=10^(-15); %mean cell volume
	t_end=70; %growth time
	m=100;
	lambda=10^(-10); % rate of conversion of mass into area
	M0=10^(-17); % initial mass of membrane lipids
	K=10^(-10); %rate of production of the membrane lipids
	a=10^(-4); %length of the side of a chamber
	n=100;
	tspan = linspace(0,t_end,m);
	
	r0=(cellv/pi*3/4)^(1/3); %radius of the cell
	alpha=(3*lambda)/(32*pi^2); % constant which defines the rate at which the radius of the cell grows
	
	x=zeros(n);
	y=zeros(n);
	
	x=linspace(-a,+a,n);
	y=linspace(-a,+a,n);
	
	figure(1)
	clf
	axis([-a +a -a +a])
	hold on
	
	r=r0;
	k=0;

	while r<=a %below is the code to draw a series of circles representing the cell with a growing radius over the timespan until the cell touches the boundary
		
		k=k+1;
		t=tspan(k);
		r=(5*alpha*(M0*t+0.5*K*t^2)+r0^5)^(1/5); %radius which changes with time
		c(k)= 2*pi*r; %circumference of the cell

                rectangle('Position',[-r,-r,r*2,2r], 'Curvature',[1], 'Edgecolor',[0,0,1], Facecolour, [0,0,1]) %plots the circle with the new radius

		pause(0.5)
		radius(k)=r;

	end;

        hold all
        l=0;
        circ=c(k);
        crv=1 %curvature (represented by the fracture of the rectangle's side which is curved

        while crv>0.15 %limit of curvature % the bit below codes for a cycle which draws the new cell outline
              
		k=k+1;
		t=tspan(k);
		r=(5*alpha*(M0*t+0.5*K*t^2)+r0^5)^(1/5);
		cn(k)= 2*pi*r; %new cell circumference

                l=((cn(k)-circ/4); %expansion of the line of coherence between the cell membrane and the chamber 
                crv=1-2*l/a;

                rectangle('Position',[-a,-a,a*2,2a], 'Curvature',[crv], 'Edgecolor',[0,0,1], Facecolour, [0,0,1]) %plots a rectangle with the side equal to the side of the chamber and curvature determined by the growth of the cell's circumference due to lipid synthesis and integration of lipids into the cell membrane.

		pause(0.5)
		radius(k)=r;

       end;

Simulation

2D simulation of the Cell Growth model.

References

Svetina S, Zeks B. 2002. Shape Behaviour of Lipid Vesicles as the Basis of Some Cellular Processes. The Anatomical Record 268:215-225 Shih YL, Rothfield L. 2006. The Bacterial Cytoskeleton. Microbiology and Molecular Biology Reviews p.729-754. Fanelli D 2008. Thermodynsmics of vesicle growth and instability. Phydsical Review E 78, 051406