Tổng quan luận án

Bệnh Alzheimer (AD) là dạng bệnh thoái hóa thần kinh phổ biến nhất ở người lớn tuổi, chiếm khoảng 60–80% tổng số ca mất trí nhớ. Bệnh đặc trưng bởi sự suy giảm dần trí nhớ, nhận thức, đi kèm sự xuất hiện của các mảng lão suy ngoại bào hình thành từ peptide Amyloid beta (Aβ) và các đám rối sợi thần kinh nội bào từ protein tau siêu phosphoryl hóa. Giả thuyết tầng amyloid (amyloid cascade) chỉ ra rằng sự tích tụ ngoại bào quá mức của peptide Aβ là sự kiện khởi phát chuỗi bệnh lý dẫn đến phá hủy tế bào thần kinh và khớp thần kinh (synapse).

Peptide Aβ được tạo ra từ quá trình phân giải protein tiền chất amyloid (APP) bởi các enzyme β-secretase và γ-secretase, trong đó hai biến thể chính là Aβ1-40 và Aβ1-42. Mặc dù nồng độ Aβ1-40 trong cơ thể cao gấp khoảng 9 lần Aβ1-42, Aβ1-42 lại là thành phần chủ yếu trong các mảng bám thần kinh do có tốc độ kết tụ nhanh hơn và độc tính cao hơn đối với tế bào thần kinh.

Khoảng trống nghiên cứu then chốt xuất phát từ sự thay đổi trong quan điểm bệnh học: các nghiên cứu gần đây chứng minh rằng các dạng kết tụ trung gian hòa tan như oligomer bậc thấp (đặc biệt là trimer, tetramer, hexamer, dodecamer của Aβ1-42) có độc tính tế bào thần kinh cao hơn hẳn so với sợi hoàn chỉnh (mature fibrils). Tuy nhiên, do đặc tính cấu trúc mất trật tự nội tại (intrinsically disordered), tốc độ kết hợp cao và thời gian sống ngắn của các oligomer, các phương pháp phân tích thực nghiệm như tinh thể học tia X hay cộng hưởng từ hạt nhân rất khó xác định cấu trúc ba chiều chi tiết của chúng. Do đó, phương pháp mô phỏng động lực học phân tử (MD) là công cụ cần thiết để tiếp cận cấu trúc và động lực học của các oligomer này.

Luận án tiến sĩ chuyên ngành Vật lý kỹ thuật (mã số: 9520401) của tác giả Nguyễn Hoàng Linh, thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh năm 2022, dưới sự hướng dẫn của GS.TSKH. Mai Xuân Lý và PGS.TS. Huỳnh Quang Linh, tập trung giải quyết các mục tiêu nghiên cứu cụ thể được đánh số như sau:

  1. Xác định ảnh hưởng của cấu trúc ban đầu lên mô phỏng oligomer bậc thấp như trimer Aβ1-42.
  2. Xác định các tính chất hoá lý của trimer, tetramer Aβ1-42 bằng mô phỏng động lực học phân tử.
  3. Đánh giá và so sánh ảnh hưởng của các dạng kết tụ dodecamer Aβ1-42, sợi hoàn chỉnh Aβ1-42 lên màng tế bào thần kinh.

Đối tượng nghiên cứu của luận án là peptide Aβ1-42 ở các trạng thái tập hợp: trimer, tetramer, dodecamer và sợi hoàn chỉnh Aβ1-42 (chuỗi 42 amino acid đầy đủ độ dài), cùng mô hình màng tế bào thần kinh đa thành phần ở trạng thái bệnh lý (DPPC, POPC, SM, Cholesterol, GM1). Phạm vi nghiên cứu tập trung vào không gian cấu hình phân tử và các đặc tính động lực học ở thang thời gian từ nano giây đến micro giây thông qua mô hình tính toán số.

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

Luận án điểm lại và kế thừa hệ thống các công trình nghiên cứu lý thuyết và thực nghiệm trong và ngoài nước:

  • Về độc tính của các dạng oligomer Aβ: Jana và cộng sự (2016), Ono và cộng sự (2017) chỉ ra tetramer Aβ có thể là dạng oligomer mang độc tính cao nhất đối với tế bào thần kinh. Drew và cộng sự (2020) phát hiện tetramer hoặc oligomer bậc lớn hơn có khả năng tương tác vượt qua màng tế bào thần kinh, mở đường cho ion Ca2+ thâm nhập. Economou và cộng sự (2007) xác định hexamer và dodecamer Aβ1-42 là hạt nhân chính của quá trình sinh sợi. Bernstein và cộng sự (2009) chỉ ra dodecamer Aβ1-42 là dạng oligomer bền vững, hòa tan, mang độc tính thần kinh cơ bản.
  • Về cấu trúc sợi hoàn chỉnh và mô hình amyloid: Các cấu trúc sợi Aβ hoàn chỉnh từ thực nghiệm được xác định bởi Lührs (2005), Petkova (2002), Meinhardt (2009), Paravastu (2008), Colvin (2016). Các công trình này chỉ ra cấu trúc "cross beta", trong đó các dải β định hướng vuông góc với trục sợi, có sự tham gia của cầu muối Asp23-Lys28. Cấu trúc sợi thực nghiệm PDB ID 2NAO/2NAQ được sử dụng làm hệ quy chiếu cấu trúc trong luận án.
  • Về cấu trúc dạng thùng β (β-barrel): Eisenberg và cộng sự (2012) chứng minh các đoạn Aβ dài 11 residue (24–31, 25–35, 26–36) có thể tạo thành cấu trúc thùng β hình trụ rỗng. Xi và cộng sự (2019) xây dựng mô hình thùng β nhân tạo chứa vùng β2-turn-β3 (residue 27–42) trong trimer và tetramer Aβ1-42.
  • Về mô phỏng tương tác màng lipid: Brown và cộng sự (2018) mô phỏng tetramer Aβ1-42 và tương tác với màng nhưng chỉ khảo sát tại một nhiệt độ cố định. Drolle và cộng sự (2017) đề xuất mô hình màng đa thành phần phức hợp bệnh lý Disease 1 (D1) gồm 5 loại lipid, phản ánh sát thực tế màng tế bào thần kinh người bệnh AD hơn các mô hình lipid đơn thành phần (như POPC hoặc DOPC đơn thuần).

Vị trí của luận án: Luận án chọn giải quyết khoảng trống về việc nghiên cứu cấu trúc oligomer Aβ1-42 bậc thấp không bị thiên lệch bởi cấu hình ban đầu; khảo sát sự hình thành tự phát của motif thùng β mà không cần áp đặt cấu trúc dựng sẵn; và so sánh cơ chế phá hủy màng tế bào thần kinh giữa oligomer (dodecamer) và sợi hoàn chỉnh trên mô hình màng lipid phức hợp D1.

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

Cơ sở lý thuyết về cấu trúc protein

Luận án áp dụng các nguyên lý hóa lý và sinh học cấu trúc của protein:

  • Cấu trúc bậc một của Aβ1-42 gồm 42 amino acid: 1DAEFR5 6HDSGY10 11EVHHQ15 16KLVFF20 21AEDVG25 26SNKGA30 31IIGLM35 36VGGVV40 41IA, chia làm ba vùng: đầu N (1–16), vùng trung tâm kị nước (17–29) và đuôi C (30–42).
  • Các góc nhị diện nhánh chính: phi ($\phi$, hình thành từ C-N-C$\alpha$-C), psi ($\psi$, hình thành từ N-C$\alpha$-C-N), và omega ($\omega$, hình thành từ C$\alpha$-C-N-C$\alpha$).
  • Cấu trúc bậc hai: xoắn $\alpha$-helix, dải $\beta$ (song song và đối song), turn (phân loại $\alpha, \beta, \gamma, \delta, \pi$), loop, beta-hairpin.
  • Các liên kết ổn định: liên kết hydro (khoảng cách donor-acceptor < 3.5 Å, góc donor-H-acceptor > 135°), cầu muối (khoảng cách tâm nhánh phụ của Asp/Glu và Arg/Lys/His mang điện tích trái dấu < 4 Å), tương tác tĩnh điện, tương tác van der Waals (vdW), tương tác kị nước (sử dụng thang đo hydropathy của Kyte và Doolittle).

Phương pháp nghiên cứu và công cụ mô phỏng

Luận án kết hợp đa phương pháp mô phỏng tính toán định lượng:

Hạng mục phương pháp Công cụ / Thuật toán / Mô hình sử dụng Mục đích áp dụng trong luận án
Mô hình hạt thô (Coarse-grained) Trường lực UNRES (UNited RESidue), mô hình hóa chuỗi C$\alpha$ cùng các hạt đại diện cho nhánh phụ (SC) và nhóm peptide (p). Lấy mẫu không gian pha cấu hình rộng lớn cho tetramer Aβ1-42 với thang thời gian hiệu dụng lớn hơn 3–4 bậc độ lớn.
Mô hình tất cả nguyên tử (All-atom) Trường lực CHARMM36m, AMBER99SB-ILDN, OPLS-AA/L kết hợp mô hình nước tường minh TIP3P. Tối ưu hóa cấu trúc, đánh giá độ ổn định và tính chất hóa lý chi tiết của trimer, tetramer, dodecamer và sợi Aβ1-42.
Kỹ thuật lấy mẫu nâng cao Động lực học phân tử trao đổi bản sao (REMD - Replica Exchange Molecular Dynamics). Vượt qua các rào cản thế năng, tránh bị mắc kẹt tại cực tiểu cục bộ trong không gian năng lượng tự do.
Mô phỏng Docking phân tử Chương trình ClusPro, PIPER scoring function (tính điểm dựa trên năng lượng đẩy, hút vdW, tĩnh điện và khử dung môi DARS). Tạo cấu hình ban đầu ngẫu nhiên cho trimer và dodecamer Aβ1-42 từ việc kết hợp monomer và dimer.
Tích phân chuyển động & Ràng buộc Thuật toán Verlet, Verlet vận tốc, Leap-frog; thuật toán ràng buộc LINCS; lưới hạt Ewald (PME). Giải phương trình chuyển động cơ học cổ điển, duy trì liên kết holonomic, xử lý tương tác tĩnh điện tầm xa với điều kiện biên tuần hoàn (PBC).
Mô hình màng tế bào thần kinh Mô hình màng D1 (Disease 1): DPPC, POPC, SM (Sphingomyelin), CHL1 (Cholesterol), GM1 (Ganglioside GM1). Tái tạo màng lipid kép không đồng nhất của bệnh nhân Alzheimer để khảo sát tương tác với protein.
Công cụ phân tích cấu trúc & Nhiệt động WHAM (Weighted Histogram Analysis Method), dPCA (dihedral Principal Component Analysis), DSSP, STRIDE. Xây dựng bề mặt năng lượng tự do (PMF), phân tích thành phần cấu trúc bậc hai, cụm cấu hình (clustering).
Đại lượng vật lý đo lường RMSD, RMSF, SASA, hSASA (diện tích kị nước tiếp xúc dung môi), CSS (tiết diện tán xạ), môđun nén diện tích, môđun uốn màng. Định lượng mức độ co cụm, độ lệch hình học, hình dạng phỏng cầu, diện tích tiếp xúc dung môi và tính đàn hồi màng.

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

Chương Tên chương trong luận án Nội dung trọng tâm
Chương 1 Giới thiệu Trình bày bối cảnh bệnh lý AD, tính cấp thiết của việc nghiên cứu cấu trúc oligomer Aβ1-42, lý do sử dụng mô phỏng MD, mục tiêu và bố cục của luận án.
Chương 2 Tổng quan Tổng quan chi tiết về bệnh Alzheimer, cơ chế cắt APP (lộ trình amyloidogenic và non-amyloidogenic), các trạng thái cấu trúc của Aβ (monomer, oligomer, tiền sợi, sợi hoàn chỉnh) và các giả thuyết bệnh học (amyloid, tau).
Chương 3 Cơ sở lý thuyết và phương pháp nghiên cứu Trình bày cơ sở lý thuyết cấu trúc protein, phân loại amino acid, trường lực hạt thô UNRES, các trường lực tất cả nguyên tử (AMBER, CHARMM, OPLS), thuật toán MD, REMD, docking và các công cụ phân tích số liệu.
Chương 4 Kết quả và thảo luận Công bố và phân tích kết quả mô phỏng cấu trúc trimer Aβ1-42, tetramer Aβ1-42, tương tác của dodecamer và sợi Aβ1-42 với màng tế bào thần kinh D1.
Chương 5 Kết luận Tổng kết các phát hiện chính của luận án và đề xuất hướng phát triển tiếp theo.

Chương 1: Giới thiệu

Chương này tổng kết gánh nặng y tế của bệnh Alzheimer và sự cấp thiết trong phát triển hiểu biết cấu trúc vi mô. Đặt vấn đề về việc các oligomer Aβ1-42 hòa tan bậc thấp là tác nhân chính gây độc tế bào thần kinh, nhưng thực nghiệm chưa thể giải cấu trúc nguyên tử do tính mất trật tự và đời sống ngắn. Xác lập 3 mục tiêu nghiên cứu và định hướng sử dụng mô phỏng tính toán làm phương pháp tiếp cận chủ đạo.

Chương 2: Tổng quan

Chương 2 phân tích cơ chế phân tử của AD:

  • Quá trình xử lý APP: Lộ trình nonamyloidogenic thông qua $\alpha$-secretase tạo ra sAPP-$\alpha$ và đoạn p3; lộ trình amyloidogenic thông qua $\beta$-secretase (tạo sAPP-$\beta$ và $\beta$CTF) và $\gamma$-secretase cắt tại các vị trí khác nhau trong màng tạo ra các biến thể Aβ (Aβ38, Aβ40, Aβ42, Aβ43).
  • So sánh Aβ1-40 và Aβ1-42: Khác biệt ở 2 amino acid kị nước ở đuôi C (Ile41 và Ala42) làm tăng đáng kể tốc độ tự tập hợp và độc tính thần kinh của Aβ1-42.
  • Vai trò của oligomer so với sợi: Phân tích các dữ liệu lâm sàng cho thấy việc loại bỏ mảng bám sợi không ngăn chặn được tiến triển AD, trong khi nồng độ oligomer hòa tan tương quan thuận rõ rệt với mức độ suy giảm nhận thức và tổn thương synapse.
  • Cơ chế đầu độc màng: Khái quát giả thuyết oligomer tạo lỗ thủng/kênh dẫn ion không kiểm soát trên màng tế bào thần kinh, gây mất cân bằng nội môi ion Ca2+, stress oxy hóa và làm chết neuron.

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

Chương 3 trình bày toán học và vật lý của các phương pháp mô phỏng:

  • Xây dựng hàm thế năng trong UNRES và trường lực all-atom (thành phần liên kết, góc hóa trị, góc nhị diện, góc improper, tương tác tĩnh điện Coulomb và tương tác van der Waals Lennard-Jones 12-6).
  • Quy trình mô phỏng docking tạo phức hợp ban đầu.
  • Thiết lập hệ thống REMD cho trimer và tetramer: phân bố nhiệt độ, kiểm tra sự trao đổi bản sao và đánh giá sự hội tụ thông qua nhiệt dung riêng ($C_v$), entropy cấu hình, và phân tích chuyển dời trạng thái.
  • Xây dựng hệ màng D1 phức hợp: DPPC/POPC/SM/Cholesterol/GM1 với tỷ lệ sinh học, cân bằng hệ trong hộp nước chứa ion sinh lý.

Chương 4: Kết quả và thảo luận

1. Cấu trúc và tính chất hóa lý của Trimer Aβ1-42

  • Khảo sát hai thiết lập mô phỏng REMD all-atom với trường lực CHARMM36m và nước TIP3P:
    • Mô phỏng 1 (bắt đầu từ docking ngẫu nhiên monomer vào dimer): Trimer hình thành các cụm cấu trúc có motif dạng thùng $\beta$ ($\beta$-barrel) nhỏ xuất hiện tự phát. Đây là lần đầu tiên cấu trúc thùng $\beta$ tự phát được quan sát trực tiếp từ mô phỏng mà không cần áp đặt cấu hình ban đầu. Trimer tạo ra các lỗ xuyên tâm có bán kính khoảng 1.1 Å, đủ lớn để phân tử nước và ion Ca2+ thâm nhập.
    • Mô phỏng 2 (bắt đầu từ cấu hình tách ra từ sợi hoàn chỉnh 2NAQ): Không xuất hiện motif thùng $\beta$; đuôi C của các chuỗi Aβ1-42 duy trì bó cấu trúc ổn định, kế thừa từ trật tự của sợi thực nghiệm.
  • So sánh tính chất: Trong cả hai mô phỏng, trimer Aβ1-42 có cấu trúc bó chặt và kém trật tự hơn nhiều so với sợi hoàn chỉnh. Các amino acid mạch thơm ở trimer tiếp xúc với dung môi nhiều hơn (hSASA cao hơn), phù hợp với các dữ liệu thực nghiệm phổ Raman.

2. Cấu trúc và tính chất hóa lý của Tetramer Aβ1-42

  • Kết hợp mô phỏng REMD hạt thô UNRES (thang thời gian micro giây) và MD tất cả nguyên tử (AMBER99SB-ILDN, OPLS-AA/L):
    • Tetramer Aβ1-42 thể hiện đặc tính đa hình (polymorphism) rõ rệt, giải thích sự tồn tại của nhiều lộ trình kết tụ khác nhau.
    • Cấu hình tetramer phân hóa thành hai nhóm chuỗi: "chuỗi ngoài" (linh động, mất trật tự) và "chuỗi lõi" (bó chặt, hình thành các dải $\beta$).
    • Tetramer có dạng hình phỏng cầu dẹt (oblate ellipsoid), khác biệt căn bản với cấu trúc dạng thanh kéo dài giả một chiều của sợi hoàn chỉnh.
    • Phân tích tương tác dung môi - chất hòa tan chỉ ra rằng tương tác với nước giữ cho tetramer bó chặt hơn nhưng bên trong "ít khô" hơn (chứa nhiều phân tử nước nội tại hơn) so với lõi kị nước của sợi hoàn chỉnh.

3. Tương tác giữa Dodecamer, Sợi hoàn chỉnh Aβ1-42 và Màng tế bào thần kinh

  • Mô phỏng MD all-atom hệ màng D1 tương tác với dodecamer (12 chuỗi Aβ1-42 dạng vòng khuyên) và đoạn sợi hoàn chỉnh (12 chuỗi):
    • Dodecamer tiếp cận và làm biến dạng bề mặt màng lipid kép mạnh hơn rõ rệt so với sợi hoàn chỉnh.
    • Tương tác protein - màng bị chi phối bởi lực đẩy tĩnh điện giữa các điện tích âm của peptide Aβ và đầu phân tử lipid GM1 (Monosialotetrahexosylganglioside).
    • Về tính chất cơ học màng: Dodecamer Aβ1-42 làm tăng đáng kể môđun nén diện tích của màng, nhưng ít ảnh hưởng đến môđun uốn; sợi hoàn chỉnh tương tác yếu nên làm thay đổi không đáng kể các tính chất đàn hồi này.
    • Sự biến dạng cục bộ và tạo các vết nứt bề mặt màng bởi dodecamer giải thích cơ chế phân tử ở các sự kiện ban đầu dẫn tới phá vỡ tính thấm màng tế bào thần kinh.

Chương 5: Kết luận và hướng phát triển của luận án

Tổng kết toàn bộ các luận điểm đạt được về mặt cấu trúc và động lực học của oligomer Aβ1-42; khẳng định tính hợp lý của phương pháp mô phỏng tính toán; nêu bật triển vọng ứng dụng các mô hình cấu trúc này để xây dựng oligomer bậc cao hơn và phục vụ sàng lọc phân tử ức chế kết tụ amyloid.

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

Luận án đã đạt được các đóng góp mới về mặt khoa học lý thuyết và tính toán sinh lý phân tử:

  1. Phát hiện sự hình thành tự phát của motif thùng $\beta$ trong trimer Aβ1-42: Lần đầu tiên qua mô phỏng động lực học phân tử (bắt đầu từ cấu hình docking không áp đặt trước), cấu trúc dạng thùng $\beta$ nhỏ và lỗ rỗng kích thước 1.1 Å có khả năng cho nước và ion Ca2+ đi qua được ghi nhận, cung cấp bằng chứng cấu trúc cho giả thuyết tạo kênh ion của oligomer.
  2. Xác định tính chất đa hình và cấu trúc phân tử của tetramer Aβ1-42: Chứng minh tetramer tồn tại ở dạng phỏng cầu dẹt, có sự phân tách giữa chuỗi lõi và chuỗi ngoài; làm rõ cơ chế tương tác với dung môi nước là nguyên nhân khiến tetramer có cấu trúc phân tử và sự phân bố độ khô bên trong khác biệt hoàn toàn với sợi hoàn chỉnh.
  3. Làm sáng tỏ cơ chế tương tác ban đầu giữa dodecamer Aβ1-42 và màng tế bào thần kinh bệnh lý: Chứng minh dodecamer làm rối loạn trật tự lipid, làm tăng môđun nén diện tích và gây biến dạng màng mạnh hơn sợi hoàn chỉnh; phát hiện vai trò chi phối của lực đẩy tĩnh điện giữa peptide Aβ và ganglioside GM1 trong động học bám màng.
  4. Đề xuất giả thuyết về sự khác biệt cấu trúc giữa oligomer và sợi: Cung cấp cơ sở lý thuyết cho thấy oligomer hòa tan và sợi hoàn chỉnh có thể thuộc các lộ trình cấu trúc độc lập thay vì chỉ là sự chuyển tiếp tuyến tính đơn thuần.

Kiến nghị và đề xuất: Sử dụng các mô hình cấu trúc trimer, tetramer và dodecamer đã được tối ưu hóa trong luận án làm cấu trúc bia (target) để thiết kế các phối tử, peptide mô phỏng hoặc hợp chất ức chế giai đoạn kết tụ sớm của Aβ1-42.

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

  • Hạn chế: Do giới hạn tài nguyên tính toán và thang thời gian mô phỏng (tính bằng micro giây), quá trình mô phỏng chưa quan sát được toàn bộ sự kiện đục thủng màng hoàn chỉnh để hình thành kênh dẫn ion xuyên màng hoàn chỉnh, mà chỉ dừng lại ở việc nắm bắt các sự kiện biến dạng và rối loạn cấu trúc màng ban đầu.
  • Hướng nghiên cứu tiếp:
    • Mở rộng mô phỏng động lực học phân tử cho các oligomer bậc cao hơn (như hexamer, octamer, 24-mer).
    • Tăng thang thời gian mô phỏng hoặc sử dụng các kỹ thuật mô phỏng gia tốc để theo dõi tiến trình hình thành hoàn chỉnh lỗ xuyên màng sinh học.
    • Khảo sát ảnh hưởng của các đột biến di truyền (như đột biến Iowa, Arctic...) lên cấu trúc oligomer và tương tác màng.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho:

  • Nghiên cứu sinh, học viên cao học, giảng viên ngành Vật lý sinh học, Vật lý tính toán, Hóa tin học, và Sinh học cấu trúc: Cung cấp phương pháp luận kết hợp giữa trường lực hạt thô UNRES và trường lực tất cả nguyên tử (CHARMM36m, AMBER99SB-ILDN, OPLS-AA/L), kỹ thuật REMD và xử lý dữ liệu bằng WHAM, dPCA.
  • Các nhà nghiên cứu y dược học và bệnh thoái hóa thần kinh: Cung cấp mô hình tọa độ nguyên tử 3 chiều của các oligomer độc tính Aβ1-42 (trimer, tetramer, dodecamer) và cơ chế tương tác với màng lipid đa thành phần chứa GM1, phục vụ công tác sàng lọc ảo và phát triển thuốc điều trị Alzheimer nhắm đích oligomer.

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

1. Vì sao cấu hình ban đầu lại tạo ra sự khác biệt lớn trong mô phỏng trimer Aβ1-42?

Trong mô phỏng 1 (bắt đầu từ docking ngẫu nhiên monomer-dimer), hệ đại diện cho các oligomer mới hình thành từ dung dịch, các chuỗi có độ tự do cao và tự sắp xếp tạo thành motif thùng $\beta$ tự phát. Trong mô phỏng 2 (bắt đầu từ cấu trúc tách ra từ sợi 2NAQ), hệ đại diện cho trạng thái oligomer đã ổn định trong cân bằng sợi; các liên kết hydro và cấu trúc đuôi C bền vững kế thừa từ sợi hoàn chỉnh ngăn cản chuỗi chuyển dịch sang motif thùng $\beta$ trong thang thời gian mô phỏng.

2. Cấu trúc tetramer Aβ1-42 khác biệt căn bản như thế nào so với sợi amyloid hoàn chỉnh?

Tetramer Aβ1-42 có dạng hình phỏng cầu dẹt với cấu trúc mất trật tự chiếm ưu thế, gồm các chuỗi ngoài linh hoạt và chuỗi lõi, tiếp xúc với dung môi nhiều hơn và bên trong ít khô hơn (chứa nhiều phân tử nước). Ngược lại, sợi hoàn chỉnh có cấu trúc kéo dài giả một chiều dạng "cross beta", các lớp $\beta$ xếp song song có trật tự cao và lõi kị nước bên trong bị khử dung môi gần như hoàn toàn.

3. Lipid Ganglioside GM1 đóng vai trò gì trong tương tác giữa Aβ1-42 và màng tế bào thần kinh?

Trong mô hình màng D1, tương tác tĩnh điện giữa các gốc mang điện tích âm của peptide Aβ1-42 và đầu phân tử lipid GM1 tạo ra lực đẩy tĩnh điện chi phối pha tiếp cận ban đầu của protein lên bề mặt màng, ảnh hưởng trực tiếp đến sự phân bố lipid và độ sâu các vết nứt bề mặt màng.

4. Ưu thế của việc sử dụng trường lực hạt thô UNRES trước khi mô phỏng all-atom cho tetramer là gì?

Trường lực UNRES gộp các nhóm nguyên tử thành hạt đại diện (C$\alpha$, SC, p) dựa trên thế năng lực trung bình (PMF) vật lý, cho phép mở rộng không gian lấy mẫu cấu hình hiệu dụng nhanh hơn 3–4 bậc độ lớn so với all-atom. Điều này giúp hệ vượt qua các rào cản năng lượng để tìm kiếm các cụm cấu hình ổn định của tetramer mà không đòi hỏi tài nguyên tính toán quá lớn, sau đó các cấu trúc đại diện được đưa vào mô phỏng all-atom để chuẩn hóa chi tiết nguyên tử.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Hoàng Linh đã khảo sát cấu trúc ba chiều và động lực học phân tử của các dạng kết tụ oligomer bậc thấp của peptide Aβ1-42 (trimer, tetramer, dodecamer) và tương tác của chúng với màng tế bào thần kinh đa thành phần. Bằng việc kết hợp hiệu quả các phương pháp mô phỏng hạt thô UNRES, mô phỏng tất cả nguyên tử và kỹ thuật trao đổi bản sao REMD, nghiên cứu đã phát hiện sự xuất hiện tự phát của motif thùng $\beta$ ở trimer, làm rõ tính đa hình dạng phỏng cầu của tetramer và chứng minh dodecamer làm biến dạng màng tế bào thần kinh mạnh hơn sợi hoàn chỉnh. Các phát hiện này cung cấp cơ sở phân tử quan trọng ủng hộ giả thuyết về độc tính vượt trội của oligomer trong cơ chế bệnh sinh Alzheimer, đồng thời đóng góp bộ dữ liệu cấu trúc phục vụ các nghiên cứu định hướng thuốc điều trị trong tương lai.