You are here

Unsolved

The problem hasn't been solved

cluster error

Hi all,

I tried clustering decoy pdbs from a docking run using the following command

--------------
/opt/rosetta/rosetta_source/bin/cluster.linuxgccrelease -database /opt/rosetta/rosetta_database -in:file:l ./list-of-pdbs.txt -nooutput > cluster.log3
-------------------

It seems to read the files. Below is the very end of the log file

.................
.
.
.
protocols.jobdist.main: Starting complexppk_moveddecoy_56890_0001 ...
core.conformation.Conformation: Found disulfide between residues 23 93
core.conformation.Conformation: Found disulfide between residues 135 209

Post Situation: 

Failure appending a residue

On my system, the following produces memory errors 3/5 of the time at the third or fourth insertion step, segmentation faults 1/5, and finishes without problems 1/5. Deleting and inserting one residue always works, but the scripts become increasingly unreliable with increasing numbers of residues deleted and inserted. If this problem is not due to a correctable mistake on my part, could someone please suggest how I might dock two domains of a protein without interference from their interdomain linker, and then restore the linker (presumably by building a loop).

Post Situation: 

positive results in total_score by symmetry docking

Hi

I got positive results for total_score after symmetry docking. At the first time, I run the symmetry docking without native structure. the rmsd was not created in the score.fasc. I guessed maybe because I did not specify the native structure. So I use one of the reliable structure obtained from protein-protein docking as a native structure. After doing this, the rmsd is appeared in the score.fasc file. But I got the positive results for total_score. I do not know what the positive number means.

Post Situation: 

Jump question

It appears that one should be able to assign the rotation and translation values of jump 1 in pose1 to jump 1 in pose2 (even if poses 1 and 2 possess different numbers of residues) with the following:
pose2.jump(1).set_rotation(pose1.jump(1).get_rotation())
pose2.jump(1).set_translation(pose1.jump(1).get_translation())

Although I find that these commands suitably alter the translation and rotation values obtained with the commands
pose2.jump(1).get_rotation()
pose2.jump(1).get_translation()

Post Situation: 

Why do non proten ligands get a very high score in kinematic loop modeling/ in general?

hi all,

I am trying to model a loop (using kinematic loop modeling) that is in close proximity to a ligand. So I included the ligand by converting the HETATMs into Rosetta atom types and including centroid (for remodel) and all-atom (for refine) parameter files via the -extra_res_cen and -extra_res_fa command lines. I used the script molfile_to_params.py to create the parameter files (both centroid and all-atom) and the ligand.pdb file.
[I made them from a .mol file using the .pdb file of the ligand using openbabel.... maybe there is a problem here..]

Post Situation: 

score to use for picking top docking solutions?

Hi,

I have done protein-protein docking using rosetta 3.2 to generate 60,000 structures. I'm a bit confused about the next steps.

1) Which score, the total score or I_sc, to use for picking solutions for further consideration? Which of these are referred to in publications?

2) one of the older webpages I found describes three methods for picking solutions (http://graylab.jhu.edu/~mdaily/tutorial/blind_main.html). Are these still followed as a general procedure by the community?

Post Situation: 

How to define specific residues in alascan ?

I have been trying to compute ala-scanning values using the dock design parser applications,
and I also want to specify some residues to be designed (like -resfile option in fixbb design).
Can somebody help me to define specific residues in alanine scanning?
and How to generate output file (.report file)?

Thanks for any help with this problem!

Post Situation: 

Dunbrack rotamer energy term...

This question is regarding the implementation of the Dunbrack rotamer energy term (fa_dun).

Based on my understanding of the implementation in Rosetta 2.3, I noticed that the fa_dun energy is comprised of two components:
1) A backbone-dependent probability of the rotamer : E = -log (p)
2) An "un-normalized" gaussian penalty parametrized using the average chi angles and the standard deviations for the rotamer: E_penalty = -log (gaussian)

Post Situation: 

FlexPepDocking with GDP-bound protein

I have a problem when doing the FlexPepDocking to a receptor with GDP bound, but it stop with an error, saying that GDP is an unrecognized residue. When I add a flag "-remember_unrecognized_res", it goes well but GDP was still missing in the output structure. I believe GDP is involved in the binding and I don't want to ingore it in the docking. And I check the database, GDP.params is in /rosetta_database/chemical/residue_type_sets/fa_standard/residue_types/nucleic.

Post Situation: 

Pages