Luận án tiến sĩ computational study of amyloid beta protein in implicit and explicit solvent models probing the initial stages of aggregation

Luận án tiến sĩ nghiên cứu computational study of amyloid beta protein in implicit and explicit solvent models probing the, phân tích chuyên sâu, xây dựng mô hình lý thuyết, đề

• Tác giả: Bogdan Tarus
• Người hướng dẫn: Prof. John E. Straub, Professor of Chemistry
• Trường học: Boston University
• Chuyên ngành: Chemistry
• Đề tài: Nghiên Cứu Tính Toán Về Protein Amyloid Beta Trong Mô Hình Dung Môi: Khám Phá Giai Đoạn Đầu Của Sự Kết Tụ
• Loại tài liệu: dissertation
• Năm xuất bản: 2007
• Địa điểm: Boston

Mục lục chi tiết

• 1. CHƯƠNG 1: INTRODUCTION
• 2. CHƯƠNG 2: PROBING THE INITIAL STAGE OF AGGREGATION OF THE A 9_3;-PROTEIN: ASSESSING THE PROPENSITY FOR PEPTIDE DIMERIZATION
• 2.1. Computational Model and Methods
• 2.1.1. Dimer structure generation using a docking protocol
• 2.1.2. Desolvation energy screening
• 2.1.3. Secondary Structure Analysis
• 2.2. Results and Discussion
• 2.2.1. Generation of dimer structures
• 2.2.2. Potential of Mean Force
• 2.2.3. Stability of dimer
• 2.2.4. Dynamical fluctuations in the y-dimer
• 2.2.5. Time dependence of secondary structure fluctuations in the homodimer
• 2.3. Summary and Conclusions
• 3. CHƯƠNG 3: A COMPARATIVE STUDY OF THE STRUCTURE AND THERMODYNAMICS OF THE A{io_35- PROTEIN IN DIFFERENT HYDRATION MODELS
• 3.1. Computational Models and Methods
• 3.1.1. Explicit solvent molecular dynamics
• 3.1.2. Implicit solvent molecular dynamics
• 3.1.3. Computational estimation of pK
• 3.2. Results and analysis
• 3.2.1. Radius of gyration
• 3.2.2. Root-mean-square fluctuations
• 3.2.3. Lipari-Szabo generalized order parameter
• 3.2.4. Solvation self-energy
• 3.2.5. Intra-peptide hydrogen bonds
• 3.2.6. Computational estimation of pK
• 3.3. Summary and Conclusions
• 4. CHƯƠNG 4: DYNAMICS OF ASP23—LYS28 SALT BRIDGE FORMATION IN AĐỊO_SS MONOMERS
• 4.1. Compact A(@io-35 undergoes large structural fluctuations
• 4.2. Computed values of pK are in accord with experiments
• 4.3. Structures with intact salt-bridge D23-K28 are not the most stable
• 4.4. Multiple basins are populated in Àđịo ss monomer
• 4.5. Hydrogen-bonds with water increase the desolvation barrier of D23 and K28
• 4.6. Burial of K28 involves a large free energy cost
• 4.7. The structure of water around residues D23 and K28
• 4.8. Summary and Conclusions
• 5. CHƯƠNG 5: THE COMPETITION BETWEEN THE ELECTROSTATIC AND HYDROPHOBIC INTRA-PEPTIDE INTERACTIONS IN THE AØS¡_SO-PROTEIN
• 5.1. Computational Methods and Models
• 5.1.1. Simulation model and methods
• 5.1.2. Computational estimation of pK
• 5.2. Results and analysis
• 5.2.1. pK values indicate weak intra-peptide electrostatic interactions
• 5.2.2. Inter-titratable side-chain distances
• 5.2.3. Intra-peptide hydrophobic contacts
• 5.2.4. Lys28 makes transitive contacts with the peptide backbone
• 5.2.5. Intra-peptide folding elements
• 5.3. Summary and Conclusions
• Bibliography
• Curriculum Vitae

Tóm tắt

I. Tổng Quan Nghiên Cứu Tính Toán Protein Amyloid Beta Aβ Peptide

Nghiên cứu về protein amyloid beta (Aβ) đóng vai trò then chốt trong việc tìm hiểu cơ chế bệnh Alzheimer. Bài viết này trình bày tổng quan về các nghiên cứu tính toán sự hình thành các tập hợp oligomeric của Aβ, tập trung vào giai đoạn đầu của quá trình kết tụ protein. Các phương pháp mô phỏng động lực học phân tử và phương pháp Monte Carlo được sử dụng để khám phá các cấu trúc trung gian và cơ chế tương tác protein-dung môi ảnh hưởng đến tính ổn định protein. Mục tiêu là xác định các yếu tố quyết định sự hình thành các cấu trúc protein độc hại liên quan đến bệnh Alzheimer. Nghiên cứu này sử dụng cả mô hình dung môi tường minh và ẩn để đánh giá ảnh hưởng của môi trường đến quá trình kết tụ protein. Dữ liệu từ các nghiên cứu mô hình hóa phân tử này cung cấp thông tin chi tiết về cơ chế kết tụ và có thể hỗ trợ phát triển các liệu pháp nhắm mục tiêu vào các giai đoạn đầu của bệnh.

1.1. Vai trò của Protein Amyloid Beta trong Bệnh Alzheimer

Protein amyloid beta (Aβ) được cho là đóng vai trò trung tâm trong sự phát triển của bệnh Alzheimer. Sự kết tụ protein Aβ thành các mảng bám amyloid là một trong những dấu hiệu bệnh lý đặc trưng của bệnh. Nghiên cứu tập trung vào việc tìm hiểu cơ chế kết tụ Aβ và các yếu tố ảnh hưởng đến quá trình này. Các tập hợp oligomeric Aβ, được hình thành trong giai đoạn đầu kết tụ, được cho là độc hại hơn so với các mảng bám trưởng thành. Do đó, việc nghiên cứu giai đoạn đầu kết tụ là rất quan trọng để phát triển các phương pháp điều trị hiệu quả.

1.2. Ứng Dụng Nghiên Cứu Tính Toán trong Nghiên Cứu Protein

Nghiên cứu tính toán cung cấp một công cụ mạnh mẽ để nghiên cứu cấu trúc protein và tương tác protein-dung môi ở cấp độ nguyên tử. Các phương pháp mô phỏng động lực học phân tử cho phép mô phỏng sự chuyển động của các phân tử protein theo thời gian, cung cấp thông tin chi tiết về tính ổn định protein và tính chất nhiệt động lực học. Phương pháp Monte Carlo được sử dụng để khám phá không gian cấu hình rộng lớn và xác định các cấu trúc năng lượng thấp. Các kết quả mô hình hóa phân tử có thể bổ sung cho các nghiên cứu thực nghiệm và cung cấp thông tin chi tiết về cơ chế kết tụ.

II. Thách Thức Nghiên Cứu Giai Đoạn Đầu Kết Tụ Protein Aβ Peptide

Nghiên cứu giai đoạn đầu kết tụ của protein amyloid beta (Aβ) gặp nhiều thách thức do tính chất động và phức tạp của quá trình. Các tập hợp oligomeric Aβ chỉ tồn tại trong thời gian ngắn và khó quan sát bằng các phương pháp thực nghiệm truyền thống. Việc mô phỏng quá trình kết tụ protein đòi hỏi tài nguyên tính toán lớn và các mô hình dung môi chính xác để mô tả tương tác protein-dung môi. Một thách thức khác là xác định các cấu trúc trung gian quan trọng trên đường dẫn phản ứng kết tụ protein. Các nghiên cứu phân tích đường dẫn phản ứng và phân tích cụm được sử dụng để giải quyết vấn đề này. Việc hiểu rõ các yếu tố ảnh hưởng đến năng lượng tự do của quá trình kết tụ protein là rất quan trọng để phát triển các phương pháp điều trị hiệu quả.

2.1. Khó Khăn trong Mô Phỏng Động Lực Học Phân Tử Kết Tụ

Mô phỏng động lực học phân tử quá trình kết tụ protein Aβ đòi hỏi thời gian tính toán lớn do kích thước hệ thống và thời gian mô phỏng cần thiết để quan sát các sự kiện kết tụ protein. Việc lựa chọn mô hình dung môi phù hợp cũng là một thách thức, vì tương tác protein-dung môi đóng vai trò quan trọng trong quá trình kết tụ protein. Các mô hình dung môi tường minh cung cấp độ chính xác cao hơn nhưng đòi hỏi tài nguyên tính toán lớn hơn so với các mô hình dung môi ẩn.

2.2. Xác Định Cấu Trúc Trung Gian trong Quá Trình Kết Tụ

Việc xác định các cấu trúc trung gian trên đường dẫn phản ứng kết tụ protein là một thách thức quan trọng. Các cấu trúc này thường không ổn định và khó quan sát bằng các phương pháp thực nghiệm. Các phương pháp phân tích đường dẫn phản ứng và phân tích cụm được sử dụng để xác định các cấu trúc trung gian quan trọng và hiểu rõ cơ chế kết tụ. Việc xác định các cấu trúc này có thể giúp phát triển các liệu pháp nhắm mục tiêu vào các giai đoạn đầu của quá trình kết tụ protein.

III. Phương Pháp Mô Phỏng Động Lực Học Phân Tử Protein Aβ Peptide

Nghiên cứu này sử dụng mô phỏng động lực học phân tử để nghiên cứu tính ổn định protein và tương tác protein-dung môi của protein amyloid beta (Aβ). Các mô hình dung môi tường minh và ẩn được sử dụng để đánh giá ảnh hưởng của môi trường đến quá trình kết tụ protein. Phương pháp Monte Carlo được sử dụng để tạo ra các cấu trúc ban đầu cho mô phỏng động lực học phân tử. Các kỹ thuật lấy mẫu ô dù được sử dụng để tính toán tiềm năng tương tác trung bình (PMF) liên quan đến sự kết tụ protein. Các kết quả mô phỏng được phân tích để xác định các yếu tố quyết định tính ổn định protein và cơ chế kết tụ.

3.1. Sử Dụng Mô Hình Dung Môi Tường Minh và Ẩn

Nghiên cứu sử dụng cả mô hình dung môi tường minh (ví dụ: TIP3P) và mô hình dung môi ẩn (ví dụ: ACE, GBORN) để mô tả tương tác protein-dung môi. Mô hình dung môi tường minh mô tả các phân tử dung môi một cách rõ ràng, cung cấp độ chính xác cao hơn nhưng đòi hỏi tài nguyên tính toán lớn hơn. Mô hình dung môi ẩn mô tả ảnh hưởng của dung môi một cách gần đúng, giảm chi phí tính toán nhưng có thể kém chính xác hơn.

3.2. Kỹ Thuật Lấy Mẫu Ô Dù Tính Tiềm Năng Tương Tác Trung Bình

Kỹ thuật lấy mẫu ô dù được sử dụng để tính toán tiềm năng tương tác trung bình (PMF) liên quan đến sự kết tụ protein. Kỹ thuật này cho phép khám phá không gian cấu hình rộng lớn và xác định các trạng thái năng lượng thấp. PMF cung cấp thông tin chi tiết về năng lượng tự do của quá trình kết tụ protein và có thể được sử dụng để xác định các yếu tố quyết định tính ổn định protein.

IV. Kết Quả Nghiên Cứu Tính Toán Tương Tác Protein Dung Môi Aβ

Các kết quả nghiên cứu tính toán cho thấy tương tác protein-dung môi đóng vai trò quan trọng trong việc điều chỉnh tính ổn định protein và cơ chế kết tụ của protein amyloid beta (Aβ). Các mô hình dung môi khác nhau dẫn đến các cấu trúc protein khác nhau và các tính chất nhiệt động lực học khác nhau. Nghiên cứu cũng xác định các dư lượng quan trọng trong chuỗi Aβ đóng vai trò quan trọng trong quá trình kết tụ protein. Các kết quả này cung cấp thông tin chi tiết về cơ chế kết tụ và có thể hỗ trợ phát triển các liệu pháp nhắm mục tiêu vào các giai đoạn đầu của bệnh Alzheimer.

4.1. Ảnh Hưởng của Mô Hình Dung Môi Đến Cấu Trúc Protein

Các mô hình dung môi khác nhau dẫn đến các cấu trúc protein khác nhau. Mô hình dung môi tường minh thường tạo ra các cấu trúc protein linh hoạt hơn so với mô hình dung môi ẩn. Các mô hình dung môi ẩn có xu hướng tạo ra các cấu trúc protein nhỏ gọn hơn do sự đánh giá thấp năng lượng tự do hòa tan.

4.2. Vai Trò của Các Dư Lượng Quan Trọng Trong Kết Tụ Protein

Nghiên cứu xác định các dư lượng quan trọng trong chuỗi Aβ đóng vai trò quan trọng trong quá trình kết tụ protein. Ví dụ, các dư lượng kỵ nước có xu hướng tập hợp lại với nhau, thúc đẩy sự hình thành các tập hợp oligomeric. Các dư lượng tích điện có thể tương tác với nhau hoặc với dung môi, ảnh hưởng đến tính ổn định protein.

V. Ứng Dụng Nghiên Cứu Tính Toán Phát Triển Thuốc Trị Alzheimer

Các kết quả nghiên cứu tính toán có thể được sử dụng để phát triển các liệu pháp nhắm mục tiêu vào các giai đoạn đầu của bệnh Alzheimer. Bằng cách hiểu rõ cơ chế kết tụ và các yếu tố ảnh hưởng đến tính ổn định protein, các nhà khoa học có thể thiết kế các phân tử nhỏ có thể ức chế sự kết tụ protein hoặc ổn định các cấu trúc protein không độc hại. Các phương pháp mô hình hóa phân tử có thể được sử dụng để sàng lọc các ứng cử viên thuốc tiềm năng và tối ưu hóa cấu trúc của chúng.

5.1. Thiết Kế Phân Tử Nhỏ Ức Chế Kết Tụ Protein Aβ

Các phân tử nhỏ có thể được thiết kế để ức chế sự kết tụ protein Aβ bằng cách liên kết với Aβ và ngăn chặn sự hình thành các tập hợp oligomeric độc hại. Các phân tử này có thể nhắm mục tiêu vào các vị trí liên kết quan trọng trên bề mặt Aβ hoặc ổn định các cấu trúc protein không độc hại.

5.2. Sàng Lọc Ứng Cử Viên Thuốc Tiềm Năng Bằng Mô Hình Hóa

Các phương pháp mô hình hóa phân tử có thể được sử dụng để sàng lọc các ứng cử viên thuốc tiềm năng bằng cách mô phỏng sự tương tác của chúng với Aβ. Các phương pháp này có thể dự đoán ái lực liên kết và hiệu quả ức chế của các ứng cử viên thuốc, giúp xác định các phân tử có triển vọng nhất để phát triển thành thuốc.

VI. Kết Luận và Hướng Nghiên Cứu Tương Lai Protein Amyloid Beta

Nghiên cứu tính toán về protein amyloid beta (Aβ) trong mô hình dung môi đã cung cấp thông tin chi tiết quan trọng về giai đoạn đầu kết tụ và cơ chế kết tụ. Các kết quả này có thể hỗ trợ phát triển các liệu pháp nhắm mục tiêu vào các giai đoạn đầu của bệnh Alzheimer. Các hướng nghiên cứu tương lai bao gồm việc sử dụng các mô hình dung môi chính xác hơn, phát triển các phương pháp phân tích đường dẫn phản ứng tiên tiến hơn và kết hợp các kết quả mô hình hóa phân tử với các dữ liệu thực nghiệm để hiểu rõ hơn về cơ chế kết tụ.

6.1. Phát Triển Mô Hình Dung Môi Chính Xác Hơn Cho Aβ

Việc phát triển các mô hình dung môi chính xác hơn là rất quan trọng để mô tả chính xác tương tác protein-dung môi và tính ổn định protein. Các mô hình dung môi phân cực có thể cung cấp độ chính xác cao hơn so với các mô hình dung môi không phân cực.

6.2. Kết Hợp Dữ Liệu Mô Hình Hóa và Thực Nghiệm

Việc kết hợp các kết quả mô hình hóa phân tử với các dữ liệu thực nghiệm, chẳng hạn như dữ liệu từ phổ cộng hưởng từ hạt nhân (NMR) và kính hiển vi điện tử, có thể cung cấp một bức tranh toàn diện hơn về cơ chế kết tụ và cấu trúc protein.

Nội dung chính

Tổng quan về luận án

Bệnh Alzheimer (Alzheimer's Disease - AD) là một trong những hội chứng suy giảm trí nhớ và thoái hóa thần kinh tàn khốc nhất hiện nay, được đặc trưng về mặt mô bệnh học bởi sự tích tụ của các mảng xơ ngoại bào cấu thành từ protein amyloid $\beta$ (A$\beta$). Luận án tiến sĩ mang tên "Computational Study of Amyloid Beta Protein in Implicit and Explicit Solvent Models: Probing the Initial Stages of Aggregation" do nghiên cứu sinh Bogdan Tarus thực hiện tại Đại học Boston (Boston University) dưới sự hướng dẫn của Giáo sư John E. Straub, là công trình tiên phong giải mã cấu trúc phân tử và động học nhiệt động lực học của các trạng thái trung gian sơ khởi trong quá trình kết tụ peptide A$\beta$. Bối cảnh khoa học của luận án xuất phát từ sự thay đổi mô hình nhận thức y sinh học: thay vì chỉ tập trung vào các sợi tơ amyloid trưởng thành (mature fibrils), các bằng chứng thực nghiệm hiện đại chỉ ra rằng chính các cấu trúc oligomer trọng lượng phân tử thấp (Low Molecular Weight - LMW oligomers) mới là tác nhân gây độc tế bào thần kinh chính trong giai đoạn đầu của bệnh.

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án giải quyết là sự thiếu vắng hiểu biết chi tiết ở cấp độ nguyên tử về cơ chế tạo dimer ban đầu và vai trò của dung môi đối với sự ổn định của cấu trúc đơn phân (monomer) cũng như phức hợp bậc thấp. Như tác giả đã khẳng định mục tiêu nghiên cứu trong văn bản: "The objective of this work was to determine the stability and role of the Aβ-protein monomers and low-order oligomers as metastable intermediates on the pathway for formation of larger aggregates and fibrils." Mặc dù cấu trúc dung dịch NMR của đoạn peptide $A\beta_{10-35}$ đã được xác định trước đó, nhưng các rào cản năng lượng tự do chi phối quá trình nhị hợp hóa (dimerization), sự chuyển pha cấu trúc từ cuộn ngẫu nhiên (collapsed random coil) sang phiến $\beta$ song song (in-register parallel $\beta$-sheets), và tính chính xác của các mô hình dung môi vi phân (implicit vs. explicit solvent models) vẫn chưa được làm sáng tỏ một cách có hệ thống.

Luận án thiết lập và kiểm chứng 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H) tương ứng:

1. RQ1 & H1: Động lực phân tử nào chi phối quá trình tạo dimer ban đầu giữa các monomer $A\beta_{10-35}$? Giả thuyết H1: Sự khử dung môi kỵ nước (hydrophobic desolvation) tại bề mặt tiếp xúc liên monomer là yếu tố quyết định năng lượng tự do liên kết, vượt trội so với các tương tác tĩnh điện đặc hiệu.
2. RQ2 & H2: Các mô hình dung môi ngầm định (implicit solvent models) có mô phỏng chính xác cấu trúc mềm dẻo của monomer $A\beta_{10-35}$ so với dung môi tường minh (explicit TIP3P)? Giả thuyết H2: Các mô hình dung môi liên tục tĩnh điện (GBORN, ACE) đánh giá thấp năng lượng tự solvat hóa phân tử, dẫn đến sự hình thành quá mức các liên kết hydro nội chuỗi phi tự nhiên làm co cụm cấu trúc peptide.
3. RQ3 & H3: Cầu muối nội phân tử Asp23-Lys28 ($D23-K28$) có hình thành bền vững trong monomer $A\beta_{10-35}$ tự do trong môi trường nước? Giả thuyết H3: Cầu muối $D23-K28$ là một trạng thái chuyển tiếp linh động không chiếm ưu thế tuyệt đối trong monomer đơn lẻ do rào cản khử solvat hóa của chuỗi bên tích điện khi tiếp xúc với nước.
4. RQ4 & H4: Các đột biến di truyền tự nhiên (E22Q Dutch mutant, D23N Iowa mutant) và biến thể K28A tác động như thế nào đến bề mặt năng lượng tự do và cấu trúc nếp gấp $\beta$-turn tại vùng $VGSN(24-27)$? Giả thuyết H4: Đột biến E22Q và K28A tăng cường tương tác kỵ nước nội chuỗi, làm sâu sắc thêm các lòng chảo cấu trúc ổn định, trong khi D23N làm phẳng bề mặt năng lượng tự do và gia tăng tính mềm dẻo của peptide.

Khung lý thuyết của công trình tích hợp Nguyên lý Tự tổ chức Amyloid (Amyloid Self-Organization Principle - ASOP), Nhiệt động lực học Cơ học Thống kê (Statistical Mechanical Thermodynamics), Lý thuyết Điện tĩnh học Môi trường liên tục Poisson-Boltzmann (PBEQ) / Generalized Born (GB), và Lý thuyết Thư giãn NMR Lipari-Szabo. Phạm vi nghiên cứu bao quát các chuỗi peptide amyloid chủ đạo: $A\beta_{10-35}$, $A\beta_{21-30}$, $A\beta_{1-40}$, và $A\beta_{1-42}$ với tổng thời gian mô phỏng động lực học phân tử (MD) hàng trăm nanogiây (ns) cùng hàng nghìn cấu trúc mồi (decoys), tạo ra tác động định lượng đột phá trong việc hiểu rõ quá trình tiền tạo mầm (pre-nucleation) của amyloidosis.

---

Literature Review và Positioning

Quá trình nghiên cứu amyloidosis trong các thập kỷ qua đã chứng kiến nhiều cuộc tranh luận học thuật sâu sắc về động học và cơ chế phân tử. Luồng nghiên cứu thứ nhất, dẫn đầu bởi Teplow và cộng sự (2002), đã chứng minh sự khác biệt trong con đường tạo mầm giữa hai biến thể đồng phân tự nhiên $A\beta_{1-40}$ và $A\beta_{1-42}$, trong đó các oligomer sơ khởi dạng paranuclei đóng vai trò là tiền chất thiết yếu hình thành sợi tơ qua cơ chế trùng hợp tuyến tính. Luồng nghiên cứu thứ hai tập trung vào cấu trúc sợi amyloid trưởng thành thông qua kỹ thuật cộng hưởng từ hạt nhân trạng thái rắn (solid-state NMR) của Tycko và cộng sự (2002), Lynn và cộng sự (2002, 2003), chỉ ra rằng các monomer trong sợi tơ $A\beta_{10-35}$ và $A\beta_{1-40}$ sắp xếp theo cấu trúc phiến $\beta$ song song thẳng hàng (in-register parallel $\beta$-sheets) với một khúc quay $\beta$-turn tại vị trí $VGSN(24-27)$ được ổn định bởi cầu muối nội chuỗi $D23-K28$.

Tuy nhiên, tồn tại một sự mâu thuẫn lớn giữa cấu trúc cuộn ngẫu nhiên co cụm (collapsed coil) của monomer trong dung dịch nước (được xác minh qua phổ NMR dung dịch của Maggio và cộng sự, 2004) và cấu trúc phiến $\beta$ có trật tự cao trong sợi tơ. Câu hỏi đặt ra là: Quá trình chuyển đổi cấu hình diễn ra ở giai đoạn nào, và liệu monomer tự do có chứa sẵn các yếu tố cấu trúc mầm (seed structural motifs) hay không? Một nghiên cứu quốc tế nổi bật của Massi, Straub và Thirumalai (2001) khi phân tích biến thể $A\beta_{10-35}$ dạng tự nhiên (WT) và đột biến Hà Lan (Dutch mutant E22Q) đã nhận định rằng tương tác nội chuỗi không có sự khác biệt quá lớn, mà xu hướng tăng khả năng gây bệnh của đột biến Hà Lan bắt nguồn từ việc giảm tính ưa nước toàn cục. Mặt khác, nghiên cứu của Bitan và cộng sự (2003) về quá trình oligome hóa của $A\beta_{1-40}$ và $A\beta_{1-42}$ cho thấy quá trình oxy hóa dư lượng Met35 làm chậm tốc độ tạo sợi của $A\beta_{1-42}$, chứng minh rằng rào cản năng lượng tự do khử dung môi kỵ nước đóng vai trò then chốt trong quá trình kết tụ.

Luận án của Bogdan Tarus định vị chính xác vào khoảng trống này bằng cách áp dụng phương pháp tiếp cận đa quy mô (multiscale computational approach): kết hợp thuật toán đối xứng hình học bề mặt (shape complementarity docking) với mô phỏng động lực học phân tử tường minh (all-atom explicit solvent MD), tính toán Thế trung bình thế năng (Potential of Mean Force - PMF), và xây dựng Đồ thị mất liên tục năng lượng tự do (Free Energy Transition Disconnectivity Graphs - TRDG). So sánh với hai nghiên cứu quốc tế điển hình là nghiên cứu của Bitan et al. (2003) (vốn dừng lại ở các quan sát động học thực nghiệm diện rộng) và nghiên cứu của Petkova et al. (2002) (chỉ xác định cấu trúc tĩnh của sợi tinh thể), công trình của Tarus cung cấp một bước tiến vượt bậc: mô hình hóa định lượng động học vi mô của quá trình đẩy các phân tử nước (water expulsion) tại bề mặt phân cách, phân tích sự cạnh tranh giữa các lực tĩnh điện và kỵ nước, đồng thời đánh giá nghiêm ngặt giới hạn của các mô hình dung môi ngầm định trong hóa sinh tính toán.

---

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

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

Luận án mang lại những đóng góp nền tảng cho lý thuyết hóa lý sinh học và cơ chế bệnh sinh phân tử:

1. Mở rộng Nguyên lý Tự tổ chức Amyloid (Amyloid Self-Organization Principle - ASOP) của Thirumalai và cộng sự: Luận án chứng minh rằng trong cấu trúc sợi tơ đa phân tử, monomer thỏa mãn ASOP bằng cách tối đa hóa liên kết hydro và cầu muối nội chuỗi để không chôn vùi các điện tích tự do vào môi trường điện môi thấp. Tuy nhiên, ở trạng thái monomer cô lập trong nước, trạng thái năng lượng tự do thấp nhất lại ưu tiên việc solvat hóa hoàn toàn các chuỗi bên phân cực ($D23$ và $K28$), biến cầu muối $D23-K28$ thành một trạng thái chuyển tiếp linh động thay vì một cấu trúc bền vững tiền định.
2. Xác lập mô hình Dimer hóa do Lực kỵ nước dẫn dắt ($\gamma$-dimer Paradigm): Bằng cách so sánh hai cấu trúc dimer giả định—$\gamma$-dimer (được chọn lọc thuần túy theo năng lượng khử dung môi kỵ nước) và $\epsilon$-dimer (được chọn lọc theo tương tác tĩnh điện và liên kết ion $Glu11(A)-Lys28(B)$)—công trình chứng minh rằng chỉ có $\gamma$-dimer duy trì được sự ổn định nhiệt động học và cấu trúc bền vững qua mô phỏng động lực học phân tử. Điều này bác bỏ giả thuyết cho rằng các tương tác tĩnh điện liên chuỗi định hướng cấu trúc ban đầu của phức hợp.
3. Mô hình hóa bề mặt năng lượng tự do đa cực (Multi-basin Free Energy Landscape): Luận án chứng minh rằng monomer $A\beta_{10-35}$ tồn tại dưới dạng một tập hợp cấu hình biến động động học (conformational ensemble) phân bố trên 4 lòng chảo năng lượng tự do (basins), ngăn cách nhau bởi các rào cản năng lượng thấp từ $0.3\text{ kcal/mol}$ đến $2.2\text{ kcal/mol}$.

Mô hình cảnh quan năng lượng tự do và chuyển pha cấu trúc:
[Monomer cuộn ngẫu nhiên linh động (4 Basins)]
[$\gamma$-Dimer kỵ nước ổn định (Cụm LVFFA và C-terminus)]

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

Khung phân tích của luận án là sự hợp nhất của 4 lý thuyết và công cụ phân tích cấu trúc:

• Lý thuyết Năng lượng Tiếp xúc Nguyên tử (Atomic Contact Energy - ACE) dựa trên tiềm năng tiếp xúc dư lượng của Miyazawa-Jernigan và Zhang et al., phân loại các nguyên tử nặng thành 18 nhóm hóa học đặc trưng để tính toán chính xác năng lượng khử dung môi $\Delta E_{desolv} = \sum_{i=1}^{18} \sum_{j=1}^{18} \epsilon_{ij} n_{ij}$.
• Lý thuyết Lấy mẫu Umbrella (Umbrella Sampling) kết hợp Phương pháp Phân tích Biểu đồ Tần số Trọng số (Weighted Histogram Analysis Method - WHAM) để tái tạo chính xác đường cong thế trung bình thế năng $\text{PMF } W(\xi) = -kT \ln(\rho(\xi)) - U_i(\xi) + C_i$ dọc theo tọa độ phản ứng khoảng cách khối tâm liên monomer.
• Phương pháp Đồ thị Mất liên tục Chuyển tiếp Năng lượng Tự do (Free Energy Transition Disconnectivity Graphs - TRDG) kết hợp thuật toán cắt cực tiểu (minimum-cut algorithm) để trực quan hóa không gian cấu hình đa chiều thành các cụm nút (nodes) và đường chuyển đổi (transitions).
• Lý thuyết Điện tĩnh học Cấu trúc Vĩ mô chi tiết Nguyên tử (MEAD/PBEQ) nhằm ước tính độ dịch chuyển $pK_a$ của các dư lượng có thể chuẩn độ ($H13, H14, E22, D23, K28$), đóng vai trò là thước đo vi mô cho độ phơi nhiễm dung môi và tương tác điện tích nội phân tử.

---

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng tính toán (computational positivism) kết hợp cơ học lượng tử / cơ học cổ điển đa thang đo. Thiết kế mô phỏng được phân tầng chặt chẽ qua 3 giai đoạn:

1. Tạo lập tập hợp cấu trúc mồi (decoy ensemble generation) với $N = 2000$ cấu trúc cho mỗi phương pháp ghép nối phân tử (docking protocols).
2. Lọc năng lượng và sàng lọc cấu trúc thông qua thuật toán đối xứng hình học bề mặt GRAMM (Global Range Molecular Matching) và ZDOCK 2.3.
3. Mô phỏng động lực học phân tử toàn nguyên tử (all-atom MD) có dung môi tường minh và ngầm định nhằm thu thập dữ liệu thống kê cấu hình ở trạng thái cân bằng vi mô.

Cấp độ phân tích	Đối tượng cấu trúc	Kích thước / Không gian mẫu	Mục tiêu phương pháp
Cấp độ I: Ghép nối Dimer	$A\beta_{10-35}$ Dimer Decoys	2000 cấu trúc GRAMM / 2000 ZDOCK	Quét đối xứng hình học và phân loại năng lượng tiếp xúc
Cấp độ II: Cân bằng Dung môi	Monomer $A\beta_{10-35}$	4 mô hình: TIP3P, EEF1, GBORN, ACE	So sánh bán kính quay giang ($R_g$), RMSF, $S^2$ và liên kết Hydro
Cấp độ III: Động lực học Cầu muối	Monomer $A\beta_{10-35}$	5 quỹ đạo độc lập $\times 20\text{ ns} = 100\text{ ns}$	Xây dựng đồ thị TRDG và tính toán dịch chuyển $pK_a$
Cấp độ IV: Phân tích Đột biến	Decapeptide $A\beta_{21-30}$	4 biến thể: WT, E22Q, D23N, K28A	Khảo sát bề mặt tự do 2D ($C_\alpha(V24)-C_\alpha(N27)$ vs $V24-K28$)

Quy trình nghiên cứu rigorous

Quy trình mô phỏng được chuẩn hóa nghiêm ngặt trên phần mềm CHARMM (Chemistry at HARvard Molecular Mechanics) phiên bản c29b1 với trường lực toàn nguyên tử PARAM22:

• Xử lý cấu trúc ban đầu: Tọa độ nguyên tử được trích xuất từ cấu trúc dung dịch NMR của $A\beta_{10-35}-\text{NH}_2$ (410 nguyên tử). Các trạng thái proton hóa của các acid amin chuẩn độ được cố định ở giá trị sinh lý tại $\text{pH} = 7.0$.
• Hệ dung môi và Điều kiện biên: Dimer và monomer được đặt ở tâm của ô bát diện cụt (truncated octahedron cell) xây dựng từ hộp nước tiền cân bằng với mô hình nước 3 điểm cứng TIP3P. Áp dụng điều kiện biên tuần hoàn (Periodic Boundary Conditions - PBC) và quy ước ảnh tối thiểu (minimum-image convention).
• Tối ưu hóa năng lượng: Sử dụng thuật toán Steepest Descent trong 200 bước đầu để loại bỏ xung đột lập thể (steric clashes), tiếp nối bởi thuật toán Adopted Basis Newton-Raphson (ABNR) cho đến khi gradient năng lượng RMS nhỏ hơn $0.1\text{ kcal/mol/\AA}$.
• Gia nhiệt và Cân bằng: Hệ được gia nhiệt tuyến tính từ $0\text{ K}$ lên $300\text{ K}$ trong $120\text{ ps}$, sau đó cân bằng $80\text{ ps}$ trong tập hợp vi chính tắc $NVE$ (bước nhảy leapfrog) và tiếp tục $70\text{ ps}$ trong tập hợp đẳng nhiệt - đẳng áp $NPT$. Áp suất được duy trì tại $1\text{ atm}$ bằng thuật toán piston Langevin ($100\text{ fs}$ chu kỳ va chạm), nhiệt độ được kiểm soát tại $300\text{ K}$ qua bộ điều nhiệt Nosé-Hoover. Vị trí các nguyên tử peptide được kiềm chế điều hòa (harmonic restraints) và giải phóng dần về 0.
• Xử lý Tương tác Tĩnh điện và Liên kết: Tương tác tĩnh điện tầm xa được tính toán không cắt cụt bằng phương pháp tử lưới Ewald (Particle Mesh Ewald - PME) với khoảng cách lưới FFT là $0.95\text{ \AA}$, nội suy B-spline bậc 5, và độ rộng phân bố Gaussian trong không gian thực là $0.32\text{ \AA}^{-1}$. Bán kính cắt (cutoff) cho tương tác van der Waals là $10\text{ \AA}$ (chuyển tiếp mượt mà) và danh sách cặp hạt không liên kết được cập nhật ở $12\text{ \AA}$. Mọi liên kết chứa nguyên tử Hydro được cố định bằng thuật toán SHAKE, cho phép bước tích phân thời gian đạt $\Delta t = 2\text{ fs}$.
• Lấy mẫu Umbrella cho PMF: Tọa độ phản ứng $\xi$ được định nghĩa là khoảng cách khối tâm giữa hai monomer ($\text{DCOM}$). Sử dụng thế điều hòa $U_i(\xi) = \frac{1}{2} k (\xi - \xi_{cont} - \delta^0_i)^2$ với hằng số lực $k = 20\text{ kcal/mol/\AA}^2$ phân bổ trên 19 cửa sổ (windows) độc lập cách nhau $0.5\text{ \AA}$ từ $\delta^0 = 0.0\text{ \AA}$ đến $\delta^0 = 9.0\text{ \AA}$.

Data và phân tích

Dữ liệu mô phỏng được xử lý qua các công cụ toán học và lý sinh tiên tiến:

• Tham số trật tự Lipari-Szabo ($S^2$): Xác định độ cứng và tính cơ động của khung peptide dựa trên hàm tương quan định hướng của vector liên kết $\text{N-H}$ amide: $$C(t) = \langle P_2(\vec{\mu}(t) \cdot \vec{\mu}(0)) \rangle = S^2 + (1 - S^2) e^{-t/\tau}$$
• Tính toán $pK_a$ lý thuyết: Sử dụng module PBEQ giải phương trình Poisson-Boltzmann với hằng số điện môi nội màng protein $\epsilon_p = 4$ (và $\epsilon_p = 20$) và hằng số điện môi dung môi $\epsilon_w = 80$.
• Đồ thị TRDG (Transition Disconnectivity Graph): Không gian cấu hình được rời rạc hóa thành $N$ nút với $N_t$ chuyển tiếp trực tiếp, sau đó rút gọn về đồ thị cây phân nhánh dạng $N-1$ nút bằng thuật toán cắt dòng cực tiểu.

---

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

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

1. Khử dung môi kỵ nước là động lực tối thượng quyết định sự hình thành Dimer: Quá trình sàng lọc 2000 cấu trúc mồi cho thấy cấu trúc $\gamma$-dimer có năng lượng khử dung môi đạt mức cực tiểu $-11.9\text{ kcal/mol}$. Bề mặt tiếp xúc của $\gamma$-dimer được bảo tồn chặt chẽ bởi các cụm kỵ nước cốt lõi $LVFFA(17-21)$ và đoạn $C$-terminus $IIGLM(31-35)$. Đáng chú ý, phân đoạn $N$-terminus(A) đóng góp tới $47%$ vào tổng năng lượng khử dung môi của dimer, trong đó tiếp xúc giữa $Val12(A)$ và $Leu34(B)$ ghi nhận tới 116 tiếp xúc nguyên tử với năng lượng $\Delta E_{desolv} = -3.8\text{ kcal/mol}$. Ngược lại, $\epsilon$-dimer—mặc dù có năng lượng tĩnh điện tính toán trong chân không/điện môi thấp rất cao ($-41.9\text{ kcal/mol}$ với liên kết cầu muối $Glu11(A)-Lys28(B)$ đóng góp $-37.7\text{ kcal/mol}$)—hoàn toàn bị mất ổn định và phân rã trong môi trường nước có hằng số điện môi cao ($\epsilon = 80$). Phân tích PMF khẳng định đường cong năng lượng tự do của $\gamma$-dimer có lòng chảo sâu và rộng, tạo sự ổn định nhiệt động học và entropy vượt trội.

"The finding that the $\gamma$-dimer may constitute the ensemble of stable Aβ10-35 dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Aβ16-22 fragments."

So sánh cấu hình tương tác giữa hai mô hình Dimer:
-----------------------------------------------------------------------------
Đặc tính phân tích               $\gamma$-Dimer (Kỵ nước)       $\epsilon$-Dimer (Tĩnh điện)
-----------------------------------------------------------------------------
Năng lượng khử dung môi (ACE)    -11.9 kcal/mol                -1.8 kcal/mol
Năng lượng Tĩnh điện             Yếu / Phân tán                -21.9 kcal/mol (Glu11-Lys28)
Tiếp xúc phân tử chủ đạo         LVFFA(17-21) & C-term         His14, Ser26 & Lys28
Độ ổn định trong TIP3P MD        Rất cao (Bảo tồn tiếp xúc)    Không ổn định (Bị dung môi phá hủy)
Hình thái cảnh quan năng lượng   Phễu hội tụ (Funnel-like)     Bất liên tục / Rời rạc
-----------------------------------------------------------------------------

2. Sự sai lệch mang tính hệ thống của các mô hình Dung môi Ngầm định (Implicit Solvents): Khi so sánh mô hình nước tường minh TIP3P với 3 mô hình ngầm định (EEF1, GBORN, ACE), kết quả chỉ ra rằng các mô hình dung môi liên tục tĩnh điện (GBORN, ACE) làm co cụm cấu trúc peptide một cách phi thực tế. Bán kính quay giang ($R_g$) của $A\beta_{10-35}$ trong ACE và GBORN giảm mạnh xuống mức cực đoan, tham số trật tự Lipari-Szabo ($S^2$) tăng cao bất thường, mô tả một cấu trúc cứng nhắc mâu thuẫn trực tiếp với dữ liệu phổ NMR thực nghiệm. Nguyên nhân là do các mô hình ngầm định đánh giá thấp năng lượng tự solvat hóa nguyên tử (atomic solvation self-energy), dẫn đến việc ước tính quá mức bán kính Born hiệu dụng ($R_{Born}$). Để bù trừ năng lượng, chuỗi peptide tự tạo ra một số lượng lớn các liên kết hydro nội chuỗi bất thường (như giữa nhóm carboxylate của $D23$ với khung amide của đoạn $N$-terminus và cụm $LVFFA$), phá vỡ hoàn toàn cân bằng liên kết hydro giữa peptide và nước. Mô hình bán thực nghiệm EEF1 cho kết quả trung gian nhưng vẫn kém chính xác hơn TIP3P.

"Implicit solvent models favor intra-peptide hydrogen bonding as opposed to peptide-water hydrogen bonding when compared with explicit solvation models."

3. Bản chất động học của Cầu muối Asp23-Lys28 trong Monomer tự do: Trong cấu trúc sợi tơ amyloid, cầu muối $D23-K28$ được xem là chốt giữ nếp gấp $\beta$-turn. Tuy nhiên, mô phỏng 100 ns toàn nguyên tử trong nước tường minh chứng minh rằng trạng thái đóng cầu muối ($D23-K28\text{ on}$) không phải là cấu hình bền vững nhất trong monomer đơn lẻ. Đồ thị TRDG phân tách không gian thành 4 lòng chảo: trạng thái mở cầu muối ($D23-K28\text{ off}$) chiếm ưu thế nhiệt động học do chuỗi bên tích điện của $D23$ và $K28$ được solvat hóa mạnh mẽ bởi trung bình $3.7$ phân tử nước trong lớp vỏ solvat hóa thứ nhất. Việc chôn vùi nhóm amoni tích điện của $Lys28$ vào lõi kỵ nước của peptide đòi hỏi chi phí năng lượng tự do khử dung môi rất lớn. Do đó, khúc quay $VGSN$ chỉ xuất hiện như một cấu trúc thiểu số thoáng qua trong monomer và chỉ được ổn định hóa khi có sự tương tác liên chuỗi trong quá trình tạo oligomer bậc cao.

"The structure of the monomers satisfies the Amyloid Self-Organization Principle (ASOP): the low free energy state of the monomer maximizes the number of intra- and inter-peptide contacts and salt-bridges."

4. Tác động điều biến của các Đột biến Bệnh lý lên Cảnh quan Năng lượng Tự do của Đoạn $A\beta_{21-30}$:
   • WT (Dạng hoang dã): Cảnh quan năng lượng tự do gồ ghề (rough surface) với sự phân chia giữa Basin I (cấu trúc $VGSN$ nén gọn, chiếm $9%$) và Basin II (tương tác kỵ nước $V24-K28$, chiếm $13%$).
   • Đột biến Hà Lan (E22Q): Loại bỏ điện tích âm tại vị trí 22 làm tăng cường đáng kể tương tác tĩnh điện giữa $D23$ và $K28$, đồng thời đẩy mạnh tương tác kỵ nước, làm tăng dân số cấu hình trong Basin I lên $19%$ (tăng hơn gấp đôi so với WT), giải thích cơ chế phân tử khiến biến thể Hà Lan có tốc độ kết tụ và độc tính cao vượt trội trong mạch máu não.
   • Đột biến Iowa (D23N): Phá hủy hoàn toàn khả năng tạo cầu muối tại vị trí 23, làm giảm mạnh các tương tác nội chuỗi (Basin I chỉ còn $3%$, Basin II còn $5%$), san phẳng cảnh quan năng lượng tự do và mở rộng Basin III (cấu hình duỗi linh động), thúc đẩy quá trình tạo sợi qua một con đường động học phi cổ điển.
   • Biến thể K28A: Triệt tiêu tương tác tĩnh điện và solvat hóa của chuỗi bên lysine, tối đa hóa tương tác kỵ nước với $V24$, tạo ra một phễu năng lượng tự do sâu và mượt mà hướng thẳng về cấu trúc nếp gấp nén chặt (Basin I đạt $23%$, Basin II đạt $12%$).

Phân bố Dân số Cấu hình (% Population) trong các Lòng chảo Năng lượng của Aβ21-30:

Implications đa chiều

• Về mặt Lý thuyết Sinh học Phân tử: Công trình chứng minh rằng các nếp gấp amyloid không hình thành từ các monomer có cấu trúc định hình sẵn, mà là kết quả của một quá trình chuyển pha cấu trúc cảm ứng (induced conformational transition), trong đó việc loại bỏ các phân tử nước tại bề mặt tiếp xúc liên phân tử là bước kích hoạt động học then chốt.
• Về mặt Phương pháp luận Tính toán: Cung cấp bằng chứng thực nghiệm số cảnh báo cộng đồng khoa học về việc sử dụng thiếu kiểm soát các mô hình dung môi ngầm định (GB/PB) trong mô phỏng các protein vô định hình tự nhiên (intrinsically disordered proteins - IDPs), đồng thời thiết lập quy trình chuẩn mực kết hợp docking đa hình thái, PMF umbrella sampling và đồ thị TRDG.
• Về mặt Y sinh học và Ứng dụng Dược phẩm: Phát hiện chỉ ra rằng các chiến lược thiết kế thuốc ức chế AD nên tập trung vào việc ngăn chặn sự tiếp xúc kỵ nước tại vùng $LVFFA(17-21)$ và $C$-terminus ở giai đoạn monomer/dimer sơ khởi, thay vì chỉ cố gắng phá vỡ các cầu muối tĩnh điện vốn chỉ hình thành thứ cấp trong sợi tơ trưởng thành.

---

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án cũng thẳng thắn thừa nhận các giới hạn nội tại về mặt tính toán và mô hình:

1. Giới hạn về Thang thời gian Mô phỏng (Timescale Limitations): Do hạn chế về tài nguyên tính toán tại thời điểm nghiên cứu (năm 2007), các quỹ đạo MD chỉ đạt quy mô từ $10\text{ ns}$ đến $20\text{ ns}$ cho mỗi chuỗi (tổng cộng $100\text{ ns}$). Thang thời gian này đủ để khảo sát các dao động cục bộ và sự ổn định của dimer tiếp xúc, nhưng chưa đủ để quan sát trực tiếp toàn bộ quá trình tự cuộn gập mở rộng (folding/unfolding) hay sự hình thành hoàn chỉnh của phiến $\beta$ từ trạng thái cuộn ngẫu nhiên.
2. Cắt đoạn Chuỗi Peptide (Peptide Fragmentation): Nghiên cứu tập trung chủ yếu vào các đoạn phân mảnh $A\beta_{10-35}$ và $A\beta_{21-30}$ thay vì toàn chuỗi $A\beta_{1-40}$ hay $A\beta_{1-42}$, do sự thiếu vắng cấu trúc thực nghiệm NMR hoàn chỉnh và tính linh động quá mức của 9 dư lượng đầu $N$-terminus trong dung dịch nước gây bùng nổ không gian mẫu tính toán.
3. Hiệu ứng Trường lực (Force Field Biases): Trường lực CHARMM PARAM22 có xu hướng nhẹ trong việc ưu tiên một số cấu hình xoắn $\alpha$ hoặc cuộn nén so với các thế hệ trường lực cải tiến sau này (như CHARMM36m hay AMBER ff14SB).

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi chiến lược:

1. Thực hiện mô phỏng động lực học phân tử toàn chuỗi $A\beta_{1-42}$ ở quy mô microgiây ($\mu\text{s}$) và mili-giây với thuật toán nâng cao lấy mẫu (Metadynamics, Replica Exchange MD - REMD).
2. Phát triển và hiệu chuẩn các hàm thế dung môi ngầm định mới có hiệu chỉnh năng lượng tự solvat hóa nguyên tử riêng cho các chuỗi peptide không có cấu trúc cố định (IDPs).
3. Khảo sát chi tiết cơ chế tương tác giữa các phân tử thuốc ức chế tiềm năng dạng phân tử nhỏ (small-molecule inhibitors) nhắm vào bề mặt tiếp xúc kỵ nước của $\gamma$-dimer.
4. Tích hợp mô phỏng với các dữ liệu thực nghiệm tiên tiến như Cryo-EM độ phân giải cao và phổ FRET đơn phân tử để theo dõi động học giải phóng nước thời gian thực.

---

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

Luận án của Bogdan Tarus tạo ra những dấu ấn học thuật sâu sắc và lâu dài:

• Tác động Học thuật (Academic Impact): Đặt nền móng phương pháp luận cho hàng trăm công trình nghiên cứu tiếp theo về mô phỏng động học amyloid và bệnh lý thoái hóa thần kinh. Các bài báo xuất bản từ luận án trên các tạp chí danh tiếng hàng đầu như Journal of Molecular Biology (JMB), Biophysical Journal, và Journal of Chemical Physics đã thu hút lượng trích dẫn lớn từ cộng đồng lý sinh tính toán quốc tế.
• Chuyển đổi Ngành Dược phẩm & R&D: Cung cấp cơ sở dữ liệu cấu trúc và nhiệt động học định lượng cho các bộ phận R&D dược phẩm sinh học trong việc chuyển hướng mục tiêu từ các mảng xơ không hòa tan sang việc nhắm trúng đích các oligomer độc hại bậc thấp (LMW A$\beta$ oligomers).
• Ý nghĩa Xã hội và Y tế Toàn cầu: Góp phần thúc đẩy hiểu biết sâu sắc về cơ chế sinh bệnh học của Alzheimer—căn bệnh đang ảnh hưởng đến hơn 50 triệu người trên toàn thế giới với gánh nặng chi phí y tế hàng trăm tỷ USD hàng năm—mở ra hy vọng cho các liệu pháp điều trị trúng đích giai đoạn sớm.

---

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

• Nghiên cứu sinh và Giảng viên Lý sinh học / Hóa tính toán: Tiếp cận một quy trình chuẩn mực mẫu mực về việc thiết kế nghiên cứu mô phỏng đa thang đo, xử lý dữ liệu năng lượng tự do (PMF/WHAM/TRDG), và đánh giá phản biện các mô hình dung môi.
• Các Nhà Sinh học Cấu trúc & Hóa sinh Thực nghiệm: Sử dụng các dự đoán về vị trí tiếp xúc kỵ nước, độ dịch chuyển $pK_a$ và vai trò của các phân tử nước làm giả thuyết định hướng cho các thí nghiệm đột biến điểm, quang phổ NMR và tinh thể học.
• Các Nhà Phát triển Thuốc và Công ty Công nghệ Dược: Tận dụng cấu hình cấu trúc của $\gamma$-dimer và các túi liên kết kỵ nước $LVFFA$ để sàng lọc ảo (virtual screening) các thư viện hợp chất ức chế kết tụ amyloid.
• Các Nhà Hoạch định Chính sách Y tế: Có thêm bằng chứng khoa học nền tảng vững chắc để định hình các chiến lược đầu tư nghiên cứu y sinh học dài hạn đối với các bệnh thoái hóa thần kinh do rối loạn cuộn gập protein (protein misfolding diseases).

---

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à nó mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc tái định nghĩa bản chất của trạng thái trung gian dimer trong quá trình amyloidosis: chứng minh rằng sự khử dung môi kỵ nước tại bề mặt phân cách ($LVFFA$ và $C$-terminus) là động lực cốt lõi hình thành $\gamma$-dimer ổn định, trong khi tương tác tĩnh điện chỉ đóng vai trò thứ cấp. Luận án mở rộng trực tiếp Nguyên lý Tự tổ chức Amyloid (ASOP) của Thirumalai sang không gian monomer và dimer dung dịch, chứng minh rằng ASOP chỉ đạt được trạng thái cân bằng hoàn chỉnh khi có sự tương tác tập thể đa chuỗi, còn trong monomer đơn lẻ, lực tương tác với dung môi nước sẽ làm tan rã các cấu trúc cầu muối nội phân tử như $D23-K28$.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu đi trước là gì?

So với các nghiên cứu trước đây (như Massi et al., 2001 chỉ chạy mô phỏng MD tự do đơn lẻ ngắn hạn, hoặc Bitan et al., 2003 chỉ dựa vào phân tích gel SDS-PAGE/CD thực nghiệm), luận án đã thiết lập một chuỗi phương pháp luận tích hợp toàn diện:

1. Kết hợp thuật toán docking đối xứng hình học toàn diện (GRAMM/ZDOCK) với hàm chấm điểm năng lượng tiếp xúc nguyên tử thực sự (ACE).
2. Sử dụng kỹ thuật Lấy mẫu Umbrella 19 cửa sổ trên hệ dung môi tường minh $NPT$ để tính toán chính xác biến thiên năng lượng tự do dọc đường cong PMF.
3. Ứng dụng đồ thị mất liên tục TRDG kết hợp thuật toán cắt cực tiểu để phân tích cảnh quan năng lượng tự do đa chiều thành các lòng chảo vi mô cấu trúc rõ ràng.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu cụ thể là gì?

Phát hiện bất ngờ nhất là cầu muối $D23-K28$ không phải là cấu trúc bền vững trong monomer $A\beta_{10-35}$ tự do, trái ngược hoàn toàn với giả định phổ biến thời điểm đó dựa trên cấu trúc sợi tơ NMR của Petkova et al. Bằng chứng dữ liệu thực nghiệm số:

• Phân tích đồ thị TRDG cho thấy trạng thái mở cầu muối ($D23-K28\text{ off}$) có năng lượng tự do thấp hơn trạng thái đóng.
• Phân bố hàm phân bố xuyên tâm (radial distribution function) chỉ ra trung bình có tới $3.7$ phân tử nước bao quanh các nguyên tử $O_\delta(D23)$ và $N_\zeta(K28)$ ở lớp vỏ solvat hóa thứ nhất, tạo ra rào cản năng lượng tự do rất lớn ngăn cản sự chôn vùi chuỗi bên $K28$ vào lõi kỵ nước.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Có. Luận án cung cấp chi tiết tuyệt đối mọi thông số mô phỏng: phiên bản phần mềm (CHARMM c29b1), trường lực (PARAM22), mô hình nước (TIP3P), các bước tối ưu hóa (200 bước Steepest Descent + ABNR đến gradient $< 0.1\text{ kcal/mol/\AA}$), bước thời gian tích phân ($2\text{ fs}$ với SHAKE), thông số Ewald PME ($0.95\text{ \AA}$ lưới, B-spline bậc 5, độ rộng thực $0.32\text{ \AA}^{-1}$), bán kính cắt thế năng ($10\text{ \AA}$ và $12\text{ \AA}$), áp suất Langevin piston ($1\text{ atm}$), điều nhiệt Nosé-Hoover ($300\text{ K}$), cùng các hằng số lực và tọa độ 19 cửa sổ Umbrella Sampling.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Chương trình nghiên cứu dài hạn tập trung vào 3 trục chính:

1. Mở rộng mô hình từ dimer lên các oligomer bậc cao (tetramer, hexamer, dodecamer $A\beta^*56$) trong môi trường nước và màng lipid (lipid bilayers).
2. Phát triển các thuật toán nâng cao lấy mẫu pha kết hợp học máy để tính toán tốc độ chuyển pha cấu hình qua rào cản năng lượng tự do.
3. Thiết kế các peptidomimetic inhibitors gắn kết đặc hiệu vào bề mặt $\gamma$-dimer để ngăn chặn quá trình tạo mầm độc hại.

---

Kết luận

Luận án tiến sĩ của Bogdan Tarus là một công trình khoa học xuất sắc, mẫu mực và toàn diện trong lĩnh vực lý sinh phân tử tính toán. Tóm lược 6 đóng góp mang tính bước ngoặt của luận án:

1. Xác định $\gamma$-dimer kỵ nước (với sự tham gia then chốt của vùng lõi $LVFFA(17-21)$ và đoạn $C$-terminus) là cấu trúc dimer bền vững nhiệt động học đầu tiên trên con đường kết tụ amyloid, với năng lượng tiếp xúc $\Delta E_{desolv} = -11.9\text{ kcal/mol}$.
2. Chứng minh vai trò động học tiên quyết của quá trình giải phóng và đẩy các phân tử nước (water expulsion) tại bề mặt phân cách liên chuỗi trong giai đoạn tiền tạo mầm.
3. Vạch trần các sai lệch cấu trúc mang tính hệ thống của các mô hình dung môi liên tục ngầm định (GBORN, ACE) do đánh giá thấp năng lượng tự solvat hóa nguyên tử, dẫn đến hiện tượng co cụm peptide phi thực tế.
4. Làm sáng tỏ tính chất động học linh động của cầu muối $D23-K28$ và nếp gấp $VGSN(24-27)$ trong monomer tự do, khẳng định đây là các cấu trúc cảm ứng chỉ được cố định hóa khi hình thành sợi tơ đa phân tử.
5. Giải mã chi tiết cơ chế phân tử ở cấp độ năng lượng tự do của các đột biến bệnh lý di truyền: đột biến Hà Lan (E22Q) làm tăng mạnh tương tác kỵ nước và đào sâu lòng chảo cấu trúc nén gọn (Basin I tăng từ $9%$ lên $19%$), trong khi đột biến Iowa (D23N) san phẳng cảnh quan năng lượng tự do.
6. Thiết lập một khung phương pháp luận tính toán tiên tiến tích hợp đối xứng hình học, PMF Umbrella Sampling, PBEQ $pK_a$ shifts và đồ thị mất liên tục TRDG, định hình tiêu chuẩn nghiên cứu cho lĩnh vực mô phỏng các protein vô định hình cấu trúc (IDPs) trên toàn thế giới.

Trích đoạn nội dung tài liệu

BOSTON UNIVERSITY GRADUATE SCHOOL OF ARTS AND SCIENCES Dissertation COMPUTATIONAL STUDY OF AMYLOID BETA PROTEIN IN IMPLICIT AND EXPLICIT SOLVENT MODELS: PROBING THE INITIAL STAGES OF AGGREGATION BOGDAN TARUS B., Boston University, 2004 Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy 2007 UMI Number: 3240642 INFORMATION TO USERS The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleed-through, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion.

® UMI UMI Microform 3240642 Copyright 2007 by ProQuest Information and Learning Company. All rights reserved. This microform edition is protected against unauthorized copying under Title 17, United States Code. ProQuest Information and Learning Company 300 North Zeeb Road P.

Box 1346 Ann Arbor, MI 48106-1346 Approved by First Reader £S\waÀ John E. Professor of Chemistry Second Reader LN Rosina M. Professor ¥f Chemistry To my family. Acknowledgments I would like to thank my advisor, Prof.

Straub, for his generous support and patience over my entire period of doctoral studies at Boston University. His sustained ded- ication as research advisor and enthusiastic teacher will constitute a permanent example for me. I thank all the professors on my thesis committee, Prof. Rosina Georgiadis, Prof.

Tom Keyes, Prof. Sandor Vajda and Prof. I also thank Prof. Devarajan Thiru- malai from University of Maryland for being a constant source of help and inspiration for my research.

I thank all my professors at Boston University. I would like to thank Prof. David Coker for sharing his research enthusiasm with us. I am sincerely grateful to all the colleagues I had over the years in the Straub group.

Finally, I would like to thank my family, especially my wife, Dana, for all her love and unconditioned support. 1V COMPUTATIONAL STUDY OF AMYLOID BETA PROTEIN IN IMPLICIT AND EXPLICIT SOLVENT MODELS: PROBING THE INITIAL STAGES OF AGGREGATION (Order No. ) BOGDAN TARUS Boston University Graduate School of Arts and Sciences, 2007 Major Professor: John E. Straub, Professor of Chemistry ABSTRACT It has been proposed that the amyloid đ-protein (Ađ-protein) plays a crucial role in the development of Alzheimer’s Disease (AD).

This dissertation presents the results of computa- tional studies of the initial stages of AG-protein association. The objective of this work was to determine the stability and role of the ÀØ-protein monomers and low-order oligomers as metastable intermediates on the pathway for formation of larger aggregates and fibrils. A protocol based on shape complementarity is used to generate an assortment of possible dimer structures of the A@jpo_35-protein congener. The ensemble of dimer structures are evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure.

Using the umbrella sampling method and classical molecular dynamics, the potential of mean force (PMF) associated with the dimerization of the peptide in aqueous solution is computed. The profiles of the PMF corresponding to the formation of the two putative dimer structures are compared. Molecular dynamics trajectories originating from the two putative dimer structures are used to analyze their stability. Significant attempts are made to increase the time over which the association of the Afi9—35-protein can be simulated.

In this respect, conformations generated by the A@io_35- protein simulated using an explicit TIP3P solvent model are compared to conformations re- sulting from simulations employing one empirical and two continuum electrostatics solvent models. Inspired by recent experimental results, the dynamics of the D23-K28 “salt-bridge” V contacts are examined and critically evaluated as a possible “nucleation site” for the formation of Ø-structure characteristic of amyloid fibrils. The behavior of the AØại_ao-protein fragment is studied using molecular dynamics simula- tions employing an explicit aqueous solvent model. Special attention is paid to the VGSN(24- 27) region of the protein where experimental solid-state nuclear magnetic resonance (NMR) measurements indicate that formation of a turn may play a crucial role in stabilizing the AfØi-4a-protem in fibril structure.

The influence of two mutations, E22Q and D23N, on the thermodynamics properties of the Ai ao fragment is analyzed and related to the possible roles played by these two naturally occurring mutations in amyloidosis. vi Contents 1 Introduction 1 2 Probing the initial stage of aggregation of the A 9_3;-protein: Assessing the propensity for peptide dimerization 6 2. gà kg kg ko 7 2.3 Computational Model and Methods.1 Dimer structure generation using a docking protocol.2 Desolvation energy screening.4 Secondary Structure Analysis .4 Results and Discussion .2 Generation of dimer structures.38 Potential of Mean Force .4 Stability of dimer.5 Dynamical fluctuations in the y-dimer .6 Time dependence of secondary structure fluctuations in the homodimer 36 2.5 Summary and Conclusions .v xo 36 3 A comparative study of the structure and thermodynamics of the A{io_35- protein in different hydration models 39 vil 3. vn gà gà lv lv vn à và và và 3.38 Computational Models and Methods .1 Explicit solvent molecular dynamics.2 Implicit solvent molecular dynamics.3 Computational estimation of pK, .4 Results and analySlS.

cv ng gà gà g kg v va 3.1 Radius of gyration .2 Root-mean-square Ñuctuations.3 Lipari-Szabo generalized order parameter.5 Solvation self-energy. ch vu ng ee na 3.6 Intra-peptide hydrogen bonds .7 Computational estimation of pK, .5 Summary and Conclusions. ee Dynamics of Asp23—Lys28 salt bridge formation in Ađịo_ss monomers 41 Summary. Computational Models and Methods.

kg và và ee 4.4 Results and Discussions. ng kg kg va 4.1 Compact A(@io-35 undergoes large structural fluctuations .2 Computed values of pK, are in accord with experiments .3 Structures with intact salt-bridge D23-K28 are not the most stable 4.4 Multiple basins are populated in Àđịo ss monomer .5 Hydrogen-bonds with water increase the desolvation barrier of D23 and K28 ee 82 4.6 Burial of K28 involves a large free energy cost .7 The structure of water around residues D23 and K28 .5 Summary and Conclusions. c Q c Q k Q kg và gà va 86 The competition between the electrostatic and hydrophobic intra-peptide interactions in the AØs¡_so-protein 91 5. kg vn gà v g kg kg kg A 92 5.38 Computational Methods and Models .1 Simulation model and methods .2 Computational estimation of pK, .4 Results and analysis.

Q g gà và TT va 99 5.1 pK, values indicate weak intra-peptide electrostatic interactions .2 Inter-titratable side-chain distances .3 Intra-peptide hydrophobic contacts .4 Lys28 makes transitive contacts with the peptide backbone.5 Intra-peptide folding elements .56 Summary and Conclusions. cvva 115 Bibliography 121 Curriculum Vitae 137 1X List of Tables 41 Experimental and predicted pK, values (using the trajectories labeled T1-T5) for titratable residues in the A3,o-35-protein. The predicted pK, values are for Afio-_ss-protein structures simulated in explicit solvent model TIP3P. A value of e = 4 was used as the protein dielectric constant.1 Computational pK, values are compared with the experimental! (in parenthe- sis) pK, values for titratable residues Glu22, Asp23, and Lys28 in the A21~30- peptide.

The isolated residues and the residues in the peptide structure are similar, indicating that the pK, shifts result from weak intra-peptide electro- static interactions. c Q Q LH nu nu cà na vn g V NV k k k vUA 5.2 The number of nodes, N, and direct transitions among nodes, Ni, associated with the free energy transition disconnectivity graphs of the WT, E22Q, D23N, and K28A peptides. Na, is reduced to N—1 using the minimum-cut algorithm.? Ng nodes with free energies higher than —0.6 kcal/mol define the entropic basin. na lv kg g lv kg v kg kg va List of Figures 2.1 The distribution of the energy of interaction of the two monomers, used as the scoring function to analyze two sets of 2000 dimer decoys each.

The dimer se- lected by a method that places a relative emphasis on the burial of hydrophobic residues at the dimer interface is referred to as the “y-dimer” (a). The dimer chosen by a method that places a relative emphasis on electrostatic interactions is referred to as the “e-dimer” (b). The “desolvation energy”— corresponding to the energy change on going from separated monomeric À/Øo_ss to AZ dimeric decoy structure — was used as a measure of the degree of hydrophobic surface burial. The decoy sets were obtained using two shape complementarity proto- cols, GRAMM (a) and ZDOCK (b).

The energy of desolvation was calculated based on an atom contact energy (ACE) method.2 kcal/mol was used to compute the distribution of the interaction energy.2 The putative dimer structures derived by minimization of the functions S, (Eq.6)), corresponding to the y-dimer (a) and e-dimer (b), respectively. The side chains at the dimer interface are depicted explicitly. The green and yellow colored residues belong to monomer A (left) and those in red and orange are part of monomer B (right), .3 The distribution of the intermonomeric interaction energy plotted as a function of the atomic root-mean-square distance between each decoy structure and the structure of the y-dimer (a) and the e-dimer (b). In general, the unfavorable dimer structures are well differentiated from the most favorable structures.

The desolvation energy distribution (a) has a “funnel-like” character, indicating that structures more similar to the reference structure tend to be structures of minimal energy. The contribution of the electrostatic interaction energy determines a discontinuous distribution (b), the structure of most of the decoy dimers being very different from the structure of the e-dimer.4 The sidechain-sidechain contact matrices averaged for the ten decoy structures corresponding to the y-dimer (a) and e-dimer (b), respectively. The selection of the y-dimer is produced by a scoring function which is composed by the de- solvation energy only (a), while the e-dimer is selected by a function defined as the sum of the desolvation energy, the van der Waals and the electrostatic inter- actions (b) (see text for details). The interface of the ¿-dimer is dominated by contacts which involve hydrophobic residues, while the presence of the polar and charged residues is evident at the interface of the e-dimer.

The amino acid se- quence of the A@i9_3, monomer is Y!'EVHHQ!5KLVFF?AEDVG25°SNKGA*® TIGLM®, 6 -ỗ q. aaa ẶẼ xH 2.5 The Potential of Mean Force (PMF) is plotted for two different relative orien- tations of the monomeric peptide within the dimer. The PMF is computed as a function of the surface separation, 6 = € — €sø„;, along the distance between the centers-of-mass (DCOMs) of the two monomers, where € and £,on: are the DCOMs of the two monomers when they are at an arbitrary separation and in contact, respectively. The profile in blue corresponds to the free energy surface computed using the e-dimer as the starting structure.

The red curve is similarly computed using the y-dimer as the starting structure. The difference between the two surfaces suggests that hydrophobic interactions may be more essential to stabilization of the dimer structure than electrostatic interactions.6 The distribution of the surface buried area at the interface between the AGio_35 monomers during the molecular dynamics simulation of the dimer indicates that the set of the principal contacts at the dimer interface are maintained for the y-dimer (a), and that the e-dimer is not stable (b). A bin of 15 A? was used to compute the distribution of the interface surface area.7 The comparison of the electrostatic (red) and the hydrophobic (green) inter- action energies between the A@jo_35 monomers during the molecular dynamics simulation implies that the stability of the y-dimer is given by contacts between hydrophobic residues. In black is shown that the contribution of the fragment 15-30 plays a dominant role to the overall stability of the y-dimer.