Date: Thu Aug 19 2021 - 06:34:02 CDT

Can someone please explain to me this script for developing a mutated model
in VMD.

"package require psfgen

proc mutate_res {pos mut prot} {
# in pdb file $prot.pdb, mutate the residue at position $pos to a
# residue of type $mut
#set pos 2
#set mut TRP
#set prot tnc

set file_in $prot.pdb
set file_tmp $prot.temp.pdb
set file_out_pdb $prot.mut.pdb
set file_out_psf $prot.mut.psf

set str "all and not (resid $pos and sidechain)"
mol new tnc.pdb type pdb
set sel [atomselect top $str]
$sel writepdb $file_tmp

package require psfgen
topology top_all27_prot_lipid.inp
pdbalias residue HIS HSE
pdbalias atom ILE CD1 CD
segment P {
  pdb $file_tmp
  mutate $pos $mut}
  coordpdb $file_tmp P
  writepdb $file_out_pdb
  writepsf $file_out_psf
  mol delete top
  mol new $file_out_pdb type pdb

The information contained in this electronic communication is intended 
solely for the individual(s) or entity to which it is addressed. It may 
contain proprietary, confidential and/or legally privileged information. 
Any review, retransmission, dissemination, printing, copying or other use 
of, or taking any action in reliance on the contents of this information by 
person(s) or entities other than the intended recipient is strictly 
prohibited and may be unlawful. If you have received this communication in 
error, please notify us by responding to this email or telephone and 
immediately and permanently delete all copies of this message and any 
attachments from your system(s). The contents of this message do not 
necessarily represent the views or policies of BITS Pilani.