Я выполняю моделирование белка на AQP2, анализируя влияние аминокислотных мутаций на проницаемость воды. Чтобы проверить событие проникновения, у меня хорошо работает этот сценарий проникновения на языке tcl (permeation.tcl).
Но я хочу знать результаты проникновения, происходящие в каждой поре AQP2 отдельно. Поскольку я выполнил одиночное моделирование, но с разными мутациями в каждом мономере.
PS: моделирование выполняется на системе, содержащей водяной бокс и липидный бислой, в который встроен AQP2 без градиента концентрации.
Вот мой сценарий ; Кто-нибудь может помочь изменить этот скрипт в соответствии с требованием, которое я упомянул ранее?