An Insulin Model. James K. Peterson. June 14, Department of Biological Sciences and Department of Mathematical Sciences Clemson University
|
|
- Mervyn Collins
- 6 years ago
- Views:
Transcription
1 An Insulin Model James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University June 14, 2017
2 Outline 1 The Background for the Diabetes Model 2 The Diabetes Detection Model 3 Fitting the Data 4 Gradient Descent Implementation 5 Adding Line Search
3 Abstract This lecture discusses a model of diabetes detection.
4 The Background for the Diabetes Model We are now going to discuss a very nice model of diabetes detection which was presented in the classic text on applied mathematical modeling by Braun, Differential Equations and Their Applications.
5 The Background for the Diabetes Model We are now going to discuss a very nice model of diabetes detection which was presented in the classic text on applied mathematical modeling by Braun, Differential Equations and Their Applications. In diabetes there is too much sugar in the blood and the urine. This is a metabolic disease and if a person has it, they are not able to use up all the sugars, starches and various carbohydrates because they don t have enough insulin.
6 The Background for the Diabetes Model We are now going to discuss a very nice model of diabetes detection which was presented in the classic text on applied mathematical modeling by Braun, Differential Equations and Their Applications. In diabetes there is too much sugar in the blood and the urine. This is a metabolic disease and if a person has it, they are not able to use up all the sugars, starches and various carbohydrates because they don t have enough insulin. Diabetes can be diagnosed by a glucose tolerance test (GTT). If you are given this test, you do an overnight fast and then you are given a large dose of sugar in a form that appears in the bloodstream. This sugar is called glucose. Measurements are made over about five hours or so of the concentration of glucose in the blood. These measurements are then used in the diagnosis of diabetes. It has always been difficult to interpret these results as a means of diagnosing whether a person has diabetes or not. Hence, different physicians interpreting the same data can come up with a different diagnosis, which is a pretty unacceptable state of affairs!
7 The Background for the Diabetes Model We are now going to discuss a criterion developed in the 1960 s by doctors at the Mayo Clinic and the University of Minnesota that was fairly reliable. It showcases a lot of our modeling in this course and will give you another example of how we use our tools. We start with a simple model of the blood glucose regulatory system.
8 The Background for the Diabetes Model We are now going to discuss a criterion developed in the 1960 s by doctors at the Mayo Clinic and the University of Minnesota that was fairly reliable. It showcases a lot of our modeling in this course and will give you another example of how we use our tools. We start with a simple model of the blood glucose regulatory system. Glucose plays an important role in vertebrate metabolism because it is a source of energy. For each person, there is an optimal blood glucose concentration and large deviations from this leads to severe problems including death.
9 The Background for the Diabetes Model We are now going to discuss a criterion developed in the 1960 s by doctors at the Mayo Clinic and the University of Minnesota that was fairly reliable. It showcases a lot of our modeling in this course and will give you another example of how we use our tools. We start with a simple model of the blood glucose regulatory system. Glucose plays an important role in vertebrate metabolism because it is a source of energy. For each person, there is an optimal blood glucose concentration and large deviations from this leads to severe problems including death. Blood glucose levels are autoregulated via standard forward and backward interactions like we see in many biological systems. An example is the signal that is used to activate the creation of a protein which we discussed earlier.
10 The Background for the Diabetes Model The signaling molecules are typically either bound to another molecule in the cell or are free. The equilibrium concentration of free signal is due to the fact that the rate at which signaling molecules bind equals the rate at which they split apart from their binding substrate. When an external message comes into the cell called a trigger, it induces a change in this careful balance which temporarily upgrades or degrades the equilibrium signal concentration. This then influence the protein concentration rate.
11 The Background for the Diabetes Model The signaling molecules are typically either bound to another molecule in the cell or are free. The equilibrium concentration of free signal is due to the fact that the rate at which signaling molecules bind equals the rate at which they split apart from their binding substrate. When an external message comes into the cell called a trigger, it induces a change in this careful balance which temporarily upgrades or degrades the equilibrium signal concentration. This then influence the protein concentration rate. Blood glucose concentrations work like this too, although the details differ. The blood glucose concentration is influenced by a variety of signaling molecules just like the protein creation rates can be.
12 The Background for the Diabetes Model The signaling molecules are typically either bound to another molecule in the cell or are free. The equilibrium concentration of free signal is due to the fact that the rate at which signaling molecules bind equals the rate at which they split apart from their binding substrate. When an external message comes into the cell called a trigger, it induces a change in this careful balance which temporarily upgrades or degrades the equilibrium signal concentration. This then influence the protein concentration rate. Blood glucose concentrations work like this too, although the details differ. The blood glucose concentration is influenced by a variety of signaling molecules just like the protein creation rates can be. The hormone that decreases blood glucose concentration is insulin.
13 The Background for the Diabetes Model Insulin is a hormone secreted by the β cells of the pancreas. After we eat carbohydrates, our gastrointestinal tract sends a signal to the pancreas to secrete insulin. Also, the glucose in our blood directly stimulates the β cells to secrete insulin. We think insulin helps cells pull in the glucose needed for metabolic activity by attaching itself to membrane walls that are normally impenetrable.
14 The Background for the Diabetes Model Insulin is a hormone secreted by the β cells of the pancreas. After we eat carbohydrates, our gastrointestinal tract sends a signal to the pancreas to secrete insulin. Also, the glucose in our blood directly stimulates the β cells to secrete insulin. We think insulin helps cells pull in the glucose needed for metabolic activity by attaching itself to membrane walls that are normally impenetrable. This attachment increases the ability of glucose to pass through to the inside of the cell where it can be used as fuel. So, if there is not enough enough insulin, cells don t have enough energy for their needs.
15 The Background for the Diabetes Model Insulin is a hormone secreted by the β cells of the pancreas. After we eat carbohydrates, our gastrointestinal tract sends a signal to the pancreas to secrete insulin. Also, the glucose in our blood directly stimulates the β cells to secrete insulin. We think insulin helps cells pull in the glucose needed for metabolic activity by attaching itself to membrane walls that are normally impenetrable. This attachment increases the ability of glucose to pass through to the inside of the cell where it can be used as fuel. So, if there is not enough enough insulin, cells don t have enough energy for their needs. The other hormones we will focus on all tend to change blood glucose concentrations also.
16 The Background for the Diabetes Model Glucagon is a hormone secreted by the α cells of the pancreas. Excess glucose is stored in the liver in the form of Glycogen. There is the usual equilibrium amount of storage caused by the rate of glycogen formation being equal to the rate of the reverse reaction that moves glycogen back to glucose. Hence the glycogen serves as a reservoir for glucose and when the body needs glucose, the rate balance is tipped towards conversion back to glucose to release needed glucose to the cells. The hormone glucagon increases the rate of the reaction that converts glycogen back to glucose and so serves an important regulatory function. Hypoglycemia (low blood sugar) and fasting tend to increase the secretion of the hormone glucagon. On the other hand, if the blood glucose levels increase, this tends to suppress glucagon secretion; i.e. we have another back and forth regulatory tool.
17 The Background for the Diabetes Model Epinephrine also called adrenalin is a hormone secreted by the adrenal medulla. It is part of an emergency mechanism to quickly increase the blood glucose concentration in times of extremely low blood sugar levels. Hence, epinephrine also increases the rate at which glycogen converts to glucose. It also directly inhibits how much glucose is able to be pulled into muscle tissue because muscles use a lot of energy and this energy is needed elsewhere more urgently. It also acts on the pancreas directly to inhibit insulin production which keeps glucose in the blood. There is also another way to increase glucose by converting lactate into glucose in the liver. Epinephrine increases this rate also so the liver can pump this extra glucose back into the blood stream.
18 The Background for the Diabetes Model Glucocorticoids are hormones like cortisol which are secreted by the adrenal cortex which influence how carbohydrates are metabolized which is turn increase glucose if the the metabolic rate goes up.
19 The Background for the Diabetes Model Glucocorticoids are hormones like cortisol which are secreted by the adrenal cortex which influence how carbohydrates are metabolized which is turn increase glucose if the the metabolic rate goes up. Thyroxin is a hormone secreted by the thyroid gland and it helps the liver form glucose from sources which are not carbohydrates such as glycerol, lactate and amino acids. So another way to up glucose!
20 The Background for the Diabetes Model Glucocorticoids are hormones like cortisol which are secreted by the adrenal cortex which influence how carbohydrates are metabolized which is turn increase glucose if the the metabolic rate goes up. Thyroxin is a hormone secreted by the thyroid gland and it helps the liver form glucose from sources which are not carbohydrates such as glycerol, lactate and amino acids. So another way to up glucose! Somatotrophin is called the growth hormone and it is secreted by the anterior pituitary gland. This hormone directly affect blood glucose levels (i.e. an increase in Somatotrophin increases blood glucose levels and vice versa) but it also inhibits the effect of insulin on muscle and fat cell s permeability which diminishes insulin s ability to help those cells pull glucose out of the blood stream. These actions can therefore increase blood glucose levels.
21 The Background for the Diabetes Model Now net hormone concentration is the sum of insulin plus the others. Let H denote this net hormone concentration.
22 The Background for the Diabetes Model Now net hormone concentration is the sum of insulin plus the others. Let H denote this net hormone concentration. At normal conditions, call this concentration H 0. There have been studies performed that show that under close to normal conditions, the interaction of the one hormone insulin with blood glucose completely dominates the net hormonal activity. That is normal blood sugar levels primarily depend on insulin-glucose interactions.
23 The Background for the Diabetes Model Now net hormone concentration is the sum of insulin plus the others. Let H denote this net hormone concentration. At normal conditions, call this concentration H 0. There have been studies performed that show that under close to normal conditions, the interaction of the one hormone insulin with blood glucose completely dominates the net hormonal activity. That is normal blood sugar levels primarily depend on insulin-glucose interactions. So if insulin increases from normal levels, it increases net hormonal concentration to H 0 + H and decreases glucose blood concentration.
24 The Background for the Diabetes Model On the other hand, if other hormones such as cortisol increased from base levels, this will make blood glucose levels go up.
25 The Background for the Diabetes Model On the other hand, if other hormones such as cortisol increased from base levels, this will make blood glucose levels go up. Since insulin dominates all activity at normal conditions, we can think of this increase in cortisol as a decrease in insulin with a resulting drop in blood glucose levels. A decrease in insulin from normal levels corresponds to a drop in net hormone concentration to H 0 H.
26 The Background for the Diabetes Model On the other hand, if other hormones such as cortisol increased from base levels, this will make blood glucose levels go up. Since insulin dominates all activity at normal conditions, we can think of this increase in cortisol as a decrease in insulin with a resulting drop in blood glucose levels. A decrease in insulin from normal levels corresponds to a drop in net hormone concentration to H 0 H. Now let G denote blood glucose level. Hence, in our model an increase in H means a drop in G and a decrease in H means an increase in G! Note our lumping of all the hormone activity into a single net activity is very much like how we modeled food fish and predator fish in the predator prey model.
27 The Background for the Diabetes Model The idea of our model for diagnosing diabetes from the GTT is to find a simple dynamical model of this complicated blood glucose regulatory system in which the values of two parameters would give a nice criterion for distinguishing normal individuals from those with mild diabetes or those who are pre diabetic. Here is what we will do. We describe the model as G (t) = F 1 (G, H) + J(t) H (t) = F 2 (G, H) where the function J is the external rate at which blood glucose concentration is being increased.
28 The Background for the Diabetes Model The idea of our model for diagnosing diabetes from the GTT is to find a simple dynamical model of this complicated blood glucose regulatory system in which the values of two parameters would give a nice criterion for distinguishing normal individuals from those with mild diabetes or those who are pre diabetic. Here is what we will do. We describe the model as G (t) = F 1 (G, H) + J(t) H (t) = F 2 (G, H) where the function J is the external rate at which blood glucose concentration is being increased. There are two nonlinear interaction functions F 1 and F 2 because we know G and H have complicated interactions.
29 The Background for the Diabetes Model Let s assume G and H have achieved optimal values G 0 and H 0 by the time the fasting patient has arrived at the hospital.
30 The Background for the Diabetes Model Let s assume G and H have achieved optimal values G 0 and H 0 by the time the fasting patient has arrived at the hospital. Hence, we don t expect to have any contribution to G (0) and H (0); i.e. F 1 (G 0, H 0 ) = 0 and F 2 (G 0, H 0 ) = 0.
31 The Background for the Diabetes Model Let s assume G and H have achieved optimal values G 0 and H 0 by the time the fasting patient has arrived at the hospital. Hence, we don t expect to have any contribution to G (0) and H (0); i.e. F 1 (G 0, H 0 ) = 0 and F 2 (G 0, H 0 ) = 0. We are interested in the deviation of G and H from their optimal values G 0 and H 0, so let g = G G 0 and h = H H 0. We can then write G = G 0 + g and H = H 0 + h.
32 The Background for the Diabetes Model The model can then be rewritten as (G 0 + g) (t) = F 1 (G 0 + g, H 0 + h) + J(t) (H 0 + h) (t) = F 2 (G 0 + g, H 0 + h)
33 The Background for the Diabetes Model The model can then be rewritten as (G 0 + g) (t) = F 1 (G 0 + g, H 0 + h) + J(t) (H 0 + h) (t) = F 2 (G 0 + g, H 0 + h) or g (t) = F 1 (G 0 + g, H 0 + h) + J(t) h (t) = F 2 (G 0 + g, H 0 + h)
34 The Background for the Diabetes Model We can expand our functions using tangent plane approximations: F 1 (G 0 + g, H 0 + h) = F 1 (G 0, H 0 ) + F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h + E F1 F 2 (G 0 + g, H 0 + h) = F 2 (G 0, H 0 ) + F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h + E F2
35 The Background for the Diabetes Model We can expand our functions using tangent plane approximations: F 1 (G 0 + g, H 0 + h) = F 1 (G 0, H 0 ) + F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h + E F1 F 2 (G 0 + g, H 0 + h) = F 2 (G 0, H 0 ) + F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h + E F2 but the terms F 1 (G 0, H 0 ) = 0 and F 1 (G 0, H 0 ) = 0, so we can simplify to F 1 (G 0 + g, H 0 + h) = F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h + E F1 F 2 (G 0 + g, H 0 + h) = F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h + E F2
36 The Background for the Diabetes Model It seems reasonable to assume that since we are so close to ordinary operating conditions, the errors E F1 and E F2 will be negligible. Thus our model approximation is g (t) = F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h + J(t) h (t) = F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h
37 The Background for the Diabetes Model It seems reasonable to assume that since we are so close to ordinary operating conditions, the errors E F1 and E F2 will be negligible. Thus our model approximation is g (t) = F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h + J(t) h (t) = F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h Next, we can reason out the algebraic signs of the four partial derivatives.
38 The Diabetes Detection Model Our model approximation is g (t) = F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h h (t) = F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h
39 The Diabetes Detection Model Our model approximation is g (t) = F 1 g (G 0, H 0 ) g + F 1 h (G 0, H 0 ) h h (t) = F 2 g (G 0, H 0 ) g + F 2 h (G 0, H 0 ) h Next, we can reason out the algebraic signs are F 1 g (G 0, H 0 ) = F 1 h (G 0, H 0 ) = F 2 g (G 0, H 0 ) = + F 2 h (G 0, H 0 ) = The arguments for these algebraic signs come from our understanding of the physiological processes that are going on here.
40 The Diabetes Detection Model Let s look at a small positive deviation g from the optimal value G 0 while letting the net hormone concentration be fixed at H 0. Then our model approximation is g (t) = F 1 g (G 0, H 0 ) g since we do not have an external input here. At a state where we have an increase in blood sugar levels over optimal, i.e. g > 0, the other hormones such as cortisol and glucagon will try to regulate the blood sugar level down by increasing their concentrations and for example storing more sugar into glycogen.
41 The Diabetes Detection Model Let s look at a small positive deviation g from the optimal value G 0 while letting the net hormone concentration be fixed at H 0. Then our model approximation is g (t) = F 1 g (G 0, H 0 ) g since we do not have an external input here. At a state where we have an increase in blood sugar levels over optimal, i.e. g > 0, the other hormones such as cortisol and glucagon will try to regulate the blood sugar level down by increasing their concentrations and for example storing more sugar into glycogen. Hence, the term F1 g (G 0, H 0 ) should be negative as here g is negative as g should be decreasing.
42 The Diabetes Detection Model Let s look at a small positive deviation g from the optimal value G 0 while letting the net hormone concentration be fixed at H 0. Then our model approximation is g (t) = F 1 g (G 0, H 0 ) g since we do not have an external input here. At a state where we have an increase in blood sugar levels over optimal, i.e. g > 0, the other hormones such as cortisol and glucagon will try to regulate the blood sugar level down by increasing their concentrations and for example storing more sugar into glycogen. Hence, the term F1 g (G 0, H 0 ) should be negative as here g is negative as g should be decreasing. So we model this as F1 g (G 0, H 0 ) = m 1 for some positive number m 1.
43 The Diabetes Detection Model Now consider a positive change in h from the optimal level while keeping at the optimal level G 0. Then the model is g (t) = F 1 h (G 0, H 0 ) h
44 The Diabetes Detection Model Now consider a positive change in h from the optimal level while keeping at the optimal level G 0. Then the model is g (t) = F 1 h (G 0, H 0 ) h and since h > 0, this means the net hormone concentration is up which we interpret as insulin above normal. This means blood sugar levels go down which implies g is negative again.
45 The Diabetes Detection Model Now consider a positive change in h from the optimal level while keeping at the optimal level G 0. Then the model is g (t) = F 1 h (G 0, H 0 ) h and since h > 0, this means the net hormone concentration is up which we interpret as insulin above normal. This means blood sugar levels go down which implies g is negative again. Thus, F1 h (G 0, H 0 ) must be negative which means we model it as F 1 h (G 0, H 0 ) = m 2 for some positive m 2.
46 The Diabetes Detection Model If we have a small positive deviation g from the optimal value G 0 while letting the net hormone concentration be fixed at H 0, we have h (t) = F 2 g (G 0, H 0 ) g.
47 The Diabetes Detection Model If we have a small positive deviation g from the optimal value G 0 while letting the net hormone concentration be fixed at H 0, we have h (t) = F 2 g (G 0, H 0 ) g. Again, since g is positive, this means we are above normal blood sugar levels which implies mechanisms are activated to bring the level down. Hence h > 0 as we have increasing net hormone levels.
48 The Diabetes Detection Model If we have a small positive deviation g from the optimal value G 0 while letting the net hormone concentration be fixed at H 0, we have h (t) = F 2 g (G 0, H 0 ) g. Again, since g is positive, this means we are above normal blood sugar levels which implies mechanisms are activated to bring the level down. Hence h > 0 as we have increasing net hormone levels. Thus, we must have F2 g (G 0, H 0 ) = m 3 for some positive m 3.
49 The Diabetes Detection Model Finally, if we have a positive deviation h from optimal while blood sugar levels are optimal, the model is h (t) = F 2 h (G 0, H 0 ) h.
50 The Diabetes Detection Model Finally, if we have a positive deviation h from optimal while blood sugar levels are optimal, the model is h (t) = F 2 h (G 0, H 0 ) h. Since h is positive, we have the concentrations of the hormones that pull glucose out of the blood stream are above optimal. This means that too much sugar is being removed as so the regulatory mechanisms will act to stop this action implying h < 0.
51 The Diabetes Detection Model Finally, if we have a positive deviation h from optimal while blood sugar levels are optimal, the model is h (t) = F 2 h (G 0, H 0 ) h. Since h is positive, we have the concentrations of the hormones that pull glucose out of the blood stream are above optimal. This means that too much sugar is being removed as so the regulatory mechanisms will act to stop this action implying h < 0. This tells us F2 g (G 0, H 0 ) = m 4 for some positive constant m 4.
52 The Diabetes Detection Model Hence, the four partial derivatives at the optimal points can be defined by four positive numbers m 1, m 2, m 3 and m 4 as follows: F 1 g (G 0, H 0 ) = m 1 F 1 h (G 0, H 0 ) = m 2 F 2 g (G 0, H 0 ) = +m 3 F 2 h (G 0, H 0 ) = m 4
53 The Diabetes Detection Model Hence, the four partial derivatives at the optimal points can be defined by four positive numbers m 1, m 2, m 3 and m 4 as follows: F 1 g (G 0, H 0 ) = m 1 F 1 h (G 0, H 0 ) = m 2 F 2 g (G 0, H 0 ) = +m 3 F 2 h (G 0, H 0 ) = m 4 Our model dynamics are thus approximated by g (t) = m 1 g m 2 h + J(t) h (t) = m 3 g m 4 h
54 The Diabetes Detection Model Hence, the four partial derivatives at the optimal points can be defined by four positive numbers m 1, m 2, m 3 and m 4 as follows: F 1 g (G 0, H 0 ) = m 1 F 1 h (G 0, H 0 ) = m 2 F 2 g (G 0, H 0 ) = +m 3 F 2 h (G 0, H 0 ) = m 4 Our model dynamics are thus approximated by This implies g (t) = m 1 g m 2 h + J(t) h (t) = m 3 g m 4 h g (t) = m 1 g m 2 h + J (t)
55 The Diabetes Detection Model Now plug in the formula for h to get g (t) = m 1 g m 2 (m 3 g m 4 h) + J (t) = m 1 g m 2 m 3 g + m 2 m 4 h + J (t).
56 The Diabetes Detection Model Now plug in the formula for h to get g (t) = m 1 g m 2 (m 3 g m 4 h) + J (t) = m 1 g m 2 m 3 g + m 2 m 4 h + J (t). But we can use the g equation to solve for h. This gives m 2 h = g (t) m 1 g + J(t)
57 The Diabetes Detection Model Now plug in the formula for h to get g (t) = m 1 g m 2 (m 3 g m 4 h) + J (t) = m 1 g m 2 m 3 g + m 2 m 4 h + J (t). But we can use the g equation to solve for h. This gives m 2 h = g (t) m 1 g + J(t) which leads to g (t) = m 1 g m 2 m 3 g + m 4 ( g (t) m 1 g + J(t)) + J (t) = (m 1 + m 4 ) g (m 1 m 4 + m 2 m 3 )g + m 4 J(t) + J (t).
58 The Diabetes Detection Model Now plug in the formula for h to get g (t) = m 1 g m 2 (m 3 g m 4 h) + J (t) = m 1 g m 2 m 3 g + m 2 m 4 h + J (t). But we can use the g equation to solve for h. This gives m 2 h = g (t) m 1 g + J(t) which leads to g (t) = m 1 g m 2 m 3 g + m 4 ( g (t) m 1 g + J(t)) + J (t) = (m 1 + m 4 ) g (m 1 m 4 + m 2 m 3 )g + m 4 J(t) + J (t). So our final model is g (t) + (m 1 + m 4 ) g + (m 1 m 4 + m 2 m 3 )g = m 4 J(t) + J (t).
59 The Diabetes Detection Model Let α = (m 1 + m 4 )/2 and ω 2 = m 1 m 4 + m 2 m 3 and we can rewrite as where S(t) = m 4 J(t) + J (t). g (t) + 2α g + ω 2 g = S(t).
60 The Diabetes Detection Model Let α = (m 1 + m 4 )/2 and ω 2 = m 1 m 4 + m 2 m 3 and we can rewrite as where S(t) = m 4 J(t) + J (t). g (t) + 2α g + ω 2 g = S(t). Now the right hand side here is zero except for the very short time interval when the glucose load is being ingested. Hence, we can simply search for the solution to the homogeneous model g (t) + 2α g + ω 2 g = 0.
61 The Diabetes Detection Model Let α = (m 1 + m 4 )/2 and ω 2 = m 1 m 4 + m 2 m 3 and we can rewrite as where S(t) = m 4 J(t) + J (t). g (t) + 2α g + ω 2 g = S(t). Now the right hand side here is zero except for the very short time interval when the glucose load is being ingested. Hence, we can simply search for the solution to the homogeneous model g (t) + 2α g + ω 2 g = 0. The roots of the characteristic equation here are r = 2α ± 4α 2 4ω 2 2 = α ± α 2 ω 2.
62 The Diabetes Detection Model The most interesting case is if we have complex roots. In that case, α 2 ω 2 < 0. Let Ω 2 = α 2 ω 2. Then, the general phase shifted solution has the form g = Re αt cos(ωt δ) which implies G = G 0 + Re αt cos(ωt δ).
63 The Diabetes Detection Model The most interesting case is if we have complex roots. In that case, α 2 ω 2 < 0. Let Ω 2 = α 2 ω 2. Then, the general phase shifted solution has the form g = Re αt cos(ωt δ) which implies G = G 0 + Re αt cos(ωt δ). Hence, our model has five unknowns to find: G 0, R, α, Ω and δ, The easiest way to do this is to measure G 0, the patient s initial blood glucose concentration, when the patient arrives. Then measure the blood glucose concentration N more times giving the data pairs (t 1, G 1 ), (t 2, G 2 ) and so on out to (t N, G N ).
64 The Diabetes Detection Model The most interesting case is if we have complex roots. In that case, α 2 ω 2 < 0. Let Ω 2 = α 2 ω 2. Then, the general phase shifted solution has the form g = Re αt cos(ωt δ) which implies G = G 0 + Re αt cos(ωt δ). Hence, our model has five unknowns to find: G 0, R, α, Ω and δ, The easiest way to do this is to measure G 0, the patient s initial blood glucose concentration, when the patient arrives. Then measure the blood glucose concentration N more times giving the data pairs (t 1, G 1 ), (t 2, G 2 ) and so on out to (t N, G N ). Then form the least squares error function E = N ( Gi G 0 Re αt i cos(ωt i δ) ) 2 i=1
65 Fitting the Data Here is some typical glucose versus time data. Time Glucose Level We will now try to find the parameter values which minimize the nonlinear least squares problem we have here. This appears to be a simple problem, but you will see all numerical optimization problems are actually fairly difficult. Our problem is to find the free parameters G o, R, α, Ω and δ which minimize E(G 0, R, α, Ω, δ) = N ( Gi G 0 Re αt i cos(ωt i δ) ) 2 i=1 For convenience, let X = [G 0, R, α, Ω, δ] and f i (X ) = G i G 0 Re αt i cos(ωt i δ); then we can rewrite the error function as E(X ) = N i=1 f i 2 (X ). Gradient descent requires the gradient of this error function. This is just a messy calculation;
66 Fitting the Data E G o = 2 N i=1 E N R = 2 i=1 N E α = 2 i=1 E N Ω = 2 E δ = 2 i=1 N i=1 f i (X ) f i G o f i (X ) f i R f i (X ) f i α f i (X ) f i Ω f i (X ) f i δ
67 Fitting the Data where the f i partials are given by f i G o = 1 f i R = e αt i cos(ωt i δ) f i α = t i Re αt i cos(ωt i δ), f i Ω = t i Re αt i sin(ωt i δ) f i δ = Re αt i sin(ωt i δ) Now gradient descent works like this: from a given point X 0, we find the gradient of E, E(X 0 ). We know that E decreases in the direction of the negative gradient so caculating E(X 0 λ E(X 0 )) for some positive number λ has the potential to decrease E. But how big or small should λ be? Choosing λ is a difficult problem in general.
68 Fitting the Data Now any vector D for which the value of E goes down is called a descent vector. Now suppose we are at the point X 0 and we want to know how much of the descent vector D to use. Note, if we use the amount ξ of the descent vector at X 0, we compute the new error value E(X 0 ξd(x 0 )). Let g(ξ) = E(X 0 ξd(x 0 )). We see g(0) = E(X 0 ) and given a first choice of ξ = λ, we have g(λ) = E(X 0 λd(x 0 )). Next, let Y = X 0 ξ D(X 0 ). Then, using the chain rule, we can calculate the derivative of g at 0. First, we have g (ξ) = < (E), D(X 0 ) > and using the normalized gradient of E as the descent vector, we find g (0) = (E).
69 Fitting the Data Now let s approximate g using a quadratic model. Since we are trying for a minimum, in general we try to take a step in the direction of the negative gradient which makes the error function go down. Then, we have g(0) = E(X 0 ) is less than g(λ) = E(X 0 λd(x 0 )) and the directional derivative gives g (0) = (E) < 0. Hence, if we approximate g by a simple quadratic model, g(ξ) = A + Bξ + Cξ 2, this model will have a unique minimizer and we can use the value of ξ where the minimum occurs as our next choice of descent step. This technique is called a Line Search Method and it is quite useful. To summarize, we fit our g model and find g(0) = E(X 0 ) g (0) = (E) g(λ) = A + Bλ + Cλ 2 = C = E(X 0 λd(x 0 )) E(X 0 ) + (E) λ λ 2 The minimum of this quadratic occurs at λ = B 2C and this will give us our next descent direction X 0 λ D(X 0 ).
70 Gradient Descent Implementation Let s get started on how to find the optimal parameters numerically. Along the way, we will show you how hard this is. We start with a minimal implementation. We have already discussed some root finding codes in so we have seen code kind of similar. But this will be a little different and it is good to have you see a bit about it. What is complicated here is that we have lots of functions that depend on the data we are trying to fit. So the number of functions depends on the size of the data set which makes it harder to set up.
71 Gradient Descent Implementation 1 f u n c t i o n [ E r r o r, G0, R, alpha, Omega, d e l t a, normgrad, update ] = DiabetesGradOne ( I n i t i a l, lambda, m a x i t e r, data ) % I n i t i a l g u e s s f o r G0, R, alpha, Omega, d e l t a % data = c o l l e c t i o n o f ( time, g l u c o s e ) p a i r s % maxiter = max number of i t s ; lambda = how much of descent to use % s e t u p l e a s t s q u a r e s e r r o r f u n c t i o n 6 N = l e n g t h ( d a t a ) ; E = 0. 0 ; Time = d a t a ( :, 1 ) ; G = data ( :, 2 ) ; % I n i t i a l = [ i n i t i a l G, i n i t i a l R, i n i t i a l alpha, i n i t i a l Omega, i n i t i a l d e l t a ] g = I n i t i a l ( 1 ) ; r = I n i t i a l ( 2 ) ; a = I n i t i a l ( 3 ) ; o = I n i t i a l ( 4 ) ; d = I n i t i a l ( 5 ) ; 11 f t, e q u i l G, r, a, o, d ) e q u i l G + r exp( a ˆ2 t ). ( c o s ( o t d ) ) ; E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) f o r i = 1 : m a x i t e r E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) ; % c a l c u l a t e e r r o r i f ( i == 1) 16 E r r o r ( 1 ) = E ; end E r r o r = [ E r r o r ; E ] ; [ grade, normgrad ] = E r r o r G r a d i e n t ( Time, G, g, r, a, o, d ) ; % f i n d g r a d o f E i f ( normgrad >1)% f i n d d e s c e n t d i r e c t i o n D 21 Descent = grade / normgrad ; e l s e Descent = grade ; end Del = [ g ; r ; a ; o ; d] lambda Descent ; 26 g = Del ( 1 ) ; r = Del ( 2 ) ; a = Del ( 3 ) ; o = Del ( 4 ) ; d = Del ( 5 ) ; end G0 = g ; R = r ; a l p h a = a ; Omega = o ; d e l t a = d ; E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) update = [ G0 ; R ; a l p h a ; Omega ; d e l t a ] ; 31 end
72 Gradient Descent Implementation Note inside this function, we call another function to calculate the gradient of the norm. This is given below. The inputs are these: f u n c t i o n [ grade, normgrad ] = E r r o r G r a d i e n t ( Time, G, g, r, a, o, d ) % Time i s the time data ; % G i s the glucose data 4 % g = c u r r e n t e q u i l G ; % r = c u r r e n t R % a = c u r r e n t a l p h a ; % o = c u r r e n t Omega ; % d = c u r r e n t d e l t a
73 Gradient Descent Implementation f u n c t i o n [ grade, normgrad ] = E r r o r G r a d i e n t ( Time, G, g, r, a, o, d ) 2 f e r r o r t, G, g, r, a, o, d ) G g r exp( a ˆ2 t ). ( c o s ( o t d ) ) ; g e r r o r t, r, a, o, d ) r exp( a ˆ2 t ). ( c o s ( o t d ) ) ; h e r r o r t, r, a, o, d ) r exp( a ˆ2 t ). ( s i n ( o t d ) ) ; N = l e n g t h ( Time ) ; sum = 0 ; 7 f o r k =1:N sum = sum + f e r r o r ( Time ( k ),G( k ), g, r, a, o, d ) ; end pepequilg = 2 sum ; sum = 0 ; 12 f o r k =1:N sum = sum + f e r r o r ( Time ( k ),G( k ), g, r, a, o, d ) g e r r o r ( Time ( k ), r, a, o, d ) / r ; end pepr = 2 sum ; sum = 0 ; 17 f o r k =1:N sum = sum + f e r r o r ( Time ( k ),G( k ), g, r, a, o, d ) Time ( k ) 2 a Time ( k ) g e r r o r ( Time ( k ), r, a, o, d ) ; end pepa = 2 sum ; sum = 0 ; 22 f o r k =1:N sum = sum + f e r r o r ( Time ( k ),G( k ), g, r, a, o, d ) Time ( k ) h e r r o r ( Time ( k ), r, a, o, d ) ; end pepo = 2 sum ; sum = 0 ; 27 f o r k =1:N sum = sum + f e r r o r ( Time ( k ),G( k ), g, r, a, o, d ) h e r r o r ( Time ( k ), r, a, o, d ) ; end pepd = 2 sum ; grade = [ pepequilg ; pepr ; pepa ; pepo ; pepd ] ; normgrad = norm ( grade ) ; 32 end
74 Gradient Descent Implementation We also need code for the error calculations which is given here. f u n c t i o n E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) % 3 % T = Time % G = G l o c u o s e v a l u e s % f = n o n l i n e a r i n s u l i n model % g, a, r, o, d = p a r a m e t e r s i n d i a b e t e s n o n l i n e a r model % N = s i z e o f d ata 8 N = l e n g t h (G) ; E = 0. 0 ; % c a l c u l a t e e r r o r f u n c t i o n f o r i = 1 :N E = E +( G( i ) f ( Time ( i ), g, r, a, o, d ) ) ˆ 2 ; 13 end
75 Gradient Descent Implementation After 100, 000 iterations we still do not have a good fit. Note we start with a small constant λ = 5.0e 4 here. Try it yourself. If you let this value be larger, the optimization spins out of control. Also, we have not said how we chose our initial values. Let s look at some run time results using this code. Data = [ 0, 9 5 ; 1, ; 2, ; 3, ] ; Time = Data ( :, 1 ) ; G = Data ( :, 2 ) ; f t, e q u i l G, r, a, o, d ) e q u i l G + r exp( a ˆ2 t ). ( c o s ( o t d ) ) ; time = l i n s p a c e ( 0, 3, 4 1 ) ; R I n i t i a l = ; G O I n i t i a l = 95 + R I n i t i a l ; A I n i t i a l = s q r t ( l o g ( 1 7 / 5 ) ) ; O I n i t i a l = p i ; d I n i t i a l = p i ; I n i t i a l = [ G O I n i t i a l ; R I n i t i a l ; A I n i t i a l ; O I n i t i a l ; d I n i t i a l ] ; [ E r r o r, G0, R, a lpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, 5. 0 e 4,20000, Data, 0 ) ; I n i t i a l E = E = o c t a v e :16> [ E r r o r, G0, R, alpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, 5. 0 e 4,40000, Data, 0 ) ; I n i t i a l E = E = o c t a v e :145> [ E r r o r, G0, R, alpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, 5. 0 e 4,100000, Data, 0 ) ; I n i t i a l E = E =
76 Gradient Descent Implementation We actually looked at the data on a sheet of paper and did some rough calculations to try for some decent values. We will leave that to you to figure out. If the initial values are poorly chosen, gradient descent optimization is a great way to generate really bad values! So be warned. You will have to exercise due diligence to find a sweet starting spot. We can see how we did by looking the resulting curve fit
77 Adding Line Search Now let s add line search and see if it gets better. We will also try scaling the data so all the variables in question are roughly the same size. For us, a good choice is to scale the G 0 and the R value by 50, although we could try other choices. We have already discussed line search for our problem, but here it is again in a quick nutshell. If we are minimizing a function of M variables, say f (X ), then if we are at the point X 0, we can look at the slice of this function if we move out from the base point X 0 is the direction of the negative gradient, (f (X 0 ) = (f 0 ). Define a function of the single variable ξ as g(ξ) = f (X 0 + ξ (f 0 )). Then, we can try to approximate g as a quadratic, g(ξ) = A + Bξ + Cξ 2. Of course, the actual function might not be approximated nicely by such a quadratic, but it is worth a shot! Once we fit the parameters A, B and C, we see this quadratic model is minimized at λ = B 2C. The code now adds the line search code which is contained in the block
78 Adding Line Search f u n c t i o n [ l a m b d a s t a r, D eloptimal ] = L i n e S e a r c h ( Del, EStart, E F u l l S t e p, Base, Descent, normgrad, lambda ) % % Del i s the update due to the step lambda % E s t a r t i s E a t t h e base v a l u e 5 % EFullStep i s E at the given lambda value % Base i s t h e c u r r e n t p a r a m e t e r v e c t o r % Descent i s t h e c u r r e n t d e s c e n t v e c t o r % normgrad o f g r a d o f E a t t h e base v a l u e % lambda i s t h e c u r r e n t s t e p 10 % % we have enough i n f o r m a t i o n to do a l i n e s e a r c h h e r e A = EStart ; BCheck = normgrad ; C = ( EFullStep A BCheck lambda ) /( lambda ˆ2) ; 15 lambdastar = BCheck /(2 C) ; i f (C<0 lambdastar < 0) % we a r e g o i n g to a maximum on t h e l i n e s e a r c h ; r e j e c t DelOptimal = Del ; e l s e 20 % we have a minimum on the l i n e search DelOptimal= Base lambdastar Descent ; end end
79 Adding Line Search f u n c t i o n [ E r r o r, G0, R, alpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, lambda, m a x i t e r, data, d o l i n e s e a r c h ) 2 % % I n i t i a l g u e s s f o r G0, R, alpha, Omega, d e l t a % data = c o l l e c t i o n o f ( time, g l u c o s e ) p a i r s % t o l = s t o p t o l e r a n c e % m a x i t e r = maximum number o f i t e r a t i o n s to use 7 % lambda = how much o f t h e d e s c e n t v e c t o r to use % % s e t u p l e a s t s q u a r e s e r r o r f u n c t i o n N = l e n g t h ( d a t a ) ; E = 0. 0 ; 12 Time = data ( :, 1 ) ; G = data ( :, 2 ) ; % I n i t i a l = [ i n i t i a l G, i n i t i a l R, i n i t i a l alpha, i n i t i a l Omega, i n i t i a l d e l t a ] g = I n i t i a l ( 1 ) ; r = I n i t i a l ( 2 ) ; 17 a = I n i t i a l ( 3 ) ; o = I n i t i a l ( 4 ) ; d = I n i t i a l ( 5 ) ; f t, e q u i l G, r, a, o, d ) e q u i l G + r exp( a ˆ2 t ). ( c o s ( o t d ) ) ; I n i t i a l E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) 22 % c a l c u l a t e e r r o r The rest of the code is on the next slide.
80 Adding Line Search f o r i = 1 : m a x i t e r E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) ; 3 i f ( i == 1) E r r o r ( 1 ) = E ; end E r r o r = [ E r r o r ; E ] ; [ grade, normgrad ] = E r r o r G r a d i e n t ( Time, G, g, r, a, o, d ) ; % C a l c u l a t e E r r o r g r a d i e n t 8 % f i n d d e s c e n t d i r e c t i o n i f ( normgrad >1) Descent = grade / normgrad ; e l s e Descent = grade ; 13 end Base = [ g ; r ; a ; o ; d ] ; Del = Base lambda Descent ; newg = Del ( 1 ) ; newr = Del ( 2 ) ; newa = Del ( 3 ) ; newo = Del ( 4 ) ; newd = Del ( 5 ) ; E S t a r t = E ; 18 E F u l l S t e p = D i a b e t e s E r r o r ( f, newg, newr, newa, newo, newd, Time, G) ; i f ( d o l i n e s e a r c h ==1) [ l a m b d a s t a r, D e l O p t i m a l ] = L i n e S e a r c h ( Del, EStart, E F u l l S t e p, Base, Descent, normgrad, lambda ) ; g = DelOptimal ( 1 ) ; r = DelOptimal ( 2 ) ; a = DelOptimal ( 3 ) ; o = DelOptimal ( 4 ) ; d = DelOptimal ( 5 ) ; EOptimalStep = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) ; 23 e l s e g = Del ( 1 ) ; r = Del ( 2 ) ; a = Del ( 3 ) ; o = Del ( 4 ) ; d = Del ( 5 ) ; end end G0 = g ; R = r ; a l p h a = a ; Omega = o ; d e l t a = d ; 28 E = D i a b e t e s E r r o r ( f, g, r, a, o, d, Time, G) update = [ G0 ; R ; a l p h a ; Omega ; d e l t a ] ; end
81 Adding Line Search We can see how this is working by letting some of the temporary calculations print. Here are two iterations of line search printing out the A, B and C and the relevant energy values. Our initial values don t matter much here as we are just checking out the line search algorithm. Data = [ 0, 9 5 ; 1, ; 2, ; 3, ] ; Time = Data ( :, 1 ) ; G = Data ( :, 2 ) ; f t, e q u i l G, r, a, o, d ) e q u i l G + r exp( a ˆ2 t ). ( c o s ( o t d ) ) ; time = l i n s p a c e ( 0, 3, 4 1 ) ; R I n i t i a l = 85 17/21; G O I n i t i a l = 95 + R I n i t i a l ; A I n i t i a l = s q r t ( l o g ( 1 7 / 4 ) ) ; O I n i t i a l = p i ; d I n i t i a l = p i ; I n i t i a l = [ G O I n i t i a l ; R I n i t i a l ; A I n i t i a l ; O I n i t i a l ; d I n i t i a l ] ; [ E r r o r, G0, R, a lpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, 5. 0 e 4,2, Data, 1 ) ; I n i t i a l E = E F u l l S t e p = A = BCheck = C = e+05 lambdastar = e 04 E F u l l S t e p = A = BCheck = C = e+05 lambdastar = e 04 E =
82 Adding Line Search Now let s remove those prints and let it run for awhile. We are using the original data here to try to find the fit. o c t a v e :90> [ E r r o r, G0, R, alpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, 5. 0 e 4,20000, Data, 1 ) ; I n i t i a l E = E = o c t a v e :93> [ E r r o r, G0, R, alpha, Omega, d e l t a, normgrad, update ] = D i a b e t e s G r a d ( I n i t i a l, 5. 0 e 4,30000, Data, 1 ) ; I n i t i a l E = E = We have success! The line search got the job done in 30, 000 iterations while the attempt using just gradient descent without line search failed. but remember, we do additional processing at each step.
83 Adding Line Search We show the resulting curve fit here
84 Adding Line Search The five parameter values that make this error a minimum have been found using line search. Numerous experiments have shown if T 0 = 2π/Ω then whent 0 < 4 hours, the patient is normal and if T 0 is much larger than that, the patient has mild diabetes. Finally, the optimal values of the parameters are update = % G % R % a l p h a % Omega % d e l t a The critical value of 2π/Ω = here which is less than 4 so this patient is normal!
85 Adding Line Search The five parameter values that make this error a minimum have been found using line search. Numerous experiments have shown if T 0 = 2π/Ω then whent 0 < 4 hours, the patient is normal and if T 0 is much larger than that, the patient has mild diabetes. Finally, the optimal values of the parameters are update = % G % R % a l p h a % Omega % d e l t a The critical value of 2π/Ω = here which is less than 4 so this patient is normal!
86 Adding Line Search The five parameter values that make this error a minimum have been found using line search. Numerous experiments have shown if T 0 = 2π/Ω then whent 0 < 4 hours, the patient is normal and if T 0 is much larger than that, the patient has mild diabetes. Note how much of our material from this course was involved in building this model! Finally, the optimal values of the parameters are update = % G % R % a l p h a % Omega % d e l t a The critical value of 2π/Ω = here which is less than 4 so this patient is normal!
Mathematical Modeling of Diabetes
Mathematical Modeling of Diabetes Dr. Shyam Kamal Department of Systems Design and Informatics Kyushu Institute of Technology Japan kamal@ces.kyutech.ac.jp 20 May 2016 Overview Role of Glucose Homeostasis
More informationProject for Math. 224 DETECTION OF DIABETES
Project for Math. 224 DETECTION OF DIABETES Diabetes is a disease of metabolism which is characterized by too much sugar in the blood and urine. Because of the lack of insulin (a hormone), the patient
More informationChapter 3: Linear & Non-Linear Interaction Models
Chapter 3: 155/226 Chapter develops the models above to examine models which involve interacting species or quantities. Models lead to simultaneous differential equations for coupled quantites due to the
More informationEndocrine System. Regulating Blood Sugar. Thursday, December 14, 17
Endocrine System Regulating Blood Sugar Stress results in nervous and hormonal responses. The adrenal glands are located above each kidney. Involved in stress response. Stress Upsets Homeostasis Stress
More informationHuman Biochemistry. Hormones
Human Biochemistry Hormones THE ENDOCRINE SYSTEM THE ENDOCRINE SYSTEM THE ENDOCRINE SYSTEM The ENDOCRINE SYSTEM = the organ system that regulates internal environment conditions by secreting hormones into
More informationThe Chemostat: Stability at Steady States. Chapter 5: Linear & Non-Linear Interaction Models. So, in dimensional form, α 1 > 1 corresponds to
Introduction & Simple Models Logistic Growth Models The Chemostat: Stability at Steady States 1 So, in dimensional form, α 1 > 1 corresponds to K max < V F. As K max is max bacterial repro rate with unlimited
More informationAdditional Case Study: Glands and Hormones
Student Worksheet Additional Case Study: Glands and Hormones LSM 8.5-2 This activity can be done individually or in pairs. Prepare the pieces ahead of time. Materials For each student (or pair): one copy
More informationHormones and Homeostasis
Hormones and Homeostasis The endocrine system is a system of organs that releases chemical message molecules, called hormones, into the blood. Unlike the nervous system whose action helps the body react
More informationOverview. Physiology 1. The Gastrointestinal Tract. Guyton section XI
Overview Physiology 1 The Gastrointestinal Tract Guyton section XI Basic functions of the GI tract Digestion Secretion Absorption Motility Basic functions of the GI tract Digestion: : Dissolving and breaking
More informationChapter 20 Endocrine System
Chapter 20 Endocrine System The endocrine system consists of glands and tissues that secrete Hormones are chemicals that affect other glands or tissues, many times far away from the site of hormone production
More informationLESSON 2.4 WORKBOOK. Part two: Glucose homeostasis in the blood Un-Storing energy
DEFINITIONS OF TERMS Fasting A state of abstinence from all food or drinks that provide calories. For a complete list of defined terms, see the Glossary. LESSON 2.4 WORKBOOK Part two: Glucose homeostasis
More informationEB Education Revision Guide. How to work with Homeostasis: Part 2 Blood Glucose Regulation
EB Education Revision Guide How to work with Homeostasis: Part 2 Blood Glucose Regulation Blood Glucose Regulation a) Why your body regulates glucose levels What you need to know about Homeostasis: Part
More information9.3 Stress Response and Blood Sugar
9.3 Stress Response and Blood Sugar Regulate Stress Response Regulate Blood Sugar Stress Response Involves hormone pathways that regulate metabolism, heart, rate and breathing The Adrenal Glands a pair
More informationChapter 41. Lecture 14. Animal Hormones. Dr. Chris Faulkes
Chapter 41 Lecture 14 Animal Hormones Dr. Chris Faulkes Animal Hormones Aims: To appreciate the variety and roles of hormones in the body To understand the basic types of hormones To understand how hormones
More informationModule 2 Endocrine System
Module 2 Endocrine System Student Name: 1 Lesson 1 Lesson 2 Lesson 3 Lesson 4 Lesson 5 Total Marks Total Possible Marks 10 8 14 21 16 69 Your Mark Teacher Comments: 2 (10 marks) Lesson 1: Structure and
More informationBiology 30. Morinville Community High School. Unit 2: Endocrine System. Name:
Biology 30 Morinville Community High School Unit 2: Endocrine System Name: 2 Endocrine System Unit Outline Chapter 13 text p. 434-471 Key Concept A: The endocrine system and nervous system both mediate
More informationCell Processes: Osmosis
www.mathbench.umd.edu./modules-au Osmosis May 2015 page 1 Cell Processes: Osmosis URL: http://mathbench.org.au/cellular-processes/osmosis/ Learning Outcomes After completing this module you should be able
More informationwww.lessonplansinc.com Topic: Feedback Mechanism Lab Summary: Students will model how the endocrine system uses feedback mechanisms to maintain homeostasis. Goals & Objectives: Students will be able to
More informationHomeostatic Regulation
Homeostatic Regulation A hormone is :a Water-soluble hormones: Composed of amino acids and bind a receptor protein on the of the target cell. This starts a signal cascade inside the cell and the signal
More informationSession 3 - The Type 2 Diabetes Cycle Transcript
Session 3 - The Type 2 Diabetes Cycle Transcript Essentially type 2 diabetes needs to treated with medications due to our modern lifestyle. Unfortunately it has been shown that only a very small minority
More informationRespiration Cellular Respiration Understand the relationship between glucose breakdown and ATP when you burn glucose with the help of oxygen, it
Respiration Cellular Respiration Understand the relationship between glucose breakdown and ATP when you burn glucose with the help of oxygen, it traps chemical energy into ATP Energy found in glucose stores
More informationUnit 4 Homeostasis. The term homeostasis refers to the body s attempt. Your body systems must to maintain a stable internal environment -
Unit 4 Homeostasis The term homeostasis refers to the body s attempt Your body systems must to maintain a stable internal environment - The body is trying to maintain, through a series of monitored adjustments.
More information5.0 HORMONAL CONTROL OF CARBOHYDRATE METABOLISM
5.0 HORMONAL CONTROL OF CARBOHYDRATE METABOLISM Introduction: Variety of hormones and other molecules regulate the carbohydrates metabolism. Some of these have already been cited in previous sections.
More informationEndocrine System. Chemical Control
Endocrine System Chemical Control Endocrine System - the system that secretes hormones in the body - hormones can last for minutes or for hours - a major gland, once called the master gland, is the pituitary
More informationChapter 4 Reading Guide
Chapter 4 Reading Guide 1. How many covalent bonds does carbon (C) form? 2. What is the chemical formula for glucose? 3. List the major dietary monosaccharides and disaccharides. What are the components
More informationChapter 26. Hormones and the Endocrine System. Lecture by Edward J. Zalisko
Chapter 26 Hormones and the Endocrine System PowerPoint Lectures for Biology: Concepts & Connections, Sixth Edition Campbell, Reece, Taylor, Simon, and Dickey Copyright 2009 Pearson Education, Inc. Lecture
More informationEndocrine System Notes
Endocrine System Notes is the tendency to maintain a stable internal environment. - parts of the body that secrete hormones directly into the body. - parts of the body that make secretions which travel
More informationChemical Regulation. Chapter 26. Testosterone and Male Aggression: Is There a Link? THE NATURE OF CHEMICAL REGULATION
Chapter 6 Chemical Regulation PowerPoint Lectures for Biology: Concepts and Connections, Fifth Edition Campbell, Reece, Taylor, and Simon Testosterone and Male Aggression: Is There a Link? Among male animals,
More informationMoh Tarek. Razi Kittaneh. Jaqen H ghar
14 Moh Tarek Razi Kittaneh Jaqen H ghar Naif Karadsheh Gluconeogenesis is making glucose from non-carbohydrates precursors. Although Gluconeogenesis looks like Glycolysis in many steps, it is not the simple
More informationCh. 39 Endocrine System
Ch. 39 Endocrine System The Endocrine System The endocrine system is made up of glands that release their products into the bloodstream These products deliver messages throughout the body Hormones Hormones
More informationWhat systems are involved in homeostatic regulation (give an example)?
1 UNIVERSITY OF PNG SCHOOL OF MEDICINE AND HEALTH SCIENCES DIVISION OF BASIC MEDICAL SCIENCES DISCIPLINE OF BIOCHEMISTRY AND MOLECULAR BIOLOGY GLUCOSE HOMEOSTASIS (Diabetes Mellitus Part 1): An Overview
More informationConcentration of drug [A]
Pharmacology Semester 1 page 1 of 5 PHARMACODYNAMICS 1 Receptor occupancy - mass action The interaction of a drug with a receptor is reversible due to interactions via weak bonds (not covalent). [A] +
More informationEndocrine System. Chapter 20. Endocrine Glands and Hormones. The Endocrine System. Endocrine glands
Chapter 20 Endocrine System Endocrine Glands and Hormones The endocrine system consists of glands and tissues that secrete hormones Hormones are chemicals that affect other glands or tissues, many times
More informationGlucoregulation 1 of 27 Boardworks Ltd 2012
Glucoregulation 1 of 27 Boardworks Ltd 2012 2 of 27 Boardworks Ltd 2012 Glucose 3 of 27 Boardworks Ltd 2012 Glucose is a type of sugar. It is the basic fuel for aerobic respiration, which provides the
More informationHypoglycemia. When recognized early, hypoglycemia can be treated successfully.
Hypoglycemia Introduction Hypoglycemia is a condition that causes blood sugar level to drop dangerously low. It mostly shows up in diabetic patients who take insulin. When recognized early, hypoglycemia
More informationControl of Glucose Metabolism
Glucose Metabolism Control of Glucose Metabolism The pancreas is both an exocrine and endocrine gland. It secretes digestive enzymes into the duodenum (exocrine) and 3 specific hormones into the bloodstream
More informationTestosterone and other male hormones seem to be related to aggressive behavior in some species
Testosterone and Male Aggression Testosterone and other male hormones seem to be related to aggressive behavior in some species In the fish species Oreochromis mossambicus, elevated levels have been found
More informationHormonal regulation of. Physiology Department Medical School, University of Sumatera Utara
Hormonal regulation of nutrient metabolism Physiology Department Medical School, University of Sumatera Utara Homeostasis & Controls Successful compensation Homeostasis reestablished Failure to compensate
More informationEndocrine System. Endocrine vs. Exocrine. Bio 250 Human Anatomy & Physiology
Endocrine System Bio 250 Human Anatomy & Physiology Endocrine vs. Exocrine Endocrine glands secrete their products called hormones into body fluids (the internal environment) Exocrine glands secrete their
More informationnumber Done by Corrected by Doctor
number 20 Done by Corrected by Rana Ghassan Doctor Only 4 questions in the mid-term exam are based on the 4 lectures to be given by Dr Faisal. Dr Faisal will give us 10 lectures, the first 4 are included
More informationStanford Youth Diabetes Coaches Program Instructor Guide Class #1: What is Diabetes? What is a Diabetes Coach? Sample
Note to Instructors: YOU SHOULD HAVE ENOUGH COPIES OF THE QUIZ AND THE HOMEWORK TO PASS OUT TO EACH STUDENT. Be sure to use the NOTES view in Powerpoint for what to cover during class. It is important
More informationEnergy stores in different organs for a 155 lb male, in Calories
Energy stores in different organs for a 155 lb male, in Calories Organ Glucose/ Glycogen Triacyl Glycerols* Liver 400 450 400 Brain 8 0 0 Mobile Proteins Muscle 1,200 450 24,000 Adipose Tissue 80 135,000
More informationChapter 20. Endocrine System Chemical signals coordinate body functions Chemical signals coordinate body functions. !
26.1 Chemical signals coordinate body functions Chapter 20 Endocrine System! Hormones Chemical signals Secreted by endocrine glands Usually carried in the blood Cause specific changes in target cells Secretory
More informationENERGY FROM INGESTED NUTREINTS MAY BE USED IMMEDIATELY OR STORED
QUIZ/TEST REVIEW NOTES SECTION 1 SHORT TERM METABOLISM [METABOLISM] Learning Objectives: Identify primary energy stores of the body Differentiate the metabolic processes of the fed and fasted states Explain
More informationI How the human body is put together and organized
LECTURE OUTLINE, Chapter 3: The Remarkable Body I How the human body is put together and organized The human body is composed of billions of cells that need energy, water, and. Cells are organized into
More informationEndocrine System Hormones (Ch. 45)
Endocrine System Hormones (Ch. 45) Regulation Why are hormones needed? chemical messages from one body part to another communication needed to coordinate whole body daily homeostasis & regulation of large
More informationSlim & Save Re-feed Programme
Slim & Save Re-feed Programme Before you start your re-feed plan its important to be aware of the way a ketogenic diet works in respect of weight loss and weight gain and the part glycogen plays in this
More informationEndocrine System Hormones. AP Biology
Endocrine System Hormones 2007-2008 Regulation Why are hormones needed? u chemical messages from one body part to another u communication needed to coordinate whole body u daily homeostasis & regulation
More informationIntroduction Causes Of Out Of Control Hunger...3. Dehydration...3. Hyperthyroidism...4. Your Hormones Are off...4. You Skip Meals...
Table of Contents Introduction...2 10 Causes Of Out Of Control Hunger...3 Dehydration...3 Hyperthyroidism...4 Your Hormones Are off...4 You Skip Meals...4 You re Not Sleeping...5 You Rush Through Meals...5
More informationUNIVERSITY OF PNG SCHOOL OF MEDICINE AND HEALTH SCIENCES DIVISION OF BASIC MEDICAL SCIENCES DISCIPLINE OF BIOCHEMISTRY AND MOLECULAR BIOLOGY
1 UNIVERSITY OF PNG SCHOOL OF MEDICINE AND HEALTH SCIENCES DIVISION OF BASIC MEDICAL SCIENCES DISCIPLINE OF BIOCHEMISTRY AND MOLECULAR BIOLOGY GLUCOSE HOMEOSTASIS An Overview WHAT IS HOMEOSTASIS? Homeostasis
More information6. The diagram below represents an interaction between parts of an organism.
Endocrine Review 1. Base your answer to the following question on the diagram below and on your knowledge of biology. Each arrow in the diagram represents a different hormone released by the pituitary
More informationCh45: Endocrine System
Ch45: Endocrine System Endocrine System Homeostasis is the tendency to maintain a stable internal environment. Function = coordinate and control the body with hormones to maintain homeostasis Works with
More informationCampbell's Biology: Concepts and Connections, 7e (Reece et al.) Chapter 26 Hormones and the Endocrine System Multiple-Choice Questions
Campbell's Biology: Concepts and Connections, 7e (Reece et al.) Chapter 26 Hormones and the Endocrine System 26.1 Multiple-Choice Questions 1) Hormones are chemicals produced by the endocrine system that
More informationCarbohydrates. Three Types of Carbs Fiber. Does not impact the blood sugar. Cleaners Controllers May be eaten on a Ketogenic diet
The Keto Challenge Carbohydrates Carbohydrates Three Types of Carbs STARCHES (complex carbs) Take a little longer to digest b Slower rise in blood glucose levels Trigger the Fat Truck (Slower, but still
More information3 The Endocrine System
CHAPTER 25 3 The Endocrine System SECTION Communication and Control BEFORE YOU READ After you read this section, you should be able to answer these questions: Why is the endocrine system important? How
More informationExercise Physiology: Theory and Application to Fitness and Performance By Scott Powers & Edward Howley
Exercise Physiology: Theory and Application to Fitness and Performance By Scott Powers & Edward Howley Ch 5 Cell Signaling and the Hormonal Responses to Exercise Summary Created by Dan Hechler Class Lecture
More informationWhat is the endocrine system?
Name: Class: Date: What is the endocrine system? Biology Gr11A Physiological regulation in mammals Your body depends on the endocrine system to stay healthy. This system has special organs called glands
More informationComputer Science 101 Project 2: Predator Prey Model
Computer Science 101 Project 2: Predator Prey Model Real-life situations usually are complicated and difficult to model exactly because of the large number of variables present in real systems. Computer
More informationCarbohydrate Metabolism
Chapter 34 Carbohydrate Metabolism Carbohydrate metabolism is important for both plants and animals. Introduction to General, Organic, and Biochemistry, 10e John Wiley & Sons, Inc Morris Hein, Scott Pattison,
More informationPancreas. Endocrine pancreas - Islets of Langerhans A or alpha cells glucagon B or beta cells insulin Delta cells somatostatin
Endocrine System Pancreas Endocrine pancreas - Islets of Langerhans A or alpha cells glucagon B or beta cells insulin Delta cells somatostatin Glucagon & Metabolism Produced by beta cells of Islets Primary
More information1. In this process substances move from the nephron to the circulatory system. A. filtration B. reabsorption C. secretion D.
3rd lecture test 3:30 Sp 08 1. In this process substances move from the nephron to the circulatory system. A. filtration B. reabsorption C. secretion D. excretion 2. Where does secretion take place? A.
More information3. The function of that hormone. In other words, what change does that hormone facilitate.
Slide 2 The endocrine operates to regulate internal functions. It does so, via the use of hormones, or chemical messengers. Hormones travel in the blood from the site of production to distant target cells
More informationUNDERSTANDING THE DIGESTIVE SYSTEM
Slide 1 UNDERSTANDING THE DIGESTIVE SYSTEM Slide 2 Understanding Digestion The key to solving food and digestive issues starts with an understanding of how the digestive and intestinal system works If
More informationHORMONES AND CELL SIGNALLING
HORMONES AND CELL SIGNALLING TYPES OF CELL JUNCTIONS CHEMICAL SIGNALS AND MODES OF ACTION Endocrine system produces chemical messages = hormones that are transported from endocrine gland to target cell
More informationPlant organ systems, regulation and human endocrine systems: 1. Which of the following conditions would NOT cause the stomata to close?
AP Bio Practice Test Plant organ systems, regulation and human endocrine systems: 1. Which of the following conditions would NOT cause the stomata to close? a. Drought b. High temperatures c. Excessive
More informationMODELING GLUCOSE-INSULIN METABOLIC SYSTEM AND INSULIN SECRETORY ULTRADIAN OSCILLATIONS WITH EXPLICIT TIME DELAYS. Yang Kuang
MODELING GLUCOSE-INSULIN METABOLIC SYSTEM AND INSULIN SECRETORY ULTRADIAN OSCILLATIONS WITH EXPLICIT TIME DELAYS Yang Kuang (joint work with Jiaxu Li and Clinton C. Mason) Department of Mathematics and
More informationLecture 2: Glycogen metabolism (Chapter 15)
Lecture 2: Glycogen metabolism (Chapter 15) First. Fig. 15.1 Review: Animals use glycogen for ENERGY STORAGE. Glycogen is a highly-branched polymer of glucose units: Basic structure is similar to that
More informationThe Endocrine System. By Carl Antiado, Marco Habana, Cameron Loewen, and Nikki Ooka
The Endocrine System By Carl Antiado, Marco Habana, Cameron Loewen, and Nikki Ooka Standards 2.C.2 : Organisms respond to changes in their external environments. 4.A.4: Organisms exhibit complex properties
More informationStandard B-3: The student will demonstrate an understanding of the flow of energy within and between living systems.
B-3.1 Summarize the overall process by which photosynthesis converts solar energy into chemical energy and interpret the chemical equation for the process. Taxonomy Level: 2.4-B and 2.1-B Understand Conceptual
More informationMBios 478: Systems Biology and Bayesian Networks, 27 [Dr. Wyrick] Slide #1. Lecture 27: Systems Biology and Bayesian Networks
MBios 478: Systems Biology and Bayesian Networks, 27 [Dr. Wyrick] Slide #1 Lecture 27: Systems Biology and Bayesian Networks Systems Biology and Regulatory Networks o Definitions o Network motifs o Examples
More informationThyroid Gland. Chapter 18 Part 2. Thyroid gland. Thyroid Gland. Thyroid Gland. Parathyroid Gland. Adrenal Gland. Pancreas
Thyroid Gland Chapter 18 Part 2 Synthesis and function of Thyroid hormone Calcitonin and Calcium regulation Parathyroid Gland PTH and Calcium regulation Adrenal Gland The corticosteroids Pancreas Regulation
More informationControl of Blood Sugar Levels
Why? Control of Sugar Levels What hormones are involved in the homeostasis of blood sugar? All living things use as a source of energy. In vertebrates it is critical that the levels of in the blood are
More informationEndocrine System. True or False
True or False Endocrine System Some overweight people actually eat very little; their excess weight is caused by glands A woman who becomes anxious and depressed just before menstruating is suffering from
More informationPhysiology 12. Overview. The Gastrointestinal Tract. Germann Ch 19
Physiology 12 The Gastrointestinal Tract Germann Ch 19 Overview 1 Basic functions of the GI tract Digestion Secretion Absorption Motility Basic functions of the GI tract Digestion: : Dissolving and breaking
More informationType 1 Diabetes 2/23/2015. Endocrine System Hormones. Living with Type 1 Diabetes
Endocrine System Hormones 2007-2008 Living with Type 1 Diabetes Type 1 Diabetes results from the autoimmune destruction of the insulin- producing beta-cells in the pancreas. The lack of insulin leads to
More informationEquine Metabolic Syndrome and Obesity
Equine Metabolic Syndrome and Obesity By Marielle St-Laurent, DVM - Posted: June 21, 2012 to HorseJournals.com Years ago, veterinarians recognized that obese horses develop a different metabolism than
More information18. PANCREATIC FUNCTION AND METABOLISM. Pancreatic secretions ISLETS OF LANGERHANS. Insulin
18. PANCREATIC FUNCTION AND METABOLISM ISLETS OF LANGERHANS Some pancreatic functions have already been discussed in the digestion section. In this one, the emphasis will be placed on the endocrine function
More informationName Group. 8. Growth 9. Insulin 10. Iodine 11. Manufactures 12. Nervous system 13. Neuron 14. Neurotransmitter
THE NERVOUS AND ENDOCRINE SYSTEM Name Group 0. Translate all these words 1. Adrenaline 2. Brain 3. Cerebellum 4. Cerebrum 5. Dendrite 6. Dwarfism 7. Goiter 8. Growth 9. Insulin 10. Iodine 11. Manufactures
More informationCoreModels Glucose-Insulin Model A CORE LEARNING GOALS ACTIVITY FOR SCIENCE AND MATHEMATICS
CoreModels Glucose-Insulin Model A CORE LEARNING GOALS ACTIVITY FOR SCIENCE AND MATHEMATICS Summary Students use two prebuilt STELLA models to investigate the relationship between the amount of glucose
More informationHomeostasis. Endocrine System Nervous System
Homeostasis Endocrine System Nervous System 2004-2005 Regulation Why are hormones needed? chemical messages from one body part to another communication needed to coordinate whole body homeostasis & regulation
More informationHormones. Regulation. Endocrine System. What Do Endocrine Glands Release (Secrete)? Endocrine System
Endocrine System Function: Endocrine System - For glands to release products into bloodstream to deliver messages throughout the body. Regulation Is the control & coordination of other life functions to
More informationThe Cell Membrane and Homeostasis What is the cell membrane? A quick review A. Cell Membrane and Cell Transport. Unit 2: Cells and Cell Transport
Unit 2: Cells and Cell Transport Cell Membrane and Cell Transport Name: Directions: Go to https://shimkoscience.weebly.com/ and on the Biology page, find the document labelled Cell Membrane and Cell Transport
More informationEndocrine System. Chapter 9
Endocrine System Chapter 9 Endocrine Organs Hormones Chemical messengers that are released from one tissue and transported through blood to a target tissue. Chemical classification: amino acids, steroids,
More informationThe Endocrine System 2
The Endocrine System 2 Continuing on from the previous instalment, we will now look at the adrenal glands, the pancreas and the gonads as parts of the endocrine system. Adrenal Glands The adrenal glands
More informationCHAPTER 50 Endocrine Systems. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
CHAPTER 50 Endocrine Systems Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Endocrine system All the endocrine glands and other organs with hormonesecreting
More informationUnit 1 Matter & Energy for Life
Unit 1 Matter & Energy for Life Chapter 2 Interaction of Cell Structures Biology 2201 Primary Membrane Function: Homeostasis Section 2.2 Conditions in the cell must remain more or less constant under many
More informationDiabetes and Related Emergencies. *** CME Version *** Aaron J. Katz, AEMT-P, CIC
Diabetes and Related Emergencies *** CME Version *** Aaron J. Katz, AEMT-P, CIC www.es26medic.net Agenda BLS Level review of normal physiology of glucose metabolism What happens when normal becomes abnormal
More information5. Maintaining the internal environment. Homeostasis
5. Maintaining the internal environment Homeostasis Blood and tissue fluid derived from blood, flow around or close to all cells in the body. Blood and tissue fluid form the internal environment of the
More informationTarget cells only respond to specific signals Different target cells have different kinds of receptors in their membranes or cytoplasm
1 2 3 4 5 6 BIO 1102 Lec. 9 (Guided): Chapter 13 The Endocrine System The Endocrine System Controls many body functions through chemical signals Hormones Carried in blood stream Affect target cells Composed
More informationCopyright Strength Works International Publishing. All rights are reserved.
Copyright 2016. Strength Works International Publishing. All rights are reserved. BEACH WEEK Whether you have a photo shoot, a trip to the beach, a pool party, a reunion or any event that you want to turn
More informationUNDERSTANDING THE DIGESTIVE SYSTEM
UNDERSTANDING THE DIGESTIVE SYSTEM Understanding Digestion The key to solving food and digestive issues starts with an understanding of how the digestive and intestinal system works If you know what is
More informationUnit 2: Cellular Chemistry, Structure, and Physiology Module 2: Cellular Chemistry
Unit 2: Cellular Chemistry, Structure, and Physiology Module 2: Cellular Chemistry NC Essential Standard: 1.2.1 Explain how cells use buffers to regulate cell ph 4.1.1 Compare the structure and functions
More informationCIE Biology A-level Topic 14: Homeostasis
CIE Biology A-level Topic 14: Homeostasis Notes Communication is essential for the survival of organism as all living organisms must be able to detect and respond to changes in both their internal and
More informationNOTES: ENDOCRINE SYSTEM (CH 9)
NOTES: ENDOCRINE SYSTEM (CH 9) Endocrine System *The endocrine system consists of a range of glands and tissues throughout the body Functions of the Endocrine System: 1) Maintain balance within body (homeostasis)
More informationNOTES 11.5: ENDOCRINE SYSTEM. Pages
NOTES 11.5: ENDOCRINE SYSTEM Pages 1031-1042 ENDOCRINE SYSTEM Communication system that controls metabolism, growth, and development with hormones Maintains homeostasis Hormones: chemical messengers released
More informationRev. date Kaiser Foundation Health Plan of Washington
PE3620000-01-17 Rev. date 2014013 2017 Kaiser Foundation Health Plan of Washington Gestational diabetes Information to help you stay healthy during your pregnancy What is gestational diabetes? How gestational
More informationThe Endocrine System
The Endocrine System The nervous system allows the body to respond to various stimuli in a quick manner and this allows for homeostasis. The endocrine system, using hormones also allows the body to respond
More informationEndocrine System. The Endocrine Glands
Endocrine System Working together with the nervous system, the endocrine system helps maintain homeostasis in the body. Where the nervous system uses electric signals, the endocrine system uses chemical
More informationBODY CONTROL SYSTEMS
BODY CONTROL SYSTEMS THE ENDOCRINE SYSTEM - 1 of the 2 chemical control systems of the human body - function of the endocrine system: regulate body functions = maintain homeostasis ie. physical and mental
More informationCHY2026: General Biochemistry UNIT 7& 8: CARBOHYDRATE METABOLISM
CHY2026: General Biochemistry UNIT 7& 8: CARBOHYDRATE METABOLISM Metabolism Bioenergetics is the transfer and utilization of energy in biological systems The direction and extent to which a chemical reaction
More information