Tổng quan nghiên cứu

Bệnh đau mắt đỏ (conjunctivitis) là một trong những bệnh lý nhiễm trùng nhãn khoa phổ biến nhất trên toàn cầu, với khả năng bùng phát thành các đợt dịch lớn tại các khu vực đô thị có mật độ dân cư đông đúc. Về mặt kinh tế y tế, gánh nặng từ căn bệnh này là vô cùng lớn; điển hình như tại Mỹ, chi phí điều trị đau mắt đỏ do nhiễm khuẩn đã tiêu tốn khoảng 589 triệu USD trong một năm, hay đợt dịch kéo dài 13 tuần tại đảo Réunion của Pháp đã gây thiệt hại hơn 3,3 triệu euro cho công tác kiểm soát dịch bệnh. Nguyên nhân hàng đầu gây ra các vụ dịch đau mắt đỏ nghiêm trọng là vi-rút Human Adenovirus (HAdV), một chi thuộc họ Adenoviridae với hạt vi-rút đường kính khoảng 90 nm, chứa hệ gen ADN sợi đôi kích thước xấp xỉ 36 kb. Trong số 103 chủng HAdV đã được công nhận trên thế giới, chủng HAdV-3 thuộc loài B là tác nhân lưu hành chủ yếu tại Việt Nam.

Trước thực trạng các nghiên cứu trong nước về HAdV hầu hết chỉ dừng lại ở mức độ định danh phân loại dịch tễ học cơ bản, việc nghiên cứu sâu về hệ gen và cấu trúc không gian của các protein vỏ capsid là yêu cầu cấp thiết để đánh giá độc lực và nguy cơ xuất hiện các biến chủng nguy hiểm. Đề tài luận văn thạc sĩ chuyên ngành Di truyền học thực hiện tại Đại học Quốc gia Hà Nội năm 2019 tập trung giải mã toàn bộ trình tự nucleotide của 3 gen cấu trúc chính gồm Fiber, Penton và Hexon từ chủng HAdV-3 phân lập tại Bệnh viện Mắt Trung ương Hà Nội. Thông qua việc phân tích tin sinh học và mô phỏng cấu trúc protein 3D, nghiên cứu hướng đến mục tiêu xác định chính xác các đột biến thay thế axit amin, đánh giá ảnh hưởng của biến đổi không gian đến khả năng gắn thụ thể tế bào chủ và cơ chế trốn thoát miễn dịch. Đây là công trình đầu tiên tại Việt Nam cung cấp bức tranh phân tử toàn diện về 3 protein vỏ capsid của HAdV-3 với độ chính xác giải trình tự đạt 99,99%, xác lập tỷ lệ biến đổi nucleotide từ 1,47% đến 1,56%, tạo cơ sở dữ liệu mốc quan trọng phục vụ công tác giám sát dịch tễ học phân tử và phát triển các biện pháp can thiệp y sinh trong tương lai.

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 được xây dựng trên nền tảng lý thuyết cấu trúc capsid vi-rút đối xứng 20 mặt (icosahedral capsid theory) và mô hình tương tác phân tử giữa phối tử vi-rút và thụ thể tế bào chủ (ligand-receptor interaction model). Vỏ capsid của HAdV được cấu tạo từ 240 khối trimer hexon, 12 khối pentamer penton và 12 phức hợp trimer fiber nhô ra ngoài. Quá trình xâm nhiễm tế bào của các chủng HAdV loài B (bao gồm HAdV-3) sử dụng cơ chế đặc thù: đầu tương tác fiber (fiber knob) nhận diện và liên kết trực tiếp với domain SCR1 và SCR2 của thụ thể CD46 trên màng tế bào người, kết hợp với phức hợp Dodecahedron-fiber tạo bởi penton để xâm nhập nội bào.

Khung phân tích protein dựa trên các nguyên lý hóa lý cơ bản của 20 axit amin tiêu chuẩn, phân thành 4 nhóm chức năng theo cực tính và điện tích. Hai thông số định lượng cốt lõi được áp dụng để đánh giá biến đổi phân tử là điểm đẳng điện pI (phản ánh trạng thái ion hóa và điện tích ròng ở các mức pH khác nhau) và chỉ số kỵ nước Hydropathy Index (chỉ số H theo thang đo Kyte-Doolittle, trong đó giá trị dương thể hiện tính kỵ nước mạnh và giá trị âm thể hiện tính ưa nước). Ngoài ra, nghiên cứu áp dụng mô hình mô phỏng tương đồng (Homology Modeling), dựa trên định đề cấu trúc không gian bậc ba của protein có tính bảo thủ tiến hóa cao hơn trình tự chuỗi polypeptide bậc một, cho phép tái hiện chính xác cấu hình 3D của protein từ dữ liệu ngân hàng Protein Data Bank (PDB).

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

Nghiên cứu sử dụng 1 chủng HAdV-3 đại diện thu thập từ mẫu dịch nước mắt của bệnh nhân đau mắt đỏ cấp tính tại Bệnh viện Mắt Trung ương Hà Nội, được bảo quản nghiêm ngặt ở nhiệt độ âm 20 độ C. Mẫu nghiên cứu được chọn lọc theo phương pháp chọn mẫu có chủ đích và được định danh chính xác kiểu gen HAdV-3 dựa trên trình tự vùng siêu biến HVR-7 của gen hexon.

Để giải mã hoàn chỉnh các gen có kích thước lớn, nhóm nghiên cứu đã thiết kế hệ thống mồi bắt cặp đặc hiệu gồm 2 cặp mồi cho gen fiber, 4 cặp mồi cho gen penton và 3 cặp mồi cho gen hexon dựa trên các vùng trình tự bảo thủ qua công cụ ClustalW. Kỹ thuật phản ứng chuỗi polymerase (PCR) được tiến hành với enzyme đọc sửa thế hệ mới Phusion Hotstart II DNA Polymerase nhằm triệt tiêu tối đa sai số khuếch đại, trải qua chu trình nhiệt 40 chu kỳ với tốc độ kéo dài chuẩn 30 giây trên 1 kb. Sản phẩm khuếch đại được kiểm tra độ tinh sạch và kích thước bằng điện di trên gel agarose nồng độ 1%, nhuộm Ethidium Bromide 0,5 µg/mL dưới ánh sáng tử ngoại 260-280 nm, sau đó tinh sạch bằng bộ kit chuyên dụng MEGAquick-spin Total fragment DNA purification.

Dữ liệu trình tự được xác định bằng phương pháp giải trình tự Sanger sử dụng hóa chất BigDye Terminator v3.1 thực hiện tại 1stBASE Malaysia. Phương pháp Sanger được ưu tiên lựa chọn thay vì giải trình tự thế hệ mới NGS do tính tối ưu về chi phí cho 3 gen đơn lẻ và đạt độ chính xác tiêu chuẩn vàng lên tới 99,99%. Dữ liệu thô được xử lý tín hiệu bằng phần mềm Bioedit, dịch mã chuỗi bằng SnapGene và đối chiếu với trình tự tham chiếu quốc tế NCBI mang mã số NC_011203.1. Cấu trúc không gian 3D được mô phỏng qua máy chủ SWISS-MODEL với các cấu trúc khuôn PDB tương đồng cao (mã 1H7Z cho fiber, 4aqq.A cho penton và 5TX1 cho hexon), sau đó trực quan hóa và phân tích tương tác trên phần mềm đồ họa phân tử pyMOL. Toàn bộ quy trình thực nghiệm và phân tích dữ liệu được hoàn tất trong mốc thời gian từ năm 2018 đến cuối năm 2019.

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

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

Quá trình giải trình tự và phân tích di truyền đã làm sáng tỏ cấu trúc phân tử của cả 3 gen protein vỏ thuộc chủng HAdV-3 lưu hành tại Việt Nam với nhiều đột biến đặc trưng:

Thứ nhất, đối với gen fiber, nghiên cứu đã giải mã thành công vùng mã hóa hoàn chỉnh dài 960 bp, tương ứng chuỗi polypeptide 320 axit amin. So sánh với trình tự tham chiếu NCBI NC_011203.1 phát hiện 15 đột biến nucleotide (tỷ lệ biến đổi ADN đạt 1,56%), trong đó có 10 đột biến nhầm nghĩa dẫn đến thay thế axit amin (tỷ lệ biến đổi axit amin là 3,13% trên toàn chuỗi protein). Tỷ lệ đột biến thay thế axit amin chiếm tới 66,67% tổng số biến đổi nucleotide của gen fiber.

Thứ hai, đối với gen penton, tổng chiều dài vùng mã hóa thu được là 1635 bp, mã hóa cho chuỗi protein gồm 545 axit amin. Kết quả ghi nhận 24 biến đổi nucleotide trên toàn chuỗi (tỷ lệ đột biến ADN là 1,47%), tạo ra 5 đột biến thay thế axit amin gồm I32L, T158I, I178T, D326N và A462T (tỷ lệ biến đổi axit amin đạt 0,92%, chiếm 20,83% tổng số đột biến nucleotide).

Thứ ba, đối với gen hexon, đoạn mã hóa hoàn chỉnh đạt chiều dài 2835 bp, tương ứng 945 axit amin. Phát hiện tổng cộng 42 vị trí biến đổi nucleotide (tỷ lệ biến đổi ADN là 1,48%), trong đó có 11 đột biến làm biến đổi axit amin gồm G42A, G141R, E299G, N302D, N411D, T418R, T429A, A439D, P440T, T445A và M884L (tỷ lệ biến đổi axit amin là 1,16%, chiếm 26,19% tổng số đột biến nucleotide). Đáng chú ý, có tới 9 trên 11 axit amin biến đổi (tương đương 81,82%) định vị tập trung ở mặt ngoài tiếp xúc của vỏ capsid.

Thảo luận kết quả

Các kết quả phân tích cấu trúc 3D trên phần mềm pyMOL đã làm sáng tỏ mối liên hệ trực tiếp giữa các đột biến phân tử và độc lực sinh học của vi-rút. Trên cấu trúc protein fiber, 4 đột biến thay thế axit amin (Q150E, S207L, H246D, M272T) định vị ngay tại vùng đầu tương tác tiếp giáp trực tiếp với domain SCR1 và SCR2 của thụ thể CD46 trên tế bào người. Đột biến S207L làm tăng vọt chỉ số kỵ nước từ âm 0,8 lên dương 3,8, trong khi H246D làm điểm đẳng điện thay đổi cực lớn từ tính kiềm nhẹ pI 7,64 xuống tính axit mạnh pI 2,87. Sự tái cấu trúc điện tích và độ phân cực này làm gia tăng đáng kể lực liên kết phối tử - thụ thể, giúp vi-rút bám dính nhanh và xâm nhập mạnh mẽ hơn vào biểu mô kết mạc. Toàn bộ các biến thiên hóa lý này có thể được biểu diễn trực quan qua bảng ma trận đối chiếu điểm đẳng điện pI kết hợp biểu đồ cột phân bố chỉ số kỵ nước H giữa chủng Việt Nam và chủng tham chiếu thế giới.

Trên protein penton, đột biến T158I nằm ở mặt trên làm tăng mạnh tính kỵ nước từ âm 0,7 lên dương 4,5, hỗ trợ liên kết bền vững với phần đuôi fiber. Đột biến I178T nằm ở mặt bên tiếp xúc với hexon, làm thay đổi tương tác trong cấu trúc nhóm sáu GOS liên kết các nhóm chín GON trên bề mặt capsid 20 mặt. Đặc biệt, biến đổi I32L tại đầu N tương tác với protein lõi IIIa đóng vai trò ổn định quá trình cuộn gập hạt nano dodecahedron, thúc đẩy tổng hợp cấu trúc DF giải phóng ồ ạt trong chu kỳ tan của tế bào.

Trên protein hexon, việc 9 trong số 11 biến đổi axit amin phân bố tại các quai siêu biến HVR bề mặt ngoài (điển hình như G141R, E299G, T418R, A439D với sự đảo lộn mạnh mẽ về điện tích và độ phân cực) chứng minh giả thuyết vi-rút HAdV-3 tại Việt Nam đã trải qua áp lực chọn lọc tự nhiên mạnh. Cấu hình kháng nguyên bề mặt bị biến đổi giúp vi-rút vô hiệu hóa sự nhận diện của kháng thể dịch thể IgA và né tránh phản ứng đáp ứng miễn dịch trung gian qua tế bào lympho T CD4+, lý giải vì sao các đợt dịch đau mắt đỏ do chủng này gây ra thường có tốc độ lây lan rất nhanh và tái nhiễm dễ dàng trong cộng đồng.

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

Dựa trên các phát hiện phân tử của đề tài, 4 nhóm giải pháp chiến lược được đề xuất nhằm nâng cao hiệu quả phòng ngừa và điều trị dịch bệnh:

Thứ nhất, thiết lập mạng lưới giám sát dịch tễ học phân tử Adenovirus quốc gia: Cục Y tế Dự phòng phối hợp cùng Viện Vệ sinh Dịch tễ Trung ương tiến hành lấy mẫu và giải trình tự gen định kỳ đối với 100% các ổ dịch đau mắt đỏ bùng phát tại các bệnh viện mắt tuyến tỉnh và trung ương. Mục tiêu xác định biến chủng mới trong vòng 48 giờ kể từ khi thu thập mẫu bệnh phẩm, triển khai liên tục từ năm 2024 đến năm 2026.

Thứ hai, phát triển bộ sinh phẩm Real-time PCR chẩn đoán đa mồi đặc hiệu cao: Giao Viện Hàn lâm Khoa học và Công nghệ Việt Nam phối hợp với các doanh nghiệp công nghệ sinh học thiết kế mẫu dò huỳnh quang nhắm trúng các vùng siêu biến HVR của gen hexon và các đột biến đặc trưng của gen fiber chủng HAdV-3 Việt Nam. Bộ sinh phẩm cần đạt độ nhạy trên 98%, thời gian trả kết quả dưới 3 giờ, hoàn thành thử nghiệm lâm sàng trong giai đoạn 2024 - 2025.

Thứ ba, khai thác cấu trúc phân tử để phát triển vector chuyển gen và thuốc ức chế bám dính: Các phòng thí nghiệm sinh học phân tử trọng điểm cần tận dụng dữ liệu cấu trúc phức hợp Dodecahedron-fiber (DF) từ đột biến I32L và T158I để nghiên cứu các dòng vector nano vận chuyển gen trị liệu vào tế bào người với hiệu suất tăng tối thiểu 25%, triển khai trong kế hoạch trung hạn 2025 - 2028.

Thứ tư, chuẩn hóa quy trình điều trị và kiểm soát nhiễm khuẩn nhãn khoa: Bộ Y tế kết hợp Hội Nhãn khoa Việt Nam ban hành hướng dẫn cập nhật về điều trị đau mắt đỏ vi-rút, kiểm soát chặt chẽ việc lạm dụng kháng sinh nhỏ mắt vốn không có tác dụng với HAdV, hướng đến mục tiêu kéo giảm tỷ lệ lây nhiễm chéo tại các cơ sở khám chữa bệnh xuống dưới mức 2%.

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

Nội dung và dữ liệu thực nghiệm của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng trọng tâm:

Thứ nhất, các nhà nghiên cứu và giảng viên chuyên ngành Sinh học phân tử, Di truyền học và Vi sinh học: Tài liệu cung cấp quy trình chuẩn mực từ thiết kế mồi nhiều phân đoạn, tối ưu hóa phản ứng PCR với enzyme đọc sửa, đến kỹ thuật xử lý tín hiệu giải trình tự Sanger và mô phỏng không gian 3D bằng SWISS-MODEL và pyMOL phục vụ công tác giảng dạy và đề tài nghiên cứu chuyên sâu.

Thứ hai, các bác sĩ chuyên khoa Nhãn khoa và cán bộ dịch tễ học lâm sàng: Luận văn cung cấp bằng chứng rõ ràng về cơ chế xâm nhiễm qua thụ thể CD46 và hiện tượng biến đổi kháng nguyên bề mặt hexon, giúp nhân viên y tế hiểu rõ căn nguyên sinh học của các đợt bùng phát dịch đau mắt đỏ, từ đó tối ưu hóa phác đồ chăm sóc triệu chứng và cách ly bệnh nhân.

Thứ ba, kỹ sư và chuyên viên phát triển sản phẩm tại các công ty công nghệ sinh học - dược phẩm: Nguồn số liệu về các vị trí đột biến điểm SNP và cấu trúc vùng bám thụ thể là dữ liệu đầu vào chuẩn xác để thiết kế kít chẩn đoán phân tử in vitro, bào chế thuốc nhỏ mắt peptide ức chế thụ thể và phát triển vector Adenovirus tái tổ hợp an toàn.

Thứ tư, học viên cao học, nghiên cứu sinh và sinh viên khối ngành Khoa học Sự sống và Y Dược: Luận văn là hình mẫu điển hình về phương pháp luận nghiên cứu kết hợp giữa kỹ thuật sinh học phân tử ướt (wet-lab) và công cụ tin sinh học khô (dry-lab), hỗ trợ học viên nâng cao kỹ năng xây dựng luận văn khoa học chuẩn mực.

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

Tại sao nghiên cứu lại chọn phương pháp giải trình tự Sanger thay vì công nghệ Next Generation Sequencing (NGS)? Phương pháp Sanger là lựa chọn tối ưu về mặt kinh tế và kỹ thuật khi giải mã 3 đoạn gen xác định có tổng độ dài hơn 5400 bp. Sanger được công nhận là tiêu chuẩn vàng với độ chính xác đạt 99,99%, cho tín hiệu huỳnh quang rõ nét, triệt tiêu nhiễu nền và hoàn toàn không gặp lỗi lắp ráp dữ liệu như các hệ thống NGS khi phân tích số lượng gen hạn chế.

Các đột biến trên protein Fiber làm thay đổi khả năng gây bệnh của vi-rút như thế nào? Bốn đột biến gồm Q150E, S207L, H246D và M272T nằm ngay tại vùng đầu tương tác tiếp xúc trực tiếp thụ thể CD46 của tế bào người. Đặc biệt, biến đổi S207L làm tăng chỉ số kỵ nước từ âm 0,8 lên dương 3,8 và H246D hạ điểm đẳng điện pI từ 7,64 xuống 2,87 (chuyển sang tính axit mạnh), giúp tăng ái lực bám dính vào màng tế bào kết mạc.

Đột biến I32L ở protein Penton có vai trò gì trong vòng đời của HAdV-3? Đột biến I32L nằm ở đầu N của penton là vị trí tham gia liên kết trực tiếp với hexon và protein IIIa. Đột biến này tạo thuận lợi cho sự cuộn gập cấu trúc và thúc đẩy hình thành các phức hợp đa diện Dodecahedron-fiber (DF) với số lượng vượt trội, giúp vi-rút dễ dàng phá hủy cấu trúc màng để xâm nhập và giải phóng hạt vi-rút mới.

Cơ chế nào giúp HAdV-3 tại Việt Nam trốn thoát hệ miễn dịch của người bệnh? Có 9 trên 11 đột biến axit amin của hexon (chiếm tỷ lệ 81,82%) phân bố tại các quai siêu biến HVR nằm ở bề mặt ngoài cùng của vỏ capsid. Các biến đổi phân cực và điện tích này làm biến dạng các vị trí gắn kháng thể, khiến kháng thể dịch thể IgA và tế bào lympho T CD4+ mất khả năng nhận diện và trung hòa vi-rút hiệu quả.

Nghiên cứu đóng góp gì cho sự phát triển của công nghệ sinh học y tế tại Việt Nam? Nghiên cứu đã xác lập bộ cơ sở dữ liệu phân tử đầu tiên tại Việt Nam với 81 điểm biến đổi nucleotide trên cả 3 gen vỏ capsid HAdV-3. Dữ liệu này là tiền đề then chốt để thiết kế mồi chẩn đoán nhanh, đồng thời định hướng nghiên cứu phát triển các hạt nano vi-rút làm vector chuyển gen trị liệu trong y học tái tạo.

Kết luận

Công trình nghiên cứu đã mang lại những đóng góp khoa học quan trọng và toàn diện trong lĩnh vực di truyền học vi sinh vật tại Việt Nam:

  • Giải mã thành công và chính xác 100% trình tự nucleotide của 3 gen vỏ capsid chủng HAdV-3 gây bệnh đau mắt đỏ tại Việt Nam với chiều dài 960 bp (gen fiber), 1635 bp (gen penton) và 2835 bp (gen hexon).
  • Xác định tổng cộng 81 biến đổi nucleotide và 26 đột biến thay thế axit amin trên toàn bộ 3 protein vỏ, thiết lập tỷ lệ đột biến axit amin lần lượt là 3,13% ở fiber, 0,92% ở penton và 1,16% ở hexon so với chủng tham chiếu quốc tế.
  • Làm sáng tỏ cơ chế gia tăng độc lực thông qua 4 đột biến tại đầu tương tác fiber giúp tăng cường liên kết với thụ thể tế bào CD46 và đột biến I32L ở đầu N penton giúp ổn định cấu trúc hạt nano Dodecahedron.
  • Chứng minh khả năng trốn thoát miễn dịch tự nhiên của chủng HAdV-3 Việt Nam nhờ 9 đột biến định vị tại các quai siêu biến HVR trên bề mặt ngoài của protein hexon.
  • Xây dựng thành công mô hình cấu trúc không gian 3D của cả 3 protein vỏ, cung cấp nguồn tư liệu chuẩn mực phục vụ định hướng thiết kế vector chuyển gen và sinh phẩm y tế thế hệ mới.

Trong giai đoạn 2024 - 2026, các hướng nghiên cứu tiếp theo cần tập trung mở rộng giải mã toàn bộ hệ gen vi-rút (36 kb) và tiến hành đo đạc thực nghiệm ái lực liên kết thụ thể CD46 trên mô hình nuôi cấy tế bào sống. Các nhà khoa học, bác sĩ nhãn khoa và cơ quan quản lý y tế hãy chủ động khai thác nguồn dữ liệu phân tử quý giá này để xây dựng các giải pháp giám sát dịch bệnh chủ động và phát triển các chế phẩm bảo vệ sức khỏe nhãn khoa cho cộng đồng.