near-perfect adaptation in e. coli chemotaxis signal transduction network
DESCRIPTION
Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network. Yang Yang & Sima Setayeshgar. Jan, 2007. E. coli. static.howstuffworks.com/gif/cell-ecoli.gif. 1-3 microns long 1 micron in diameter 4-6 flagella Small genome (4288 genes). www.hatetank.dk. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/1.jpg)
Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network
Yang Yang & Sima Setayeshgar
Jan, 2007
![Page 2: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/2.jpg)
E. coli
•1-3 microns long
•1 micron in diameter
•4-6 flagella
• Small genome (4288 genes)
static.howstuffworks.com/gif/cell-ecoli.gif
www.hatetank.dk
![Page 3: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/3.jpg)
Bacterial Chemotaxis
http://www.rowland.harvard.edu/labs/bacteria/index_movies.html
Increasing attractants or Decreasing repellents
Run Tumble
![Page 4: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/4.jpg)
Chemotaxis signal transduction network in E. coli
well-characterized model system
CheW: Coupling CheA to MCPs CheB: CheBp demethylate MCPs CheR: Methylate MCPs CheY: flagella motor regulator protein CheZ: Dephosphoryte CheYp; CheA: Histidine kinase
![Page 5: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/5.jpg)
ligand binding
Methylation
Phosphorylation
)()( )(7/7~5/5
)( CheRLTCheRTL pnkmkkmk
pn
ppnmkmk
ppn
pnckck
pn
CheBTLCheBTL
CheRTLCheRTL
)(14~1
)(
)(14~1
)(
)()(
)()(
PCheBCheB
PCheYCheZCheZCheY
CheBCheRTCheBCheRTL
CheYCheRTCheYCheRTL
ADPCheRTLATPCheRTL
kmbp
kmyp
pnkb
np
nky
np
npkk
n
)()()(
)()()(
)()()()( 9~7
T3 T4T2
T4pT2p T3p
LT3 LT4
LT4p
LT2
LT3pLT2p
phosphorylation
methylation
Lig
and
bind
ing
Full realistic model
Chemical reactions
![Page 6: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/6.jpg)
Perfect adaptation
Steven M., et al. Journal of bacteriology. 1983
This property allows the system to compensate for the presence of continued stimulation and to be ready to respond to further stimuli
![Page 7: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/7.jpg)
Robustness
U. Alon et al. Nature,1999
inputReaction rates
proteins outpu
t
Tau-Mu Yi* et al. Biophysics,2000
Can be achieved
by
![Page 8: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/8.jpg)
Motivation
shedding light on biochemical steps and feedback mechanisms underlying robustness
shed light on values of unknown or partially known paramters
SOLVE: we develop a novel method for elucidating regions in parameter space allowing perfect adaptation.
QUESITON: basis of robustness of perfect adaptation?
![Page 9: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/9.jpg)
START with a fine-tuned model of chemotaxis network that:
reproduces key features of experiments
is NOT robust
AUGMENT the model explicitly with the requirements that:
steady state value of CheYp
values of reaction rate constants, are independent of the external stimulus, s, thereby
achieving robustness of perfect adaptation.
s
k
F
u
skuFdt
ud
0);;(
: state variables
: reaction kinetics
: reaction constants
: external stimulus
Algorithm
![Page 10: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/10.jpg)
Decretizing s
into H points
0
||
0);;(
ds
kdds
du
skuFdt
ud
N
02
|2
|
0);;(
)1(
11
11
s
kks
uu
skuFdt
ud
sjss
jm
jm
j
jN
jN
jjj
jlowj
Augmented system
The steady state concentration of proteins in the network must satisfy:
The steady state concentration of CheYp must satisfy:
At the same time, the reaction rate constants must be independent of stimulus:
: allows for near-perfect adaptation
= CheYp
0ds
kd
0);;( skuFdt
ud
N
N
u
ds
du
||
![Page 11: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/11.jpg)
Implementation
Newton-Raphson, to solve for the steady state of augmented system:
multidimensional root finding method Efficient way of converging to a root with a sufficiently good initial guess.
x
f(x) 1
2 f(x)1
x
2
Works well unfortunate case fortunate case
x
f(x) 1
32
0
||
0);;(
ds
kdds
du
skuFdt
ud
N
![Page 12: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/12.jpg)
Dsode (stiff ODE solver), to verify Time dependent behavior of proteins for different ranges of external stimulus by solving: 0);;( skuF
dt
ud
Implementation
![Page 13: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/13.jpg)
Working progress
Exploring the parameter spaces of E. coli chemotaxis signaling transduction network
Exploring the unknown parameter ranges of chemotaxis signaling transduction network of species with multiple CheYs
![Page 14: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/14.jpg)
E .coli
Relative change of CheYp:• less than 5%• less than 3%• less than 1%• pairwise trajectory
Pairwise result: 3D surface result:
parameter spaces of E. coli
![Page 15: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/15.jpg)
E. coli
Relative change of CheYp:• less than 5%• less than 3%• less than 1%• pairwise trajectory
Pairwise result: 3D surface result:
parameter spaces of E. coli
![Page 16: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/16.jpg)
Pairwise result: 3D surface result:
Relative change of CheYp:• less than 5%• less than 3%• less than 1%• pairwise
trajectory
E .coli
parameter spaces of E. coli
![Page 17: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/17.jpg)
1%
k1c : 0.17 s-1 1 s-1
k8 : 15 s-1 12.7 s-1
Violating and restoring perfect adaptation 9%
k1c : 0.17 s-1 1 s-1(1,15)
(1,12.7)
At 250s, giving step stimulus from 0 to 1e-6M
![Page 18: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/18.jpg)
Consistency with recent work done by Bernardo A. mello and Yuhai Tu
They list a series of conditions which allow near-perfect adaptations
They are a active-ependent model which the receptors are either in active or inactive state
Our parameter space remarkable shows the same consistency with their predictions about the relationships of the parameter values although we are using a active-independent model.
![Page 19: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/19.jpg)
1. The timescale for ligand binding is much shorter than the methylation and phosphorylation timescale. This condition allows us to neglect ligand-binding/unbinding kinetics.
2. The association rates between the receptor and the methylation/demethylation enzymes, CheR and CheB-P, are linearly related to the activity of the receptor and are zero for n = 4 and n = 0, respectively: and : The dissociation rates of the enzyme receptor bound states are independent of λ.
3. The receptor activities of the nonmethylated and the maximally methylated receptors are independent of l: P0v = P0o, P4v = P4o.
4. The ratios between the CheR catalytic rate kR n and the CheB-P catalytic rate of the next methylation level kB n+1 are the same for all methylation states n: kB n+1 / kR n = const:
5. The phosphate transfer rates from CheA to CheB or CheY are proportional to CheA autophosphorylation rate:
6. The explicit dependence on [TFn] distribution can be
removed from the expression
this condition can only be strictly satisfied when
nRn PPK 4 0PPK n
Bn
nPYnn
PBn PkPk ;
4
0
2 ][)][][
(n
FnnB
F
R
F
TPK
B
K
R
B
F
R
F
KB
KR ][][
List of conditions:
![Page 20: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/20.jpg)
Condition 2
nancnRn kkPPK 4
ckbkakPPK nanamnnBn 2
0 )(
![Page 21: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/21.jpg)
Condition 3 P0v = P0o, P4v =
P4o
![Page 22: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/22.jpg)
Condition 4 kB n+1 / kR n =
constan
nc
nm kkk
)1()1(
![Page 23: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/23.jpg)
Condition 5 n
PYnn
PBn PkPk ;
kb /k8-13 and ky/k8-13 are linearly related
*The parameter value are normalized to the literature value( Peter A. S., John S.P. and Hans G.O. , A model of excitation and adaptation in bacterial chemotaxis, biochemistry 1997) while the inset is not since the literature value is zero for k11.
![Page 24: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/24.jpg)
Condition 6 B
F
R
F
KB
KR ][][
![Page 25: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/25.jpg)
kb /k8-13 and ky/k8-13 are linearly related
*The parameter value are normalized to the literature value( Peter A. S., John S.P. and Hans G.O. , A model of excitation and adaptation in bacterial chemotaxis, biochemistry 1997) while the inset is not since the literature value is zero for k11.
![Page 26: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/26.jpg)
Two CheY system
•Rhodobacter sphaeroides, Caulobacter crescentus have multiple CheYs while lack of CheZ protein.•Similar chemotaxis behaviors.
![Page 27: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/27.jpg)
Two CheY system
Our work:Reproduce the key feature of chemotaxis behavior in two CheY system by replacing CheZ with CheY2.
![Page 28: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/28.jpg)
Parameter spaces of two CheY system
Introducing [CheY2] and CheY2 (de-)phosphorylation rates.Exploring the parameter values which can give perfect adaptation.
Two CheY
Relative change of CheYp:• less than 5%• less than 3%• less than 1%
Other parameter value were set as the literature value except Kb= 1e+6 M-1s-1 instead of 8e+5 M-1s-1.
![Page 29: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/29.jpg)
parameter spaces comparison of two and single CheY case
The parameter space for single CheY case seems more restrict than the two CheY case
![Page 30: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/30.jpg)
Conclusions
Successful implementation of the augmented model of the chemotaxis signal transduction network in E. coli that explicitly takes into account robust perfect adaption.
Preliminary results on projections of robustness manifolds in parameter space of E. coli and two CheY system
Complete construction of manifolds in parameter space, allowing insight into parameter dependence giving rise to robustness
Work in progress
![Page 31: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/31.jpg)
Future workApplying the method to other cellular signal
transduction networks exhibiting robust homeostasis, such as phototransduction
Signal flow in visual transduction, Leon Lagnado and Denis Baylor,Neuron,1992
![Page 32: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/32.jpg)
Thanks and comments!
![Page 33: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/33.jpg)
![Page 34: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/34.jpg)
Physics limitation in signal sensing
25 years ago Berg and Purcell had showed that the physics limitation of the single celled organism.The derivation is mainly assumed a perfect measurement device and they determined the relative measurement accuracy is :
But for multiple and noninterating receptors shaped as a ring, the formula is derived by Willam and Sima recently as:
With know parameter value, we can get the actual physics limit to measurements of CheYp concentration corresponds to
cDac
C 1 : diffusion constant
: device size: average concentration: sampling time
210 )2
1(
1
a
g
mbcDc
C
C
a
D
0g
b
m : receptor numbers: single receptor size: geometric factor of order unity
c
M
c
C 035.3
![Page 35: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/35.jpg)
*Computational models of chemotaxis signal transduction network
Activity dependent model
Barkai & Leibler (1997)
The concept of robustness in biochemical networks introduced, showing how it may arise in bacterial chemotaxis through activity-dependent kinetics. The chemoreceptor is either in active state or inactive state.Simulations show that precision of adaptation is a robust property, while adaptation time is not, and that adaptation time is inversely proportional to receptor-complex activity.It did not show how the parameter space will change which is very important for understanding the robustness mechanism.
Activity independent model
Spiro et al. (1997)
Simplified three-methylation-state model, fine-tuned by trial and error, simulates ramp, step and saturation responses to aspartate.Although it can not achieve the robust perfect adaptation, but it is a more realistic model without assuming any of the activity dependent parameter values.And our work is start from implementing this fine-tuned model.
![Page 36: Near-Perfect Adaptation in E. coli Chemotaxis Signal Transduction Network](https://reader035.vdocument.in/reader035/viewer/2022062309/56814bfc550346895db8f682/html5/thumbnails/36.jpg)
k1c/km1, k2c/km2, k3c/km3, k4c/km4 are linearly related:
*The parameter value are normalized to the literature value( Peter A. S., John S.P. and Hans G.O. , A model of excitation and adaptation in bacterial chemotaxis, biochemistry 1997).