Member Site Forums Rosetta 3 Rosetta 3 – Applications Dimer of dimers

  • This topic has 4 replies, 2 voices, and was last updated 13 years ago by Anonymous.
Viewing 4 reply threads
  • Author
    • #873


        I want to use Fold and Dock on a small domain of a protein that forms a dimer-of-dimers. It is a 50 residue domain that is known to form helixes.

        I am not able to get a reasonable quaternary structure for this protein using the below symm files. Any help would be deeply appreciated. Thanks!

        D2.symm file:

        set_dof BASEJUMP angle_x(0:360) angle_y(0:360) angle_z(0:360)
        set_dof JUMP2 angle_z(0:90)

      • #5437

          Your symmetry definition file does not appear to allow for translational moves. Typically, the set_dof statements would look like this:

          set_dof BASEJUMP x(50) angle_x(0:360) angle_y(0:360) angle_z(0:360)
          set_dof JUMP2 z(50) angle_z(0:90.0)

          The x(50) and z(50) items allow for translational optimization and specify that the coordinates of subunits are placed 50 A from the origin at the beginning of each simulation.

          Then nstruct 9999 is probably on the low side given the number of degrees of freedom in the system.

        • #5448

            If I add the x(50) and z(50), the monomers don’t come together. In fact, x(1) and z(1) also does not bring the monomers together. But I will run more nstruct to see if it changes my results.

          • #5490

              No changes in the results. The lowest energy structure forms dimers but not dimer-of-dimers. Any other suggestions that I can try?

            • #5501

                Its hard to say without looking at the result in detail. Is there any sampling of closed conformations at all? Does it converge into a particular dimer structure? This behavior would be expected if the interactions between dimers is not sufficiently large to converge to a particular dimer-dimer orientation.

                The x(1) z(1) does not make much sense because the subunits would end up on top of each other at the beginning of the simulation. You could play around with alternative ways of sliding subunits together. In the script you can control the slide behavior ( -slide_type (RANDOM|SEQUENTIAL|ORDERED_SEQUENTIAL).

                If there is a good convergence for a dimer you can try symmetric docking of these dimers followed by a symmetric relax.

                You could play around wit

            Viewing 4 reply threads
            • You must be logged in to reply to this topic.