Luận văn thạc sĩ: Nghiên cứu mô phỏng tĩnh điện trong sinh học

Luận văn thạc sĩ HUS nghiên cứu mô phỏng tĩnh điện trong sinh học, phân tích các bài toán và ứng dụng thực tiễn trong lĩnh vực này.

• Tác giả: Nguyễn Thị Trang
• Người hướng dẫn: TS. Nguyễn Thế Toàn
• Trường học: Trường Đại học Khoa học Tự nhiên
• Chuyên ngành: Vật lí lí thuyết và vật lí toán
• Đề tài: Nghiên cứu mô phỏng của 1 số bài toán về tĩnh điện trong sinh học
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2015
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• MỞ ĐẦU
• 0.1. Lý do chọn đề tài
• 0.2. Phương pháp nghiên cứu
• 0.3. Bố cục luận văn
• 1. CHƯƠNG 1: TỔNG QUAN VỀ CẤU TRÚC VÀ MỘT SỐ TÍNH CHẤT VẬT LÝ CỦA PHÂN TỬ ADN
• 1.1. Cấu trúc hóa học của phân tử ADN
• 1.2. Cấu trúc không gian của phân tử ADN
• 1.2.1. Cấu trúc sơ cấp
• 1.2.2. Hiện tượng cuộn xoắn và các cấu trúc bậc cao của ADN
• 1.2.3. ADN được xem như một polymer tích điện, tính chất đàn hồi của ADN
• 1.2.3.1. Định nghĩa độ dài quán tính của dây polymer
• 1.2.3.2. Đóng góp của hiệu ứng tĩnh điện vào độ dài quán tính của ADN
• 2. CHƯƠNG 2: PHƯƠNG TRÌNH POISSON – BOLTZMANN VÀ LÝ THUYẾT TĨNH ĐIỆN DEBYE – HUCKEL
• 2.1. Phương trình Poisson – Boltzmann
• 2.1.1. Phương trình Poisson – Boltzmann; Phương trình trường trung bình tự hợp
• 2.1.2. Tuyến tính hóa phương trình Poisson – Boltzmann (PB); Phương trình Debye – Huckel (DH)
• 2.2. Phương trình Poisson - Boltzmann cho thế năng quanh một hình trụ tích điện
• 2.2.1. Nghiệm của phương trình Debye – Huckel cho thế năng quanh hình trụ tích điện
• 2.2.2. Hiện tượng cô đọng “Manning”
• 3. CHƯƠNG 3: PHƯƠNG PHÁP MÔ PHỎNG MONTE CARLO
• 3.1. Phương pháp mô phỏng Monte Carlo
• 3.2. Điều kiện cân bằng chi tiết và thuật toán Metropolis
• 3.2.1. Điều kiện cân bằng chi tiết
• 3.2.2. Thuật toán Metropolis
• 4. CHƯƠNG 4: NHỮNG KẾT QUẢ MÔ PHỎNG THU ĐƯỢC CHO DÂY POLYMER TÍCH ĐIỆN
• 4.1. Mô hình dây polymer tích điện
• 4.2. Những kết quả mô phỏng thu được
• 4.2.1. Sự phụ thuộc một cách tỷ lệ của độ dài đầu – cuối vào bán kính chắn
• 4.2.2. Những kết quả thu được tại giới hạn bán kính chắn lớn
• 4.2.3. Hàm tương quan góc của các liên kết
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về nghiên cứu mô phỏng tĩnh điện trong sinh học

Nghiên cứu mô phỏng tĩnh điện trong sinh học là một lĩnh vực quan trọng, giúp hiểu rõ hơn về các tương tác điện trong các hệ sinh học. Tĩnh điện có vai trò quyết định trong cấu trúc và hoạt động của các phân tử sinh học như ADN và protein. Việc mô phỏng tĩnh điện cho phép các nhà khoa học dự đoán hành vi của các phân tử trong môi trường sinh học, từ đó phát triển các ứng dụng trong y học và công nghệ sinh học.

1.1. Khái niệm về tĩnh điện trong sinh học

Tĩnh điện trong sinh học đề cập đến các lực tương tác giữa các phân tử tích điện. Các phân tử như ADN và protein thường mang điện tích, ảnh hưởng đến cách chúng tương tác với nhau và với môi trường xung quanh.

1.2. Tầm quan trọng của mô phỏng tĩnh điện

Mô phỏng tĩnh điện giúp hiểu rõ hơn về các quá trình sinh học như sao chép ADN, phiên mã và tương tác protein. Điều này có thể dẫn đến những phát hiện mới trong nghiên cứu bệnh học và phát triển thuốc.

II. Vấn đề và thách thức trong nghiên cứu tĩnh điện

Mặc dù tĩnh điện có vai trò quan trọng, nhưng việc nghiên cứu và mô phỏng các tương tác tĩnh điện trong sinh học gặp nhiều thách thức. Các yếu tố như độ chính xác của mô hình, điều kiện môi trường và sự phức tạp của các hệ sinh học đều ảnh hưởng đến kết quả nghiên cứu.

2.1. Độ chính xác của mô hình mô phỏng

Một trong những thách thức lớn nhất là đảm bảo độ chính xác của các mô hình mô phỏng. Các mô hình cần phản ánh đúng các tương tác tĩnh điện trong môi trường sinh học thực tế.

2.2. Ảnh hưởng của điều kiện môi trường

Điều kiện môi trường như pH, nhiệt độ và nồng độ ion có thể ảnh hưởng lớn đến các tương tác tĩnh điện. Việc kiểm soát và mô phỏng các điều kiện này là rất quan trọng.

III. Phương pháp mô phỏng tĩnh điện hiệu quả trong sinh học

Có nhiều phương pháp mô phỏng tĩnh điện được sử dụng trong nghiên cứu sinh học, bao gồm phương pháp Monte Carlo và phương pháp động lực học phân tử. Những phương pháp này giúp mô phỏng các tương tác tĩnh điện một cách chính xác và hiệu quả.

3.1. Phương pháp Monte Carlo

Phương pháp Monte Carlo là một kỹ thuật mô phỏng dựa trên việc sử dụng các mẫu ngẫu nhiên để dự đoán hành vi của các hệ thống phức tạp. Phương pháp này rất hữu ích trong việc nghiên cứu tĩnh điện.

3.2. Động lực học phân tử

Động lực học phân tử cho phép mô phỏng chuyển động của các phân tử trong thời gian thực. Phương pháp này giúp hiểu rõ hơn về các tương tác tĩnh điện trong các hệ sinh học.

IV. Ứng dụng thực tiễn của mô phỏng tĩnh điện trong sinh học

Mô phỏng tĩnh điện có nhiều ứng dụng thực tiễn trong nghiên cứu sinh học, từ việc phát triển thuốc đến hiểu rõ hơn về các bệnh lý. Các nghiên cứu này không chỉ giúp cải thiện hiểu biết về sinh học mà còn có thể dẫn đến các phương pháp điều trị mới.

4.1. Phát triển thuốc mới

Mô phỏng tĩnh điện giúp xác định cách các phân tử thuốc tương tác với các mục tiêu sinh học, từ đó phát triển các loại thuốc hiệu quả hơn.

4.2. Nghiên cứu bệnh lý

Các nghiên cứu về tĩnh điện có thể giúp hiểu rõ hơn về cơ chế gây bệnh, từ đó tìm ra các phương pháp điều trị hiệu quả hơn cho các bệnh lý liên quan đến ADN và protein.

V. Kết luận và tương lai của nghiên cứu tĩnh điện trong sinh học

Nghiên cứu mô phỏng tĩnh điện trong sinh học đang ngày càng trở nên quan trọng. Với sự phát triển của công nghệ và các phương pháp mô phỏng mới, tương lai của lĩnh vực này hứa hẹn sẽ mang lại nhiều khám phá mới và ứng dụng thực tiễn trong y học và công nghệ sinh học.

5.1. Xu hướng nghiên cứu trong tương lai

Các nghiên cứu trong tương lai sẽ tập trung vào việc cải thiện độ chính xác của các mô hình mô phỏng và mở rộng ứng dụng của chúng trong các lĩnh vực khác nhau của sinh học.

5.2. Tác động của công nghệ mới

Sự phát triển của công nghệ tính toán và mô phỏng sẽ mở ra nhiều cơ hội mới cho nghiên cứu tĩnh điện, từ đó giúp giải quyết các vấn đề phức tạp trong sinh học.

Nội dung chính

Tổng quan nghiên cứu

Phân tử axit deoxyribonucleic (ADN) là trung tâm lưu trữ thông tin di truyền của mọi dạng sống, mang khoảng 3 tỷ cặp nucleotide ở tế bào người với tổng chiều dài xấp xỉ 1 m nhưng được đóng gói kỳ diệu trong nhân tế bào chỉ vài micromet. Quá trình cuộn xoắn đa cấp này giúp kích thước phân tử giảm xấp xỉ 1.000.000 lần so với trạng thái duỗi thẳng. Khi tồn tại trong môi trường nước sinh lý, các gốc phosphate mất proton khiến ADN trở thành một polymer tích điện âm cực mạnh với mật độ điện tích tuyến tính đạt khoảng -1e/1,7 Å và mật độ điện tích bề mặt hình trụ xấp xỉ 1e/1 nm2. Đây là một trong những mật độ điện tích cao nhất từng được ghi nhận trong các hệ thống sinh học, làm cho tương tác tĩnh điện Coulomb trở thành yếu tố quyết định độ cứng, độ đàn hồi và năng lượng tự do đóng gói của sợi ADN.

Vấn đề cốt lõi đặt ra trong vật lý polymer và lý sinh học là xác định quy luật phụ thuộc của độ dài quán tính tĩnh điện ($l_e$) và khoảng cách đầu - cuối ($R_{ee}$) vào bán kính chắn Debye - Hückel ($r_s$) của dung môi. Suốt nhiều thập kỷ, cộng đồng khoa học tồn tại tranh cãi gay gắt giữa lý thuyết OSF của Odijk, Skolnick và Fixman ($l_e$ tỷ lệ với $r_s^2$) với các lý thuyết trường trung bình của nhóm Joanny ($l_e$ tỷ lệ với $r_s^1$). Đề tài được thực hiện bởi học viên Nguyễn Thị Trang dưới sự hướng dẫn của Tiến sĩ Nguyễn Thế Toàn tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội năm 2015, với sự hỗ trợ từ Quỹ NAFOSTED qua đề tài số 103.75. Nghiên cứu giải quyết bài toán bằng cách thực hiện mô phỏng Monte Carlo quy mô lớn lên đến 4.096 monomer, mở rộng gấp 4 lần giới hạn các nghiên cứu trước đó để loại bỏ triệt để hiệu ứng kích thước hữu hạn, từ đó làm sáng tỏ bản chất tương tác tĩnh điện trong vật lý sinh học.

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 xây dựng trên nền tảng phương trình trường trung bình tự hợp Poisson - Boltzmann (PB), mô tả thế năng tĩnh điện của các ion linh động quanh đại phân tử sinh học. Trong điều kiện thế năng tĩnh điện nhỏ hơn năng lượng nhiệt ($k_B T$), phương trình PB được tuyến tính hóa thành phương trình Debye - Hückel (DH), cho nghiệm dạng thế năng Yukawa bị chắn $V(r) = (e^2 / \varepsilon r) \exp(-r/r_s)$, trong đó $\varepsilon \approx 78$ là hằng số điện môi của nước tại nhiệt độ phòng 298 K. Bán kính chắn $r_s$ trong dung dịch muối đơn trị 1 mM đạt khoảng 100 Å, phản ánh mức độ suy giảm tương tác tĩnh điện theo hàm mũ khi nồng độ muối tăng.

Đối với các hệ tích điện cao, lý thuyết cô đọng Manning chỉ ra rằng khi tham số $\xi_m = l_B / b > 1$ (với độ dài Bjerrum $l_B \approx 7$ Å trong nước và $b$ là khoảng cách giữa các điện tích), các phản ion sẽ ngưng tụ lên bề mặt polymer để đưa mật độ điện tích hiệu dụng về ngưỡng tới hạn $e/l_B$. Độ dài quán tính tổng cộng của polymer bao gồm độ dài tự nhiên $l_0$ và đóng góp tĩnh điện $l_e$. Khung lý thuyết OSF và mô hình mở rộng Khokhlov - Khachaturian (KK) cho polymer dẻo dựa trên khái niệm chuỗi khối cầu tĩnh điện (blobs) với kích thước $\xi \approx (D k_B T l_0^2 / \eta_0^2)^{1/3}$, qua đó tái khẳng định đóng góp tĩnh điện tỷ lệ thuận với bình phương bán kính chắn $r_s^2$.

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

Mô phỏng sử dụng phương pháp Monte Carlo cổ điển kết hợp thuật toán Metropolis để lấy mẫu phân bố thống kê Boltzmann trong không gian pha cấu hình. Cỡ mẫu nghiên cứu bao gồm 7 hệ polymer tích điện với số lượng monomer $N$ lần lượt là 64, 128, 256, 512, 1024, 2048 và 4096 hạt. Độ dài liên kết giữa các monomer liền kề được cố định bằng độ dài Bjerrum $l = l_B \approx 7$ Å, tương ứng mật độ điện tích tuyến tính ngưỡng Manning.

Phương pháp chọn mẫu cấu hình áp dụng kỹ thuật kết hợp: 30% số bước thử dùng phép dịch chuyển quay khối toàn cục (pivot move) và 70% số bước thử dùng phép quay cục bộ (flip move quanh trục nối hai monomer lân cận). Tỷ lệ này tối ưu hóa không gian pha, duy trì xác suất chấp nhận cấu hình đạt khoảng 50%. Mỗi bộ tham số ($N, r_s$) được thực hiện từ 10.000 đến 20.000 bước Monte Carlo, trong đó 512 bước đầu tiên được loại bỏ để hệ đạt trạng thái cân bằng nhiệt hoàn toàn từ cấu hình que hoặc cầu Gaussian ban đầu.

Lý do lựa chọn phương pháp phân tích số trực tiếp là để khắc phục việc cắt cụt thế năng (cut-off), vốn thường gây sai số lớn ở bán kính chắn rộng. Toàn bộ tương tác Coulomb giữa mọi cặp hạt được tính toán đầy đủ với độ phức tạp tính toán $O(N^2)$ mỗi bước và thời gian hồi phục tỷ lệ với $N^3$, được tăng tốc song song hóa qua thư viện SimEngine viết bằng C++, OpenMP và OpenCL trên CPU đa lõi và GPU tính toán hiệu năng cao.

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

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

Mô phỏng mang lại ba kết quả định lượng quan trọng về cấu trúc không gian của dây polymer tích điện:

Thứ nhất, ở giới hạn chắn cực mạnh ($r_s \to 0$), chuỗi polymer tích điện hoàn toàn mất tương tác tĩnh điện tầm xa và hoạt động như một chuỗi Gaussian lý tưởng, trong đó bình phương khoảng cách đầu - cuối tuân theo hệ thức $R_{ee}^2 = l_B^2 (N - 1)$, hoàn toàn không phụ thuộc vào bán kính chắn.

Thứ hai, ở dải bán kính chắn rộng ($r_s / l_B \gg 1$), hệ số tỷ lệ $\alpha = \partial \ln(R_{ee}^2) / \partial \ln(r_s)$ biến thiên phụ thuộc rõ rệt vào số lượng monomer $N$. Với các chuỗi ngắn $N \le 512$, hệ số $\alpha$ duy trì ở mức thấp xấp xỉ 1. Tuy nhiên, khi tăng quy mô chuỗi lên $N = 2048$ và $N = 4096$, đường cong thực nghiệm của $\alpha$ tiệm cận chính xác với đường dự đoán của lý thuyết OSF ($l_e \sim r_s^2$), vượt xa mức dự đoán tuyến tính ($l_e \sim r_s$) của nhóm Joanny.

Thứ ba, hàm tương quan góc liên kết (BACF) dọc theo thân chuỗi với $N = 512$ và $r_s = 50 l_B$ cho thấy tương quan góc giảm dần theo khoảng cách và bảo toàn độ cứng cục bộ trên quy mô lớn hơn nhiều so với bán kính chắn, phù hợp với mật độ điện tích chuẩn hóa hội tụ tại các giá trị $N$ lớn.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự bất đồng giữa các nghiên cứu trước đây nằm ở hiệu ứng kích thước hữu hạn (finite-size effects). Khi các mô phỏng máy tính trước đó chỉ dừng lại ở quy mô $N \le 1024$ monomer, chiều dài toàn phần của chuỗi chưa đủ lớn so với độ dài quán tính tĩnh điện $l_e$, làm cho hiệu ứng tự do ở hai đầu chuỗi chiếm ưu thế và tạo ra cảm giác giả tạo rằng $l_e$ chỉ tỷ lệ bậc 1 với $r_s$. Khi nghiên cứu này nâng kích thước chuỗi lên mức kỷ lục 4.096 hạt, hiệu ứng biên bị triệt tiêu hoàn toàn, qua đó khẳng định lý thuyết OSF và mô hình blob KK là hoàn toàn chính xác.

Dữ liệu mô phỏng được thể hiện rõ ràng qua đồ thị logarit kép giữa $R_{ee}^2 / l_B^2$ theo $r_s / l_B$, minh họa sự tách biệt rõ rệt giữa các nhánh kích thước hạt khác nhau, kèm theo bảng đối chiếu độ dốc $\alpha$ tiệm cận giá trị lý thuyết khi $r_s$ tiến tới vô cùng. Điều này mang ý nghĩa to lớn trong việc giải thích cơ chế tương tác tĩnh điện điều hòa độ uốn cong của ADN khi kết hợp với các protein histon trong nhân tế bào sinh vật.

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

Nhằm phát triển sâu rộng các ứng dụng của mô phỏng tĩnh điện sinh học, luận văn đề xuất 4 khuyến nghị hành động cụ thể:

Thứ nhất, nâng cấp thuật toán mô phỏng song song trên kiến trúc GPGPU tiên tiến để nâng quy mô chuỗi lên $N \ge 10.000$ monomer trong vòng 12 tháng tới, do nhóm nghiên cứu vật lý lý thuyết và tin sinh học chủ trì, nhằm mô phỏng trọn vẹn cấu trúc một vòng chromatin hoàn chỉnh.

Thứ hai, tích hợp thêm các tương tác tầm ngắn thực tế như tương tác thể tích loại trừ Lennard - Jones và liên kết hydro vào mô hình hạt thô trong lộ trình 6 tháng tiếp theo, nhằm giảm sai số năng lượng tự do xuống dưới 3% khi mô tả các dung dịch sinh học phức tạp.

Thứ ba, chuyển giao và chuẩn hóa thư viện SimEngine thành gói phần mềm mô phỏng nguồn mở cho cộng đồng nghiên cứu trong nước trước quý IV/2026, hỗ trợ nâng cao 50% hiệu suất xử lý dữ liệu cho các bài toán đại phân tử sinh học.

Thứ tư, khuyến nghị các phòng thí nghiệm sinh học thực nghiệm sử dụng kỹ thuật bẫy quang học (optical tweezers) thiết lập các dải nồng độ muối monovalent từ 0,1 mM đến 100 mM để kiểm chứng thực nghiệm tính co giãn của ADN dựa trên các điểm uốn lý thuyết đã được mô phỏng.

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

Nội dung và phương pháp của luận văn mang lại giá trị học thuật và ứng dụng thiết thực cho 4 nhóm đối tượng:

1. Học viên cao học và nghiên cứu sinh ngành Vật lý lý thuyết, Vật lý toán và Lý sinh: Tài liệu cung cấp cơ sở toán học chặt chẽ về phương trình Poisson - Boltzmann, lý thuyết Debye - Hückel và kỹ thuật lập trình mô phỏng Monte Carlo nâng cao.
2. Các nhà nghiên cứu Sinh học phân tử và Công nghệ gen: Giúp hiểu rõ cơ chế vật lý chi phối quá trình đóng gói ADN, cấu trúc cuộn xoắn nucleosome và các tương tác tĩnh điện giữa axit nucleic với protein chức năng.
3. Kỹ sư Khoa học vật liệu và Công nghệ nano: Ứng dụng quy luật tĩnh điện của polymer để thiết kế các hệ dẫn thuốc thông minh, vật liệu sinh học tự lắp ghép và polymer dẫn điện có cấu trúc kiểm soát.
4. Chuyên gia phát triển phần mềm tính toán khoa học: Tham khảo mô hình tối ưu hóa tính toán nhiều hạt không cắt cụt thế năng bằng C++, OpenMP và OpenCL trên phần cứng hiệu năng cao.

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

1. Độ dài quán tính tĩnh điện của ADN thay đổi như thế nào khi nồng độ muối trong môi trường tăng lên? Khi nồng độ muối tăng, mật độ ion linh động tăng khiến bán kính chắn Debye - Hückel $r_s$ giảm tỷ lệ nghịch với căn bậc hai nồng độ ion. Do đó, lực đẩy tĩnh điện giữa các monomer bị suy giảm mạnh, làm độ dài quán tính tĩnh điện $l_e$ giảm theo quy luật bình phương bán kính chắn ($l_e \sim r_s^2$), khiến chuỗi ADN trở nên mềm dẻo hơn.

2. Tại sao nghiên cứu kết hợp 30% phép dịch chuyển pivot và 70% phép dịch chuyển flip trong mô phỏng? Phép dịch chuyển pivot giúp xoay các đoạn chuỗi lớn ngẫu nhiên để thay đổi nhanh cấu hình toàn cục, trong khi phép flip giúp tinh chỉnh các góc quay monomer cục bộ quanh trục nối. Tỷ lệ 30% pivot và 70% flip giúp cân bằng tối ưu giữa tốc độ quét không gian pha và xác suất chấp nhận cấu hình đạt mức lý tưởng khoảng 50%.

3. Hiện tượng cô đọng Manning ảnh hưởng như thế nào đến mật độ điện tích của chuỗi polymer? Khi mật độ điện tích của polymer vượt quá ngưỡng tới hạn khiến tham số Manning $\xi_m > 1$, lực hút tĩnh điện sẽ giữ các phản ion ngưng tụ quanh thân chuỗi. Quá trình này làm trung hòa bớt điện tích bề mặt và đưa mật độ điện tích hiệu dụng về hằng số $e/l_B$ bất kể mật độ điện tích ban đầu cao đến mức nào.

4. Điểm khác biệt mấu chốt giữa lý thuyết OSF và lý thuyết Joanny là gì? Lý thuyết OSF dự đoán độ dài quán tính tĩnh điện tỷ lệ với bình phương bán kính chắn ($l_e \sim r_s^2$), áp dụng tốt cho chuỗi polymer cứng và bán dẻo. Ngược lại, lý thuyết Joanny dự đoán tỷ lệ tuyến tính ($l_e \sim r_s^1$). Kết quả mô phỏng quy mô 4.096 monomer trong luận văn đã chứng minh lý thuyết OSF là hoàn toàn chính xác.

5. Vì sao các nghiên cứu trước đây lại đưa ra kết luận sai lệch về quy luật phụ thuộc của $l_e$ vào $r_s$? Các nghiên cứu trước đây bị giới hạn năng lực phần cứng nên chỉ mô phỏng chuỗi polymer ngắn dưới 1.024 monomer. Khi chuỗi quá ngắn so với độ dài tĩnh điện, hiệu ứng tự do ở hai đầu chuỗi làm biến dạng cấu trúc toàn cục, dẫn đến số mũ đo được bị kéo thấp xuống mức xấp xỉ 1 thay vì giá trị thực là 2.

Kết luận

• Công trình khẳng định vai trò quyết định của tương tác tĩnh điện Coulomb chắn Debye - Hückel đối với độ đàn hồi và cấu trúc không gian của phân tử ADN trong dung môi nước.
• Dữ liệu mô phỏng quy mô lớn xác nhận tính chính xác tuyệt đối của lý thuyết OSF và mô hình chuỗi blob Khokhlov - Khachaturian với quy luật tỷ lệ $l_e \sim r_s^2$.
• Nghiên cứu giải quyết dứt điểm cuộc tranh luận học thuật kéo dài nhiều năm nhờ mở rộng quy mô hệ lên đến 4.096 monomer, loại bỏ hoàn toàn các sai số do hiệu ứng kích thước hữu hạn.
• Xây dựng thành công quy trình mô phỏng Monte Carlo tối ưu hóa toàn diện trên CPU đa lõi và GPU qua thư viện SimEngine.
• Cung cấp nền tảng định lượng vững chắc cho việc giải mã cơ chế đóng gói thông tin di truyền và phát triển các vật liệu polymer sinh học thế hệ mới.

Trong 12 đến 24 tháng tới, các nhóm nghiên cứu nên tiếp tục mở rộng mô hình sang các hệ sinh học đa thành phần chứa ion đa hóa trị và protein liên kết. Hãy ứng dụng ngay các phát hiện và khung phương pháp mô phỏng của luận văn này để nâng cao độ chính xác cho các công trình nghiên cứu lý sinh và vật liệu nano của bạn.

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

Chương 1 của luận văn, chúng tôi sẽ giới thiệu những kiến thức cơ bản nhất về cấu trúc hóa học, cấu trúc không gian và một số tính chất vật lý của phân tử ADN chịu ảnh hưởng của các hiệu ững tĩnh điện.1 Cấu trúc hóa học của phân tử ADN Về thành phần hóa học, ADN là một loại axit hữu cơ có chứa các nguyên tố chủ yếu là: cacbon (C), hiđro (H), oxi (O), photpho (P), nitơ (N). Về cấu trúc hóa học, ADN là đại phân tử có khối lượng phân tử lớn, chiều dài có thể đạt tới hàng trăm micromet, khối lượng phân tử có thể đạt từ 4 đến 8 triệu, một số có thể đạt tới 16 triệu đơn vị cacbon. ADN là hợp chất cao phân tử cấu tạo theo nguyên tắc đa phân, trong đó mỗi đơn phân là các nucleotide gọi tắt là các Nu. Trong cơ thể người, số nucleotide có thể lên tới 3 tỉ nu với chiều dài phân tử ADN xấp xỉ 1m.

Về thành phần hóa học, mỗi nucleotide có 3 thành phần cơ bản: một phân tử đường deoxiribose, một gốc phosphate (H3PO4) và một trong bốn loại base 5 LUAN VAN CHAT LUONG download : add luanvanchat@agmail. Trong đó hai loại base Cytosine (C) và Thymine (T) có cấu trúc vòng đơn cacbon nitrogen, kích thước nhỏ (Pyrimidine), hai loại base Adenin (A) và Guanin (G) có cấu trúc vòng kép, kích thước lớn hơn (Purine) (Hình 1. Do mỗi nucleotide chỉ khác nhau ở thành phần nucleobase nên tên gọi của các nu là tên gọi của nucleobase mà nó mang.1 Cấu trúc 4 loại base của phân tử ADN (Nguồn:http://pixshark.com) Cấu trúc không gian phổ biến của ADN là cấu trúc chuỗi xoắn kép gồm hai mạch đơn (sẽ được chúng tôi giới thiệu chi tiết trong nội dung tiếp theo). Với cấu trúc 6 LUAN VAN CHAT LUONG download : add [e-mail] này, trên mỗi mạch đơn các đơn phân liên kết với nhau bằng liên kết hóa trị, là liên kết hình thành giữa đường C5H10O4 của nucleotide này với phân tử H3PO4 của nucleotide bên cạnh (hay còn được gọi là liên kết photphodieste) hình thành nên chuỗi polinucleotide (Hình 1.2 Cấu trúc hóa học của ADN (Nguồn: http://commons.org) Liên kết photphodieste là liên kết bền, đảm bảo thông tin di truyền trên mỗi mạch đơn ổn định ngay cả khi ADN tái bản và phiên mã.

Đồng thời chính liên kết hóa trị này cũng làm nên tính ổn định của phân tử ADN. Do đó chuỗi các liên kết giữa phân tử đường và gốc phosphate được xem là chuỗi “xương sống” của phân tử ADN, có ý nghĩa quan trọng trong hoạt động di truyển của các cơ thể sống. Từ chuỗi xương sống này, các nucleobase trên nhánh đơn gắn vào nhóm phân tử đường của nhánh đó như “chân rết”, trên hai nhánh đơn các nuleobase liên kết với nhau thông qua các liên kết hiđro. Hai liên kết này hình thành nên độ cứng và cấu trúc không gian của phân tử ADN.

7 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.2 Cấu trúc không gian của phân tử ADN 1.1 Cấu trúc sơ cấp Cấu trúc sơ cấp của phân tử ADN được biết đến là cấu trúc chuỗi xoắn kép gồm hai sợi đơn - chuỗi polinucleotide. Có nhiều dạng xoắn hai mạch đơn của phân tử ADN, tuy nhiên cấu trúc phổ biến 1nm nhất là cấu trúc xoắn dạng B – ADN theo mô hình của J.cric xây dựng năm 1953.4 hình này, hai mạch đơn của phân tử ADN liên kết với nm nhau và xoắn đều quanh một trục tưởng tượng theo chiều từ trái sang phải tương tự như chiếc thang dây xoắn với bán kính 1nm gồm nhiều vòng xoắn. Các chu kỳ xoắn được lặp đi lặp lại đều đặn, chiều cao mỗi vòng xoắn là 3.4𝑛𝑚, trung bình khoảng 10.5 cặp base nu một chu kỳ với khoảng cách các cặp nu xấp xỉ 0.