Luận án tiến sĩ vật lý: Động lực học cuốn protein trong môi trường tế bào - Bùi Phương Thúy (Viện Hàn Lâm KHCN Việt Nam)

Cấu trúc không gian của chuỗi polypeptide được duy trì bởi hệ thống liên kết vật lý đa dạng. Tương tác kỵ nước đóng vai trò then chốt trong việc đưa các chuỗi bên không phân cực vào lõi protein. Lực Van der Waals giúp cố định vị trí các nguyên tử ở khoảng cách tối ưu. Cầu hydro liên kết các nhóm amino và carboxyl dọc theo khung xương peptide. Cầu nối disulfide đồng hóa trị tăng cường độ bền vững cấu trúc cho protein ngoại bào. Tương tác tĩnh điện giữa các nhóm tích điện tạo nên các cầu muối đặc hiệu. Sự cạnh tranh giữa năng lượng liên kết và entropy cấu hình chi phối toàn bộ quá trình. Mỗi tương tác nội phân tử đóng góp một phần nhỏ vào độ ổn định tổng thể của trạng thái tự nhiên. Cấu trúc ba chiều hoàn chỉnh là kết quả của sự tối ưu hóa năng lượng toàn cục.

Sinh tổng hợp protein bắt đầu từ quá trình dịch mã thông tin di truyền tại ribosome. Phân tử mRNA mang thông tin mã hóa từ nhân tế bào ra tế bào chất. Phức hợp ribosome đọc từng bộ ba mã sao để gắn các amino acid tương ứng. Chuỗi polypeptide kéo dài dần bên trong đường hầm ribosome hẹp. Tốc độ dịch mã dao động tùy thuộc vào độ hiếm của codon và cấu trúc mRNA. Sự chuyển động của chuỗi mới sinh qua đường hầm bị hạn chế không gian đáng kể. Quá trình cuộn gập có thể bắt đầu ngay khi chuỗi peptide vẫn còn gắn với ribosome. Đây là hiện tượng cuộn gập đồng dịch mã (co-translational folding). Tốc độ tổng hợp peptide ảnh hưởng trực tiếp đến con đường hình thành cấu trúc không gian cục bộ.

Nghiên cứu ống nghiệm (in vitro) giúp làm sáng tỏ các quy luật vật lý thuần túy của chuỗi polypeptide. Môi trường ống nghiệm thường là dung dịch đệm loãng với nồng độ protein rất thấp. Thí nghiệm kinh điển của Anfinsen chứng minh thông tin cuộn gập nằm hoàn toàn trong trình tự amino acid. Protein bị biến tính hoàn toàn có thể tự tái tạo cấu trúc tự nhiên khi loại bỏ chất làm biến tính. Quá trình tự tổ chức diễn ra tự phát theo nguyên lý cực tiểu hóa năng lượng tự do. Tuy nhiên, môi trường ống nghiệm không phản ánh đầy đủ các tác động phức tạp của tế bào sống. Thiếu hụt các rào cản không gian và tương tác đại phân tử làm thay đổi đáng kể động lực học phản ứng. Việc so sánh giữa in vitro và in vivo là trọng tâm then chốt của sinh vật lý hiện đại.

Nghịch lý Levinthal chỉ ra rằng chuỗi polypeptide không thể tìm kiếm cấu trúc tự nhiên bằng cách duyệt ngẫu nhiên mọi cấu hình. Nếu tìm kiếm ngẫu nhiên, thời gian cuộn gập sẽ vượt quá tuổi của vũ trụ. Bề mặt năng lượng tự do (free energy landscape) định hướng dòng cấu hình trượt xuống đáy phễu một cách nhanh chóng. Sự giảm entropy cấu hình được bù đắp bởi sự giải phóng năng lượng liên kết nội tại. Dạng phễu tạo ra một gradient năng lượng tự do hướng tâm rõ rệt. Quá trình cuộn gập diễn ra đồng thời dọc theo nhiều lộ trình song song thay vì một con đường cố định duy nhất. Các cấu hình trung gian xuất hiện thoáng qua và liên tục điều chỉnh. Nhờ vậy, protein đạt trạng thái tự nhiên chỉ trong vài mili-giây đến vài giây.

Cơ chế tạo nhân cuộn gập mô tả sự hình thành sớm của một cụm cấu trúc cục bộ bền vững. Cụm hạt nhân này đóng vai trò mầm mống định hướng cho các phần còn lại của chuỗi peptide. Trạng thái chuyển tiếp nằm trên đỉnh rào cản năng lượng tự do giữa trạng thái duỗi và trạng thái cuộn. Tại trạng thái chuyển tiếp, một phần cấu trúc tự nhiên đã được thiết lập ổn định. Vượt qua bờ thế năng này là bước giới hạn tốc độ của toàn bộ quá trình cuộn gập hai trạng thái. Phân tích giá trị phi (Phi-value analysis) thực nghiệm cho phép xác định chính xác các gốc amino acid tham gia vào nhân cuộn. Khi nhân hình thành, năng lượng tự do giảm mạnh và kéo theo sự đóng gói nhanh chóng của toàn bộ phân tử.

Nguyên lý thất vọng tối thiểu khẳng định các tương tác trong trạng thái tự nhiên có tính tương hỗ mạnh mẽ. Xung đột năng lượng giữa các gốc amino acid được giảm thiểu đến mức tối đa qua chọn lọc tự nhiên. Tính nhất quán giúp chuỗi peptide tránh bị bẫy trong các trạng thái cấu hình giả ổn định. Các tiếp xúc không bản địa (non-native contacts) bị suy yếu nhanh chóng trong quá trình trượt dốc năng lượng. Ngược lại, các tiếp xúc bản địa liên tục được củng cố và tích lũy. Độ gồ ghề của bề mặt thế năng phụ thuộc vào trình tự chuỗi và điều kiện nhiệt độ môi trường. Sự cân bằng hoàn hảo giữa tính linh động và độ ổn định nhiệt động lực học giúp protein thực hiện chức năng sinh học tối ưu.

Quá trình thoát của chuỗi peptide qua đường hầm ribosome được mô tả hiệu quả bằng lý thuyết khuếch tán một chiều. Chuyển động Brown của chuỗi bị giới hạn dọc theo trục đối xứng của đường hầm. Lực đẩy entropy từ không gian tự do bên ngoài tế bào chất hỗ trợ kéo chuỗi ra khỏi ống hẹp. Thành đường hầm thuần đẩy tạo ra ít ma sát hơn, giúp rút ngắn thời gian thoát của chuỗi. Ngược lại, tương tác hút tĩnh điện hoặc kỵ nước giữa chuỗi peptide và thành ống tạo ra các hố bẫy thế năng. Sự hiện diện của các vị trí hút làm tăng thời gian lưu trú và làm chậm đáng kể quá trình tổng hợp hoàn chỉnh. Động lực học thoát phụ thuộc phi tuyến tính vào chiều dài chuỗi và bán kính hiệu dụng của đường hầm.

Vùng cửa ra của đường hầm ribosome tạo tiền đề cho sự cuộn gập cấu trúc bậc ba của protein. Khi các miền cấu trúc (domains) nhô dần ra ngoài, không gian thể tích tự do mở rộng đột ngột. Sự hình thành nhân cuộn gập tại miệng đường hầm tạo ra lực kéo cơ học hỗ trợ giải phóng phần đuôi còn lại. Quá trình này diễn ra tuần tự từ đầu N đến đầu C của chuỗi polypeptide. Cuộn gập đồng dịch mã giúp phân tách các miền không gian, tránh tương tác sai lệch giữa các đoạn chưa tổng hợp. Cấu trúc bề mặt ribosome tại vùng tiếp giáp cũng ảnh hưởng đến định hướng không gian của protein. Các quan sát mô phỏng khẳng định cửa ra ribosome là trạm kiểm soát chất lượng cuộn gập đầu tiên trong tế bào.

Lý thuyết hạt điều chỉnh tỷ lệ (scaled particle theory) của Minton và Zhou là công cụ toán học quan trọng trong việc định lượng thể tích loại trừ. Lý thuyết này tính toán chi phí năng lượng tự do để tạo một khoang trống chứa protein trong dung dịch đông đúc. Dạng hình học của các hạt chen chúc và phân tử protein được đơn giản hóa thành hình cầu hoặc hình trụ. Kết quả giải tích cho thấy trạng thái duỗi thẳng bị mất mát entropy nhiều hơn đáng kể so với trạng thái cuộn chặt. Do đó, hằng số cân bằng hóa học dịch chuyển mạnh mẽ về phía trạng thái cấu hình tự nhiên. Lý thuyết của Zhou mở rộng phân tích cho các chuỗi polymer mềm dẻo với nhiều bậc tự do không gian. Các dự đoán lý thuyết này phù hợp chặt chẽ với số liệu đo phổ huỳnh quang thực nghiệm.

Mặc dù thể tích loại trừ tạo lực đẩy làm bền cấu trúc cuộn, tương tác hút không đặc hiệu lại gây hiệu ứng ngược lại. Các đại phân tử chen chúc trong tế bào thực tế luôn mang điện tích hoặc vùng kỵ nước bề mặt. Lực hút yếu giữa protein và các tác nhân chèn ép tạo điều kiện cho trạng thái duỗi mở rộng tiếp xúc. Tương tác hút này làm mất ổn định trạng thái tự nhiên và hạ thấp rào cản năng lượng tự do chuyển tiếp. Động lực học cuộn gập protein trở nên phức tạp do sự giằng co giữa lực đẩy thể tích và lực hút tĩnh điện. Nếu lực hút chiếm ưu thế, thời gian đạt cấu trúc bản địa sẽ bị kéo dài đáng kể. Việc tính toán đồng thời cả hai yếu tố này là yêu cầu bắt buộc để mô tả chính xác môi trường nội bào thực tế.

Sự chen chúc của các hạt đại phân tử làm biến đổi sâu sắc chiều cao và vị trí của bờ thế năng lượng tự do. Khi mật độ hạt chèn tăng, rào cản nhiệt động lực học giữa trạng thái mở và trạng thái cuộn có xu hướng giảm. Vị trí của trạng thái chuyển tiếp dịch chuyển về phía cấu hình nhỏ gọn hơn dọc theo tọa độ phản ứng. Tốc độ cuộn gập có thể tăng nhanh ở nồng độ chất chen chúc vừa phải do hiệu ứng nén không gian. Tuy nhiên, độ nhớt dung dịch tăng cao ở mật độ lớn lại cản trở quá trình khuếch tán nội phân tử. Sự cạnh tranh giữa động học tiếp xúc và độ nhớt tạo ra một điểm cực trị tối ưu cho tốc độ phản ứng. Phân tích bờ thế năng cung cấp cái nhìn định lượng về hiệu suất biến đổi cấu trúc protein trong tế bào sống.

Các rào cản năng lượng cao thường khiến mô phỏng MD truyền thống bị mắc kẹt tại các cực tiểu cục bộ. Phương pháp lấy mẫu ô (umbrella sampling) khắc phục hạn chế này bằng cách áp dụng các thế năng giữ nhân tạo dọc theo tọa độ phản ứng. Không gian pha được chia thành nhiều cửa sổ lấy mẫu độc lập nhằm bao phủ trọn vẹn lộ trình chuyển trạng thái. Phương pháp phân tích biểu đồ có trọng số (WHAM - Weighted Histogram Analysis Method) tổng hợp dữ liệu xác suất từ tất cả các cửa sổ. Kỹ thuật này loại bỏ ảnh hưởng của thế năng cưỡng bức và tái thiết lập bề mặt năng lượng tự do không sai lệch. WHAM cho phép xác định chính xác năng lượng tự do tuyệt đối của trạng thái chuyển tiếp và trạng thái tự nhiên.

Mô hình Go đơn giản hóa cấu trúc protein bằng cách gán mỗi amino acid thành một hạt tương tác duy nhất. Trường lực Go giả định năng lượng hấp dẫn chỉ tồn tại giữa các cặp hạt có tiếp xúc trong cấu trúc tự nhiên. Cách tiếp cận này loại trừ hiện tượng thất vọng năng lượng, giúp tập trung phân tích vai trò của cấu trúc liên kết hình học. Mô hình hạt thô cắt giảm đáng kể bậc tự do tính toán, cho phép khảo sát các hệ phức tạp chứa hàng trăm hạt chen chúc. Nghiên cứu trên các protein tiêu chuẩn như GB1 và Ubiquitin chứng minh mô hình Go tái tạo chính xác cơ chế cuộn gập hai trạng thái. Phương pháp này cung cấp cầu nối trực tiếp giữa lý thuyết vật lý thống kê và dữ liệu cấu trúc tinh thể học.

Mô phỏng môi trường tế bào đòi hỏi thiết lập điều kiện biên tuần hoàn cùng mật độ hạt chèn thể tích xác định. Hộp mô phỏng chứa chuỗi protein trung tâm được bao quanh bởi các hạt đại phân tử hình cầu với bán kính khác nhau. Tương tác giữa các hạt được mô tả qua thế năng Lennard-Jones hoặc thế năng thuần đẩy Weeks-Chandler-Andersen (WCA). Việc kiểm soát áp suất và nhiệt độ thông qua các bộ điều nhiệt (thermostat) đảm bảo hệ thống duy trì trạng thái cân bằng vĩ mô. Sự phân bố không gian của các hạt chen chúc phản ánh chính xác độ phức tạp của tế bào chất vi khuẩn hoặc sinh vật nhân thực. Dữ liệu quỹ đạo thu được cung cấp cái nhìn trực quan về sự biến dạng liên tục của bề mặt năng lượng tự do.

Phức hợp chaperone GroEL-GroES hoạt động theo cơ chế lồng nano bảo vệ (nano-cage mechanism). Khoang rỗng bên trong GroEL tạo ra một môi trường biệt lập vật lý hoàn toàn cho một phân tử protein đơn lẻ. Buồng giam cầm nano này ngăn ngừa tuyệt đối nguy cơ tương tác va chạm giữa các chuỗi peptide với nhau. Chu kỳ thủy phân ATP thúc đẩy sự thay đổi hình thể bề mặt bên trong khoang từ kỵ nước sang ưa nước. Sự đảo chiều tính chất hóa lý này ép chuỗi protein phải chôn giấu các gốc kỵ nước vào tâm để hoàn tất quá trình cuộn gập. Sau khi gấp nếp chuẩn xác, nắp GroES mở ra và giải phóng protein chức năng vào tế bào chất. Mô hình lồng chaperone là minh chứng điển hình cho việc kiểm soát năng lượng tự do bằng rào cản cơ học.

Hiểu biết sâu sắc về động lực học cuộn gập protein mở ra những triển vọng lớn cho y học phân tử hiện đại. Nhiều căn bệnh hiểm nghèo như Alzheimer, Parkinson và tiểu đường tuýp 2 bắt nguồn trực tiếp từ sự kết tụ sợi amyloid độc hại. Việc giải mã bề mặt năng lượng tự do giúp xác định chính xác các trạng thái trung gian gây độc tế bào. Các nhà khoa học phát triển các phân tử nhỏ (chaperone hóa học) nhằm ổn định cấu trúc tự nhiên hoặc ngăn chặn tiếp xúc kỵ nước bất thường. Ngoài ra, việc tối ưu hóa động học cuộn gập trong tế bào giúp nâng cao hiệu suất biểu hiện protein tái tổ hợp trong công nghiệp sinh dược phẩm. Nghiên cứu vật lý lý thuyết tạo nền tảng vững chắc cho việc thiết kế thuốc định hướng cấu trúc trong tương lai.