AsymFoldandDockRbTrialMover.cc

Go to the documentation of this file.

// -*- mode:c++;tab-width:2;indent-tabs-mode:t;show-trailing-whitespace:t;rm-trailing-spaces:t -*-
// vi: set ts=2 noet:
//
// This file is part of the Rosetta software suite and is made available under license.
// The Rosetta software is developed by the contributing members of the Rosetta Commons consortium.
// (C) 199x-2009 Rosetta Commons participating institutions and developers.
// For more information, see http://www.rosettacommons.org/.

/// @file MinMover.cc
/// @brief
/// @author Ingemar Andre

// Unit headers
#include <protocols/simple_moves/asym_fold_and_dock/AsymFoldandDockRbTrialMover.hh>
#include <protocols/rigid/RigidBodyMover.hh>
// AUTO-REMOVED #include <core/pose/Pose.hh>
#include <core/scoring/ScoreFunction.hh>
#include <core/scoring/ScoreFunctionFactory.hh>
#include <core/kinematics/MoveMap.hh>
#include <protocols/moves/MonteCarlo.hh>

// Package headers
// options
#include <basic/options/option.hh>
#include <basic/options/keys/fold_and_dock.OptionKeys.gen.hh>

// ObjexxFCL Headers

// C++ Headers

// Utility Headers
#include <basic/Tracer.hh>

#include <utility/vector1.hh>


namespace protocols {
namespace simple_moves {
namespace asym_fold_and_dock {

static basic::Tracer TR("protocols.simple_moves.symmetry.AsymFoldandDockRbTrialMover");

AsymFoldandDockRbTrialMover::AsymFoldandDockRbTrialMover() :
  protocols::moves::Mover( "AsymFoldandDockRbTrialMover" ),
  smooth_move_(false), rot_mag_(8.0), trans_mag_(3.0), rigid_body_cycles_(50), mc_filter_(true)
{
  scorefxn_ = core::scoring::ScoreFunctionFactory::create_score_function( "interchain_cen" );
}

AsymFoldandDockRbTrialMover::AsymFoldandDockRbTrialMover(
  core::scoring::ScoreFunctionCOP scorefxn
) :
  protocols::moves::Mover( "AsymFoldandDockRbTrialMover" ),
  scorefxn_(scorefxn), smooth_move_(false), rot_mag_(8.0), trans_mag_(3.0), rigid_body_cycles_(50), mc_filter_(true)
{}

AsymFoldandDockRbTrialMover::AsymFoldandDockRbTrialMover(
  core::scoring::ScoreFunctionCOP scorefxn,
  bool smooth_move
) :
  protocols::moves::Mover( "AsymFoldandDockRbTrialMover" ),
  scorefxn_(scorefxn), smooth_move_(smooth_move), rot_mag_(8.0), trans_mag_(3.0), rigid_body_cycles_(50), mc_filter_(true)
{}

AsymFoldandDockRbTrialMover::AsymFoldandDockRbTrialMover(
  core::scoring::ScoreFunctionCOP scorefxn,
  bool smooth_move,
  core::Real rot_mag,
  core::Real trans_mag
) :
  protocols::moves::Mover( "AsymFoldandDockRbTrialMover" ),
  scorefxn_( scorefxn), smooth_move_(smooth_move), rot_mag_(rot_mag), trans_mag_(trans_mag)
{}

void
AsymFoldandDockRbTrialMover::apply( core::pose::Pose & pose )
{

  using namespace basic::options;
  using namespace protocols::moves;

  TR.Debug << "Rb move applied..." << std::endl;

  core::Real trans_mag_smooth = 0.1;
  core::Real rot_mag_smooth = 1.0;

  // Docking options
  if ( smooth_move_ ) {
    if ( option[ OptionKeys::fold_and_dock::trans_mag_smooth ].user() ) {
      trans_mag_smooth = option[ OptionKeys::fold_and_dock::trans_mag_smooth ];
    }
    if ( option[ OptionKeys::fold_and_dock::rot_mag_smooth ].user() ) {
      rot_mag_smooth = option[ OptionKeys::fold_and_dock::rot_mag_smooth ];
    }
  }

  // overrriding constructor or default values
  if ( option[ OptionKeys::fold_and_dock::rb_rot_magnitude ].user() )
  {
    rot_mag_ = ( option[ OptionKeys::fold_and_dock::rb_rot_magnitude ] );
  }

  // overrriding constructor or default values
  if ( option[ OptionKeys::fold_and_dock::rb_trans_magnitude ].user() )
  {
    trans_mag_ = ( option[ OptionKeys::fold_and_dock::rb_trans_magnitude ] );
  }

  core::Real rot_mag_trial = smooth_move_ ? rot_mag_smooth : rot_mag_;
  core::Real trans_mag_trial = smooth_move_ ? trans_mag_smooth : trans_mag_;

  // overrriding constructor or default values
  if ( option[ OptionKeys::fold_and_dock::rigid_body_cycles ].user() )
  {
    rigid_body_cycles_ = ( option[ OptionKeys::fold_and_dock::rigid_body_cycles ] );
  }

  if ( option[ OptionKeys::fold_and_dock::rigid_body_disable_mc ].user() )
  {
    mc_filter_ = false;
  }

  // Setup Monte Carlo object
  protocols::moves::MonteCarloOP monteCarlo_ = new protocols::moves::MonteCarlo(pose, *scorefxn_, 2.0 );


  // Setup the movemap
  core::kinematics::MoveMapOP movemap = new core::kinematics::MoveMap();
  movemap->set_bb( false );
  movemap->set_chi( false );
  movemap->set_jump( true );

  //set up mover for docking
  protocols::rigid::RigidBodyPerturbMover rb_perturb =
    protocols::rigid::RigidBodyPerturbMover(
      pose, *movemap, rot_mag_trial, trans_mag_trial,
      rigid::partner_downstream, false );

  for ( Size i = 1; i <= rigid_body_cycles_; ++i ) {
    rb_perturb.apply( pose );
    if ( mc_filter_ ) monteCarlo_->boltzmann( pose );
  }
  //monteCarlo_->recover_low(pose);
}

std::string
AsymFoldandDockRbTrialMover::get_name() const {
  return "AsymFoldandDockRbTrialMover";
}


} // asym_fold_and_dock
} // moves
} // protocols