**The Michaelis-Menten Equation**

In 1903, French physical chemist Victor Henri found that enzyme reactions were initiated by a bond (more generally, a binding interaction) between the enzyme and the substrate. His work was taken up by German biochemist Leonor Michaelis and Canadian physician Maud Menten, who investigated the kinetics of an enzymatic reaction mechanism, invertase, that catalyzes the hydrolysis of sucrose into glucose and fructose. In 1913, they proposed a mathematical model of the reaction. This reaction is known as The Michaelis-Menten Equation.

Les us consider an Enzyme E is binding with a substrate S. It assumes the formation of an enzyme-substrate complex ES . After that the product P is released.

- It assumes that the ES complex is in rapid equilibrium with free enzyme E
- Breakdown of ES to form products is assumed to be slower than

(1) formation of ES and

(2) breakdown of ES to re-form E and S

This may be represented schematically as,

**E**** + ****S ****↔ ****ES → ****E**** + ****P**

Vo = k2 [ES]

Rate of ES formation = k_{1} [E][S]

= k_{1} ([Etotal] – [ES]) [S]

Rate of ES breakdown = k_{-1 }[ES] + k_{2} [ES]

**At steady state assumption:**

Rate of ES formation = Rate of ES breakdow** **

k1 ([Etotal] – [ES]) [S] = k_{-1} [ES] + k_{2} [ES]

k1 [Etotal][S] – k1[ES][S] = ( k-1 + k2 )[ES]

k1 [Etotal][S] = (k1[S] + k-1 + k2 )[ES]We know, V_{o} = k_{2 } [ES] Vo = Vmax when [Etotal] = [ES] (at saturation). Therefore Vmax = k2 [Etotal]

**Enzyme Kinetics: ****Michaelis-Menton** **Equation**

## 1.Characteristics of K_{m}:

- K
_{m}—the Michaelis constant—is characteristic of an enzyme and its particular substrate, and reflects the affinity of the enzyme for that substrate. - K
_{m}does not vary with the concentration of enzyme - Small K
_{m}:

A numerically small (low) K_{m}reflects a high affinity of the enzyme for substrate. - Large K
_{m}:

A numerically large (high) K_{m}reflects a low affinity of enzyme for substrat

## 2. Relationship of velocity to enzyme concentration:

- The rate of the reaction is directly proportional to the enzyme concentration at all substrate concentrations.For example, if the enzyme concentration is halved, the initial rate of the reaction (v
_{o}), as well as that of Vmax, are reduced to half that of the original.

MATHEMATICAL INVERSE FUNCTION (EQUATION) FOR ENZYME KINETICS

ABSTRACT

The most significant application of Lineweaver–Burk plot lies in the determination of Michaelis constant (Km), substrate concentration [S] at which the velocity (v) of the enzyme catalyzed biochemical reaction attain half of it’s maximum (Vmax ÷ 2). For practical purposes, Km is the concentration of substrate which allows the enzyme velocity to achieve half of it’s maximum Vmax (Vmax ÷ 2). Most points on the plot are found far to the right of the y-axis. Through the reverse the mathematical steps and get inverse of substrate concentration (1÷S) back from some output value, say inverse of respective velocity (1÷v), one should undo each step in reverse order. It means that, one should subtract the inverse of maximum velocity (1÷Vmax) from inverse of respective velocity (1÷v) and then multiply the result by Vmax/Km . This is going to yield the equation correspond to: 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) . The 1÷S and 1÷v for given enzyme catalyzed biochemical reaction deserve symmetry, that is to say the symmetry between a Lineweaver Burk Plot (the real function) and the inverse function for enzyme kinetics of present attempt. The co-ordinates of the point of intersection of both the equations 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) and 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) correspond to: ( 1/(Vmax-Km ) , 1/(Vmax-Km ) ) .

Key Words: Baramati Constant; 1/(Vmax-Km ) ; mathematical approach

1. Seema Karna Dongare 2. Manali Rameshrao Shinde 3.Vitthalrao Bhimasha Khyade

Address:

1 and 2: Student, T. Y. B. Sci (Microbiology – 2018 – 2019), and Science Association, Shardabai Pawar Mahila Mahavidyalaya, Shardanagar Tal. Baramati Dist. Pune – 413115 (India).

3 (A): Sericulture Unit, Malegaon Sheti Farm, Agricultural Development Trust Baramati, Shardanagar, (Malegaon Khurd) Post Box No – 35, Baramati, Pune 413 115, Maharashtra, India

[email protected]

3 (B). Head, Department of Zoology, Shardabai Pawar Mahila Mahavidyalaya, Shardanagar Tal. Baramati Dist. Pune – 413115 (India). [email protected]

3 (C). Dr. APIS”, Shrikrupa Residence, Teachers Society, Malegaon Colony (Baramati) Dist. Pune – 413115. [email protected]

INTRODUCTION

The velocity (v) of biochemical reaction catalyzed by the enzyme vary according to the status of factors like: temperature; pH; enzyme concentration; substrate concentration; activators and inhibitors. There is no linear response of velocity (v) of enzyme catalyzed biochemical reaction to the substrate concentration. This may be due to saturable nature of enzyme catalyzed biochemical reactions. If the initial velocity (v) or rate of the enzyme catalyzed biochemical reaction is measured over a range of concentration of substrate [S], it appears to increase. That is to say, initial velocity (v) of the enzyme catalyzed biochemical reaction get increase according to the increase in the concentration of substrate [S]. This tendency of increase in initial velocity (v) of the enzyme catalyzed biochemical reaction according to the increase in the concentration of substrate [S] is observed upto certain level of the concentration of substrate [S]. At this concentration of substrate [S], the enzyme exhibit saturation and exert the initial velocity (v) of the enzyme catalyzed biochemical reaction to attain it’s maximum (Vmax). Hans Lineweaver and Dean Burk (1934) suggested the double reciprocal plot for presenting the data on substrate [S] and velocity (v) of the enzyme catalyzed biochemical reaction. In enzyme kinetics, double reciprocal plot suggested by Hans Lineweaver and Dean Burk is well esteemed graphical presentation of the data on substrate concentration [S] and velocity (v) of the enzyme catalyzed biochemical reaction recognized as, the “Lineweaver–Burk plot”. This plot deserve wide applicability. The most significant application of Lineweaver–Burk plot lies in the determination of substrate concentration [S] at which the velocity (v) of the enzyme catalyzed biochemical reaction attain half of it’s maximum (Vmax ÷ 2). According to Lineweaver and Dean Burk (1934), substrate concentration [S] at which the velocity (v) of the enzyme catalyzed biochemical reaction attain half of it’s maximum (Vmax ÷ 2) is called as Km (Michaelis constant). For practical purposes, Km is the concentration of substrate which allows the enzyme velocity to achieve half of it’s maximum Vmax (Vmax ÷ 2). An enzyme with a high Km has a low affinity for its substrate, and requires a greater concentration of substrate to achieve maximum Vmax (Keith J. Laidler, 1997).

The Lineweaver–Burk plot was widely used to determine important terms in enzyme kinetics, such as Km and Vmax, before the wide availability of powerful computers and non-linear regression software. The y-intercept of such a graph is equivalent to the inverse of Vmax; the x-intercept of the graph represents −1/Km. It also gives a quick, visual impression of the different forms of enzyme inhibition. Accordingly, the mathematical equation suggested by Lineweaver and Dean Burk (1934) can be written as: 1/(V ) =(Km )/Vmax X 1/S+1/(Vmax )

This Lineweaver–Burk plot deserve wide applicability. It is useful for the determination of Km , the most significant factor in enzyme kinetics. The intercept on y – axis of Lineweaver–Burk plot is the inverse of Vmax that is (1/ Vmax ). The intercept on X – axis of Lineweaver–Burk plot is the inverse of – Km that is (−1/Km ). the x-intercept of the graph represents −1/Km. This Lineweaver–Burk plot (double reciprocal plot) also gives a quick, visual impression of the different forms of enzyme activation and inhibition.

Hayakawa, et al (2006) reported the distortion of the error structure through the Lineweaver–Burk plot (double reciprocal plot). It is therefore, the Lineweaver–Burk plot (double reciprocal plot) appears to be for the determination of enzyme kinetic parameters. Although it is still used for representation of kinetic data, (Hayakawa, et al ,2006) non-linear regression or alternative linear forms of the Michaelis–Menten equation such as the Hanes-Woolf plot or Eadie–Hofstee plot are generally used for the calculation of parameters (Greco and Hakala, 1979).

Dick (2011) explained type of enzyme inhibition. Accordingly, the Lineweaver–Burk plot (double reciprocal plot) can distinguish different types of inhibitors, such as competitive, non-competitive and uncompetitive inhibitors. The competitive inhibitors have the same y-intercept (1/ Vmax ) as uninhibited enzyme (since Vmax is unaffected by competitive inhibitors the inverse of Vmax also doesn’t change) but there are different slopes and x-intercepts between the two data sets. The non-competitive inhibition produces plots with the same x-intercept (−1/Km ) as uninhibited enzyme (Km is unaffected) but different slopes and y-intercepts. Uncompetitive inhibition causes different intercepts on both the y- and x-axes (Berg, et al, 2002). John E. Dowd and Douglas Briggs (1965) reviewed the literature on “Estimates of Michaelis – Menten kinetic constants through the use of different linear transformation” and listed some problems with Lineweaver–Burk plot (double reciprocal plot). Accordingly, Lineweaver–Burk plot (double reciprocal plot) is classically used in older texts. It seems to be prone to the error. The y-axis takes the reciprocal of the rate of reaction – in turn increasing any small errors in measurement. Also, most points on the plot are found far to the right of the y-axis. This may be due to limiting solubility. It is not allowing for large values of 1/[S] and hence no small values for [S]), calling for a large extrapolation back to obtain x- and y-intercepts.

MATHEMATICAL INVERSE FUNCTION FOR ENZYME KINETICS

According to Hall, Arthur Graham and Frink, Fred Goodrich (1909), in mathematics, the equation (mathematical function) is represented in two forms: Real form and Inverse form. The inverse function may also called as anti – function. The inverse function is a function that “reverses” real form of function. If the function f applied to an input x gives a result of y, then applying its inverse function g to y gives the result x, and vice versa, i.e., f(x) = y if and only if g(y) = x (Scheinerman, Edward R., 2013). For the purpose of fortified understanding of the concept of “Inverse Function”, mathematics is citing simple example. Let us have a look on this example. Consider the real-valued function of a real variable given by f(x) = 5x – 7 or y = 5x – 7 (In “y = m x + c” format. This is because, the original Lineweaver Burk plot deserve this format). Let us take a number x, multiply it by 5, then subtract 7 from the result, to reverse this and get x back from some output value, say y, one should undo each step in reverse order. It means that one should add 7 to y and then divide the result by 5. In functional notation this inverse function would be given by: g(x) = (1÷5)x + (7÷5) or y = (1÷5)x + (7÷5).

Let us now have directly the equation for Lineweaver Burk plot: 1/(V ) = (Km )/Vmax X 1/S+1/(Vmax ) . The mathematical equation of Lineweaver Burk plot, itself is self explainatory. Inverse of substrate concentration (1÷S) is to be multiplied with the slope, that is Km÷Vmax. Finally, it is necessary to add the inverse of maximum velocity (1÷Vmax) to get the value of inverse of respective velocity (1÷v) of enzyme catalyzed biochemical reaction. Of course, this is not the new idea of the present attempt. It has already been expected in the Lineweaver Burk plot of enzyme kinetics. The efforts of present attempt are concerned with getting the inverse function (equation) for Lineweaver Burk plot of enzyme kinetics.

The mathematical attempt for obtaining the inverse function for “y = m x + c” is concerned with getting the value of “X” back. To reverse the mathematical steps and get inverse of substrate concentration (1÷S) back from some output value, say inverse of respective velocity (1÷v), one should undo each step in reverse order. It means that, one should subtract the inverse of maximum velocity (1÷Vmax) from inverse of respective velocity (1÷v) and and then multiply the result by Vmax/Km . This is going to yield the equation correspond to: 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) .

For graphical presentation of equations of regular Lineweaver Burk plot : 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ;

it is necessary to take the inverse of substrate concentration (1÷S) on X – axis and inverse of respective velocity (1÷v) of enzyme catalyzed biochemical reaction on Y – axis. This is going to serve the purpose in expectation of Lineweaver Burk plot: 1/(V ) = (Km )/Vmax X 1/S+1/(Vmax ) .

For graphical presentation of equation : 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) ;

it is necessary to take the inverse of respective velocity (1÷v) of enzyme catalyzed biochemical reaction on X – axis and the inverse of substrate concentration (1÷S) on Y – axis. This is going to serve the purpose in expectation of Inverse function for Lineweaver Burk plot. That is to say plotting the readings of (1÷S) against (1÷v) and plotting the readings of (1÷v) against (1÷S) (Fig.1).

PROPERTIES OF MATHEMATICAL INVERSE FUNCTION FOR ENZYME KINETICS:

(A). The mathematical equation for regular Lineweaver Burk plot is explaining binary relation between the inverse of substrate concentration (1÷S) and inverse of respective velocity (1÷v) of enzyme catalyzed biochemical reaction. It is as good as a special type of mathematical binary relation. It means, many of the properties of inverse function for enzyme kinetics correspond to properties of converse relations.

(B). Theoretically, mathematical equation for regular Lineweaver Burk plot for a given enzyme is with unique inverse function. That is to say, for each concentration of substrate, there is a unique value for velocity of enzyme catalyzed biochemical reaction. This follows since the inverse function must be the converse relation.

(C). If the 1÷S and 1÷v for given enzyme catalyzed biochemical reaction are inverses of each other, then the domain of 1÷S is equal to the range of 1÷v and the range of 1÷v is equal to the domain of 1÷S.

(D).The 1÷S and 1÷v for given enzyme catalyzed biochemical reaction deserve symmetry, then there is a symmetry between a function and its inverse. This statement is a consequence of the implication that for 1÷v to be invertible it must be bijective.

(E).The co-ordinates of the point of intersection of both the equations 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) and 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) correspond to: ( 1/(Vmax-Km ) , 1/(Vmax-Km ) ) .

(F). The 1÷S and 1÷v for given enzyme catalyzed biochemical reaction may deserve “one-to-one”. Even if a function 1÷S is not one-to-one relation with 1÷v ; it is possible to define a inverse through restricting the domain.

(G). The inverse function in enzyme kinetics is a function that reverses real form of mathematical equation for regular Lineweaver Burk plot. If the equation for (1÷v), applied to input (1÷ S), gives a result of (1÷v). Then applying its inverse function (1÷S) to (1÷v) gives the result 1÷S, and vice versa.

(H). If ( 1/(S ) , 1/(v ) ) is a point on the graph of the original equation of Lineweaver–Burk plot, then the point ( 1/(S ) , 1/(v ) ) must be a point on the graph of the inverse function of enzyme kinetics. The Lineweaver–Burk plot and it’s inverse function are mirror images of each other with respect to the line y=x .

APPLICATIONS OF MATHEMATICAL INVERSE FUNCTION (EQUATION) FOR ENZYME KINETICS:

Applications of Interfacial enzyme catalyzed reactions are ubiquitous both in vivo and in technical applications. Analysis of kinetics of enzyme catalyzed reactions remains controversial. Michaelis–Menten equation expect reaching a steady-state condition in the opposite experimental limit, where the concentration of enzyme far exceeded the molar concentration of accessible surface sites.In particular, it is unclear whether conventional Michaelis–Menten theory. It requires a large excess of substrate. Extensive experimental study of the enzymatic hydrolysis of insoluble cellulose indeed proving that the conventional approach of Michaelis–Menten equation had a limited applicability. The interfacial enzyme catalysis may reach a steady-state condition in the opposite experimental limit. In such situation, the concentration of enzyme far exceeded the molar concentration of accessible surface sites. Under such condition, an “inverse function or equation for enzyme kinetics” deserve applicability. Here it is essential to mention that, “inverse function or equation for enzyme kinetics” provides a general tool for kinetic analyses of interfacial enzyme reactions. It may establish theory provides a bridge to the accumulated understanding of steady-state enzyme kinetics. Itreveals the density of enzyme attack sites on the substrate surface as probed by one specific enzyme.

The regular Michaelis–Menten equation ( 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ) explains the influence of concentration of substrate [S] on velocity of enzyme catalyzed biochemical reaction. It’s inverse function may explain the influence of velocity of enzyme catalyzed biochemical reaction (v) on concentration of substrate [S] . That is to say, inverse function of the regular Michaelis–Menten equation and ( 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) ) is going to explain the role of product of enzyme catalyzed biochemical reaction in controlling the rate of reaction. The regular Michaelis–Menten equation ( 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ) is demonstrating the Km (Michaelis constant), the substrate concentration [S] at which the velocity (v) of the enzyme catalyzed biochemical reaction attain half of it’s maximum (Vmax ÷ 2). And … and … the inverse function is demonstrating the [1÷ (Vmax. – Km)], (Baramati Constant), point on both, the regular Michaelis–Menten equation ( 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ) and it’s inverse function ( 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ). At this point [1÷ (Vmax. – Km)], (Baramati Constant), both the equations are equal with each other. This point [1÷ (Vmax. – Km)], (Baramati Constant) is going to serve the saturation of enzyme molecules and the substarte molecules in enzyme catalyzed biochemical reaction.

CONCLUSION:

Through the reverse the mathematical steps and to get inverse of substrate concentration (1÷S) back from some output value, say inverse of respective velocity (1÷v), one should undo each step in reverse order. It means that, one should subtract the inverse of maximum velocity (1÷Vmax) from inverse of respective velocity (1÷v) and and then multiply the result by Vmax/Km . This is going to yield the equation correspond to: 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) . The 1÷S and 1÷v for given enzyme catalyzed biochemical reaction deserve symmetry, then there is a symmetry between a function and its inverse. The co-ordinates of the point of intersection of both (Real form of Lineweaver Burk plot and it’ Inverse form) the equations 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) and 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) correspond to: ( 1/(Vmax-Km ) , 1/(Vmax-Km ) ) . The inverse function of the regular Michaelis–Menten equation and ( 1/(S ) = (Vmax )/Km X 1/v – 1/(Km ) ) is going to explain the role of product of enzyme catalyzed biochemical reaction in controlling the rate of reaction. The [1÷ (Vmax. – Km)], (Baramati Constant), point on both, the regular Michaelis–Menten equation ( 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ) and it’s inverse function ( 1/(V ) = (Km )/( Vmax) X 1/S+1/(Vmax ) ). At this point [1÷ (Vmax. – Km)], (Baramati Constant), both the equations are equal with each other. This point [1÷ (Vmax. – Km)], (Baramati Constant) is going to serve the saturation of enzyme molecules and the substarte molecules in enzyme catalyzed biochemical reaction. The attempt on the inverse function for enzyme kinetics of present attempt may open a new chapter to classify the enzymes on the basis of mathematical approach.

ACKNOWLEDGEMENTS: Expertise support from Department of Microbiology, Shardabai Pawar Mahila Mahavidyalaya, Shardanagar, Tal. Baramati, Dist. Pune deserve appreciations and exert a grand salutary influence. The fifth virtual science congress 2018 of I S C A availed excellent opportunity of presentation of the work.

REFERENCES:

Lineweaver, H; Burk, D. (1934). “The determination of enzyme dissociation constants”. Journal of the American Chemical Society. 56 (3): 658–666. doi:10.1021/ja01318a036.

Hayakawa, K.; Guo, L.; Terentyeva, E.A.; Li, X.K.; Kimura, H.; Hirano, M.; Yoshikawa, K.; Nagamine, T.; et al. (2006). “Determination of specific activities and kinetic constants of biotinidase and lipoamidase in LEW rat and Lactobacillus casei (Shirota)”. Journal of Chromatography B. 844 (2): 240–50. doi:10.1016/j.jchromb.2006.07.006. PMID 16876490.

Greco, W.R.; Hakala, M.T. (1979). “Evaluation of methods for estimating the dissociation constant of tight binding enzyme inhibitors” (PDF). Journal of Biological Chemistry. 254 (23): 12104–12109. PMID 500698.

Dowd, John E.; Riggs, Douglas S. (1965). “A comparison of estimates of Michaelis–Menten kinetic constants from various linear transformations” (pdf). Journal of Biological Chemistry. 240 (2): 863–869.

Keith J. Laidler (1997). New Beer in an Old Bottle: Eduard Buchner and the Growth of Biochemical Knowledge, pp. 127–133, ed. A. Cornish-Bowden, Universitat de València, Valencia, Spain, 1997 http://bip.cnrs-mrs.fr/bip10/newbeer/laidler.pdf

Dick RM (2011). “Chapter 2. Pharmacodynamics: The Study of Drug Action”. In Ouellette R, Joyce JA. Pharmacology for Nurse Anesthesiology. Jones & Bartlett Learning. ISBN 978-0-7637-8607-6.

Berg, Jeremy M.; Tymoczko, John L.; Stryer, Lubert (2002). “The Michaelis-Menten Model Accounts for the Kinetic Properties of Many Enzymes”. Biochemistry. 5th edition. W H Freeman.

John E. Dowd and Douglas Briggs (1965). Estimates of Michaelis – Menten kinetic constants from various linear transformation. The Journal of Biological Chemistry Vol. 240 No. 2 February, 1965. http://www.jbc.org/content/240/2/863.full.pdf

Hall, Arthur Graham; Frink, Fred Goodrich (January 1909). “Chapter II. The Acute Angle [14] Inverse trigonometric functions”. Written at Ann Arbor, Michigan, USA. Trigonometry. Part I: Plane Trigonometry. New York, USA: Henry Holt and Company / Norwood Press / J. S. Cushing Co. – Berwick & Smith Co., Norwood, Massachusetts, USA. p. 15. Retrieved 2017-08-12. https://en.wikipedia.org/wiki/Inverse_function

Scheinerman, Edward R. (2013). Mathematics: A Discrete Introduction. Brooks/Cole. p. 173. ISBN 978-0840049421.

Hall, Arthur Graham; Frink, Fred Goodrich (January 1909). “Chapter II. The Acute Angle [14] Inverse trigonometric functions”. Written at Ann Arbor, Michigan, USA. Trigonometry. Part I: Plane Trigonometry. New York, USA: Henry Holt and Company / Norwood Press / J. S. Cushing Co. – Berwick & Smith Co., Norwood, Massachusetts, USA. p. 15. • • Keisler, Howard Jerome. “Differentiation” (PDF). Retrieved 2015-01-24. §2.4

Scheinerman, Edward R. (2013). Mathematics: A Discrete Introduction. Brooks/Cole. p. 173. ISBN 978-0840049421.

Korn, Grandino Arthur; Korn, Theresa M. (2000) [1961]. “21.2.-4. Inverse Trigonometric Functions”. Mathematical handbook for scientists and engineers: Definitions, theorems, and formulars for reference and review (3 ed.). Mineola, New York, USA: Dover Publications, Inc. p. 811;. ISBN 978-0-486-41147-7.

Oldham, Keith B.; Myland, Jan C.; Spanier, Jerome (2009) [1987]. An Atlas of Functions: with Equator, the Atlas Function Calculator (2 ed.). Springer Science+Business Media, LLC. doi:10.1007/978-0-387-48807-3. ISBN 978-0-387-48806-6. LCCN 2008937525.

Briggs, William; Cochran, Lyle (2011). Calculus / Early Transcendentals Single Variable. Addison-Wesley. ISBN 978-0-321-66414-3.

Devlin, Keith J. (2004). Sets, Functions, and Logic / An Introduction to Abstract Mathematics (3rd ed.). Chapman & Hall / CRC Mathematics. ISBN 978-1-58488-449-1.

Fletcher, Peter; Patty, C. Wayne (1988). Foundations of Higher Mathematics. PWS-Kent. ISBN 0-87150-164-3.

Lay, Steven R. (2006). Analysis / With an Introduction to Proof (4 ed.). Pearson / Prentice Hall. ISBN 978-0-13-148101-5.

Smith, Douglas; Eggen, Maurice; St. Andre, Richard (2006). A Transition to Advanced Mathematics (6 ed.). Thompson Brooks/Cole. ISBN 978-0-534-39900-9.

Thomas, Jr., George Brinton (1972). Calculus and Analytic Geometry Part 1: Functions of One Variable and Analytic Geometry (Alternate ed.). Addison-Wesley.

Wolf, Robert S. (1998). Proof, Logic, and Conjecture / The Mathematician’s Toolbox. W. H. Freeman and Co. ISBN 978-0-7167-3050-7.

Stryer L, Berg JM, Tymoczko JL (2002). Biochemistry (5th ed.). San Francisco: W.H.Freeman. ISBN 0-7167-4955-6.

Schomburg I, Chang A, Placzek S, Söhngen C, Rother M, Lang M, Munaretto C, Ulas S, Stelzer M, Grote A, Scheer M, Schomburg D (January 2013). “BRENDA in 2013: integrated reactions, kinetic data, enzyme function data, improved disease classification: new options and contents in BRENDA”. Nucleic Acids Research. 41 (Database issue): D764–72. doi:10.1093/nar/gks1049. PMC 3531171. PMID 23203881

Radzicka A, Wolfenden R (January 1995). “A proficient enzyme”. Science. 267 (5194): 90–931. Bibcode:1995Sci…267…90R. doi:10.1126/science.7809611. PMID 7809611.

Callahan BP, Miller BG (December 2007). “OMP decarboxylase—An enigma persists”. Bioorganic Chemistry. 35 (6): 465–9. doi:10.1016/j.bioorg.2007.07.004. PMID 17889251.

Michaelis L, Menten M (1913). “Die Kinetik der Invertinwirkung” [The Kinetics of Invertase Action]. Biochem. Z. (in German). 49: 333–369.; Michaelis L, Menten ML, Johnson KA, Goody RS (2011). “The original Michaelis constant: translation of the 1913 Michaelis-Menten paper”. Biochemistry. 50 (39): 8264– 9. doi:10.1021/bi201284u. PMC 3381512. PMID 21888353.

Briggs GE, Haldane JB (1925). “A Note on the Kinetics of Enzyme Action”. The Biochemical Journal. 19 (2): 339–339. doi:10.1042/bj0190338. PMC 1259181. PMID 16743508.

Ellis RJ (October 2001). “Macromolecular crowding: obvious but underappreciated”. Trends in Biochemical Sciences. 26 (10): 597–604. doi:10.1016/S0968- 0004(01)01938-7. PMID 11590012.

Kopelman R (September 1988). “Fractal reaction kinetics”. Science. 241 (4873): 1620

Bibcode:1988Sci…241.1620K. doi:10.1126/science.241.4873.1620. PMID 17820893.

Segel, L.A.; Slemrod, M. (1989). “The quasi-steady-state assumption: A case study in perturbation”. ThermochimActa31 (3): 446–477. doi:10.1137/1031091.

Leskovac, V. (2003). Comprehensive enzyme kinetics. New York: Kluwer Academic/Plenum Pub. ISBN 978-0-306-46712-7.

Greco, W.R.; Hakala, M.T. (1979). “Evaluation of methods for estimating the dissociation constant of tight binding enzyme inhibitors,”.J BiolChem254 (23): 12104–12109. PMID 500698.

Hayakawa, K.; Guo, L.; Terentyeva, E.A.; Li, X.K.; Kimura, H.; Hirano, M.; Yoshikawa,K.; Nagamine, T. et al. (2006). “Determination of specific activities and kinetic constants of biotinidase and lipoamidase in LEW rat and Lactobacillus casei (Shirota)”. J Chromatogr B AnalytTechnol Biomed Life Sci844 (2): 240–50.doi:10.1016/j.jchromb.2006.07.006. PMID 16876490.

Vitthalrao B. Khyade; Vivekanand V. Khyade and Sunanda V. Khyade (2017). The Baramati Quotient for the accuracy in calculation of Km value for Enzyme Kinetics. International Academic Journal of Innovative Research Vol. 4, No. 1, 2017, pp. 1-22.ISSN 2454-390X