Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Enzyme acetyl-CoA carboxylase (ACC) (EC 6.2) là một enzyme dị phân tử phụ thuộc biotin, giữ vai trò xúc tác trung tâm trong quá trình điều hòa sinh tổng hợp và oxy hóa axit béo, đồng thời duy trì cân bằng nội môi năng lượng ở người. ACC xúc tác phản ứng carboxyl hóa acetyl-CoA để tạo thành malonyl-CoA – tiền chất trực tiếp cho quá trình tổng hợp axit béo và là chất ức chế tự nhiên đối với carnitine palmitoyltransferase I (CPT1), qua đó kìm hãm con đường $\beta$-oxy hóa axit béo.

Ở cơ thể người, ACC tồn tại dưới hai dạng đồng phân chính gồm ACC1 và ACC2 với sự khác biệt về khối lượng phân tử, phân bố mô và chức năng tế bào:

  • ACC1: Khối lượng phân tử 265 kDa, gồm 2346 axit amin, khu trú chủ yếu tại bào tương của các mô giàu lipid như gan và mô mỡ, chịu trách nhiệm xúc tác sinh tổng hợp axit béo mới (de novo lipogenesis - DNL).
  • ACC2: Khối lượng phân tử 280 kDa, gồm 2458 axit amin (chứa thêm đoạn peptide 114 axit amin ở đầu tận cùng N so với ACC1), phân bố chủ yếu tại màng ngoài ty lạp thể của mô cơ tim và cơ xương, tham gia kìm hãm quá trình $\beta$-oxy hóa thông qua việc điều hòa nồng độ malonyl-CoA ức chế CPT1.

Khi cơ chế chuyển hóa bị rối loạn, lượng axit béo dư thừa tích tụ sẽ dẫn tới hội chứng chuyển hóa (Metabolic Syndrome - MetS). MetS đặc trưng bởi rối loạn lipid máu, béo bụng (tiêu chuẩn chẩn đoán khi vòng eo > 94 cm ở nam và > 80 cm ở nữ), gan nhiễm mỡ không do rượu, kháng insulin, làm gia tăng tỷ lệ mắc đái tháo đường týp 2, tăng huyết áp và biến cố tim mạch. Theo số liệu của Tổ chức Y tế Thế giới (WHO), MetS gây ra khoảng 1,9 triệu ca tử vong mỗi năm tại các quốc gia có thu nhập thấp và trung bình; tỷ lệ hiện mắc tại Nam Phi đạt 21,8% (15,6% ở nam và 24,8% ở nữ), và tỷ lệ mắc ở thanh thiếu niên béo phì lên tới gần 30%. Người mắc MetS có nguy cơ mắc đái tháo đường týp 2 tăng gấp 5 lần, nguy cơ phát triển bệnh tim mạch tăng gấp 2 lần trong vòng 5–10 năm, nguy cơ đột quỵ tăng 2–4 lần, nhồi máu cơ tim tăng 3–4 lần và nguy cơ tử vong tăng gấp 2 lần so với người bình thường.

Việc ức chế đồng thời hoặc chọn lọc ACC1 (nhằm giảm sinh tổng hợp axit béo) và ACC2 (nhằm kích thích quá trình oxy hóa axit béo) được xác định là mục tiêu điều trị trúng đích đối với MetS. Tuy nhiên, các nghiên cứu tính toán và thực nghiệm trên thế giới giai đoạn 2015–2021 chủ yếu khảo sát trên mô hình enzyme của động vật (chuột), vi khuẩn hoặc nấm men. Luận án chỉ ra khoảng trống nghiên cứu thực tế:

  • Chưa có hệ thống mô hình hóa đầy đủ về tương tác phân tử giữa các nhóm chất ức chế với cấu trúc enzyme ACC1 và ACC2 ở người trên cả ba vùng chức năng: biotin carboxylase (BC), protein vận chuyển biotin carboxyl (BCCP) và carboxyltransferase (CT).
  • Thiếu các mô hình định lượng cấu trúc – hoạt tính (QSAR) kết hợp các công cụ học máy phi tuyến và logic mờ nhằm làm sáng tỏ các yếu tố cấu trúc và nhóm chức quyết định hiệu quả ức chế ACC ở người.

Mục tiêu nghiên cứu

  • Mục tiêu chính: Ứng dụng các phương pháp hóa tính toán, nghiên cứu mô hình ức chế ACC ở người và xác định các yếu tố cấu trúc ảnh hưởng đến hoạt tính ức chế ACC của các nhóm chất ức chế, góp phần định hướng, hỗ trợ điều trị MetS với các tác nhân ức chế hiệu quả hơn.
  • Mục tiêu cụ thể (được đánh số theo văn bản):
    1. Khảo sát các mô hình tương tác giữa ACC1 và ACC2 ở người với các nhóm hợp chất ức chế đã có $IC_{50}$ thực nghiệm. Xác định các tâm axit amin quan trọng, vị trí gắn kết thuận lợi, hốc phản ứng, vùng gắn kết thuận lợi trên 2 enzyme. Đánh giá loại liên kết, tương tác chủ yếu giữa ligand với các phân tử axit amin trên enzyme.
    2. Khảo sát mối liên hệ định lượng cấu trúc – hoạt tính (QSAR), xác định dạng cấu trúc/nhóm chức chủ yếu có ảnh hưởng quan trọng đến hoạt tính ức chế enzyme.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu:
    • Cấu trúc 3D của hai đồng phân enzyme ACC1 và ACC2 ở người, được thu thập từ Ngân hàng Dữ liệu Protein (RCSB PDB), bao gồm các mã PDB: 2YL2 (chuỗi axit amin 78–617), 4ASI (chuỗi axit amin 1493–2260) cho ACC1; 3JRX (chuỗi axit amin 217–775), 3FF6 (chuỗi axit amin 1693–2450) và 4JRX cho ACC2.
    • Tập hợp 21 hợp chất ức chế ACC1 và 29 hợp chất ức chế ACC2 đã có giá trị hoạt tính thực nghiệm $\log(IC_{50})$ được công bố để phục vụ khảo sát tương tác enzyme – ligand.
    • Tập hợp 130 hợp chất thơm (hương phương) có hoạt tính ức chế enzyme thực nghiệm từ cơ sở dữ liệu mở của Đại học California – San Diego (UCSD) để xây dựng và kiểm định mô hình QSAR.
  • Phạm vi nghiên cứu: Nghiên cứu lý thuyết và thực nghiệm tính toán hóa tin (in silico), sử dụng các kỹ thuật docking phân tử, mô phỏng động lực học phân tử (MD) trong lồng dung môi nước ở $37^\circ\text{C}$, tính toán cơ học lượng tử/cơ học phân tử (QM/MM) theo phương pháp ONIOM, hồi quy tuyến tính đa biến (MLR), giải thuật di truyền (GA), mạng nơron nhân tạo (ANN) và hệ logic mờ (Fuzzy Logic - FL).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước được điểm lại trong văn bản

Tác giả đã hệ thống hóa các công trình nghiên cứu trong và ngoài nước liên quan đến chất ức chế ACC:

  • Harwood và cộng sự: Thí nghiệm trên chuột, chuột hamster và khỉ Rhesus với axit 5-(tetradecyloxy)-2-furancarboxylic (TOFA). Ở môi trường nội bào, TOFA chuyển hóa thành dạng acyl-CoA thioester ức chế ACC, làm giảm tổng hợp axit béo, giảm nồng độ cholesterol và triglyceride huyết tương. Đồng thời, nghiên cứu chất ức chế không chọn lọc CP-640186 (Pfizer) chứng minh khả năng kìm hãm phản ứng chuyển vị nhóm carboxyl, làm giảm malonyl-CoA trong tế bào gan và tăng cường oxy hóa axit béo toàn cơ thể.
  • Matthew và cộng sự (2014): Tổng kết các nhóm chất ức chế ACC thực nghiệm trên nấm men và động vật:
    • Dẫn xuất bipiperidylcarboxamides (CP-640186): Gắn vào tâm hoạt động carboxyltransferase (CT), che chắn vị trí chuyển nhóm carboxyl và tương tác một phần với tâm gắn ATP của chuỗi biotin carboxylase (BC).
    • Nhóm polyketide (Soraphen A): Dựa trên nghiên cứu của Hoffmann-LaRoche; Shen và cộng sự xác định Soraphen A gắn vào tâm BC của nấm men tại vị trí cách tâm gắn ATP khoảng $25\text{ \AA}$.
    • Nhóm aryloxyphenoxypropionate và cyclohexanedione: Gồm haloxyfop, diclofop và sethoxydim. Zhang và cộng sự đã đồng kết tinh haloxyfop và diclofop với tâm CT của ACC nấm men, chỉ ra vị trí gắn kết lân cận nhóm sulfhydryl của coenzyme A.
    • Nhóm chuỗi axit béo acyl-CoA tổng hợp: TOFA, MEDICA 16, ESP-55016, S2E, chuyển hóa thành dạng acyl-CoA thioester gây ức chế CT nhưng độ đặc hiệu chưa cao do tác động chéo lên các enzyme phụ thuộc acyl-CoA khác.
    • Nhóm ức chế cạnh tranh cơ chất: Waldrop và cộng sự phát triển chloroacetylate biotin (CABI), liên kết cộng hóa trị giữa CoA và biotin thông qua cầu nối acyl giữa nhóm thiol của CoA và $1'\text{-N}$ của biotin, ức chế cạnh tranh quá trình hình thành malonyl-CoA.
    • Dẫn xuất acylsulfonamides: Nghiên cứu của Fujirubio và Ajinomoto cho thấy khả năng ức chế ACC ở dải nồng độ $10\text{--}100\ \mu\text{M}$.
    • Các dẫn xuất khác: Các cấu trúc ức chế do Taisho, Sanofi-Aventis, Boehringer Ingelheim, AstraZeneca, Takeda, Pfizer (spirochromanone, spirocyclic salicylamide, spirolactam), Haselkorn, Nimbus (thienopyrimidine), Amgen (oxadiazole) đề xuất trên mô hình chuột và nấm men.
  • Jones và cộng sự (2017): Chứng minh việc ức chế ACC giúp giảm DNL và tăng tốc độ chuyển hóa tế bào, kìm hãm sự phát triển của khối u mỡ.
  • Huang và cộng sự (2015, 2016, 2017): Tổng hợp hai nhóm hợp chất mang khung quinoline có hoạt tính ức chế ACC; hợp chất 7a đạt $IC_{50}$ trên ACC1 là $189\text{ nM}$ và trên ACC2 là $172\text{ nM}$. Tác giả sử dụng docking phân tử xác định kiểu gắn kết và phát triển dẫn xuất tối ưu PP-7a điều trị rối loạn chuyển hóa trên mô hình chuột.
  • Tống Thị Thu Cúc: Áp dụng mô phỏng tương tác protein – ligand kết hợp phương pháp QM/MM nghiên cứu cơ chế thủy phân acetylcholine của acetylcholinesterase, xác định vị trí gắn kết của donepezil, galantamin, tacrin, neostigmin, physostigmin và rivastigmin.

Khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết các hạn chế tồn đọng từ các nghiên cứu trước:

  1. Chuyển trọng tâm khảo sát từ mô hình enzyme trên nấm men/chuột sang cấu trúc tinh thể enzyme ACC1 và ACC2 hoàn chỉnh ở người, bao quát cả 3 vùng cấu trúc chức năng (BC, BCCP, CT).
  2. Định vị chính xác các hốc phản ứng, các tâm gắn kết axit amin then chốt, xác định năng lượng gắn kết và phân tích tương tác năng lượng QM/MM đối với các nhóm ligand đại diện.
  3. Ứng dụng mô phỏng động lực học phân tử trong môi trường lỏng TIP3P ở nhiệt độ sinh lý $37^\circ\text{C}$ để đánh giá độ bền liên kết phức hợp enzyme – ligand.
  4. Xây dựng mô hình QSAR tích hợp mạng nơron nhân tạo và logic mờ nhằm lượng hóa mức độ và xu hướng tác động của từng biến cấu trúc lên hoạt tính sinh học.

Cơ sở lý thuyết và phương pháp nghiên cứu

Lý thuyết, khái niệm và khung phân tích

  • Cấu trúc phân vùng của ACC: Gồm 3 vùng chức năng BC, BCCP và CT.
  • Lý thuyết hóa học lượng tử: Sử dụng phương pháp bán kinh nghiệm PM3 và lý thuyết phiếm hàm mật độ (DFT) ở mức lý thuyết B3LYP/6-31G(d,p) để tối ưu hóa cấu trúc hình học và tính toán các thông số orbital biên HOMO, LUMO, moment lưỡng cực, năng lượng gắn kết.
  • Lý thuyết docking phân tử: Sử dụng các thuật toán tìm kiếm không gian cấu hình gồm gia nhiệt mô phỏng (simulated annealing), thuật giải di truyền (genetic algorithm - GA), thuật giải di truyền Lamarckian (Lamarckian GA) trong AutoDock và hàm tính điểm lai hóa (hybrid scoring function) trong AutoDock Vina.
  • Phương pháp QM/MM (ONIOM): Phân tách hệ thống phức hợp thành 2 lớp tính toán: lớp cao (vùng tương tác trực tiếp giữa ligand và các gốc axit amin tạo liên kết) tính bằng cơ học lượng tử PM3; lớp thấp (các gốc axit amin còn lại trong bán kính quét $5\text{ \AA}$) tính bằng cơ học phân tử với trường lực UFF.
  • Lý thuyết xử lý thống kê cổ điển: Phân tích hồi quy tuyến tính đa biến (MLR), hệ số tương quan ($r$), hệ số xác định ($R^2$), khả năng dự đoán mô hình ($Q^2$), sai số tuyệt đối trung bình (MAE), kiểm định ý nghĩa thống kê $t$ (mức ý nghĩa $p < 0{,}05$), kiểm định Fisher ($F$), kiểm định Mann – Whitney và kiểm định Kolmogorov – Smirnov để so sánh hai mẫu độc lập.
  • Lý thuyết trí tuệ nhân tạo và thống kê tin sinh:
    • Mạng nơron nhân tạo (ANN): Mô hình mạng nhiều lớp sử dụng thuật toán lan truyền ngược (backpropagation training), hàm kích hoạt phi tuyến logistic, cấu trúc mạng ANN 4-6-1.
    • Phương pháp logic mờ (Fuzzy Logic - FL): Khái niệm tập mờ, hàm liên thuộc (membership function) $\mu(x) \in [0, 1]$, hệ quy tắc điều kiện IF-THEN để mô tả quan hệ phi tuyến giữa các thông số cấu trúc đầu vào và hoạt tính $\log(IC_{50})$ đầu ra.
    • Giải thuật di truyền (GA): Tối ưu hóa tổ hợp biến cấu trúc dựa trên cơ chế tiến hóa sinh học (chọn lọc, đột biến, trao đổi chéo).

Phương pháp thực tế, công cụ và quy trình phân tích

Quy trình tính toán thực tế trong luận án được triển khai qua các bước:

Nhóm nội dung Công cụ / Phần mềm Phương pháp / Kỹ thuật chi tiết
Tối ưu hóa cấu trúc ligand Gaussian 09w, Hyperchem 8.10 Bán kinh nghiệm PM3, DFT B3LYP/6-31G(d,p)
Docking phân tử & Sàng lọc ảo AutoDock 4, AutoDock Vina, Raccoon, Cygwin Sàng lọc phối tử trên hốc phản ứng cục bộ và phủ toàn bộ enzyme; định dạng pdbqt
Phân tích tương tác không gian Discovery Studio Visualizer 2017 R2 Client Xác định vị trí gắn kết, liên kết hydro, tương tác kỵ nước, khoảng cách tương tác
Động lực học phân tử (MD) UCSF Chimera 1.10 Mô phỏng trong lồng dung môi nước TIP3P ở nhiệt độ cơ thể ($37^\circ\text{C}$)
Tính toán năng lượng tương tác Gaussian 09w (ONIOM) Khảo sát QM(PM3)/MM(UFF) trong bán kính quét $5\text{ \AA}$
Trích xuất thông số cấu trúc PaDEL-Descriptor 2.18, Hyperchem 8.10 17 thông số hóa lượng tử + 16018 thông số 2D/3D
Sàng lọc biến & Tối ưu hóa Statgraphics Centurion XVI, phần mềm ĐH Jadavpur Phân tích MLR ($p < 0{,}05$), giải thuật di truyền GA nhiều chu kỳ
Xây dựng mô hình phi tuyến R-console 3.2.3, Fuzzy Logic Toolbox ANN lan truyền ngược 5000 vòng lặp (85% train, 15% test; $r_{\text{test}} > 0{,}90$, $\text{SSE} < 9$, $\text{MSE}_{\text{test}} < 0{,}15$), chuyển đổi hệ mờ ANFIS

Nguồn tư liệu và dữ liệu

  • Cơ sở dữ liệu cấu trúc protein RCSB PDB: ACC1 (2YL2, 4ASI) và ACC2 (3FF6, 3JRX, 4JRX).
  • Bộ dữ liệu hoạt tính sinh học của 21 chất ức chế ACC1 và 29 chất ức chế ACC2 từ công trình của Matthew và cộng sự (2014).
  • Bộ dữ liệu mở gồm 130 hợp chất thơm có hoạt tính ức chế thực nghiệm của Đại học California – San Diego (UCSD).

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 hệ thống hóa cơ sở sinh học phân tử của ACC và MetS, đồng thời tổng hợp các chất ức chế ACC đã công bố. Văn bản cung cấp các đặc điểm cấu trúc của ACC1 và ACC2 ở người:

Đặc điểm cấu trúc Vị trí trên ACC1 Axit amin / Phối tử trên ACC1 Vị trí trên ACC2 Axit amin / Phối tử trên ACC2
Tâm gắn kim loại 424–424 Axit glutamic ($\text{Mn}^{2+}$ 1) 567–567 Axit glutamic ($\text{Mn}^{2+}$ 1)
Tâm gắn kim loại 437–437 Axit glutamic ($\text{Mn}^{2+}$ 1) 580–580 Axit glutamic ($\text{Mn}^{2+}$ 1)
Tâm gắn kim loại 437–437 Axit glutamic ($\text{Mn}^{2+}$ 2) 580–580 Axit glutamic ($\text{Mn}^{2+}$ 2)
Tâm gắn kim loại 439–439 Axit aspartic ($\text{Mn}^{2+}$ 2) 582–582 Axit aspartic ($\text{Mn}^{2+}$ 2)
Tâm hoạt động 441–441 Phức hợp Arginine + 6 axit amin (LEU365, THR345, ASN362, PRO363, GLU137, TRP132) 584–584 Phức hợp Arginine + 6 axit amin (GLU146, PRO373, VAL377, LEU375, GLN376, THR355)
Tâm gắn Coenzyme A 1823–1823 Arginine (Coenzyme A) 1934–1934 Arginine (Coenzyme A)
Tâm gắn Coenzyme A 2127–2127 Lysine (Coenzyme A) 2238–2238 Lysine (Coenzyme A)
Tâm gắn Coenzyme A 2129–2129 Arginine (Coenzyme A) 2240–2240 Arginine (Coenzyme A)
Tâm gắn nucleotide 315–320 Glycine, Lysine (ATP) 458–463 Glycine, Lysine (ATP)

Chương này cũng phân loại các nhóm chất ức chế: dẫn xuất bipiperidylcarboxamides (CP-640186), polyketide (Soraphen A), aryloxyphenoxypropionate (haloxyfop, diclofop), cyclohexanedione (sethoxydim), chuỗi acyl-CoA tổng hợp (TOFA, MEDICA 16, ESP-55016, S2E), chất ức chế cạnh tranh (CABI-CoA), acylsulfonamides, cùng các cấu trúc đề xuất bởi AstraZeneca, Takeda, Pfizer (spirochromanone, spirocyclic salicylamide, spirolactam), Boehringer Ingelheim, Haselkorn, Nimbus và Amgen.

Chương 2: Thực nghiệm tính toán

Chương 2 trình bày chi tiết quy trình tính toán hai nhóm đối tượng:

  • Nhóm 1: Khảo sát tương tác docking, MD và QM/MM trên 21 hợp chất ức chế ACC1 và 29 hợp chất ức chế ACC2 có giá trị thực nghiệm $\log(IC_{50})$.

Dữ liệu phối tử ức chế ACC1 và ACC2 được tác giả thống kê như sau:

Tên hợp chất $\log(IC_{50})$ trên ACC1 $\log(IC_{50})$ trên ACC2 Tên hợp chất $\log(IC_{50})$ trên ACC1 $\log(IC_{50})$ trên ACC2
Taisho1 2,00 1,36 Takeda3 2,08 1,26
Taisho2 1,52 1,52 Takeda4 2,20 0,93
Taisho3 1,87 1,87 Takeda5 2,11 1,36
Taisho4 1,38 1,38 Takeda6 1,32 0,69
Taisho5 2,28 1,98 Pfizer3 — 1,04
Sanofi1 — 2,80 Pfizer4 3,40 3,40
Sanofi2 2,28 1,48 Pfizer5 2,07 1,88
Astra1 — 3,40 Pfizer6 1,00 0,60
Astra2 — 3,56 Boe1 — 1,70
Astra3 — 2,78 Boe2 — 2,26
Astra4 — 2,32 Nimbus1 1,00 1,00
Astra5 4,11 3,11 Amgen1 3,08 1,57
Astra6 3,08 2,04 Amgen2 2,28 1,73
Astra7 3,00 2,00 Takeda1 — 0,73
Takeda2 1,50 —
  • Nhóm 2: 130 hợp chất thơm từ nguồn dữ liệu UCSD để xây dựng mô hình QSAR thông qua các bước: trích xuất 17 thông số hóa lượng tử và 16018 thông số PaDEL, sàng lọc MLR ($p < 0{,}05$), tối ưu biến bằng giải thuật di truyền (GA), huấn luyện mạng ANN lan truyền ngược 5000 chu kỳ trên tập huấn luyện (85%) và tập kiểm tra (15%), thiết lập hệ logic mờ (ANN 4-6-1) và trích xuất quy tắc mờ để dự đoán xu hướng ảnh hưởng.

Chương 3: Kết quả - Bàn luận

Chương 3 trình bày các phát hiện về tương tác phối tử – enzyme và mô hình QSAR:

  • Tương tác enzyme – phối tử:
    • Xác định các hốc phản ứng không gian gồm các hốc A đến E trên ACC1 và các hốc F đến I trên ACC2.
    • Phân tích tương quan giữa năng lượng gắn kết (Binding Energy - BE) và hoạt tính thực nghiệm $\log(IC_{50})$.
    • Khảo sát mô phỏng động lực học phân tử (MD) trong môi trường cơ thể ($37^\circ\text{C}$, TIP3P) xác nhận sự ổn định cấu dạng của các phức hợp tiêu biểu:
      • Trên ACC1: Phức hợp với Pfizer6 đạt $\text{BE} = -9{,}30\text{ kcal/mol}$; phức hợp với Astra5 đạt $\text{BE} = -6{,}40\text{ kcal/mol}$.
      • Trên ACC2: Phức hợp với Pfizer6 đạt $\text{BE} = -10{,}90\text{ kcal/mol}$; phức hợp với Astra1 đạt $\text{BE} = -7{,}27\text{ kcal/mol}$.
    • Đánh giá sự biến đổi liên kết và tương tác axit amin tại các hốc phản ứng trước và sau mô phỏng động lực học phân tử.
    • Kết quả tính toán ONIOM QM(PM3)/MM(UFF) làm rõ sự phân bố mật độ điện tử và năng lượng tương tác tại tâm gắn kết của các ligand Soraphen A và CP-640186 (các dạng A, B, C, D).
  • Mô hình QSAR và logic mờ:
    • Chọn lọc được 10 thông số cấu trúc có tương quan thống kê chặt chẽ nhất ($p < 0{,}05$) với hoạt tính ức chế, sắp xếp theo thứ tự độ lớn trọng số tuyệt đối $|W|$ giảm dần.
    • Xây dựng thành công mô hình hệ mờ lai hóa mạng nơron cho ACC1 (theo phương pháp PM3) và ACC2 (theo phương pháp DFT).
    • So sánh đối chiếu giữa $\log(IC_{50})$ thực nghiệm và $\log(IC_{50})$ tính toán qua mô hình ANFIS/FL cho thấy sự tương thích thống kê.
    • Xác định các khoảng giá trị tối ưu của từng biến cấu trúc đối với ACC1 và ACC2, từ đó thiết lập quy tắc dự đoán xu hướng nâng cao hoạt tính ức chế enzyme.

Chương 4: Kết luận – Hướng phát triển

  • Tóm lược các vị trí gắn kết thuận lợi, các hốc phản ứng và các tương tác axit amin chính trên ACC1 và ACC2 ở người.
  • Đánh giá vai trò của mô hình QSAR lai hóa nơron – mờ trong việc lượng hóa các yếu tố cấu trúc ảnh hưởng đến khả năng ức chế ACC1 và ACC2.
  • Định hướng tiếp tục mở rộng mô hình tính toán cho các cấu trúc phối tử mới và kết hợp kiểm chứng thực nghiệm.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và khoa học

  1. Thiết lập hệ thống mô hình tương tác phân tử 3D hoàn chỉnh giữa hai đồng phân enzyme ACC1 và ACC2 ở người với các nhóm chất ức chế đã có số liệu thực nghiệm, bao quát đầy đủ cả ba vùng chức năng BC, BCCP và CT.
  2. Xác định các hốc phản ứng đặc trưng (hốc A–E trên ACC1 và F–I trên ACC2), định vị các tâm axit amin then chốt tham gia hình thành liên kết hydro và tương tác kỵ nước, đồng thời đánh giá sự thay đổi cấu dạng phức hợp trong điều kiện mô phỏng nhiệt độ sinh lý ($37^\circ\text{C}$).
  3. Ứng dụng phương pháp ONIOM QM/MM để tính toán năng lượng tương tác điện tử chi tiết tại tâm gắn kết, làm rõ cơ chế tương tác ở mức độ cơ học lượng tử.
  4. Xây dựng quy trình tích hợp giữa giải thuật di truyền (GA), mạng nơron nhân tạo (ANN) và logic mờ (FL) để xử lý dữ liệu QSAR phi tuyến, cho phép mô hình hóa mối quan hệ giữa cấu trúc phân tử và hoạt tính ức chế sinh học $\log(IC_{50})$.

Đóng góp mới về mặt thực tiễn và giải pháp đề xuất

  1. Xác định năng lượng gắn kết và mức độ ổn định của các phối tử tiêu biểu, trong đó Pfizer6 thể hiện ái lực liên kết cao trên cả ACC1 ($\text{BE} = -9{,}30\text{ kcal/mol}$) và ACC2 ($\text{BE} = -10{,}90\text{ kcal/mol}$).
  2. Sàng lọc và thiết lập bảng trọng số $|W|$ của 10 thông số cấu trúc quyết định hoạt tính ức chế enzyme.
  3. Xác định các khoảng giá trị tối ưu cho từng biến cấu trúc trên ACC1 và ACC2, cung cấp tiêu chí định lượng phục vụ định hướng thiết kế và sàng lọc ảo các dẫn xuất thuốc mới điều trị hội chứng chuyển hóa MetS.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  • Nghiên cứu được thực hiện hoàn toàn dựa trên các phương pháp hóa tin và mô phỏng tính toán (in silico).
  • Dữ liệu hoạt tính sinh học phục vụ mô phỏng docking và xây dựng QSAR phụ thuộc vào các nguồn dữ liệu thực nghiệm đã công bố (tập 21 chất trên ACC1, 29 chất trên ACC2 và 130 chất vòng thơm từ UCSD), chưa tự thực hiện các phép đo sinh hóa thực nghiệm trực tiếp trên mẫu sinh phẩm người.

Hướng nghiên cứu tiếp theo

  • Dựa trên các khoảng giá trị cấu trúc tối ưu đã xác định từ mô hình logic mờ, tiến hành thiết kế cấu trúc các phân tử ức chế mới có độ chọn lọc cao trên ACC1 hoặc ACC2.
  • Thực hiện tổng hợp hóa học các hợp chất được dự đoán có hoạt tính cao và tiến hành khảo sát thực nghiệm sinh học in vitro, in vivo để kiểm chứng hiệu quả ức chế enzyme ACC ở người và khả năng cải thiện hội chứng chuyển hóa.

Giá trị tham khảo

  • Đối tượng thụ hưởng:
    • Nghiên cứu sinh, học viên cao học và giảng viên thuộc các chuyên ngành Hóa lý thuyết và Hóa lý, Hóa dược, Tin sinh học và Dược học tính toán.
    • Các nhà nghiên cứu phát triển thuốc điều trị rối loạn chuyển hóa, đái tháo đường và bệnh lý tim mạch.
  • Nội dung đáng tham khảo:
    • Dữ liệu chi tiết về vị trí các hốc phản ứng và các gốc axit amin tương tác tại tâm hoạt động BC và CT của ACC1 (2YL2, 4ASI) và ACC2 (3FF6, 3JRX, 4JRX) ở người.
    • Quy trình kỹ thuật kết hợp AutoDock Vina, UCSF Chimera (mô phỏng MD trong lồng dung môi nước TIP3P ở $37^\circ\text{C}$) và Gaussian 09w (ONIOM QM/MM).
    • Phương pháp xây dựng mô hình QSAR phi tuyến kết hợp giải thuật di truyền (GA), mạng nơron nhân tạo (ANN 4-6-1) và chuyển đổi sang hệ logic mờ (FL/ANFIS) trên môi trường R để biểu diễn mối quan hệ cấu trúc – hoạt tính.

Câu hỏi thường gặp

1. Đồng phân ACC1 và ACC2 ở người có điểm gì khác nhau về cấu trúc và vai trò sinh học trong điều trị MetS?

Theo văn bản luận án, ACC1 có khối lượng phân tử 265 kDa (2346 axit amin), phân bố chủ yếu ở bào tương tế bào gan và mô mỡ, xúc tác tạo malonyl-CoA để sinh tổng hợp axit béo mới (DNL). ACC2 có khối lượng phân tử 280 kDa (2458 axit amin, nhiều hơn ACC1 114 axit amin ở đầu N), khu trú tại màng ty lạp thể mô tim và cơ xương, điều hòa ức chế CPT1 để kìm hãm $\beta$-oxy hóa axit béo. Trong điều trị trúng đích MetS, ức chế ACC1 giúp kìm hãm sinh tổng hợp lipid, còn ức chế ACC2 giúp thúc đẩy quá trình oxy hóa axit béo dư thừa.

2. Luận án sử dụng những mã cấu trúc PDB nào từ Ngân hàng Dữ liệu Protein để khảo sát docking trên ACC1 và ACC2?

Tác giả sử dụng 4 mã cấu trúc PDB từ RCSB PDB:

  • Đối với ACC1: Cấu trúc 2YL2 (chuỗi axit amin 78–617) và 4ASI (chuỗi axit amin 1493–2260).
  • Đối với ACC2: Cấu trúc 3JRX (chuỗi axit amin 217–775), 3FF6 (chuỗi axit amin 1693–2450) cùng mã bổ trợ 4JRX.

3. Quy trình xây dựng mô hình QSAR trong luận án được tiến hành qua các bước cơ bản nào?

Quy trình QSAR trên 130 hợp chất thơm gồm:

  1. Tối ưu hóa cấu trúc phân tử bằng phương pháp PM3 và DFT B3LYP/6-31G(d,p).
  2. Tính toán 17 thông số hóa lượng tử (Hyperchem 8.10) và 16018 thông số mô tả 2D/3D (PaDEL-Descriptor 2.18).
  3. Sàng lọc biến bằng hồi quy tuyến tính đa biến (MLR) với tiêu chuẩn ý nghĩa thống kê $p < 0{,}05$.
  4. Tối ưu tổ hợp biến bằng giải thuật di truyền (GA).
  5. Huấn luyện mạng nơron nhân tạo (ANN) lan truyền ngược 5000 chu kỳ với điều kiện dừng $r_{\text{test}} > 0{,}90$, $\text{SSE} < 9$, $\text{MSE}_{\text{test}} < 0{,}15$.
  6. Thiết lập hệ logic mờ (mô hình ANN 4-6-1), so sánh lựa chọn mô hình và trích xuất khoảng giá trị tối ưu của các biến cấu trúc.

4. Kết quả mô phỏng phức hợp trong môi trường nhiệt độ cơ thể ($37^\circ\text{C}$) ghi nhận năng lượng gắn kết (BE) của các ligand tiêu biểu là bao nhiêu?

Mô phỏng động lực học phân tử trong lồng dung môi nước TIP3P ở $37^\circ\text{C}$ xác định:

  • Trên ACC1: Hợp chất Pfizer6 đạt $\text{BE} = -9{,}30\text{ kcal/mol}$; hợp chất Astra5 đạt $\text{BE} = -6{,}40\text{ kcal/mol}$.
  • Trên ACC2: Hợp chất Pfizer6 đạt $\text{BE} = -10{,}90\text{ kcal/mol}$; hợp chất Astra1 đạt $\text{BE} = -7{,}27\text{ kcal/mol}$.

5. Phương pháp logic mờ (Fuzzy Logic) mang lại lợi điểm gì trong việc phân tích mối quan hệ cấu trúc – hoạt tính (QSAR) của luận án?

Logic mờ cho phép mô hình hóa các mối quan hệ phi tuyến phức tạp giữa tập dữ liệu thông số cấu trúc đầu vào và hoạt tính $\log(IC_{50})$ đầu ra thông qua các hàm liên thuộc và hệ quy tắc mệnh đề điều kiện IF-THEN. Phương pháp này cung cấp biểu đồ bề mặt biểu diễn mức độ và xu hướng ảnh hưởng của từng biến cấu trúc, hỗ trợ xác định các khoảng giá trị tối ưu để định hướng thiết kế các phân tử mới có hoạt tính cao hơn.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Trương Công Minh đã hoàn thành nghiên cứu mô hình ức chế enzyme ACC1 và ACC2 ở người bằng việc kết hợp các phương pháp hóa tin gồm docking phân tử, mô phỏng động lực học phân tử ở $37^\circ\text{C}$, tính toán cơ học lượng tử/cơ học phân tử ONIOM QM/MM và mô hình hóa QSAR phi tuyến. Công trình đã làm sáng tỏ các hốc phản ứng, các tâm gắn kết axit amin then chốt, cũng như năng lượng tương tác của các nhóm phối tử tiêu biểu trên cả ba vùng chức năng BC, BCCP và CT của enzyme. Đồng thời, mô hình lai hóa mạng nơron nhân tạo và logic mờ đã xác lập được 10 thông số cấu trúc quan trọng cùng các khoảng giá trị tối ưu, tạo cơ sở dữ liệu định lượng phục vụ công tác nghiên cứu và định hướng thiết kế các chất ức chế ACC trong điều trị hội chứng chuyển hóa axit béo.