Starting phenix.real_space_refine on Fri Dec 8 09:04:27 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hpt_34943/12_2023/8hpt_34943_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hpt_34943/12_2023/8hpt_34943.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.39 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hpt_34943/12_2023/8hpt_34943.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hpt_34943/12_2023/8hpt_34943.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hpt_34943/12_2023/8hpt_34943_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hpt_34943/12_2023/8hpt_34943_updated.pdb" } resolution = 3.39 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } 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 Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- 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.000 sd= 0.024 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 Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 48 5.16 5 C 5051 2.51 5 N 1343 2.21 5 O 1455 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A PHE 75": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 297": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C GLU 226": "OE1" <-> "OE2" Residue "C PHE 253": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ASP 322": "OD1" <-> "OD2" Residue "C PHE 335": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H ASP 90": "OD1" <-> "OD2" Residue "H GLU 220": "OE1" <-> "OE2" Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 7897 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 1872 Number of conformers: 1 Conformer: "" Number of residues, atoms: 255, 1872 Classifications: {'peptide': 255} Incomplete info: {'backbone_only': 3, 'truncation_to_alanine': 29} Link IDs: {'PTRANS': 11, 'TRANS': 243} Chain breaks: 2 Unresolved non-hydrogen bonds: 137 Unresolved non-hydrogen angles: 180 Unresolved non-hydrogen dihedrals: 118 Unresolved non-hydrogen chiralities: 16 Planarities with less than four sites: {'TYR:plan': 2, 'ASN:plan1': 2, 'TRP:plan': 2, 'ASP:plan': 1, 'PHE:plan': 2, 'GLU:plan': 1, 'ARG:plan': 3} Unresolved non-hydrogen planarities: 73 Chain: "D" Number of atoms: 62 Number of conformers: 1 Conformer: "" Number of residues, atoms: 6, 62 Classifications: {'peptide': 6} Modifications used: {'COO': 1, 'PEPT-D': 1} Link IDs: {'PTRANS': 1, 'TRANS': 4} Unresolved non-hydrogen bonds: 3 Unresolved non-hydrogen angles: 6 Planarities with less than four sites: {'ZAL:plan-1': 1, 'ALC:plan-1': 1, 'MEA:plan-2': 1} Unresolved non-hydrogen planarities: 3 Chain: "B" Number of atoms: 1458 Number of conformers: 1 Conformer: "" Number of residues, atoms: 206, 1458 Classifications: {'peptide': 206} Modifications used: {'COO': 1} Incomplete info: {'backbone_only': 2, 'truncation_to_alanine': 48} Link IDs: {'PTRANS': 4, 'TRANS': 201} Chain breaks: 3 Unresolved non-hydrogen bonds: 182 Unresolved non-hydrogen angles: 220 Unresolved non-hydrogen dihedrals: 151 Unresolved non-hydrogen chiralities: 10 Planarities with less than four sites: {'GLN:plan1': 2, 'ASP:plan': 7, 'TYR:plan': 1, 'ASN:plan1': 3, 'HIS:plan': 2, 'GLU:plan': 7, 'ARG:plan': 2} Unresolved non-hydrogen planarities: 95 Chain: "C" Number of atoms: 2436 Number of conformers: 1 Conformer: "" Number of residues, atoms: 338, 2436 Classifications: {'peptide': 338} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 51} Link IDs: {'PTRANS': 5, 'TRANS': 332} Unresolved non-hydrogen bonds: 165 Unresolved non-hydrogen angles: 205 Unresolved non-hydrogen dihedrals: 125 Unresolved non-hydrogen chiralities: 13 Planarities with less than four sites: {'GLN:plan1': 5, 'ASN:plan1': 1, 'ASP:plan': 8, 'ASN%COO:plan1': 1, 'GLU:plan': 5, 'ARG:plan': 7} Unresolved non-hydrogen planarities: 105 Chain: "G" Number of atoms: 395 Number of conformers: 1 Conformer: "" Number of residues, atoms: 56, 395 Classifications: {'peptide': 56} Incomplete info: {'truncation_to_alanine': 9} Link IDs: {'PTRANS': 4, 'TRANS': 51} Unresolved non-hydrogen bonds: 34 Unresolved non-hydrogen angles: 41 Unresolved non-hydrogen dihedrals: 27 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'GLU:plan': 2, 'ARG:plan': 2, 'ASP:plan': 1} Unresolved non-hydrogen planarities: 21 Chain: "H" Number of atoms: 1674 Number of conformers: 1 Conformer: "" Number of residues, atoms: 234, 1674 Classifications: {'peptide': 234} Incomplete info: {'truncation_to_alanine': 38} Link IDs: {'PCIS': 1, 'PTRANS': 9, 'TRANS': 223} Chain breaks: 1 Unresolved non-hydrogen bonds: 126 Unresolved non-hydrogen angles: 150 Unresolved non-hydrogen dihedrals: 102 Unresolved non-hydrogen chiralities: 7 Planarities with less than four sites: {'GLN:plan1': 4, 'GLU:plan': 6, 'ARG:plan': 4, 'ASN:plan1': 1, 'ASP:plan': 2} Unresolved non-hydrogen planarities: 69 Time building chain proxies: 4.91, per 1000 atoms: 0.62 Number of scatterers: 7897 At special positions: 0 Unit cell: (89.0706, 126.025, 128.868, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 48 16.00 O 1455 8.00 N 1343 7.00 C 5051 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS H 22 " - pdb=" SG CYS H 96 " distance=2.04 Simple disulfide: pdb=" SG CYS H 159 " - pdb=" SG CYS H 229 " distance=2.03 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 Time building additional restraints: 3.46 Conformation dependent library (CDL) restraints added in 1.8 seconds 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Skipped Input residue name is d-peptide pdb=" CB DAR D 406 " Input volumes are d-peptide like pdb=" CB ZAL D 404 " Number of C-beta restraints generated: 2032 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 27 helices and 15 sheets defined 34.4% alpha, 21.1% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.63 Creating SS restraints... Processing helix chain 'A' and resid 37 through 64 Proline residue: A 53 - end of helix removed outlier: 4.507A pdb=" N PHE A 64 " --> pdb=" O TRP A 60 " (cutoff:3.500A) Processing helix chain 'A' and resid 70 through 97 removed outlier: 4.309A pdb=" N TRP A 74 " --> pdb=" O ASN A 71 " (cutoff:3.500A) removed outlier: 3.594A pdb=" N LEU A 87 " --> pdb=" O LEU A 84 " (cutoff:3.500A) Proline residue: A 90 - end of helix removed outlier: 3.550A pdb=" N THR A 94 " --> pdb=" O VAL A 91 " (cutoff:3.500A) removed outlier: 3.628A pdb=" N THR A 95 " --> pdb=" O LEU A 92 " (cutoff:3.500A) removed outlier: 3.932A pdb=" N VAL A 96 " --> pdb=" O PHE A 93 " (cutoff:3.500A) Processing helix chain 'A' and resid 107 through 109 No H-bonds generated for 'chain 'A' and resid 107 through 109' Processing helix chain 'A' and resid 112 through 139 Processing helix chain 'A' and resid 143 through 146 No H-bonds generated for 'chain 'A' and resid 143 through 146' Processing helix chain 'A' and resid 150 through 173 Proline residue: A 170 - end of helix Processing helix chain 'A' and resid 202 through 232 removed outlier: 5.591A pdb=" N LEU A 214 " --> pdb=" O VAL A 210 " (cutoff:3.500A) Proline residue: A 215 - end of helix removed outlier: 4.208A pdb=" N SER A 232 " --> pdb=" O LEU A 228 " (cutoff:3.500A) Processing helix chain 'A' and resid 234 through 236 No H-bonds generated for 'chain 'A' and resid 234 through 236' Processing helix chain 'A' and resid 240 through 265 Proline residue: A 258 - end of helix Processing helix chain 'A' and resid 277 through 304 removed outlier: 4.681A pdb=" N SER A 284 " --> pdb=" O GLU A 280 " (cutoff:3.500A) removed outlier: 4.865A pdb=" N LEU A 285 " --> pdb=" O LYS A 281 " (cutoff:3.500A) Proline residue: A 298 - end of helix removed outlier: 3.970A pdb=" N TYR A 301 " --> pdb=" O ASN A 297 " (cutoff:3.500A) removed outlier: 3.614A pdb=" N VAL A 302 " --> pdb=" O PRO A 298 " (cutoff:3.500A) Processing helix chain 'A' and resid 306 through 314 Processing helix chain 'B' and resid 9 through 31 Processing helix chain 'B' and resid 46 through 52 Processing helix chain 'B' and resid 209 through 216 removed outlier: 4.616A pdb=" N HIS B 214 " --> pdb=" O LYS B 210 " (cutoff:3.500A) removed outlier: 6.025A pdb=" N CYS B 215 " --> pdb=" O LYS B 211 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N PHE B 216 " --> pdb=" O TRP B 212 " (cutoff:3.500A) Processing helix chain 'B' and resid 246 through 255 Processing helix chain 'B' and resid 272 through 281 Processing helix chain 'B' and resid 297 through 309 Processing helix chain 'B' and resid 331 through 350 Processing helix chain 'C' and resid 4 through 24 Processing helix chain 'C' and resid 30 through 35 removed outlier: 3.747A pdb=" N ASN C 35 " --> pdb=" O SER C 31 " (cutoff:3.500A) Processing helix chain 'G' and resid 8 through 24 removed outlier: 3.523A pdb=" N GLN G 11 " --> pdb=" O SER G 8 " (cutoff:3.500A) removed outlier: 3.569A pdb=" N ASN G 24 " --> pdb=" O MET G 21 " (cutoff:3.500A) Processing helix chain 'G' and resid 30 through 43 Processing helix chain 'G' and resid 45 through 47 No H-bonds generated for 'chain 'G' and resid 45 through 47' Processing helix chain 'H' and resid 29 through 31 No H-bonds generated for 'chain 'H' and resid 29 through 31' Processing helix chain 'H' and resid 62 through 64 No H-bonds generated for 'chain 'H' and resid 62 through 64' Processing helix chain 'H' and resid 88 through 90 No H-bonds generated for 'chain 'H' and resid 88 through 90' Processing helix chain 'H' and resid 221 through 223 No H-bonds generated for 'chain 'H' and resid 221 through 223' Processing sheet with id= A, first strand: chain 'B' and resid 319 through 323 removed outlier: 6.805A pdb=" N ILE B 265 " --> pdb=" O TYR B 320 " (cutoff:3.500A) removed outlier: 8.359A pdb=" N HIS B 322 " --> pdb=" O ILE B 265 " (cutoff:3.500A) removed outlier: 6.653A pdb=" N LEU B 267 " --> pdb=" O HIS B 322 " (cutoff:3.500A) removed outlier: 4.116A pdb=" N ALA B 221 " --> pdb=" O LYS B 35 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'C' and resid 100 through 105 removed outlier: 6.956A pdb=" N GLY C 115 " --> pdb=" O MET C 101 " (cutoff:3.500A) removed outlier: 4.747A pdb=" N CYS C 103 " --> pdb=" O ALA C 113 " (cutoff:3.500A) removed outlier: 6.703A pdb=" N ALA C 113 " --> pdb=" O CYS C 103 " (cutoff:3.500A) removed outlier: 5.250A pdb=" N TYR C 105 " --> pdb=" O TYR C 111 " (cutoff:3.500A) removed outlier: 7.359A pdb=" N TYR C 111 " --> pdb=" O TYR C 105 " (cutoff:3.500A) removed outlier: 6.188A pdb=" N ASN C 125 " --> pdb=" O VAL C 135 " (cutoff:3.500A) removed outlier: 6.145A pdb=" N VAL C 135 " --> pdb=" O ASN C 125 " (cutoff:3.500A) Processing sheet with id= C, first strand: chain 'C' and resid 146 through 151 removed outlier: 6.569A pdb=" N SER C 160 " --> pdb=" O SER C 147 " (cutoff:3.500A) removed outlier: 4.542A pdb=" N CYS C 149 " --> pdb=" O VAL C 158 " (cutoff:3.500A) removed outlier: 6.791A pdb=" N VAL C 158 " --> pdb=" O CYS C 149 " (cutoff:3.500A) removed outlier: 4.972A pdb=" N PHE C 151 " --> pdb=" O GLN C 156 " (cutoff:3.500A) removed outlier: 7.432A pdb=" N GLN C 156 " --> pdb=" O PHE C 151 " (cutoff:3.500A) removed outlier: 3.673A pdb=" N ILE C 157 " --> pdb=" O TRP C 169 " (cutoff:3.500A) removed outlier: 7.033A pdb=" N ASP C 170 " --> pdb=" O GLN C 176 " (cutoff:3.500A) removed outlier: 6.817A pdb=" N GLN C 176 " --> pdb=" O ASP C 170 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'C' and resid 200 through 202 removed outlier: 3.589A pdb=" N GLN C 220 " --> pdb=" O LEU C 210 " (cutoff:3.500A) Processing sheet with id= E, first strand: chain 'C' and resid 229 through 234 removed outlier: 7.074A pdb=" N GLY C 244 " --> pdb=" O ASN C 230 " (cutoff:3.500A) removed outlier: 4.799A pdb=" N ILE C 232 " --> pdb=" O ALA C 242 " (cutoff:3.500A) removed outlier: 6.913A pdb=" N ALA C 242 " --> pdb=" O ILE C 232 " (cutoff:3.500A) removed outlier: 5.369A pdb=" N PHE C 234 " --> pdb=" O ALA C 240 " (cutoff:3.500A) removed outlier: 7.357A pdb=" N ALA C 240 " --> pdb=" O PHE C 234 " (cutoff:3.500A) removed outlier: 4.612A pdb=" N THR C 249 " --> pdb=" O SER C 245 " (cutoff:3.500A) removed outlier: 3.554A pdb=" N MET C 262 " --> pdb=" O LEU C 252 " (cutoff:3.500A) removed outlier: 5.820A pdb=" N ASP C 254 " --> pdb=" O GLU C 260 " (cutoff:3.500A) removed outlier: 5.985A pdb=" N GLU C 260 " --> pdb=" O ASP C 254 " (cutoff:3.500A) Processing sheet with id= F, first strand: chain 'C' and resid 273 through 278 removed outlier: 6.416A pdb=" N GLY C 288 " --> pdb=" O THR C 274 " (cutoff:3.500A) removed outlier: 4.533A pdb=" N VAL C 276 " --> pdb=" O LEU C 286 " (cutoff:3.500A) removed outlier: 6.953A pdb=" N LEU C 286 " --> pdb=" O VAL C 276 " (cutoff:3.500A) removed outlier: 4.581A pdb=" N PHE C 278 " --> pdb=" O LEU C 284 " (cutoff:3.500A) removed outlier: 6.716A pdb=" N LEU C 284 " --> pdb=" O PHE C 278 " (cutoff:3.500A) removed outlier: 4.501A pdb=" N ASN C 293 " --> pdb=" O TYR C 289 " (cutoff:3.500A) removed outlier: 5.629A pdb=" N ASP C 298 " --> pdb=" O ARG C 304 " (cutoff:3.500A) removed outlier: 6.262A pdb=" N ARG C 304 " --> pdb=" O ASP C 298 " (cutoff:3.500A) Processing sheet with id= G, first strand: chain 'C' and resid 47 through 51 removed outlier: 6.726A pdb=" N ILE C 338 " --> pdb=" O ARG C 48 " (cutoff:3.500A) removed outlier: 4.458A pdb=" N THR C 50 " --> pdb=" O LEU C 336 " (cutoff:3.500A) removed outlier: 6.747A pdb=" N LEU C 336 " --> pdb=" O THR C 50 " (cutoff:3.500A) Processing sheet with id= H, first strand: chain 'C' and resid 58 through 63 removed outlier: 7.118A pdb=" N ALA C 73 " --> pdb=" O TYR C 59 " (cutoff:3.500A) removed outlier: 4.558A pdb=" N MET C 61 " --> pdb=" O VAL C 71 " (cutoff:3.500A) removed outlier: 6.621A pdb=" N VAL C 71 " --> pdb=" O MET C 61 " (cutoff:3.500A) Processing sheet with id= I, first strand: chain 'C' and resid 78 through 83 removed outlier: 3.802A pdb=" N HIS C 91 " --> pdb=" O ILE C 81 " (cutoff:3.500A) removed outlier: 6.571A pdb=" N ASP C 83 " --> pdb=" O LYS C 89 " (cutoff:3.500A) removed outlier: 6.036A pdb=" N LYS C 89 " --> pdb=" O ASP C 83 " (cutoff:3.500A) Processing sheet with id= J, first strand: chain 'H' and resid 3 through 7 Processing sheet with id= K, first strand: chain 'H' and resid 10 through 12 Processing sheet with id= L, first strand: chain 'H' and resid 93 through 98 removed outlier: 3.592A pdb=" N VAL H 97 " --> pdb=" O HIS H 35 " (cutoff:3.500A) removed outlier: 3.818A pdb=" N HIS H 35 " --> pdb=" O VAL H 97 " (cutoff:3.500A) removed outlier: 3.561A pdb=" N ALA H 49 " --> pdb=" O TRP H 36 " (cutoff:3.500A) removed outlier: 5.409A pdb=" N ARG H 38 " --> pdb=" O TRP H 47 " (cutoff:3.500A) removed outlier: 4.846A pdb=" N TRP H 47 " --> pdb=" O ARG H 38 " (cutoff:3.500A) Processing sheet with id= M, first strand: chain 'H' and resid 146 through 149 removed outlier: 6.667A pdb=" N LYS H 244 " --> pdb=" O VAL H 147 " (cutoff:3.500A) removed outlier: 7.335A pdb=" N VAL H 149 " --> pdb=" O LYS H 244 " (cutoff:3.500A) removed outlier: 5.801A pdb=" N GLU H 246 " --> pdb=" O VAL H 149 " (cutoff:3.500A) No H-bonds generated for sheet with id= M Processing sheet with id= N, first strand: chain 'H' and resid 155 through 160 removed outlier: 3.523A pdb=" N ALA H 211 " --> pdb=" O SER H 208 " (cutoff:3.500A) Processing sheet with id= O, first strand: chain 'H' and resid 225 through 231 removed outlier: 6.214A pdb=" N LEU H 178 " --> pdb=" O LEU H 187 " (cutoff:3.500A) removed outlier: 5.547A pdb=" N LEU H 187 " --> pdb=" O LEU H 178 " (cutoff:3.500A) 360 hydrogen bonds defined for protein. 960 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: 2.15 Time building geometry restraints manager: 3.60 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.19 - 1.32: 1277 1.32 - 1.44: 2223 1.44 - 1.56: 4498 1.56 - 1.69: 1 1.69 - 1.81: 65 Bond restraints: 8064 Sorted by residual: bond pdb=" C PRO A 298 " pdb=" O PRO A 298 " ideal model delta sigma weight residual 1.234 1.193 0.041 1.00e-02 1.00e+04 1.72e+01 bond pdb=" N LEU A 276 " pdb=" CA LEU A 276 " ideal model delta sigma weight residual 1.456 1.494 -0.038 1.18e-02 7.18e+03 1.05e+01 bond pdb=" C PRO A 270 " pdb=" N PRO A 271 " ideal model delta sigma weight residual 1.333 1.361 -0.028 1.01e-02 9.80e+03 7.45e+00 bond pdb=" C PRO A 170 " pdb=" O PRO A 170 " ideal model delta sigma weight residual 1.237 1.205 0.032 1.26e-02 6.30e+03 6.57e+00 bond pdb=" N ASP H 223 " pdb=" CA ASP H 223 " ideal model delta sigma weight residual 1.455 1.487 -0.032 1.33e-02 5.65e+03 5.81e+00 ... (remaining 8059 not shown) Histogram of bond angle deviations from ideal: 99.51 - 106.41: 186 106.41 - 113.31: 4416 113.31 - 120.21: 2714 120.21 - 127.11: 3610 127.11 - 134.01: 85 Bond angle restraints: 11011 Sorted by residual: angle pdb=" N PRO A 170 " pdb=" CA PRO A 170 " pdb=" C PRO A 170 " ideal model delta sigma weight residual 113.53 108.21 5.32 1.39e+00 5.18e-01 1.47e+01 angle pdb=" N ILE C 270 " pdb=" CA ILE C 270 " pdb=" C ILE C 270 " ideal model delta sigma weight residual 113.53 109.85 3.68 9.80e-01 1.04e+00 1.41e+01 angle pdb=" O ILE A 169 " pdb=" C ILE A 169 " pdb=" N PRO A 170 " ideal model delta sigma weight residual 121.10 117.10 4.00 1.14e+00 7.69e-01 1.23e+01 angle pdb=" CB MET H 83 " pdb=" CG MET H 83 " pdb=" SD MET H 83 " ideal model delta sigma weight residual 112.70 122.36 -9.66 3.00e+00 1.11e-01 1.04e+01 angle pdb=" CA ASP H 223 " pdb=" CB ASP H 223 " pdb=" CG ASP H 223 " ideal model delta sigma weight residual 112.60 115.78 -3.18 1.00e+00 1.00e+00 1.01e+01 ... (remaining 11006 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 26.57: 4427 26.57 - 53.13: 229 53.13 - 79.70: 15 79.70 - 106.27: 2 106.27 - 132.83: 2 Dihedral angle restraints: 4675 sinusoidal: 1525 harmonic: 3150 Sorted by residual: dihedral pdb=" CD DAR D 406 " pdb=" NE DAR D 406 " pdb=" CZ DAR D 406 " pdb=" NH1 DAR D 406 " ideal model delta sinusoidal sigma weight residual 0.00 -132.83 132.83 1 1.00e+01 1.00e-02 1.61e+02 dihedral pdb=" CA GLY A 51 " pdb=" C GLY A 51 " pdb=" N VAL A 52 " pdb=" CA VAL A 52 " ideal model delta harmonic sigma weight residual -180.00 -158.56 -21.44 0 5.00e+00 4.00e-02 1.84e+01 dihedral pdb=" CG ALC D 405 " pdb=" CD1 ALC D 405 " pdb=" CE1 ALC D 405 " pdb=" CZ ALC D 405 " ideal model delta sinusoidal sigma weight residual -55.84 61.55 -117.39 1 3.00e+01 1.11e-03 1.56e+01 ... (remaining 4672 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.047: 1028 0.047 - 0.095: 219 0.095 - 0.142: 56 0.142 - 0.190: 2 0.190 - 0.237: 2 Chirality restraints: 1307 Sorted by residual: chirality pdb=" CB ILE A 169 " pdb=" CA ILE A 169 " pdb=" CG1 ILE A 169 " pdb=" CG2 ILE A 169 " both_signs ideal model delta sigma weight residual False 2.64 2.41 0.24 2.00e-01 2.50e+01 1.41e+00 chirality pdb=" CA DAR D 406 " pdb=" N DAR D 406 " pdb=" C DAR D 406 " pdb=" CB DAR D 406 " both_signs ideal model delta sigma weight residual False -2.51 -2.72 0.21 2.00e-01 2.50e+01 1.10e+00 chirality pdb=" CB ILE A 266 " pdb=" CA ILE A 266 " pdb=" CG1 ILE A 266 " pdb=" CG2 ILE A 266 " both_signs ideal model delta sigma weight residual False 2.64 2.46 0.18 2.00e-01 2.50e+01 8.45e-01 ... (remaining 1304 not shown) Planarity restraints: 1398 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CD DAR D 406 " 0.801 9.50e-02 1.11e+02 3.59e-01 7.84e+01 pdb=" NE DAR D 406 " -0.047 2.00e-02 2.50e+03 pdb=" CZ DAR D 406 " -0.010 2.00e-02 2.50e+03 pdb=" NH1 DAR D 406 " 0.024 2.00e-02 2.50e+03 pdb=" NH2 DAR D 406 " -0.002 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C GLN H 13 " -0.058 5.00e-02 4.00e+02 8.77e-02 1.23e+01 pdb=" N PRO H 14 " 0.152 5.00e-02 4.00e+02 pdb=" CA PRO H 14 " -0.046 5.00e-02 4.00e+02 pdb=" CD PRO H 14 " -0.047 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C ASP H 74 " 0.057 5.00e-02 4.00e+02 8.49e-02 1.15e+01 pdb=" N PRO H 75 " -0.147 5.00e-02 4.00e+02 pdb=" CA PRO H 75 " 0.043 5.00e-02 4.00e+02 pdb=" CD PRO H 75 " 0.047 5.00e-02 4.00e+02 ... (remaining 1395 not shown) Histogram of nonbonded interaction distances: 2.21 - 2.75: 843 2.75 - 3.29: 7662 3.29 - 3.82: 12807 3.82 - 4.36: 13887 4.36 - 4.90: 25620 Nonbonded interactions: 60819 Sorted by model distance: nonbonded pdb=" O ARG H 67 " pdb=" OG1 THR H 84 " model vdw 2.212 2.440 nonbonded pdb=" O SER B 310 " pdb=" OG SER B 310 " model vdw 2.226 2.440 nonbonded pdb=" NH1 ARG H 180 " pdb=" O GLU H 222 " model vdw 2.306 2.520 nonbonded pdb=" OG SER C 245 " pdb=" OD1 ASP C 247 " model vdw 2.325 2.440 nonbonded pdb=" OG1 THR C 128 " pdb=" O ASN C 132 " model vdw 2.351 2.440 ... (remaining 60814 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 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: 4.070 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.010 Extract box with map and model: 10.090 Check model and map are aligned: 0.130 Set scattering table: 0.080 Process input model: 26.300 Find NCS groups from input model: 0.240 Set up NCS constraints: 0.030 Set refine NCS operators: 0.000 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:2.720 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 43.670 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7842 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.003 0.045 8064 Z= 0.173 Angle : 0.598 9.662 11011 Z= 0.315 Chirality : 0.042 0.237 1307 Planarity : 0.011 0.359 1398 Dihedral : 14.890 132.831 2637 Min Nonbonded Distance : 2.212 Molprobity Statistics. All-atom Clashscore : 6.51 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.37 % Favored : 96.63 % Rotamer: Outliers : 0.40 % Allowed : 8.30 % Favored : 91.30 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.29 (0.27), residues: 1067 helix: 1.87 (0.28), residues: 364 sheet: -0.21 (0.32), residues: 291 loop : -1.00 (0.30), residues: 412 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.031 0.001 TRP A 268 HIS 0.003 0.001 HIS B 189 PHE 0.025 0.001 PHE H 68 TYR 0.017 0.001 TYR B 297 ARG 0.006 0.001 ARG A 148 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 97 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 94 time to evaluate : 0.910 Fit side-chains outliers start: 3 outliers final: 1 residues processed: 96 average time/residue: 0.2122 time to fit residues: 27.7263 Evaluate side-chains 93 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 92 time to evaluate : 0.813 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0750 time to fit residues: 1.3968 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 89 optimal weight: 0.9990 chunk 80 optimal weight: 1.9990 chunk 44 optimal weight: 0.2980 chunk 27 optimal weight: 4.9990 chunk 54 optimal weight: 0.7980 chunk 42 optimal weight: 1.9990 chunk 83 optimal weight: 4.9990 chunk 32 optimal weight: 5.9990 chunk 50 optimal weight: 1.9990 chunk 61 optimal weight: 0.9990 chunk 96 optimal weight: 5.9990 overall best weight: 1.0186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 283 ASN ** C 32 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 237 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** H 179 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7982 moved from start: 0.1155 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.044 8064 Z= 0.321 Angle : 0.619 7.397 11011 Z= 0.320 Chirality : 0.044 0.179 1307 Planarity : 0.005 0.058 1398 Dihedral : 7.520 179.734 1169 Min Nonbonded Distance : 2.174 Molprobity Statistics. All-atom Clashscore : 5.66 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.12 % Favored : 95.88 % Rotamer: Outliers : 1.61 % Allowed : 14.06 % Favored : 84.34 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.12 (0.27), residues: 1067 helix: 1.57 (0.29), residues: 362 sheet: -0.18 (0.31), residues: 302 loop : -1.03 (0.31), residues: 403 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A 74 HIS 0.007 0.002 HIS C 142 PHE 0.038 0.002 PHE H 68 TYR 0.017 0.002 TYR H 190 ARG 0.005 0.001 ARG A 237 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 117 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 105 time to evaluate : 0.875 Fit side-chains outliers start: 12 outliers final: 10 residues processed: 113 average time/residue: 0.1930 time to fit residues: 30.3173 Evaluate side-chains 109 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 99 time to evaluate : 0.940 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.0724 time to fit residues: 2.5321 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 53 optimal weight: 0.0470 chunk 29 optimal weight: 0.0670 chunk 80 optimal weight: 1.9990 chunk 65 optimal weight: 0.8980 chunk 26 optimal weight: 0.9990 chunk 96 optimal weight: 4.9990 chunk 104 optimal weight: 1.9990 chunk 85 optimal weight: 0.9980 chunk 95 optimal weight: 0.4980 chunk 32 optimal weight: 2.9990 chunk 77 optimal weight: 0.6980 overall best weight: 0.4416 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 32 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7932 moved from start: 0.1277 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 8064 Z= 0.176 Angle : 0.528 7.354 11011 Z= 0.271 Chirality : 0.041 0.169 1307 Planarity : 0.004 0.049 1398 Dihedral : 7.186 178.940 1169 Min Nonbonded Distance : 2.196 Molprobity Statistics. All-atom Clashscore : 5.33 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.56 % Favored : 96.44 % Rotamer: Outliers : 1.20 % Allowed : 17.14 % Favored : 81.66 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.11 (0.27), residues: 1067 helix: 1.53 (0.29), residues: 362 sheet: -0.19 (0.31), residues: 301 loop : -1.01 (0.31), residues: 404 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A 74 HIS 0.003 0.001 HIS C 142 PHE 0.030 0.001 PHE H 68 TYR 0.011 0.001 TYR H 190 ARG 0.004 0.000 ARG A 237 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 111 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 102 time to evaluate : 0.843 Fit side-chains outliers start: 9 outliers final: 7 residues processed: 106 average time/residue: 0.1991 time to fit residues: 29.3104 Evaluate side-chains 101 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 94 time to evaluate : 0.843 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0783 time to fit residues: 2.2096 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 95 optimal weight: 0.7980 chunk 72 optimal weight: 0.9980 chunk 50 optimal weight: 0.8980 chunk 10 optimal weight: 0.2980 chunk 46 optimal weight: 0.9980 chunk 64 optimal weight: 0.9980 chunk 96 optimal weight: 1.9990 chunk 102 optimal weight: 1.9990 chunk 91 optimal weight: 3.9990 chunk 27 optimal weight: 1.9990 chunk 85 optimal weight: 0.8980 overall best weight: 0.7780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 194 ASN ** C 32 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 237 ASN C 239 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7987 moved from start: 0.1672 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.036 8064 Z= 0.251 Angle : 0.568 8.180 11011 Z= 0.291 Chirality : 0.043 0.153 1307 Planarity : 0.004 0.045 1398 Dihedral : 7.210 179.013 1169 Min Nonbonded Distance : 2.185 Molprobity Statistics. All-atom Clashscore : 5.46 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.03 % Favored : 95.97 % Rotamer: Outliers : 1.74 % Allowed : 19.81 % Favored : 78.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.11 (0.27), residues: 1067 helix: 1.43 (0.29), residues: 358 sheet: -0.08 (0.31), residues: 305 loop : -0.98 (0.31), residues: 404 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP A 74 HIS 0.004 0.001 HIS C 142 PHE 0.033 0.002 PHE H 68 TYR 0.014 0.001 TYR H 190 ARG 0.003 0.000 ARG H 18 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 118 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 105 time to evaluate : 0.923 Fit side-chains outliers start: 13 outliers final: 7 residues processed: 111 average time/residue: 0.2044 time to fit residues: 31.2922 Evaluate side-chains 109 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 102 time to evaluate : 0.941 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0788 time to fit residues: 2.2707 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 58 optimal weight: 0.7980 chunk 1 optimal weight: 0.3980 chunk 76 optimal weight: 0.7980 chunk 42 optimal weight: 0.8980 chunk 87 optimal weight: 0.9990 chunk 70 optimal weight: 0.6980 chunk 0 optimal weight: 0.5980 chunk 52 optimal weight: 2.9990 chunk 91 optimal weight: 2.9990 chunk 25 optimal weight: 1.9990 chunk 34 optimal weight: 1.9990 overall best weight: 0.6580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 183 HIS C 237 ASN C 239 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7978 moved from start: 0.1800 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 8064 Z= 0.220 Angle : 0.541 7.481 11011 Z= 0.278 Chirality : 0.042 0.174 1307 Planarity : 0.004 0.043 1398 Dihedral : 7.127 179.135 1169 Min Nonbonded Distance : 2.191 Molprobity Statistics. All-atom Clashscore : 5.46 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.75 % Favored : 96.25 % Rotamer: Outliers : 1.07 % Allowed : 21.02 % Favored : 77.91 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.12 (0.27), residues: 1067 helix: 1.45 (0.29), residues: 358 sheet: -0.14 (0.32), residues: 299 loop : -0.91 (0.31), residues: 410 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A 74 HIS 0.003 0.001 HIS C 142 PHE 0.030 0.001 PHE H 68 TYR 0.012 0.001 TYR H 190 ARG 0.004 0.000 ARG H 18 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 114 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 106 time to evaluate : 0.867 Fit side-chains outliers start: 8 outliers final: 3 residues processed: 111 average time/residue: 0.2021 time to fit residues: 30.5353 Evaluate side-chains 102 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 99 time to evaluate : 0.815 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0695 time to fit residues: 1.4754 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 92 optimal weight: 1.9990 chunk 20 optimal weight: 0.5980 chunk 60 optimal weight: 0.9990 chunk 25 optimal weight: 0.8980 chunk 102 optimal weight: 3.9990 chunk 85 optimal weight: 0.9990 chunk 47 optimal weight: 3.9990 chunk 8 optimal weight: 1.9990 chunk 33 optimal weight: 0.0470 chunk 53 optimal weight: 2.9990 chunk 98 optimal weight: 0.9980 overall best weight: 0.7080 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 55 ASN A 77 ASN C 183 HIS C 237 ASN C 239 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7993 moved from start: 0.2015 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.035 8064 Z= 0.234 Angle : 0.556 7.498 11011 Z= 0.284 Chirality : 0.042 0.179 1307 Planarity : 0.004 0.042 1398 Dihedral : 7.113 179.159 1169 Min Nonbonded Distance : 2.167 Molprobity Statistics. All-atom Clashscore : 5.53 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.31 % Favored : 95.69 % Rotamer: Outliers : 1.61 % Allowed : 21.69 % Favored : 76.71 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.11 (0.27), residues: 1067 helix: 1.41 (0.29), residues: 358 sheet: -0.22 (0.32), residues: 300 loop : -0.82 (0.31), residues: 409 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP A 74 HIS 0.007 0.001 HIS C 183 PHE 0.034 0.002 PHE H 68 TYR 0.016 0.001 TYR B 297 ARG 0.004 0.000 ARG H 18 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 113 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 101 time to evaluate : 0.826 Fit side-chains outliers start: 12 outliers final: 8 residues processed: 110 average time/residue: 0.1966 time to fit residues: 29.5837 Evaluate side-chains 104 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 96 time to evaluate : 0.805 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0799 time to fit residues: 2.2558 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 11 optimal weight: 0.9990 chunk 58 optimal weight: 0.5980 chunk 74 optimal weight: 1.9990 chunk 86 optimal weight: 1.9990 chunk 57 optimal weight: 0.7980 chunk 102 optimal weight: 3.9990 chunk 63 optimal weight: 0.7980 chunk 62 optimal weight: 3.9990 chunk 47 optimal weight: 0.4980 chunk 40 optimal weight: 0.7980 chunk 61 optimal weight: 0.9980 overall best weight: 0.6980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 77 ASN B 205 GLN C 183 HIS C 237 ASN C 239 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7993 moved from start: 0.2127 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.042 8064 Z= 0.231 Angle : 0.553 7.501 11011 Z= 0.284 Chirality : 0.042 0.194 1307 Planarity : 0.004 0.042 1398 Dihedral : 7.108 179.232 1169 Min Nonbonded Distance : 2.182 Molprobity Statistics. All-atom Clashscore : 6.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.12 % Favored : 95.88 % Rotamer: Outliers : 0.94 % Allowed : 22.76 % Favored : 76.31 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.09 (0.27), residues: 1067 helix: 1.38 (0.29), residues: 358 sheet: -0.26 (0.31), residues: 305 loop : -0.79 (0.32), residues: 404 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.001 TRP A 74 HIS 0.003 0.001 HIS C 142 PHE 0.031 0.002 PHE H 68 TYR 0.011 0.001 TYR H 190 ARG 0.007 0.000 ARG H 18 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 98 time to evaluate : 0.959 Fit side-chains outliers start: 7 outliers final: 4 residues processed: 102 average time/residue: 0.2189 time to fit residues: 30.5104 Evaluate side-chains 99 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 95 time to evaluate : 1.005 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0739 time to fit residues: 1.9668 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 30 optimal weight: 0.9990 chunk 20 optimal weight: 2.9990 chunk 19 optimal weight: 0.7980 chunk 64 optimal weight: 0.9980 chunk 69 optimal weight: 0.4980 chunk 50 optimal weight: 0.6980 chunk 9 optimal weight: 1.9990 chunk 80 optimal weight: 0.6980 chunk 92 optimal weight: 0.9980 chunk 97 optimal weight: 3.9990 chunk 89 optimal weight: 0.4980 overall best weight: 0.6380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 77 ASN B 194 ASN C 183 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7989 moved from start: 0.2222 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 8064 Z= 0.217 Angle : 0.549 7.502 11011 Z= 0.281 Chirality : 0.042 0.211 1307 Planarity : 0.004 0.041 1398 Dihedral : 7.079 179.363 1169 Min Nonbonded Distance : 2.190 Molprobity Statistics. All-atom Clashscore : 5.99 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.94 % Favored : 96.06 % Rotamer: Outliers : 0.94 % Allowed : 23.03 % Favored : 76.04 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.02 (0.27), residues: 1067 helix: 1.27 (0.29), residues: 364 sheet: -0.27 (0.31), residues: 305 loop : -0.84 (0.32), residues: 398 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP A 74 HIS 0.006 0.001 HIS C 183 PHE 0.037 0.001 PHE H 68 TYR 0.016 0.001 TYR B 297 ARG 0.008 0.000 ARG H 18 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 106 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 99 time to evaluate : 0.903 Fit side-chains outliers start: 7 outliers final: 4 residues processed: 102 average time/residue: 0.2001 time to fit residues: 28.1300 Evaluate side-chains 103 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 99 time to evaluate : 0.811 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0769 time to fit residues: 1.7373 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 95 optimal weight: 0.8980 chunk 97 optimal weight: 0.8980 chunk 57 optimal weight: 0.6980 chunk 41 optimal weight: 0.0070 chunk 74 optimal weight: 1.9990 chunk 29 optimal weight: 0.3980 chunk 86 optimal weight: 0.9990 chunk 90 optimal weight: 0.4980 chunk 94 optimal weight: 1.9990 chunk 62 optimal weight: 1.9990 chunk 100 optimal weight: 1.9990 overall best weight: 0.4998 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 77 ASN B 194 ASN B 270 ASN C 183 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7970 moved from start: 0.2254 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.037 8064 Z= 0.187 Angle : 0.539 9.331 11011 Z= 0.273 Chirality : 0.042 0.199 1307 Planarity : 0.004 0.040 1398 Dihedral : 7.026 179.657 1169 Min Nonbonded Distance : 2.193 Molprobity Statistics. All-atom Clashscore : 5.66 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.84 % Favored : 96.16 % Rotamer: Outliers : 0.13 % Allowed : 23.83 % Favored : 76.04 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.86 % Cis-general : 0.00 % 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.05 (0.27), residues: 1067 helix: 1.28 (0.29), residues: 366 sheet: -0.21 (0.32), residues: 296 loop : -0.85 (0.31), residues: 405 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP A 74 HIS 0.003 0.001 HIS C 142 PHE 0.022 0.001 PHE H 68 TYR 0.010 0.001 TYR H 190 ARG 0.009 0.000 ARG H 18 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2134 Ramachandran restraints generated. 1067 Oldfield, 0 Emsley, 1067 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 99 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 98 time to evaluate : 0.915 Fit side-chains outliers start: 1 outliers final: 0 residues processed: 99 average time/residue: 0.2065 time to fit residues: 28.1821 Evaluate side-chains 97 residues out of total 926 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 97 time to evaluate : 0.916 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.2049 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/iotbx/cli_parser.py", line 931, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/programs/real_space_refine.py", line 191, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 312, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 737, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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: 55.8 > 50: distance: 77 - 145: 19.760 distance: 86 - 135: 26.615 distance: 89 - 132: 13.316 distance: 100 - 103: 14.409 distance: 103 - 104: 41.993 distance: 104 - 105: 27.465 distance: 105 - 106: 28.542 distance: 105 - 107: 16.704 distance: 107 - 108: 16.742 distance: 108 - 109: 15.400 distance: 108 - 111: 32.532 distance: 109 - 110: 19.789 distance: 109 - 112: 4.754 distance: 112 - 113: 7.666 distance: 113 - 114: 10.409 distance: 113 - 116: 8.952 distance: 114 - 115: 19.628 distance: 114 - 119: 10.355 distance: 116 - 117: 21.951 distance: 116 - 118: 23.533 distance: 119 - 120: 15.353 distance: 120 - 121: 19.170 distance: 120 - 123: 19.786 distance: 121 - 122: 33.096 distance: 121 - 126: 33.275 distance: 123 - 124: 25.602 distance: 123 - 125: 30.707 distance: 126 - 127: 22.627 distance: 126 - 226: 14.095 distance: 127 - 128: 20.644 distance: 127 - 130: 23.469 distance: 128 - 129: 23.059 distance: 128 - 132: 16.690 distance: 129 - 223: 22.431 distance: 130 - 131: 18.388 distance: 132 - 133: 12.939 distance: 133 - 134: 38.877 distance: 133 - 136: 18.798 distance: 134 - 135: 13.953 distance: 134 - 137: 24.463 distance: 137 - 138: 16.694 distance: 137 - 214: 11.907 distance: 138 - 139: 31.046 distance: 138 - 141: 15.502 distance: 139 - 140: 42.755 distance: 139 - 145: 23.404 distance: 140 - 211: 25.556 distance: 141 - 142: 43.324 distance: 142 - 143: 10.810 distance: 142 - 144: 27.250 distance: 145 - 146: 5.113 distance: 146 - 147: 27.806 distance: 146 - 149: 9.390 distance: 147 - 148: 11.924 distance: 147 - 159: 33.388 distance: 149 - 150: 43.064 distance: 150 - 151: 24.822 distance: 150 - 152: 18.048 distance: 151 - 153: 31.988 distance: 152 - 154: 14.509 distance: 152 - 155: 9.367 distance: 153 - 154: 5.742 distance: 154 - 156: 19.688 distance: 155 - 157: 18.612 distance: 156 - 158: 7.219 distance: 157 - 158: 17.051 distance: 159 - 160: 40.208 distance: 160 - 161: 34.389 distance: 160 - 163: 9.351 distance: 161 - 162: 26.134 distance: 161 - 167: 32.042 distance: 163 - 164: 34.006 distance: 164 - 165: 28.687 distance: 164 - 166: 29.593 distance: 167 - 168: 23.954 distance: 168 - 169: 21.394 distance: 168 - 171: 26.249 distance: 169 - 170: 18.295 distance: 169 - 175: 5.566 distance: 171 - 172: 17.704 distance: 171 - 173: 18.649 distance: 172 - 174: 26.364