Hi, I am using secondary restraints in the phenix.refine. I got the following error message. hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 141:143) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 141:143) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 160:162) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 160:162) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 242:244) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 242:244) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 360:363) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 360:363) and (altloc 'A' or altloc ' ') and name O" => 4 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 387:389) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 387:389) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 455:457) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 455:457) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 143:152) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 143:152) and (altloc 'A' or altloc ' ') and name O" => 10 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 156:165) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 156:165) and (altloc 'A' or altloc ' ') and name O" => 10 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 190:193) and (altloc 'A' or altloc ' ') and name H" => 1 donors "(chain 'B' and resseq 190:193) and (altloc 'A' or altloc ' ') and name O" => 4 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 222:227) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 222:227) and (altloc 'A' or altloc ' ') and name O" => 6 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 243:251) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 243:251) and (altloc 'A' or altloc ' ') and name O" => 9 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 253:256) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 253:256) and (altloc 'A' or altloc ' ') and name O" => 4 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 278:282) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 278:282) and (altloc 'A' or altloc ' ') and name O" => 5 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 301:305) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 301:305) and (altloc 'A' or altloc ' ') and name O" => 5 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 311:321) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 311:321) and (altloc 'A' or altloc ' ') and name O" => 11 acceptors Traceback (most recent call last): File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/phenix/phenix/refinement/command_line.py", line 71, in run overwrite = command_line_interpreter.command_line.options.overwrite) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/phenix/phenix/refinement/driver.py", line 121, in __init__ h_bond_table = self.sec_str.get_bonds_table(log=log) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 1023, in get_bonds_table log=log) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 444, in hydrogen_bonds_from_selections ss_type="helix") File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 652, in _donors_and_acceptors donor_isel.insert(k, n_atoms) IndexError: Index out of range. How do I fix it? Thanks. -- Jianghai
On Mon, Jun 21, 2010 at 11:55 AM, Jianghai Zhu
I am using secondary restraints in the phenix.refine. I got the following error message. hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 141:143) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 141:143) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors . . .
Traceback (most recent call last):
. . .
File
"/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 652, in _donors_and_acceptors donor_isel.insert(k, n_atoms) IndexError: Index out of range.
The code will try to switch intelligently between restraining H-O and N-O distances, depending on whether or not hydrogen atoms are present in the PDB file - but it does so globally, and *any* hydrogen atoms are enough to trigger use of the H-O restraints. The warning messages are telling you that hydrogens are missing on those residues; this could happen if you added residues in Coot, but didn't re-run phenix.ready_set to fully hydrogenate them. The IndexError is just a bug, probably due to some special case present in your model. If you intended to refine with hydrogens, you should run phenix.ready_set first; if not, I'd recommend deleting whatever hydrogen atoms are present now, and/or add this parameter: h_bond_restraints.substitute_n_for_h=True which will tell it to only restrain N-O distances. I'd still like to figure out why it crashes, however - could you please send me the PDB file privately? (I don't need data for this, just the model.) thanks, Nat
participants (2)
-
Jianghai Zhu
-
Nathaniel Echols