Động lực học phân tử (molecular dynamics - MD) là phương pháp mô phỏng máy tính dựa trên cơ học thống kê và cơ học cổ điển nhằm nghiên cứu sự tiến triển theo thời gian của một hệ thống nhiều hạt (nguyên tử, phân tử hoặc hạt thô). Bằng cách tích phân số các phương trình chuyển động của Newton cho toàn bộ nguyên tử trong hệ, mô phỏng MD cho phép quan sát chi tiết động học phân tử ở độ phân giải không gian cấp độ angstrom (0.1 nanomet) và độ phân giải thời gian cấp độ femtosecond (1 fs = giây).
Cơ sở lý thuyết và phương trình chuyển động
Trong mô phỏng động lực học phân tử cổ điển, mỗi nguyên tử có khối lượng được coi là một chất điểm tại vị trí . Lực tương tác tác dụng lên nguyên tử tại thời điểm được xác định bởi đạo hàm riêng bậc nhất của thế năng toàn phần :
\mathbf{F}_i = - abla_i V(\mathbf{r}_1, \mathbf{r}_2, \dots, \mathbf{r}_N) = m_i rac{d^2 \mathbf{r}_i}{dt^2}
Hệ phương trình vi phân cấp 2 liên kết này không thể giải chính xác bằng giải tích đối với hệ nhiều hạt (), do đó cần sử dụng các thuật toán tích phân số như thuật toán Verlet, Velocity Verlet hoặc Leap-frog với bước thời gian tích phân (time step) điển hình từ 1 đến 2 femtoseconds nhằm đảm bảo tính bảo toàn năng lượng và độ ổn định số học.
Trường lực trong mô phỏng phân tử
Thế năng toàn phần được mô tả thông qua trường lực (force field) kinh nghiệm, chia thành hai nhóm đóng góp chính: tương tác liên kết (bonded interactions) và tương tác phi liên kết (non-bonded interactions):
V(\mathbf{r}) = \sum_{ ext{bonds}} K_b (r - r_0)^2 + \sum_{ ext{angles}} K_ heta ( heta - heta_0)^2 + \sum_{ ext{dihedrals}} V_n [1 + \cos(n\phi - \gamma)] + \sum_{i < j} \left[ 4 arepsilon_{ij} \left( \left( rac{\sigma_{ij}}{r_{ij}} ight)^{12} - \left( rac{\sigma_{ij}}{r_{ij}} ight)^6 ight) + rac{q_i q_j}{4\pi arepsilon_0 r_{ij}} ight]
- Tương tác liên kết: Gồm dao động dãn dài liên kết (harmonic bond stretching), góc hóa trị (angle bending) và năng lượng xoay quanh trục liên kết (torsional dihedral angles).
- Tương tác phi liên kết: Gồm thế tương tác Van der Waals biểu diễn qua hàm thế Lennard-Jones 12-6 và tương tác tĩnh điện Coulomb giữa các điện tích điểm cục bộ .
Các bộ trường lực phổ biến trong mô phỏng phân tử sinh học bao gồm AMBER, CHARMM, GROMOS và OPLS-AA, được tối ưu hóa dựa trên tính toán hóa học lượng tử ab initio kết hợp thực nghiệm phổ học.
Tập hợp thống kê và kiểm soát nhiệt động
Để tái lập chính xác các điều kiện nhiệt động lực học trong phòng thí nghiệm, mô phỏng MD áp dụng các tập hợp thống kê khác nhau:
- Tập hợp vi chính tắc (NVE): Bảo toàn số hạt , thể tích và tổng năng lượng .
- Tập hợp chính tắc (NVT): Giữ cố định nhiệt độ thông qua các bộ điều nhiệt (thermostat) như Berendsen, Langevin dynamics hoặc Nosé-Hoover thermostat.
- Tập hợp đẳng nhiệt - đẳng áp (NPT): Giữ cố định nhiệt độ và áp suất (thường ở 300 K và 1 bar) bằng cách kết hợp bộ điều áp (barostat) như Berendsen hoặc Parrinello-Rahman barostat.
Điều kiện biên và xử lý tương tác tầm xa
Để triệt tiêu hiệu ứng bề mặt nhân tạo của hộp mô phỏng hữu hạn, điều kiện biên tuần hoàn (Periodic Boundary Conditions - PBC) được thiết lập bằng cách sao chép vô hạn hộp mô phỏng trong không gian ba chiều. Tương tác tĩnh điện tầm xa được tính toán hiệu quả và chính xác bằng thuật toán tổng Ewald lưới hạt (Particle Mesh Ewald - PME), giảm độ phức tạp tính toán từ xuống .
Ứng dụng và phương pháp lấy mẫu nâng cao
Động lực học phân tử là công cụ then chốt trong sinh học cấu trúc, hóa dược và khoa học vật liệu: từ giải mã cơ chế gấp cuộn protein (protein folding), vận chuyển ion qua kênh màng tế bào, đến tính toán năng lượng liên kết tự do của thuốc với protein đích (free energy perturbation - FEP, MM-PBSA). Khi nghiên cứu các quá trình chuyển pha hoặc thay đổi cấu hình diễn ra ở thang thời gian mili-giây đến giây vượt quá khả năng của MD thông thường, các kỹ thuật lấy mẫu nâng cao (enhanced sampling) như Metadynamics, Umbrella Sampling và Replica Exchange Molecular Dynamics (REMD) được ứng dụng rộng rãi.