Cực tiểu hoá thay đổi điều hoà (Regulatory On/Off Minimization

ROOM)

4.2.3.1.Giới thiệu về ROOM

Tối thiểu hoá thay đổi điều hoà (ROOM) là một mô hình dùng để dự đoán hành động của mạng trao đổi chất khi gen bị loại bỏ. Nó dựa vào việc tối thiểu hoá số lƣợng của các thay đổi dòng đáng kể so với dạng nguyên thuỷ. ROOM có cùng ý tƣởng nhƣ MOMA nhƣng sử dụng một độ đo khoảng cách khác so với MOMA dùng khoảng cách Euclid. Ngoài ra ROOM còn có khả năng xác định chính xác các đƣờng dẫn khác cho các phản ứng liên kết với các gen bị loại bỏ, điều này càng làm tăng tính tin cậy sinh học của ROOM. ROOM cũng tốt hơn MOMA trong việc dự đoán các dòng trong tế bào và khả năng sát thân của E. coli khi bị đột biến loại bỏ gen.

Một thách thức khác là việc dự đoán tính sát thân và kiểu hình của sinh vật sau khi xóa bỏ gen. Nhƣ đã trình bày trong nghiên cứu của Segre (2002) [4], giả thiết hành động tối ƣu của sinh vật bị xoá bỏ gen có thể không phản ánh chính xác các hành động trong thực tế. Thay vào đó, Segre gợi ý rằng các sinh vật đó điều chỉnh bằng cách tối thiểu hoá các thay đổi phân phối dòng của nó nhƣ cách tiếp cận của MOMA. Phƣơng pháp này cực tiểu theo tiêu chuẩn Euclid của sai khác dòng giữa các mạng trao đổi chất của sinh vật nguyên thuỷ và các dạng đột biến. MOMA đƣợc ghi nhận dự đoán thành công tính sát thân và giá trị dòng trao đổi chất của mạng đột biến của E. coli.

FBA giả thiết về tốc độ phát triển tối ƣu còn MOMA sử dụng độ đo Euclid nhƣng đều không phản ánh chính xác các điều chỉnh của mạng trao đổi chất sau đột biến. Nhƣợc điểm này là do độ đo Euclid có thể ngăn cản các thay đổi lớn của các dòng đơn. Các thay đổi lớn có thể cần để chuyển hƣớng dòng sang các đƣờng phản ứng khác và có thể thấy đƣợc trong thực nghiệm.

Chúng tôi đề nghị thay thế độ đo khoảng cách trong MOMA bằng cực tiểu hoá trạng thái lân cận. Phƣơng pháp mới của chúng tôi, ROOM đơn giản cực tiểu hoá số lƣợng các thay đổi dòng đáng kể. Đặc biệt, ROOM tìm kiếm phân phối dòng cho một giống đột biến mà thoả mãn các ràng buộc ma trận hệ số cân bằng, nhiệt động lực học và hiệu suất dòng, trong khi cực tiểu số lƣợng thay đổi dòng đáng kể so với dạng nguyên thuỷ. ROOM cải tiến dự đoán dòng trao đổi

chất tế bào của E. coli so với MOMA. ROOM cũng dự đoán tốt hơn MOMA về tính sát thân của men bia khi bị xoá bỏ gen ở mức độ lớn.

Kinh nghiệm tăng cƣờng của khoảng cách ROOM dựa trên ý tƣởng:

 Các thay đổi điều hoà gen để thay đổi dòng sau khi gen bị xoá bỏ đƣợc cực tiểu hoá bởi tế bào khi cố gắng cực tiểu hoá chi phí thích nghi.

 Những thay đổi điều hoà này có thể miêu tả bởi các giá trị on/off. ROOM hoàn toàn mô tả các thay đổi dòng bằng cách xác định các thay đổi dòng đáng kể.

4.2.3.2.Phƣơng pháp ROOM

ROOM tìm kiếm một phân phối dòng thoả mãn cùng ràng buộc nhƣ FBA trong khi cực tiểu hoá số lƣợng thay đổi dòng đáng kể. Chúng ta chỉ tính các thay đổi dòng đáng kể vì luôn tồn tại các nhiễu cố hữu trong các hệ thống sinh học và cũng là để giảm thời gian xử lý. Chúng ta dùng MILP nhƣ sau:

 Minimize yi i1 m  subject to Sv = 0 vmin v vmax vj 0,j  A vi yi(vmax,iwiu)wiu (*) vi yi(vmin,i wil) wil (**) yi {0,1} (***) wiu wiwi  wil wiwi              

với mỗi dòng i, 1<= i <= m, yi = 1 nếu vi có thay đổi dòng đáng kể và yi = 0 trong trƣờng hợp còn lại, wu

and wl là các ngƣỡng để xác định thay đổi dòng, với

δ, ε là các khoảng dung sai tƣơng đối và tuyệt đối tƣơng ứng với w và A trong MOMA.

Thật vậy, với yi = 1, bất đẳng thức (*) và (**) không tạo ra ràng buộc mới

với vi, trong khi nếu yi = 0 thì đó là các ràng buộc mới. Kích thƣớc của δ, ε ảnh hƣởng đến thời gian chạy của thuật toán MILP; chúng tôi chọn các giá trị nhỏ nhất mà vẫn có kết quả trong thời gian có thể chấp nhận. (Đặc biệt, chúng tôi

chọn δ = 0.03, ε =0.001 để dự đoán dòng và δ = 0.1, ε =0.01 để dự đoán tính sát thân). Sự lựa chọn các tham số ảnh hƣởng đến kết quả phân phối dòng của giải thuật bởi nó cho phép cộng (δ) và nhân (ε) một lƣợng nhỏ dòng để điều chỉnh các đƣờng thay thế mà không mất chi phí. Thời gian chạy với bài toán MILP chỉ là vài giây.

Khi bỏ qua các ràng buộc số nguyên trong (***) thành 0<=yi <=1 thì bài toán trở thành biến thể của quy hoạch tuyến tính. Với biến thể này, các tham số δ, ε đƣợc cho bằng 0. Dự đoán dựa trên biến thể LP của ROOM cũng khá chính xác dù kém hơn so với ROOM ban đầu. Công thức MILP dùng để dự đoán dòng của E. coli và dự đoán tính sát thân của men, trong khi LP dùng để dự đoán tính sát thân và tốc độ phát triển của E. coli. Điều thú vị là công thức LP của ROOM đôi khi tƣơng tự với một biến thể của MOMA.

4.2.3.3.So sánh MOMA và ROOM

Dự đoán dòng

Các số liệu dòng trao đổi chất của E. coli bị xoá gen có trong báo cáo của Emmerling [17]. Bộ sƣu tập 17 dòng nội bào trong mạng trao đổi chất carbon trung tâm của E. coli đƣợc xác định theo kinh nghiệm bằng cách kết hợp quang phổ NMR trong các thí nghiệm đánh dấu bằng 13C (đồng vị carbon) và các số liệu sinh lý học (xem thêm ở tài liệu tham khảo). Các dòng của cả dạng nguyên thuỷ (JM101) và dạng đột biến pyk (PB25) đều đƣợc đo với 3 điều kiện phát

triển khác nhau để cho ra sáu tập kết quả. Hai điều kiện hạn chế glucose đƣợc dùng (nồng độ thấp / cao) và một điều kiện hạn chế nitrogen (nitơ).

Để mô hình hoá dạng đột biến pyk (PB25, pykA::kan pykF::cat), hai phản ứng từ Phosphoenolpyruvate (PEP) tới Pyruvate (PYR) đƣợc đặt bằng không.

Trong báo cáo của Segre (2002) [4] đã so sánh dự đoán của FBA và MOMA với tập hợp các dòng thực nghiệm này. Với dạng hoang dã FBA dự đoán tốt trong cả 3 trƣờng hợp, trong khi với dạng đột biến pyk kết quả dự đoán tồi cho cả 2 trƣờng hợp môi trƣờng glucose thấp và nitrogen thấp. MOMA dự đoán cho dạng đột biến tốt hơn nhiều so với FBA.

Dự đoán của ROOM cho dạng đột biến pyk chính xác hơn đáng kể so với FBA và MOMA. Với cả 2 thí nghiệm glucose thấp và nitrogen hạn chế, độ tƣơng quan của ROOM với dữ liệu dòng thực nghiệm cao hơn 12% so với MOMA (Hình 4.10A). Hình 4.10 B,C so sánh độ chính xác dự đoán của MOMA và ROOM.

Hình 4.9 So sánh FBA, MOMA và ROOM [8].

Hình (A) là độ tƣơng quan giữa dự đoán của FBA, MOMA và ROOM với các dữ liệu thực nghiệm cho một dòng không nhân trong 3 điều kiện phát triển khác nhau. Hình (B,C) thể hiện so sánh MOMA và ROOM trong điều kiện glucose thấp. Các dòng đƣợc thể hiện bằng phần trăm glucose sử dụng.

Ngoài ra, với điều kiện glucose cao, ROOM tìm thấy một nghiệm khả thi mà chỉ các dòng trên đƣờng thay thế của các phản ứng bị bỏ là thay đổi đáng kể (đƣờng v16-v14v15 trong hình 4.11). Thật vậy, trong môi trƣờng glucose cao, dữ liệu thực nghiệm cho thấy hơn 3 dòng có chênh lệch lớn nhất giữa dạng hoang dã và dạng đột biến (Emmerling [17]). MOMA dự đoán thay đổi dòng nhỏ hơn trên đƣờng này, từ đó cho thấy phƣơng pháp cực tiểu hóa khoảng cách Euclid có thể không cho phép các thay đổi lớn.

Hình 4.10 Giản đồ miêu tả mạng trao đổi chất carbon trung tâm của E. coli [8].

Các dòng thực nghiệm gán nhãn từ v1–v24 và đƣợc đo theo hƣớng mũi tên. Chú ý rằng trong một số trƣờng hợp các dòng theo hƣớng ngƣợc lại đƣợc thể hiện bằng các giá trị âm. Mũi tên đậm thể hiện các đƣờng thay thế ngắn (v19– v18–v21) đƣợc dự đoán bởi ROOM với gen pyk bị loại bỏ (đánh dấu X) trong

điều kiện glucose trung bình. Các thay đổi lớn của flux trên đƣờng này đƣợc xác định qua thực nghiệm.

Một đánh giá rộng bằng mô phỏng việc xoá bỏ gen của Saccharomyces cerevisiae đã đƣợc trình bày bởi Forster et al. (2003). FBA đƣợc dùng để tính tốc độ phát triển tối đa cho 555 gen xoá bỏ. Một gen đƣợc cho là sát thân nếu tốc độ phát triển tối đa của dạng đột biến tƣơng ứng nhỏ hơn 5% dạng hoang dã. FBA dự đoán chính xác tính sát thân của 89.6% số gen. Chúng tôi kiểm tra cả MOMA và ROOM trên mô hình này để nghiên cứu một tập hợp các gen có liên quan với nhau.

Các gen đƣợc dự đoán là sát thân bởi FBA cũng đƣợc dự đoán bởi MOMA và ROOM vì tốc độ phát triển tối đa tính đƣợc bởi FBA chỉ có thể giảm với MOMA và ROOM. Cùng chạy MOMA với các gen có thể dự đoán bởi FBA cho kết quả tồi hơn đáng kể hơn FBA. Đặc biệt, ngoài các gen đó thì MOMA dự đoán nhầm 50 gen là sát thân (hình 4.12). Ngƣợc lại, kết quả của ROOM trên các gen đó tƣơng đƣơng hoàn toàn với FBA.

Hình 4.11 Danh sách 50 gen MOMA dự đoán nhầm.

Chúng tôi chuyển sang kiểm tra 50 gen mà ROOM và FBA đều dự đoán đúng còn MOMA dự đoán sai. Nhƣ thực nghiệm E.coli trong môi trƣờng

glucose cao, ROOM tìm ra một đƣờng thay thế ngắn nhất cho các phản ứng bị loại bỏ. Dự đoán sai của MOMA khiến nó đƣa ra một lời giải sai.

4.2.3.4.Mở rộng của ROOM

Như đã trình bày ở trên, các con đường khác đôi được sử dụng bởi các giống đột biến để chuyển hướng dòng. Do đó các mô hình dựa trên ràng buộc cần phải chứa cả các đường thay thế hay các đặc tính cấu trúc mạng như là khoảng cách từ phản ứng bị xoá bỏ. Để làm điều này, một biến thể của ROOM có thể gắn các giá trị chi phí khác nhau cho các thay đổi dòng đáng kể

của các phản ứng khác nhau. Đặc biệt, hàm mục tiêu có thể được thay đổi để

cực tiểu 

m i ciyi

1 , với ci là chi phí để thay đổi dòng i. Ví dụ, để tính cho các đường thay thế, mỗi phản ứng trên một đường thay thế được gán một chi phí tỷ lệ với độ dài đường. Các phản ứng còn lại sẽ được gán chi phí cao hơn.

Thực tế, chúng tôi thấy rằng khi gắn các chi phí lớn hơn cho các thay đổi dòng đáng kể không sử dụng biến thể của ROOM, các kết quả dự đoán có tốt hơn một chút. Dự đoán này cho thấy để điều hoà một gen không phải nguyên thuỷ

CHƢƠNG 5 - BÀI TOÁN CỤ THỂ VÀ PHƢƠNG PHÁP TIẾP CẬN