You are here
I have a question regarding Pyrosetta. Currently, I want to use pyrosetta to score some protein structures but with missing atoms in sidechain. So my protein structures contrain full-atoms for backbone but missing one or two atoms per residue on sidechain. I used pose_from_pdb to read my protein structure (pdb). But it seems pysoretta will fill in the missing atoms according to rotamer library. Is there a way to prevent it so I can just score my structures with missing atoms on sidechains?
Thanks a lot for the help,
A pose can have a constraint set, a modifiable residue set, an unnoticed `pose.pdb_info().obsolete() == True` case, virtual residues, energies and a few more. Therefore when saving to disk none of the following truly work as they save coordinates:
pose.dump_pdb('foo.pdb') pose.dump_scored_pdb('foo.pdb', pyrosetta.get_fa_scorefxn()) pose.dump_cif('foo.cif') pose.dump_mmtf('foo.mmtf') pyrosetta.io.poses_to_silent([pose], 'foo.silent')
So I was hoping for an all-in-one
I have been trying to load 'HYP.params' for the hydroxyproline residue. According to previous answers in the forum there are 2 options: add the params file to the residue_types folder and editing the residue_types.txt or loading the params file via command lines in pyrosetta. Neither of these options work for me and raise the following error:
I'm providing the molfile_to_params.py with the enzyme-ligand complex of my molecule but unfortunetly it gives the "ValueError: can only read V2000 format files" error. Previously i gave it the individual ligand file and it processed with no problems. But when i try to process this file it doesn't seem to work. I suspect the problem is about the chemical table file formats and something in my code is from V3000 but i'm a begginer so i don't know what.
I will attach the text files of both the ligand .mol file and the complex .mol file.