Trouble enforcing a disulfide in loop building

Member Site Forums Rosetta 3 Rosetta 3 – Applications Trouble enforcing a disulfide in loop building

Viewing 1 reply thread
  • Author
    Posts
    • #2205
      Anonymous

        Hello,

        Like some others on this forum, I’m having some trouble getting disulfides to work. I pass

        -in:fix_disulf ./disulf

        where disulf is a text file containing the residue numbers of the cysteines I am trying to bridge like so:

        269 280

        With an empty line at the end, which others have mentioned may be necessary. Unfortunately, this has no discernable effect on the loop in question–I have done pretty extensive testing, and I never see a disulfide form. Should I try to encode this disulfide as a constraint, or is there something I’m clearly doing wrong?

        EDIT:

        Given that -in:fix_disulf has no effect by my lights, can any of you suggest a benchmark or standard to test whether this command is working properly in my installation? I’ve had success with other projects, so it would have to be a pretty specific issue.

        Thanks,
        CH

      • #11019
        Anonymous

          -in:fix_disulf (note the underscore and not a colon between fix and disulf) is handled at the application level, rather than at a general Rosetta level – depending on your protocol it may not even be read. Which protocol are you using?

          If your protocol does not obey -in:fix_disulf, then (atom pair) constraints may be the way to go – many protocols, even if they don’t obey -in:fix_disulf, will automatically detect disulfides if the two cysteines are close enough.

          For an example of -in:fix_disulf usage, you can take a look at Rosetta/demos/protocol_capture/2012/replica_docking/, which makes use of it.

        • #11026
          Anonymous

            Looking back at what I’ve done, I was passing -in:fix_disulf with correct punctuation, I just reproduced it incorrectly in my first post.

            I’m using a homology modelling protocol to do loop building on a membrane protein, here’s the contents of the options file:

            -in:file:fullatom
            -s ./hTRPM8_full_c-term_threaded.pdb
            -loops:loop_file ./test.loops ##loop file
            -loops:extended true
            -loops:build_initial
            -loops:idealize_after_loop_close
            -loops:relax fastrelax
            -loops:fast
            -out:file:fullatom
            -use_input_sc
            -ex1 #Include extra chi1 rotamers
            -ex2
            -loops:frag_sizes 9 3 1
            -loops:frag_files ./aat000_09_05.200_v1_3 ./aat000_03_05.200_v1_3 none
            -loops:remodel quick_ccd
            -loops:refine refine_kic
            -in:file:psipred_ss2 ./t000_.psipred_ss2
            -in:file:spanfile ./modellerrebuild.span
            -in:file:lipofile ./hTRP.lips4
            -in:fix_disulf ./disulf
            -score:weights /usr/local/apps/rosetta/rosetta-3.5/rosetta_database/scoring/weights/membrane_highres_Menv_smooth.wts
            -overwrite

            Any idea what the problem might be? I’ve begun to try using an atom pair constraint, we’ll see how that goes.

            EDIT: For anyone with similar problems, I seem to have successfully generated disulfides by passing this constraint file:

            AtomPair SG 269 SG 280 HARMONIC 1.800 2.05

            with these options:

            -cst_fa_file disulf.cst
            -cst_fa_weight 1000

        Viewing 1 reply thread
        • You must be logged in to reply to this topic.