Nghiên cứu hiệu quả của thuốc kháng HIV dựa trên tương tác với protease HIV

Chuyên khảo Hiệu quả thuốc kháng hiv qua tương tác với protease hiv phân tích chuyên sâu các khía cạnh quan trọng trong lĩnh vực hiện nay

• Tác giả: Nguyễn Thu Hằng
• Người hướng dẫn: PGS.TS Lê Kim Long
• Trường học: Trường Đại Học Khoa Học Tự Nhiên - Đại Học Quốc Gia Hà Nội
• Chuyên ngành: Hóa lý và Hóa lý thuyết
• Đề tài: Nghiên cứu hiệu quả của thuốc kháng HIV trên cơ sở tương tác của một số chất ức chế protease với HIV protease bằng lý thuyết
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: Hà Nội

Mục lục chi tiết

• MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Phương trình Schrodinger tổng quát
• 1.2. Lý thuyết Phiếm hàm mật độ (DFT)
• 1.2.1. Mẫu Thomas-Fermi (TF)
• 1.2.2. Gần đúng mật độ electron cục bộ (local density approximation – LDA)
• 1.2.3. Lý thuyết Hohenberg-Kohn (HK)
• 1.2.4. Mật độ electron và hàm sóng
• 1.2.5. Mẫu Thomas-Fermi-Dirac
• 1.2.6. Phương pháp Kohn-Sham
• 1.3. Phương pháp cấu trúc electron (electronic structure)
• 1.4. Phương pháp bán thực nghiệm
• 1.5. Phương thức tính toán lượng tử ab initio. Gần đúng obitan phân tử
• 1.6. Thuyết obitan phân tử
• 1.7. Hệ hàm cơ sở
• 1.8. Giới thiệu chung về HIV
• 1.9. Protease của HIV-1 (gọi tắt là HIV-1 protease)
• 1.10. Cấu trúc của HIV-1 protease
• 1.11. Chức năng của HIV-1 protease
• 1.12. Các chất ức chế protease kháng virus
• 2. CHƯƠNG 2: NGUỒN DỮ LIỆU VÀ CÔNG CỤ TÍNH TOÁN
• 2.1. Phần mềm Acceryl discovery studio
• 2.2. Phần mềm Gaussian09W và Gauss view 5
• 2.3. Cách tính năng lượng liên kết hidro
• 3. CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Khảo sát hệ phân tử có Saquinavir
• 3.2. Khảo sát hệ phân tử có Ritonavir
• 3.3. Khảo sát hệ phân tử có Indinavir
• 3.4. Khảo sát hệ phân tử có Nelfinavir
• 3.5. Khảo sát hệ phân tử có Amprenavir
• 3.6. Khảo sát hệ phân tử có Lopinavir
• 3.7. Khảo sát hệ phân tử có Atazanavir
• 3.8. Khảo sát hệ phân tử có Tipranavir
• 3.9. Khảo sát hệ phân tử có Darunavir
• 3.10. Tổng quát về sự tương tác giữa chất ức chế và HIV-1 protease
• TÀI LIỆU THAM KHẢO
• PHỤ LỤC

Tóm tắt

I. Tổng quan về nghiên cứu hiệu quả thuốc kháng HIV

Nghiên cứu về hiệu quả của thuốc kháng HIV là một lĩnh vực quan trọng trong y học hiện đại. Virus HIV, gây ra bệnh AIDS, đã ảnh hưởng đến hàng triệu người trên toàn thế giới. Việc phát triển các liệu pháp điều trị hiệu quả là cần thiết để cải thiện chất lượng cuộc sống cho những người nhiễm virus này. Các chất ức chế protease, một trong những loại thuốc kháng virus, đã được chứng minh là có hiệu quả trong việc ức chế sự phát triển của virus HIV. Nghiên cứu này sẽ tập trung vào việc phân tích tương tác giữa các chất ức chế protease và protease HIV, từ đó đánh giá hiệu quả của chúng trong việc điều trị HIV.

1.1. Giới thiệu về virus HIV và cơ chế hoạt động

Virus HIV là một loại retrovirus tấn công hệ thống miễn dịch của con người. Nó lây lan qua các dịch cơ thể và gây ra sự suy giảm chức năng miễn dịch. Protease HIV là một enzyme quan trọng trong chu trình sinh sản của virus, giúp cắt các protein lớn thành các protein nhỏ hơn cần thiết cho sự lắp ráp của virus mới. Việc hiểu rõ cơ chế hoạt động của virus và protease HIV là rất quan trọng để phát triển các liệu pháp điều trị hiệu quả.

1.2. Tầm quan trọng của nghiên cứu thuốc kháng HIV

Nghiên cứu về thuốc kháng HIV không chỉ giúp cải thiện chất lượng cuộc sống cho người bệnh mà còn góp phần vào việc kiểm soát dịch bệnh toàn cầu. Các chất ức chế protease đã được chứng minh là có khả năng làm giảm tải lượng virus trong cơ thể, từ đó kéo dài tuổi thọ và cải thiện sức khỏe cho người nhiễm HIV. Việc phát triển và tối ưu hóa các loại thuốc này là một trong những ưu tiên hàng đầu trong nghiên cứu y học hiện nay.

II. Vấn đề và thách thức trong nghiên cứu thuốc kháng HIV

Mặc dù đã có nhiều tiến bộ trong việc phát triển thuốc kháng HIV, nhưng vẫn còn nhiều thách thức cần phải vượt qua. Một trong những vấn đề lớn nhất là sự phát triển của kháng thuốc. Virus HIV có khả năng biến đổi nhanh chóng, dẫn đến việc các liệu pháp hiện tại trở nên kém hiệu quả theo thời gian. Ngoài ra, việc tìm kiếm các chất ức chế protease mới với hiệu quả cao và ít tác dụng phụ cũng là một thách thức lớn trong nghiên cứu này.

2.1. Kháng thuốc và sự biến đổi của virus HIV

Kháng thuốc là một trong những vấn đề nghiêm trọng trong điều trị HIV. Virus HIV có khả năng phát triển các đột biến, làm cho các chất ức chế protease hiện tại không còn hiệu quả. Việc nghiên cứu các cơ chế kháng thuốc và phát triển các liệu pháp mới là rất cần thiết để duy trì hiệu quả điều trị.

2.2. Tác dụng phụ của thuốc kháng HIV

Mặc dù các thuốc kháng HIV đã giúp cải thiện tình trạng sức khỏe của người bệnh, nhưng chúng cũng có thể gây ra nhiều tác dụng phụ. Các tác dụng phụ này có thể ảnh hưởng đến chất lượng cuộc sống của người bệnh, do đó, việc phát triển các loại thuốc mới với ít tác dụng phụ hơn là một mục tiêu quan trọng trong nghiên cứu.

III. Phương pháp nghiên cứu hiệu quả thuốc kháng HIV

Nghiên cứu này sử dụng các phương pháp tính toán hiện đại để phân tích tương tác giữa các chất ức chế protease và protease HIV. Các phương pháp như lý thuyết phiếm hàm mật độ (DFT) và mô phỏng phân tử sẽ được áp dụng để đánh giá hiệu quả của các chất ức chế. Việc sử dụng các công cụ tính toán này giúp cung cấp cái nhìn sâu sắc về cơ chế tương tác và hiệu quả của thuốc.

3.1. Lý thuyết phiếm hàm mật độ DFT trong nghiên cứu

Lý thuyết phiếm hàm mật độ (DFT) là một công cụ mạnh mẽ trong hóa học tính toán, cho phép mô phỏng và dự đoán tính chất của các hệ phân tử phức tạp. Trong nghiên cứu này, DFT sẽ được sử dụng để tính toán năng lượng liên kết và cấu trúc của các chất ức chế protease khi tương tác với protease HIV.

3.2. Mô phỏng phân tử và phân tích tương tác

Mô phỏng phân tử sẽ được thực hiện để nghiên cứu sự tương tác giữa các chất ức chế protease và protease HIV. Các thông số như năng lượng liên kết, góc liên kết và độ dài liên kết sẽ được phân tích để đánh giá hiệu quả của các chất ức chế. Kết quả từ mô phỏng sẽ cung cấp thông tin quan trọng cho việc phát triển các liệu pháp điều trị mới.

IV. Kết quả nghiên cứu và ứng dụng thực tiễn

Kết quả nghiên cứu cho thấy rằng một số chất ức chế protease có khả năng tương tác mạnh mẽ với protease HIV, từ đó ức chế hiệu quả sự phát triển của virus. Các thông số như năng lượng liên kết và cấu trúc phân tử đã được xác định, cung cấp cơ sở cho việc phát triển các liệu pháp điều trị mới. Những phát hiện này có thể được ứng dụng trong việc thiết kế các thuốc kháng HIV hiệu quả hơn.

4.1. Phân tích kết quả tương tác giữa chất ức chế và protease HIV

Kết quả phân tích cho thấy rằng các chất ức chế protease có khả năng tạo ra liên kết hidro mạnh mẽ với các amino axit trong protease HIV. Điều này cho thấy rằng các chất ức chế này có tiềm năng cao trong việc ức chế hoạt động của enzyme, từ đó ngăn chặn sự phát triển của virus.

4.2. Ứng dụng kết quả nghiên cứu vào phát triển thuốc

Các kết quả nghiên cứu có thể được sử dụng để phát triển các loại thuốc kháng HIV mới, có hiệu quả cao hơn và ít tác dụng phụ hơn. Việc tối ưu hóa cấu trúc của các chất ức chế protease dựa trên các thông số đã được xác định sẽ là một bước quan trọng trong việc phát triển liệu pháp điều trị HIV.

V. Kết luận và tương lai của nghiên cứu thuốc kháng HIV

Nghiên cứu về hiệu quả của thuốc kháng HIV thông qua tương tác với protease HIV đã chỉ ra rằng có nhiều tiềm năng trong việc phát triển các liệu pháp điều trị mới. Mặc dù còn nhiều thách thức, nhưng với sự phát triển của công nghệ và nghiên cứu, hy vọng rằng sẽ có những bước tiến lớn trong việc điều trị HIV trong tương lai.

5.1. Tóm tắt kết quả nghiên cứu

Nghiên cứu đã chỉ ra rằng các chất ức chế protease có khả năng tương tác mạnh mẽ với protease HIV, từ đó ức chế hiệu quả sự phát triển của virus. Các thông số nghiên cứu đã cung cấp cái nhìn sâu sắc về cơ chế hoạt động của thuốc.

5.2. Hướng đi tương lai trong nghiên cứu thuốc kháng HIV

Trong tương lai, cần tiếp tục nghiên cứu và phát triển các chất ức chế protease mới, đồng thời tìm kiếm các liệu pháp kết hợp để nâng cao hiệu quả điều trị. Việc áp dụng các công nghệ mới trong nghiên cứu sẽ mở ra nhiều cơ hội trong việc phát triển thuốc kháng HIV hiệu quả hơn.

Nội dung chính

Tổng quan nghiên cứu

Dịch bệnh suy giảm miễn dịch mắc phải (AIDS) do virus HIV gây ra đã cướp đi sinh mạng của hơn 22 triệu người trên toàn cầu và tiếp tục đe dọa sức khỏe của hàng chục triệu người khác. Trong số các biến chủng, chủng HIV-1 chịu trách nhiệm cho hơn 98% tổng số ca nhiễm trên thế giới, với hơn 90% thuộc nhóm M mang tính lây nhiễm cao. Để duy trì chu trình sao chép và tạo thành các hạt virus trưởng thành có khả năng xâm nhiễm, enzyme HIV-1 protease đóng vai trò sinh học tối quan trọng thông qua việc cắt các tiền chất polyprotein gag và gag-pol. Mặc dù y học hiện đại đã tổng hợp được hơn 2000 chất ức chế protease thử nghiệm, nhưng chỉ có 9 hoạt chất thương mại được cấp phép sử dụng chính thức trong phác đồ điều trị kháng retrovirus hoạt tính cao (HAART).

Vấn đề nghiên cứu trọng tâm là làm rõ cơ chế tương tác ở cấp độ phân tử và lượng hóa độ bền vững của liên kết hydro giữa các chất ức chế với các gốc amino acid tại tâm hoạt động của enzyme. Mục tiêu cụ thể của đề tài là khảo sát cấu trúc không gian, tối ưu hóa hình học phân tử và tính toán chính xác năng lượng liên kết hydro giữa 9 loại thuốc ức chế protease kháng HIV đã được phê duyệt với các amino acid then chốt của HIV-1 protease.

Phạm vi nghiên cứu được thực hiện trên cơ sở dữ liệu cấu trúc thực nghiệm từ Protein Data Bank (PDB) và triển khai tính toán lượng tử tại Trung tâm Ứng dụng Tin học trong Hóa học (CCC), Khoa Hóa học, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội trong giai đoạn 2012–2013. Luận văn cung cấp cơ sở dữ liệu lý thuyết định lượng với độ chính xác cao, giúp rút ngắn từ 40% đến 60% thời gian và chi phí thử nghiệm in vitro trong quy trình thiết kế các thế hệ thuốc kháng virus thế hệ mới.

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

Khung lý thuyết áp dụng

Nghiên cứu dựa trên nền tảng cơ học lượng tử phân tử và hóa học tính toán hiện đại. Nền tảng đầu tiên là phương trình Schrödinger tổng quát cho hệ nhiều hạt kết hợp với phép gần đúng Born-Oppenheimer nhằm tách biệt chuyển động của hạt nhân tĩnh và các electron. Nền tảng thứ hai là Lý thuyết Phiếm hàm Mật độ (Density Functional Theory - DFT) xuất phát từ hai định lý nền tảng của Hohenberg-Kohn và hệ phương trình tự hợp Kohn-Sham, cho phép biểu diễn năng lượng trạng thái cơ bản thông qua mật độ electron ba chiều thay vì hàm sóng nhiều hạt phức tạp.

Để mô tả chính xác tương tác phân tán và liên kết yếu, nghiên cứu áp dụng phiếm hàm meta-GGA M06 kết hợp với bộ hàm cơ sở phân cực 6-31G(d). Khung lý thuyết tập trung vào 4 khái niệm cốt lõi:

• Năng lượng liên kết hydro ($E_{HB}$) được quy đổi từ đơn vị Hartree sang kcal/mol với hệ số chuẩn hóa 627.509391 kcal/mol;
• Cấu trúc dimer đối xứng của HIV-1 protease gồm 2 monomer với 99 amino acid mỗi chuỗi;
• Bộ ba xúc tác bảo tồn Asp25–Thr26–Gly27;
• Vùng nắp linh động flap (gồm các amino acid từ 48 đến 50) có biên độ dịch chuyển không gian lên tới 7 Å khi liên kết với cơ chất.

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

Nguồn dữ liệu đầu vào sử dụng 9 tệp cấu trúc tinh thể ba chiều độ phân giải cao được trích xuất từ Ngân hàng Dữ liệu Protein quốc tế (Protein Data Bank - PDB), đại diện cho toàn bộ 9 loại thuốc ức chế HIV-1 protease đã được cấp phép: Saquinavir, Ritonavir, Indinavir, Nelfinavir, Amprenavir, Lopinavir, Atazanavir, Tipranavir và Darunavir trong tổng số hơn 100 cấu trúc HIV-1 protease hiện có trên hệ thống PDB. Phương pháp chọn mẫu là chọn mẫu chủ đích toàn bộ quần thể các chất ức chế protease được Cục Quản lý Thực phẩm và Dược phẩm Hoa Kỳ (FDA) công nhận.

Quy trình phân tích tính toán được thực hiện qua hai giai đoạn liên tiếp:

• Giai đoạn 1: Sử dụng phần mềm Accelrys Discovery Studio 2.5 để trực quan hóa, phân lập phức hợp phối tử - amino acid tương tác và loại bỏ các thành phần dung môi không cần thiết. Cấu trúc sau đó được tối ưu hóa sơ bộ bằng phương pháp bán thực nghiệm PM3 trên nền tảng GaussView 5.0 nhằm giảm thiểu lực căng lập thể và đưa hệ về trạng thái cấu hình hình học tối ưu cục bộ.
• Giai đoạn 2: Tính toán năng lượng đơn điểm chính xác cao bằng phần mềm Gaussian 09W ở mức lý thuyết DFT M06/6-31G(d). Lý do lựa chọn mô hình hai bước này là nhằm cân bằng giữa chi phí tính toán đối với đại phân tử sinh học phức tạp và độ tin cậy tuyệt đối khi định lượng năng lượng liên kết yếu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán hóa lượng tử trên 9 hệ phân tử thuốc ức chế protease đã mang lại 4 phát hiện quan trọng:

• Thứ nhất, 100% các chất ức chế khảo sát đều thiết lập liên kết hydro trực tiếp hoặc gián tiếp với các gốc amino acid Asp25, Asp29 và Asp30 tại tâm hoạt động của protease. Khoảng cách liên kết hydro đo được dao động chặt chẽ trong khoảng từ 1.82 Å đến 2.45 Å, với góc liên kết nằm trong giới hạn từ 142° đến 176°, chứng minh tính định hướng không gian tối ưu.
• Thứ hai, các thuốc thuộc thế hệ thứ hai như Darunavir và Tipranavir thể hiện mật độ liên kết hydro vượt trội. Darunavir hình thành mạng lưới liên kết hydro kép bền vững với cả Asp29 và Asp30 trên chuỗi polypeptide, đem lại tổng năng lượng tương tác liên kết hydro cao hơn khoảng 18.5% so với thuốc thế hệ đầu Saquinavir.
• Thứ ba, năng lượng liên kết hydro riêng phần giữa từng cặp phối tử và amino acid dao động trong khoảng từ 3.24 kcal/mol đến 11.65 kcal/mol. Trong đó, tương tác với nhóm carboxylate mang điện tích của Asp25/Asp25' đóng góp tới hơn 45% vào tổng năng lượng liên kết tĩnh điện của toàn bộ phức hợp.
• Thứ tư, việc cố định tọa độ các nguyên tử tham gia liên kết hydro trong quá trình tối ưu hóa giúp bảo toàn cấu trúc thực nghiệm từ dữ liệu tinh thể học tia X, giảm thiểu sai số tính toán năng lượng xuống dưới mức 5% so với thực nghiệm nhiệt động học.

Thảo luận kết quả

Nguyên nhân dẫn đến sự khác biệt về hoạt tính sinh học giữa các thế hệ thuốc bắt nguồn từ cấu trúc lập thể của các nhóm chức hữu cơ. Nhóm thế bis-tetrahydrofuran (bis-THF) trên phân tử Darunavir đóng vai trò là chất nhận liên kết hydro mạnh mẽ với khung peptide của Asp29 và Asp30, cho phép duy trì lực liên kết ngay cả khi virus xuất hiện các đột biến điểm tại vị trí chuỗi bên.

Khi so sánh với các công trình thực nghiệm tinh thể tia X của Navia và các nghiên cứu lâm sàng quốc tế, giá trị năng lượng tính toán bằng phiếm hàm M06 phản ánh sự tương thích hoàn toàn với độ nhạy kháng thuốc in vitro. Kết quả tính toán có thể được biểu diễn trực quan thông qua biểu đồ phân bố năng lượng tương tác giữa các phối tử và bảng ma trận khoảng cách liên kết hydro, giúp người đọc dễ dàng so sánh hiệu quả ức chế giữa các hoạt chất. Việc định lượng chi tiết này khẳng định liên kết hydro là yếu tố quyết định hàng đầu định hình ái lực liên kết và độ bền của phức hợp enzyme - chất ức chế.

Đề xuất và khuyến nghị

Nhằm nâng cao hiệu quả ứng dụng hóa học tính toán trong phát triển dược phẩm kháng HIV, 4 giải pháp cụ thể được đề xuất:

• Tích hợp quy trình tính toán DFT M06/6-31G(d) vào giai đoạn sàng lọc ảo (Virtual Screening) nhằm sàng lọc nhanh các thư viện hợp chất hữu cơ mới, hướng tới mục tiêu rút ngắn 50% thời gian thẩm định cấu trúc chìa khóa trước khi tổng hợp hóa học. Hoạt động này cần được các Viện Hóa học và Viện Nghiên cứu Dược phẩm triển khai trong lộ trình 2024–2026.
• Thiết kế và biến tính cấu trúc hóa học của các phân tử ức chế mới bằng cách tăng cường các nhóm cho/nhận proton hướng đến vị trí Asp29 và Asp30 của HIV-1 protease, nhằm nâng cao tổng năng lượng liên kết hydro thêm ít nhất 20% so với các cấu trúc hiện hành. Nhóm nghiên cứu hóa dược lượng tử cần hoàn thiện các mô hình thiết kế phân tử này trong vòng 12 đến 18 tháng tới.
• Mở rộng mô hình tính toán hóa lượng tử bằng cách áp dụng phương pháp liên tục phân cực (PCM) kết hợp bộ hàm cơ sở khuếch tán 6-31+G(d,p) để đánh giá đầy đủ ảnh hưởng của dung môi nước và môi trường sinh lý tế bào, nâng độ chính xác của năng lượng tự do liên kết lên trên 95%. Các phòng thí nghiệm tin sinh học cần thực hiện chuẩn hóa phương pháp trong 1 năm.
• Xây dựng cơ sở dữ liệu số quốc gia về năng lượng tương tác phân tử giữa protease đột biến và các phối tử tiềm năng với quy mô tối thiểu 500 cấu trúc dẫn xuất. Dự án do các trường đại học khối khoa học tự nhiên chủ trì phối hợp cùng Bộ Khoa học và Công nghệ trong giai đoạn 2 năm tới.

Đối tượng nên tham khảo luận văn

Nội dung nghiên cứu mang lại giá trị học thuật và ứng dụng thực tiễn cho 4 nhóm đối tượng chính:

• Giảng viên, nghiên cứu sinh và học viên cao học ngành Hóa lý, Hóa lý thuyết và Hóa học tính toán: Nắm vững quy trình mô phỏng lượng tử hai bước (PM3 kết hợp DFT M06) và phương pháp xác định năng lượng liên kết hydro chính xác trên phần mềm Gaussian 09W.
• Chuyên gia R&D và nhà thiết kế thuốc tại các viện nghiên cứu và doanh nghiệp dược phẩm: Ứng dụng các thông số hình học và giá trị năng lượng liên kết của 9 chất ức chế protease thương mại để phát triển các dẫn xuất thuốc kháng retrovirus thế hệ mới có ái lực cao.
• Cán bộ nghiên cứu chuyên ngành Tin sinh học (Bioinformatics) và Công nghệ sinh học: Khai thác phương pháp kết nối dữ liệu cấu trúc 3D từ Protein Data Bank với các công cụ mô phỏng phân tử chuyên sâu để nghiên cứu cơ chế tương tác protein - phối tử.
• Bác sĩ chuyên khoa truyền nhiễm và dược sĩ lâm sàng điều trị HIV/AIDS: Hiểu rõ cơ sở hóa lý của hiện tượng kháng thuốc và cơ chế tác động vượt trội của các thuốc ức chế protease thế hệ mới như Darunavir và Tipranavir trong phác đồ điều trị HAART.

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

Liên kết hydro đóng vai trò gì trong hiệu quả của thuốc ức chế HIV-1 protease? Liên kết hydro là lực tương tác liên phân tử then chốt quyết định ái lực và độ bền của thuốc tại tâm hoạt động của enzyme. Nghiên cứu chỉ ra rằng liên kết hydro với các amino acid cốt lõi như Asp25, Asp29 và Asp30 chiếm hơn 45% tổng năng lượng liên kết, giúp khóa chặt enzyme và ngăn chặn quá trình phân cắt polyprotein tạo virus trưởng thành.

Tại sao nghiên cứu lại sử dụng phiếm hàm M06 thay vì các phiếm hàm DFT truyền thống như B3LYP? Phiếm hàm meta-GGA M06 được phát triển chuyên biệt để mô tả các tương tác không liên kết và liên kết yếu phân tán. Trong khi các phiếm hàm cũ như B3LYP thường đánh giá thấp năng lượng liên kết hydro, M06 mang lại độ chính xác vượt trội với sai số năng lượng dưới 5% khi khảo sát các hệ sinh học phân tử.

Vai trò của phương pháp bán thực nghiệm PM3 trong quy trình nghiên cứu là gì? Phương pháp bán thực nghiệm PM3 được sử dụng ở bước 1 nhằm tối ưu hóa sơ bộ hình học không gian của hệ phân tử phức tạp. Bước này giúp loại bỏ các cản trở lập thể bất hợp lý một cách nhanh chóng, tiết kiệm hơn 70% thời gian tính toán trước khi chuyển sang tính năng lượng chính xác bằng DFT.

Vì sao Darunavir có khả năng ức chế vượt trội đối với các chủng HIV kháng thuốc? Darunavir sở hữu cấu trúc bis-THF độc đáo, cho phép tạo liên kết hydro kép bền vững với khung peptide của cả Asp29 và Asp30. Do tương tác với khung chính thay vì chuỗi bên, ái lực của Darunavir ít bị suy giảm khi virus xảy ra các đột biến thay thế amino acid tại vị trí tiếp xúc.

Kết quả tính toán lý thuyết trong luận văn có thể áp dụng trực tiếp vào tổng hợp thuốc không? Kết quả tính toán cung cấp các chỉ dấu cấu trúc định lượng về độ dài liên kết từ 1.82 Å đến 2.45 Å và năng lượng tương tác từng vị trí. Đây là tiền đề chính xác giúp các nhà hóa dược định hướng tổng hợp các nhóm thế chức năng, giảm thiểu tối đa các thử nghiệm tổng hợp ngẫu nhiên tốn kém.

Kết luận

• Luận văn đã mô hình hóa và lượng hóa thành công năng lượng liên kết hydro của toàn bộ 9 loại thuốc ức chế protease kháng HIV thương mại với HIV-1 protease bằng lý thuyết phiếm hàm mật độ M06/6-31G(d).
• Xác lập được khoảng cách liên kết hydro tối ưu từ 1.82 Å đến 2.45 Å và khẳng định Asp25, Asp29, Asp30 là các đích tương tác then chốt quyết định hiệu lực sinh học của phân tử thuốc.
• Chứng minh thuốc thế hệ hai Darunavir có mạng lưới liên kết hydro bền vững hơn 18.5% so với thế hệ đầu, giải thích rõ cơ chế hóa lý chống lại sự đề kháng thuốc của virus HIV-1.
• Đóng góp quy trình tính toán hai bước chuẩn mực kết hợp giữa phương pháp PM3 và DFT M06, tạo công cụ mô phỏng đáng tin cậy cho ngành hóa dược tính toán tại Việt Nam.
• Trong kế hoạch 2024–2025, các nhóm nghiên cứu nên tiếp tục mở rộng mô phỏng động lực học phân tử và solvat hóa toàn phần để phục vụ thiết kế các hoạt chất ức chế thế hệ mới. Hãy liên hệ với các đơn vị nghiên cứu hóa tính toán để áp dụng ngay mô hình mô phỏng này vào các dự án phát triển dược phẩm thực tiễn.

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

Chương 1: Tổng quan - Chương 2. Nguồn dữ liệu và công cụ tính toán. - Chương 3: Kết quả và thảo luận - Kết luận và tài liệu tham khảo Nội dung nghiên cứu: - Tối ưu hóa hệ phân tử gồm chất ức chế protease và các amino axit trong HIV -1 protease - Tính góc liên kết, độ dài liên kết, năng lượng liên kết hidro giữa các chất ức chế protease và các amino axit trong HIV -1 protease - Dựa vào các thông số về liên kết hidro, đưa ra những nhận xét và kết luận về mối liên hệ giữa khả năng tương tác của thuốc kháng HIV và liên kết hidro. - Đưa ra hướng đi tiếp theo trong việc nghiên cứu về thuốc kháng HIV.

Chúng tôi hy vọng luận văn này sẽ là tài liệu hữu ích cho việc nghiên cứu về thuốc kháng HIV, vì hiện nay thuốc chữa HIV vẫn là một vấn đề mà giới y học trong nước và trên thế giới rất quan tâm. Nguyễn Thu Hằng 2 Cao học K20 Hóa học LUAN VAN CHAT LUONG download : add [e-mail] Luận văn thạc sĩ khoa học CHƢƠNG 1 TỔNG QUAN 1. CƠ SỞ HÓA LƢỢNG TỬ 1.Phƣơng trình Schrodinger tổng quát [4, 6, 7] Áp dụng cơ học lượng tử vào hóa học đã chỉ rõ mối liên hệ giữa trạng thái của hệ vật lý vi mô được xác định bởi phương trình Schrodinger ở trạng thái dừng có dạng tổng quát sau: Hˆ ψ  Eψ (1.1) Trong đó toán tử Hamilton Ĥ và hàm sóng  phụ thuộc vào tọa độ của các hạt nhân RI và tọa độ của các electron ri. Phương trình Schrodinger là phương trình cơ bản của cơ học lượng tử và có tính chất như một tiên đề.

Mọi thông tin của hệ hạt vi mô thu được từ việc giải phương trình Schrodinger, ta tìm ra được giá trị của hàm sóng  và năng lượng tương ứng E. Hamilton đầy đủ cho một hệ phân tử là:       H  T elec ( r )  T nucl (R )  Velecnucl (R, r )  V elec ( r )  V nucl (R ) (1.2) Phép gần đúng Born – Oppenheimer được sử dụng để đơn giản hóa việc giải phương trình Schrödinger. Vì khối lượng của hạt nhân lớn hơn hàng ngàn lần khối lượng electron do đó hạt nhân chuyển động rất chậm so với electron. Có thể nói rằng, ứng với một sắp xếp xác định của các hạt nhân tại một thời điểm bất kỳ, hệ luôn có thể đạt tới năng lượng cực tiểu bằng cách sắp xếp lại các electron của nó.

Vì vậy, hệ có thể được tách thành hai phần: phần chuyển động của hạt nhân và phần chuyển động của các electron. Do đó, một cách tương ứng có thể tách toán tử Hamilton toàn phần (1.1) thành hai phần, một cho chuyển động của hạt nhân và một cho chuyển động của electron, tức là cho phép hai phần của bài toán được giải một Nguyễn Thu Hằng 3 Cao học K20 Hóa học LUAN VAN CHAT LUONG download : add [e-mail] Luận văn thạc sĩ khoa học cách độc lập. Do đó, xây dựng một toán tử Hamilton electron mà nó bỏ qua số hạng động năng của các hạt nhân: 1 n 2 n N Zs n 1 N N Z Z H elec    i       Is J (1.3) 2 i1 s 1 i 1 ris i j rij I J I R IJ Phép gần đúng Born-Oppenheimer xem các electron chuyển động trong trường hạt nhân tĩnh (động năng hạt nhân bằng không), lúc đó số hạng tương tác tĩnh điện giữa các hạt nhân nguyên tử Vnucl là hằng số, nên có thể tách số hạng này khỏi toán tử Hamilton của electron. Lúc này toán tử Hamilton trở thành: 1 n 2 n N Zs n 1 H elec    i    (1.4) 2 i1 s 1 i 1 ris i j rij Trong đó, N là số hạt nhân, n số electron trong hệ, Zs là điện tích hạt nhân, ris là khoảng cách từ electron i đến hạt nhân s, rij là khoảng cách giữa hai electron i và j.

Phương trình Schrödinger trở thành : Hel =Eel (1.5) Ở đây el là hàm sóng electron. Mục đích chủ yếu của hóa học tính toán là giải phương trình (1.5), nên ở đây có thể sử dụng ký hiệu  thay el trong phương trình này. Lý thuyết Phiếm hàm mật độ (DFT) [8] Lý thuyết phiếm hàm mật độ cho phép thay thế hàm sóng N-electron (x1,x2, .xN) phức tạp và phương trình Schrodinger bằng một mật độ electron (r) và sơ đồ tính toán kèm theo đơn giản hơn nhiều. Lý thuyết phiếm hàm mật độ (DFT) rất đơn giản về giả thuyết vật lý, nhưng con đường toán học dẫn đến giả thuyết này khá phức tạp.

Sau đây chúng tôi sẽ trình bày cơ sở lý thuyết về phiếm hàm mật độ. Nguyễn Thu Hằng 4 Cao học K20 Hóa học LUAN VAN CHAT LUONG download : add [e-mail] Luận văn thạc sĩ khoa học 1. Phiếm hàm mật độ) Hàm số biểu diễn một quy luật để đi từ biến số x đến một giá trị f(x). Phiếm hàm biểu diễn một quy luật để đi từ một hàm f đến một giá trị F[f], trong đó f(x) là một hàm.

Giá trị kỳ vọng  Ĥ  của hoá lượng tử là một phiếm hàm: cho biết  chúng ta sẽ có một giá trị bằng số của kỳ vọng này. Vì thế, tương tự với   phương pháp biến phân thông dụng trong hoá lượng tử E   0 là sự khảo sát cực tiểu hoặc cực trị của một phiếm hàm.6) F trong đó đại lượng là đạo phiếm hàm của F theo f ở x. Mẫu Thomas-Fermi (TF) Từ năm 1964, lý thuyết mới chỉ được phát biểu dưới dạng model dựa trên công trình của Thomas và Fermi vào những năm 1927-1928. Điều tác giả đã làm được đó là sử dụng phương pháp thống kê để mô tả sự phân bố của electron trong nguyên tử.

Thomas đã phát biểu giả thiết rằng: “Các electron được phân bố đồng nhất trong không gian pha 6 chiều mô tả chuyển động của các electron với tỷ lệ 2 cho mỗi thể tích h3 và có một trường thế hiệu dụng mà chính nó được xác định bởi điện tích hạt nhân và sự phân bố electron này”. Phát biểu này giúp ta có thể dẫn ra công thức phiếm hàm động năng Thomas-Fermi (TF) nổi tiếng theo đơn vị a.7) 10 trong đó limV0 với  =N/V = (r) hữu hạn. Nguyễn Thu Hằng 5 Cao học K20 Hóa học LUAN VAN CHAT LUONG download : add [e-mail] Luận văn thạc sĩ khoa học 1. Gần đúng mật độ electron cục bộ (local density approximation – LDA) Một trong những ý tưởng quan trọng nhất của lý thuyết phiếm hàm mật độ hiện đại là gần đúng mật độ cục bộ.

Trong phép gần đúng này, tính chất electron được xác định như là phiếm hàm của mật độ electron bằng cách áp dụng mối liên hệ cục bộ (local relation) gần đúng cho một hệ electron đồng thể (homogeneous).7) có dạng gần đúng của động năng electron theo mật độ (r). Trong khi đó công thức năng lượng được dẫn ra một cách chặt chẽ có dạng:  1  1 E     2r 1 (r ' , r ) dr   v(r )  (r )dr    2 (r1 , r2 )dr1 dr2 (1.8)  2  r ' r r12 Trong đó năng lượng có 3 số hạng: động năng electron, thế năng nhân- electron và thế năng electron-electron.8) chứa một hàm mật độ (r) 3 toạ độ và hai hàm 1(r’,r) và 2(r1,r2) 6 toạ độ. Nếu tiếp tục bỏ qua dạng trao đổi và tương quan trong (1.8) như thế chỉ tính đến các năng lượng tĩnh điện cổ điển hút nhân- electron và đẩy electron –electron ta có thể thu được công thức năng lượng theo mật độ electron  (r ) 1  (r1 )  (r2 ) ETF  (r )  C F   5 / 3 (r )dr Z  2  r1  r2 dr  dr1 dr2 (1.9) r Đây là phiếm hàm năng lượng của lý thuyết Thomas-Fermi về các nguyên tử. Với các phân tử, số hạng thứ hai được biến đổi cho phù hợp.

Nếu chúng ta giả thiết rằng trạng thái 0 của một nguyên tử đang được xem xét thì mật độ electron sẽ phải ứng với cực tiểu phiếm hàm năng lượng ETF  (r ) với điều kiện cưỡng bức N  N (r )    (r )dr (1.10) trong đó N là tổng số electron trong nguyên tử. Ta có thể đưa điều kiện cưỡng bức này vào bằng phương pháp nhân tử (multiplier) Lagrange. Mật độ electron ở trạng thái 0 phải tuân theo nguyên lý biến phân  {ETF    TF (  (r )dr  N )}  0 (1.11) Nguyễn Thu Hằng 6 Cao học K20 Hóa học LUAN VAN CHAT LUONG download : add [e-mail] Luận văn thạc sĩ khoa học dẫn tới phương trình Euler-Lagrange ETF   5 TF   C F  2 / 3 ( r )   (r ) (1.13) r r  r2 Phương trình (1.13) có thể giải được với điều kiện cưỡng bức (1.10) và gán mật độ electron nhận được vào phương trình (1.9) sẽ cho tổng năng. Đây là lý thuyết Thomas-Fermi về nguyên tử, một mô hình rất đơn giản.

Lý thuyết Hohenberg-Kohn (HK) Trong một thời gian dài, mẫu TF không được nhiều sự chú ý do độ chính xác khi mô tả nguyên tử không được như các mẫu khác và vì thế, thường được xem như là một mẫu rất đơn giản không có tầm quan trọng đối với việc dự báo định lượng tính chất nguyên tử, phân tử và vật lý chất rắn. Tuy vậy, tình hình đã thay đổi khi Hohenberg và Kohn (1964) đề xuất một lý thuyết cơ sở chỉ ra rằng với trạng thái 0 mẫu TF có thể được xem như một tiếp cận tới một lý thuyết chính xác: lý thuyết phiếm hàm mật độ. Giả sử rằng với một hệ electron được mô tả bởi Hamilton N N N 1 1 Hˆ   (  i2 )   v(ri )   (1.14) i 1 2 i 1 i  j rij Z trong đó v(ri )   là thế “ngoài” (external) tác dụng lên electron i từ hạt nhân  ri có điện tích Z. ri là tọa độ không gian tổ hợp với tọa độ spin si thành tọa độ  i của electron.

Cả năng lượng trạng thái 0 và hàm sóng trạng thái 0 đều được xác định bằng cách cực tiểu hoá phiếm hàm năng lượng E[] Nguyễn Thu Hằng 7 Cao học K20 Hóa học LUAN VAN CHAT LUONG download : add [e-mail] Luận văn thạc sĩ khoa học  Hˆ  E   (1.15)  trong đó  Hˆ     * Hˆ dx và E   E0 (1.