Difference between revisions of "Team:Aalto-Helsinki/Modeling micelle"

m (added pic of overlapping)
m
Line 103: Line 103:
 
If we assume that our enzymes have a density that is common to enzymes, our task becomes easier. Based on this paper, the average density of proteins is 1.37 g/ml. Because we want to calculate the volume, and in effect, their radius, we invert this value, thus getting 0.73 ml/g.</p>
 
If we assume that our enzymes have a density that is common to enzymes, our task becomes easier. Based on this paper, the average density of proteins is 1.37 g/ml. Because we want to calculate the volume, and in effect, their radius, we invert this value, thus getting 0.73 ml/g.</p>
  
<p>Using some clever calculation, we get that the relationship between mass and volume for proteins is \[V(nm^3)=\frac{0.73\tfrac{nm^3}{g}10^{21}\tfrac{nm^3}{cm^3}}{6.023\cdot 10^-3 \tfrac{nm^3}{Da}} M(Da).\] We can calculate the radius by \[r_{min} = \left( \frac{3V}{4\pi} \right)^{1/3}\] if we know the volume of a sphere, so we get that \( R_{min}=0.066\cdot M(Da) \).</p>
+
<p>Using some clever calculation, we get that the relationship between mass and volume for proteins is  
 +
\[V(nm^3)=\frac{0.73\tfrac{nm^3}{g}10^{21}\tfrac{nm^3}{cm^3}}{6.023\cdot 10^-3 \tfrac{nm^3}{Da}} M(Da).\]  
 +
We can calculate the radius by  
 +
\[r_{min} = \left( \frac{3V}{4\pi} \right)^{1/3}\]  
 +
if we know the volume of a sphere, so we get that \( R_{min}=0.066\cdot M(Da) \).</p>
  
 
<p>The mass of CAR is <a href = "http://www.uniprot.org/uniprot/B2HN69">127 797 DA</a> and the mass of ADO is <a href="http://www.rcsb.org/pdb/explore/explore.do?structureId=4KVS">27 569.15 Da</a>. Now the radius would be 3,5 nm for CAR and 2 nm for ADO. The mass of Gfp is 26 890 Da, which makes its radius roughly 2 nm.</p>
 
<p>The mass of CAR is <a href = "http://www.uniprot.org/uniprot/B2HN69">127 797 DA</a> and the mass of ADO is <a href="http://www.rcsb.org/pdb/explore/explore.do?structureId=4KVS">27 569.15 Da</a>. Now the radius would be 3,5 nm for CAR and 2 nm for ADO. The mass of Gfp is 26 890 Da, which makes its radius roughly 2 nm.</p>
Line 125: Line 129:
 
<h2 id="adocar">Calculations for Ado and Car</h2>
 
<h2 id="adocar">Calculations for Ado and Car</h2>
  
<p>We can estimate how many amphiphilic proteins we can theoretically fit in one micelle by calculating how big solid angles they take with attached enzymes. The easiest way to estimate the solid angles is to think the amphiphilic proteins linked with enzymes as cones. We can calculate the solid angle \( \Omega \) for these by \[ \Omega = 2\pi \left( 1-\cos(\theta) \right), \] where \( \theta \) is half of the apex angle. So for CAR we get \[ \Omega_{cone-CAR} = 2\pi \left( 1-\cos\left( \arctan\left(\frac{3.5}{14.1}\right)\right) \right) \approx 0.185 \text{ rad} \] and for ADO \[  \Omega_{cone-ADO} = 2\pi \left( 1-\cos\left( \arctan\left(\frac{2}{9.8}\right)\right) \right) \approx 0.127 \text{ rad}.\]</p>
+
<p>We can estimate how many amphiphilic proteins we can theoretically fit in one micelle by calculating how big solid angles they take with attached enzymes. The easiest way to estimate the solid angles is to think the amphiphilic proteins linked with enzymes as cones. We can calculate the solid angle \( \Omega \) for these by  
 +
\[ \Omega = 2\pi \left( 1-\cos(\theta) \right), \]  
 +
where \( \theta \) is half of the apex angle. So for CAR we get  
 +
\[ \Omega_{cone\text{-}CAR} = 2\pi \left( 1-\cos\left( \arctan\left(\frac{3.5}{14.1}\right)\right) \right) \approx 0.185 \text{ rad} \]  
 +
and for ADO  
 +
\[  \Omega_{cone\text{-}ADO} = 2\pi \left( 1-\cos\left( \arctan\left(\frac{2}{9.8}\right)\right) \right) \approx 0.127 \text{ rad}.\]</p>
  
 
<p style="color:gray">--picture of this cone-like structure? is it needed or can this be understood without it?--</p>
 
<p style="color:gray">--picture of this cone-like structure? is it needed or can this be understood without it?--</p>
Line 135: Line 144:
 
<p style="color:gray">--picture of pyramid structure? is it needed or can this be understood without it?--</p>
 
<p style="color:gray">--picture of pyramid structure? is it needed or can this be understood without it?--</p>
  
<p>The solid angle \( \Omega\) for this kind of structure can be calculated by \[\Omega = 4\arcsin\left( \sin\left(\theta\right) ^2 \right),\] where \( \theta\) is again half of the apex angle. This yields us \[\Omega_{pyramid-CAR} = 4\arcsin\left( \sin\left(\arctan\left( \frac{3.5}{14.1} \right) \right) ^2 \right) \approx 0.232 \text{ rad}\] and \[\Omega_{pyramid-ADO} = 4\arcsin\left( \sin\left(\arctan\left( \frac{2}{9.8} \right) \right) ^2 \right) \approx 0.16 \text{ rad}.\] </p>
+
<p>The solid angle \( \Omega\) for this kind of structure can be calculated by  
 +
\[\Omega = 4\arcsin\left( \sin\left(\theta\right) ^2 \right),\] where \( \theta\) is again half of the apex angle. This yields us \[\Omega_{pyramid\text{-}CAR} = 4\arcsin\left( \sin\left(\arctan\left( \frac{3.5}{14.1} \right) \right) ^2 \right) \approx 0.232 \text{ rad}\]  
 +
and  
 +
\[\Omega_{pyramid\text{-}ADO} = 4\arcsin\left( \sin\left(\arctan\left( \frac{2}{9.8} \right) \right) ^2 \right) \approx 0.16 \text{ rad}.\] </p>
  
 
<p>By this method of calculation we could get at most 32 of both fusion proteins in one micelle.</p>
 
<p>By this method of calculation we could get at most 32 of both fusion proteins in one micelle.</p>
Line 149: Line 161:
 
</figure>
 
</figure>
  
<p>We can think that this structure consists of a single cone whose centre is the centre of CAR fusion protein and the side goes along ADO fusion protein. We can further take some of the empty areas into account by thinking this as pyramid instead of cone. This yields us \[\Omega_{CAR\&ADO} = 4\arcsin\left( \sin\left( \arccos\left( \frac{14.1^2+9.8^2-5.5^2}{2\cdot14.1\cdot 9.8} \right) \right)  ^2 \right) \approx 0.3336 \text{ rad}.\]</p>
+
<p>We can think that this structure consists of a single cone whose centre is the centre of CAR fusion protein and the side goes along ADO fusion protein. We can further take some of the empty areas into account by thinking this as pyramid instead of cone. This yields us  
 +
\[\Omega_{CAR\&ADO} = 4\arcsin\left( \sin\left( \arccos\left( \frac{14.1^2+9.8^2-5.5^2}{2\cdot14.1\cdot 9.8} \right) \right)  ^2 \right) \approx 0.3336 \text{ rad}.\]</p>
  
 
<p>This means that there fits about 37 of these pyramid structures in one micelle, so 37 CAR-enzymes. For ADO we can approximate that there is about twice as many of them than CAR fusion proteins (this is justified in infinite field so we approximate with it here), so the amount of ADO would be 74 and the whole amount of fusion proteins in this micelle 111. Since there is probably even more efficient way of packing these proteins in one micelle, the real upper bound might be even larger.</p>
 
<p>This means that there fits about 37 of these pyramid structures in one micelle, so 37 CAR-enzymes. For ADO we can approximate that there is about twice as many of them than CAR fusion proteins (this is justified in infinite field so we approximate with it here), so the amount of ADO would be 74 and the whole amount of fusion proteins in this micelle 111. Since there is probably even more efficient way of packing these proteins in one micelle, the real upper bound might be even larger.</p>
Line 155: Line 168:
 
<h2 id="gfp">Calculations for Gfp</h2>
 
<h2 id="gfp">Calculations for Gfp</h2>
  
<p>For comparison we calculated how big micelles we could possibly get with green fluorescent protein. Since the Gfp is same size as Ado, we can use values from previous calculations. With cone-approximation we get \[\frac{4\pi}{ \Omega_{cone-ADO}} \approx 98\] of these fusion proteins in one micelle, and with pyramid-approximation \[\frac{4\pi}{ \Omega_{pyramid-ADO}} \approx 78\] fusion proteins. The real value is thus probably somewhere between them.
+
<p>For comparison we calculated how big micelles we could possibly get with green fluorescent protein. Since the Gfp is same size as Ado, we can use values from previous calculations. With cone-approximation we get  
 +
\[\frac{4\pi}{ \Omega_{cone\text{-}ADO}} \approx 98\]  
 +
of these fusion proteins in one micelle, and with pyramid-approximation  
 +
\[\frac{4\pi}{ \Omega_{pyramid\text{-}ADO}} \approx 78\] fusion proteins. The real value is thus probably somewhere between them.
 
</p>
 
</p>
  

Revision as of 07:13, 5 August 2015

Introduction

--Picture of the pathway here, CAR, ADO and butyraldehyde highlighted to clarify what we are talking about.--

The product of second to last enzyme of our pathway, butyraldehyde, is toxic to the cell. Because of that and about 15 naturally occurring butyraldehyde-eating enzymes in the cell it is essential for the propane production that Butyraldehyde goes swiftly to the enzyme we want it to go, ADO. As the solution to this our team wanted to put CAR and ADO close together in a micelle so that butyraldehyde would go with more probability to ADO than to any other enzyme.

We have made a model of effectiveness of having enzymes close together, but our team also wanted to know if the micelle structure was possible at the first place. We know (references as links for this statement!) that it is possible to form the micelle without any proteins at the end and with green fluorecent protein (Gfp), but could CAR and ADO be part of this kind of structure?

Geometrical approach

2d simplification of micelle

Micelle structure

The micelle is formed by amphiphilic proteins that have both hydrophilic and hydrophobic parts. At the end of hydrophilic part there is short protein, a linker that attaches CAR or ADO to the amphiphilic part.

Image form structure
prediction software

Amphiphilic proteins are 10 nm long, 5 nm for both hydrophilic and hydrophobic parts. (Here where we got amphiphilic proteins sizes.) The linker (here link for more info about this. Structure and such, does lab have that somewhere?) consists of eight amino acids, for which the maximum lengths are 3,8Å. From this we can calculate that at most the length of one linker is 2,8 nm. If the linker would form α-helical structure, then the length for one peptide would be about 1,5 Å so the one linker would be 1,2 nm long. (we need some source for the Å-lengths) However, we can estimate that the linkers are straight, since when running the structure in peptide structure prediction software doesn't yield strong folding or helical structure. CAR uses two of these linkers and ADO one.

One problem we are facing here is that we need some sort of approximations for the enzymes’ radii. Since we don’t know the exact three-dimensional structure of the proteins, we approximated the enzymes as perfect spheres. If we assume that our enzymes have a density that is common to enzymes, our task becomes easier. Based on this paper, the average density of proteins is 1.37 g/ml. Because we want to calculate the volume, and in effect, their radius, we invert this value, thus getting 0.73 ml/g.

Using some clever calculation, we get that the relationship between mass and volume for proteins is \[V(nm^3)=\frac{0.73\tfrac{nm^3}{g}10^{21}\tfrac{nm^3}{cm^3}}{6.023\cdot 10^-3 \tfrac{nm^3}{Da}} M(Da).\] We can calculate the radius by \[r_{min} = \left( \frac{3V}{4\pi} \right)^{1/3}\] if we know the volume of a sphere, so we get that \( R_{min}=0.066\cdot M(Da) \).

The mass of CAR is 127 797 DA and the mass of ADO is 27 569.15 Da. Now the radius would be 3,5 nm for CAR and 2 nm for ADO. The mass of Gfp is 26 890 Da, which makes its radius roughly 2 nm.

Car fusion protein
Ado fusion protein

Calculations for Ado and Car

We can estimate how many amphiphilic proteins we can theoretically fit in one micelle by calculating how big solid angles they take with attached enzymes. The easiest way to estimate the solid angles is to think the amphiphilic proteins linked with enzymes as cones. We can calculate the solid angle \( \Omega \) for these by \[ \Omega = 2\pi \left( 1-\cos(\theta) \right), \] where \( \theta \) is half of the apex angle. So for CAR we get \[ \Omega_{cone\text{-}CAR} = 2\pi \left( 1-\cos\left( \arctan\left(\frac{3.5}{14.1}\right)\right) \right) \approx 0.185 \text{ rad} \] and for ADO \[ \Omega_{cone\text{-}ADO} = 2\pi \left( 1-\cos\left( \arctan\left(\frac{2}{9.8}\right)\right) \right) \approx 0.127 \text{ rad}.\]

--picture of this cone-like structure? is it needed or can this be understood without it?--

This means that there would be at most 40 of both ADO and CAR fusion proteins in one micelle by this method of calculation.

However, when approximating enzymes with spheres it is not possible for them to fill the whole surface of the micelle; there will always be gaps. This is why it might be better to approximate the solid angle these complexes of amphiphilic proteins and enzymes by using pyramides instead of cones.

--picture of pyramid structure? is it needed or can this be understood without it?--

The solid angle \( \Omega\) for this kind of structure can be calculated by \[\Omega = 4\arcsin\left( \sin\left(\theta\right) ^2 \right),\] where \( \theta\) is again half of the apex angle. This yields us \[\Omega_{pyramid\text{-}CAR} = 4\arcsin\left( \sin\left(\arctan\left( \frac{3.5}{14.1} \right) \right) ^2 \right) \approx 0.232 \text{ rad}\] and \[\Omega_{pyramid\text{-}ADO} = 4\arcsin\left( \sin\left(\arctan\left( \frac{2}{9.8} \right) \right) ^2 \right) \approx 0.16 \text{ rad}.\]

By this method of calculation we could get at most 32 of both fusion proteins in one micelle.

Above isn’t a perfect arrangement of these fusion proteins either, but the problem is too hard in ball surface with two different sizes of enzymes. The real maximum value if we think the problem this way is somewhere between the ones obtained, so somewhere between 64 and 80 fusion proteins in a micelle.

The previous calculations have not taken into account that CAR and ADO might overlap because CAR has two linkers when ADO has just one. We don’t know the ideal structure of the overlapping, but we can estimate it by the structure shown below.

Estimation of the overlapping micelle structure.

We can think that this structure consists of a single cone whose centre is the centre of CAR fusion protein and the side goes along ADO fusion protein. We can further take some of the empty areas into account by thinking this as pyramid instead of cone. This yields us \[\Omega_{CAR\&ADO} = 4\arcsin\left( \sin\left( \arccos\left( \frac{14.1^2+9.8^2-5.5^2}{2\cdot14.1\cdot 9.8} \right) \right) ^2 \right) \approx 0.3336 \text{ rad}.\]

This means that there fits about 37 of these pyramid structures in one micelle, so 37 CAR-enzymes. For ADO we can approximate that there is about twice as many of them than CAR fusion proteins (this is justified in infinite field so we approximate with it here), so the amount of ADO would be 74 and the whole amount of fusion proteins in this micelle 111. Since there is probably even more efficient way of packing these proteins in one micelle, the real upper bound might be even larger.

Calculations for Gfp

For comparison we calculated how big micelles we could possibly get with green fluorescent protein. Since the Gfp is same size as Ado, we can use values from previous calculations. With cone-approximation we get \[\frac{4\pi}{ \Omega_{cone\text{-}ADO}} \approx 98\] of these fusion proteins in one micelle, and with pyramid-approximation \[\frac{4\pi}{ \Omega_{pyramid\text{-}ADO}} \approx 78\] fusion proteins. The real value is thus probably somewhere between them.

Discussion

The goal of this modeling approach was to understand if it was possible to form micelles which have CAR or ADO at the end of the amphiphilic proteins. The main thing was to prove that the proteins aren’t too big to have impact to micelle modeling, and since we already knew that this arrangement works with green fluorescent protein it was only natural to compare these two. Because based on our calculations the green fluorescent protein micelles have upper bound of building blocks somewhere between 78 and 98 and the micelles with CAR and ADO somewhere between 64 and 111 we can say that geometrically it is possible for CAR and ADO to be part of a micelle structure.

There are some assumptions of the model that might have some effect on its accuracy. We have assumed that the amphiphilic protein could be approximated by just its length, and it has no width that could have any effect on our calculations. It is also to be noted that we didn’t even aim to be accurate in assembly of ADO and CAR in ball surface. The best possible formation is very hard to find in this situation and there wasn’t any need to be that accurate in our calculations. Further, since we don't know what shapes the enzymes are, we have estimated them as spheres.

Even though our model seems to prove that the formation of these micelles is possible, there are lots of things we couldn’t take into account that might have effects on micelle formation and make it impossible.