Microcompartments are proteinaceous capsules that have only recently been discovered to exist in a wide range of bacteria. They contain enzymes required for a particular metabolic process. The carboxysome is a particularly well-studied example of a specialised microcompartment, and is shown here as a schematic diagram (S. Frank et al, 2013).

The Pdu microcompartment that we are using in our project is composed of several proteins encoded in an operon consisting of the genes pduA, -B, -T, -U, -N, -J, -K. These are expressed in various stoichiometries to form different polyhedral shapes with a diameter of up to 100-150 nm. The faces of the polyhedron are formed by the hexagonal shell proteins PduA, PduB, and PduJ, while PduN is thought to form the vertices (Joshua B. Parsons et al, 2010).

We are expressing the pdu-ABTUNJK codon in E. coli in the pUNI vector, shown below, which we received from the University of Dundee iGEM team.









Additionally, we are expressing pduABUTNJK in P. putida, which to our knowledge has not been done before. For this purpose, we transferred the ABTUNJK sequence into the pBBR1MCS vector, shown below:



Oxford iGEM 2014

Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information

Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information Information

The aim of this model was to predict a theoretical maximum number of enzyme molecules that can be packed into a single microcompartment. To get a first estimate, without taking into consideration whether this volume of protein would interrupt the biological processes in the cell, we approached this problem volumetrically.

Due to the complexity of the enzyme movements and their interactions, I simplified their structures by approximating them as ellipsoids with axes lengths calculated through modelling the monomers and predicting the structures of the FdhA tetramer and DcmA hexamer respectively.





Oxford iGEM 2014

The ellipsoid packing problem

Once treated as ellipsoids, the problem was then reduced to the classical ‘sand packing’ problem. Because the dimensions of these proteins was substantially smaller than the icosahedron (by approximately a factor of 20 in every dimension), I assumed that the geometry of the container i.e. the microcompartment, was not significant.

Another assumption made in these calculations was that the enzymes could be treated as homogenous. They are of very similar dimensions, varying by no more than 20-30% on any axis, and also have very similar sphericities- the key variable in determining the packing efficiency of the molecules. Sphericity is defined as:

φ = sphericity

V_p = volume of ellipsoid (nm^3)

A_p = surface area of ellipsoid (nm^2)

For ellipsoids, a surface area approximation was used:

a = length of axis 1 (nm)

b = length of axis 2 (nm)

c = length of axis 3 (nm)

After calculating the sphericities of the enzymes, the porosity of the system could then be determined through empirical data from literature. Because the DcmA and FdhA sphericities were very similar (0.953 and 0.981 respectively), we considered the system to be composed of a homogenous spheroid species of porosity 0.973 i.e. the weighted average of the two species.

The image in the top right is a visualization of the ellipsoids packing in the microcompartment

The graph on the right shows the empirical ellipsoid packing efficiencies as defined by sphericity

Results

Because enzyme-enzyme interactions are very difficult to predict, and it is difficult to assign an analogous friction factor to their movements, I took the average of the predicted porosities across a range of different friction factors and alternative methods giving a predicted porosity of 37.3%.

Combining this information with the relative expression rates of the two enzymes, we predict that:





We should note that these numbers are a first approximation and do not take into consideration the strain placed on the cell involved in expressing this amount of protein. This model is based entirely on volumetric and geometric laws rather than biological ones and should be considered a first assumption and theoretical maximum given optimum conditions.