{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "*This notebook contains material from [PyRosetta](https://RosettaCommons.github.io/PyRosetta.notebooks);\n", "content is available [on Github](https://github.com/RosettaCommons/PyRosetta.notebooks.git).*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "< [RosettaCarbohydrates: Trees, Selectors and Movers](http://nbviewer.jupyter.org/github/RosettaCommons/PyRosetta.notebooks/blob/master/notebooks/13.01-Glycan-Trees-Selectors-and-Movers.ipynb) | [Contents](toc.ipynb) | [Index](index.ipynb) | [RNA in PyRosetta](http://nbviewer.jupyter.org/github/RosettaCommons/PyRosetta.notebooks/blob/master/notebooks/14.00-RNA-Basics.ipynb) >

\"Open" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# RosettaCarbohydrates: Modeling and Design\n", "Keywords: carbohydrate, glycan, glucose, mannose, sugar, design, prediction\n", "\n", "## Overview\n", "Here, you will learn how to model glycans and design optimal glycosylation positions in a protein.\n", "\n", "We will be using the RosettaCarbohydrate framework to build and model glycans. The `GlycanModeler`, which is our main method for modeling glycans, will be published in 2020. We will be using some custom glycan options to load pdbs. \n", "First, one needs the `-include_sugars` option, which will tell Rosetta to load sugars and add the sugar_bb energy term to a default scorefunction. This scoreterm is like rama for the sugar dihedrals which connect each sugar residue. \n", "\n", "\t\t-include_sugars\n", "\n", "\n", "When loading structures from the PDB that include glycans, we use these options. This includes an option to write out the structures in pdb format instead of the Rosetta format (which is actually better). Again, this is included in the config/flags files you will be using.\n", "\n", "\t\t-maintain_links\n", "\t\t-auto_detect_glycan_connections\n", "\t\t-alternate_3_letter_codes pdb_sugar\n", "\t\t-write_glycan_pdb_codes\n", "\n", "\n", "More information on working with glycans can be found at this page: [Working With Glycans](https://www.rosettacommons.org/docs/latest/application_documentation/carbohydrates/WorkingWithGlycans)\n", "\n", "## Algorithm\n", " \n", "The `GlycanModeler` essentially builds glycans from the root (The first residue of the Tree) out to the trees in a way that simulates a tree growing. It uses a notion of a 'layer' where the layer is defined as the number of residues to the glycan root (with the glycan root being layer 0). Within modeling, all glycan residues other than the ones being optimized are 'virtualized'. In Rosetta, the term 'Virtual' means that these residues are present, but not scored. (It should be noted that it is now possible to turn any residues Virtual and back to Real using two movers: `ConvertVirtualToRealMover` and `ConvertRealToVirtualMover`. )\n", "\n", "Within the modeling application, sampling of glycan DOFs is done through the `GlycanSampler`. The sampler attempts to sample the large amount of DOFs available to a glycan tree. The GlycanSampler is a `WeightedRandomSampler`, which is a container of highly specific sampling strategies, where each strategy is weighted by a particular probability. At each apply, the mover selects one of these samplers using the probability set to it. This is the same way the SnugDock algorithm for antibody modeling works. \n", "\n", "Sampling is always scaled with the number of glycan residues that you are modeling, so run-time will increase proportionally as well. \n", "If you are modeling a huge viral particle with lots of glycans, one can use quench mode, which will optimize each glycan individually. \n", "Tpyically for these cases, multiple rounds of glycan modeling is desired. \n", "\n", "\n", "### GlycanSampler Major components\n", "\n", "Some of these components were covered in the previous tutorial.\n", "\n", "1. __Glycan Conformers__\n", "\n", "\tThese conformers have been generated through an in-depth bioinformatic analysis of the PDB using adaptive kernal density estimates and are unique for each linkage type including glycan residues connected to ASN residues. A conformer is a specific conformation of all of the backbone dihedrals of a particular glycan linkage. Essentialy glycan 'fragments' for a particular type of linkage.\n", "\n", "\n", "2. __SugarBB Sampling__ \n", "\n", "\tThis sampling is done through turning the `sugar_bb` energy term into a set of probabilities using the -log(e) function. This allows us to sample on the QM derived torsonal potentials during modeling. \n", "\n", "\n", "3. __Random Sampling and Shear Moves__\n", "\n", "\tWe sample random torsions at +/- 15 , +/- 45, +/- 90 degrees, each at decreasing probabilities at a 4:2:1 ratio of sampling Small,Medium,Large. \n", "\tShear sampling is done where torsions are set for two residues in order to reduce downsteam effects and allow 'flipping' of the glycan torsions.\n", "\n", "\n", "4. __Minimization__\n", "\t\n", "\tWe Minimize Sugar residues by randomly selecting a residue from what is set to model, and selecting all residues out to the tree that are not virtualized. This reduces computational time that would otherwise restrict the total number of glycan residues we could model at once.\n", " \n", "\n", "5. __Packing__\n", "\n", "\tOf the residues set to optimize, we chooses a random residue and pack that residue and all residues out to the tree that are not virtualized. We pack the sugar residues (OH and constituents) and any neighboring protein sidechains. TaskOperations may be set to allow design of protein residues during this. We do packing this way to once again reduce total computational time.\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "!pip install pyrosettacolabsetup\n", "import pyrosettacolabsetup; pyrosettacolabsetup.install_pyrosetta()\n", "import pyrosetta; pyrosetta.init()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Make sure you are in the directory with the pdb files:**\n", "\n", "`cd google_drive/MyDrive/student-notebooks/`" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# General Setup and Inputs\n", "\n", "You will be using a few different inputs. We will be designing in glycosylation spots in order to block antibody binding at a highly curved epitope, and we will be loading a human structure from the PDB that has internal glycans. \n", "\n", "\n", "## Notes for Tutorial Shortening\n", "\n", "\n", "Typically, the value of `-glycan_sampler_rounds` is set to 25 (which typically is enough) and nstruct is about 5-10k per input structure. You may increase glycan_sampler_rounds to 100 and then decrease output to 1-2500 nstruct in order to have the same level of sampling, which will result in very good models as well. Since this is denovo modeling of glycans, more nstruct is almost always better. For some tutorials, we may decrease this value below our optimal value in order to shorten the length of the tutorial.\n", "\n", "\n", "## General Notes\n", "\n", "We will use a flags file for all common options in this tutorial. Note that instead of passing this flag on init, you can instead put it into your working directory or a particular place in your home directory and rename it common. \n", " \n", "See this page for more info on using rosetta with custom config files: \n", "\n", "All tutorials have generated output in output_files and their approximate time to finish on a single (core i7) processor.\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "PyRosetta-4 2019 [Rosetta PyRosetta4.Release.python36.mac 2019.39+release.93456a567a8125cafdf7f8cb44400bc20b570d81 2019-09-26T14:24:44] retrieved from: http://www.pyrosetta.org\n", "(C) Copyright Rosetta Commons Member Institutions. Created in JHU by Sergey Lyskov and PyRosetta Team.\n", "\u001b[0mcore.init: \u001b[0mChecking for fconfig files in pwd and ./rosetta/flags\n", "\u001b[0mcore.init: \u001b[0mReading fconfig.../Users/jadolfbr/.rosetta/flags/common\n", "\u001b[0mcore.init: \u001b[0m\n", "\u001b[0mcore.init: \u001b[0m\n", "\u001b[0mcore.init: \u001b[0mRosetta version: PyRosetta4.Release.python36.mac r233 2019.39+release.93456a567a8 93456a567a8125cafdf7f8cb44400bc20b570d81 http://www.pyrosetta.org 2019-09-26T14:24:44\n", "\u001b[0mcore.init: \u001b[0mcommand: PyRosetta @inputs/glycans/common_glycans @inputs/glycans/pdb_flags @inputs/glycans/map_flags -database /Users/jadolfbr/Library/Python/3.6/lib/python/site-packages/pyrosetta-2019.39+release.93456a567a8-py3.6-macosx-10.6-intel.egg/pyrosetta/database\n", "\u001b[0mbasic.random.init_random_generator: \u001b[0m'RNG device' seed mode, using '/dev/urandom', seed=2011870376 seed_offset=0 real_seed=2011870376\n", "\u001b[0mbasic.random.init_random_generator: \u001b[0mRandomGenerator:init: Normal mode, seed=2011870376 RG_type=mt19937\n" ] } ], "source": [ "#Python\n", "from pyrosetta import *\n", "from pyrosetta.rosetta import *\n", "from pyrosetta.teaching import *\n", "init('@inputs/glycans/common_glycans @inputs/glycans/pdb_flags @inputs/glycans/map_flags')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial\n", "\n", "GlycanModeling is done through the RosettaScripts interface. Each tutorial has you copying a base XML and adding/modifying specific components to achieve a goal. ALL of these movers are available as components in PyRosetta - however, setup is much more difficult and time consuming. So for now, we will rely on the RS interface, but \n", "\n", "## Tutorial A: Epitope Blocking, De-novo Glycan Modeling\n", "\n", "Here, we will start with the antigen known as Bee Hyaluronidase, from PDB ID 2J88. The PDB file has an antibody bound to it as a HIGHLY immunogenic site. We would like to block this in order to begin to use this enzyme for therapy as Hyaluronidase can be effective in breaking down sugars in the extracellular matrix, allowing certain larger drugs to get to regions of interest. The antibody is renumbered into the AHo numbering scheme that we used in the RAbD tutorial, and it has been relaxed with constraints into the Rosetta energy function. \n", "\n", "We will be designing in at least one optimal glycan at the most immunogenic site.\n", "Note that a prototocol called SugarCoat is in development that will scan regions of interest for potential ideal glycosylation, however, one can certainly do this manually as we do below. \n", "\n", "### A 1. Designing in a Glycosylation Site: \n", "\n", "CreateGlycanSequonMover` and `CreateSequenceMotifMover`\n", "\n", "A sugar glycosylation site is known as a `Sequon`. The glycan sequon is made up of three protein residues which are recognized by the GlycosylTransferase Enzyme during translation in the ER. This enzyme adds the root of nascent glycan onto a protein. In this case, we use the sequon for ASN glycosylation. The sequon is as follows: `N[^P][S/T]`. The `[^P]` notation means that any residue other than P can be there. The `[S/T]` notation means that either S or T is recognized. This notation can be used to directly create Motifs in proteins using the `CreateSequenceMotifMover` and associated `SequenceMotifTaskOperation`. Documentation for these is available here:\n", "\n", "\n", "\t\n", "\n", "The create GlycanSequonMover can also be used for glycosylation of different AA than ASN.\n", "\n", "#### A1.1 Design using a typical sequon\n", "\n", "Before we begin, take a look at the complex either by PyMol or use the PyMolMover as you have previously. \n", "The file is `inputs/glycans/2j88_complex.pdb`\n", " \n", "Where can we introduce a glycan to block binding?\n", " \n", " \n", "Where do you think the optimal glycan position would be for this particular antibody? Take a look at the xml. Is this the position we are targeting? Typically, we may want to allow some backbone movement in our sequon. The full glycan scanning protocol can be found in an input file, simple_glycan_scanner_manual.xml, where we relax the motif residues with constraints, add the sequon, and then relax again, comparing the energy between them to get the full energetic contributions of the sequon on the structure. In order to reduce the run time in these tutorials, we will be removing this going forward.\n", "\n", "The XML syntax is below: \n", " \n", " \n", " \n", "Go ahead and run the xml (`inputs/glycans/tutA11.xml`) using what you have learned previously, or run the code below (about 15 seconds)\n", "\n", "Note that the xml uses `SimpleMetrics` to output a variety of metrics that are in pose.scores or output into a scorefile.\n", "\n", "Here is the full XML. We will only use part of it in code." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```xml\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\t\t\n", "\t\t\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n", "\n", "```" ] }, { "cell_type": "code", "execution_count": 39, "metadata": {}, "outputs": [], "source": [ "from rosetta.protocols.carbohydrates import *\n", "from rosetta.core.select.residue_selector import *\n", "from rosetta.core.simple_metrics.metrics import *\n", "from rosetta.core.simple_metrics.composite_metrics import *\n", "from rosetta.core.simple_metrics.per_residue_metrics import *" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.import_pose.import_pose: \u001b[0mFile 'inputs/glycans/2j88_complex.pdb' automatically determined to be of type PDB\n", "\u001b[0mcore.io.pdb.pdb_reader: \u001b[0mParsing 660 .pdb records with unknown format to search for Rosetta-specific comments.\n", "\u001b[0mcore.io.util: \u001b[0mAutomatic glycan connection is activated.\n", "\u001b[0mcore.io.util: \u001b[0mStart reordering residues.\n", "\u001b[0mcore.io.util: \u001b[0mCorrected glycan residue order (internal numbering): []\n", "\u001b[0mcore.io.util: \u001b[0m\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 56\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 57\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 359\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 360\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 434\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 435\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 452\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 453\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 461\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 462\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 478\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 479\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 616\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 617\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 633\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 634\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 13 299\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 13 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 299 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 13 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 299 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 175 187\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 175 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 187 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 175 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 187 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 342 406\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 342 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 406 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 342 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 406 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 437 480\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 437 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 480 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 437 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 480 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 512 587\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 512 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 587 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 512 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 587 CYD\n", "\u001b[0mcore.import_pose.import_pose: \u001b[0mFile 'inputs/glycans/2j88_antigen.pdb' automatically determined to be of type PDB\n", "\u001b[0mcore.io.pdb.pdb_reader: \u001b[0mParsing 0 .pdb records with unknown format to search for Rosetta-specific comments.\n", "\u001b[0mcore.io.util: \u001b[0mAutomatic glycan connection is activated.\n", "\u001b[0mcore.io.util: \u001b[0mStart reordering residues.\n", "\u001b[0mcore.io.util: \u001b[0mCorrected glycan residue order (internal numbering): []\n", "\u001b[0mcore.io.util: \u001b[0m\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 56\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 57\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 13 299\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 13 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 299 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 13 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 299 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 175 187\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 175 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 187 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 175 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 187 CYD\n" ] } ], "source": [ "pose_complex = pose_from_pdb(\"inputs/glycans/2j88_complex.pdb\")\n", "pose_complex_original = pose_complex.clone()\n", "\n", "pose_antigen = pose_from_pdb(\"inputs/glycans/2j88_antigen.pdb\")\n", "pose_antigen_original = pose_antigen.clone()" ] }, { "cell_type": "code", "execution_count": 47, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.carbohydrates.CreateGlycanSequonMover: \u001b[0mSequons created\n", "\u001b[0mprotocols.carbohydrates.CreateGlycanSequonMover: \u001b[0mApplying...\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mSCOREFUNCTION: \u001b[32mref2015\u001b[0m\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mThe -include_sugars flag was used with no sugar_bb weight set in the weights file. Setting sugar_bb weight to 1.0 by default.\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/elec_cp_reps.dat\n", "\u001b[0mcore.scoring.elec.util: \u001b[0mRead 40 countpair representative atoms\n", "\u001b[0mcore.pack.dunbrack.RotamerLibrary: \u001b[0mshapovalov_lib_fixes_enable option is true.\n", "\u001b[0mcore.pack.dunbrack.RotamerLibrary: \u001b[0mshapovalov_lib::shap_dun10_smooth_level of 1( aka lowest_smooth ) got activated.\n", "\u001b[0mcore.pack.dunbrack.RotamerLibrary: \u001b[0mBinary rotamer library selected: /Users/jadolfbr/Library/Python/3.6/lib/python/site-packages/pyrosetta-2019.39+release.93456a567a8-py3.6-macosx-10.6-intel.egg/pyrosetta/database/rotamer/shapovalov/StpDwn_0-0-0/Dunbrack10.lib.bin\n", "\u001b[0mcore.pack.dunbrack.RotamerLibrary: \u001b[0mUsing Dunbrack library binary file '/Users/jadolfbr/Library/Python/3.6/lib/python/site-packages/pyrosetta-2019.39+release.93456a567a8-py3.6-macosx-10.6-intel.egg/pyrosetta/database/rotamer/shapovalov/StpDwn_0-0-0/Dunbrack10.lib.bin'.\n", "\u001b[0mcore.pack.dunbrack.RotamerLibrary: \u001b[0mDunbrack 2010 library took 0.246233 seconds to load from binary\n", "\u001b[0mprotocols.calc_taskop_movers.CreateSequenceMotifMover: \u001b[0mFocus [129, 130, 131]Applying...\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 393 rotamers at 9 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating PDInteractionGraph\n" ] } ], "source": [ "start = pose_antigen.pdb_info().pdb2pose(\"A\", 143)\n", "end = pose_antigen.pdb_info().pdb2pose(\"A\", 145)\n", "select = ResidueIndexSelector()\n", "select.set_index(start)\n", "\n", "motif = ResidueIndexSelector()\n", "motif.set_index_range(start, end)\n", "\n", "sequon_creator = CreateGlycanSequonMover(select)\n", "sequon_creator.apply(pose_antigen)" ] }, { "cell_type": "code", "execution_count": 53, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "original FWD\n", "designed NWS\n" ] } ], "source": [ "seq_metric = SequenceMetric(motif)\n", "print(\"original\", seq_metric.calculate(pose_antigen_original))\n", "print(\"designed\", seq_metric.calculate(pose_antigen))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Did we design in the sequon? (`N[^P][S/T]`)\n", "\n", "What do our scores look like?\n", "Note we are packing neighbors here, so our score does not totally say that we actually reduced the score." ] }, { "cell_type": "code", "execution_count": 54, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mSCOREFUNCTION: \u001b[32mref2015\u001b[0m\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mThe -include_sugars flag was used with no sugar_bb weight set in the weights file. Setting sugar_bb weight to 1.0 by default.\n", "original -643.5306907221014\n", "designed -654.908075182865\n" ] } ], "source": [ "score = get_score_function()\n", "print(\"original\", score.score(pose_antigen_original))\n", "print(\"designed\", score.score(pose_antigen))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Note we are packing neighbors here, so our score does not totally say that we actually reduced the score." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### A1.2 Design using the `N[^P][T]` motif\n", "\n", "\n", "This motif has been shown to have higher occupancy of the glycosation site with glycans in the resulting protein. Glycosylation is not 100% in some cases at some positions for (currently) unknown reasons, but it has been shown that a THR at the +2 site results in higher occupancy. If we were creating a drug, we can use chromatography during protein isolation to choose peaks which include our glycan. \n", "\n", "Here, we are using the [-] notation as to not actually design the second position. We will use what is in the native protein here. Again, you may run the XML (`inputs/glycans/tutA12.xml`) within PyRosetta replacing start and end with: `start=143A end=145A`, or run the code below. \n", "\n", "The key XML line is here:\n", " \n" ] }, { "cell_type": "code", "execution_count": 56, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mSCOREFUNCTION: \u001b[32mref2015\u001b[0m\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mThe -include_sugars flag was used with no sugar_bb weight set in the weights file. Setting sugar_bb weight to 1.0 by default.\n", "\u001b[0mprotocols.calc_taskop_movers.CreateSequenceMotifMover: \u001b[0mFocus [129, 130, 131]Applying...\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 108 rotamers at 9 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating PDInteractionGraph\n" ] } ], "source": [ "from rosetta.protocols.calc_taskop_movers import *\n", "\n", "general_motif_creator = CreateSequenceMotifMover(select)\n", "general_motif_creator.set_motif(\"N[-]T\")\n", "pose_antigen_d1 = pose_antigen.clone()\n", "\n", "general_motif_creator.apply(pose_antigen)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Was the sequon successfully designed?" ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "original FWD\n", "designed1 NWS\n", "designed NWT\n" ] } ], "source": [ "print(\"original\", seq_metric.calculate(pose_antigen_original))\n", "print(\"designed1\", seq_metric.calculate(pose_antigen_d1))\n", "print(\"designed2\", seq_metric.calculate(pose_antigen))" ] }, { "cell_type": "code", "execution_count": 58, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "original -643.5306907221014\n", "designed1 -654.9080751828645\n", "designed2 -653.9549454739722\n" ] } ], "source": [ "print(\"original\", score.score(pose_antigen_original))\n", "print(\"designed1\", score.score(pose_antigen_d1))\n", "print(\"designed2\", score.score(pose_antigen)) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### A2. Adding a man5 glycan: \n", "\n", "__SimpleGlycosylateMover__\n", "\n", "Now, we will expand on our first tutorial by glycosylating afterward. We will use the common name for a man5 sugar, which is a high-mannose branching glycan of 7 sugar residues (and 5 mannoses). You can use a few common names to make glycosylation easier, or an IUPAC string, or a file that has the IUPAC string in the first name of the file. Common names include man5,man7,man9 and a few others. You can find these in\n", "\n", "\tdatabase/chemical/carbohydrates/common_glycans\n", "\n", "The IUPAC nomenclature of the man5 is as follows: \n", "\n", "\n", "\t\ta-D-Manp-(1->3)-[a-D-Manp-(1->3)-[a-D-Manp-(1->6)]-a-D-Manp-(1->6)]\n", "\t\t -b-D-Manp-(1->4)-b-D-GlcpNAc-(1->4)-b-D-GlcpNAc-\n", "\n", "\n", "\n", "More information on IUPAC nomenclature of sugar trees is here: . There is also a very detailed README in the common glycan directory for your reference.\n", "\n", "Note that within the `SimpleGlycosylateMover` you may also give multiple glycans using the `glycans` option, which will randomly choose a glycan tree to use for glycosylation from the list given. Glycosylation is not deterministic in that you always get a man5 at a particular position and is influenced by a great deal of structural biology that is not yet fully determined. For now, since we are aiming to create a drug and purifying our result, using a man5 is sufficient. Once again, either run the xml (inputs/glycans/tutA2.xml) as before or run the code below. This takes about 15 seconds. \n", "\n", "XML syntax is as follows:\n", "\n", "\t\n", "\n" ] }, { "cell_type": "code", "execution_count": 59, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mN-glycan_core N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mglycan_core N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mman3 N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mman5 man5.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mman9 man9.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mbisected_N-glycan_core bisected_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mbisected_fucosylated_N-glycan_core bisected_fucosylated_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_1_O-glycan core_1_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_2_O-glycan core_2_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_3_O-glycan core_3_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_4_O-glycan core_4_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_5_O-glycan core_5_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_6_O-glycan core_6_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_7_O-glycan core_7_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mcore_8_O-glycan core_8_O-glycan.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mfucosylated_N-glycan_core fucosylated_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mhigh-mannose_N-glycan_core high-mannose_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mhybrid_N-glycan_core hybrid_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mhybrid_bisected_N-glycan_core hybrid_bisected_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mhybrid_bisected_fucosylated_N-glycan_core hybrid_bisected_fucosylated_N-glycan_core.iupac\n", "\u001b[0mcore.chemical.carbohydrates.database_io: \u001b[0mhybrid_fucosylated_N-glycan_core hybrid_fucosylated_N-glycan_core.iupac\n", "\u001b[0mcore.pose.carbohydrates.util: \u001b[0mDelete to residue is not carbohydrate and not a branch point. Nothing to be done.\n", "\u001b[0mprotocols.carbohydrates.SimpleGlycosylateMover: \u001b[0mGlycosylating at 129 : 129 a-D-Manp-(1->3)-[a-D-Manp-(1->3)-[a-D-Manp-(1->6)]-a-D-Manp-(1->6)]-b-D-Manp-(1->4)-b-D-GlcpNAc-(1->4)-b-D-GlcpNAc-\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mappending residue by a chemical bond in the foldtree: 320 ->4)-beta-D-Glcp:2-AcNH anchor: ND2 129 root: C1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mappending residue by a chemical bond in the foldtree: 324 ->3)-alpha-D-Manp:->6)-branch anchor: O6 322 root: C1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mappending residue by a chemical bond in the foldtree: 326 ->4)-alpha-D-Manp:non-reducing_end anchor: O6 324 root: C1\n", "\u001b[0mcore.pose.carbohydrates.util: \u001b[0mGlycosylated pose with a(n) a-D-Manp-(1->3)-[a-D-Manp-(1->3)-[a-D-Manp-(1->6)]-a-D-Manp-(1->6)]-b-D-Manp-(1->4)-b-D-GlcpNAc-(1->4)-b-D-GlcpNAc--ND2ASN129 bond.\n", "\u001b[0mcore.conformation.carbohydrates.GlycanTreeSet: \u001b[0mSetting up Glycan Trees\n", "\u001b[0mcore.conformation.carbohydrates.GlycanTreeSet: \u001b[0mFound 1 glycan trees.\n" ] } ], "source": [ "glycosylate = SimpleGlycosylateMover()\n", "glycosylate.set_glycosylation(\"man5\")\n", "glycosylate.set_residue_selector(select)\n", "glycosylate.apply(pose_antigen)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Take a look at the structure using the PyMol Mover or output it. What does the structure look like? The SimpleGlycosylateMover simply adds in glycans at ideal values for rings and backbones. \n", "Lets optimize this and model some glycans!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### A3. Modeling glycans\n", "\n", "`GlycanResidueSelector` and the `GlycanModeler`\n", "\n", "We will run the previous tutorials in a single rosetta script where we end with modeling the glycan residues. We use a very short run time and nstruct, so results will not be as clean as they would otherwise, but this should give you an idea of how all this works. Typically, we would model different positions of potential glycosylations, but here to save time, we will simply continue to build and model the glycan position we started with. Output files have been provided for you if you wish to use these. We will not be giving the mover a residue selector as it uses all glycans by default, but you can use the `GlycanResidueSelector` to choose specific trees or even glycan residues within those trees to model as we have seen in the previous tutorial.\n", "\n", "This takes about 380 seconds to run. You can either use the XML (`inputs/glycans/tutA3.xml`) or the code below. Output 10 structures.\n", "\n", "XML Syntax:\n", "\n", "\t\n" ] }, { "cell_type": "code", "execution_count": 61, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mSCOREFUNCTION: \u001b[32mref2015\u001b[0m\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mThe -include_sugars flag was used with no sugar_bb weight set in the weights file. Setting sugar_bb weight to 1.0 by default.\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRandomizing glycan torsions\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mSmallest glycan layer: 0\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mLargest glycan layer: 4\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mNtrees to model 1\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mStarting Score: 7529.89\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mstarting round 1\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mGoing in the forward direction\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mModeling up to max end: 4\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mVirtualizing new foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning the GlycanSampler on layer [ start -> end (including) ]: 0 1\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mApplying normal protocol\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 2 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mSetting up Tree-based packing.\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -291.629\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mTotal Rounds = 50 ( 2 glycans * 25 )\n", "\u001b[0mcore.scoring.ScoringManager: \u001b[0mshould be setting up for sampling...\n", "\u001b[0mcore.scoring.carbohydrates.CHIEnergyFunction: \u001b[0mSetting up chi sampling\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 4 which has a population of 0.0169\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.0124\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 321\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 320\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 2 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 7 which has a population of 0.0058\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.0124\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 30\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 4 which has a population of 0.0169\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 320\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 2 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 40\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.3218\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.9876\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 50\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -291.629\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -649.102\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 77/80\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 327 rotamers at 3 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 78/80\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 327 rotamers at 3 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -650.195\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -650.195\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 215; accepts= 0.8698; energy_drop/trial= -1.66087\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mDe-Virtualizing current foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mVirtualizing new foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning the GlycanSampler on layer [ start -> end (including) ]: 1 2\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 2 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -469.11\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mTotal Rounds = 24 ( 2 glycans * 12 )\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.0124\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.0124\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.9876\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 321\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 2 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.9876\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 321\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 2 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 322\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-beta-D-Manp-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 321\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 2 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -469.11\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -642.762\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 72/80\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 336 rotamers at 3 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 75/80\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 336 rotamers at 3 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -643.293\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -644.723\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 189; accepts= 0.8677; energy_drop/trial= -0.92160\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning full window protocol sampling\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 3 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -644.723\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 24\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.4755\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.4755\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 320\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 3 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 321\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 2 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 322\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 322\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -644.723\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -648.181\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 110/120\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 337 rotamers at 4 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 106/120\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 337 rotamers at 4 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -645.986\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -648.769\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 269; accepts= 0.8699; energy_drop/trial= -0.00470\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mDe-Virtualizing current foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mVirtualizing new foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning the GlycanSampler on layer [ start -> end (including) ]: 2 3\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 3 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: 1719.85\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 24\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 323\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 323\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 323\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: 1719.85\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -630.877\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 106/120\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 427 rotamers at 4 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 104/120\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 427 rotamers at 4 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -631.385\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -632.623\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 269; accepts= 0.8625; energy_drop/trial= -8.74065\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning full window protocol sampling\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 5 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -632.623\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 36\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.9876\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.0124\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 320\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 5 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-beta-D-Manp-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.4755\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-beta-D-Manp-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-beta-D-Manp-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 30\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 320\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 5 residues\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.9876\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.0124\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -632.623\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -634.56\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 183/200\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 591 rotamers at 6 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 189/200\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 591 rotamers at 6 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mcore.conformation.Conformation: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Branch 1does not have enough heavy atoms about its connection to define a torsion angle!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -633.67\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -634.869\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 441; accepts= 0.8912; energy_drop/trial= -0.00232\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mDe-Virtualizing current foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mVirtualizing new foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning the GlycanSampler on layer [ start -> end (including) ]: 3 4\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 4 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: 2843.92\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 36\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mcore.scoring.ScoringManager: \u001b[0mshould be setting up for sampling...\n", "\u001b[0mcore.scoring.carbohydrates.OmegaPreferencesFunction: \u001b[0mSetting up omega preference sampling\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 323\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 30\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 325\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: 2843.92\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -451.244\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 160/160\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 752 rotamers at 5 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 160/160\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 752 rotamers at 5 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -452.757\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -452.757\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 357; accepts= 0.9636; energy_drop/trial= -9.23223\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning full window protocol sampling\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 7 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -452.757\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 48\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-beta-D-Manp-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 323\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-beta-D-Manp-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 326\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 322\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 5 residues\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.4755\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 3 which has a population of 0.1613\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 30\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?->4)-beta-D-GlcpNAc- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.9876\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 7 which has a population of 0.0058\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 40\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 324\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 3 residues\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -452.757\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -630.598\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 253/280\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 241/280\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -630.233\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -630.645\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 613; accepts= 0.8630; energy_drop/trial= -0.27731\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mDe-Virtualizing current foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mVirtualizing new foliage layer\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning the GlycanSampler on layer [ start -> end (including) ]: 4 4\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 2 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -630.645\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 48\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 325\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 30\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mLast mover failed. Continueing!\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 325\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 325\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 40\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 2 which has a population of 0.4401\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing from carbohydrate root: 325\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeMinMover: \u001b[0mMinimizing 1 residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -630.645\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -636.196\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 80/80\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 498 rotamers at 3 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.simple_moves.BackboneMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m no movable positions in ShearMover!\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 80/80\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 498 rotamers at 3 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -636.196\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -636.196\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 209; accepts= 0.9139; energy_drop/trial= -0.01769\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning full window protocol sampling\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0minitializing objects\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mModeling 7 glycan residues\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mInitialized\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mstarting energy: -636.196\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mCustom total rounds: 24\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->3)-alpha-D-Manp-(?->6)-beta-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 10\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1078 rotamers at 7 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0mRound: 20\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling ->4)-beta-D-GlcpNAc-(?ASN linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 6 which has a population of 0.0063\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->3)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 1\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling alpha-D-Manp-(?->6)-alpha-D-Manp- linkage\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mSampling conformer 1 which has a population of 0.5599\n", "\u001b[0mprotocols.carbohydrates.LinkageConformerMover: \u001b[0mComplete\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy start: -636.196\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final: -636.558\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 265/280\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 238/280\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy final post min: -636.289\n", "\u001b[0mprotocols.carbohydrates.GlycanSampler: \u001b[0menergy recov post min: -636.824\n", "\u001b[0mprotocols.moves.TrialCounter: \u001b[0m unk trials= 589; accepts= 0.8879; energy_drop/trial= 0.00060\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mround 1 -636.527\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0maccepted 1\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0m\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mround 1 -636.527 final\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0m\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mFinal Score: -636.527\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mRunning final min-pack.\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 247/280\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mShear accepts: 238/280\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mmin accept: 1\n", "\u001b[0mcore.pack.task: \u001b[0mPacker task: initialize from command line()\n", "\u001b[0mcore.pack.rotamer_set.RotamerSet_: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m including current in order to get at least 1 rotamer !!!!!! 129 ASN:N-glycosylated\n", "\u001b[0mcore.pack.pack_rotamers: \u001b[0mbuilt 1079 rotamers at 8 positions.\n", "\u001b[0mcore.pack.interaction_graph.interaction_graph_factory: \u001b[0mInstantiating DensePDInteractionGraph\n", "\u001b[0mprotocols.carbohydrates.util: \u001b[0mpack accepted: 1\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mStart- : 7529.89\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mPre - : -636.527\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mPost - : -635.756\n", "\u001b[0mprotocols.carbohydrates.GlycanTreeModeler: \u001b[0mFinal- : -636.711\n" ] } ], "source": [ "#Here, we setup some defaults that will be set in master shortly. \n", "# Note that the name of the GTM will change to simply GlycanModeler within the next week or two.\n", "\n", "modeler = GlycanTreeModeler()\n", "modeler.set_hybrid_protocol(True)\n", "modeler.set_use_shear(True)\n", "modeler.set_use_gaussian_sampling(True)\n", "modeler.apply(pose_antigen)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Run this code for a total of 10 nstruct. Compare energies and take a look at the lowest one. Scores and pdbs are also available in `expected_outputs/glycans` with prefix `tutA3`" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "How does it look? Load the native (complexed structure) into pymol as well. Would this glycan block this particular antibody? Where else could we place a glycan?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Tutorial B: Using Glycan Density\n", "\n", "In this tutorial we will load a pdb directly into Rosetta with sugars already present. The config for this has been provided for you. We will exclusively be using the XML interface to PyRosetta for these tutorials, however, using the knowledge gained in the Density and Symmetry tutorials, all of this is available code-wise as well, but the XML interface greatly simplifies this task.\n", "\n", "\n", "The glycan tree that we will be working with is 5 residues long. I use coot to look at density maps. \n", "Density maps were generated by downloading the cif file from the PDB and using PHENIX maps and default `maps.params`. This command was used to generate them. This is covered in the Working-With-Density tutorial:\n", "\n", "\t\tphenix.maps inputs/glycans/4do4.pdb inputs/glycans/4do4-sf.cif \n", "\n", "The density map generated is too large to be distributed with the rest of the tuorial, so I have uploaded it to Google Drive for you to download. Please download and place it in your working directory, if it's not already there.\n", "\n", "The tutorial will use a density map that is small enough to fit into the repo, but you are welcome to use 4d04 instead with the downloaded map.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### B1. Calculating Density Fit\n", "\n", "\n", "Although a structure may be solved with high resolution, not all solved residues may fit the density well. A structure from the PDB is still a model afterall, informed through experimentation. This is especially true of glycan residues, which are fairly mobile. Crystal contacts of neighboring proteins help to reduce the movemment of glycans and may help to induce a state that can be solved more easily given high-resolution density. In this tutorial, we will be using Rosetta to determine how well a residue fits into the given density. There are methods to do this in the coot program, but we want to be able to do this for any structure in a streamlined way - especially if we need to calculate RMSDs on only well-fitting glycan residues. The methods we will be employing in Rosetta are based on Frank Dimaio's work with Rosetta density. \n", "\n", "To do this, we will once again be employing the SimpleMetric system. In this case, we use the `PerResidueDensityFitMetric`, which is a PerResidueRealMetric. This type of SimpleMetric calculates a particular value for each residue given a residue selector. Very useful here. \n", "\n", "We will also be employing the DensityFitResidueSelector, which uses the metric. Since this is a fairly slow metric, we will use in-built functionality for using our calculated values from the metric, which are stored in the pose. We will then use the SelectedResidueCountMetric to determine how many residues have great fit. In later tutorials, we will be using the RMSDMetric with this selector in order to calculate RMSD on well-fitting glycan residues. \n", "\n", "\tResidues higher than .8 are great fit to density.\n", "\tResidues between .6 - .8 are good fit to density\n", "\tResidues below .4 fit to density are BAD fits\n", "\n", "\n", "The XML is `inputs/glycans/tutB1.xml`\n", "\n", "The key XML lines are:\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\n", "Run the xml and while it is running, take a look at the XML (runtime is about 80 seconds). It is fairly complicated and we will be building on it during the rest of these tutorials. Note that we first define the density metric, and then we use it within the selector. At the bottom, we add these to our set of native_metrics. What other metrics are we using?\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "branch=\"200A\"\n", "map_file=\"inputs/1jnd_2mFo-DFc_map.ccp4\"\n", "symmdef=\"inputs/1jnd_crys.symm\"\n", "\n", "xml = f'''\n", "\n", " \n", " \n", "\n", " NEEDED FOR CACHING density fit info\n", " \n", " \n", " \n", " \n", " \n", " \n", " END Density Fit Setup\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n", " \n", "\n", " \n", "\n", " \n", " \n", " \n", " \n", " \n", "\n", " \n", " \n", " \n", " \n", " \n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n", "'''" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.chemical.GlobalResidueTypeSet: \u001b[0mFinished initializing fa_standard residue type set. Created 1251 residue types\n", "\u001b[0mcore.chemical.GlobalResidueTypeSet: \u001b[0mTotal time to initialize 1.1842 seconds.\n", "\u001b[0mcore.import_pose.import_pose: \u001b[0mFile 'inputs/1jnd_refined.pdb.gz' automatically determined to be of type PDB\n", "\u001b[0mcore.io.pdb.pdb_reader: \u001b[0mParsing 411 .pdb records with unknown format to search for Rosetta-specific comments.\n", "\u001b[0mcore.io.util: \u001b[0mAutomatic glycan connection is activated.\n", "\u001b[0mcore.io.util: \u001b[0mStart reordering residues.\n", "\u001b[0mcore.io.util: \u001b[0mCorrected glycan residue order (internal numbering): [401, 402, 403, 404]\n", "\u001b[0mcore.io.util: \u001b[0m\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m missing density found at residue (rosetta number) 140\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Residue connection id changed when creating a new residue at seqpos 141\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m ResConnID info stored on the connected residue (residue 0) is now out of date!\n", "\u001b[0mcore.conformation.Residue: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Connection atom name (in src): C\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Glc401 has an unfavorable ring conformation; the coordinates for this input structure may have been poorly assigned.\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Glc402 has an unfavorable ring conformation; the coordinates for this input structure may have been poorly assigned.\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Man403 has an unfavorable ring conformation; the coordinates for this input structure may have been poorly assigned.\n", "\u001b[0mcore.io.pose_from_sfr.PoseFromSFRBuilder: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Man404 has an unfavorable ring conformation; the coordinates for this input structure may have been poorly assigned.\n", "\u001b[0mcore.chemical.AtomICoor: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m IcoorAtomID::atom_id(): Cannot get atom_id for POLYMER_LOWER of residue ->4)-beta-D-Glcp:2-AcNH 401. Returning BOGUS ID instead.\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 5 32\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 5 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 32 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 5 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 32 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mFound disulfide between residues 302 385\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 302 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 385 CYS\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 302 CYD\n", "\u001b[0mcore.conformation.Conformation: \u001b[0mcurrent variant for 385 CYD\n", "\u001b[0mcore.conformation.carbohydrates.GlycanTreeSet: \u001b[0mSetting up Glycan Trees\n", "\u001b[0mcore.conformation.carbohydrates.GlycanTreeSet: \u001b[0mFound 1 glycan trees.\n" ] } ], "source": [ "pose = pose_from_pdb(\"inputs/1jnd_refined.pdb.gz\")" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mGenerating XML Schema for rosetta_scripts...\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0m...done\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mInitializing schema validator...\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0m...done\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mValidating input script...\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0m...done\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mParsed script:\n", "\n", "\t\n", "\t\n", "\t\t\n", "\t\n", "\t\n", "\t\t\n", "\t\n", "\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\n", "\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\n", "\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\n", "\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\t\n", "\t\n", "\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mSCOREFUNCTION: \u001b[32mref2015\u001b[0m\n", "\u001b[0mcore.scoring.etable: \u001b[0mStarting energy table calculation\n", "\u001b[0mcore.scoring.etable: \u001b[0msmooth_etable: changing atr/rep split to bottom of energy well\n", "\u001b[0mcore.scoring.etable: \u001b[0msmooth_etable: spline smoothing lj etables (maxdis = 6)\n", "\u001b[0mcore.scoring.etable: \u001b[0msmooth_etable: spline smoothing solvation etables (max_dis = 6)\n", "\u001b[0mcore.scoring.etable: \u001b[0mFinished calculating energy tables.\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/hbonds/ref2015_params/HBPoly1D.csv\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/hbonds/ref2015_params/HBFadeIntervals.csv\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/hbonds/ref2015_params/HBEval.csv\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/hbonds/ref2015_params/DonStrength.csv\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/hbonds/ref2015_params/AccStrength.csv\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/rama/fd/all.ramaProb\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/rama/fd/prepro.ramaProb\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/omega/omega_ppdep.all.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/omega/omega_ppdep.gly.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/omega/omega_ppdep.pro.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/omega/omega_ppdep.valile.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/P_AA_pp/P_AA\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/P_AA_pp/P_AA_n\n", "\u001b[0mcore.scoring.P_AA: \u001b[0mshapovalov_lib::shap_p_aa_pp_smooth_level of 1( aka low_smooth ) got activated.\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/P_AA_pp/shapovalov/10deg/kappa131/a20.prop\n", "\u001b[0mcore.scoring.ScoreFunctionFactory: \u001b[0mThe -include_sugars flag was used with no sugar_bb weight set in the weights file. Setting sugar_bb weight to 1.0 by default.\n", "\u001b[0mcore.scoring.ScoringManager: \u001b[0mCreating CHI Energy Function.\n", "\u001b[0mcore.scoring.ScoringManager: \u001b[0mCreating carbohydrate omega preferences function.\n", "\u001b[0mcore.scoring.etable: \u001b[0mStarting energy table calculation\n", "\u001b[0mcore.scoring.etable: \u001b[0msmooth_etable: changing atr/rep split to bottom of energy well\n", "\u001b[0mcore.scoring.etable: \u001b[0msmooth_etable: spline smoothing lj etables (maxdis = 6)\n", "\u001b[0mcore.scoring.etable: \u001b[0msmooth_etable: spline smoothing solvation etables (max_dis = 6)\n", "\u001b[0mcore.scoring.etable: \u001b[0mFinished calculating energy tables.\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/PairEPotential/pdb_pair_stats_fine\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/InterchainPotential/interchain_env_log.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/InterchainPotential/interchain_pair_log.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/EnvPairPotential/env_log.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/EnvPairPotential/cbeta_den.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/EnvPairPotential/pair_log.txt\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/EnvPairPotential/cenpack_log.txt\n", "\u001b[0mcore.scoring.ramachandran: \u001b[0mshapovalov_lib::shap_rama_smooth_level of 4( aka highest_smooth ) got activated.\n", "\u001b[0mbasic.io.database: \u001b[0mDatabase file opened: scoring/score_functions/rama/shapovalov/kappa25/all.ramaProb\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector tree\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector tree\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector fits8\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector fits6\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector fits4\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector layer01\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector tree\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector tree\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector root\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector tree\n", "\u001b[0mcore.select.residue_selector.util: \u001b[0mFound residue selector root\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Setting weight of master jump ( jump-id=1 ) to 1.0 (was undefined)\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Setting jump_group JUMPGROUP1: [master 1] 6:0 7:0 2:0 3:0 4:0 5:0 8:0 9:0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Setting weight of master jump ( jump-id=10 ) to 1.0 (was undefined)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Setting jump_group JUMPGROUP2: [master 10] 20:0 22:0 12:0 14:0 16:0 18:0 24:0 26:0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0msymmetry name: 1jnd_P_32_2_1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mnumber of subunits: 9\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mnumber of interfaces: 1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mscore subunit number: VRT_0_0_0_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0manchor the subunits at residue: 1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_0_0_0_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: 1 0 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0 1 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 73.8756 49.0221 11.025\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_0_0_0_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: 1 0 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0 1 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 72.8756 48.0221 10.025\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_1_1_0_n1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0.866025 -0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 26.9528 39.4671 -18.9707\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_1_1_0_n1_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0.866025 -0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 28.3188 39.1011 -19.9707\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_2_1_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 -0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0.866025 -0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 58.6891 3.60826 41.0206\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_2_1_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 -0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0.866025 -0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 58.3231 4.97428 40.0206\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_3_0_n1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 58.6891 -3.60826 -11.025\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_3_0_n1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 58.3231 -4.97428 -10.025\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_3_1_0_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 111.862 88.4892 -11.025\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_3_1_0_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: 0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 111.496 87.1232 -10.025\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_4_2_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 -0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 80.1253 52.6304 48.9664\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_4_2_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 -0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 81.4913 52.9964 49.9664\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_4_2_1_n1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 -0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 80.1253 52.6304 -41.0206\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_4_2_1_n1_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: -0.5 -0.866025 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0.866025 0.5 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 81.4913 52.9964 -40.0206\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_5_0_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: 1 0 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0 -1 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 20.7031 43.0754 18.9707\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_5_0_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: 1 0 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0 -1 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 19.7031 44.0754 19.9707\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_5_1_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: 1 0 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0 -1 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 127.048 43.0754 18.9707\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0m Virtual coordinate system VRT_5_1_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx: 1 0 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my: -0 -1 0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0morigin: 126.048 44.0754 19.9707\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mDof for jump: 1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx 0:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my 0:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mz 0:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx_angle 1:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my_angle 1:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mz_angle 1:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mDof for jump: 10\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx 1:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my 1:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mz 1:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mx_angle 0:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0my_angle 0:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mz_angle 0:0,0:0,0 n2c\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_0_0_0_0_to_com VRT_0_0_0_0 VRT_0_0_0_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_0_0_0_0_to_subunit VRT_0_0_0_0_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_1_1_0_n1 VRT_0_0_0_0 VRT_1_1_0_n1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_1_1_0_n1_to_com VRT_1_1_0_n1 VRT_1_1_0_n1_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_1_1_0_n1_to_subunit VRT_1_1_0_n1_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_2_1_1_0 VRT_0_0_0_0 VRT_2_1_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_2_1_1_0_to_com VRT_2_1_1_0 VRT_2_1_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_2_1_1_0_to_subunit VRT_2_1_1_0_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_3_0_n1_0 VRT_0_0_0_0 VRT_3_0_n1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_3_0_n1_0_to_com VRT_3_0_n1_0 VRT_3_0_n1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_3_0_n1_0_to_subunit VRT_3_0_n1_0_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_3_1_0_0 VRT_0_0_0_0 VRT_3_1_0_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_3_1_0_0_to_com VRT_3_1_0_0 VRT_3_1_0_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_3_1_0_0_to_subunit VRT_3_1_0_0_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_4_2_1_0 VRT_0_0_0_0 VRT_4_2_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_4_2_1_0_to_com VRT_4_2_1_0 VRT_4_2_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_4_2_1_0_to_subunit VRT_4_2_1_0_base SUBUNIT\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_4_2_1_n1 VRT_0_0_0_0 VRT_4_2_1_n1\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_4_2_1_n1_to_com VRT_4_2_1_n1 VRT_4_2_1_n1_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_4_2_1_n1_to_subunit VRT_4_2_1_n1_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_5_0_1_0 VRT_0_0_0_0 VRT_5_0_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_5_0_1_0_to_com VRT_5_0_1_0 VRT_5_0_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_5_0_1_0_to_subunit VRT_5_0_1_0_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_5_1_1_0 VRT_0_0_0_0 VRT_5_1_1_0\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_5_1_1_0_to_com VRT_5_1_1_0 VRT_5_1_1_0_base\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mJump JUMP_5_1_1_0_to_subunit VRT_5_1_1_0_base SUBUNIT\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mInclude subunit:\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mOutput subunit:\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mSlideType: RANDOM\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mSlideCriteriaType: CEN_DOCK_SCORE\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mSlideCriteriaVal: AUTOMATIC\n", "\u001b[0mcore.conformation.symmetry.SymmData: \u001b[0mSlideOrder: none\n", "\u001b[0mprotocols.simple_moves.symmetry.SetupForSymmetryMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m Symmetry is being set globally, which can cause conflicts with certain protocols, but is necessary for backwards compatibility.\n", "\u001b[0mprotocols.simple_moves.symmetry.SetupForSymmetryMover: \u001b[0m\u001b[1m[ WARNING ]\u001b[0m If you are experiencing unexpected behaviour related to symmetry, try setting set_global_symmetry_at_parsetime=0 in .\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"setup_symm\" of type SetupForSymmetry\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"loaddens\" of type LoadDensityMap\n", "\u001b[0mprotocols.electron_density.util: \u001b[0mParsing SetupForDensityScoringMover----\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"setupdens\" of type SetupForDensityScoring\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric PerResidueDensityFitMetric.\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"native_metrics\" of type RunSimpleMetrics\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric PerResidueGlycanLayerMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResiduesPyMOLMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResiduesPyMOLMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResiduesMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResiduesMetric.\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"selections\" of type RunSimpleMetrics\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResidueCountMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResidueCountMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResidueCountMetric.\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric SelectedResidueCountMetric.\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"counts\" of type RunSimpleMetrics\n", "\u001b[0mcore.simple_metrics.util: \u001b[0mAdded simple metric TimingProfileMetric.\n", "\u001b[0mprotocols.rosetta_scripts.RosettaScriptsParser: \u001b[0mDefined mover named \"timings\" of type RunSimpleMetrics\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0mParsedProtocol mover with the following movers and filters\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0madded mover \"loaddens\" with filter \"true_filter\"\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0madded mover \"setupdens\" with filter \"true_filter\"\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0madded mover \"selections\" with filter \"true_filter\"\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0madded mover \"native_metrics\" with filter \"true_filter\"\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0madded mover \"counts\" with filter \"true_filter\"\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0madded mover \"timings\" with filter \"true_filter\"\n" ] } ], "source": [ "mover = pyrosetta.rosetta.protocols.rosetta_scripts.XmlObjects.create_from_string(xml).get_mover(\"ParsedProtocol\")" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0m=======================BEGIN MOVER LoadDensityMap - loaddens=======================\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0mLoading Density Map\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0mLoading density mapinputs/1jnd_2mFo-DFc_map.ccp4\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m Setting resolution to AUTO\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m atom mask to 3.2A\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m CA mask to 6A\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m Read density map'inputs/1jnd_2mFo-DFc_map.ccp4'\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m extent: 218 x 229 x 265\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m origin: -66 x 82 x 209\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m altorigin: 0 x 0 x 0\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m grid: 360 x 360 x 288\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m celldim: 106.345 x 106.345 x 89.987\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m cellangles: 90 x 90 x 120\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0m voxel vol.: 0.0236128\n", "\u001b[0mcore.scoring.electron_density.ElectronDensity: \u001b[0mEffective resolution = 0.937365\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0m=======================BEGIN MOVER SetupForDensityScoring - setupdens=======================\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0m=======================BEGIN MOVER RunSimpleMetrics - selections=======================\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: PerResidueGlycanLayerMetric - calculating glycan_layer\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResiduesPyMOLMetric - calculating glycans_pymol_selection\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResiduesPyMOLMetric - calculating branch_pymol_selection\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResiduesMetric - calculating glycans_selection\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResiduesMetric - calculating branch_selection\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0m=======================BEGIN MOVER RunSimpleMetrics - native_metrics=======================\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: PerResidueDensityFitMetric - calculating res_density_fit\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0m=======================BEGIN MOVER RunSimpleMetrics - counts=======================\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResidueCountMetric - calculating tree_size_selection_count\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResidueCountMetric - calculating layer01_selection_count\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResidueCountMetric - calculating fit6_selection_count\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3100 A 0.85068\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0mGood Fit\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3101 A 0.834764\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0mGood Fit\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3102 A 0.761029\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0mGood Fit\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3103 A 0.649696\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0mGood Fit\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: SelectedResidueCountMetric - calculating fit8_selection_count\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3100 A 0.85068\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0mGood Fit\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3101 A 0.834764\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0mGood Fit\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3102 A 0.761029\n", "\u001b[0mcore.select.residue_selector.DensityFitResidueSelector: \u001b[0m3103 A 0.649696\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0m=======================BEGIN MOVER RunSimpleMetrics - timings=======================\n", "\u001b[0mprotocols.analysis.simple_metrics.RunSimpleMetricsMover: \u001b[0mRunning: TimingProfileMetric - calculating timing_profile\n", "\u001b[0mprotocols.rosetta_scripts.ParsedProtocol: \u001b[0msetting status to success\n" ] } ], "source": [ "mover.apply(pose)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Ok, take a look at the (tutB1) scorefile or pose.scores. You can use the scorefile.py script to output as tabs if you would like. How many residues have great fit to density (hint, look for `fit6_selection_count` and `fit8_selection_count` data terms)? Are there any residues that fit poorly into the density?" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'branch_pymol_selection': 'select rosetta_sele, (chain A and resid 200)', 'branch_selection': '200A', 'glycans_pymol_selection': 'select rosetta_sele, (chain A and resid 3100,3101,3102,3103)', 'glycans_selection': '3100A,3101A,3102A,3103A', 'fit6_selection_count': 4.0, 'fit8_selection_count': 2.0, 'glycan_layer_3100A': 0.0, 'glycan_layer_3101A': 1.0, 'glycan_layer_3102A': 2.0, 'glycan_layer_3103A': 3.0, 'layer01_selection_count': 2.0, 'native_res_density_fit_3100A': 0.850679730589326, 'native_res_density_fit_3101A': 0.8347641236660035, 'native_res_density_fit_3102A': 0.7610287167541959, 'native_res_density_fit_3103A': 0.6496957888062153, 'timing_profile': 1.0479055643081665, 'tree_size_selection_count': 4.0}\n" ] } ], "source": [ "print(pose.scores)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "4.0\n" ] } ], "source": [ "print(pose.scores[\"fit6_selection_count\"])" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2.0\n" ] } ], "source": [ "print(pose.scores[\"fit8_selection_count\"])" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "3100A 0.850679730589326\n", "3101A 0.8347641236660035\n", "3102A 0.7610287167541959\n", "3103A 0.6496957888062153\n" ] } ], "source": [ "import re\n", "\n", "\n", "for term in pose.scores:\n", " if re.search(\"native_res_density_fit\", term):\n", " print(term.split(\"_\")[-1], pose.scores[term])\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### B2. Refinement into density\n", "\t\t\n", "Here, we will be doing a short refinement protocol into the density, with its crystal symmetry. This is a short protocol, but will work for our purposes. For a much longer (albeit very similar) refinement protocol of the glycan and whole protein, see Frenz et al (referenced in the Working With Glycans tutorial). The full protocol used in this paper is included in the input files as cryoem_glycan_refinement.xml. Take a look and see how it compares to what we are doing here. As usual, output files are available. Runtime for all 10 structures is about 2 hours. \n", "\n", "The XML is `inputs/glycans/tutB2.xml`. Since this takes so long to run, we will skip running it here. \n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Take a look at the scorefile either using the scorefile.py script or through pandas: `expected_outputs/glycans/tutB2_score.sc`\n", "\n", "Are the density fit scores higher? How different are the RMSDs of the glycan residues? Take a look at the structure of the lowest energy - how different does it look? Are any new contacts created? Were we able to improve the density fit for some of those residues?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### B3. Denovo building into Density\n", "\t\n", "In this tutorial, we will be once again loading our crystal structure with density and symmetry. However, we will be randomizing the bb torsions and building the glycan out from scratch. In reality, we would have some idea of what glycan we are building and we would glycosylate the protein with the chemical motif we have figured out from means such as mass spec. We would then model the glycan to solve the crystal structure. With the new PackerPalette machinery in Rosetta and the ability to design glycans, we could actually build a protocol to sample chemical motifs of the glycans we are building out into the density, however, since this a very very large combinatorial problem, we should have some idea of what exists in the structure. \n", "\n", "\n", "We will first rebuild the glycan tree using the density as a guide, and then refine it further using what we learned in the previous tutorial. Note that like tutorial B2, this one takes a good long while (4 hours)\n", "\n", "The XML for this is `inputs/glycans/tutB3.xml`. Once again, we will not be running this here.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Take a look at the scorefile and PDBs. (prefix tutB3_)\n", "\n", "How are our RMSDs? Were we able to do enough sampling to get close to the native structure? Are the energies acceptable? Are there parts of the glycan that are closer to native than others? Why might this be? Is an nstruct of 10 enough??" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Thank you for doing this tutorial! I hope you learned a lot and are ready to work with these crazy carbohydrates! Cheers!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "< [RosettaCarbohydrates: Trees, Selectors and Movers](http://nbviewer.jupyter.org/github/RosettaCommons/PyRosetta.notebooks/blob/master/notebooks/13.01-Glycan-Trees-Selectors-and-Movers.ipynb) | [Contents](toc.ipynb) | [Index](index.ipynb) | [RNA in PyRosetta](http://nbviewer.jupyter.org/github/RosettaCommons/PyRosetta.notebooks/blob/master/notebooks/14.00-RNA-Basics.ipynb) >

\"Open" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.0" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 1 }