<PE_scan.pl><pot_energy.pl>Thanks for your help. I got it to run ok. I was a bit taken back by the numbers until I discovered the energy is in KJ/mol. I'm not sure if anyone is interested but, I have attached 2 perl scripts (sorry, I don't know python) that generates the python script for chimera then analyzes the output.Make PE_scan.pl (pot_energy.pl) and executable and type PE_scan.pl (and pot_energy.pl) for instructions.
One requirement is that the atoms (involved in the dihedral angle) must have unique names in the structure file.
~>PE_scan.pl C1 C2 C3 C4 > PEscan.pyopen your molecule in chimera and on the command line type: open PEscan.py
An output file is generated named: reply_log.txt
An angle (degrees) vs. energy (kJ/mole) table is produced upon running pot_energy.pl.
~>pot_energy.pl reply_log.txt >energy
If you have xmgrace installed:
~>xmgrace energy
will plot it.Hopefully I have caught all of the bugs, but one never knows.
On Tue, Aug 26, 2014 at 1:51 PM, Eric Pettersen <pett@cgl.ucsf.edu> wrote:
On Aug 26, 2014, at 11:49 AM, Eric Pettersen <pett@cgl.ucsf.EDU> wrote:dihed = angle(":241,a@CA,CB,CG,ND")Not that it matter too much, but I meant :241.a, not :241,a in the example.--Eric
---- Milo
===================================================
National Magnetic Resonance Facility at Madison
An NIH-Supported Resource Center
W. Milo Westler, Ph.D.
NMRFAM Director
Senior Scientist
and
Adjunct Professor
Department of Biochemistry
University of Wisconsin-Madison
DeLuca Biochemistry Laboratories
433 Babcock Drive
Rm B160D
Madison, WI USA 53706-1544
EMAIL: milo@nmrfam.wisc.edu
PHONE: (608)-263-9599
FAX: (608)-263-1722
======================================================================= ========