ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Brief Report
Revised

A biochemical logarithmic sensor with broad dynamic range

[version 3; peer review: 2 approved]
PUBLISHED 18 Apr 2018
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Sensory perception often scales logarithmically with the input level. Similarly, the output response of biochemical systems sometimes scales logarithmically with the input signal that drives the system. How biochemical systems achieve logarithmic sensing remains an open puzzle. This article shows how a biochemical logarithmic sensor can be constructed from the most basic principles of chemical reactions. Assuming that reactions follow the classic Michaelis-Menten kinetics of mass action or the more generalized and commonly observed Hill equation response, the summed output of several simple reactions with different sensitivities to the input will often give an aggregate output response that logarithmically transforms the input. The logarithmic response is robust to stochastic fluctuations in parameter values. This model emphasizes the simplicity and robustness by which aggregate chemical circuits composed of sloppy components can achieve precise response characteristics. Both natural and synthetic designs gain from the power of this aggregate approach.

Keywords

Biochemical circuit, Hill equation, synthetic biology, systems biology, aggregation, robustness

Revised Amendments from Version 2

This version adds two new figures to show how the logarithmic sensor responds for various parameter combinations.

See the author's detailed response to the review by Sudin Bhattacharya
See the author's detailed response to the review by Marco Archetti

Introduction

I present a simple biochemical circuit that logarithmically transforms input signals. This circuit adds the outputs of several reactions that follow standard mass action Michaelis-Menten kinetics. Alternatively, the biochemical kinetics may follow the commonly observed Hill equation response, which includes Michaelis-Menten kinetics as a special case. This sensor has high dynamic range, responding logarithmically across many orders of magnitude. The high dynamic range is achieved by adding together reactions with different sensitivity ranges. The aggregate nature of this circuit provides robustness to parameter variations. Aggregate sensor design may explain the commonly observed high dynamic range of logarithmic biological responses and may also provide a useful tool for synthetic biology.

Results and discussion

Many biochemical reactions and cellular responses transform an input, x, into an output, y, according to the Hill equation

y=xkck+xk,
in which c is the value of the input x that yields one-half of the maximal output response, k is the Hill coefficient that determines the shape of the response, and y is normalized to be between 0 and 11. Simple mass action kinetics often follows the Hill equation with k = 1, which corresponds to classical Michaelis-Menten kinetics2. For example, the input may drive production of the output, and the output may spontaneously transform back to a prior state.

The output of the Hill equation scales approximately logarithmically with its input through the middle part of its response range, because y is roughly linear with respect to log x. Prior studies have emphasized that a Hill equation response can act as a logarithmic sensor3,4. However, a single Hill equation response provides a logarithmic sensor with limited dynamic range (Figure 1A).

8c31b728-6073-4217-adae-fa49f9c7e7d7_figure1.gif

Figure 1. The logarithmic sensor response, y, from the sum of Hill equations given in Equation 1, versus the logarithm of the input, x, on a log2 scale.

Each unit on the input scale corresponds to a doubling of the input. The range of 13 doublings is approximately four orders of magnitude, 213 ≈ 104. A logarithmic sensor responds linearly with respect to the logarithm of the input. The solid blue lines show the response, y, from Equation 1. The dashed gold lines are linear fits to the response. (A) A single Hill equation with c = 1 and k = 2. A linear response with low sensitivity occurs over a few doublings at small input levels. The following plots all use n = 7 and additional parameters as described from the summed Hill equation in Equation 1, in which ci = bi. (B) Response with k = 1 and b = 4. (C) Response with k = 2 and b = 4. (D) Response with k = 1 and b = 8. (E) Response with k = 2 and b = 8. (F) Same average parameters values of k = 2 and ci = bi for b = 4 as in the plot above, but with random parameter fluctuations around those average values. For each of the 7 Hill equations in the sum given in Equation 1, each parameter was obtained by a different random number drawn from a normal distribution. For each k, the parameters were drawn from a normal distribution with mean 2 and standard deviation 2 × 0.25 = 0.5. For each ci, the parameters were drawn from a normal distribution with mean 4i and standard deviation 4i ×0.25. The response remains nearly linear in spite of the random parameter fluctuations (see Data availability).

My extended dynamic range sensor arises by adding together n Hill equations with increasing values of the half-maximal response, ci, as

y=i=0n1xkcik+xk.(1)

For example, if ci = bi, then each reaction in the sum has an increasing input value at which its maximal sensitivity occurs. Figure 1B–E shows that simple combinations of k and b create a logarithmic sensor, in which the output is linearly related to the logarithm of the input. The logarithmic relation holds robustly when the parameters k and b of the individual reactions vary stochastically (Figure 1F). Figure 2 and Figure 3 show the response for various parameter combinations.

8c31b728-6073-4217-adae-fa49f9c7e7d7_figure2.gif

Figure 2. The logarithmic sensor response as in Figure 1, with the rows from top to bottom showing n = 1,2,3,4,5,6 and the columns from left to right showing b = 2,4,8,16.

Note that the scale along the x-axis is increased to 214 compared with 212 in Figure 1. All plots use the Michaelis-Menten relation with Hill equation coefficient k = 1 (see Data availability).

8c31b728-6073-4217-adae-fa49f9c7e7d7_figure3.gif

Figure 3. The logarithmic sensor response as in Figure 2, with Hill equation coefficient k = 2 (see Data availability).

Several biochemical circuit responses and many aspects of perception scale logarithmically5. A robust and generic pattern of this sort seems likely to depend on a robust and generic underlying design. In the search for a generic circuit, my biochemical logarithmic sensor has three advantages over prior designs. First, prior models depended on particular molecular assumptions about biochemical kinetics or reaction pathways3,4. My design requires only Michaelis-Menten or Hill equation responses, which are very widely observed in biochemical and cellular systems1. Second, prior models focused on single input-output processes, which have relatively narrow dynamic range. My aggregate design provides a logarithmic response over a vastly greater range. Third, the prior models’ responses are easily perturbed by parameter fluctuations. My design performs robustly with respect to broad fluctuations in parameters. The robustness of my logarithmic sensor emphasizes the potential to achieve precise response characteristics from underlying sloppy components when using an aggregate design6,7.

Data availability

Mathematica source code to produce Figure 1 can be found at: https://doi.org/10.5281/zenodo.12176588

License: CC BY 4.0

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 16 Feb 2018
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Frank SA. A biochemical logarithmic sensor with broad dynamic range [version 3; peer review: 2 approved]. F1000Research 2018, 7:200 (https://doi.org/10.12688/f1000research.14016.3)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 2
VERSION 2
PUBLISHED 23 Mar 2018
Revised
Views
20
Cite
Reviewer Report 11 Apr 2018
Marco Archetti, School of Biological Sciences, University of East Anglia, Norwich, UK 
Approved
VIEWS 20
This is a brilliant short paper: a clever idea, explained very clearly and that immediately makes sense (but that is by no means obvious at first sight). The revision and the comments about the first version seem appropriate to me. 
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Archetti M. Reviewer Report For: A biochemical logarithmic sensor with broad dynamic range [version 3; peer review: 2 approved]. F1000Research 2018, 7:200 (https://doi.org/10.5256/f1000research.15653.r32535)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response (F1000Research Advisory Board Member) 18 Apr 2018
    Steven Frank, University of California, Irvine, USA
    18 Apr 2018
    Author Response F1000Research Advisory Board Member
    I thank Dr. Archetti for his suggestion to illustrate the response over various values of the parameters n and b. In the new revision of the article, I have added ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response (F1000Research Advisory Board Member) 18 Apr 2018
    Steven Frank, University of California, Irvine, USA
    18 Apr 2018
    Author Response F1000Research Advisory Board Member
    I thank Dr. Archetti for his suggestion to illustrate the response over various values of the parameters n and b. In the new revision of the article, I have added ... Continue reading
Views
15
Cite
Reviewer Report 04 Apr 2018
Sudin Bhattacharya, Department of Biomedical Engineering, Department of Pharmacology & Toxicology, Institute for Quantitative Health Science and Engineering, Michigan State University, East Lansing, MI, USA 
Approved
VIEWS 15
I am satisfied with the author's response to my comments about his article, which at this ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Bhattacharya S. Reviewer Report For: A biochemical logarithmic sensor with broad dynamic range [version 3; peer review: 2 approved]. F1000Research 2018, 7:200 (https://doi.org/10.5256/f1000research.15653.r32373)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Version 1
VERSION 1
PUBLISHED 16 Feb 2018
Views
27
Cite
Reviewer Report 19 Mar 2018
Sudin Bhattacharya, Department of Biomedical Engineering, Department of Pharmacology & Toxicology, Institute for Quantitative Health Science and Engineering, Michigan State University, East Lansing, MI, USA 
Approved with Reservations
VIEWS 27
The author presents an interesting model of a robust biochemical logarithmic switch with a broad dynamical range as a sum of multiple Hill-equation based models.
  1. My main problem with this work is the lack of a
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Bhattacharya S. Reviewer Report For: A biochemical logarithmic sensor with broad dynamic range [version 3; peer review: 2 approved]. F1000Research 2018, 7:200 (https://doi.org/10.5256/f1000research.15234.r30929)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response (F1000Research Advisory Board Member) 23 Mar 2018
    Steven Frank, University of California, Irvine, USA
    23 Mar 2018
    Author Response F1000Research Advisory Board Member
    I thank Dr. Bhattacharya for sharing his expertise on topics related to the important Hill equation (see ref. 1 of the article). Dr. Bhattacharya made three helpful comments in his ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response (F1000Research Advisory Board Member) 23 Mar 2018
    Steven Frank, University of California, Irvine, USA
    23 Mar 2018
    Author Response F1000Research Advisory Board Member
    I thank Dr. Bhattacharya for sharing his expertise on topics related to the important Hill equation (see ref. 1 of the article). Dr. Bhattacharya made three helpful comments in his ... Continue reading

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 16 Feb 2018
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.