Starting phenix.real_space_refine (version: dev) on Fri Feb 24 17:04:25 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7t9i_25758/02_2023/7t9i_25758_trim.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7t9i_25758/02_2023/7t9i_25758.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.9 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7t9i_25758/02_2023/7t9i_25758.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7t9i_25758/02_2023/7t9i_25758.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7t9i_25758/02_2023/7t9i_25758_trim.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7t9i_25758/02_2023/7t9i_25758_trim.pdb" } resolution = 2.9 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.009 sd= 0.038 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Process input model Symmetric amino acids flipped Residue "A PHE 18": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "A TYR 65": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "A GLU 77": "OE1" <-> "OE2" Residue "B TYR 45": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B ASP 56": "OD1" <-> "OD2" Residue "B GLU 63": "OE1" <-> "OE2" Residue "B TYR 74": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "N ASP 50": "OD1" <-> "OD2" Residue "N ASP 106": "OD1" <-> "OD2" Residue "R PHE 286": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "R TYR 385": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "R PHE 434": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "R TYR 601": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "R TYR 643": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "R PHE 681": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "X ASP 33": "OD1" <-> "OD2" Residue "X ASP 378": "OD1" <-> "OD2" Residue "X TYR 391": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "Y GLU 10": "OE1" <-> "OE2" Residue "Y ASP 66": "OD1" <-> "OD2" Residue "Y ASP 322": "OD1" <-> "OD2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/chem_data/mon_lib" Total number of atoms: 23529 Number of models: 1 Model: "" Number of chains: 11 Chain: "A" Number of atoms: 1205 Number of conformers: 1 Conformer: "" Number of residues, atoms: 86, 1205 Classifications: {'peptide': 86} Link IDs: {'PTRANS': 6, 'TRANS': 79} Unresolved non-hydrogen bonds: 30 Unresolved non-hydrogen angles: 36 Unresolved non-hydrogen dihedrals: 28 Planarities with less than four sites: {'PHE:plan': 2, 'GLU:plan': 1, 'GLN:plan1': 1} Unresolved non-hydrogen planarities: 20 Chain: "B" Number of atoms: 1636 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 1636 Classifications: {'peptide': 109} Modifications used: {'NH3': 1} Link IDs: {'PTRANS': 7, 'TRANS': 101} Unresolved non-hydrogen bonds: 25 Unresolved non-hydrogen angles: 30 Unresolved non-hydrogen dihedrals: 20 Unresolved non-hydrogen chiralities: 1 Planarities with less than four sites: {'GLN:plan1': 1, 'GLU:plan': 2, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 17 Chain: "N" Number of atoms: 1893 Number of conformers: 1 Conformer: "" Number of residues, atoms: 126, 1893 Classifications: {'peptide': 126} Modifications used: {'NH3': 1} Link IDs: {'CIS': 1, 'PTRANS': 5, 'TRANS': 119} Chain: "R" Number of atoms: 9416 Number of conformers: 1 Conformer: "" Number of residues, atoms: 596, 9416 Classifications: {'peptide': 596} Link IDs: {'PCIS': 3, 'PTRANS': 24, 'TRANS': 568} Chain breaks: 2 Unresolved non-hydrogen bonds: 34 Unresolved non-hydrogen angles: 41 Unresolved non-hydrogen dihedrals: 27 Unresolved non-hydrogen chiralities: 5 Planarities with less than four sites: {'GLU:plan': 1, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 9 Chain: "X" Number of atoms: 3286 Number of conformers: 1 Conformer: "" Number of residues, atoms: 221, 3286 Classifications: {'peptide': 221} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 5, 'TRANS': 215} Chain breaks: 2 Unresolved non-hydrogen bonds: 152 Unresolved non-hydrogen angles: 182 Unresolved non-hydrogen dihedrals: 122 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'GLN:plan1': 2, 'GLU:plan': 8, 'ASN:plan1': 4, 'ARG:plan': 4, 'ASP:plan': 10} Unresolved non-hydrogen planarities: 102 Chain: "Y" Number of atoms: 5110 Number of conformers: 1 Conformer: "" Number of residues, atoms: 338, 5110 Classifications: {'peptide': 338} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 5, 'TRANS': 332} Chain: "Z" Number of atoms: 871 Number of conformers: 1 Conformer: "" Number of residues, atoms: 56, 871 Classifications: {'peptide': 56} Link IDs: {'PTRANS': 4, 'TRANS': 51} Chain: "C" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen chiralities: 2 Chain: "A" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen chiralities: 1 Chain: "B" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen chiralities: 1 Chain: "R" Number of atoms: 56 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 56 Unusual residues: {'NAG': 4} Classifications: {'undetermined': 4} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen chiralities: 4 Time building chain proxies: 10.72, per 1000 atoms: 0.46 Number of scatterers: 23529 At special positions: 0 Unit cell: (103.272, 99.962, 202.572, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 100 16.00 O 2224 8.00 N 2028 7.00 C 7620 6.00 H 11557 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=16, symmetry=0 Simple disulfide: pdb=" SG CYS A 7 " - pdb=" SG CYS A 31 " distance=2.03 Simple disulfide: pdb=" SG CYS A 10 " - pdb=" SG CYS A 60 " distance=2.04 Simple disulfide: pdb=" SG CYS A 28 " - pdb=" SG CYS A 82 " distance=2.03 Simple disulfide: pdb=" SG CYS A 32 " - pdb=" SG CYS A 84 " distance=2.03 Simple disulfide: pdb=" SG CYS A 59 " - pdb=" SG CYS A 87 " distance=2.02 Simple disulfide: pdb=" SG CYS B 2 " - pdb=" SG CYS B 52 " distance=2.03 Simple disulfide: pdb=" SG CYS B 16 " - pdb=" SG CYS B 67 " distance=2.03 Simple disulfide: pdb=" SG CYS B 19 " - pdb=" SG CYS B 105 " distance=2.04 Simple disulfide: pdb=" SG CYS B 27 " - pdb=" SG CYS B 83 " distance=2.03 Simple disulfide: pdb=" SG CYS B 31 " - pdb=" SG CYS B 85 " distance=2.04 Simple disulfide: pdb=" SG CYS B 88 " - pdb=" SG CYS B 95 " distance=2.04 Simple disulfide: pdb=" SG CYS R 24 " - pdb=" SG CYS R 31 " distance=2.03 Simple disulfide: pdb=" SG CYS R 29 " - pdb=" SG CYS R 41 " distance=2.03 Simple disulfide: pdb=" SG CYS R 283 " - pdb=" SG CYS R 398 " distance=2.03 Simple disulfide: pdb=" SG CYS R 301 " - pdb=" SG CYS R 390 " distance=2.04 Simple disulfide: pdb=" SG CYS R 494 " - pdb=" SG CYS R 569 " distance=2.04 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Links applied BETA1-4 " NAG C 1 " - " NAG C 2 " NAG-ASN " NAG A 101 " - " ASN A 78 " " NAG B 201 " - " ASN B 23 " " NAG C 1 " - " ASN A 52 " " NAG R 801 " - " ASN R 99 " " NAG R 802 " - " ASN R 77 " " NAG R 803 " - " ASN R 198 " " NAG R 804 " - " ASN R 177 " Time building additional restraints: 18.26 Conformation dependent library (CDL) restraints added in 1.9 seconds 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2902 Finding SS restraints... Secondary structure from input PDB file: 35 helices and 19 sheets defined 32.2% alpha, 22.7% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.74 Creating SS restraints... Processing helix chain 'A' and resid 40 through 45 Processing helix chain 'B' and resid 43 through 48 removed outlier: 3.625A pdb=" N SER B 48 " --> pdb=" O TYR B 45 " (cutoff:3.500A) Processing helix chain 'N' and resid 28 through 32 removed outlier: 3.688A pdb=" N TYR N 32 " --> pdb=" O PHE N 29 " (cutoff:3.500A) Processing helix chain 'N' and resid 62 through 65 removed outlier: 3.507A pdb=" N LYS N 65 " --> pdb=" O GLY N 62 " (cutoff:3.500A) No H-bonds generated for 'chain 'N' and resid 62 through 65' Processing helix chain 'N' and resid 87 through 91 Processing helix chain 'R' and resid 172 through 176 removed outlier: 3.791A pdb=" N LEU R 175 " --> pdb=" O PHE R 172 " (cutoff:3.500A) removed outlier: 4.317A pdb=" N CYS R 176 " --> pdb=" O GLN R 173 " (cutoff:3.500A) No H-bonds generated for 'chain 'R' and resid 172 through 176' Processing helix chain 'R' and resid 218 through 222 removed outlier: 4.464A pdb=" N PHE R 221 " --> pdb=" O LYS R 218 " (cutoff:3.500A) Processing helix chain 'R' and resid 280 through 287 removed outlier: 3.622A pdb=" N LYS R 287 " --> pdb=" O CYS R 284 " (cutoff:3.500A) Processing helix chain 'R' and resid 294 through 304 Processing helix chain 'R' and resid 381 through 391 Processing helix chain 'R' and resid 416 through 442 removed outlier: 4.834A pdb=" N TRP R 422 " --> pdb=" O ARG R 418 " (cutoff:3.500A) Processing helix chain 'R' and resid 447 through 477 Processing helix chain 'R' and resid 480 through 491 removed outlier: 3.731A pdb=" N HIS R 484 " --> pdb=" O GLU R 480 " (cutoff:3.500A) removed outlier: 5.010A pdb=" N ILE R 486 " --> pdb=" O TYR R 482 " (cutoff:3.500A) removed outlier: 5.454A pdb=" N ASP R 487 " --> pdb=" O ASN R 483 " (cutoff:3.500A) Processing helix chain 'R' and resid 491 through 525 removed outlier: 4.163A pdb=" N ASN R 495 " --> pdb=" O GLY R 491 " (cutoff:3.500A) Processing helix chain 'R' and resid 528 through 532 Processing helix chain 'R' and resid 534 through 559 removed outlier: 3.572A pdb=" N ALA R 538 " --> pdb=" O ARG R 534 " (cutoff:3.500A) Proline residue: R 556 - end of helix removed outlier: 3.511A pdb=" N GLY R 559 " --> pdb=" O LEU R 555 " (cutoff:3.500A) Processing helix chain 'R' and resid 576 through 609 Processing helix chain 'R' and resid 620 through 649 removed outlier: 3.902A pdb=" N LEU R 629 " --> pdb=" O ARG R 625 " (cutoff:3.500A) Proline residue: R 639 - end of helix Processing helix chain 'R' and resid 655 through 666 removed outlier: 3.776A pdb=" N LYS R 660 " --> pdb=" O VAL R 656 " (cutoff:3.500A) removed outlier: 4.657A pdb=" N ILE R 661 " --> pdb=" O SER R 657 " (cutoff:3.500A) removed outlier: 3.943A pdb=" N LEU R 662 " --> pdb=" O ASN R 658 " (cutoff:3.500A) Processing helix chain 'R' and resid 668 through 679 Proline residue: R 675 - end of helix Processing helix chain 'R' and resid 682 through 697 Processing helix chain 'X' and resid 15 through 39 Processing helix chain 'X' and resid 52 through 59 Processing helix chain 'X' and resid 233 through 239 removed outlier: 3.521A pdb=" N ASN X 239 " --> pdb=" O GLN X 236 " (cutoff:3.500A) Processing helix chain 'X' and resid 265 through 279 Processing helix chain 'X' and resid 293 through 304 removed outlier: 3.528A pdb=" N GLY X 304 " --> pdb=" O LYS X 300 " (cutoff:3.500A) Processing helix chain 'X' and resid 312 through 317 Processing helix chain 'X' and resid 331 through 352 removed outlier: 3.947A pdb=" N THR X 335 " --> pdb=" O ASP X 331 " (cutoff:3.500A) removed outlier: 4.180A pdb=" N SER X 352 " --> pdb=" O ILE X 348 " (cutoff:3.500A) Processing helix chain 'X' and resid 370 through 391 Processing helix chain 'Y' and resid 4 through 25 Processing helix chain 'Y' and resid 29 through 34 Processing helix chain 'Z' and resid 8 through 24 removed outlier: 3.541A pdb=" N ASN Z 24 " --> pdb=" O LYS Z 20 " (cutoff:3.500A) Processing helix chain 'Z' and resid 29 through 45 Processing helix chain 'Z' and resid 46 through 48 No H-bonds generated for 'chain 'Z' and resid 46 through 48' Processing helix chain 'Z' and resid 55 through 59 removed outlier: 3.670A pdb=" N GLU Z 58 " --> pdb=" O PRO Z 55 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 11 through 14 Processing sheet with id=AA2, first strand: chain 'A' and resid 53 through 57 removed outlier: 4.160A pdb=" N CYS A 32 " --> pdb=" O CYS B 31 " (cutoff:3.500A) removed outlier: 4.048A pdb=" N CYS B 31 " --> pdb=" O CYS A 32 " (cutoff:3.500A) removed outlier: 10.327A pdb=" N THR B 25 " --> pdb=" O PRO A 38 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 59 through 67 removed outlier: 6.129A pdb=" N ALA A 81 " --> pdb=" O LYS A 63 " (cutoff:3.500A) removed outlier: 6.115A pdb=" N TYR A 65 " --> pdb=" O HIS A 79 " (cutoff:3.500A) removed outlier: 6.308A pdb=" N HIS A 79 " --> pdb=" O TYR A 65 " (cutoff:3.500A) removed outlier: 6.533A pdb=" N ARG A 67 " --> pdb=" O GLU A 77 " (cutoff:3.500A) removed outlier: 6.905A pdb=" N GLU A 77 " --> pdb=" O ARG A 67 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'B' and resid 51 through 63 removed outlier: 6.759A pdb=" N SER B 82 " --> pdb=" O ARG B 55 " (cutoff:3.500A) removed outlier: 6.732A pdb=" N PHE B 57 " --> pdb=" O ALA B 80 " (cutoff:3.500A) removed outlier: 6.934A pdb=" N ALA B 80 " --> pdb=" O PHE B 57 " (cutoff:3.500A) removed outlier: 6.467A pdb=" N TYR B 59 " --> pdb=" O PRO B 78 " (cutoff:3.500A) removed outlier: 6.599A pdb=" N THR B 61 " --> pdb=" O SER B 76 " (cutoff:3.500A) removed outlier: 6.461A pdb=" N SER B 76 " --> pdb=" O THR B 61 " (cutoff:3.500A) removed outlier: 6.796A pdb=" N GLU B 63 " --> pdb=" O TYR B 74 " (cutoff:3.500A) removed outlier: 7.225A pdb=" N TYR B 74 " --> pdb=" O GLU B 63 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'N' and resid 3 through 7 removed outlier: 3.535A pdb=" N GLN N 5 " --> pdb=" O ALA N 23 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'N' and resid 58 through 60 removed outlier: 3.549A pdb=" N SER N 59 " --> pdb=" O ASP N 50 " (cutoff:3.500A) removed outlier: 6.391A pdb=" N TRP N 36 " --> pdb=" O VAL N 48 " (cutoff:3.500A) removed outlier: 4.516A pdb=" N ASP N 50 " --> pdb=" O MET N 34 " (cutoff:3.500A) removed outlier: 6.752A pdb=" N MET N 34 " --> pdb=" O ASP N 50 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'R' and resid 30 through 33 removed outlier: 3.648A pdb=" N ARG R 38 " --> pdb=" O HIS R 32 " (cutoff:3.500A) removed outlier: 3.581A pdb=" N THR R 56 " --> pdb=" O PHE R 37 " (cutoff:3.500A) removed outlier: 6.662A pdb=" N LEU R 57 " --> pdb=" O TYR R 82 " (cutoff:3.500A) removed outlier: 8.222A pdb=" N SER R 84 " --> pdb=" O LEU R 57 " (cutoff:3.500A) removed outlier: 6.344A pdb=" N LEU R 59 " --> pdb=" O SER R 84 " (cutoff:3.500A) removed outlier: 6.721A pdb=" N ILE R 81 " --> pdb=" O GLU R 107 " (cutoff:3.500A) removed outlier: 8.202A pdb=" N ARG R 109 " --> pdb=" O ILE R 81 " (cutoff:3.500A) removed outlier: 6.266A pdb=" N VAL R 83 " --> pdb=" O ARG R 109 " (cutoff:3.500A) removed outlier: 6.536A pdb=" N THR R 111 " --> pdb=" O VAL R 83 " (cutoff:3.500A) removed outlier: 6.438A pdb=" N ILE R 106 " --> pdb=" O GLY R 132 " (cutoff:3.500A) removed outlier: 7.671A pdb=" N PHE R 134 " --> pdb=" O ILE R 106 " (cutoff:3.500A) removed outlier: 6.017A pdb=" N ILE R 108 " --> pdb=" O PHE R 134 " (cutoff:3.500A) removed outlier: 6.167A pdb=" N THR R 136 " --> pdb=" O ILE R 108 " (cutoff:3.500A) removed outlier: 7.988A pdb=" N ASN R 110 " --> pdb=" O THR R 136 " (cutoff:3.500A) removed outlier: 6.662A pdb=" N LEU R 131 " --> pdb=" O GLU R 157 " (cutoff:3.500A) removed outlier: 7.918A pdb=" N THR R 159 " --> pdb=" O LEU R 131 " (cutoff:3.500A) removed outlier: 6.520A pdb=" N ILE R 133 " --> pdb=" O THR R 159 " (cutoff:3.500A) removed outlier: 5.479A pdb=" N LEU R 180 " --> pdb=" O LEU R 202 " (cutoff:3.500A) removed outlier: 4.001A pdb=" N ALA R 204 " --> pdb=" O LEU R 180 " (cutoff:3.500A) removed outlier: 6.720A pdb=" N ASP R 203 " --> pdb=" O LEU R 230 " (cutoff:3.500A) removed outlier: 8.913A pdb=" N ASP R 232 " --> pdb=" O ASP R 203 " (cutoff:3.500A) removed outlier: 6.930A pdb=" N VAL R 205 " --> pdb=" O ASP R 232 " (cutoff:3.500A) removed outlier: 6.713A pdb=" N LEU R 231 " --> pdb=" O ILE R 253 " (cutoff:3.500A) removed outlier: 6.814A pdb=" N LEU R 252 " --> pdb=" O ASP R 276 " (cutoff:3.500A) removed outlier: 6.821A pdb=" N ALA R 275 " --> pdb=" O THR R 399 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'R' and resid 66 through 67 Processing sheet with id=AA9, first strand: chain 'R' and resid 97 through 98 removed outlier: 6.320A pdb=" N PHE R 97 " --> pdb=" O LYS R 123 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA9 Processing sheet with id=AB1, first strand: chain 'R' and resid 166 through 167 removed outlier: 3.573A pdb=" N ILE R 167 " --> pdb=" O SER R 191 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'X' and resid 210 through 214 removed outlier: 10.205A pdb=" N ALA X 243 " --> pdb=" O THR X 40 " (cutoff:3.500A) removed outlier: 7.508A pdb=" N ARG X 42 " --> pdb=" O ALA X 243 " (cutoff:3.500A) removed outlier: 8.266A pdb=" N ILE X 245 " --> pdb=" O ARG X 42 " (cutoff:3.500A) removed outlier: 6.536A pdb=" N LEU X 44 " --> pdb=" O ILE X 245 " (cutoff:3.500A) removed outlier: 7.377A pdb=" N VAL X 247 " --> pdb=" O LEU X 44 " (cutoff:3.500A) removed outlier: 6.559A pdb=" N LEU X 46 " --> pdb=" O VAL X 247 " (cutoff:3.500A) removed outlier: 8.069A pdb=" N ASP X 249 " --> pdb=" O LEU X 46 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'Y' and resid 46 through 51 removed outlier: 6.644A pdb=" N ILE Y 338 " --> pdb=" O ARG Y 48 " (cutoff:3.500A) removed outlier: 4.537A pdb=" N THR Y 50 " --> pdb=" O LEU Y 336 " (cutoff:3.500A) removed outlier: 6.535A pdb=" N LEU Y 336 " --> pdb=" O THR Y 50 " (cutoff:3.500A) removed outlier: 3.859A pdb=" N CYS Y 317 " --> pdb=" O GLY Y 330 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'Y' and resid 58 through 63 removed outlier: 3.905A pdb=" N ALA Y 60 " --> pdb=" O ALA Y 73 " (cutoff:3.500A) removed outlier: 3.561A pdb=" N LEU Y 70 " --> pdb=" O TRP Y 82 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'Y' and resid 58 through 63 removed outlier: 3.905A pdb=" N ALA Y 60 " --> pdb=" O ALA Y 73 " (cutoff:3.500A) removed outlier: 3.561A pdb=" N LEU Y 70 " --> pdb=" O TRP Y 82 " (cutoff:3.500A) removed outlier: 3.574A pdb=" N LEU Y 79 " --> pdb=" O ILE Y 93 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'Y' and resid 100 through 105 removed outlier: 7.027A pdb=" N GLY Y 115 " --> pdb=" O MET Y 101 " (cutoff:3.500A) removed outlier: 4.511A pdb=" N CYS Y 103 " --> pdb=" O ALA Y 113 " (cutoff:3.500A) removed outlier: 6.554A pdb=" N ALA Y 113 " --> pdb=" O CYS Y 103 " (cutoff:3.500A) removed outlier: 5.186A pdb=" N TYR Y 105 " --> pdb=" O TYR Y 111 " (cutoff:3.500A) removed outlier: 7.494A pdb=" N TYR Y 111 " --> pdb=" O TYR Y 105 " (cutoff:3.500A) removed outlier: 3.969A pdb=" N GLY Y 116 " --> pdb=" O ILE Y 120 " (cutoff:3.500A) removed outlier: 4.667A pdb=" N ILE Y 120 " --> pdb=" O GLY Y 116 " (cutoff:3.500A) removed outlier: 3.972A pdb=" N LEU Y 139 " --> pdb=" O CYS Y 121 " (cutoff:3.500A) removed outlier: 6.032A pdb=" N ASN Y 125 " --> pdb=" O VAL Y 135 " (cutoff:3.500A) removed outlier: 6.347A pdb=" N VAL Y 135 " --> pdb=" O ASN Y 125 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'Y' and resid 146 through 151 removed outlier: 3.724A pdb=" N CYS Y 148 " --> pdb=" O SER Y 160 " (cutoff:3.500A) removed outlier: 6.484A pdb=" N ASP Y 170 " --> pdb=" O GLN Y 176 " (cutoff:3.500A) removed outlier: 6.106A pdb=" N GLN Y 176 " --> pdb=" O ASP Y 170 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'Y' and resid 187 through 192 removed outlier: 7.047A pdb=" N GLY Y 202 " --> pdb=" O MET Y 188 " (cutoff:3.500A) removed outlier: 4.440A pdb=" N LEU Y 190 " --> pdb=" O VAL Y 200 " (cutoff:3.500A) removed outlier: 6.791A pdb=" N VAL Y 200 " --> pdb=" O LEU Y 190 " (cutoff:3.500A) removed outlier: 5.034A pdb=" N LEU Y 192 " --> pdb=" O LEU Y 198 " (cutoff:3.500A) removed outlier: 7.370A pdb=" N LEU Y 198 " --> pdb=" O LEU Y 192 " (cutoff:3.500A) removed outlier: 3.625A pdb=" N ALA Y 203 " --> pdb=" O SER Y 207 " (cutoff:3.500A) removed outlier: 4.523A pdb=" N SER Y 207 " --> pdb=" O ALA Y 203 " (cutoff:3.500A) removed outlier: 6.396A pdb=" N ASP Y 212 " --> pdb=" O CYS Y 218 " (cutoff:3.500A) removed outlier: 6.052A pdb=" N CYS Y 218 " --> pdb=" O ASP Y 212 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'Y' and resid 231 through 234 removed outlier: 3.735A pdb=" N ALA Y 231 " --> pdb=" O GLY Y 244 " (cutoff:3.500A) removed outlier: 6.479A pdb=" N CYS Y 250 " --> pdb=" O THR Y 263 " (cutoff:3.500A) removed outlier: 4.484A pdb=" N THR Y 263 " --> pdb=" O CYS Y 250 " (cutoff:3.500A) removed outlier: 6.612A pdb=" N LEU Y 252 " --> pdb=" O LEU Y 261 " (cutoff:3.500A) removed outlier: 3.500A pdb=" N GLN Y 259 " --> pdb=" O ASP Y 254 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'Y' and resid 273 through 278 removed outlier: 3.597A pdb=" N SER Y 275 " --> pdb=" O GLY Y 288 " (cutoff:3.500A) removed outlier: 5.953A pdb=" N ASP Y 298 " --> pdb=" O ARG Y 304 " (cutoff:3.500A) removed outlier: 6.118A pdb=" N ARG Y 304 " --> pdb=" O ASP Y 298 " (cutoff:3.500A) 506 hydrogen bonds defined for protein. 1398 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 9.48 Time building geometry restraints manager: 18.90 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.82 - 1.02: 11523 1.02 - 1.23: 57 1.23 - 1.43: 5136 1.43 - 1.63: 6939 1.63 - 1.83: 132 Bond restraints: 23787 Sorted by residual: bond pdb=" ND2 ASN A 15 " pdb="HD21 ASN A 15 " ideal model delta sigma weight residual 0.860 0.981 -0.121 2.00e-02 2.50e+03 3.63e+01 bond pdb=" ND2 ASN X 23 " pdb="HD22 ASN X 23 " ideal model delta sigma weight residual 0.860 0.980 -0.120 2.00e-02 2.50e+03 3.62e+01 bond pdb=" ND2 ASN X 23 " pdb="HD21 ASN X 23 " ideal model delta sigma weight residual 0.860 0.980 -0.120 2.00e-02 2.50e+03 3.61e+01 bond pdb=" ND2 ASN A 15 " pdb="HD22 ASN A 15 " ideal model delta sigma weight residual 0.860 0.980 -0.120 2.00e-02 2.50e+03 3.58e+01 bond pdb=" ND2 ASN A 52 " pdb="HD21 ASN A 52 " ideal model delta sigma weight residual 0.860 0.971 -0.111 2.00e-02 2.50e+03 3.11e+01 ... (remaining 23782 not shown) Histogram of bond angle deviations from ideal: 97.69 - 104.99: 207 104.99 - 112.29: 26666 112.29 - 119.59: 6478 119.59 - 126.89: 9263 126.89 - 134.19: 152 Bond angle restraints: 42766 Sorted by residual: angle pdb=" C ASN A 15 " pdb=" CA ASN A 15 " pdb=" CB ASN A 15 " ideal model delta sigma weight residual 108.86 99.99 8.87 1.44e+00 4.82e-01 3.79e+01 angle pdb=" N PHE A 74 " pdb=" CA PHE A 74 " pdb=" C PHE A 74 " ideal model delta sigma weight residual 110.41 117.56 -7.15 1.23e+00 6.61e-01 3.38e+01 angle pdb=" CA GLU X 21 " pdb=" C GLU X 21 " pdb=" N ALA X 22 " ideal model delta sigma weight residual 117.07 110.47 6.60 1.14e+00 7.69e-01 3.35e+01 angle pdb=" C PHE A 74 " pdb=" N LYS A 75 " pdb=" CA LYS A 75 " ideal model delta sigma weight residual 122.84 130.14 -7.30 1.30e+00 5.92e-01 3.16e+01 angle pdb=" C LYS X 24 " pdb=" N LYS X 25 " pdb=" CA LYS X 25 " ideal model delta sigma weight residual 120.28 127.35 -7.07 1.34e+00 5.57e-01 2.78e+01 ... (remaining 42761 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.84: 8816 17.84 - 35.68: 722 35.68 - 53.52: 114 53.52 - 71.37: 37 71.37 - 89.21: 27 Dihedral angle restraints: 9716 sinusoidal: 4378 harmonic: 5338 Sorted by residual: dihedral pdb=" CB CYS R 494 " pdb=" SG CYS R 494 " pdb=" SG CYS R 569 " pdb=" CB CYS R 569 " ideal model delta sinusoidal sigma weight residual -86.00 -35.08 -50.92 1 1.00e+01 1.00e-02 3.55e+01 dihedral pdb=" CB CYS R 301 " pdb=" SG CYS R 301 " pdb=" SG CYS R 390 " pdb=" CB CYS R 390 " ideal model delta sinusoidal sigma weight residual -86.00 -130.41 44.41 1 1.00e+01 1.00e-02 2.74e+01 dihedral pdb=" CA THR A 86 " pdb=" C THR A 86 " pdb=" N CYS A 87 " pdb=" CA CYS A 87 " ideal model delta harmonic sigma weight residual 180.00 155.09 24.91 0 5.00e+00 4.00e-02 2.48e+01 ... (remaining 9713 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.091: 1769 0.091 - 0.182: 157 0.182 - 0.273: 6 0.273 - 0.365: 1 0.365 - 0.456: 1 Chirality restraints: 1934 Sorted by residual: chirality pdb=" C1 NAG C 2 " pdb=" O4 NAG C 1 " pdb=" C2 NAG C 2 " pdb=" O5 NAG C 2 " both_signs ideal model delta sigma weight residual False -2.40 -2.26 -0.14 2.00e-02 2.50e+03 4.90e+01 chirality pdb=" C1 NAG A 101 " pdb=" ND2 ASN A 78 " pdb=" C2 NAG A 101 " pdb=" O5 NAG A 101 " both_signs ideal model delta sigma weight residual False -2.40 -2.86 0.46 2.00e-01 2.50e+01 5.19e+00 chirality pdb=" C1 NAG C 1 " pdb=" ND2 ASN A 52 " pdb=" C2 NAG C 1 " pdb=" O5 NAG C 1 " both_signs ideal model delta sigma weight residual False -2.40 -2.72 0.32 2.00e-01 2.50e+01 2.54e+00 ... (remaining 1931 not shown) Planarity restraints: 3566 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB ASN A 15 " 0.265 2.00e-02 2.50e+03 4.54e-01 3.09e+03 pdb=" CG ASN A 15 " -0.021 2.00e-02 2.50e+03 pdb=" OD1 ASN A 15 " -0.258 2.00e-02 2.50e+03 pdb=" ND2 ASN A 15 " -0.001 2.00e-02 2.50e+03 pdb="HD21 ASN A 15 " 0.749 2.00e-02 2.50e+03 pdb="HD22 ASN A 15 " -0.734 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB ASN X 23 " 0.292 2.00e-02 2.50e+03 3.04e-01 1.38e+03 pdb=" CG ASN X 23 " -0.050 2.00e-02 2.50e+03 pdb=" OD1 ASN X 23 " -0.261 2.00e-02 2.50e+03 pdb=" ND2 ASN X 23 " -0.005 2.00e-02 2.50e+03 pdb="HD21 ASN X 23 " -0.434 2.00e-02 2.50e+03 pdb="HD22 ASN X 23 " 0.457 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C2 NAG C 1 " 0.305 2.00e-02 2.50e+03 2.65e-01 8.81e+02 pdb=" C7 NAG C 1 " -0.065 2.00e-02 2.50e+03 pdb=" C8 NAG C 1 " 0.027 2.00e-02 2.50e+03 pdb=" N2 NAG C 1 " -0.464 2.00e-02 2.50e+03 pdb=" O7 NAG C 1 " 0.198 2.00e-02 2.50e+03 ... (remaining 3563 not shown) Histogram of nonbonded interaction distances: 1.67 - 2.26: 2688 2.26 - 2.84: 51425 2.84 - 3.43: 59957 3.43 - 4.01: 82640 4.01 - 4.60: 125495 Nonbonded interactions: 322205 Sorted by model distance: nonbonded pdb=" O ALA R 380 " pdb=" HG SER R 383 " model vdw 1.670 1.850 nonbonded pdb=" H LYS R 291 " pdb=" O GLU R 394 " model vdw 1.694 1.850 nonbonded pdb=" HZ1 LYS X 233 " pdb=" OD2 ASP Y 228 " model vdw 1.696 1.850 nonbonded pdb=" O LYS R 250 " pdb=" H THR R 273 " model vdw 1.702 1.850 nonbonded pdb=" O ILE B 36 " pdb=" HZ2 LYS B 39 " model vdw 1.709 1.850 ... (remaining 322200 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians S 100 5.16 5 C 7620 2.51 5 N 2028 2.21 5 O 2224 1.98 5 H 11557 0.53 5 sf(0) = scattering factor at diffraction angle 0. Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 8.640 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.350 Extract box with map and model: 29.310 Check model and map are aligned: 0.350 Process input model: 72.190 Find NCS groups from input model: 0.420 Set up NCS constraints: 0.080 Set refine NCS operators: 0.000 Set scattering table: 0.210 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:13.880 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 125.430 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6263 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.058 12230 Z= 0.342 Angle : 0.909 9.661 16612 Z= 0.526 Chirality : 0.052 0.456 1934 Planarity : 0.010 0.265 2096 Dihedral : 12.717 84.578 4362 Min Nonbonded Distance : 2.016 Molprobity Statistics. All-atom Clashscore : 3.30 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.64 % Favored : 96.36 % Rotamer Outliers : 0.62 % Cbeta Deviations : 0.07 % Peptide Plane: Cis-proline : 5.17 % Cis-general : 0.07 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.16 (0.22), residues: 1510 helix: 1.51 (0.25), residues: 428 sheet: -1.14 (0.29), residues: 305 loop : -0.52 (0.22), residues: 777 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Residue LYS 45 is missing expected H atoms. Skipping. Residue LYS 75 is missing expected H atoms. Skipping. Residue LEU 40 is missing expected H atoms. Skipping. Residue LYS 101 is missing expected H atoms. Skipping. Residue THR 21 is missing expected H atoms. Skipping. Residue LYS 262 is missing expected H atoms. Skipping. Residue ILE 295 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LEU 299 is missing expected H atoms. Skipping. Residue MET 300 is missing expected H atoms. Skipping. Residue SER 305 is missing expected H atoms. Skipping. Residue LYS 621 is missing expected H atoms. Skipping. Residue LYS 17 is missing expected H atoms. Skipping. Residue LYS 24 is missing expected H atoms. Skipping. Residue LYS 25 is missing expected H atoms. Skipping. Residue LYS 28 is missing expected H atoms. Skipping. Residue LYS 32 is missing expected H atoms. Skipping. Residue SER 54 is missing expected H atoms. Skipping. Residue THR 55 is missing expected H atoms. Skipping. Residue LYS 58 is missing expected H atoms. Skipping. Residue LYS 216 is missing expected H atoms. Skipping. Residue LYS 274 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LYS 305 is missing expected H atoms. Skipping. Residue LYS 307 is missing expected H atoms. Skipping. Evaluate side-chains 277 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 269 time to evaluate : 1.834 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 8 outliers final: 1 residues processed: 271 average time/residue: 2.0356 time to fit residues: 613.0358 Evaluate side-chains 172 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 171 time to evaluate : 1.718 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 1 residues processed: 0 time to fit residues: 2.2052 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 149 random chunks: chunk 125 optimal weight: 2.9990 chunk 112 optimal weight: 5.9990 chunk 62 optimal weight: 2.9990 chunk 38 optimal weight: 1.9990 chunk 76 optimal weight: 1.9990 chunk 60 optimal weight: 0.2980 chunk 116 optimal weight: 2.9990 chunk 45 optimal weight: 3.9990 chunk 71 optimal weight: 1.9990 chunk 86 optimal weight: 1.9990 chunk 135 optimal weight: 2.9990 overall best weight: 1.6588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 27 GLN ** A 79 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** N 123 GLN ** R 113 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** R 135 ASN R 173 GLN Y 13 GLN Y 340 ASN Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6439 moved from start: 0.2745 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.046 12230 Z= 0.301 Angle : 0.629 7.533 16612 Z= 0.343 Chirality : 0.043 0.206 1934 Planarity : 0.005 0.056 2096 Dihedral : 5.113 27.578 1711 Min Nonbonded Distance : 1.968 Molprobity Statistics. All-atom Clashscore : 7.95 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.70 % Favored : 95.30 % Rotamer Outliers : 2.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.17 % Cis-general : 0.07 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.14 (0.21), residues: 1510 helix: 1.81 (0.25), residues: 415 sheet: -1.19 (0.29), residues: 299 loop : -0.64 (0.22), residues: 796 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Residue LYS 45 is missing expected H atoms. Skipping. Residue LYS 75 is missing expected H atoms. Skipping. Residue LEU 40 is missing expected H atoms. Skipping. Residue LYS 101 is missing expected H atoms. Skipping. Residue THR 21 is missing expected H atoms. Skipping. Residue LYS 262 is missing expected H atoms. Skipping. Residue ILE 295 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LEU 299 is missing expected H atoms. Skipping. Residue MET 300 is missing expected H atoms. Skipping. Residue SER 305 is missing expected H atoms. Skipping. Residue LYS 621 is missing expected H atoms. Skipping. Residue LYS 17 is missing expected H atoms. Skipping. Residue LYS 24 is missing expected H atoms. Skipping. Residue LYS 25 is missing expected H atoms. Skipping. Residue LYS 28 is missing expected H atoms. Skipping. Residue LYS 32 is missing expected H atoms. Skipping. Residue SER 54 is missing expected H atoms. Skipping. Residue THR 55 is missing expected H atoms. Skipping. Residue LYS 58 is missing expected H atoms. Skipping. Residue LYS 216 is missing expected H atoms. Skipping. Residue LYS 274 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LYS 305 is missing expected H atoms. Skipping. Residue LYS 307 is missing expected H atoms. Skipping. Evaluate side-chains 199 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 172 time to evaluate : 2.066 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 27 outliers final: 14 residues processed: 191 average time/residue: 1.9186 time to fit residues: 410.8481 Evaluate side-chains 160 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 146 time to evaluate : 1.802 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 12 residues processed: 2 average time/residue: 0.9843 time to fit residues: 4.7922 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 149 random chunks: chunk 75 optimal weight: 2.9990 chunk 41 optimal weight: 2.9990 chunk 112 optimal weight: 5.9990 chunk 92 optimal weight: 2.9990 chunk 37 optimal weight: 5.9990 chunk 135 optimal weight: 0.8980 chunk 146 optimal weight: 3.9990 chunk 120 optimal weight: 0.9980 chunk 134 optimal weight: 1.9990 chunk 46 optimal weight: 2.9990 chunk 108 optimal weight: 3.9990 overall best weight: 1.9786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 27 GLN A 90 HIS ** R 113 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Y 9 GLN Y 225 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6561 moved from start: 0.3905 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.045 12230 Z= 0.355 Angle : 0.662 7.139 16612 Z= 0.361 Chirality : 0.045 0.276 1934 Planarity : 0.006 0.065 2096 Dihedral : 5.473 33.111 1711 Min Nonbonded Distance : 1.884 Molprobity Statistics. All-atom Clashscore : 9.18 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.70 % Favored : 94.30 % Rotamer Outliers : 3.04 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.17 % Cis-general : 0.07 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.69 (0.21), residues: 1510 helix: 1.48 (0.25), residues: 424 sheet: -1.35 (0.28), residues: 316 loop : -1.16 (0.21), residues: 770 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Residue LYS 45 is missing expected H atoms. Skipping. Residue LYS 75 is missing expected H atoms. Skipping. Residue LEU 40 is missing expected H atoms. Skipping. Residue LYS 101 is missing expected H atoms. Skipping. Residue THR 21 is missing expected H atoms. Skipping. Residue LYS 262 is missing expected H atoms. Skipping. Residue ILE 295 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LEU 299 is missing expected H atoms. Skipping. Residue MET 300 is missing expected H atoms. Skipping. Residue SER 305 is missing expected H atoms. Skipping. Residue LYS 621 is missing expected H atoms. Skipping. Residue LYS 17 is missing expected H atoms. Skipping. Residue LYS 24 is missing expected H atoms. Skipping. Residue LYS 25 is missing expected H atoms. Skipping. Residue LYS 28 is missing expected H atoms. Skipping. Residue LYS 32 is missing expected H atoms. Skipping. Residue SER 54 is missing expected H atoms. Skipping. Residue THR 55 is missing expected H atoms. Skipping. Residue LYS 58 is missing expected H atoms. Skipping. Residue LYS 216 is missing expected H atoms. Skipping. Residue LYS 274 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LYS 305 is missing expected H atoms. Skipping. Residue LYS 307 is missing expected H atoms. Skipping. Evaluate side-chains 223 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 39 poor density : 184 time to evaluate : 2.088 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 39 outliers final: 20 residues processed: 211 average time/residue: 1.7524 time to fit residues: 421.6089 Evaluate side-chains 185 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 165 time to evaluate : 1.965 Switching outliers to nearest non-outliers outliers start: 20 outliers final: 15 residues processed: 6 average time/residue: 1.4806 time to fit residues: 12.8298 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 149 random chunks: chunk 133 optimal weight: 1.9990 chunk 101 optimal weight: 1.9990 chunk 70 optimal weight: 1.9990 chunk 15 optimal weight: 0.0570 chunk 64 optimal weight: 3.9990 chunk 91 optimal weight: 2.9990 chunk 136 optimal weight: 2.9990 chunk 144 optimal weight: 0.9980 chunk 71 optimal weight: 0.8980 chunk 128 optimal weight: 1.9990 chunk 38 optimal weight: 1.9990 overall best weight: 1.1902 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: R 113 ASN R 173 GLN Y 13 GLN Y 225 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6546 moved from start: 0.4127 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.091 12230 Z= 0.246 Angle : 0.578 9.306 16612 Z= 0.311 Chirality : 0.042 0.227 1934 Planarity : 0.004 0.044 2096 Dihedral : 5.215 32.677 1711 Min Nonbonded Distance : 1.851 Molprobity Statistics. All-atom Clashscore : 8.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.64 % Favored : 95.36 % Rotamer Outliers : 2.58 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.17 % Cis-general : 0.07 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.54 (0.22), residues: 1510 helix: 1.71 (0.26), residues: 424 sheet: -1.38 (0.28), residues: 309 loop : -1.08 (0.21), residues: 777 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Residue LYS 45 is missing expected H atoms. Skipping. Residue LYS 75 is missing expected H atoms. Skipping. Residue LEU 40 is missing expected H atoms. Skipping. Residue LYS 101 is missing expected H atoms. Skipping. Residue THR 21 is missing expected H atoms. Skipping. Residue LYS 262 is missing expected H atoms. Skipping. Residue ILE 295 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LEU 299 is missing expected H atoms. Skipping. Residue MET 300 is missing expected H atoms. Skipping. Residue SER 305 is missing expected H atoms. Skipping. Residue LYS 621 is missing expected H atoms. Skipping. Residue LYS 17 is missing expected H atoms. Skipping. Residue LYS 24 is missing expected H atoms. Skipping. Residue LYS 25 is missing expected H atoms. Skipping. Residue LYS 28 is missing expected H atoms. Skipping. Residue LYS 32 is missing expected H atoms. Skipping. Residue SER 54 is missing expected H atoms. Skipping. Residue THR 55 is missing expected H atoms. Skipping. Residue LYS 58 is missing expected H atoms. Skipping. Residue LYS 216 is missing expected H atoms. Skipping. Residue LYS 274 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LYS 305 is missing expected H atoms. Skipping. Residue LYS 307 is missing expected H atoms. Skipping. Evaluate side-chains 202 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 33 poor density : 169 time to evaluate : 2.014 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 33 outliers final: 26 residues processed: 196 average time/residue: 1.8339 time to fit residues: 407.3563 Evaluate side-chains 183 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 26 poor density : 157 time to evaluate : 1.869 Switching outliers to nearest non-outliers outliers start: 26 outliers final: 22 residues processed: 5 average time/residue: 0.6772 time to fit residues: 7.0545 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 149 random chunks: chunk 119 optimal weight: 2.9990 chunk 81 optimal weight: 4.9990 chunk 2 optimal weight: 4.9990 chunk 107 optimal weight: 4.9990 chunk 59 optimal weight: 8.9990 chunk 122 optimal weight: 2.9990 chunk 99 optimal weight: 1.9990 chunk 0 optimal weight: 0.8980 chunk 73 optimal weight: 3.9990 chunk 129 optimal weight: 1.9990 chunk 36 optimal weight: 1.9990 overall best weight: 1.9788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: R 113 ASN Y 9 GLN Y 13 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6613 moved from start: 0.4658 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.048 12230 Z= 0.354 Angle : 0.647 7.182 16612 Z= 0.353 Chirality : 0.044 0.275 1934 Planarity : 0.005 0.058 2096 Dihedral : 5.499 36.623 1711 Min Nonbonded Distance : 1.859 Molprobity Statistics. All-atom Clashscore : 9.65 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.28 % Favored : 92.72 % Rotamer Outliers : 4.53 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.17 % Cis-general : 0.07 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.84 (0.21), residues: 1510 helix: 1.40 (0.25), residues: 429 sheet: -1.45 (0.29), residues: 304 loop : -1.30 (0.22), residues: 777 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3020 Ramachandran restraints generated. 1510 Oldfield, 0 Emsley, 1510 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Residue LYS 45 is missing expected H atoms. Skipping. Residue LYS 75 is missing expected H atoms. Skipping. Residue LEU 40 is missing expected H atoms. Skipping. Residue LYS 101 is missing expected H atoms. Skipping. Residue THR 21 is missing expected H atoms. Skipping. Residue LYS 262 is missing expected H atoms. Skipping. Residue ILE 295 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LEU 299 is missing expected H atoms. Skipping. Residue MET 300 is missing expected H atoms. Skipping. Residue SER 305 is missing expected H atoms. Skipping. Residue LYS 621 is missing expected H atoms. Skipping. Residue LYS 17 is missing expected H atoms. Skipping. Residue LYS 24 is missing expected H atoms. Skipping. Residue LYS 25 is missing expected H atoms. Skipping. Residue LYS 28 is missing expected H atoms. Skipping. Residue LYS 32 is missing expected H atoms. Skipping. Residue SER 54 is missing expected H atoms. Skipping. Residue THR 55 is missing expected H atoms. Skipping. Residue LYS 58 is missing expected H atoms. Skipping. Residue LYS 216 is missing expected H atoms. Skipping. Residue LYS 274 is missing expected H atoms. Skipping. Residue LEU 296 is missing expected H atoms. Skipping. Residue LYS 305 is missing expected H atoms. Skipping. Residue LYS 307 is missing expected H atoms. Skipping. Evaluate side-chains 223 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 58 poor density : 165 time to evaluate : 1.898 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 58 outliers final: 35 residues processed: 209 average time/residue: 1.8905 time to fit residues: 445.8000 Evaluate side-chains 190 residues out of total 1344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 35 poor density : 155 time to evaluate : 1.999 Switching outliers to nearest non-outliers outliers start: 35 outliers final: 31 residues processed: 5 average time/residue: 0.9970 time to fit residues: 8.5011 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/iotbx/cli_parser.py", line 865, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/phenix/phenix/programs/real_space_refine.py", line 192, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 297, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 312, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 673, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 736, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1479, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1381, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1262, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 67.816 > 50: distance: 66 - 75: 25.899 distance: 75 - 76: 35.811 distance: 75 - 81: 29.077 distance: 76 - 77: 17.453 distance: 76 - 79: 18.545 distance: 76 - 82: 22.922 distance: 77 - 78: 24.158 distance: 77 - 89: 43.030 distance: 79 - 80: 8.216 distance: 79 - 83: 10.588 distance: 79 - 84: 12.520 distance: 80 - 81: 22.701 distance: 80 - 85: 20.173 distance: 80 - 86: 15.927 distance: 81 - 87: 13.560 distance: 81 - 88: 28.985 distance: 89 - 90: 18.836 distance: 89 - 94: 14.925 distance: 90 - 93: 30.240 distance: 90 - 95: 21.460 distance: 91 - 92: 10.181 distance: 91 - 99: 18.734 distance: 93 - 96: 8.741 distance: 93 - 97: 8.763 distance: 93 - 98: 9.983 distance: 99 - 100: 38.395 distance: 99 - 105: 38.290 distance: 100 - 101: 29.551 distance: 100 - 103: 16.774 distance: 100 - 106: 19.067 distance: 101 - 102: 31.263 distance: 101 - 113: 33.838 distance: 103 - 104: 21.289 distance: 103 - 107: 19.762 distance: 103 - 108: 11.794 distance: 104 - 105: 12.401 distance: 104 - 109: 5.005 distance: 104 - 110: 5.408 distance: 105 - 111: 3.613 distance: 105 - 112: 30.796 distance: 113 - 114: 27.307 distance: 113 - 124: 36.967 distance: 114 - 115: 16.490 distance: 114 - 117: 36.349 distance: 114 - 125: 38.765 distance: 115 - 116: 29.139 distance: 115 - 133: 21.541 distance: 117 - 118: 4.621 distance: 117 - 126: 34.274 distance: 117 - 127: 8.579 distance: 118 - 119: 9.967 distance: 118 - 120: 14.684 distance: 119 - 121: 8.463 distance: 119 - 128: 13.401 distance: 120 - 122: 11.495 distance: 120 - 129: 7.864 distance: 121 - 123: 10.101 distance: 121 - 130: 6.878 distance: 122 - 123: 10.050 distance: 122 - 131: 6.671 distance: 123 - 132: 6.946 distance: 133 - 134: 4.103 distance: 133 - 140: 25.267 distance: 134 - 135: 11.257 distance: 134 - 137: 15.388 distance: 134 - 141: 9.236 distance: 135 - 136: 32.628 distance: 135 - 147: 13.822 distance: 137 - 138: 15.422 distance: 137 - 139: 6.839 distance: 137 - 142: 10.325 distance: 138 - 143: 7.627 distance: 139 - 144: 8.110 distance: 139 - 145: 16.050 distance: 139 - 146: 4.902 distance: 147 - 148: 16.113 distance: 147 - 158: 26.983 distance: 148 - 149: 8.703 distance: 148 - 151: 23.342 distance: 148 - 159: 34.719 distance: 149 - 150: 19.684 distance: 149 - 171: 21.234 distance: 151 - 152: 28.145 distance: 151 - 160: 27.280 distance: 151 - 161: 5.774 distance: 152 - 153: 7.228 distance: 152 - 162: 10.325 distance: 152 - 163: 26.119 distance: 153 - 154: 4.969 distance: 153 - 164: 15.549 distance: 153 - 165: 10.423 distance: 154 - 155: 3.855 distance: 154 - 166: 4.766 distance: 155 - 156: 3.868 distance: 155 - 157: 4.437 distance: 171 - 172: 23.561 distance: 171 - 179: 16.011 distance: 172 - 173: 10.967 distance: 172 - 175: 35.330 distance: 172 - 180: 39.701 distance: 173 - 174: 42.306 distance: 173 - 183: 50.151 distance: 175 - 176: 13.130 distance: 175 - 181: 18.355 distance: 175 - 182: 9.056 distance: 176 - 177: 29.785 distance: 176 - 178: 14.241 distance: 183 - 184: 37.735 distance: 183 - 189: 31.954 distance: 184 - 185: 11.660 distance: 184 - 187: 51.334 distance: 184 - 190: 18.141 distance: 185 - 194: 41.730 distance: 187 - 188: 11.684 distance: 187 - 191: 7.023 distance: 187 - 192: 3.100 distance: 188 - 193: 28.497 distance: 194 - 195: 13.517 distance: 194 - 205: 7.923 distance: 195 - 196: 15.520 distance: 195 - 198: 5.287 distance: 195 - 206: 13.175 distance: 196 - 197: 14.891 distance: 196 - 214: 38.644 distance: 198 - 199: 16.539 distance: 198 - 207: 24.964 distance: 198 - 208: 20.814 distance: 199 - 200: 9.240 distance: 199 - 201: 10.284 distance: 200 - 202: 5.812 distance: 200 - 209: 5.924 distance: 201 - 203: 8.708 distance: 201 - 210: 6.821 distance: 202 - 204: 3.559 distance: 204 - 213: 3.947 distance: 214 - 215: 6.713 distance: 214 - 222: 17.760 distance: 215 - 216: 14.934 distance: 215 - 218: 25.076 distance: 215 - 223: 5.614 distance: 216 - 217: 21.756 distance: 216 - 226: 3.794 distance: 218 - 219: 13.862 distance: 218 - 224: 19.484 distance: 218 - 225: 14.618 distance: 219 - 220: 13.408 distance: 219 - 221: 10.235