remarks changed Tue Sep 2 09:48:07 CDT 2003 remarks by jfd to include a fes residue !@generate_tmoc.inp topology reset end structure reset end topology {===>} @TOPPAR:topallhdg5.3.pro end ! no chain id. segment name=" " chain @TOPPAR:topallhdg5.3.pep sequence @1xxx.seq end end end ! removes the hydrogen on the cys-en ! and changes the atom type from SH1E to S. patch DISU reference=1=( resid 50 ) reference=2=( resid 53 ) end ! turn peptide from trans to cis patch CISP reference=nil=( resid 92 ) end evaluate ( $par_nonbonded = "OPLSX" ) parameter reset end parameter @TOPPAR:parallhdg5.3.pro end !coor @1xxx_xplor.pdb ! show everything is consistent (or not) !energy end ! print thres=0.05 bonds ! print thres=5. angles ! print thres=5. impropers ! print thres=140 dihe !coor init sele ( all ) end write psf output=1xxx.psf end stop ! junk not reached by script. print thres=0.05 bonds print thres=5. angles print thres=0.5 impropers print thres=140 dihe