zoukankan      html  css  js  c++  java
  • rosetta constraints

    通过设置constraints,可以实现一些键长键角以及二面角的限制:

    AtomPair P1 289 OD1 19 HARMONIC 2.10 0.02
    
    Angle O9 289  P1 289 OD1  19 CIRCULARHARMONIC 1.605703 0.02
    Angle P1 289 OD1  19  CG  19 CIRCULARHARMONIC 2.268928 0.02
    
    Dihedral O5 289  O9 289  P1 289 OD1 19 CIRCULARHARMONIC 3.141592 0.02
    Dihedral P1 289 OD1  19  CG 19  CB  19 CIRCULARHARMONIC 3.403392 0.02
    Dihedral O9 289  P1 289 OD1  19 CG  19 CIRCULARHARMONIC 6.143559 0.02

    倒数第二列是弧度值,pai代表180度。

    距离设置规则可以参考rosetta官网:

    https://www.rosettacommons.org/docs/latest/rosetta_basics/file_types/constraint-file

    在程序中引入constraints,可以使用mover或者直接在flags里面引入:

    https://www.rosettacommons.org/docs/latest/scripting_documentation/RosettaScripts/Movers/movers_pages/ConstraintSetMover

    附上自己写的做限制性对接的scripts,其中用到了constraints, cst等

    <ROSETTASCRIPTS>This protocol will simply do low-resolution followed by high-resolution docking.
    It will also report the binding energy (ddg) and buried-surface area (sasa) in the score file.
        <SCOREFXNS>
            <ScoreFunction name="ligand_soft_rep" weights="ligand_soft_rep">
            </ScoreFunction>
            <ScoreFunction name="hard_rep" weights="ligand">
            <Reweight scoretype="metalbinding_constraint" weight="1.0" />
            <Reweight scoretype="atom_pair_constraint" weight="0.333" />
            <Reweight scoretype="angle_constraint" weight="0.333" />
            <Reweight scoretype="dihedral_constraint" weight="0.333" />
            </ScoreFunction>
        </SCOREFXNS>
    
        <LIGAND_AREAS>
            <LigandArea name="docking_sidechain_X" chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>
            <LigandArea name="final_sidechain_X" chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>
            <LigandArea name="final_backbone_X" chain="X" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>
        </LIGAND_AREAS>
        <INTERFACE_BUILDERS>
            <InterfaceBuilder name="side_chain_for_docking" ligand_areas="docking_sidechain_X"/>
            <InterfaceBuilder name="side_chain_for_final" ligand_areas="final_sidechain_X"/>
            <InterfaceBuilder name="backbone" ligand_areas="final_backbone_X" extension_window="3"/>
        </INTERFACE_BUILDERS>
        <MOVEMAP_BUILDERS>
            <MoveMapBuilder name="docking" sc_interface="side_chain_for_docking" minimize_water="true"/>
            <MoveMapBuilder name="final" sc_interface="side_chain_for_final" bb_interface="backbone" minimize_water="true"/>
        </MOVEMAP_BUILDERS>
        <SCORINGGRIDS ligand_chain="X" width="8">
            <ClassicGrid grid_name="classic" weight="1.0"/>
        </SCORINGGRIDS>
        <MOVERS>
        single movers_X
            <AddOrRemoveMatchCsts name="cstadd" cst_instruction="add_new"/> add catalytic constraints
            <SetupMetalsMover name="setup_metals" metal_resnums="19,21,232,284,286,289" remove_hydrogens="false" metals_detection_LJ_multiplier="1.0" metals_distance_constraint_multiplier="1.0" metals_angle_constraint_multiplier="1.0" />
            <ConstraintSetMover name="constraints_S6P" add_constraints="true" cst_file="inputs/constraints.S6P" />
            <Transform name="transform_X" chain="X" box_size="8.0" move_distance="0.2" angle="5" cycles="800" repeats="3" temperature="5"/>
            <AddOrRemoveMatchCsts name="cstrem" cst_instruction="remove" keep_covalent="1"/> remove constraints
            <HighResDocker name="high_res_docker" cycles="6" repack_every_Nth="3" scorefxn="ligand_soft_rep" movemap_builder="docking"/>
            <FinalMinimizer name="final" scorefxn="hard_rep" movemap_builder="final"/>
            <AddOrRemoveMatchCsts name="cstfinadd" cst_instruction="add_pregenerated"/>
            <InterfaceScoreCalculator name="add_scores" chains="X" scorefxn="hard_rep"/>
            
        compound movers
            <ParsedProtocol name="low_res_dock">
                <Add mover_name="cstadd"/>
                <Add mover_name="transform_X"/>
                <Add mover_name="cstrem"/>
            </ParsedProtocol>
            <ParsedProtocol name="high_res_dock">
                <Add mover_name="high_res_docker"/>
                <Add mover_name="constraints_S6P"/>
                <Add mover_name="final"/>
                <Add mover_name="cstfinadd"/>
            </ParsedProtocol>
        </MOVERS>
        <PROTOCOLS>
            <Add mover_name="low_res_dock"/>
            <Add mover_name="high_res_dock"/>
            <Add mover_name="add_scores"/>        
        </PROTOCOLS>
    </ROSETTASCRIPTS>
    

      

     
  • 相关阅读:
    [转]HDR渲染器的实现(基于OpenGL)
    32位微机的内存管理模式
    王爽汇编检测点1.1
    标志位的说明
    通用寄存器的作用
    Debug初次使用
    昨天纠结了一天要始学习dos汇编还是win32汇编,终于想通了
    CPU对存储器的读写(二、数据总线、控制总线)
    CPU对存储器的读写(一、地址总线)
    段寄存器的引用
  • 原文地址:https://www.cnblogs.com/wq242424/p/11264028.html
Copyright © 2011-2022 走看看