Tổng quan về luận án

Bệnh Alzheimer (Alzheimer’s Disease – AD) là một hội chứng thoái hóa thần kinh tiến triển không hồi phục, được Alois Alzheimer mô tả lần đầu tiên vào năm 1906 [104]. Đây là nguyên nhân hàng đầu dẫn đến tình trạng sa sút trí tuệ nghiêm trọng ở người cao tuổi, làm suy giảm chức năng nhận thức, biến đổi hành vi và mất khả năng tự chủ trong sinh hoạt hằng ngày [101]. Theo ước tính dịch tễ học toàn cầu, năm 2018 thế giới có khoảng 50 triệu người sống chung với hội chứng sa sút trí tuệ liên quan đến AD, và con số này được dự báo sẽ tăng gấp ba lần, đạt mốc 152 triệu người vào năm 2050 [89]. Gánh nặng kinh tế - xã hội của AD đối với hệ thống y tế toàn cầu là vô cùng to lớn, với chi phí chăm sóc y tế ước tính tăng từ 950 tỷ USD năm 2015 lên tới 2.100 tỷ USD vào năm 2030 [70]. Tại Việt Nam, các thống kê dịch tễ ghi nhận hiện có khoảng 500.000 người trên 60 tuổi mắc sa sút trí tuệ [4].

Vấn đề cốt lõi của y học hiện đại nằm ở chỗ các thuốc điều trị AD đã được phê duyệt lâm sàng bao gồm các chất ức chế acetylcholinesterase (donepezil, galantamin, rivastigmin) và chất đối kháng thụ thể NMDA (memantin) chỉ mang tính điều trị triệu chứng tạm thời, không thể làm chậm hay đảo ngược quá trình thoái hóa neuron thần kinh [73], [101]. Khoảng trống nghiên cứu (research gap) mang tính then chốt hiện nay là sự thất bại của mô hình trị liệu "đơn đích" (single-target drugs) khi đối mặt với bản chất bệnh học đa yếu tố (multi-factorial pathogenesis) của AD [29], [101]. Các bằng chứng khoa học chỉ ra rằng sinh bệnh học AD là mạng lưới phức tạp gồm: sự lắng đọng ngoại bào của mảng $\beta$-amyloid ($A\beta$), sự tích tụ nội bào của các đám rối sợi thần kinh chứa protein Tau bị phosphoryl hóa quá mức [101], [165], sự suy thoái hệ thống dẫn truyền cholinergic [59], phản ứng viêm thần kinh kéo dài [80] và sự tổn thương hệ vi mạch não [62]. Do đó, chiến lược thiết kế phối tử tác động đa đích (Multi-Target-Directed Ligands – MTDL) nhắm đồng thời vào hai enzym trọng yếu: acetylcholinesterase (AChE – chịu trách nhiệm cho khiếm khuyết dẫn truyền cholinergic) và $\beta$-secretase (BACE-1 – enzym cắt mở đầu quá trình sinh tổng hợp mảng $A\beta$) đang trở thành hướng tiếp cận đột phá [97], [169], [176].

Luận án tiến sĩ Dược học của nghiên cứu sinh Trần Thái Sơn, chuyên ngành Hóa dược (Mã số: 62720403), được thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Thái Khắc Minh và PGS.TS. Trần Thành Đạo tại Đại học Y Dược Thành phố Hồ Chí Minh (năm 2022) với đề tài: "Thiết kế, tổng hợp và đánh giá tác động ức chế acetylcholinesterase và beta-amyloid của một số dẫn chất tương đồng curcumin và flavonoid". Luận án giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Liệu việc kết hợp các khung cấu trúc tương đồng curcumin và flavonoid (đặc biệt là chalcon và flavon) có thể tạo ra các phối tử mới có ái lực gắn kết đồng thời trên cả hai đích tác động AChE và BACE-1 hay không?
  • RQ2: Quy trình thiết kế in silico tích hợp (Pharmacophore, 2D-QSAR, Docking phân tử) có khả năng dự đoán chính xác hoạt tính sinh học và định hướng tối ưu hóa tổng hợp hóa học các dẫn chất phenothiazin chalcon và dẫn chất flavon hay không?
  • RQ3: Mức độ tương quan thực nghiệm in vitro giữa khả năng ức chế AChE (phương pháp Ellman) và BACE-1 (phương pháp FRET) của các dẫn chất tổng hợp so với các chất đối chiếu chuẩn (galantamin, quercetin, verubecestat) đạt giá trị như thế nào?
  • H1: Khung cấu trúc lai hóa phenothiazin chalcon có thể tiếp cận đồng thời khoang xúc tác sâu 20 Å và vùng ngoại biên PAS của AChE, đồng thời chèn vào túi kỵ nước S1/S3 của BACE-1.
  • H2: Các mô hình 2D-QSAR xây dựng bằng thuật toán bình phương tối thiểu từng phần (PLS) và các mô hình pharmacophore có độ tin cậy thống kê cao ($R^2 > 0.8$, RMSD $< 1.5\text{ \AA}$) đủ chuẩn xác để sàng lọc ảo thư viện tổ hợp.
  • H3: Các dẫn chất bán tổng hợp từ baicalein, diosmetin và tổng hợp toàn phần phenothiazin chalcon thể hiện hoạt tính ức chế kép ở nồng độ micromolar đến dưới micromolar.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa Thuyết cholinergic (Cholinergic Hypothesis) [59], Thuyết dòng thác Amyloid (Amyloid Cascade Hypothesis) [37], [181], Lý thuyết phối tử đa đích (MTDL) và Lý thuyết thiết kế thuốc với sự hỗ trợ của máy tính (CADD) [172], [173]. Luận án tạo ra bước đột phá về mặt phương pháp luận và thực nghiệm khi xây dựng thành công thư viện tổ hợp ảo, tối ưu hóa các mô hình sàng lọc in silico, tổng hợp hóa học thành công 22 dẫn chất mới (gồm chuỗi phenothiazin chalcon AC1–AC13, dẫn xuất baicalein B1–B3 và dẫn xuất diosmetin D1–D6), đồng thời xác nhận hoạt tính sinh học in vitro đạt hiệu quả ức chế cao trên cả hai enzym mục tiêu.

Literature Review và Positioning

Tổng quan tài liệu cho thấy sự phát triển của các nghiên cứu thuốc điều trị AD tập trung chủ yếu vào hai nhánh enzym bệnh học chính. Đối với đích tác động AChE, enzym xúc tác quá trình thủy phân chất dẫn truyền thần kinh acetylcholin (ACh) thành cholin và acetat tại khe synap [35]. Cấu trúc tinh thể học tia X của AChE (điển hình như mã PDB 1EVE, 1DX6, 1ACJ) đã làm sáng tỏ khoang tác động dạng hẻm sâu khoảng 20 Å và hẹp khoảng 5 Å [182], bao gồm hai vị trí then chốt: (1) Vùng xúc tác (catalytic site) nằm ở đáy hẻm gồm bộ ba xúc tác Ser200–His440–Glu327, túi oxyanion (Gly118, Gly119, Ala201), túi acyl (Trp233, Phe288, Phe290, Phe331) và vùng anion gắn kết cholin (Trp84, Glu199, Phe330); (2) Vùng anion ngoại biên (Peripheral Anionic Site – PAS) nằm ở lối vào hẻm gồm Tyr70, Asp72, Tyr121, Trp279, Tyr334 [182]. Đáng chú ý, vùng PAS không chỉ điều hòa hoạt tính xúc tác mà còn đóng vai trò như một "khuôn mẫu" thúc đẩy quá trình kết tập và hình thành các sợi $A\beta$ độc tính thần kinh [182].

Đối với đích tác động BACE-1 (aspartic protease xuyên màng), enzym này chịu trách nhiệm phân cắt protein tiền sinh chất amyloid (APP) tại vị trí $\beta$ để giải phóng phân mảnh CTF99 và $sAPP\beta$, tạo tiền đề cho phức hợp $\gamma$-secretase phân cắt tiếp tục tạo ra các monomer $A\beta_{1-40}$ và $A\beta_{1-42}$ [37], [181]. Cấu trúc BACE-1 (PDB 1FKN, 5HU1, 6EQM) có vùng gắn kết cơ chất rộng, được che phủ bởi nắp gập linh động dạng kẹp tóc $\beta$ (gọi là vòng "flap" gồm các acid amin ở vị trí 67–75) cùng bộ đôi xúc tác acid-base Asp32/Asp228 [105]. Vị trí hoạt động gồm các túi kỵ nước S1, S3 (Leu30, Phe108, Ile110, Trp115) và các túi thân nước tiếp xúc dung môi S2, S4 (Lys9, Ser10, Thr72, Arg235, Arg307) [105].

Trong kho tàng hợp chất tự nhiên, curcumin (polyphenol chính từ Curcuma longa L. với cấu trúc 1,7-diaryl-1,6-heptadien-3,5-dion) và flavonoid (với khung cơ bản $C_6-C_3-C_6$ gồm hai vòng benzen A và B nối qua nhân pyran C) đã được chứng minh sở hữu hoạt tính sinh học đa dạng: chống oxy hóa, dọn dẹp gốc tự do, tạo phức chelat kim loại, ức chế kết tập $A\beta$, chống viêm thần kinh và ức chế trực tiếp AChE, BACE-1 [67], [117], [127]. Uriarte-Pueyo và Calvo (2011) [162] đã thống kê 128 dẫn chất flavonoid có hoạt tính ức chế AChE tiềm năng, trong đó nổi bật nhất là các phân nhóm flavon và chalcon.

Trong y văn tồn tại hai trường phái tranh luận đối nghịch:

  • Trường phái 1 (Đơn đích ái lực cao): Cho rằng cần tập trung ức chế triệt để enzym BACE-1 bằng các phân tử nhỏ có ái lực cực cao ở dải picomolar đến nanomolar (như chất ức chế verubecestat). Tuy nhiên, các thử nghiệm lâm sàng giai đoạn 3 của các chất này liên tiếp thất bại do xuất hiện độc tính trên hệ thần kinh trung ương và thiếu hiệu quả cải thiện nhận thức toàn diện [98].
  • Trường phái 2 (Đa đích cân bằng - MTDL): Khẳng định tính ưu việt của việc phát triển các phân tử có khả năng ức chế đồng thời nhiều cơ chế bệnh sinh ở mức độ vi phân (micromolar đến dưới micromolar), giúp hiệp đồng tác dụng, hạn chế phản ứng bù trừ của mạng lưới sinh học và giảm thiểu độc tính [97], [101].

Luận án định vị nghiên cứu của mình theo trường phái MTDL, giải quyết trực tiếp hạn chế về tính kém bền hóa học của curcumin (dễ bị phân hủy enol trong kiềm và thoái biến quang hóa) [164] bằng cách lai hóa cấu trúc chalcon với nhân dị vòng phenothiazin cồng kềnh giàu mật độ electron, đồng thời khai thác các dẫn chất tự nhiên baicalein và diosmetin.

Khi so sánh với các công bố quốc tế:

  • Nghiên cứu của Sakata và cộng sự (2020) [137] tổng hợp các dẫn chất chalcon đơn giản ức chế AChE nhưng chưa tối ưu hóa trên BACE-1;
  • Nghiên cứu của Mphahlele và cộng sự (2020) [106] thiết kế dẫn xuất chalcon chứa quinolin tập trung vào độc tính tế bào hơn là cơ chế gắn kết trên cấu trúc tinh thể học;
  • Nghiên cứu của Shimmyo và cộng sự (2008) [144] cùng Han và cộng sự (2019) [60] khảo sát ức chế BACE-1 của flavon tự nhiên nhưng chưa xây dựng được mô hình QSAR định lượng hoàn chỉnh.

Luận án này tiến xa hơn khi thiết lập một quy trình CADD khép kín từ sàng lọc ảo tổ hợp, mô hình hóa pharmacophore, lượng hóa 2D-QSAR đa biến, mô phỏng docking phân tử đa tinh thể đến tổng hợp hóa học và đánh giá in vitro đối chiếu trực tiếp.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp quan trọng vào việc mở rộng và kết nối hai lý thuyết kinh điển trong bệnh học thần kinh: Thuyết cholinergic (được phát triển bởi Bartus và cộng sự, 1982; Francis và cộng sự, 1999) [59] và Thuyết dòng thác Amyloid (được đề xuất bởi Hardy và Higgins, 1992) [37]. Thay vì xem xét hai con đường thoái hóa này một cách độc lập, luận án củng cố mô hình lý thuyết "Tương tác hiệp đồng Cholinergic-Amyloid" thông qua vai trò cầu nối của enzym AChE. Cụ thể, vị trí ngoại biên PAS của AChE tương tác với các monomer $A\beta$ để thúc đẩy quá trình tạo sợi fibril; việc một phân tử thuốc vừa chèn vào khoang xúc tác lẫn PAS của AChE, vừa khóa hoạt tính phân cắt của BACE-1 sẽ cắt đứt đồng thời cả nguồn tạo $A\beta$ ngoại sinh lẫn quá trình kết tập mảng $A\beta$ nội sinh.

Luận án xác lập 2 mệnh đề lý thuyết (Theoretical Propositions) trong hóa dược đa đích:

  • Mệnh đề 1 (P1): Sự tích hợp dị vòng phenothiazin (với tính chất cho electron và độ cồng kềnh lập thể) vào khung nối $\alpha,\beta$-ceton không no của chalcon sẽ tạo ra phân tử có độ dài liên kết tối ưu ($>15\text{ \AA}$), cho phép một đầu tương tác tĩnh điện/$\pi$-$\pi$ với Trp84/Phe330 tại đáy khoang xúc tác AChE, trong khi đầu còn lại tạo liên kết xếp chồng thơm với Trp279/Tyr334 tại vùng PAS.
  • Mệnh đề 2 (P2): Cấu trúc liên hợp thơm mở rộng của dẫn xuất flavon (baicalein, diosmetin) sau khi được ether hóa hoặc ester hóa chọn lọc các nhóm hydroxyl tự do sẽ làm gia tăng tính thân dầu ($\text{Log}P$), tối ưu hóa khả năng lấp đầy túi kỵ nước S1/S3 của enzym BACE-1 mà không làm mất đi liên kết hydro thiết yếu với bộ đôi xúc tác Asp32/Asp228.

Sự chuyển dịch mô thức (Paradigm shift) được minh chứng rõ ràng: từ tư duy thiết kế phối tử cấu trúc cứng nhắc theo mô hình "ổ khóa - chìa khóa" truyền thống của Paul Ehrlich sang mô hình "phối tử đa chức năng linh hoạt" (network-based pharmacophore), có khả năng điều biến đồng thời nhiều thụ thể sinh học trong mạng lưới bệnh sinh phức tạp.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn 4 trụ cột kỹ thuật điện toán và thực nghiệm hóa dược tiên tiến:

[Thư viện tổ hợp ảo (QuaSAR-CombiGen)]
[Sàng lọc Pharmacophore 3D (A1 cho AChE & B1 cho BACE-1)]
[Dự đoán hoạt tính 2D-QSAR (Thuật toán PLS - 184 Descriptors)]
[Mô phỏng Docking phân tử linh hoạt (FlexX trên đa cấu trúc tinh thể PDB)]
[Tổng hợp hóa học (Claisen-Schmidt & Bán tổng hợp Ether/Ester)]
[Kiểm chứng sinh học In Vitro (Ellman Assay & FRET Assay)]

Khung phân tích này tích hợp 3 lý thuyết công cụ:

  1. Lý thuyết Pharmacophore: Mô tả sự sắp xếp không gian 3 chiều của các thuộc tính lập thể và điện tử cần thiết (nhóm cho/nhận liên kết hydro, vòng thơm, vùng kỵ nước) để đảm bảo tương tác siêu phân tử (supramolecular interaction) tối ưu với đích sinh học [77], [126].
  2. Lý thuyết QSAR 2 chiều (2D-QSAR): Toán học hóa mối quan hệ giữa cấu trúc và hoạt tính thông qua 184 thông số mô tả phân tử (phân tử lượng, diện tích bề mặt chia nhỏ, chỉ số hình dạng Kier-Hall, điện tích thành phần) [57], [103].
  3. Lý thuyết nhiệt động học gắn kết phân tử (Molecular Docking): Sử dụng thuật toán tam giác phân (triangle matching) và thuật toán xây dựng lớn dần (incremental construction) kết hợp hàm chấm điểm thực nghiệm (empirical scoring function) để định lượng năng lượng tự do gắn kết $\Delta G$ [90], [91].

Điều kiện biên (Boundary conditions): Mô hình QSAR và Docking áp dụng nghiêm ngặt cho các khung polyphenol có phân tử lượng trong khoảng $250 - 600\text{ Da}$, giá trị $\text{Log}P$ dự đoán từ 1.5 đến 5.5, thỏa mãn các tiêu chuẩn về độ lệch căn quân phương re-docking ($\text{RMSD} \le 1.5\text{ \AA}$ hoặc $\le 2.0\text{ \AA}$) để loại trừ các cấu dạng giả tạo (false positives).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ lập trường triết học thực chứng (Positivism) với bản thể luận khách quan (Objectivism) và nhận thức luận thực nghiệm lượng hóa. Thiết kế nghiên cứu là sự kết hợp chặt chẽ giữa nghiên cứu tính toán lý thuyết (In silico computational modeling) và nghiên cứu thực nghiệm thực nghiệm hóa lý - sinh học (Experimental wet-lab synthesis & bioassay).

Quy mô dữ liệu và tiêu chí chọn mẫu tinh thể học:

  • Hệ thống dữ liệu tinh thể AChE: Lựa chọn 6 cấu trúc tinh thể có độ phân giải cao từ Protein Data Bank (PDB), gồm: 1ACJ, 1DX6, 1EVE, 1W6R, 4EY6 (chuỗi A, B), 4EY7 (chuỗi A, B). Tinh thể 1DX6 (phức hợp với galantamin) và 1EVE (phức hợp với donepezil) được sử dụng làm khuôn mẫu chuẩn để phân tích tương tác vùng xúc tác và vùng PAS.
  • Hệ thống dữ liệu tinh thể BACE-1: Lựa chọn 5 cấu trúc tinh thể chất lượng cao, gồm: 1FKN, 3VEU, 4B78, 5HU0 (chuỗi A, B), 5HU1 (chuỗi A, B) và 6EQM (phức hợp với quercetin). Tinh thể 5HU1 (phức hợp với verubecestat) và 6EQM được chọn làm chuẩn đối chiếu docked poses.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được chuẩn hóa qua các giao thức nghiêm ngặt:

- QuaSAR-CombiGen       - Acepromazin base       - Đo màu Ellman (AChE)
- MOE 2008.10           - Claisen-Schmidt        - FRET huỳnh quang (BACE-1)
- PLS 2D-QSAR           - Ether hóa Williamson   - Xác định IC50 / pIC50
- FlexX Docking         - Định danh 1H/13C-NMR   - Đối chiếu Galantamin/Quercetin

  1. Thiết kế in silico và sàng lọc ảo:

    • Xây dựng thư viện tổ hợp ảo bằng công cụ QuaSAR-CombiGen trong phần mềm MOE 2008.10 bằng cách gắn các nhóm thế đa dạng vào các vị trí "port" được định vị trước trên khung chalcon, baicalein và diosmetin.
    • Xây dựng mô hình pharmacophore dựa trên phối tử bằng công cụ Pharmacophore Elucidation (MOE 2008.10). Mô hình được đánh giá thông qua các chỉ số: hệ số làm giàu EF (Enrichment Factor), độ bao phủ hoạt tính (Active coverage), chất lượng nhận diện GH (Goodness of hits) và đường cong đặc trưng hoạt động máy thu ROC (Receiver Operating Characteristic).
    • Xây dựng mô hình 2D-QSAR: Sử dụng 184 thông số mô tả phân tử 2D thuộc 7 nhóm (Physical Properties, Subdivided Surface Areas, Atom/Bond Counts, Kier-Hall & Kappa Indices, Adjacency/Distance Matrix, Pharmacophore Feature, Partial Charge Descriptors). Áp dụng thuật toán bình phương tối thiểu từng phần (Partial Least Squares – PLS) để xây dựng phương trình hồi quy tuyến tính đa biến dự đoán giá trị hoạt tính sinh học $\text{pIC}{50} = -\log(\text{IC}{50})$.
    • Mô phỏng Docking phân tử: Sử dụng phần mềm thương mại FlexX với thuật toán tam giác phân (triangle matching) và thuật toán xây dựng lớn dần (incremental construction). Giao thức kiểm định tái gắn kết (re-docking) được thực hiện trên tất cả các phối tử đồng kết tinh để đảm bảo giá trị $\text{RMSD} \le 1.5\text{ \AA}$ (đối với AChE) và $\le 2.0\text{ \AA}$ (đối với BACE-1).

  2. Quy trình tổng hợp hóa học:

    • Chuỗi phenothiazin chalcon (AC1–AC13): Xử lý nguyên liệu acepromazin maleat bằng dung dịch kiềm để giải phóng acepromazin base tự do. Thực hiện phản ứng ngưng tụ Claisen–Schmidt giữa acepromazin base và các dẫn chất benzaldehyd đa thế khác nhau trong dung môi ethanol (EtOH) với xúc tác natri ethanolat (EtONa) hoặc dung dịch natri hydroxyd (NaOH 10–60%) ở nhiệt độ phòng hoặc đun hồi lưu.
    • Chuỗi dẫn chất flavon: Bán tổng hợp từ hợp chất tự nhiên baicalein và diosmetin. Thực hiện phản ứng ether hóa Williamson bằng alkyl halide / dimethyl sulfat với xúc tác $\text{K}_2\text{CO}_3$ trong acetone/DMF khan, hoặc phản ứng ester hóa nhóm –OH bằng acid anhydrid/acyl clorid với xúc tác pyridin/DMAP.
    • Kiểm soát chất lượng và xác định cấu trúc: Sản phẩm được tinh chế bằng sắc ký cột cổ điển hoặc kết tinh phân đoạn; kiểm tra độ tinh khiết bằng sắc ký lớp mỏng (TLC). Cấu trúc hóa học được chứng minh bằng các phương pháp phổ học hiện đại: Phổ cộng hưởng từ hạt nhân proton ($^1\text{H-NMR}$), phổ cộng hưởng từ hạt nhân carbon-13 ($^{13}\text{C-NMR}$), phổ hồng ngoại (IR) và phổ khối lượng (MS).

  3. Giao thức thử nghiệm hoạt tính sinh học in vitro:

    • Thử nghiệm ức chế AChE: Tiến hành theo phương pháp đo màu quang phổ của Ellman [48]. Enzym AChE (nguồn gốc từ Electrophorus electricus) xúc tác thủy phân cơ chất acetylthiocholin iodid (ATCI) tạo thành thiocholin; thiocholin phản ứng ngay với thuốc thử Ellman (acid 5,5'-dithio-bis-nitrobenzoic – DTNB) tạo ra anion 5-thio-2-nitrobenzoat có màu vàng đặc trưng. Mật độ quang (OD) được đo tự động trên đầu đọc vi đĩa ở bước sóng $\lambda = 405\text{ nm}$. Chất đối chiếu chuẩn là galantamin.
    • Thử nghiệm ức chế BACE-1: Tiến hành bằng phương pháp truyền năng lượng cộng hưởng huỳnh quang (Fluorescence Resonance Energy Transfer – FRET) [42]. Cơ chất peptid đặc hiệu mang một đầu huỳnh quang và một đầu dập tắt huỳnh quang (quencher). Khi enzym BACE-1 tái tổ hợp phân cắt cơ chất, sự truyền năng lượng nội phân tử bị cắt đứt, phát ra tín hiệu huỳnh quang mạnh. Tín hiệu được ghi nhận trên máy đo huỳnh quang ở bước sóng kích thích $\lambda_\text{ex} = 320\text{ nm}$ và bước sóng phát xạ $\lambda_\text{em} = 405\text{ nm}$. Chất đối chiếu chuẩn là quercetin và verubecestat.

Data và phân tích

Dữ liệu thực nghiệm được xử lý bằng các công cụ thống kê chuyên sâu:

  • Phân tích hồi quy PLS trong QSAR được tối ưu hóa qua phương pháp kiểm định chéo rút từng phần tử (Leave-One-Out – LOO) với hệ số tương quan chéo $q^2_\text{LOO} > 0.6$, hệ số xác định hiệu chỉnh $R^2 > 0.8$, sai số căn quân phương (Root-Mean-Square Error – RMSE) $< 0.35$ đơn vị $\text{pIC}_{50}$.
  • Giá trị $\text{IC}_{50}$ (nồng độ ức chế 50% hoạt tính enzym) được tính toán từ đường cong ức chế liều - đáp ứng (dose-response curve) phi tuyến tính thông qua phần mềm GraphPad Prism, với độ tin cậy $95%$ Confidence Interval ($p < 0.05$).
  • Độ tin cậy phương pháp Docking được xác nhận khi giá trị RMSD giữa cấu dạng docked pose và phối tử đồng kết tinh thực nghiệm trong tinh thể 1DX6 (galantamin) đạt $0.68\text{ \AA}$, và trong tinh thể 5HU1 (verubecestat) đạt $1.12\text{ \AA}$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Thiết kế và sàng lọc thành công thư viện tổ hợp ảo: Từ hàng nghìn phân tử ảo tạo bởi QuaSAR-CombiGen, mô hình pharmacophore A1 (cho AChE) và B1 (cho BACE-1) đã sàng lọc được nhóm chất tối ưu; mô hình 2D-QSAR PLS dự đoán chính xác dải hoạt tính $\text{pIC}_{50}$ từ $4.5$ đến $7.8$.
  2. Tổng hợp thành công 22 dẫn chất chọn lọc: Gồm 13 dẫn chất phenothiazin chalcon mới (AC1–AC13) và 9 dẫn chất flavon (B1–B3, D1–D6) với hiệu suất phản ứng Claisen-Schmidt đạt $45 - 82%$, cấu trúc được khẳng định tuyệt đối bằng dữ liệu $^1\text{H-NMR}$ và $^{13}\text{C-NMR}$ (như phổ phân tích tín hiệu nhóm olefin $trans$ $J = 15.5 - 16.0\text{ Hz}$, tín hiệu đặc trưng của khung phenothiazin tại $\delta 6.8 - 7.6\text{ ppm}$, nhóm carbonyl tại $\delta 185 - 190\text{ ppm}$).
  3. Phát hiện các phân tử ức chế kép vượt trội:
    • Dẫn chất phenothiazin chalcon AC12 và AC4 thể hiện hoạt tính ức chế AChE mạnh với $\text{IC}_{50}$ ở mức micromolar thấp, tương đương với hoạt lực của chất đối chiếu galantamin.
    • Các dẫn chất này đồng thời ức chế enzym BACE-1 với giá trị $\text{IC}_{50}$ vượt trội hơn so với chất đối chiếu tự nhiên quercetin trong thử nghiệm FRET.
    • Dẫn chất flavon D5 và D6 (dẫn xuất alkyl hóa từ diosmetin) thể hiện hoạt lực ức chế chọn lọc cao trên BACE-1.
  4. Làm sáng tỏ cơ chế gắn kết siêu phân tử ở cấp độ nguyên tử:
    • Trên AChE (tinh thể 1DX6): Dẫn chất AC12 định vị hoàn hảo dọc theo hẻm xúc tác; nhân phenothiazin cồng kềnh tạo tương tác xếp chồng thơm $\pi$-$\pi$ và tương tác cation-$\pi$ mạnh mẽ với các gốc thơm Trp84 và Phe330 tại vùng đáy anion, trong khi nhóm liên kết chalcon mở rộng ra miệng ống để tạo liên kết kỵ nước và van der Waals với Trp279 và Tyr334 tại vùng ngoại biên PAS.
    • Trên BACE-1 (tinh thể 5HU1 và 6EQM): Dẫn chất AC4 và AC12 chèn sâu vào khe liên kết cơ chất, tạo liên kết hydro bền vững với nhóm carboxyl của cặp đôi xúc tác Asp32 và Asp228, đồng thời vòng thế thơm lấp đầy túi kỵ nước S1/S3, khóa chặt nắp "flap" ở cấu dạng đóng.
  5. Hiện tượng phản trực giác (Counter-intuitive Finding): Trái với quan niệm thông thường cho rằng khung phenothiazin quá cồng kềnh sẽ bị cản trở lập thể không thể chui vào hẻm hẹp 5 Å của AChE, kết quả mô phỏng docking và thử nghiệm in vitro chứng minh rằng nhờ tính linh động góc quay của cầu nối carbonyl-olefin, phân tử tự điều chỉnh cấu dạng uốn lượn (butterfly conformation), giúp gia tăng diện tích tiếp xúc bề mặt và tạo ái lực liên kết cao hơn đáng kể so với các khung chalcon phẳng đơn giản.

Implications đa chiều

  • Về mặt lý thuyết: Khẳng định tính đúng đắn của chiến lược MTDL, cung cấp bằng chứng thực nghiệm rõ nét rằng việc lai hóa cấu trúc tự nhiên (flavonoid/chalcon) với dị vòng tổng hợp (phenothiazin) là một hướng đi giàu tiềm năng để tạo ra các phối tử tác động đa cơ chế cho bệnh Alzheimer.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa "4 pha tích hợp" (In silico Combinatorial Library $\rightarrow$ Pharmacophore/QSAR Filtering $\rightarrow$ Rational Chemical Synthesis $\rightarrow$ Standardized Bioassay) có khả năng chuyển giao và áp dụng rộng rãi cho các dự án phát minh thuốc chống ung thư, kháng viêm hoặc các bệnh thoái hóa thần kinh khác tại các viện nghiên cứu và trường đại học.
  • Về mặt ứng dụng thực tiễn: Dẫn chất AC12, AC4 và D5 là các hợp chất chìa khóa (lead compounds) đầy triển vọng, có giá trị đăng ký sáng chế và tiếp tục các thử nghiệm tiền lâm sàng sâu hơn.
  • Về khả năng khái quát hóa: Khung phân tích QSAR-PLS và mô hình docking đa tinh thể hoàn toàn có thể mở rộng để sàng lọc các thư viện hợp chất tự nhiên phong phú của hệ thực vật Việt Nam nhằm tìm kiếm các hoạt chất ức chế các mục tiêu sinh học khác.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn khoa học nội tại:

  1. Giới hạn về mô phỏng động lực học phân tử (MD Simulation): Nghiên cứu dừng lại ở kỹ thuật docking phân tử tĩnh (FlexX), chưa thực hiện mô phỏng động lực học phân tử (Molecular Dynamics) trong môi trường dung môi nước ở thang thời gian $100 - 500\text{ ns}$ để đánh giá độ ổn định nhiệt động học, năng lượng tự do liên kết theo phương pháp MM-PBSA/MM-GBSA và sự thay đổi cấu dạng linh hoạt của vòng "flap" BACE-1 theo thời gian.
  2. Giới hạn về mô hình thử nghiệm sinh học: Đánh giá hoạt tính mới dừng lại ở mức độ enzym in vitro (cell-free assay), chưa tiến hành thử nghiệm trên các dòng tế bào thần kinh sống (như dòng tế bào u nguyên bào thần kinh SH-SY5Y hoặc tế bào PC12) để đánh giá khả năng bảo vệ neuron chống lại độc tính gây bởi oligomer $A\beta_{1-42}$.
  3. Giới hạn về dữ liệu dược động học và hàng rào máu não: Chưa thực hiện các thử nghiệm thấm màng nhân tạo song song (PAMPA-BBB assay) để kiểm chứng thực nghiệm khả năng vượt qua hàng rào máu não của các dẫn chất phenothiazin chalcon có phân tử lượng tương đối lớn.
  4. Quy mô tổng hợp hóa học: Số lượng dẫn chất tổng hợp được kiểm chứng in vitro (22 chất) còn ở mức khiêm tốn do hạn chế về thời gian và nguồn hóa chất tinh khiết trong điều kiện phòng thí nghiệm.

Chương trình nghiên cứu tương lai (Future Research Agenda) bao gồm 5 định hướng:

  • Hướng 1: Ứng dụng mô phỏng MD chuyên sâu ($>200\text{ ns}$) kết hợp phân tích tương tác lượng tử QM/MM để khảo sát cơ chế proton hóa của bộ đôi Asp32/Asp228 khi tương tác với dẫn chất phenothiazin.
  • Hướng 2: Thực hiện thử nghiệm in vitro trên tế bào SH-SY5Y nhằm đánh giá hoạt tính chống chết theo chương trình (anti-apoptosis), giảm stress oxy hóa nội bào và ức chế tạo mảng $A\beta$ nội sinh.
  • Hướng 3: Thử nghiệm mô hình màng PAMPA-BBB và đánh giá độc tính cấp, độc tính bán trường diễn trên chuột thực nghiệm.
  • Hướng 4: Tối ưu hóa cấu trúc dẫn xuất AC12 thông qua gắn thêm các nhóm thế ưa nước (như morpholin, piperazin) để cải thiện độ tan và nâng cao sinh khả dụng đường uống.
  • Hướng 5: Đánh giá hiệu quả cải thiện trí nhớ và khả năng học tập in vivo trên mô hình chuột chuyển gen AD (APP/PS1 transgenic mice).

Tác động và ảnh hưởng

  • Tác động học thuật: Luận án là một trong những công trình tiên phong tại Việt Nam áp dụng toàn diện và chuyên sâu các kỹ thuật CADD hiện đại trong phát minh thuốc điều trị Alzheimer. Kết quả nghiên cứu đã được công bố trên các tạp chí khoa học chuyên ngành Hóa dược uy tín, đóng góp nguồn dữ liệu cấu trúc và hoạt tính sinh học chuẩn xác cho cộng đồng khoa học quốc tế, với tiềm năng trích dẫn cao trong các nghiên cứu phát triển phối tử đa đích (MTDL).
  • Chuyển đổi công nghiệp Dược: Cung cấp quy trình tổng hợp hóa học khả thi, sử dụng nguồn nguyên liệu đầu vào sẵn có (acepromazin maleat, baicalein, diosmetin) với hiệu suất phản ứng cao, mở ra cơ hội cho các công ty dược phẩm trong nước tiếp cận hướng nghiên cứu phát triển hoạt chất mới điều trị bệnh thần kinh thay vì chỉ sản xuất thuốc generic.
  • Ý nghĩa chính sách y tế và xã hội: Trong bối cảnh tốc độ già hóa dân số tại Việt Nam đang diễn ra nhanh chóng, việc chủ động nghiên cứu và phát triển các tác nhân trị liệu AD thế hệ mới từ nguồn tài nguyên nội lực sẽ góp phần giảm bớt gánh nặng ngân sách y tế quốc gia, hạ giá thành điều trị và nâng cao chất lượng cuộc sống cho hàng triệu bệnh nhân cao tuổi và gia đình.
  • Tầm vóc quốc tế: Nghiên cứu khẳng định năng lực hội nhập khoa học công nghệ đỉnh cao của các nhà khoa học Dược học Việt Nam trong lĩnh vực thiết kế thuốc hỗ trợ bởi máy tính, sánh ngang với các nghiên cứu tiên tiến trên thế giới.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học thuật Hóa dược: Tiếp cận một tài liệu mẫu mực về phương pháp luận kết hợp in silico (Pharmacophore, 2D-QSAR PLS, FlexX Docking) với tổng hợp hữu cơ và thử nghiệm sinh học chuẩn quốc tế.
  • Bộ phận Nghiên cứu & Phát triển (R&D) tại các Tập đoàn Dược phẩm: Tiếp nhận danh mục các cấu trúc phân tử tiềm năng (đặc biệt là chuỗi dẫn chất phenothiazin chalcon AC12, AC4) đã được tối ưu hóa để làm chất khởi nguồn (lead optimization) cho các chương trình phát triển thuốc tiền lâm sàng.
  • Các nhà Hóa học Hợp chất Tự nhiên: Khai thác quy trình bán tổng hợp gia tăng giá trị sinh học từ các flavonoid tự nhiên phổ biến như baicalein (từ Scutellaria baicalensis) và diosmetin (từ Chrysanthemum morifolium).
  • Bác sĩ lâm sàng và Bệnh nhân Alzheimer: Đón nhận triển vọng về các liệu pháp điều trị đa cơ chế trong tương lai, giúp giải quyết tận gốc căn nguyên bệnh học thay vì chỉ trì hoãn triệu chứng suy giảm nhận thức.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc phát triển và chứng minh trên thực nghiệm mô hình "Khung lai hóa phân tử linh hoạt" (Hybrid Flexible Scaffold). Luận án đã mở rộng Lý thuyết phối tử đa đích (MTDL) bằng cách tích hợp thành công cấu trúc dị vòng phenothiazin giàu mật độ electron vào mạch nối chalcon liên hợp. Cấu trúc này phá vỡ rào cản cản trở lập thể truyền thống, tạo ra phân tử có khả năng gắn kết lưỡng vùng (dual-site binding) đồng thời tại vùng đáy xúc tác (Trp84/Phe330) và vùng miệng ngoại biên PAS (Trp279/Tyr334) của AChE, song song với việc ức chế vị trí xúc tác Asp32/Asp228 của BACE-1.

2. Sự đổi mới về phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế trước đây?

So với các nghiên cứu của Sakata và cộng sự (2020) [137] (chỉ dừng ở tổng hợp chalcon đơn đích trên AChE) hay Shimmyo và cộng sự (2008) [144] (chỉ khảo sát in vitro flavonoid đơn thuần), luận án đã thiết lập quy trình tích hợp 4 lớp in silico liên hoàn:

  1. Xây dựng thư viện tổ hợp ảo bằng QuaSAR-CombiGen;
  2. Sàng lọc cấu trúc qua mô hình Pharmacophore Elucidation 3D;
  3. Dự đoán lượng hóa hoạt tính bằng mô hình 2D-QSAR PLS với 184 thông số mô tả phân tử;
  4. Khảo sát ái lực gắn kết trên hệ thống đa cấu trúc tinh thể PDB bằng phần mềm FlexX trước khi tiến hành thực nghiệm tổng hợp và kiểm chứng qua 2 kỹ thuật đo màu Ellman và FRET. Quy trình này giúp tiết kiệm $80%$ thời gian và chi phí sàng lọc thực nghiệm.

3. Phát hiện bất ngờ nhất trong quá trình nghiên cứu là gì và bằng chứng thực nghiệm chứng minh?

Phát hiện bất ngờ nhất là hoạt tính ức chế AChE cực mạnh của dẫn chất phenothiazin chalcon AC12 và AC4. Ban đầu, các dự đoán lý thuyết cổ điển e ngại thể tích phân tử lớn của nhân phenothiazin sẽ ngăn cản việc thâm nhập vào khe hẹp 5 Å của AChE. Tuy nhiên, dữ liệu phổ $^1\text{H-NMR}$/ $^{13}\text{C-NMR}$ kết hợp docked poses trên tinh thể 1DX6 và thực nghiệm đo màu Ellman đã chứng minh: liên kết olefin-carbonyl có khả năng quay tạo cấu dạng "cánh bướm" gập góc, giúp nhân phenothiazin chui lọt xuống đáy khoang để tạo tương tác cation-$\pi$ mạnh với Trp84, trong khi phần còn lại của phân tử khóa chặt vùng ngoại biên PAS, đạt hiệu lực $\text{IC}_{50}$ ở mức micromolar tương đương galantamin.

4. Giao thức nhân rộng (Replication Protocol) có được mô tả đầy đủ không?

Luận án cung cấp chi tiết toàn bộ thông số thực nghiệm để tái lập:

  • In silico: Thông số cấu hình phần mềm MOE 2008.10, danh mục 184 biến mô tả, thuật toán PLS, mã PDB tinh thể (1DX6, 1EVE, 5HU1, 6EQM), thông số FlexX.
  • Tổng hợp: Nồng độ chất phản ứng, tỷ lệ mol giữa acepromazin base và benzaldehyd, loại dung môi (EtOH khan), chất xúc tác (EtONa, NaOH), nhiệt độ phản ứng ($25 - 80^\circ\text{C}$), quy trình sắc ký tinh chế và dữ liệu phổ định danh ($^1\text{H-NMR}$, $^{13}\text{C-NMR}$, MS, IR).
  • In vitro: Nồng độ cơ chất ATCI, thuốc thử DTNB, đệm phosphat pH 8.0, bước sóng đo quang 405 nm cho thử nghiệm Ellman; trình tự peptid FRET, bước sóng kích thích 320 nm, phát xạ 405 nm cho thử nghiệm BACE-1.

5. Lộ trình nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào?

  • Năm 1–2: Mô phỏng động lực học phân tử MD 500 ns, tính toán năng lượng liên kết MM-GBSA và tổng hợp mở rộng 50 dẫn xuất thế hệ 2 nhằm tối ưu hóa độ tan.
  • Năm 3–4: Thử nghiệm độc tính tế bào, khả năng bảo vệ neuron trên dòng tế bào SH-SY5Y và đo độ thấm màng PAMPA-BBB.
  • Năm 5–6: Đánh giá dược động học (ADMET), xác định sinh khả dụng và thử nghiệm hiệu quả in vivo trên chuột cống trắng gây mô hình sa sút trí tuệ bằng scopolamin.
  • Năm 7–8: Đánh giá tiền lâm sàng toàn diện trên chuột biến đổi gen mang đột biến APP/PS1; hoàn thiện hồ sơ đăng ký sáng chế quốc tế.
  • Năm 9–10: Tiến hành thử nghiệm lâm sàng Phase I trên người tình nguyện khỏe mạnh để đánh giá độ an toàn và dung nạp của ứng viên thuốc dẫn đầu.

Kết luận

  1. Xây dựng thành công hệ thống công cụ In Silico hoàn chỉnh: Phát triển thư viện tổ hợp ảo, mô hình pharmacophore 3D và phương trình hồi quy 2D-QSAR (sử dụng thuật toán bình phương tối thiểu từng phần PLS với 184 thông số mô tả) có độ tin cậy thống kê cao ($R^2 > 0.8$, $q^2_\text{LOO} > 0.6$), dự đoán chuẩn xác hoạt tính ức chế AChE và BACE-1.
  2. Tổng hợp và xác định cấu trúc 22 dẫn chất mới: Tổng hợp thành công 13 dẫn chất phenothiazin chalcon (AC1–AC13) bằng phản ứng ngưng tụ Claisen–Schmidt từ acepromazin base và 9 dẫn chất flavon (B1–B3, D1–D6) thông qua bán tổng hợp ether hóa/ester hóa từ baicalein và diosmetin; cấu trúc hóa học được chứng minh sáng tỏ qua phổ $^1\text{H-NMR}$, $^{13}\text{C-NMR}$, MS và IR.
  3. Chứng minh hoạt tính ức chế kép đột phá In Vitro: Xác nhận thực nghiệm bằng phương pháp Ellman và FRET cho thấy các dẫn chất AC12, AC4 sở hữu hoạt tính ức chế mạnh mẽ trên cả AChE và BACE-1 ở mức micromolar, mở ra nhóm hoạt chất MTDL tiềm năng vượt trội so với các chất đối chiếu chuẩn galantamin và quercetin.
  4. Làm sáng tỏ cơ chế tương tác ở cấp độ phân tử: Dữ liệu Docking trên phần mềm FlexX với đa cấu trúc tinh thể PDB (1DX6, 5HU1, 6EQM) đã giải mã thành công cơ chế gắn kết lưỡng vùng của phenothiazin chalcon tại đáy xúc tác và lối vào PAS của AChE, cùng tương tác ức chế bộ đôi xúc tác Asp32/Asp228 của BACE-1.
  5. Thúc đẩy bước chuyển dịch mô thức nghiên cứu Hóa dược: Đóng góp bằng chứng thực nghiệm vững chắc củng cố lý thuyết phối tử đa đích (MTDL), xóa bỏ giới hạn của các phương pháp đơn đích truyền thống trong điều trị bệnh thoái hóa thần kinh.
  6. Mở ra 3 dòng nghiên cứu khoa học mới:
    • Hướng nghiên cứu tối ưu hóa cấu trúc dẫn xuất phenothiazin chalcon hướng thần kinh;
    • Hướng nghiên cứu bán tổng hợp nâng cao giá trị dược lý từ nguồn flavonoid tự nhiên Việt Nam;
    • Hướng ứng dụng quy trình sàng lọc ảo CADD tích hợp trong phát minh thuốc điều trị các bệnh lý phức tạp. Công trình khẳng định tầm vóc học thuật xuất sắc và để lại giá trị khoa học bền vững cho ngành Dược học Việt Nam và quốc tế.