Bootstrap Example

Overview

The goal of the model was to determine what parameters were pertinent to the speed and accuracy of our design. We wanted to our model to speed up the experimentation process in the lab and give us a good idea of the characteristics of our construct. We also wanted to determine the optimal protein to biomarker ratio to meet our time constraint of 5 mins, since our survey indicated users prefer a fast test.

Reaction

The reaction that we are trying to optimize is shown below. Each half of our split protein is attached to a part of horseradish peroxidase (HRP). The halves of our protein form a complex with the biomarker. This brings the two parts of split HRP (sHRP) closer together and primed to reconstitute. After sHRP is reconstituted, it is functional and oxidizes tetramethyl benzidine with hydrogen peroxide to form a blue product for our colorimetric output.

The reaction between the two halves of the protein complex and the biomarker, tyrosol/farnesol, were modeled as following two-step mass action kinetics (the biomarker first binds with on half to form a complex and then the complex binds with the other half). We assumed that once both halves of the protein bound together and formed a complex, the reconstitution of the split-HRP was instantaneous. Or in other words, the formation of the protein complex was the rate limiting step. The reconstituted HRP then oxidizes TMB into a blue compound for our colorimetric response. The reconstituted HRP was modeled as following Michaelis Menten kinetics.

The initial values for the parameters are listed in the table below.

Constant	Value	Source
Kd1	0.2 nM	[1]
Kd2	26 µM	[1]
Kd3	12 nM	[1]
Kd1	100 fM	[1]
Kd (avg)	6.5 µM
k_cat	435	[2]
Km	90 µM	[2]

The dissociation constants from the FRB, FKBP, and rapamycin reaction were used for the split protein reconstitution kinetics because of their similar interactions compared to our product.

The dissociation constants were then averaged to simplify our analysis later. It was determined that using the real dissociation constants versus using the average of them had very little effect on the color output over time graph. In the figure below, a simulation of the reaction was run with the real and averaged dissociation constants:

As you can see, there is little difference between the two curves, and the calculated standard square error (SSE) between them was around 1E-5. This led us to conclude that matlab cannot be as precise with such small numbers and that the use of real or averaged k_d didn’t matter that much.

Protein Concentrations

In modeling methods we wanted to implement in the future to better characterize our product, we would need to determine certain reaction constants which relies on known reagent concentrations, including protein concentrations. Since poly-histidine tags were not included in our genetic construct and comparing the colorimetric output between our sHRP and wild type HRP was not viable (since the difference in output could be attributed to different concentrations or activity levels), we had to figure out a different way to determine our protein concentrations once we lysed the cell.

Our model offered us one solution. It was observed in simulations of the reaction at fixed protein concentrations but varying biomarker concentrations, that the concentration of reconstituted sHRP at equilibrium was the greatest (around 80% of initial split protein concentration) when the protein halves and biomarker concentration were the same.

This was an important observation to make because it allows us to “titrate” our lysate to determine the protein concentration.

Sensitivity Analysis

In our sensitivity analysis, we adjusted four parameters, the dissociation constant (averaged), the Michaelis Menten constants for split-HRP (k_cat and k_m), and the protein concentration. The purpose of the sensitivity analysis was to determine which of these parameters affected the time in which all the TMB was reacted the greatest. Each parameter except one was held constant at the initial values listed in the table above. The varied parameter was swept across 7 orders of magnitude in either direction.

We can see clearly that protein concentration and k_cat affect the time that it takes TMB to fully react the greatest, changing it by around 15 and 10 orders of magnitude respectively. Each lower the time considerably as they increase. K_m and k_d on the other hand don’t change the time as much or only do once they wander out of the domain of possible values.

In conclusion, according to our sensitivity analysis, it seems it would be crucial to determine our sHRP’s k_cat value and adjust our protein concentrations in order to meet our time constraint.

Future Modeling

Now that we have determined the important parameters, we know what to focus on once our protein construct is synthesized. We can find k_cat of our sHRP by producing a Lineweaver Burk plot. Then, we can run multiple simulations at different protein to biomarker concentrations (at both infected and healthy concentration levels) to determine the optimal ratio that performs under 5 mins. but still has a large enough difference in complete reaction time to distinguish between a positive and negative test.

Reciever Operating Characteristic (ROC)

The ROC is an important property of diagnostic tests in industry, plotting the true positive rate (sensitivity) over the false-positive rate (1 – specificity). The integral of this curve gives you the accuracy of your diagnostic test.

Whilst waiting for our wet lab crew to construct the actual protein construct, we were able to model its behavior and develop a plausible ROC for our product using the protein concentration estimation method described above.

References

Banaszynski, L. A., Liu, C. W., & Wandless, T. J. (2005). Characterization of the FKBP·Rapamycin·FRB Ternary Complex. J. Am. Chem. Soc. 127(13), p. 4721.
Brenda:https://www.brenda-enzymes.org/enzyme.php?ecno=1.11.1.7