Xây dựng mô hình toán học dự đoán khả năng thẩm thấu qua da dựa trên cấu trúc phân tử​

VN
U

ĐẠI HỌC QUỐC GIA HÀ NỘI

ha
r

ma
c

y,

KHOA Y DƢỢC

ne
an
dP

Người thực hiện: ĐỖ TRẦN LONG

ici

XÂY DỰNG MƠ HÌNH TỐN HỌC

ho
ol
o

fM

ed

DỰ ĐỐN KHẢ NĂNG THẨM THẤU QUA DA
DỰA TRÊN CẤU TRÚC PHÂN TỬ

Sc

KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC NGÀNH DƢỢC HỌC

Khóa: QH.2013.Y

Co
p

yri

gh
t@

Người hướng dẫn: TS. LÊ THỊ THU HƢỜNG

Hà Nội - 2018


VN
U

LỜI CẢM ƠN

y,

Lời đầu tiên, em xin được gửi lời cảm ơn tới Ban giám đốc Đại học
Quốc Gia Hà Nội, các phòng ban đã cho phép em được thực hiện khóa luận
tốt nghiệp.

ha
r

ma
c

Em xin được cảm ơn các thầy giáo, cô giáo khoa Y Dược Đại học Quốc
Gia Hà Nội đã truyền đạt cho em những kiến thức, kĩ năng để trở thành một
dược sĩ.

ici

ne
an
dP

Bằng cả tấm lòng, em xin gửi lời cảm ơn chân thành tới TS. Lê Thị
Thu Hƣờng (Bộ môn Dược liệu – Khoa Y dược Đại học Quốc Gia Hà Nội),
người đã dìu dắt em vào con đường nghiên cứu khoa học; đã tận tình chỉ bảo,
truyền đạt kinh nghiệm nghiên cứu, kinh nghiệm sống cho em; đã khơng quản
ngại khó khăn, dồn tâm huyết, sức lực để tôi luyện, rèn giũa cho em kĩ năng
và đạo đức nghiên cứu dù em còn nhiều thiếu sót, nhiều lần vấp ngã. Cơ là
tấm gương sáng cho thế hệ học trò chúng em học tập.

fM

ed

Cuối cùng, con xin cảm ơn bố mẹ đã sinh thành, nuôi dạy con nên
người; cảm ơn mẹ đáng kính đã ln bên con, ủng hộ con; cảm ơn chị gái và
bạn bè đã luôn dõi theo em trên mỗi chặng đường.

Co
p

yri

gh
t@

Sc

ho
ol
o

Em xin chân thành cảm ơn!

Hà Nội, ngày 9 tháng 5 năm 2018
Sinh Viên

Đỗ Trần Long


VN

U

DANH MỤC KÝ HIỆU VIẾT TẮT

Tên Tiếng Anh

ANN

Artificial
Networks

Tên Tiếng Việt

y,

Tên viết tắt

ma
c

Neural Mạng nơ-ron nhân tạo

In vitro release testing

Kow

Hệ số phân chia octanolnước

Kp

Hệ số thẩm thấu qua da

Living skin equivalent Mơ hình tương đương
models
với da người sống

MLR

Multiple
Regression

fM

ed

ici

LSEs

PCA

Principal
Analysis

gh
t@

Q2
Q2ext

yri

QSAR
R2

Đường kính ngồi

Component Phân tích thành phần
chính

Polydimethylsiloxan

Sc

PDMS
Q

Linear Phương pháp hồi quy
tuyến tính đa biến

Outside diameter

ho
ol
o

OD

Co
p

Thử nghiệm giải phóng
thuốc in vitro

ne
an
dP

IVRT

ha
r

Cơ sở dữ liệu

CSDL

Hệ số tương quan chéo
Hệ số xác định cho tập
kiểm tra
Quantitative Structure – Tương quan định lượng
Activity Relationship
cấu trúc – tác dụng
Hệ số xác định


Test set

Tập kiểm tra

Tr

Training set

Tập huấn luyện

Tham số phân tử

Co
p

yri

gh
t@

Sc

ho
ol
o

fM

ed

ici

ne
an

dP

ha
r

ma
c

y,

TSPT

VN
U

Te


VN
U

DANH MỤC HÌNH VẼ

Hình 1.1. Thiết bị Franz và chi tiết các bộ phận ............................................... 2

y,

Hình 1.2. Cấu tạo gồm hai phần chính của thiết bị Franz................................. 3

ma

c

Hình 2.1. Các bước xây dựng mơ hình QSAR ............................................... 20

Co
p

yri

gh
t@

Sc

ho
ol
o

fM

ed

ici

ne
an
dP

ha
r

Hình 3.1. Miền ứng dụng của mơ hình QSAR đối với các chất ..................... 26


VN
U

DANH MỤC BẢNG

Bảng 1.1. Thuật tốn tính hệ số thấm (Kp) được công bố bởi Flynn (1990) .. 10

y,

Bảng 1.2. Các nghiên cứu về khả năng thẩm thấu qua da của các hợp chất

ma
c

được công bố sau 1990.................................................................................... 11
Bảng 2.1. Giá trị log Kp các hợp chất xây dựng mô hình QSAR ................... 18

Co
p

yri

gh
t@

Sc

ho
ol
o

fM

ed

ici

ne
an
dP

ha
r

Bảng 3.1. Thơng số thống kê của các mơ hình QSAR log Kp ........................ 24


VN
U

MỤC LỤC
DANH MỤC KÝ HIỆU VIẾT TẮT

y,

DANH MỤC HÌNH VẼ

ma
c

DANH MỤC BẢNG

ĐẶT VẤN ĐỀ .................................................................................................. 1

ha
r

CHƢƠNG 1 - TỔNG QUAN.......................................................................... 2
1.1. Thiết bị Franz .......................................................................................... 2

ne
an
dP

1.1.1. Vai trò của thiết bị Franz trong thí nghiệm đo khả năng thẩm thấu
qua da.......................................................................................................... 2
1.1.2. Cấu tạo của thiết bị Franz ................................................................ 3
1.1.3. Cách hoạt động của thiết bị Franz ................................................... 4
1.2. Tính thấm qua da..................................................................................... 4

ici

1.2.1. Khó khăn khi nghiên cứu khả năng thẩm thấu qua da ..................... 4

ed

1.2.2. Các phương pháp đánh giá khả năng thẩm thấu của các chất hóa
học qua da ................................................................................................... 5

fM

1.2.3. Mối quan hệ định lượng cấu trúc-tác dụng (QSAR) và khả năng
thẩm thấu qua da ........................................................................................ 7

ho
ol
o

1.2.4. Quá trình phát triển của QSAR trong nghiên cứu tính thấm qua da...... 8
1.3. Tổng quan về mối quan hệ định lượng cấu trúc-tác dụng của hợp chất
(QSAR) ........................................................................................................ 11
1.3.1. Khái niệm và nguyên lý chung của QSAR ...................................... 11

Sc

1.3.2. Xây dựng và đánh giá mơ hình QSAR ........................................... 13

gh
t@

CHƢƠNG 2 - NGUYÊN LIỆU VÀ PHƢƠNG PHÁP NGHIÊN CỨU ... 18
2.1. Nguyên liệu ........................................................................................... 18
2.1.1. Cơ sở dữ liệu ................................................................................... 18
2.1.2. Phần mềm sử dụng .......................................................................... 19

Co

p

yri

2.2. Phương pháp nghiên cứu ..................................................................... 19
2.2.1 Tính tốn tham số mơ tả phân tử ..................................................... 19
2.2.2. Phân chia tập huấn luyện / Tập kiểm tra........................................ 20


VN
U

2.2.3. Xây dựng mơ hình QSAR ................................................................ 20
2.2.4. Đánh giá mơ hình ........................................................................... 22
CHƢƠNG 3 – KẾT QUẢ VÀ BÀN LUẬN ................................................ 23

y,

3.1. Kết quả ................................................................................................. 23

ma
c

3.1.1. Mơ hình tốn học thu được ............................................................. 23
3.1.2. Đánh giá mơ hình ........................................................................... 24

ha
r

3.1.3. Đánh giá mơ hình theo 5 tiêu chí của OECD ................................. 25

3.2. Bàn luận .............................................................................................. 27

ne
an
dP

3.2.1. Về phương pháp .............................................................................. 27
3.2.2. Về kết quả của thu được.................................................................. 28
3.2.3. So sánh với các mơ hình khác trên thế giới .................................... 30
KẾT LUẬN VÀ KIẾN NGHỊ ...................................................................... 31

Co
p

yri

gh
t@

Sc

ho
ol
o

fM

ed

ici

TÀI LIỆU THAM KHẢO ............................................................................ 32


ne
an
dP

ha
r

ma
c

y,

VN
U

ĐẶT VẤN ĐỀ
Da là cơ quan có diện tích lớn nhất của cơ thể con người (bao phủ 1,52
2,3m diện tích bề mặt và chiếm khoảng 16% trọng lượng cơ thể). Nó có cấu
tạo phức tạp nhiều lớp, bao gồm: biểu bì, bì và lớp mỡ dưới da cùng các thành
phần phụ (lơng, móng, tuyến bã, tuyến mồ hơi) hình thành hệ da. Da người
đóng vai trị như một lớp bảo vệ cơ thể khỏi môi trường xung quanh, với chức
năng chính là điều chỉnh sự xâm nhập của các chất lạ vào cơ thể. Vai trò là
một hàng rào bảo vệ này của da đã thu hút sự quan tâm rất lớn trong giới khoa
học vì thơng qua đó, các nhà khoa học có thể phát triển nhiều nghiên cứu có
tính ứng dụng cao [14].

ed

ici

Khả năng thẩm thấu của các hợp chất qua da đã được khai thác từ khá
lâu trước đây, ví dụ như phát triển thuốc bơi ngồi da, miếng dán ngồi da....
Ngồi ra, tính thấm qua da cũng dự đốn các độc tính có thể có của hoá chất
hay mỹ phẩm tiếp xúc với da người. Để hỗ trợ cho việc phát triển các thuốc
dùng ngoài da cũng như nghiên cứu độc chất và mỹ phẩm, nhiều mơ hình da
người đã được đề xuất, trong đó mơ hình thiết bị Franz được sử dụng phổ
biến nhất.

gh
t@

Sc

ho
ol
o

fM

Tuy nhiên mơ hình này triển khai khơng đơn giản, địi hỏi tốn kém thời
gian và chi phí cũng như hố chất máy móc. Do đó từ hơn một thập kỷ qua,
giới khoa học sử dụng cơng cụ tính tốn để dự đoán độ hấp thu qua da của
hoá chất. Một trong các cách tiếp cận được khai thác nhiều là sử dụng mơ
hình liên quan định lượng giữa cấu trúc hố học với độ thẩm thấu qua da. Các
mơ hình toán học biểu diễn quan hệ định lượng cấu trúc-tác dụng QSAR
(Quantitative structure-activity relationship) là giải pháp ứng dụng hiệu quả

để giải quyết vấn đề này. Qua đó, mục tiêu chung của nghiên cứu này là “Xây
dựng mơ hình tốn học dự đoán khả năng thẩm thấu qua da dựa trên cấu trúc
phân tử”. Với mục tiêu cụ thể:

Co
p

yri

1. Xây dựng mơ hình QSAR nhằm dự đốn khả năng thẩm thấu qua da
thơng qua của các hợp chất.
2. Phân tích và đánh giá mơ hình tìm được theo 5 tiêu chí của tố chức Hợp
tác và Phát triền Kinh tế thế giới (OECD).

1


VN
U

CHƢƠNG 1 - TỔNG QUAN
1.1. Thiết bị Franz

ma
c

y,

1.1.1. Vai trò của thiết bị Franz trong thí nghiệm đo khả năng thẩm thấu
qua da

ne
an
dP

ha
r

Thử nghiệm đo lường sự giải phóng thuốc từ một dạng thuốc được cho
là điều cơ bản để phát triển các sản phẩm thuốc. Xác định cụ thể kỹ thuật
được sử dụng bởi chính dạng liều lượng và đường phân phối mong muốn của
thuốc. Đối với các dạng thuốc rắn, thử nghiệm về độ tan của thuốc đã được sử
dụng trong 50 năm qua. Gần đây, đối với sản phẩm thuốc bán rắn, thử nghiệm
giải phóng thuốc in vitro (IVRT) là một phương tiện để đánh giá tính giải
phóng của thuốc. Tỷ lệ giải phóng in vitro có thể phản ánh hiệu quả kết hợp
của một số thông số vật lý và hóa học, bao gồm độ hịa tan và kích thước hạt
và các đặc tính lưu biến của dạng liều [18].

Co
p

yri

gh
t@

Sc

ho
ol

o

fM

ed

ici

Phương pháp IVRT phổ biến nhất sử dụng thiết kế buồng mở như hệ
thống thiết bị khuếch tán Franz và có thể được sử dụng với màng tổng hợp,
cấu trúc mô hoặc mẫu sinh học, chẳng hạn như da người [18].

Hình 1.1. Thiết bị Franz và chi tiết các bộ phận

2


VN
U

1.1.2. Cấu tạo của thiết bị Franz

ici

ne
an
dP

ha
r

ma
c

y,

Thiết bị Franz là các thiết bị giải phóng riêng biệt được làm từ hai phần
thủy tinh borosilicate. Phần trên có thể được gọi là nắp thiết bị, đỉnh của thiết
bị, buồng cấp, hoặc khoang cấp. Phần dưới thường được gọi là phần thân của
thiết bị, hay cịn được gọi là buồng thụ thể.

ed

Hình 1.2. Cấu tạo gồm hai phần chính của thiết bị Franz

fM

Buồng trên hay buồng cấp được nối với buồng dưới hay còn là buồng
thụ thể bằng một khớp. Màng ngăn qua đó thẩm thấu hoặc vận chuyển chất
được đặt ở giữa khớp và giữ bằng một chiếc kẹp.

Sc

ho
ol
o

Các thiết bị có vỏ bọc có thể được nối với nguồn nước để kiểm sốt
nhiệt độ. Các thiết bị khơng có vỏ bọc có thể được đặt trong máy sấy khơ, bồn
nước hoặc mơi trường được kiểm sốt nhiệt độ nếu cần. Mặc dù các thiết bị

có vỏ đắt hơn nhưng lại được sử dụng phổ biến hơn bởi vì chúng có thể dễ
dàng kiểm soát nhiệt độ.

gh
t@

Các thiết bị Franz thường sử dụng các máy khuấy V-Series. Nhưng nếu
chỉ sử dụng một thiết bị Franz đơn lẻ trong thí nghiệm in vitro thì cũng có thể
sử dụng máy khuấy thí nghiệm thơng thường.

Co
p

yri

Tay lấy mẫu, ống dẫn, ống thốt nước nóng, và lưu lượng được tùy
chỉnh của các thiết bị Franz tiêu chuẩn được làm từ ống OD 6mm có độ dày
thành 1,0mm. Dung sai của ống được nhà sản xuất ống đưa ra, đường kính
bên trong của ống là khoảng 4mm. Các xi lanh với có kích thước kim tiêm

3


VN
U

OD 2mm có thể dễ dàng được sử dụng để lấy mẫu thiết bị Franz. Lượng lấy
mẫu tiêu chuẩn của thiết bị Franz hiện chưa có quy chuẩn rõ ràng nên được
lấy một lượng càng nhỏ càng tốt [14].

y,

1.1.3. Cách hoạt động của thiết bị Franz

ho
ol
o

fM

ed

ici

ne
an
dP

ha
r

ma
c

Các sản phẩm được thử nghiệm được đặt ở buồng cấp. Màng ngăn tách
buồng cấp chứa mẫu thử nghiệm với buồng thụ thể chứa đầy các chất trung
gian. Muối Phosphate Buffered Saline (PBS) thường được lựa chọn đầu tiên
làm chất trung gian, mặc dù nó có thể khơng khả thi khi sử dụng trong một số
trường hợp [18]. Sự thẩm thấu của thuốc bán rắn qua màng được theo dõi
bằng cách lấy mẫu từ buồng thụ thể theo thời gian và kiểm tra. Tại các thời

điểm xác định trước, một lượng chất trung gian được lấy ra từ buồng thụ thể
để phân tích, thường là bằng HPLC. Chất trung gian sẽ được bù vào với mỗi
lần lấy mẫu như vậy. Buồng dưới chứa chất lỏng từ đó các mẫu được lấy theo
các khoảng thời gian đều đặn để phân tích. Thử nghiệm này nhằm xác định
lượng thẩm thấu qua màng tại mỗi thời điểm. Buồng được duy trì ở nhiệt độ
khơng đổi 37oC. Tùy thuộc vào phương tiện, khả năng thấm của từng chất
nhất định được xác định thơng qua phân tích thiết bị Franz có thể thay đổi
đáng kể (có lẽ từ 10 đến 50 lần). Kết quả thử nghiệm trên thiết bị Franz khơng
nhất thiết phải dự đốn hiệu quả của một cấu trúc; mà có thể là xác định liệu
cấu trúc đó có tác nhân hoạt động qua da hay khơng [24].
1.2. Tính thấm qua da

Sc

1.2.1. Khó khăn khi nghiên cứu khả năng thẩm thấu qua da

Co
p

yri

gh
t@

Rào cản chính đối với khả năng thẩm thấu qua da là lớp sừng, và sự
xâm nhập của một hợp chất được xác định chủ yếu bởi cấu trúc hóa học của
nó. Ngồi ra có nhiều yếu tố khác đóng vai trị trong việc ảnh hưởng đến sự
thẩm thấu qua da bao gồm (i) cách sử dụng hoặc tiếp xúc với các hợp chất
trên bề mặt da, (ii) các cách thức làm thay đổi tính chất rào cản của tầng lớp
sừng, (iii) quy trình xảy ra trong các mô, (iv) các yếu tố sinh học. Lớp sừng

nằm trên bề mặt ngoài của da, là một lớp tế bào chết chứa đầy keratin được
bao quanh bởi một hệ thống ngoại bào giàu lipid. Trong cơ thể người, nó bao
gồm 10 và 25 lớp tế bào chết, kéo dài, các tế bào sừng ở trong một hệ thống
4


ma
c

y,

VN
U

lipid kép. Lớp sừng chỉ dày 7–16 µm ở hầu hết các vùng của cơ thể nhưng
dày 400–600 µm ở lòng bàn tay và bàn chân. Tầng lớp sừng bao gồm 40%
protein trong đó 80% là keratin. Chủng loại và lượng lipid có trong lớp sừng
phụ thuộc vào từng vùng của cơ thể có thể làm ảnh hưởng lớn đến tính thấm
qua da [13].

ne
an
dP

ha
r

Nhiều nỗ lực đã được thực hiện để phát triển các phương pháp tiếp cận
thực nghiệm đo lường sự hấp thu qua da. Hầu hết các dữ liệu hiện nay đã thu
được bằng các thí nghiệm in vitro, trong khi dữ liệu in vivo thu được bằng các

xét nghiệm sinh học .Về tỷ lệ chuyển hóa qua da của các loại thuốc khác nhau
thì thường được nghiên cứu nhiều nhất bằng cách tiến hành các thí nghiệm
trong in vitro, đơn giản hơn nhiều so với các phép đo trên lâm sàng.
Mức độ thẩm thấu có thể được đo bằng :

Sự chuyển hóa thuốc trong động vật hoặc người thử nghiệm.
Sử dụng da được cắt bỏ từ người hoặc động vật.
Sử dụng mơ hình màng tổng hợp trong thí nghiệm in vitro.

ici

i.
ii.
iii.

ho
ol
o

fM

ed

Da người có từ nhiều nguồn khác nhau như từ phẫu thuật thẩm mỹ và
phẩu thuật cắt bớt, đã được sử dụng trong thí nghiệm in vitro. Da từ các loài
động vật bao gồm lợn, chuột, lợn, rắn, thỏ và các loài khác được lựa chọn
thay thế phù hợp cho da người. Trong số các loại, da của lợn và chuột thường
được sử dụng phổ biến nhất vì có nhiều điểm chung với da người. Mặc dù da
động vật cho thấy sự thẩm thấu các hóa chất trong nghiên cứu rất tốt, nhưng
nó khơng có sự phức tạp như lớp sừng ở da người.

gh
t@

Sc

1.2.2. Các phương pháp đánh giá khả năng thẩm thấu của các chất hóa
học qua da

Co
p

yri

Các phương pháp ex vivo và in vitro để thử nghiệm tính thẩm thấu của
da được thiết lập và công bố rộng rãi [2, 3]. Chúng được sử dụng phổ biến
trong nhiều lĩnh vực nghiên cứu. Da người từ nhiều nguồn khác nhau, bao
gồm phẫu thuật thẩm mỹ và phẩu thuật cắt bỏ đã được sử dụng để đánh giá ex
vivo về khả năng thẩm thấu qua da [2, 3, 8]. Việc sử dụng các loại da động vật
khác nhau cũng là một thành phần được chấp nhận phổ biến trong các nghiên
cứu thẩm thấu qua da trong thí nghiệm in vitro. Da từ các loài động vật bao
5


y,

VN
U

gồm lợn, chuột, khỉ, rắn và những loài khác đã được coi như là một sự thay

thế thích hợp cho da người [1, 4]. Trong đó, da từ lợn và chuột thường được
sử dụng rộng rãi nhất, vì chúng có các rào cản tương tự để thẩm thấu các phân
tử gần giống với da người [1, 4, 17].

ne
an
dP

ha
r

ma
c

Mặc dù da động vật có thể cung cấp các đặc tính rào cản và tạo ra các
đặc điểm thẩm thấu tương tự như da người, nhưng chúng khơng thể mang đầy
đủ tính chất phức tạp của da người nói chung, và đặc tính rào cản của lớp
sừng nói riêng. Vì vậy, một số nhà nghiên cứu đã nỗ lực phát triển da nhân
tạo có đặc điểm tương đương với da người để giảm bớt các vấn đề như vậy.
Một số loại da nhân tạo đã được sử dụng thành công trong cấy ghép da và
trong điều trị phẫu thuật bỏng [21]. Da nhân tạo có các vùng ưa nước và kỵ
nước và do đó mang đầy đủ các đặc điểm của lớp sừng trên da người.

ho
ol
o

fM

ed

ici

Các mơ hình tương đương với da người (LSEs) cũng được sử dụng để
đánh giá khả năng thẩm thấu qua da. Chúng bao gồm các màng da như: biểu
bì tái tạo được phát triển trong ni cấy mô và được sử dụng làm phương
pháp thay thế cho các mơ động vật. Nói chung, các mơ này yếu hơn và có ít
hơn đáng kể các rào cản như ở da người khi chúng được xây dựng với hàm
lượng Ceramide 1 và 3 thấp hơn bình thường ở da người. Một số nhà nghiên
cứu đã chỉ ra rằng chúng có thể cung cấp các đặc điểm thẩm thấu tương tự
cho da người nhưng trong nhiều trường hợp, chúng vượt quá mức ước lượng
thấm qua da [10, 11, 22, 28].

Co
p

yri

gh
t@

Sc

Khi khơng có da người hoặc da động vật, hoặc khi thực hiện thí nghiệm
với một số lượng lớn, đặc biệt là đối với các thí nghiệm sàng lọc, màng tổng
hợp thường được sử dụng thay thế. Các màng nhân tạo hay được sử dụng nhất
là polydimethylsiloxan (PDMS) và cellulose axetat. Màng PDMS là rào cản
đồng nhất đã được chứng minh giống với các đặc tính rào cản của lớp sừng và
đã sử dụng phổ biến trong các thí nghiệm [5, 36]. Màng cellulose acetate cũng
được sử dụng tương tự trong các thí nghiệm để mơ phỏng tính thẩm thấu của

tế bào. Cấu trúc xốp của nó được so sánh như các đường nước và nang trong
da mà qua đó các phân tử có thể thẩm thấu [7]. Tuy nhiên, các màng này đã
được chứng minh thường làm tăng đáng kể khả năng thẩm thấu của các chất
nên việc sử dụng chúng đã bị giới hạn. Có rất nhiều các cách thức thử nghiệm
6


ne
an
dP

ha
r

ma
c

y,

VN
U

tính rào cản để xác định khả năng thẩm thấu qua da. Để đánh giá mối quan hệ
giữa các phương pháp khác nhau, khả năng thẩm thấu tương đối thông qua
các màng của cùng một hợp chất sẽ được so sánh. Trong thực tế, có ít dữ liệu
so sánh cho mục đích này tồn tại. Để phá vỡ vấn đề này, mối quan hệ định
lượng giữa cấu trúc-tác dụng dựa trên các dữ liệu có thể được mang ra so
sánh. Các mối quan hệ định lượng cấu trúc-tác dụng (QSARs) thống kê khả
năng thẩm thấu qua da của một loạt các hóa chất ngoại sinh với các thơng số
hóa lý đã biết khi được xác định bằng thực nghiệm. Tuy nhiên, có một số

phản đối do những hạn chế của việc phát triển mơ hình thống kê, một số
thống kê khi xây dựng cần được cân nhắc [25].
1.2.3. Mối quan hệ định lượng cấu trúc-tác dụng (QSAR) và khả năng
thẩm thấu qua da

ho
ol
o

fM

ed

ici

Mơ hình QSAR thống kê sự liên hệ giữa hoạt tính sinh học của một loạt
các hợp chất với tính chất hóa lý và hoặc cấu trúc của chúng [20]. Như vậy,
nó cung cấp các mơ hình tiên đốn để ước tính hoạt động sinh học được
nghiên cứu. Trải dài hơn một thế kỷ, phần lớn công việc trong lĩnh vực này đã
tập trung chủ yếu vào việc tối ưu hóa hoạt động của thuốc và dự đốn độc
tính. Với mục đích của luận văn này để nghiên cứu dự đoán khả năng thẩm
thấu qua da của các chất [25].

Sc

Khả năng thẩm thấu của một hợp chất qua da được xác định bằng nhiều
phương pháp. Trong đó, đầu tiên có định luật Fick bao gồm lưu lượng và hệ
số thấm như sau [25]:

gh

t@

J: Lưu lượng thẩm thấu
Dm: Hệ số khuyếch tán

yri

Kp: Hệ số thẩm thấu

Co
p

∆Cm: Chênh lệch nồng độ giữa 2 bên của màng
t: Độ dày của màng

7


VN
U

1.2.4. Quá trình phát triển của QSAR trong nghiên cứu tính thấm qua da.
 Trước 1990:

ed

ici

ne
an

dP

ha
r

ma
c

y,

Những nghiên cứu liên quan đến tính thấm của các hợp chất qua da đối
với tính chất hóa lý của chúng đã được thực hiện trong ít nhất ba mươi năm.
Tuy nhiên, đặc biệt trong các nghiên cứu trước đây, hiếm khi các mơ hình dự
báo định lượng (tức là QSPR) được thực hiện. Phần lớn các nghiên cứu dựa
trên việc phân tích các chuỗi phân tử tương đồng hoặc có liên quan chặt chẽ.
Thơng thường chỉ có một số lượng nhỏ các hợp chất được đánh giá. Nhiều
nghiên cứu trong số này cho thấy mối quan hệ tuyến tính với khả năng kỵ
nước, tức là tăng tính kỵ nước có liên quan đến việc tăng cường độ thẩm thấu
của da [31, 32]. Một số nghiên cứu thì cho thấy mối quan hệ dạng parabol của
độ thẩm thấu qua da với tính kỵ nước, đặc biệt là khi lựa chọn các hợp chất
bao gồm những chất có tính kỵ nước cao [32]. Qua những nghiên cứu này,
cho thấy rằng có thể áp dụng được các mơ hình QSAR cho độ thẩm thấu qua
da, mỗi mơ hình được thiết lập riêng cho một nhóm chất hóa học cụ thể, hoặc
một loạt các hợp chất.

fM

 Tập dữ liệu Flynn (1990) và các phân tích QSAR sau đó:

Co

p

yri

gh
t@

Sc

ho
ol
o

Nỗ lực để phát triển cách tiếp cận mơ hình QSAR tồn diện hơn và phổ
biến hơn về dự đốn khả năng thẩm thấu qua da được tạo điều kiện thuận lợi
nhờ công bố bởi Flynn vào năm 1990 của 97 hệ số thấm cho 95 hợp chất bằng
thí nghiệm in vitro qua da người. Nó cung cấp một cơ sở dữ liệu lớn đầu tiên
và là lớn nhất cho đến gần đây. Tại thời điểm này, cần lưu ý rằng trong khi
các dữ liệu này được đo bằng cách sử dụng da người trong thí nghiệm in
vitro, kho dữ liệu này được biên tập từ 15 nguồn tài liệu khác nhau. Do đó,
khơng thể tránh khỏi rằng những dữ liệu này chứa rủi ro sai lệch giữa các thử
nghiệm là cao, do sự biến đổi giữa các phịng thí nghiệm, đặc biệt là sự biến
đổi phát sinh từ việc sử dụng da từ các nguồn khác nhau và vị trí của phần da
đó trên cơ thể. Về phát triển QSAR cho dự đốn tính thấm của da, cơ sở dữ
liệu của Flynn là một mốc quan trọng [25].
Trong nghiên cứu ban đầu của Flynn (1990) một số thuật toán để dự
đoán khả năng thẩm thấu qua da được đề xuất [6]. Cách tiếp cận được thực

8

Co
p

yri

gh
t@

Sc

ho
ol
o

fM

ed

ici

ne
an
dP

ha
r

ma
c

y,

VN
U

hiện dựa trên nguyên lý tính thấm của da phần lớn do chức năng phân vùng
giữa các lớp nước và khơng chứa nước, và do đó nó có thể được mô tả bởi hệ
số phân chia octanol-nước (Kow) và cũng bị ảnh hưởng bởi kích thước phân
tử. Flynn đã trình bày một cách tiếp cận đơn giản với các thuật tốn cho
những hợp chất có trọng lượng phân tử thấp và cao. Các thuật toán chỉ đơn
giản chỉ ra rằng các hợp chất rất ưa nước và kỵ nước có độ thẩm thấu da thấp
và cao tương ứng, và log Kow khác nhau phụ thuộc vào mơ hình QSAR có thể
được sử dụng để dự đốn tính thấm của da đối với các hợp chất có trọng
lượng phân tử cao và thấp. Khơng có biện pháp thống kê phù hợp nào được
cung cấp tại thời điểm này [6, 25].

9


Các hợp chất có trọng
lượng phân tử cao (>150
đvC)

Log Kp = -3

Log Kp = -5

ma
c

y,

Các hợp chất có trọng
lượng phân tử thấp
(<150 đvC)
Log Kow < 0,5

Log Kp = log Kow - 3,5

0,5 ≤ Log Kow ≤ 3,5

Log Kp = log Kow - 5,5

Log Kp = -0,5

ne
an
dP

Log Kow > 3

ha
r

0,5 ≤ Log Kow ≤ 3,0

VN
U

Bảng 1.1. Thuật tốn tính hệ số thấm (Kp) được công bố bởi Flynn (1990)

Log Kow > 3,5

Log Kp = -1,5

ed

ici

Các dữ liệu được công bố bởi Flynn (1990) đã được phân tích sau đó
bởi rất nhiều nhà nghiên cứu được trình bày trong bảng 1.2 [6, 25, 27]. Potts
and Guy (1992) đã chứng minh việc sử dụng log Kow kết hợp với trọng lượng
phân tử (MW) hoặc thể tích phân tử để dự đốn độ thấm của da và so sánh với
Flynn (1990). Đối với 93 hợp chất, công thức sau được công bố [25, 27]:

fM

Log Kp = 0.71 log Kow – 0.0061 MW – 6.3

Trong đó:

ho
ol
o

(n = 93; r2 = 0.67; s : khơng cơng bố; F : không công bố)

n là số quan sát

Sc

r là hệ số tương quan

s là sai số chuẩn của ước tính

Co
p

yri

gh
t@

F là số liệu thống kê của Fisher

10


R2

Mơ hình

Potts và Guy

93

0,67

Tuyến tính

Pugh và Hadgraft

24

0,935

Tuyến tính

Lien và Gao

22

0,96

Barratt

60

0,90

Potts và Guy

37

0,94

Abraham

46

0,96

Hostynek and Magee

20

Roy và cộng sự

60

Abraham và cộng sự

53

Gute và cộng sự

60

Minghetti và cộng sự

14

Dearden và cộng sự

91

Poulin và Krishnan

47

ma
c

Tuyến tính

ha
r

Tuyến tính
Tuyến tính
Tuyến tính
Tuyến tính

0,64

Tuyến tính

0,96

Tuyến tính

0,67

Tuyến tính

0,749

Tuyến tính

0,83

Tuyến tính

0,95

Mơ hình dựa trên Kp

39

0,94

Tuyến tính

38

0,92

Tuyến tính

37

0,97

Lơgic mờ

54

0,95

Lơgic mờ

Degim và cộng sự

38

0,997

ANN

Moody và MacPherson

65

0,89

Tuyến tính

Majumdar và cộng sự

62

0,934

Tuyến tính

ed

ici

0,80

fM

Tên tác giả

y,

Số chất

ne
an
dP

VN
U

Bảng 1.2. Các nghiên cứu về khả năng thẩm thấu qua da của các hợp chất
được công bố sau 1990

Ghafourian và Fooladi
Gonzalez và Helguera
Pannier và cộng sự

Sc

ho
ol
o

Pannier và cộng sự

gh
t@

1.3. Tổng quan về mối quan hệ định lƣợng cấu trúc-tác dụng của hợp
chất (QSAR)
1.3.1. Khái niệm và nguyên lý chung của QSAR

Co
p

yri

QSAR là một hàm số biểu thị mối quan hệ định lượng giữa cấu trúc
phân tử và hoạt tính sinh học hoặc hóa học của hợp chất có dạng tổng quát
như sau:

11


VN
U

Y = f1(x1) + f2(x2) + ... + fn(xn)

ne
an
dP

ha
r

ma
c

y,

Trong đó, Y là biến đáp ứng sinh học hoặc hóa học thường thu được từ
thực nghiệm và thể hiện bằng giá trị log Kp với Kp là hệ số thẩm thấu qua da.
Các biến x1, x2,…xn là các đại lượng được sử dụng để mô tả cấu trúc phân tử
hợp chất gọi là các tham số phân tử (molecular descriptors). Đến nay đã có
hàng ngàn tham số phân tử (TSPT) được sử dụng trong nghiên cứu xây dựng
mơ hình QSAR. Các TSPT này thu được bằng thực nghiệm hoặc tính tốn
bằng lý thuyết bởi các phần mềm chuyên dụng như Dapt, Dragon, Padell,
MOE, Hyperchem.
Các mơ hình QSAR được xây dựng trên nguyên lý cơ bản là cấu trúc
của mỗi phân tử chứa các yếu tố quyết định đặc tính sinh học, vật lý, hóa lý
của nó; do đó, các hợp chất có cấu trúc giống nhau thường có hoạt tính tương
tự nhau [26].

ici

Từ năm 1868, Crum-Brown và Fraser đã đưa ra nhận định tác dụng
sinh học là hàm số của cấu trúc hóa học:

ed

  f (C)

ho
ol
o

fM

Đến năm 1893, Richet đã rằng sự khác nhau về tác dụng sinh học là do
sự thay đổi về tính chất hóa học hay lý hóa [33]:
  f (C)

Đây là hai tư tưởng nền tảng để thiết lập các mơ hình QSAR sau này.

Sc

Năm 1900, Meyer và Overton tìm ra mối quan hệ tuyến tính giữa mức
độ thân dầu và hoạt tính của thuốc gây mê [29].

gh
t@

Và đến năm 1935, phương trình Hammett ra đời, được coi là phương
trình đầu tiên biểu diễn mối quan hệ định lượng hoạt tính-cấu trúc [29]:
Log 

yri

Với K, Ko là hằng số acid,  là hằng số Hammett (thơng số hóa lý đặc

Co
p

trưng cho khả năng hút hoặc đấy điện tử của nhóm thế).

12


ma
c

Log (1/C) = -k12 + k2 + k3 +k4Es + k5

y,

VN
U

Từ những năm 60 của thế kỉ XX, mơ hình QSAR thực sự được nghiên
cứu bởi Corwin Herman Hansch. Mô hình QSAR Hansch là phương trình hồi
quy nhiều tham số biểu thị mối quan hệ giữa hoạt tính sinh học và các thơng
số hóa lý:

Trong đó: C biểu thị đáp ứng sinh học;  là hằng số thế kỵ nước;  là

ha
r

hằng số thế Hammett; Es là hằng số thế Taft và k1, k2, k3, k4, k5 là các hệ số
phân tích hồi quy

ne

an
dP

Từ đó đến nay, số lượng nghiên cứu về xây dựng các mơ hình QSAR
được cơng bố tăng đột biến và các phương pháp xây dựng, ứng dụng mơ hình
QSAR ngày càng đa dạng.
1.3.2. Xây dựng và đánh giá mơ hình QSAR

Mơ hình QSAR được xây dựng và đánh giá gồm 5 bước như sau:

ici

Bước 1: Chuẩn bị cơ sở dữ liệu

ho
ol
o

fM

ed

Cơ sở dữ liệu (CSDL) để xây dựng các mơ hình QSAR thường là hoạt
tính, độc tính, đặc tính dược động học của các hợp chất được công bố trong
các nghiên cứu khoa học tạo thành dữ liệu tồn cầu (global database). CSDL
cho một mơ hình QSAR cụ thể sẽ bao gồm cấu trúc phân tử của dãy chất và
hoạt tính của chúng, thường gặp nhất là hoạt tính thu được bằng thực nghiệm
in vitro [15].

gh

t@

Sc

CSDL ban đầu được chia thành tập huấn luyện (training set) và tập
kiểm nghiệm (test set). Tập huấn luyện chiếm khoảng 80-90% CSDL ban đầu
và được sử dụng để xây dựng mơ hình QSAR [19]. Tập kiểm nghiệm được sử
dụng để đánh giá mơ hình.
Bước 2: Tính tốn tham số phân tử

Co
p

yri

Cấu trúc phân tử hợp chất được vẽ bằng các phần mềm như Chemdraw,
hoặc bằng chức năng tích hợp trong các phần mềm xây dựng mơ hình QSAR
chun dụng. Từ các cấu trúc này, tham số phân tử (TSPT) mô tả thông tin
cấu trúc 1D, 2D, 3D của hợp chất được tính tốn bằng phần mềm máy tính.

13


ma
c

y,

VN
U

Một số phần mềm thường được sử dụng để tính toán các TSPT hiện nay là
Dragon, Padell, MOE, Hyperchem,... Cùng với sự phát triển của QSAR, ngày
càng nhiều nhóm TSPT được tính tốn để mơ tả chi tiết, chính xác hơn cấu
trúc của hợp chất. Phiên bản cập nhật Dragon 7.0 (2017) có thể tính tốn được
đến 5270 TSPT.
Các TSPT có thể được phân loại thành 3 nhóm như sau:

ne
an
dP

ha
r

Tham số phân tử 1D: Những TSPT có thể thu được từ công thức tổng
quát của hợp chất (khối lượng phân tử, số nhóm thế, số nguyên tử C...hoặc là
những TSPT tính tốn từ sự tổng hợp các ngun tử thành phần như chỉ số
khúc xạ, độ phân cực. Trong đó hệ số phân bố n-octal/nước là tham số thường
sử dụng để đánh giá sinh khả dụng, khả năng thấm qua màng tế bào của hợp
chất.

ici

Tham số phân tử 2D: Những TSPT để mơ tả hình dạng, trật tự liên kết
của các nguyên tử trong cấu trúc của hợp chất. Ví dụ số liên kết -, số vòng
thơm, số liên kết H.

ho
ol

o

fM

ed

Tham số phân tử 3D: Những TSPT mô tả cấu trúc không gian hợp chất,
được xác định bằng hệ trục tọa độ nguyên tử thành phần, bao gồm các TSPT
như mơ men lưỡng cực (dipole moment), diện tích bề mặt cực (polar surface
area), thể tích van der Waals...
Bước 3: Xử lí số liệu

gh
t@

Sc

Để xây dựng được một mơ hình QSAR một cách thuận tiện và đảm bảo
tính tin cậy, dữ liệu ban đầu cần được xử lí sơ bộ. Cơng việc này có thể bao
gồm loại bỏ các giá trị nằm ngồi (outlier), chuẩn hóa đơn vị, chuyển đổi hàm
thức (logarit hóa, mũ hóa, nghịch đảo,...). Những cơng việc này giúp loại bỏ
các TSPT gây nhiễu khi xây dựng mô hình.
Bước 4: Xây dựng mơ hình QSAR

Co
p

yri

Ngày nay, nhiều phương pháp xây dựng mơ hình QSAR được phát

triển và ứng dụng phổ biến. Tuy nhiên, khơng có phương pháp nào được coi
là hiệu quả nhất. Một phương pháp có thể phù hợp để xây dựng mơ hình
QSAR này nhưng có thể sẽ không hiệu quả khi xây dựng các mô hình QSAR

14


VN
U

khác vì mỗi CSDL, mỗi mối quan hệ cấu trúc-hoạt tính sẽ có sự phụ thuộc
khác nhau giữa biến TSPT và biến đáp ứng [15].

ma
c

y,

Có hai dạng mơ hình QSAR chính: một là mơ hình QSAR định lượng,
biến đầu ra là biến liên tục; hai là mơ hình QSAR định tính, biến đầu ra là giá
trị nhị phân có hoặc khơng, âm tính hoặc dương tính (cịn gọi là mơ hình
QSAR phân loại).

ne
an
dP

ha
r

Đối với dạng thứ nhất, các phương pháp xây dựng mơ hình thường áp
dụng là hồi quy đa biến tuyến tính - Multiple Linear Regression (MLR), bình
phương tối thiểu từng phần - Partial Least Squares (PLS), phân tích thành
phần chính - Principal Component Analysis (PCA). Các phương pháp thường
dùng để xây dựng mơ hình QSAR phân loại là mạng nơ-ron nhân tạo Artificial Neural Networks (ANN), phân tích cụm (Cluster Analysis), thuật
toán k láng giềng gần nhất (k-nearest neighbors algorithm).

fM

ed

ici

Các phương pháp này có thể có nguyên tắc lựa chọn biến giống hoặc
khác nhau tuy nhiên đều nhằm mục đích chung là TSPT nào có ảnh hưởng
nhiều nhất đến hoạt tính của phân tử sẽ được ưu tiên đưa vào mơ hình trước;
nhờ đó, giá trị biến đáp ứng dự đoán được sẽ gần với giá trị thực nghiệm nhất.
Bước 5: Đánh giá mơ hình QSAR

ho
ol
o

Đánh giá mơ hình là giai đoạn quan trọng liên quan đến khả năng ứng
dụng của mơ hình. Đánh giá mơ hình thường sử dụng 5 nguyên tắc của
OECD [23].

Sc

Nguyên tắc 1: Có đích xác định (defined end-points).

Co
p

yri

gh
t@

Đích là những giá trị thực nghiệm về hoạt tính sinh học thu được từ
CSDL. Nếu toàn bộ CSDL đều sử dụng một protocol để xác định giá trị của
đích, thì mơ hình xây dựng được có độ tin cậy cao, có thể dùng để dự đốn
hoạt tính các hợp chất mới (nếu thực nghiệm dùng theo quy trình này). Cịn
nếu CSDL sử dụng nhiều protocol khác nhau thì mơ hình sẽ khơng có độ tin
cậy cao.

15


VN
U

Ngun tắc 2: Các thuật tốn sử dụng được mơ tả rõ ràng.

Các thuật toán sử dụng để xây dựng mơ hình cần được mơ tả rõ ràng,
có khả năng ứng dụng để xây dựng các mơ hình khác.

y,

Ngun tắc 3: Có miền cấu trúc ứng dụng xác định.

ne
an
dP

ha
r

ma
c

Miền cấu trúc ứng dụng là khoảng không gian cấu trúc được xác định
bởi các hợp chất trong tập huấn luyện để xây dựng mơ hình. Những hợp chất
thuộc tập huấn luyện nếu có cấu trúc nằm ngồi miền cấu trúc ứng dụng sẽ
ảnh hưởng đến độ chính xác của mơ hình. Những hợp chất thuộc tập dự đốn,
nếu có cấu trúc nằm ngồi miền này sẽ cho kết quả dự đốn khơng chính xác.
Các phương pháp xác định miền ứng dụng thường gặp như: phương pháp đòn
bẩy, phương pháp ba lân cận gần nhất,...
Nguyên tắc 4: Có độ khớp, độ ổn định, và khả năng dự đoán tốt.
Độ khớp hay độ tuyến tính được đánh giá thơng qua giá hệ số xác định
R , R càng cao thì mức độ khớp (tuyến tính) của mơ hình với các giá trị thực
nghiệm càng tốt. R2> 0,6 thì mơ hình mới có ý nghĩa. Cơng thức tính hệ số
xác định như sau:
2

fM

ed

ici

2

ho
ol
o

R2 = 1-

^

∑

y i  yi

∑

( yi  y)
^

Trong đó n là số các quan sát của tập huấn luyện; yi, y i , lần lượt là giá
trị thực tế, giá trị dự đoán của biến phụ thuộc ở quan sát thứ i; y là giá trị thực

Sc

tế trung bình của biến phụ thuộc.

Co
p

yri

gh
t@

Độ ổn định được đánh giá thông qua hệ số tương quan chéo Q2, Q2
được đánh giá dựa trên phương pháp tham số chéo (cross validation leave one
out) bằng cách lần lượt loại một quan sát ra khỏi tập huấn luyện, giữ nguyên
các biến đã lựa chọn và tính tốn lại các thơng số của mơ hình. Q2 càng gần 1,
tính tổng qt hóa của mơ hình càng cao, mơ hình càng ổn định. Thơng
thường, u cầu Q2 > 0,5 thì mơ hình mới bền vững. Cơng thức hệ số tương
quan chéo tính như sau:

16


y i  yi / i

∑

VN
U

Q2 = 1 -

^

∑

( yi  y)

y,

Trong đó n là số các quan sát tập huấn luyện; yi, là giá trị thực tế của
^

ma
c

quan sát thứ i; yi / i giá trị dự đoán của biến phụ thuộc ở quan sát thứ i sử dụng
mô hình đã loại biến i; y là giá trị thực tế trung bình của biến phụ thuộc.

2
Qext

ne
an
dP

ha
r

Khả năng dự đốn ngoại được đánh giá thông qua giá trị Q2ext, Q2ext
càng lớn mơ hình càng cho thấy độ tuyến tính của tập kiểm tra và do vậy
chứng tỏ khả năng dự đốn của mơ hình cao.

=1-

^

∑

y i  yi

∑

( yi  y)

^

Trong đó n là số các quan sát tập kiểm tra; yi, y , lần lượt là giá trị thực

ed

trung bình của biến phụ thuộc.

ici

tế, giá trị dự đốn của biến phụ thuộc ở quan sát thứ i; y là giá trị thực tế

fM

Nguyên tắc 5: Giải thích được cơ chế (nếu có thể).

Co
p

yri

gh

t@

Sc

ho
ol
o

Mơ hình cần được giải thích về vai trị của các biến trong mơ hình, qua
đó giúp định hướng thiết kế các hợp chất mới. Tuy nhiên việc giải thích các
biến này nhiều khi gặp phải khó khăn vì cơ chế của chúng phức tạp địi hỏi
những hiểu biết chuyên sâu về hóa lượng tử. Bên cạnh đó có nhiều biến chưa
có các tài liệu tham khảo, do đó chưa đủ cơ sở để giải thích.

17