跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
自然。作者手稿;2018年4月16日PMC发布。
以最终编辑形式发布为:
2017年10月25日在线发布。 数字对象标识:10.1038/自然24299
预防性维修识别码:项目编号5901684
NIHMSID公司:NIHMS956215
PMID:29072300

全球资源分配策略控制大肠杆菌

关联数据

补充资料
数据可用性声明

摘要

系统生物学的一个重大挑战是从潜在的分子相互作用出发,预测生命系统对扰动的动力学响应。长期以来,营养环境的变化一直被用于研究微生物的调节和适应现象1它们仍然是一个积极调查的主题411虽然人们对控制关键代谢过程对应用扰动的调节的分子相互作用知之甚多1217,对于自下而上的预测建模,它们的量化不够。在这里,我们开发了一种自顶向下的方法,扩展了最近建立的粗粒度蛋白质组分配模型15,1820从稳态增长到动力学状态。仅使用基本调控过程的定性知识并施加通量平衡条件,我们导出了一个独立于无法访问的动力学参数的细菌生长转变定量模型。由此产生的通量控制调节模型准确地预测了基因表达和生物量积累的时间进程,以响应碳向上转移和向下转移(例如,双向转移),而无需调整参数。正如模型预测和定量蛋白质组学验证的那样,由于严格的蛋白质合成分配策略,细胞在营养物质转移时表现出次优的恢复动力学,这并不是为了缓解特定的代谢瓶颈。我们的方法不依赖于动力学参数,因此指出了一个理论框架,用于描述广泛的此类动力学过程,而无需详细了解潜在的生化反应。

我们专注于可共同利用的碳基质之间的生长转变21或那些密切相关且易于切换的,避免长时间生长停滞或异质生长特征8,9,11甚至稳态行为的特征也很差。我们首先通过种植野生型植物来研究营养物质的迁移大肠杆菌K-12菌株(NCM3722)在最小培养基中,以琥珀酸盐为唯一的碳底物。在指数增长期间,第二个共同利用的碳基质21添加(葡萄糖酸盐)(图1和方法)。600 nm(OD)处光密度的时间过程600纳米)细胞培养(细胞质量M的代表,扩展数据图1a),在整个过渡期间进行了测量(图1a); 转变后不久,增长速度加快。核糖体含量(以RNA丰度报告16)众所周知,其增长速度超过细胞质量,而分解代谢蛋白含量(通过在含有异丙基β-D-1-硫代吡喃半乳糖苷、IPTG的培养基中诱导LacZ表达来报告15)花了几个小时才适应。最终,所有这些数量都以相同的速度增长,这是平衡增长所必需的。瞬时增长率,λ=dlnM(M)/d日,慢慢收敛到其最终值λ(f)经过几个小时(图1b、和扩展数据图2a时间更长)。生物质合成通量,J型=天M(M)/d日,反应迅速:在升档后30分钟内达到最终增加速度(图1c). 新RNA和LacZ的合成也观察到类似的快速反应(图1d,e),在升档后一个质量内接近最终值(虚线的斜率)的表达水平(数据的斜率)加倍。

保存图片、插图等的外部文件。对象名称为nihms956215f11.jpg
营养升档和降档

大肠杆菌在最小培养基中生长的K-12(和诱导LacZ的IPTG)通过添加第二种营养素或耗尽营养素而发生移位。a–e,添加第二种营养素后上移:20 mM琥珀酸盐(持续存在),20 mM葡萄糖酸盐(仅在= 0).f–j营养素耗尽后降档:20 mM琥珀酸(持续存在),0.03%葡萄糖(约耗尽= 0).a、 (f)、光密度(OD600纳米,红色圆圈),RNA丰度(每培养体积的RNA数量,蓝色倒三角形),以及每培养体积LacZ活性,均与=0(绿色三角形),绘制时间。虚线的斜率表示最终增长率(升档,λ(f)=0.98小时−1; 降档,λ(f)=0.45小时−1).b、 克,瞬时增长率λ()(红色钻石)。红色虚线表示λ(f).c、 小时,OD的时间导数600纳米(生物量通量,红色方块)。虚线的斜率表示λ(f).d、 e、i、j、绝对LacZ活性和RNA丰度与OD的关系图600纳米绿色和蓝色线条的斜率表示最终表达水平(定义为总蛋白的分数,参见补充注释2中的补充方程(8)); 橙色线的斜率表示转移后的预测表达水平。数据具有高度可重复性,请参阅扩展数据图2在所有面板中,黑色曲线是FCR模型的无参数拟合预测,按照方法中的概述获得。

接下来,我们描述了由于培养基中营养成分的消耗而导致的营养物质向下转移的动力学,如莫诺德的双倍生长研究中所述1野生型细胞在饱和的琥珀酸盐和少量葡萄糖上生长(图1f):葡萄糖与琥珀酸共同利用21,细胞首先使用两种碳底物,并在葡萄糖耗尽后切换到只使用琥珀酸。此时,瞬时增长率λ()急剧下降(图1g). 在恢复过程中,RNA和LacZ(图1i,j)分别显示表达水平下降和增加(数据斜率),最终达到最终水平(虚线斜率)。升档和降档都是高度可重复的(扩展数据图2b)在其他八对衬底上也观察到了类似的行为(扩展数据图3)和其他四对基板用于降档(扩展数据图4),涵盖了广泛的初始和最终增长率组合。

中显示的恢复过程图1受监管互动的控制。尽管有丰富的分子知识,特别是关于四磷酸鸟苷(ppGpp)对核糖体合成的调节13,14,16,17以及cAMP-cAMP受体蛋白(cAMP-CRP)信号系统对碳分解代谢的调节12,15,缺乏对这些调控过程动态的定量描述。不包括基因调控动力学的模型无法一致地描述不同养分转移后的不同类型的恢复动力学(补充说明1扩展数据图5). 在这里,我们介绍了一个粗粒度动力学模型,该模型定量地纳入了分解代谢和蛋白质合成的调控,尽管缺乏描述分子过程的动力学参数。我们假设生长恢复的动力学受到许多关键内部代谢物(包括酮酸和氨基酸)丰度的限制11,22,我们将其粗粒化为一个“中央前体”池,因为这些代谢物通过快速可逆反应连接起来23。如所示图2,这个前驱物池被碳流入填满J型C类(方程式(1)in图2)在蛋白质合成过程中被核糖体耗尽。因为核糖体的翻译活性,σ(),数学上定义为平均翻译率24,由前体库设置,我们假设它在粗粒度时间尺度上进行调整,以便蛋白质合成的流量J型R(右)平衡碳流入J型C类(方程式(2a)和(2b)in图2).

保存图片、插图等的外部文件。对象名为nihms956215f12.jpg
FCR模型

流量平衡和蛋白质合成分配。碳底物S1和S2通过摄取丰富的Cat1和Cat2蛋白导入M(M)类别1M(M)第2类按费率k个1k个2分别是。这些底物外部浓度的变化导致碳流入的变化J型C类(红色,方程式(1)),提供了一个中央前体池(正方形)。这些前体被核糖体(Rb)消耗用于蛋白质合成,其通量J型R(右)(橙色,方程式(2b),带α作为从碳到蛋白质的转换因子)必须在粗粒度时间尺度上平衡碳流入的变化。对于给定数量的核糖体(丰度)达到这种蛋白质合成通量M(M)Rb公司)需要改变翻译活动σ如方程式(2a)所定义。分子上σ归因于前体库的变化16,24部分蛋白质合成通量分配给分解代谢和核糖体蛋白质的表达,M(M)Rb公司M(M)类别1、2(方程式(3a)和(4a))。该分配由全球监管职能部门决定χ卢比()和χ猫(),由前体池通过ppGpp进行分子设置(蓝色箭头)13,16和cAMP–CRP系统(绿色箭头)15以及特定于基层的法规小时1,2(补充说明3).b条、监管职能的形式χ^Rb公司(σ)χ^(σ)(黑线)可以从已知的稳态测量值中确定(扩展数据图1c、d)导出于补充说明2.2d–f模型的中心简化假设是,生长过渡期间调节函数的时滞依赖性仅通过σ,如方程式(3b)和(4b)中明确表示的。

蛋白质合成的分配由基因调节功能决定:前体库的增加上调了核糖体的合成(蓝色,方程式(3a))图2)通过函数χRb公司(),由ppGpp信号控制13,14,16,17同时,分解代谢蛋白的合成(绿色,方程式(4a)图2)被压抑通过χ()并由cAMP–CRP信号系统控制12,15因此,这两种调节功能都依赖于前体库,前体库也设置了翻译活动σ.这些观察结果使我们形成了我们的关键方法,它关闭了模型方程并使其独立于分子细节:我们使用翻译活性σ()设置监管职能χRb公司()和χ()通过稳态关系,χ^Rb公司(σ)χ^(σ)。后者如所示图2b并从前面描述的“增长规律”中获得15,18(扩展数据图1c–g). 这在中的方程式(3b)和(4b)中明确说明图2相当于准静态假设(补充注释3和4). 通过构造,得到的粗粒度模型并没有试图解决非常快速的动力学过程,例如,向富介质迁移期间的核糖体振荡25然而,它描述了细胞生长和基因表达在较慢的时间尺度上的动力学,例如从大约10分钟到达到新的稳定状态,前提是生长过渡不涉及长时间的生长停滞8,这将使准静态近似无效。

我们将这种基因调控模型称为通量控制调控(FCR),因为关键的动态变量σ()通过调节功能控制基因表达χRb公司χ(方程式(3a)和(4a)in图2),由通量平衡条件设置(公式(2a)和(2b)图2). 组合中显示的所有方程式图2(详见补充注释4.1),我们得到了描述位移后动力学的单一非线性微分方程(> 0):

d日σd日=σ()[μ(f)χ^(σ())σ()χ^Rb公司(σ())],
(5)

哪里χ^χ^Rb公司监管职能是否在图2b.方程(5)可以用解析法求解(补充说明4.2). 这种逻辑式方程的时间刻度由速率设定μ(f),表示变换后介质中的最大碳吸收速率;它的值可以通过与稳态增长率的关系来确定λ(f)在换班后介质中,μ(f)λ(f)/(1 −λ(f)/λC类),其中λC类=1.17小时−1是生长定律中已知的应变特定常数(扩展数据图1d补充说明2.3). 群众的时间进程M(M)(),核糖体M(M)Rb公司(),以及引入的所有其他数量图2可以从FCR模型中分析得出(补充注释4.4和方法). 对于许多宏观观察,第二个时间刻度由λ(f)方程式(5)由于生长介导的遗传蛋白质组的缓慢稀释,已经松弛。

解决方案还需要指定初始条件σ(0),转换后直接进行的翻译活动。对于升档,这取决于升档时摄取系统Cat2(用于新添加的底物)的丰度,因为预先表达的Cat2导致碳流入的突然增加导致σ(0). 许多吸收系统在缺乏同源底物的情况下不表达,因此生物量通量J型()和σ()是连续的,并且方程式(5)能够在没有任何可调参数的情况下描述整个升档动力学(图1a-e扩展数据图3,黑色面板)。

在某些情况下,我们观察到生物合成通量的突然跳跃J型(t) 升档瞬间(例如,请参见,扩展数据图6)可能是由于Cat2在移前生长期间的表达。这种预表达式可以通过改变初始条件来适应FCR模型σ(0) (扩展数据图6,黑线)。通过在移前生长期间综合预表达Cat2,流量突变的幅度J型()在转变的时刻可以系统地改变,即使瞬时增长率暂时超过转变后状态的增长率(扩展数据图7a、b). 对于这些情况,FCR模型定量地预测了通量跳跃的幅度,如扩展数据图7d.

相同的FCR模型定量地捕获了共同利用的碳基质之间非常不同的降档动态(图1f–j扩展数据图4,黑线)。在这里,使用具有既定Michaelis常数的Monod动力学来模拟碳底物消耗(补充说明4.2a),同样没有可调参数。底物耗尽后通量和生长速率下降的幅度由Cat2的表达水平决定,与仅底物S2相比,Cat2在两种底物上的生长过程中下调(补充注释4.2d和补充方程(57)). FCR模型进一步扩展,以捕捉分级使用的碳基质之间的生长变化,特别是聚合葡萄糖-乳糖二糖1(图3). 这里,诱导剂排除26防止乳糖的摄入和紫胶只要葡萄糖存在并被吸收,操纵子就会被激活。在不调用如何精确删除诱导子排除的详细信息的情况下,我们引入了紫胶操纵子作为孤立拟合参数(补充说明3.2g)并发现FCR模型也准确地描述了葡萄糖-乳糖二糖(图3,黑线)。

保存图片、插图等的外部文件。对象名为nihms956215f13.jpg
分级使用的碳基质之间的双氧水迁移

,大肠杆菌K-12生长在0.03%葡萄糖和0.2%乳糖上。最初,只使用葡萄糖1.大约时间=0的葡萄糖被耗尽,乳糖吸收和代谢被激活。红色虚线表示仅在乳糖上的生长速率,λ(f)=0.95小时−1.b条、瞬时增长率λ()葡萄糖耗尽后突然下降,随后迅速恢复到λ(f)(水平虚线)。c(c),葡萄糖消耗后,生物量通量以增加的速率增加,然后达到最终速率。红色虚线的斜率表示λ(f),橙色虚线的斜率表示μ(f) d日,乳糖分解代谢酶LacZ(绿色三角形)在转换前受到严格抑制(98 U ml−1OD处600纳米),然后在转移后迅速增加,然后稳定在绿色虚线斜率(1557 U ml)指示的最终表达水平−1OD处600纳米,仅从乳糖的稳态生长中获得)。带斜率的橙色虚线μ(f)/λ(f)×1557单位ml−1外径600纳米=7866单位毫升−1外径600纳米,是根据补充注释2的补充方程(18)黑色曲线是使用以下时间点对FCR理论的完整预测紫胶作为单个拟合参数的操纵子激活(参见补充说明3.2g).

考虑到模型的简单性,特别是缺乏自由参数,模型输出与实测观测值的定量一致性是显著的。如中所述扩展数据图8,升档和降档显示出不同的恢复动力学,尽管受相同的基本监管策略控制。在升档时,蛋白质合成的重新分配很快就会解决(扩展数据图8a,虚线),分解代谢蛋白和核糖体丰度向其最终状态的生长介导的松弛是主要的恢复过程。在降档时,蛋白质合成的主动重新分配阶段跨越整个恢复(扩展数据图8b,虚线)。由于适应能力不同,升档主要受时间尺度的影响λ(f)−1,而降档主要是μ(f)−1(补充说明6扩展数据图9),解释了为什么图3(μ(f)=5.1小时−1)比升档快五倍图1(λ(f)=0.98小时−1). 在这两种转变中,尽管调节函数是从它们派生出来的,但它们并没有沿着稳态线滑动(补充说明2; 比较中的黑色线条和彩色线条扩展数据图8c).

考虑到克服代谢瓶颈所需的所需蛋白质数量很少,营养素转换后的完全恢复速度惊人地慢;例如,用于升档的Cat2和用于降档的Cat1不超过整个蛋白质组的百分之几。因此,如果按照最新的最佳生长转换策略理论,细胞只将蛋白质合成导向有限数量的瓶颈蛋白质,那么它们只需要加倍时间的百分之几就可以合成27,28相反,根据FCR模型,细胞采用了一种调节策略,其特点是蛋白质组的协调全球重塑远远超出了对少数特定蛋白质的调节,但包括例如大量与培养基中特定糖无关的合成代谢和分解代谢蛋白质。事实上,FCR模型对蛋白质组重塑动力学进行了定量预测:根据碳限制下的稳态反应对蛋白质组进行分类19,我们将在碳限制下增加、减少或不改变的蛋白质分为C↑、C↓和C“部门”。因为分解代谢蛋白属于C↑,核糖体属于C↓19这些蛋白质组部门的动态,由调节功能分配χC↑χC↓.,由监管职能部门在FCR模型中规定χ()和χRb公司()根据中的关系式(6a)和(6b)扩展数据图8d,所有常数均由稳态蛋白质组分配固定(补充说明5). 扇区C↑和C↓的动态,表现为这些扇区中所有蛋白质的总质量分数的变化,表示为φC↑()和φC↓(),然后由方程式(7a)和(7b)规定,如所示扩展数据图8d升档和降档产生的扇区动态定性地描述在扩展数据图8e和f分别是。

我们通过扩展定量质谱方法来测试这些预测19到动力学状态。图3a显示了数百种蛋白质在向上和向下移动过程中丰度的时间依赖性变化图1我们观察到,蛋白质的动态反应符合其稳态分类:碳降档后,C↑蛋白质增加(绿色),C↓蛋白质减少(红色),而升档时则相反。C的大多数蛋白质对养分转移几乎没有反应。蛋白质丰度的变化主要取决于蛋白质表达和生长介导的稀释,因为在过渡过程中降解作用微不足道(扩展数据图10). 我们估计了属于这三个部分中每一部分的蛋白质的总丰度(见方法),并绘制了它们随时间变化的曲线(图4b、c). 在整个过渡过程中,它们的动力学与模型预测(线)密切相关。同样的方法也被用于表征属于各种主要代谢途径或官能团的蛋白质的总丰度(图4d和方法)。这些组在升档和降档时的动力学也基本上遵循了模型预测,但运动蛋白除外,运动蛋白在营养降档时明显表现出延迟。

保存图片、插图等的外部文件。对象名为nihms956215f14.jpg
恢复期间的蛋白质组组成

,单个蛋白质在图1,相对于预换档稳定状态(色标:对数2). 蛋白质根据其对碳限制降低生长速率的稳态反应(见方法)进行分组:C↑(增加;绿色)、C↓(减少;蓝色)和C(无响应;黑色)19.b条,c(c),蛋白质部分的绝对蛋白质组分数φC↑,φC↓φC−(符号)和非参数预测(线),假设C↑、C↓和C扇区中的所有蛋白质遵循简单的全球分配规则(扩展数据图9补充说明5).d日,升档(实线和符号)和降档(虚线和开放符号)期间单个路径或蛋白质组的质量分数。颜色表示部门隶属关系。所有生物合成途径都遵循假设共同调控的理论预测。运动性表达式(右下角)是一个明显的例外,与表达式的预期增加值(虚线)相比,从降档恢复过程中的延迟增加了2小时。TCA,三羧酸。

模型预测与不带可调参数的蛋白质组数据的定量匹配表明,许多蛋白质组确实是共同调节的,由模型中假设的单一动态控制变量在整个生长转变过程中控制。在我们的模型中,这个动态变量是平移活动σ(),它本身是中央代谢物池的代理(例如酮酸和氨基酸),如所示图2a因为在分子上,控制是通过信使cAMP来实现的12,15和ppGpp13,14,16,17,我们的结果表明,这两条信号通路通过中央代谢物池紧密耦合15,22这是全球监管控制的拟议驱动力。值得注意的是,即使是已知不受cAMP或ppGpp直接调控的途径,例如氨基酸合成,也会遵循碳转移后预测的全球反应,这可能是由于对转录或翻译资源的竞争16,17.

尽管这些信号具有中央调节作用,但生长过渡动力学对这些信号通路的动力学细节并不敏感,这一点可以从FCR模型的预测能力中得到证明,FCR模型并未明确包含信号。相反,动力学是由蛋白质合成和生长介导的稀释的全球重新分配的较慢时间尺度控制的。观察到的资源分配策略保持蛋白质组的全局协调,而不是按顺序优先处理各种瓶颈的表达27,28通过将代谢瓶颈限制在中心前体,这种保守的控制策略可能更加稳健,而中心前体推动全球监管控制(图2a). 定量了解蛋白质组重建动力学是预测细胞动力学的先决条件,例如遗传回路的动力学或应力生长的动力学,使此处建立的框架成为理解其他动力学现象的广泛光谱的基础。

方法

没有使用统计方法预先确定样本量。这些实验不是随机化的,研究人员在实验和结果评估期间也没有对分配盲视。

生长培养基

本研究中使用的所有生长培养基均基于N-C最小培养基29,包含K2SO公司4(1克),K2高性能操作4·3小时2O(17.7克),千赫2人事军官4(4.7克),硫酸镁4·7小时2O(0.1克)和NaCl(2.5克)/升。用20mM NH补充培养基4氯作为氮源,各种碳源。必要时向培养基中添加1 mM IPTG,以完全诱导本机紫胶加入操纵子或0.1 mM IPTG诱导P拉卡-O1启动子驱动XylR用于可滴定表达dctA公司在应变NQ537中。的表达式聚氨基甲酸酯在指定浓度下添加3-甲基苄基醇(3MBA)诱导后,启动子被调节剂XylR激活。扩展数据图1h,i包括所有使用的碳源介质支持的稳态增长率。

培养程序

每个实验分三个步骤进行:种子培养、预培养和实验培养。实验前,将从−80°C甘油原料中提取的细胞接种在LB琼脂平板上。选取一个菌落,在37°C和250 r.p.m的摇晃水浴中的新鲜LB上生长。选择与实验培养物相同的预培养基。选择用洗过的种子培养物接种预培养物,使隔夜培养的预培养物在实验早上仍呈指数增长。细胞保存至少10倍(生长速度低于0.4小时时至少5倍−1)在37°C下的预培养中,在250 r.p.m下摇晃。选择早上接种到预加热的实验培养物中,以便在测量生长之前,细胞在37°C下至少再重复3次,并在250 r.p.m下摇动。

对于较小的培养体积,在20 mm×150 mm的玻璃试管中培养5 ml实验培养物。对于较大的容量,25 ml(100 ml)生长在125 ml(500 ml)带挡板的Erlenmeyer烧瓶中。在每个时间点,提取200μl样品,并在ThermoScific Genesys 20分光光度计中使用1 mm路径长度的Starna Sub-Micro试管在600 nm处测量光密度。

高精度光密度数据如所示图1使用流动细胞装置获得。将25 ml体积的实验培养物在37°C和250 r.p.m.的温度下在空气培养箱中的125 ml带挡板的Erlenmayer烧瓶中摇晃。使用蠕动泵使培养基循环通过流动细胞比色皿。外径600纳米使用LED光源(592 nm)和由微控制器控制的光电探测器进行测量。所有仪器和管道在空气培养箱内保持在37°C。

菌株

本研究中使用的所有菌株均来自野生型大肠杆菌K-12菌株NCM372230和总结如下扩展数据图1j.可滴定dctA公司如前所述构建表达菌株15简而言之,包含公里基因和聚氨基甲酸酯从质粒pKDPu中PCR扩增启动子并整合到大肠杆菌上游应变NQ351dctA公司(相对于的平移起点−1到−182 bpdctA公司),通过使用λ红色系统31.结果km-Pu-dctA将上述菌株中的P(P)拉卡-O1-木马R通过P1转导,产生菌株NQ530。构建可滴定的lacZ公司表达菌株,km-Pu-lacZ公司在里面补充方程(40)通过P1转导将其转入菌株NQ386,得到菌株NQ537。应变NQ381和NQ399已在前面进行了描述16.A型dctA公司敲除子通过P1转导途径从JW3496-1转移到NCM3722,产生菌株NQ1324。

β-半乳糖苷酶定量

β-半乳糖苷酶的定量基于Miller测定法32并如前所述使用15.

总RNA定量

RNA定量方法基于上述方法33按说明进行修改15.

总蛋白质定量

使用商业微BCA测定法对总蛋白进行定量。总之,通过离心收集1 ml样品,并将其储存在−80°C下,直至分析。为了进行分析,再次清洗细胞,将其重新悬浮在N-C培养基中,并以1:10稀释。接下来,在60°C的水浴中培养200μl稀释样品和200μl microBCA工作试剂。1小时后,将样品放在冰上,通过离心去除细胞碎片,并记录592 nm处的吸光度。每个分析都使用牛白蛋白血清标准物进行单独校准。

干质量量化

将150ml细胞培养物在带挡板的250ml烧瓶中生长至OD600纳米在37°C的水浴振动器中约为0.5。随后收集50 ml细胞培养物的三份样品,在冰水浴中摇晃冷却,并通过离心浓缩。每一步都记录上清液,并从OD中减去上清液600纳米正在读取。将所得颗粒转移到1.5 ml试管中,并在80°C下干燥过夜。用高精度天平测定颗粒的最终质量。

生长过渡质谱协议

我们使用定量质谱分析了升档和降档的动力学序列,如图1在大约两倍的时间内,我们为每个动力学转变收集了两个班前和六个对数间隔的班后样本,每个样本都包含OD下1ml液体培养物的等效生物量600纳米=1,17000离心收集持续2分钟。为了使用质谱法实现轻/重相对定量,我们准备了一个15N稳定同位素标记参比物,包含稳态葡萄糖和琥珀酸-动物培养基培养物的1:1(w/w)混合物15全日空航空公司4氯。

在纯水中重新悬浮后,100μl混合液15含1 ml OD液体培养物当量的N参考600纳米=将1个细胞生物量添加到之前储存在−20°C下的每个样品细胞颗粒中。随后的样品处理,包括TCA沉淀、半胱氨酸还原、烷基化、胰蛋白酶消化和脱盐程序,如前所述进行24.

在与Sciex 5600 TripleTOF系统耦合的Eksigent NanoLC Ultra系统上分离色氨酸肽,并在数据依赖的鸟枪模式下进行分析(250 ms MS1扫描,然后40个150 ms MS2扫描周期)24每个样品注射两次,得到两个技术复制品。

如前所述24,供应商WIFF格式的仪器数据文件使用Sciex软件转换为轮廓和质心MZML格式。使用跨蛋白质管道34,以中心为中心的mzML文件被转换为mzXML并使用X!串联35违背习俗大肠杆菌数据库(来源于UniProt生物83333)。使用SpectraST将MS2搜索结果合并到原始光谱库和一致光谱库中。

MS1强度封套包含14N(轻)-和15使用最小二乘傅里叶变换卷积算法拟合N(重)标记离子19使用isodist实现36在Perl程序屠杀中实现37.相对14N个/15N比率被校正为不平衡15N参考移液管通过累积轻(样品)和重(参考)X的商!串联光谱计数大肠杆菌以每个样本为基础的肽。绝对蛋白质表达(即特定蛋白质的总蛋白质质量分数)通过X!串联肽-光谱匹配。

对于每次升档和降档,“绝对表达式数据”的时间序列可以表达式矩阵的形式表示E类m、 n个,带有时间序列(n个=1…8),以及蛋白质的数量,,作为行。

以前分类为在碳限制下增加、减少或无反应的蛋白质19分为C↑、C↓和C部门。24个以前未检测到的蛋白质19根据其线性响应的最近距离(升档和降档的平均值)和属于每个扇区的蛋白质的平均线性响应进行分类。

数据如所示图4a使用相对表达数据,这些数据是通过将“绝对表达数据”标准化为预移位表达数据(一个数据点的平均值=−15分钟,以及换档前的第二个数据点)。如果两个换档前条目都为空,则无法将数据标准化并从分析中删除。此外,为了专注于高质量的表达数据,删除了包含四个以上空条目的行(蛋白质),总共得到=647个蛋白质,提供于补充表1.数据为日志2在中转换并绘制为热图图4a.

数据如所示图4b、c使用“绝对表达数据”(即蛋白质的绝对质量分数)。单蛋白数据见补充表2.一个部门的丰富性X(X)在时间输入n个,φ十、 n个,通过蛋白质表达求和计算第页m、 n个所有条目在蛋白质的条件下是该行业的一个元素X(X),提供于补充表3蛋白质组的计算类似于蛋白质部分;所有蛋白质组的质量分数约为蛋白质组的60%。数据见补充表4蛋白质分类(部门和组)见补充表1.

脉冲标记质谱协议

细胞在含10 mM的N-C培养基中生长14全日空航空公司4氯是氮的唯一来源。在葡萄糖耗尽(降档)和添加葡萄糖酸盐(升档)时,细胞15通过添加脉冲来同位素标记N15全日空航空公司4Cl至10 mM。脉冲允许预移位蛋白质的分化(14N) 来自移位后产生的蛋白质(50%15N) ●●●●。我们收集了两个班前和六个班后样本(1个OD600纳米×ml总生物量)。A类15将N同位素富集参比物制备为稳态琥珀酸盐和葡萄糖最低限度培养基的1:1(w/w)混合物,在15全日空航空公司4Cl并如上所述添加到每个样品中。

数据采集和处理按照“生长转换质谱协议”中的描述进行。

MS1强度封套包含14N标记(预移位)离子,50%15N(位移后)离子,100%15使用isodist定量N-标记(恒定参考)离子36以符合三物种模型,包括50%15N-标记的信封。

FCR模型在特定增长转变中的应用

FCR模型基于稳态增长规律(补充方程式(9)和(12)),如所示扩展数据图1c、d并在中进行了讨论补充说明2)作为模型的一般输入,与具体的增长变化无关。线性增长定律中出现的常数(φ卢比,0,γ,λC类)修正调节功能中出现的参数χ^χ^Rb公司(补充方程式(28)和(29)),进入翻译活动的中心方程σ()英寸方程式(5)(或补充方程(43)). 此外,FCR模型需要特定于生长转移的参数:初始生长速率λ和最终增长率λ(f).最终增长率λ(f)修复参数μ(f)(出现在方程式(5)),它通过生长规律描述了换班后营养素的质量(补充方程(17)).

翻译活动的初始条件σ(0),定义见补充方程(50),由位移类型决定,并可能导致明显不同的过渡动力学。对于在共同利用的碳源之间没有预先表达的升档,初始条件σ(0)等于预设值σ由班前增长率决定,λ看见补充方程(52)对于带有预表达式的升档,初始条件σ(0)需要预表达式级别的拟合,请参见补充方程(54)对于共同使用的碳基底之间的降档,初始条件σ(0)由剩余基质的碳流入量设定,可根据初始和最终生长速率计算,λ、和λ(f)和增长规律,请参见补充方程式(56)对于分级使用的碳基板之间的降档,初始条件σ(0)需要拟合一个参数。

所有变量的移位前增长由初始增长率决定λ以及增长规律,请参见补充方程式(58)–(64).位移后增长是根据σ()(解析为补充方程(49)); 该解用于导出所有其他动态量的一般分析解,包括核糖体丰度(M(M)Rb公司(),补充方程(66))分解代谢蛋白(M(M)(),补充方程(73)),生物质合成(J型R(右)(),补充方程(68)),生物量(M(M)(),补充方程(76))、和增长率(λ(),补充方程式(77)); 看见补充说明4.4用于推导这些结果。如果碳吸收受到基质限制(补充方程式(19)和(20))当底物浓度降低到非常低的水平时,该模型只能通过数值求解(补充说明4).

数据可用性

支持本研究结果的数据可根据合理要求从通讯作者处获得。

扩展数据

扩展数据图1

保存图片、插图等的外部文件。对象名为nihms956215f1.jpg
增长规律

,每光密度OD的干质量600纳米与所研究条件下的生长速率无关,平均干重为509μg/ml×OD600纳米值(黑线)。b条,每光密度OD的蛋白质质量600纳米根据细胞中RNA和DNA在较高生长速率下的增加,显示出对生长速率的轻微依赖性(黑线:眼睛指南)38。由于所示升档和降档的增长率范围内的变化很小,因此我们认为在所示的整个换档过程中,从干质量到总蛋白质的转换是恒定的。干质量和总蛋白质数据取自参考文献。39.c(c),细胞核糖体蛋白质组部分ϕRb公司(λ)稳态增长率λ*使用RNA或蛋白质作为报告者。数据拟合成线性关系ϕRb公司=ϕRb公司,0+λ/γ具有φ卢比,0=(0.049±0.02)和γ=(11.02±0.44)小时−1.d日,分解代谢蛋白(ϕ(λ),左边)以LacZ(用1 mM IPTG诱导)作为报告物进行测量,并测量米勒单位(MU)(右轴;参见参考文献。32对于MU)。线性关系拟合为ϕLacZ公司=ϕLacZ公司,最大值(1λ/λC类)具有ϕLacZ公司,最大值=(21.8±0.5)×10λC类=(1.17±0.05)小时−1.e(电子),翻译活动σ*计算依据为ϕRb公司(λ)(c(c))作为σ=λ/ϕRb公司(λ).f、 克、监管职能χ^Rb公司(σ)χ^(σ),翻译活动的功能σ*,通过替换计算λ*由σ*从配合c(c)d日(请参见补充方程式(28)和(29)).h、 我,表格显示绘制的所有数据c(c)d日中的数据小时通过不同碳基质上的生长测量通过具有可滴定转运体表达的菌株的生长进行测量,参见j个.

扩展数据图2

保存图片、插图等的外部文件。对象名为nihms956215f2.jpg
长期动态和重复性

,升档的系列稀释试验图1(琥珀酸,添加葡萄糖酸)。上图显示了光密度OD600纳米作为生物量的测量M(M)(),中间的图表显示了生物量通量J型(),下图显示了增长率λ() =J型()/M(M)(). 升档后不久,细胞培养物在新鲜培养基中稀释4倍和16倍,与原始培养物平行生长,当稀释培养物的光密度达到0.05左右时,记录光密度。由于生长速率与细胞密度无关,原始培养物和稀释培养物的数据会崩溃,表明生长速率在较长时间内达到其最终水平(虚线)。b条,可重复性。生物质(顶部)、通量(中部)和增长率(底部)从左到右图1a–e以及中所示的降档图1f–j理论(黑线)和数据(红钻石)与图1四个独立的重复序列(紫色、橙色、绿色和蓝色)绘制在中所示数据的顶部图1.

扩展数据图3

保存图片、插图等的外部文件。对象名称为nihms956215f3.jpg
从琥珀酸和丙酮酸上升

a–h,NCM3722在20 mM琥珀酸盐(成功)或20 mM丙酮酸盐(丙酮酸)作为唯一碳基质上呈指数增长。=0秒,随后添加共同利用的碳基质:0.2%阿拉伯糖(,e(电子)); 0.2%木糖(b条,(f)); 0.2%甘油(c(c),); 0.2%葡萄糖(d日); 20 mM葡萄糖酸盐(小时). 左侧面板显示光学密度OD600纳米(生物量的测量M(M)(),红色圆圈),中间面板显示OD的导数600纳米(生物量通量的测量J型(),红色方块),右侧面板显示增长率(λ() =J型()/M(M)(),红色钻石)。理论路线M(M)(),J型()和λ()计算方法为补充方程式(76)、(70)和(77),使用初始条件补充方程(52)用于不带预表达的升档。在各碳基质上的稳态生长期间,测量了初始和最终生长速率。

扩展数据图4

保存图片、插图等的外部文件。对象名为nihms956215f4.jpg
使用共同利用的碳基板降档

a–d,NCM3722在20 mM琥珀酸盐上呈指数增长(成功)(a、 b条)或20 mM丙酮酸(丙酮酸)(c、 d日),与0.56mM葡萄糖酸盐(葡萄糖)组合(a、 c(c))或1.11 mM甘油(甘油)(b、 d日). 在附近=0,葡萄糖酸或甘油被耗尽。左侧面板显示光学密度OD600纳米(生物量的测量M(M)(t) ,红色圆圈),中间面板显示OD的导数600纳米(生物量通量的测量J型(),红色方块),右侧面板显示增长率(λ() =J型()/M(M)(),红色钻石)。理论路线M(M)(),J型()和λ()通过求解平移活动的微分方程进行计算σ() (补充方程(22))以及定义内部通量和蛋白质含量的(补充方程(19)-(21)和(22)-(24))数字上。消耗底物(甘油或葡萄糖酸盐)的摄取量计算如下补充方程(20),使用Michaelis常数K(K)=212μM(葡萄糖酸盐40)和K(K)=5.6μM(甘油41). 由于葡萄糖酸转运亲和力低,它消耗缓慢(参见c(c))我们的分解代谢模型准确地描述了这一点。在各碳基质上的稳态生长期间,测量了初始和最终生长速率。

扩展数据图5

保存图片、插图等的外部文件。对象名为nihms956215f5.jpg
无基因调控生长动力学模型的性能

a–c,该图显示了正文中讨论的三种营养素变化的结果:琥珀酸葡萄糖从图1b(); 琥珀酸-葡萄糖从图1g(b条); 葡萄糖-乳糖二糖从图3b(c(c)). 实验数据显示为红色符号。在现有的增长转型模型中,Dennis和Bremer提出的模型可以预测升档的增长转型42该模型假设通量瞬时跳至最终状态,无需调节方案(参见补充说明1查看)。该模型的预测,这里称为“瞬时模型”,在所有面板中显示为绿色虚线。尽管瞬时模型适用于升档(面板),它是为其开发的42,它无法描述降档(面板b条c(c)),因为它不能描述通量减少。瞬时模型无法通过匹配葡萄糖消耗后的流量来固定(面板中的蓝色虚线b条c(c))由于缺乏监管计划,它低估了回收率(补充说明1). 相比之下,本研究中引入的FCR模型基于生长过渡期间蛋白质合成的主动重新分配,一致预测了升档和降档的正确回收率(黑色实线)。

扩展数据图6

保存图片、插图等的外部文件。对象名称为nihms956215f6.jpg
甘露糖升档

a–c,NCM3722在0.1%甘露糖作为唯一碳基质上呈指数增长。添加=0(虚线)草酰乙酸(OAA),随后共同使用。光学密度OD600纳米(生物量的测量M(M)(),红色圆圈)(),OD的导数600纳米(生物量通量的测量J型(),红色方块)(b条),增长率(λ() =J型()/M(M)(),红色钻石)(c(c)).d日,分解代谢蛋白报告者LacZ的表达在移位后暂时受到抑制(垂直虚线)。黑色虚线表示使用无预表达式的初始条件的理论解(补充方程(52)),这与实验数据(符号)不一致。使用初始条件进行预表达式(补充方程(54))该理论(黑线)拟合了预表达水平,很好地描述了数据,包括LacZ在d。LacZ活动单位U定义为MU×OD600纳米.

扩展数据图7

保存图片、插图等的外部文件。对象名称为nihms956215f7.jpg
预表达对升档动力学的影响

,大肠杆菌NQ530在含有不同水平琥珀酸转运体(DctA)的甘露糖最小培养基中呈指数增长,通过可滴定法预先表达聚氨基甲酸酯启动子(调节子:XylR,诱导子:3MBA)。时间=0,培养基补充20 mM琥珀酸盐,DctA表达设置为所有升档的共同水平。b条、生物质M(M)()(外径600纳米),生物质通量J型()(OD的导数600纳米)和增长率λ()针对不同水平的预表达(面板中显示的3MBA浓度)。DctA预表达对转移后生长有积极影响。生物量通量J型()(中间)显示跳跃式增长,然后以最终速度指数增长λ(f)对于预表达水平较高的培养物,瞬时生长率λ()(右)在放松到最终值之前短暂超调λ(f)这种增加是暂时的,因为合成结构中的DctA水平最终会降低到预表达水平以下。当使用初始条件时,动力学理论定量地捕获了所有预表达水平的升档动力学σ(0)作为单个拟合参数(实线)。(c(c))NQ 1324(ΔdctA公司)尽管琥珀酸转运体DctA被淘汰,但升档后增长率略有增加。单用甘露糖进行稳态生长时,最终生长速率没有显著增加。NQ1324仅琥珀酸盐的稳态增长率几乎无法检测到(λ≲0.02/h)。e(电子)、理论预测的验证。生长速率Δ的跃增幅度λ(如右面板所示b、 c(c))线性依赖于预表达水平。跳跃式增长dctA公司中显示的淘汰c(c)用三角形表示。

扩展数据图8

保存图片、插图等的外部文件。对象名称为nihms956215f8.jpg
蛋白质组重塑

a、 b条、监管职能χ和蛋白质组分φ分解代谢蛋白质(相对于最大表达)、核糖体(每毫克蛋白质的毫克RNA吸光度)以及生长速度λ()和翻译活动σ()在增长转变期间图1.分解代谢蛋白的调节功能增强后不久χ()和核糖体χRb公司(),由翻译活动驱动σ()(绿色、蓝色和红色虚线),已放松至最终状态(). 因此,一个新的蛋白质组以最终比例合成,通过生长介导的稀释导致缓慢收敛到最终状态。葡萄糖消耗后碳流入不足导致χ()以及下调χRb公司()由翻译活动驱动σ()(绿色、蓝色和红色虚线)(b条).c(c)、监管职能χ()和χRb公司(),由翻译活动设置σ(),绘制与增长率λ(). 细长的黑线表示稳态增长规律(扩展数据图1c、d补充说明2). 尽管调控功能来源于生长规律(黑线),但它们在生长过渡期间仍有很大差异,因为它们受翻译活动的控制σ()而不是增长率λ().d日假设蛋白质组扇区C↑和C↓与分解代谢蛋白质和核糖体蛋白质共同调节(方程式(6a)和(6b)),则蛋白质组扇形的动力学可以描述(方程式(7a)和方程式(7b);看见补充说明5详细信息)。e(电子),(f),蛋白质组重构的图形概要b条,如方程式(6a)、(6b)、(7a)和(7b)所预测。红色边界和箭头表示动态。在升档调节功能中χ迅速放松到最终状态;在降档时,调节功能最初超调,然后缓慢放松到其最终状态(参见补充说明6供扩展讨论).

扩展数据图9

保存图片、插图等的外部文件。对象名称为nihms956215f9.jpg
双相弛豫动力学

该图显示了从FCR模型获得的生长位移的弛豫动力学。-c(c),增长缓慢λ=0.3小时−1(橙色)和从快速的换档前增长降档λ=0.95小时−1(绿色),这两个值在转换后增长率相同,λ(0)=0.3小时−1,以及同样快速的最终增长率λ(f)=0.95小时−1.d–f日,增长缓慢λ=0.2小时−1(橙色)和从快速的换档前增长降档λ=0.95小时−1(绿色),带有λ(0)=0.2小时−1和相同的缓慢增长率λ(f)=0.95小时−1.,d日,增长率空间中的恢复轨迹(左-轴),与分解代谢蛋白质部分成比例(右-轴)和核糖体部分φRb公司()(底部x-轴)。初始条件φRb公司(0)取决于班前增长率λ(顶部x-轴)。轨迹从这些点开始(φRb公司(0),λ(0),用彩色方块表示,并在最终状态结束(φRb、f,λ(f))用黑色十字标记。详见补充说明6,动力学是双相的,可以通过简单的几何结构来理解,从而产生虚线和虚线(参见补充注释6.2). 首先,与时间尺度相关的放松μ(f)−1(补充方程(116))沿着灰色虚线,在此期间,蛋白质合成比蛋白质组的最终组成含有更多分解代谢蛋白质和更少核糖体。该动力学阶段的结束由所有面板中的圆圈表示。第二,沿时间刻度上的对角线虚线缓慢移动λ(f)−1(补充方程(120)). 沿着这条路线,翻译活动,σ() =λ()/φRb公司(),已放松到最终状态(请参阅b条,e(电子)),其中σ()在“圆圈”之后放松了。因此,调节功能也放松到了最终状态,蛋白质合成包含的分解代谢蛋白质和核糖体数量与最终状态相同。在第二阶段,由于生长介导的遗传蛋白质稀释,蛋白质组逐渐适应。c(c),(f),增长率在第一阶段几乎完全适应(从=0到圆圈)用于降档(绿色、青色),但不用于升档(橙色、紫色)。这种差异是因为降档中的核糖体丰度较高,这使得分解代谢蛋白的表达增加,因此几乎完全避免了缓慢阶段,并且在向快速最终生长转变时更为突出(a–c),而不是减缓最终增长(d–f日).

扩展数据图10

保存图片、插图等的外部文件。对象名为nihms956215f10.jpg
蛋白质降解

,脉冲标记允许区分脉冲前和脉冲后蛋白质质量。我们添加了浓缩15全日空航空公司4在葡萄糖酸盐升档时或降档期间葡萄糖耗尽前几分钟,Cl进入培养基。与第三种同位素相比(峰值15N-reference culture)允许跟踪一段时间内的预移位蛋白水平。b条,示意图显示了当培养物瞬时上移时,脉冲时间内存在的总细胞蛋白(黑色)和细胞蛋白(蓝色)的水平。稳定的蛋白质水平以零斜线(蓝色实线)为特征,而降解或输出的细胞蛋白质呈现负斜线(蓝虚线、红色箭头)。c(c)采用脉冲标记法对40种质量分数最高的细胞蛋白的位移后水平进行了量化。灯(L,14N) 超重(H,15N) 相对蛋白质水平绘制为时间的函数。这些蛋白质具有不同的生物学功能、细胞定位、大小和结构。总的来说,它们占整个转换过程中检测到的总蛋白质质量的35-40%(通过汇总表中列出的它们的质量丰度进行估算补充表2对于每个条件)。除了鞭毛蛋白(fliC,红框)外,我们没有观察到大量蛋白质在2小时时间尺度上的升档或降档序列的蛋白质水平下降,我们可以自信地量化。Flagellin通过一个专用的运输系统输出到细胞外围,并可能在稳态生长期间释放到培养基中。由于脱落蛋白的收集效率与细胞中的蛋白质不同,因此预先标记的FliC的减少可作为该方法的阳性对照。

补充材料

补充信息

单击此处查看。(120万,pdf)

致谢

我们感谢E.Mateescu发起这项工作,并感谢J.Hasty和R.Young进行讨论。这项工作得到了国立卫生研究院(NIH;拨款1R01GM109069)和西蒙斯基金会(拨款330378)通过T.H.的支持,NIH(拨款1R01-GM118850)通过J.R.W.的支持,以及德国研究基金会通过卓越集群“慕尼黑纳米系统倡议”和优先项目SPP1617(拨款GE1098/6-2)的支持通过U.G。

脚注

在线内容方法以及任何其他扩展数据显示项和源数据都可以在论文的在线版本中获得;这些章节独有的参考文献仅出现在在线论文中。

补充信息可在论文的在线版本中找到。

作者贡献

D.W.E.、S.J.S.、U.G.和TH设计了这项研究。生长位移由D.W.E.和S.J.S.执行,并由S.J.S、D.W.E和T.H.进行分析。定量质谱分析由V.P.执行,由S.JS.、VP和J.R.W.D.W.E..进行分析,S.J.S..、U.G.和T.H开发了模型,所有作者都参与了论文的撰写补充信息.

作者信息

作者声明没有竞争性的经济利益。欢迎读者对该论文的在线版本发表评论。

出版商注:Springer Nature在出版的地图和机构关联中的管辖权主张方面保持中立。

工具书类

1Monod J.酶适应现象及其与遗传学和细胞分化问题的关系。增长。1947;11:223–289. [谷歌学者]
2Kjeldgaard NO,Maalöe O,Schaechter M.平衡生长期间不同生理状态之间的转换鼠伤寒沙门菌.微生物杂志。1958;19:607–616.[公共医学][谷歌学者]
三。Dennis PP,Bremer H.《细胞内核糖核酸合成的调控》大肠杆菌:对转移的分析。分子生物学杂志。1973;75:145–159.[公共医学][谷歌学者]
4Thattai M,Shraiman BI大肠杆菌.生物物理学杂志。2003;85:744–754. [PMC免费文章][公共医学][谷歌学者]
5Kussell E,Leibler S.波动环境中的表型多样性、种群增长和信息。科学。2005;309:2075–2078.[公共医学][谷歌学者]
6Acar M,Mettetal JT,van Oudenaarden A.波动环境中作为生存策略的随机切换。自然遗传学。2008;40:471–475.[公共医学][谷歌学者]
7Boulineau S等人。单细胞动力学揭示了双相转变期间的持续增长。公共科学图书馆一号。2013;8:e61686。 [PMC免费文章][公共医学][谷歌学者]
8.Kotte O、Volkmer B、Radzikowski JL、Heinemann M大肠杆菌中枢碳代谢。分子系统生物学。2014;10:736. [PMC免费文章][公共医学][谷歌学者]
9van Heerden JH等人。迷失在过渡中:糖酵解的启动产生非生长细胞的亚群。科学。2014;343:1245114.[公共医学][谷歌学者]
10Venturelli OS、Zuleta I、Murray RM、EI-Samad H。酵母代谢项目中的人口多样化促进了对环境变化的预期。《公共科学图书馆·生物》。2015;13:e1002042。 [PMC免费文章][公共医学][谷歌学者]
11Link H,Fuhrer T,Gerosa L,Zamboni N,Sauer U。饥饿和生长之间代谢转换的实时代谢组分析。自然方法。2015;12:1091–1097.[公共医学][谷歌学者]
12Kolb A、Busby S、Buc H、Garges S、Adhya S。cAMP及其受体蛋白的转录调控。生物化学年度收益。1993;62:749–795.[公共医学][谷歌学者]
13Potrykus K,Cashel M.(p)ppGpp:仍然神奇吗?微生物年鉴。2008;62:35–51.[公共医学][谷歌学者]
14Lemke JJ等人。直接监管大肠杆菌转录因子ppGpp和DksA的核糖体蛋白启动子。美国国家科学院程序。2011;108:5712–5717. [PMC免费文章][公共医学][谷歌学者]
15You C等。通过循环AMP信号调节细菌蛋白质组与代谢的协调。自然。2013;500:301–306. [PMC免费文章][公共医学][谷歌学者]
16Scott M,Klumpp S,Mateescu EM,Hwa T.核糖体合成优化调控中稳健生长规律的出现。分子系统生物学。2014;10:747. [PMC免费文章][公共医学][谷歌学者]
17.Hauryliuk V、Atkinson GC、Murakami KS、Tenson T、Gerdes K。(p)ppGpp在细菌生理学中作用的最新功能见解。《自然微生物评论》。2015;13:298–309. [PMC免费文章][公共医学][谷歌学者]
18Scott M,Gunderson CW,Mateescu EM,Zhang Z,Hwa T.细胞生长和基因表达的相互依赖性:起源和后果。科学。2010;330:1099–1102.[公共医学][谷歌学者]
19Hui S等。定量蛋白质组分析揭示了细菌全球资源分配的简单策略。分子系统生物学。2015;11:784. [PMC免费文章][公共医学][谷歌学者]
20Basan M等人大肠杆菌高效蛋白质组分配的结果。自然。2015;528:99–104. [PMC免费文章][公共医学][谷歌学者]
21.Hermsen R、Okano H、You C、Werner N、Hwa T大肠杆菌在共同利用的碳衬底上。分子系统生物学。2015;11:801. [PMC免费文章][公共医学][谷歌学者]
22Chubukov V,Gerosa L,Kochanowski K,Sauer U。微生物代谢的协调。《自然微生物评论》。2014;12:327–340.[公共医学][谷歌学者]
23Reitzer L.氮同化与全球调节大肠杆菌.微生物年鉴。2003;57:155–176.[公共医学][谷歌学者]
24Dai X等。减少翻译核糖体使大肠杆菌以在缓慢生长期间保持伸长率。自然微生物。2016;2:16231. [PMC免费文章][公共医学][谷歌学者]
25高斯K。核糖体:结构、功能和遗传学。帕克大学出版社;1980年,第693–718页。[谷歌学者]
26Postma PW、Lengeler JW、Jacobson GR.磷酸烯醇丙酮酸:细菌的碳水化合物磷酸转移酶系统。微生物评论。1993;57:543–594. [PMC免费文章][公共医学][谷歌学者]
27Pavlov MY,Ehrenberg M.基因表达的最佳控制,以实现蛋白质组对环境变化的快速适应。美国国家科学院程序。2013;110:20527–20532. [PMC免费文章][公共医学][谷歌学者]
28Giordano N、Mairet F、GouzéJL、Geiselmann J、de Jong H。细胞资源的动态分配是一个优化控制问题:对微生物生长策略的新见解。公共科学图书馆计算生物学。2016;12:e1004802。 [PMC免费文章][公共医学][谷歌学者]
29Csonka LN、Ikeda TP、Fletcher SA、Kustu S。谷氨酸的积累对于鼠伤寒沙门菌在高渗压介质中,但不诱导proU(程序)操纵子。细菌杂志。1994;176:6324–6333. [PMC免费文章][公共医学][谷歌学者]
30Soupene E等人大肠杆菌MG1655菌株:生长缺陷和基因表达的明显交叉调节。细菌杂志。2003;185:5611–5626. [PMC免费文章][公共医学][谷歌学者]
31.Datsenko KA,Wanner BL.染色体基因的一步失活大肠杆菌K-12使用PCR产品。美国国家科学院程序。2000;97:6640–6645. [PMC免费文章][公共医学][谷歌学者]
32.Miller JH公司。分子遗传学实验。冷泉港实验出版社;1972[谷歌学者]
33Benthin S,Nielsen J,Villadsen J.测定细胞RNA含量的简单可靠方法。生物技术。1991;5:39–42. [谷歌学者]
34Deutsch EW等人,《跨太平洋石油管道导游》。蛋白质组学。2010;10:1150–1159. [PMC免费文章][公共医学][谷歌学者]
35克雷格·R,比维斯·RC。串联质谱:用串联质谱匹配蛋白质。生物信息学。2004;20:1466–1467.[公共医学][谷歌学者]
36Sperling E,Bunner AE,Sykes MT,Williamson JR。使用最小二乘傅里叶变换卷积对蛋白质组质谱中同位素分布进行定量分析。分析化学。2008;80:4906–4917. [PMC免费文章][公共医学][谷歌学者]
37Sykes MT,Sperling E,Chen SS,Williamson JR。用质谱法定量核糖体蛋白自动调节网络。分析化学。2010;82:5038–5045. [PMC免费文章][公共医学][谷歌学者]
38Schaechter M,Maalöe O,Kjeldgaard NO.平衡生长过程中细胞大小和化学成分对介质和温度的依赖性鼠伤寒沙门菌.微生物杂志。1958;19:592–606.[公共医学][谷歌学者]
39Basan M等人。通过增加蛋白质合成使细菌细胞膨胀。分子系统生物学。2015;11:836. [PMC免费文章][公共医学][谷歌学者]
40Tong S,Porco A,Isturiz T,Conway T大肠杆菌gntR,gntK(全球导航卫星系统)、和gntU(全球导航卫星)葡萄糖酸代谢的主要系统GntI的基因。细菌杂志。1996;178:3260–3269. [PMC免费文章][公共医学][谷歌学者]
41Voegele RT、Sweet GD、Boos W.甘油激酶大肠杆菌通过与甘油促进剂的相互作用激活。细菌杂志。1993;175:1087–1094. [PMC免费文章][公共医学][谷歌学者]
42Dennis PP,Bremer H大肠杆菌B/r(右)。分子生物学杂志。1974;84:407–422.[公共医学][谷歌学者]