Bối cảnh và vấn đề nghiên cứu

Tăng huyết áp (THA) là một trong những thách thức y tế cộng đồng hàng đầu do tỷ lệ mắc cao và mối liên hệ mật thiết với nhiều biến chứng tim mạch nguy hiểm. Theo số liệu thống kê năm 2015 của Tổ chức Y tế Thế giới (WHO), toàn cầu có khoảng 1,13 tỷ người mắc bệnh, với tỷ lệ mắc ở nam giới (25%) cao hơn nữ giới (20%), và gây ra khoảng 10 triệu ca tử vong mỗi năm (trong đó 4,9 triệu ca do bệnh mạch vành và 3,5 triệu ca do đột quỵ). Tại Việt Nam, tỷ lệ tăng huyết áp dao động từ 18,4% đến 21,1% theo các điều tra dịch tễ học giai đoạn 2005–2018; đáng chú ý, ở nhóm người từ 40–69 tuổi, tỷ lệ người biết mình mắc bệnh chỉ đạt 67,3%, số người được điều trị chiếm 33,2% và chỉ có 12,2% kiểm soát huyết áp hiệu quả.

Hệ renin-angiotensin-aldosteron (RAA) đóng vai trò trung tâm trong cơ chế bệnh sinh của tăng huyết áp. Trong hệ này, enzym chuyển angiotensin (Angiotensin Converting Enzyme - ACE) thực hiện chuyển hóa decapeptide không hoạt động angiotensin I thành octapeptide co mạch mạnh angiotensin II, đồng thời bất hoạt chất giãn mạch bradykinin. Mặc dù các thuốc ức chế ACE kinh điển (captopril, lisinopril, ramipril...) là liệu pháp đầu tay trong điều trị, việc sử dụng lâu dài thường đi kèm các tác dụng phụ như ho khan (chiếm 5–20% bệnh nhân), phù mạch (0,1–0,5%), tăng kali máu và chống chỉ định tuyệt đối ở phụ nữ có thai.

Trước thực trạng đó, việc tìm kiếm các hoạt chất có nguồn gốc tự nhiên nhằm phát triển các tác nhân hạ huyết áp mới là hướng đi cần thiết. Nhóm hợp chất flavonoid từ thực vật với cấu trúc polyphenol đa dạng đã được chứng minh có khả năng ức chế ACE in vitro và in vivo lên tới 50%. Tuy nhiên, cơ chế tương tác ở cấp độ phân tử giữa flavonoid và enzym ACE vẫn chưa được làm sáng tỏ đầy đủ. Đề tài được triển khai nhằm thực hiện hai mục tiêu nghiên cứu cụ thể:

  1. Sàng lọc các hợp chất flavonoid có tác dụng ức chế enzym chuyển angiotensin bằng phương pháp docking phân tử (molecular docking).
  2. Dự đoán tính giống thuốc (drug-likeness), các thông số dược động học và độc tính (ADMET) của các hợp chất tốt nhất thu được sau quá trình sàng lọc.

Đối tượng nghiên cứu của khóa luận là cấu trúc không gian 3 chiều của enzym chuyển angiotensin (mã PDB: 1O86) và tập dữ liệu gồm 116 phân tử flavonoid tự nhiên thu thập từ cơ sở dữ liệu PubChem. Phạm vi nghiên cứu tập trung vào mô phỏng tính toán in silico, đánh giá ái lực liên kết, vị trí tương tác và hồ sơ an toàn sơ bộ trên máy tính.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết và mô hình sinh học

  • Cấu trúc enzym ACE: Nghiên cứu sử dụng cấu trúc tinh thể tia X của ACE tinh hoàn người (tACE) gắn phối tử đồng kết tinh lisinopril ở độ phân giải 2,0 Å (mã PDB: 1O86). Trung tâm hoạt động của enzym chứa mô típ HEXXH gắn ion $\text{Zn}^{2+}$, phối trí với hai gốc histidine (His383, His387) và một glutamate (Glu411), cùng các axit amin quan trọng như Lys511, Tyr520, Gln281, His353 và Ala354.
  • Hợp chất Flavonoid: Cấu trúc khung cơ bản 2-phenyl benzopyrone gồm ba phân nhóm chính: euflavonoid (2-phenylbenzopyran), isoflavonoid (3-benzopyran) và neoflavonoid (4-benzopyran), có khả năng tạo phức chelat với ion kẽm và hình thành các liên kết hydro, liên kết $\pi$-$\pi$ với protein đích.
  • Quy tắc 5 Lipinski (Rule of 5 - RO5): Đánh giá khả năng phát triển thành thuốc đường uống qua 4 tiêu chuẩn: Khối lượng phân tử (MW) < 500 Da; Số nhóm cho liên kết hydro (HBD) < 5; Số nhóm nhận liên kết hydro (HBA) < 10; Hệ số phân bố octanol/nước (LogP) < 5 (chấp nhận vi phạm không quá 1 tiêu chí).
  • Mô hình dự đoán ADMET: Đánh giá các chỉ số hấp thu (tính thấm qua màng Caco-2, hấp thu ở ruột người - HIA), phân bố (thể tích phân bố VDss, thấm qua hàng rào máu não - BBB, thấm vào hệ thần kinh trung ương - CNS), chuyển hóa (khả năng ức chế CYP2D6, CYP3A4), thải trừ (độ thanh thải toàn phần, cơ chất của OCT2) và độc tính (độc tính đột biến gen AMES, ức chế kênh kali tim hERG, độc tính gan, kích ứng da và liều gây chết cấp tính $\text{LD}_{50}$ trên chuột).

Phương pháp nghiên cứu và công cụ tính toán

Quy trình nghiên cứu được tiến hành tuần tự qua 3 bước trên máy tính Dell Inspiron 15 (hệ điều hành Windows 10):

Công cụ / Phần mềm Phiên bản / Nguồn Mục đích sử dụng
Protein Data Bank (PDB) Mã 1O86 Tải cấu trúc 3D của enzym ACE tinh hoàn người
PubChem CSDL trực tuyến Tải cấu trúc 3D và mã SMILES của 116 phối tử flavonoid
Avogadro Gradient liên hợp (Conjugate Gradients) Tối ưu hóa năng lượng phối tử
Discovery Studio Client Bản 2020 Chuẩn bị protein (loại nước, tách lisinopril) và biểu diễn tương tác 2D/3D
AutoDock Tools Bản 1.5.6 Gắn trường lực Kollman, tối ưu hóa hydro phân cực, tạo file pdbqt
AutoDock Vina / AutoDock Bản 4.2 Mô phỏng docking phân tử (Hộp lưới 60x70x50 Å, khoảng cách ô 0,375 Å)
UCSF Chimera Bản 1.14 Tính toán độ lệch bình phương trung bình gốc (RMSD) khi re-dock
pkCSM Công cụ online (CSDL đồ thị phân tử) Dự đoán quy tắc Lipinski và các chỉ số dược động học, độc tính ADMET

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ở lý luận từ các tài liệu trong và ngoài nước về:

  • Dịch tễ và sinh lý bệnh tăng huyết áp: Định nghĩa phân độ huyết áp theo Hội Tim mạch học Việt Nam; cơ chế gây tổn thương cơ quan đích (phì đại thất trái, xơ vữa động mạch não, suy thận mạn, xuất huyết võng mạc).
  • Hệ RAA và cấu trúc ACE: Phân biệt hai dạng đồng dạng sACE (1277 axit amin) và tACE (701 axit amin); cơ chế phân cắt angiotensin I và bất hoạt bradykinin; cấu trúc rãnh hoạt động chứa ion $\text{Zn}^{2+}$ và hai ion $\text{Cl}^-$.
  • Đặc điểm các thuốc ức chế ACE: Phân loại theo nhóm chức (chứa sulfhydryl như captopril; carboxyl như lisinopril, enalapril; phosphinyl như fosinopril) cùng các giới hạn về mặt lâm sàng.
  • Tiềm năng của flavonoid và công nghệ sàng lọc ảo in silico: Phân loại flavonoid; nguyên lý docking phân tử (các thuật toán lấy mẫu Monte Carlo, thuật toán di truyền, động lực học phân tử và hàm tính điểm năng lượng tự do Gibbs $\Delta G$); quy tắc số 5 Lipinski và các thông số ADMET.

Chương 2: Nguyên liệu, thiết bị, phương pháp nghiên cứu

Chương 2 mô tả chi tiết quy trình thực nghiệm trên máy tính:

  • Thiết lập quy trình thẩm định docking thông qua kỹ thuật re-dock lisinopril đồng kết tinh, lấy tiêu chuẩn sai lệch RMSD < 1,5 Å và năng lượng liên kết của lisinopril làm ngưỡng sàng lọc.
  • Thiết lập thông số không gian tìm kiếm: Tọa độ tâm hộp lưới $x = 41,97$ (hoặc 41,62), $z = 46,62$, kích thước 60 $\times$ 70 $\times$ 50 Å bao trọn vị trí ion $\text{Zn}^{2+}$ và các axit amin xúc tác.
  • Chuẩn hóa cấu trúc hóa học các chất bằng Avogadro, chuyển đổi định dạng SMILES và nạp vào máy chủ pkCSM để tính toán thuộc tính hóa lý và sinh học.

Chương 3: Kết quả

Chương 3 trình bày chuỗi kết quả sàng lọc qua từng màng lọc:

  1. Thẩm định quy trình docking: Tái docking lisinopril vào ACE cho giá trị RMSD = 1,5 Å, năng lượng liên kết $\Delta G = -8,4\text{ kcal/mol}$, tạo liên kết hydro với His353, Ala354, Tyr523 và tương tác $\pi$-cation/$\pi$-sigma với Lys511, His383, Glu411 và $\text{Zn}^{2+}$.
  2. Sàng lọc docking 116 chất: Có 86 hợp chất đạt năng lượng liên kết $\Delta G \le -8,4\text{ kcal/mol}$.
  3. Sàng lọc tính giống thuốc (Lipinski RO5): Trong 86 chất, có 33 hợp chất thỏa mãn tối thiểu 3/4 tiêu chí của quy tắc 5 Lipinski.
  4. Dự đoán ADMET: Sàng lọc được 16 hợp chất không độc tính, từ đó xác định được 5 hợp chất tối ưu nhất gồm: Rotenone, Amorphigenin, Quercetin, Apigenin và Hesperitin.
116 Hợp chất Flavonoid (PubChem)
86 Hợp chất liên kết tốt
33 Hợp chất giống thuốc
16 Hợp chất an toàn sơ bộ
05 HỢP CHẤT TIỀM NĂNG NHẤT
(Rotenone, Amorphigenin, Quercetin, Apigenin, Hesperitin)

Chương 4: Bàn luận

Chương 4 tập trung phân tích sâu về dược lý phân tử và phương pháp luận:

  • Phân tích tương quan cấu trúc - tác dụng (SAR): Nhóm hydroxyl liên kết trực tiếp trên vòng thơm và liên kết đôi C2=C3 ở vòng C giúp tăng khả năng tạo liên kết chelat với kẽm và tương tác $\pi$-$\pi$ tại khoang hoạt động.
  • Đánh giá cụ thể từng hoạt chất tiềm năng:
    • Rotenone (ID: 6758): Năng lượng liên kết $-9,2\text{ kcal/mol}$, hấp thu qua ruột 97,05%, Caco-2 đạt 1,415 nhưng có tính thấm qua BBB/CNS và ức chế CYP3A4, cần lưu ý nguy cơ độc tính thần kinh.
    • Amorphigenin (ID: 92207): Năng lượng liên kết $-9,2\text{ kcal/mol}$, hấp thu qua ruột 98,07%, tạo mạng lưới liên kết dày đặc (His353, Tyr520, Gln281, Glu162, His383, Phe457), không gây độc gan/tim/da.
    • Quercetin (ID: 5280343): Năng lượng $-8,5\text{ kcal/mol}$, an toàn tuyệt đối với thần kinh (không qua BBB/CNS), không ức chế CYP; tuy nhiên tính thấm Caco-2 thấp (-0,559).
    • Apigenin (ID: 5280443) & Hesperitin (ID: 72281): Năng lượng liên kết $-8,4\text{ kcal/mol}$, hấp thu qua ruột tốt (>50%), hồ sơ độc tính an toàn nhưng ái lực tương tác tại các axit amin chủ chốt thấp hơn so với 3 chất trên.
  • Bàn luận về ưu - nhược điểm của phương pháp in silico và việc sử dụng phần mềm AutoDock Vina kết hợp công cụ pkCSM.

Kết quả và đóng góp

Tổng hợp số liệu mô phỏng docking và tương tác phân tử

Nghiên cứu đã xác định 5 hợp chất flavonoid có tiềm năng ức chế ACE vượt trội, thể hiện qua các thông số gắn kết và axit amin tương tác tại vị trí hoạt động:

Tên hợp chất PubChem ID Năng lượng liên kết $\Delta G$ (kcal/mol) Các axit amin tại vị trí hoạt động tham gia tạo liên kết
Rotenone 6758 -9,2 Ala354, His353, Tyr520, Gln281, His383, Val380, Asp453, Phe457, Trp279
Amorphigenin 92207 -9,2 Ala354, Glu162, Tyr520, Val380, Leu61, Gln281, Trp279, His383, Asp453, Phe457
Quercetin 5280343 -8,5 Ala354, Tyr520, Gln281, His383, Glu162, Val380, Asp415
Apigenin 5280443 -8,4 Ser222, Met223, Lys118, Phe570, Trp59
Hesperitin 72281 -8,4 His383, Gln281, Ser284, Asp453, Ala354
Lisinopril (Đối chứng) - -8,4 Ala354, His353, His513, Tyr520, Gln281, Lys511, His383, Glu411, Tyr523, Val518, $\text{Zn}^{2+}$

Đánh giá hồ sơ Dược động học và Độc tính (ADMET)

Các chỉ số ADMET dự đoán bằng pkCSM trên 5 hợp chất tiêu biểu:

Thông số đánh giá Rotenone Amorphigenin Quercetin Apigenin Hesperitin
Thấm màng Caco-2 ($\log \text{Papp}$) 1,301 1,011 -0,559 1,301 0,746
Hấp thu ở ruột (HIA) (%) 97,05 98,07 77,20 93,00 82,10
Thấm BBB ($\log \text{BB}$) Qua một phần Qua một phần Không Thấm Không
Thấm CNS ($\log \text{PS}$) Thấm Không Không Thấm Không
Ức chế CYP2D6 / CYP3A4 Không / Có Không / Không Không / Không Không / Không Không / Không
Thanh thải toàn phần ($\log \text{mL/min/kg}$) 0,224 0,226 0,407 0,566 0,044
Cơ chất OCT2 Không Không Không Không Không
Độc tính AMES / Độc tim hERG Không / Không Không / Không Không / Không Không / Không Không / Không
Độc gan / Kích ứng da Không / Không Không / Không Không / Không Không / Không Không / Không
$\text{LD}_{50}$ cấp tính chuột ($\text{mol/kg}$) 2,640 2,663 2,471 2,450 2,042

Đóng góp học thuật của đề tài

  • Làm rõ cơ chế gắn kết ở cấp độ nguyên tử/phân tử của các phân nhóm flavonoid đối với trung tâm xúc tác của ACE, chứng minh vai trò phối trí của các nhóm -OH và vòng thơm đối với các axit amin His383, His353, Tyr520 và Ala354.
  • Cung cấp dữ liệu sàng lọc in silico toàn diện gồm 116 phối tử, thu hẹp danh sách xuống 5 phân tử đích thực sự có triển vọng, giúp tiết kiệm thời gian và kinh phí nghiên cứu thực nghiệm.
  • Phát hiện Amorphigenin là chất có khả năng liên kết mạnh nhất ($\Delta G = -9,2\text{ kcal/mol}$), cấu hình tương tác tương đồng cao với lisinopril và có hồ sơ an toàn thuận lợi, mở ra định hướng nghiên cứu mới vì trước đây chưa có công bố về tác dụng ức chế ACE của chất này.

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

Hạn chế được ghi nhận

  • Toàn bộ kết quả đều dựa trên mô phỏng tính toán in silico; các thuật toán học máy và hàm tính điểm có thể tồn tại sai số nhất định so với điều kiện sinh lý thực tế trong cơ thể người.
  • Nghiên cứu chưa tiến hành kiểm chứng hoạt tính ức chế enzym bằng các phép thử in vitro hoặc thử nghiệm tiền lâm sàng in vivo trên động vật tăng huyết áp.
  • Một số chất có ái lực gắn kết cao nhưng gặp hạn chế về tính thấm màng (như Quercetin có chỉ số Caco-2 thấp) hoặc có nguy cơ đi qua hàng rào máu não (như Rotenone, Apigenin).

Kiến nghị và hướng mở rộng

  1. Tiến hành thử nghiệm thực nghiệm in vitro đánh giá giá trị $\text{IC}_{50}$ ức chế enzym ACE và thử nghiệm in vivo đánh giá khả năng hạ huyết áp của 5 hợp chất (đặc biệt là Amorphigenin và Quercetin).
  2. Nghiên cứu công nghệ bào chế, hệ dẫn thuốc (như nano lipid, micelle) nhằm cải thiện độ tan và tính thấm qua màng ruột Caco-2 cho Quercetin và Hesperitin.
  3. Nghiên cứu sâu hơn về độc tính thần kinh lâu dài của Rotenone và dẫn xuất.
  4. Mở rộng mô hình sàng lọc in silico trên các đích phân tử tim mạch khác trong hệ RAA (như thụ thể $\text{AT}_1$, enzym Renin) và các thư viện hợp chất tự nhiên khác.

Giá trị tham khảo

  • Đối tượng hưởng lợi: Sinh viên, học viên cao học và nghiên cứu viên thuộc các ngành Dược học, Hóa dược, Công nghệ Sinh học, Tin sinh học và Y học cổ truyền.
  • Phần giá trị chuyên môn cao:
    • Quy trình thực nghiệm in silico: Hướng dẫn chi tiết các bước chuẩn bị protein 1O86, tối ưu hóa phối tử, thiết lập thông số hộp lưới docking trên AutoDock/AutoDock Vina và cách thẩm định bằng giá trị RMSD qua re-docking.
    • Dữ liệu tương tác phân tử: Danh mục các axit amin then chốt trong túi liên kết của enzym ACE (His383, His387, Glu411, Lys511, Tyr520, Gln281) là tài liệu tham khảo giá trị cho các nghiên cứu thiết kế thuốc dựa trên cấu trúc (SBDD).
    • Phương pháp luận tích hợp: Cách kết hợp liên hoàn giữa sàng lọc docking, quy tắc 5 Lipinski và công cụ pkCSM để đánh giá toàn diện một phân tử thuốc tiềm năng.

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

1. Tại sao nghiên cứu lại chọn cấu trúc tinh thể ACE có mã PDB 1O86 làm mục tiêu docking?

Mã 1O86 là cấu trúc tinh thể tia X của ACE tinh hoàn người (tACE) gắn với phối tử đồng kết tinh lisinopril ở độ phân giải cao (2,0 Å). Miền tACE có trình tự axit amin và trung tâm hoạt động gần như tương đồng hoàn toàn với miền C của sACE (vùng chịu trách nhiệm chính trong việc chuyển đổi angiotensin I thành angiotensin II và gắn kết chất ức chế), do đó đây là mô hình chuẩn để nghiên cứu thiết kế thuốc ức chế ACE.

2. Tiêu chuẩn nào được sử dụng để thẩm định quy trình docking có đáng tin cậy hay không?

Nghiên cứu sử dụng phương pháp re-dock lisinopril (tách phân tử lisinopril đồng kết tinh ra khỏi protein và dock ngược trở lại vào vị trí hoạt động). Kết quả đo độ lệch bình phương trung bình gốc (RMSD) bằng phần mềm Chimera đạt 1,5 Å (đạt tiêu chuẩn $\le 1,5\text{ Å}$), chứng minh thuật toán và các thông số hộp lưới (grid box) đã thiết lập phản ánh chính xác cấu hình tự nhiên của phức hợp.

3. Trong số các chất sàng lọc, hợp chất nào thể hiện ái lực tương tác tốt nhất với ACE và cấu trúc có đặc điểm gì nổi bật?

Amorphigenin và Rotenone đạt năng lượng liên kết âm nhất ($\Delta G = -9,2\text{ kcal/mol}$), thấp hơn đáng kể so với lisinopril ($-8,4\text{ kcal/mol}$). Trong đó, Amorphigenin thể hiện tương tác tối ưu nhất khi tạo mạng lưới liên kết với các gốc axit amin trọng yếu (His353, Tyr520, Gln281, Glu162) và liên kết $\pi$-$\pi$ rất bền vững với His383 và Phe457 nhờ cấu trúc giàu nhóm hydroxyl chui sâu vào rãnh hoạt động của enzym.

4. Tại sao Quercetin có năng lượng docking tốt ($-8,5\text{ kcal/mol}$) và an toàn nhưng lại cần cải thiện khi phát triển thuốc đường uống?

Mặc dù Quercetin hấp thu tốt ở ruột non (77,20%), không độc tính và không qua hàng rào máu não, nhưng chỉ số thấm qua tế bào biểu mô ruột Caco-2 của Quercetin lại mang giá trị âm (-0,559 $\log \text{Papp}$, tiêu chuẩn cao là $> 0,3$). Điều này cho thấy tính thẩm thấu qua màng tế bào còn hạn chế, đòi hỏi các giải pháp công nghệ bào chế để tăng sinh khả dụng đường uống.

5. Có bao nhiêu hợp chất flavonoid vượt qua toàn bộ các màng lọc từ 116 chất ban đầu?

Từ 116 chất ban đầu trên PubChem: có 86 chất đạt năng lượng liên kết $\Delta G \le -8,4\text{ kcal/mol}$; 33 chất vượt qua quy tắc 5 Lipinski; 16 chất đạt các tiêu chuẩn dược động học và không mang độc tính sơ bộ (AMES, hERG, độc gan, kích ứng da); và cuối cùng 5 hợp chất tối ưu nhất được lựa chọn để phân tích chi tiết là Rotenone, Amorphigenin, Quercetin, Apigenin và Hesperitin.

Kết luận

Khóa luận tốt nghiệp của tác giả Cao Thị Xuân Quỳnh đã xây dựng thành công quy trình sàng lọc ảo in silico dựa trên cấu trúc để đánh giá 116 hợp chất flavonoid ức chế enzym chuyển angiotensin (ACE). Kết quả nghiên cứu đã chọn lọc được 5 hoạt chất tiềm năng nhất (Rotenone, Amorphigenin, Quercetin, Apigenin và Hesperitin) thỏa mãn các tiêu chí về ái lực liên kết cao ($\Delta G$ từ $-8,4$ đến $-9,2\text{ kcal/mol}$), tuân thủ quy tắc 5 Lipinski và sở hữu hồ sơ an toàn ADMET thuận lợi. Đây là nguồn dữ liệu khoa học tin cậy, định hướng cho các nghiên cứu thực nghiệm in vitro, in vivo tiếp theo nhằm phát triển các thuốc mới hoặc sản phẩm hỗ trợ điều trị tăng huyết áp có nguồn gốc tự nhiên.