Bảng 5: Mồi khuếch đại gen fiber của HAdV-3 ở Việt Nam Tên mồi Trình tự (chiều 5’-3’) Kích thước sản phẩm khuếch đại (bp) FF1 GGGTGGAATTCTCCCAATG 2818 FR1 TGAAACCCGAGCTTCTATGA FF2 CTTCGAGACCTCCTACCCAT 1531 FR2 TTCGGATTATGATTCCCATCG
3.1.2. Kết quả khuếch đại PCR gen fiber của HAdV-3 tại Việt Nam
Hình 8 cho thấy băng sản phẩm PCR của cả 2 cặp mồi đều sáng bằng hoặc
sáng hơn marker tương ứng, điều này cho thấy lượng sản phẩm khuếch đại đã đạt chuẩn. Vị trí của các băng cũng cho thấy kích thước các sản phẩm khuếch đại đều phù hợp với tính tốn lý thuyết khi thiết kế mồi. Bên cạnh đó, khơng quan sát thấy có vệt smear hoặc sản phẩm phụ nào, chứng tỏ sản phẩm khuếch đại rất đặc hiệu và đáng tin cậy. Như vậy, các sản phẩm khuếch đại bằng phương pháp PCR mà chúng tôi thực hiện, đều đủ điều kiện để được tinh sạch trực tiếp, không cần phải thực hiện thôi gel để loại bỏ các sản phẩm không mong muốn. Sản phẩm ADN sau tinh sạch đã được gửi cho công ty 1stBASE để giải trình tự.
Hình 8: Hình ảnh điện di sản phẩm khuếch đại gen fiber
3.1.3. Kết quả giải trình tự gen fiber của HAdV-3 tại Việt Nam
Hình 9 cho thấy các đỉnh tín hiệu giải trình tự gen fiber đều rất rõ ràng và đạt
độ phân giải cao. Tín hiệu nhiễu nền gần như khơng có. Tổng độ dài đã giải được từ 2 cặp mồi FF1-FR1 và FF2-FR2 là 2707 bp. Tuy nhiên, sau khi lược bỏ những trình tự khơng thuộc vùng mã hóa thì chiều dài gen fiber thu được là 960 bp. Như vậy, độ dài chuỗi axit amin của protein fiber của HAdV-3 tại Việt Nam là 320 axit amin. Kích thước này cũng tương tự so với trình tự đã cơng bố trước đó của HAdV-3 trên cơ sở dữ liệu NCBI. Chứng tỏ, chúng tôi đã giải thành cơng trình tự của gen fiber của HAdV-3 gây bệnh đau mắt đỏ tại Việt Nam với độ tin cậy cao.
3.1.4. Kết quả so sánh trình tự mã hóa gen fiber của HAdV-3 tại Việt Nam
Sau khi thu được trình tự ADN hồn chỉnh của gen fiber. Chúng tơi đã so sánh với trình tự genome tham khảo của HAdV-3 trên cơ sở dữ liệu NCBI, mã trình tự là NC_011203.1. Thơng qua việc so sánh này, chúng tôi đã phát hiện được những biến đổi ở cả cấp độ ADN và axit amin của gen và protein fiber. Kết quả so sánh được biểu diễn trên phần mềm Bioedit (Hình 10), và được tổng hợp lại trong Bảng 6.
Bảng 6: Các biến đổi trên gen fiber của HAdV-3 tại Việt Nam so với trình tự tham chiếu
STT Vị trí Nucleotide Axit amin biến đổi NC_011203.1 HAdV-3 Việt Nam
1 31432 G A S22N 2 31435 C T S23L 3 31706 A G - 4 31751 G C - 5 31815 C G Q150E 6 31904 C T - 7 31947 C T - 8 31987 C T S207L 9 32010 G C E215Q 10 32031 G A A222T 11 32034 G C D223H 12 32103 C G H246D 13 32153 T C M272T 14 32182 T C - 15 32314 G C R316T
Hình 10: Hình ảnh kết quả so sánh trình tự ADN và axit amin của gen fiber
Như vậy, phần mã hóa gen fiber của HAdV-3 tại Việt Nam có 15 biến đổi ADN trên tổng chiều dài là 960 bp. Trong đó, 10 biến đổi ADN dẫn đến thay thế axit amin. Từ đó, chúng tơi đã tính ra tỉ lệ biến đổi trình tự ở mức độ ADN của gen fiber là 1,56%. Tương tự, tỉ lệ số lượng axit amin bị thay thế trên tổng số biến đổi nucleotide của gen fiber 66.67%. Cuối cùng, tỉ lệ biến đổi axit amin trên tổng số axit amin của protein fiber là 3,13%.
Để có cái nhìn rõ ràng hơn về vị trí của các biến đổi axit amin trên cấu trúc không gian của protein fiber, và đánh giá được chi tiết hơn về những tác động của các biến đổi này có thể mang lại, chúng tơi thực hiện bước tiếp theo là xây dựng cấu trúc 3D mô phỏng của protein này.
3.1.5. Kết quả phân tích cấu trúc ba chiều protein fiber của HAdV-3 tại Việt Nam Sử dụng chương trình SWISS-MODEL kết hợp với dữ liệu của ngân hàng dữ liệu protein (PDB), chúng tôi đã mô phỏng thành công cấu trúc 3 chiều của protein fiber của HAdV-3 gây bệnh đau mắt đỏ ở Việt Nam. Các biến đổi axit amin được thể hiện dưới dạng tinh thể (sphere) trong phần mềm pyMOL. Có một số biến đổi khơng được thể hiện lên hình là do phần mềm không xây dựng được phần cấu trúc mang biến đổi đó. Cụ thể, phần cấu trúc khơng xây dựng được là cấu trúc trục của protein fiber.
Cấu trúc ba chiều của đầu tương tác fiber HAdV-3 ở Việt Nam với thụ thể của tế bào vật chủ, dựa trên mã truy cập trong cơ sở dữ liệu PDB là 1H7Z [15], vì đây là cấu trúc có độ tương đồng trình tự lớn nhất với trình tự fiber HAdV-3 ở Việt Nam (95,36%). Tuy nhiên, cấu trúc này không thể hiện được liên kết với domain SCR1-2 của thụ thể CD46, nên chúng tôi đã tiến hành xếp chồng cấu trúc của 1H7Z (HAdV-
3) với 2O39 (HAdV-11) [40] và 3L89 (HAdV-21) [13] để đánh giá mức độ tương đồng cấu trúc của chúng (Hình 11). Kết quả cho thấy rằng cấu trúc đầu tương tác fiber của 1H7Z, 2O39 và 3L89 gần như trùng khớp lên nhau cho dù trình tự axit amin của chúng rất khác biệt. Từ đó, chúng tơi đã xây dựng cấu trúc của fiber HAdV-3 ở Việt Nam dựa trên 1H7Z và thể hiện sự liên kết với các domain SCR1-2 của thụ thể CD46 của 2O39 (Hình 12B).
Hình 11: Hình ảnh xếp chồng của đầu tương tác protein fiber HAdV
Hình 12: Đầu tương tác fiber của HAdV-3 tại Việt Nam và domain SCR1-2 của thụ thể CD46
A. So sánh trình tự axit amin protein fiber, B. Hướng nhìn từ bên ngồi vi-rút
Trong số 10 axit amin bị biến đổi, có 6 biến đổi khơng được thể hiện trên hình gồm: bốn biến đổi (E215Q, A222T, D223H, R316T) cũng nằm ở đầu tương tác của fiber nhưng ở vị trí khơng tương tác trực tiếp với thụ thể CD46, và hai biến đổi (S22N, S23L) nằm ở trục của protein fiber (phần mềm không xây dựng được cấu trúc). Bốn biến đổi còn lại (Q150E, S207L, H246D, M272T) nằm rất gần với domain SCR1-2 nên có khả năng lớn ảnh hưởng tới sự tương tác giữa đầu tương tác fiber của HAdV- 3 và thụ thể CD46 ở người. Vì vậy chúng tơi đã tổng hợp lại giá trị pI và chỉ số H của các axit amin trong 4 biến đổi này để tiện cho việc so sánh và đánh giá hậu quả (Bảng
7).
A