{ "cells": [ { "cell_type": "markdown", "id": "5cbc8640", "metadata": {}, "source": [ "# **Enzyme Kinetics** : \n", "#### An _accurate numerical solution_ of an enzymatic reaction `E + S <-> ES -> E + P` - compared to the **Michaelis-Menten** model approximation and to the alternative **Morrison** model.\n", "\n", "#### Same as experiments `enzyme_1_a` and `enzyme_1_b` but with _a lavish amount of Enzyme_ relative to the initial Substrate concentration. \n", "Experiment `enzyme_1_b` had a ten-fold increase in enzyme, relative to `enzyme_1_a`. Here, we do _another x10 increase!_\n", "\n", "#### Real-life kinetic parameters. \n", "\n", "Please refer to `enzyme_1_a` and `enzyme_1_b` for more details." ] }, { "cell_type": "markdown", "id": "604b150b-7812-4fd3-9403-69a06dd7e397", "metadata": {}, "source": [ "#### THE REACTION: \n", "the enzyme `Adenosine deaminase` with the substrate `2,6-Diamino-9-β-D-deoxyribofuranosyl-9-H-purine` (same as in experiment `enzyme_1`) \n", "satisfies the customary Michaelis-Menten assumptions that `k1_reverse >> k2_forward`\n", "\n", "However, the initial concentration values choosen below, do NOT satisfy the expected Michaelis-Menten assumptions that `[E] << [S]` \n", "\n", "Source of kinetic parameters: *page 16 of \"Analysis of Enzyme Reaction Kinetics, Vol. 1\", by F. Xavier Malcata, Wiley, 2023*" ] }, { "cell_type": "markdown", "id": "c123db4f-c802-47f0-a3d3-0b857314d8e5", "metadata": {}, "source": [ "### TAGS : \"uniform compartment\", \"chemistry\", \"numerical\", \"enzymes\"" ] }, { "cell_type": "code", "execution_count": 1, "id": "6e9d0902-6fc9-4692-ac39-0651d08902ca", "metadata": {}, "outputs": [], "source": [ "LAST_REVISED = \"Sep. 2, 2025\"\n", "LIFE123_VERSION = \"1.0.0rc6\" # Library version this experiment is based on" ] }, { "cell_type": "code", "execution_count": 2, "id": "1e0ae9a9-9d0c-4edf-a5f2-1c589419e6cf", "metadata": {}, "outputs": [], "source": [ "#import set_path # Using MyBinder? Uncomment this before running the next cell!" ] }, { "cell_type": "code", "execution_count": 3, "id": "a29db1c7", "metadata": { "tags": [] }, "outputs": [], "source": [ "#import sys\n", "#sys.path.append(\"C:/some_path/my_env_or_install\") # CHANGE to the folder containing your venv or libraries installation!\n", "# NOTE: If any of the imports below can't find a module, uncomment the lines above, or try: import set_path \n", "\n", "import pandas as pd\n", "\n", "from life123 import check_version, ChemData, UniformCompartment, ReactionEnzyme, PlotlyHelper" ] }, { "cell_type": "code", "execution_count": 4, "id": "af15ecf0-e083-4fef-b68e-abe794dcc86e", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "OK\n" ] } ], "source": [ "check_version(LIFE123_VERSION) # To check compatibility" ] }, { "cell_type": "code", "execution_count": null, "id": "3713fa4d-e9bb-4e33-8734-c1d4d23cf177", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "34d1cefc-f644-410a-9fe4-5204964742ac", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "91a729ee-ad2f-4771-bc84-49ecb9d6a851", "metadata": {}, "source": [ "# PART 1. Accurate numerical solution" ] }, { "cell_type": "code", "execution_count": 5, "id": "ac9eea69-174c-43e5-9eed-443cbc5e2ba7", "metadata": {}, "outputs": [], "source": [ "chem_data = ChemData(names=[\"P\", \"ES\"], plot_colors=[\"green\", \"red\"])" ] }, { "cell_type": "code", "execution_count": 6, "id": "90f9ec41-6e1b-4558-9400-47e35ece6efc", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
| \n", " | name | \n", "label | \n", "plot_color | \n", "
|---|---|---|---|
| 0 | \n", "P | \n", "P | \n", "green | \n", "
| 1 | \n", "ES | \n", "ES | \n", "red | \n", "
| 2 | \n", "Adenosine deaminase | \n", "E | \n", "violet | \n", "
| 3 | \n", "2,6-Diamino-9-β-D-deoxyribofuranosyl-9-H-purine | \n", "S | \n", "darkturquoise | \n", "
| \n", " | SYSTEM TIME | \n", "P | \n", "ES | \n", "E | \n", "S | \n", "step | \n", "caption | \n", "rxn0_rate_2 | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "100.000000 | \n", "20.000000 | \n", "\n", " | Set concentration | \n", "0.000000 | \n", "
| 1 | \n", "0.000010 | \n", "0.000000 | \n", "0.360000 | \n", "99.640000 | \n", "19.640000 | \n", "1 | \n", "1st reaction step | \n", "17.640000 | \n", "
| 2 | \n", "0.000011 | \n", "0.000018 | \n", "0.395523 | \n", "99.604477 | \n", "19.604459 | \n", "2 | \n", "\n", " | 19.380617 | \n", "
| 3 | \n", "0.000011 | \n", "0.000019 | \n", "0.397295 | \n", "99.602705 | \n", "19.602686 | \n", "3 | \n", "\n", " | 19.467446 | \n", "
| 4 | \n", "0.000011 | \n", "0.000019 | \n", "0.398181 | \n", "99.601819 | \n", "19.601800 | \n", "4 | \n", "\n", " | 19.510856 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 993 | \n", "0.047918 | \n", "17.775498 | \n", "2.101933 | \n", "97.898067 | \n", "0.122569 | \n", "993 | \n", "\n", " | 102.994725 | \n", "
| 994 | \n", "0.048397 | \n", "17.824842 | \n", "2.055364 | \n", "97.944636 | \n", "0.119793 | \n", "994 | \n", "\n", " | 100.712859 | \n", "
| 995 | \n", "0.048881 | \n", "17.873576 | \n", "2.009370 | \n", "97.990630 | \n", "0.117055 | \n", "995 | \n", "\n", " | 98.459109 | \n", "
| 996 | \n", "0.049370 | \n", "17.921695 | \n", "1.963952 | \n", "98.036048 | \n", "0.114353 | \n", "996 | \n", "\n", " | 96.233637 | \n", "
| 997 | \n", "0.049863 | \n", "17.969197 | \n", "1.919114 | \n", "98.080886 | \n", "0.111689 | \n", "997 | \n", "\n", " | 94.036596 | \n", "
998 rows × 8 columns
\n", "| \n", " | SYSTEM TIME | \n", "P | \n", "ES | \n", "E | \n", "S | \n", "step | \n", "caption | \n", "rxn0_rate_2 | \n", "Michaelis_rate | \n", "
|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "100.000000 | \n", "20.000000 | \n", "\n", " | Set concentration | \n", "0.000000 | \n", "3465.618861 | \n", "
| 1 | \n", "0.000010 | \n", "0.000000 | \n", "0.360000 | \n", "99.640000 | \n", "19.640000 | \n", "1 | \n", "1st reaction step | \n", "17.640000 | \n", "3447.122503 | \n", "
| 2 | \n", "0.000011 | \n", "0.000018 | \n", "0.395523 | \n", "99.604477 | \n", "19.604459 | \n", "2 | \n", "\n", " | 19.380617 | \n", "3445.270565 | \n", "
| 3 | \n", "0.000011 | \n", "0.000019 | \n", "0.397295 | \n", "99.602705 | \n", "19.602686 | \n", "3 | \n", "\n", " | 19.467446 | \n", "3445.178054 | \n", "
| 4 | \n", "0.000011 | \n", "0.000019 | \n", "0.398181 | \n", "99.601819 | \n", "19.601800 | \n", "4 | \n", "\n", " | 19.510856 | \n", "3445.131800 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 993 | \n", "0.047918 | \n", "17.775498 | \n", "2.101933 | \n", "97.898067 | \n", "0.122569 | \n", "993 | \n", "\n", " | 102.994725 | \n", "71.495562 | \n", "
| 994 | \n", "0.048397 | \n", "17.824842 | \n", "2.055364 | \n", "97.944636 | \n", "0.119793 | \n", "994 | \n", "\n", " | 100.712859 | \n", "69.899661 | \n", "
| 995 | \n", "0.048881 | \n", "17.873576 | \n", "2.009370 | \n", "97.990630 | \n", "0.117055 | \n", "995 | \n", "\n", " | 98.459109 | \n", "68.323957 | \n", "
| 996 | \n", "0.049370 | \n", "17.921695 | \n", "1.963952 | \n", "98.036048 | \n", "0.114353 | \n", "996 | \n", "\n", " | 96.233637 | \n", "66.768545 | \n", "
| 997 | \n", "0.049863 | \n", "17.969197 | \n", "1.919114 | \n", "98.080886 | \n", "0.111689 | \n", "997 | \n", "\n", " | 94.036596 | \n", "65.233510 | \n", "
998 rows × 9 columns
\n", "| \n", " | SYSTEM TIME | \n", "P | \n", "ES | \n", "E | \n", "S | \n", "step | \n", "caption | \n", "rxn0_rate_2 | \n", "Michaelis_rate | \n", "Morrison_rate | \n", "
|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "0.000000 | \n", "0.000000 | \n", "0.000000 | \n", "100.000000 | \n", "20.000000 | \n", "\n", " | Set concentration | \n", "0.000000 | \n", "3465.618861 | \n", "889.979162 | \n", "
| 1 | \n", "0.000010 | \n", "0.000000 | \n", "0.360000 | \n", "99.640000 | \n", "19.640000 | \n", "1 | \n", "1st reaction step | \n", "17.640000 | \n", "3447.122503 | \n", "889.979162 | \n", "
| 2 | \n", "0.000011 | \n", "0.000018 | \n", "0.395523 | \n", "99.604477 | \n", "19.604459 | \n", "2 | \n", "\n", " | 19.380617 | \n", "3445.270565 | \n", "889.978385 | \n", "
| 3 | \n", "0.000011 | \n", "0.000019 | \n", "0.397295 | \n", "99.602705 | \n", "19.602686 | \n", "3 | \n", "\n", " | 19.467446 | \n", "3445.178054 | \n", "889.978342 | \n", "
| 4 | \n", "0.000011 | \n", "0.000019 | \n", "0.398181 | \n", "99.601819 | \n", "19.601800 | \n", "4 | \n", "\n", " | 19.510856 | \n", "3445.131800 | \n", "889.978321 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 993 | \n", "0.047918 | \n", "17.775498 | \n", "2.101933 | \n", "97.898067 | \n", "0.122569 | \n", "993 | \n", "\n", " | 102.994725 | \n", "71.495562 | \n", "100.506658 | \n", "
| 994 | \n", "0.048397 | \n", "17.824842 | \n", "2.055364 | \n", "97.944636 | \n", "0.119793 | \n", "994 | \n", "\n", " | 100.712859 | \n", "69.899661 | \n", "98.280755 | \n", "
| 995 | \n", "0.048881 | \n", "17.873576 | \n", "2.009370 | \n", "97.990630 | \n", "0.117055 | \n", "995 | \n", "\n", " | 98.459109 | \n", "68.323957 | \n", "96.082240 | \n", "
| 996 | \n", "0.049370 | \n", "17.921695 | \n", "1.963952 | \n", "98.036048 | \n", "0.114353 | \n", "996 | \n", "\n", " | 96.233637 | \n", "66.768545 | \n", "93.911274 | \n", "
| 997 | \n", "0.049863 | \n", "17.969197 | \n", "1.919114 | \n", "98.080886 | \n", "0.111689 | \n", "997 | \n", "\n", " | 94.036596 | \n", "65.233510 | \n", "91.768007 | \n", "
998 rows × 10 columns
\n", "