Расчёт энергии сольватации димера из двух мономеров полиэтиленгликоля методом FEP
Нам дан 1,2-диметоксиэтан.Его SMILES:COCCOC. Мы разбиваем его на 2 мономера, их SMILES: СОС.
Разобьем на 2 мономера. Это диметиловый эфир. Его SMILES:COC.
Расчитаем частичные заряды на атомах этого остатка используя как заглушку метильную группу. Новый SMILES:COCC.
При подсчёте зарядов полезно использовать специальные возможности RED-vIII.4.pl.
Надо добавить в p2n файл mol_red1.p2n вот такую директиву об обнулении зарядов на метильной группе-заглушке и удалении этой группы из финального файла:
REMARK INTRA-MCC 0.0 | 9 10 11 12 | R
Построение записи для остатка:
Изменим файл aminoacids.rtp, добавив в конец запись нашего мономера:
[ EGL ]
[ atoms ]
CA CT 0.02750 1
HA1 H1 0.05010 2
HA2 H1 0.05010 3
HA3 H1 0.05010 4
O O -0.36850 5
CB CT 0.13770 6
HB1 H1 0.02000 7
HB2 H1 0.02000 8
[ bonds ]
CA HA1
CA HA2
CA HA3
CA O
O CB
CB HB1
CB HB2
CA +CB
Заряды были взяты отсюда: Mol_m1_o1-sm.mol2.
Теперь сделаем файл формата pdb для 1,2-диметоксиэтана. Получился такой файл.

© Свистунова Даша