Starting phenix.real_space_refine (version: 1.21rc1) on Fri Apr 21 19:11:31 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7lb8_23251/04_2023/7lb8_23251.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7lb8_23251/04_2023/7lb8_23251.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.4 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7lb8_23251/04_2023/7lb8_23251.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7lb8_23251/04_2023/7lb8_23251.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7lb8_23251/04_2023/7lb8_23251.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7lb8_23251/04_2023/7lb8_23251.pdb" } resolution = 3.4 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.000 sd= 0.003 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 64 5.16 5 C 6779 2.51 5 N 1785 2.21 5 O 1811 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "B PHE 79": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 89": "OE1" <-> "OE2" Residue "B TYR 156": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 167": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 260": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 272": "NH1" <-> "NH2" Residue "B ARG 276": "NH1" <-> "NH2" Residue "B PHE 350": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 385": "OD1" <-> "OD2" Residue "B GLU 423": "OE1" <-> "OE2" Residue "B PHE 447": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 474": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 504": "NH1" <-> "NH2" Residue "B TYR 517": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 521": "OD1" <-> "OD2" Residue "B ARG 543": "NH1" <-> "NH2" Residue "B PHE 589": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 603": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 635": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 651": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 653": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ASP 61": "OD1" <-> "OD2" Residue "D ARG 84": "NH1" <-> "NH2" Residue "D GLU 90": "OE1" <-> "OE2" Residue "D PHE 99": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D TYR 106": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D PHE 122": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ASP 126": "OD1" <-> "OD2" Residue "D GLU 140": "OE1" <-> "OE2" Residue "D TYR 158": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D GLU 159": "OE1" <-> "OE2" Residue "D PHE 216": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D TYR 230": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ASP 232": "OD1" <-> "OD2" Residue "D PHE 238": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ASP 241": "OD1" <-> "OD2" Residue "D ASP 247": "OD1" <-> "OD2" Residue "D PHE 274": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D PHE 284": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U ARG 20": "NH1" <-> "NH2" Residue "U PHE 83": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U TYR 111": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U ARG 119": "NH1" <-> "NH2" Residue "U GLU 130": "OE1" <-> "OE2" Residue "U ASP 165": "OD1" <-> "OD2" Residue "U ASP 179": "OD1" <-> "OD2" Residue "U TYR 214": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U TYR 245": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U PHE 263": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U TYR 265": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C GLU 67": "OE1" <-> "OE2" Residue "C PHE 83": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C GLU 104": "OE1" <-> "OE2" Residue "C GLU 126": "OE1" <-> "OE2" Residue "C ASP 146": "OD1" <-> "OD2" Residue "C ARG 155": "NH1" <-> "NH2" Residue "C ASP 165": "OD1" <-> "OD2" Residue "C ASP 179": "OD1" <-> "OD2" Residue "C ASP 185": "OD1" <-> "OD2" Residue "C ASP 207": "OD1" <-> "OD2" Residue "C TYR 214": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C GLU 234": "OE1" <-> "OE2" Residue "C GLU 242": "OE1" <-> "OE2" Residue "C PHE 263": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-4938/modules/chem_data/mon_lib" Total number of atoms: 10439 Number of models: 1 Model: "" Number of chains: 4 Chain: "B" Number of atoms: 4610 Number of conformers: 1 Conformer: "" Number of residues, atoms: 630, 4610 Classifications: {'peptide': 630} Incomplete info: {'truncation_to_alanine': 25} Link IDs: {'PTRANS': 30, 'TRANS': 599} Chain breaks: 1 Unresolved chain link angles: 3 Unresolved non-hydrogen bonds: 100 Unresolved non-hydrogen angles: 129 Unresolved non-hydrogen dihedrals: 90 Unresolved non-hydrogen chiralities: 5 Planarities with less than four sites: {'PHE:plan': 1, 'ASN:plan1': 3, 'TRP:plan': 2, 'ASP:plan': 3, 'ARG:plan': 3} Unresolved non-hydrogen planarities: 57 Chain: "D" Number of atoms: 2085 Number of conformers: 1 Conformer: "" Number of residues, atoms: 265, 2085 Classifications: {'peptide': 265} Link IDs: {'PTRANS': 20, 'TRANS': 244} Chain: "U" Number of atoms: 1896 Number of conformers: 1 Conformer: "" Number of residues, atoms: 255, 1896 Classifications: {'peptide': 255} Incomplete info: {'truncation_to_alanine': 12} Link IDs: {'PTRANS': 17, 'TRANS': 237} Unresolved chain link angles: 1 Unresolved non-hydrogen bonds: 52 Unresolved non-hydrogen angles: 63 Unresolved non-hydrogen dihedrals: 44 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'GLN:plan1': 1, 'GLU:plan': 2, 'ARG:plan': 3, 'HIS:plan': 1} Unresolved non-hydrogen planarities: 32 Chain: "C" Number of atoms: 1848 Number of conformers: 1 Conformer: "" Number of residues, atoms: 255, 1848 Classifications: {'peptide': 255} Incomplete info: {'truncation_to_alanine': 24} Link IDs: {'PTRANS': 17, 'TRANS': 237} Unresolved chain link angles: 4 Unresolved non-hydrogen bonds: 106 Unresolved non-hydrogen angles: 132 Unresolved non-hydrogen dihedrals: 92 Unresolved non-hydrogen chiralities: 4 Planarities with less than four sites: {'PHE:plan': 1, 'GLU:plan': 3, 'ARG:plan': 6, 'ASP:plan': 1, 'HIS:plan': 3} Unresolved non-hydrogen planarities: 66 Residues with excluded nonbonded symmetry interactions: 35 residue: pdb=" N ILE B 111 " occ=0.00 ... (6 atoms not shown) pdb=" CD1 ILE B 111 " occ=0.00 residue: pdb=" N PRO B 112 " occ=0.00 ... (3 atoms not shown) pdb=" CB PRO B 112 " occ=0.00 residue: pdb=" N GLY B 113 " occ=0.00 ... (2 atoms not shown) pdb=" O GLY B 113 " occ=0.00 residue: pdb=" N VAL B 135 " occ=0.00 ... (3 atoms not shown) pdb=" CB VAL B 135 " occ=0.00 residue: pdb=" N ALA B 136 " occ=0.00 ... (3 atoms not shown) pdb=" CB ALA B 136 " occ=0.00 residue: pdb=" N TRP B 137 " occ=0.00 ... (3 atoms not shown) pdb=" CB TRP B 137 " occ=0.00 residue: pdb=" N GLY B 138 " occ=0.00 ... (2 atoms not shown) pdb=" O GLY B 138 " occ=0.00 residue: pdb=" N LEU B 217 " occ=0.00 ... (6 atoms not shown) pdb=" CD2 LEU B 217 " occ=0.00 residue: pdb=" N ASP B 218 " occ=0.00 ... (3 atoms not shown) pdb=" CB ASP B 218 " occ=0.00 residue: pdb=" N ASP B 219 " occ=0.00 ... (6 atoms not shown) pdb=" OD2 ASP B 219 " occ=0.00 residue: pdb=" N TRP B 449 " occ=0.00 ... (3 atoms not shown) pdb=" CB TRP B 449 " occ=0.00 residue: pdb=" N LEU B 450 " occ=0.00 ... (3 atoms not shown) pdb=" CB LEU B 450 " occ=0.00 ... (remaining 23 not shown) Time building chain proxies: 6.08, per 1000 atoms: 0.58 Number of scatterers: 10439 At special positions: 0 Unit cell: (81.498, 95.88, 140.624, 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 64 16.00 O 1811 8.00 N 1785 7.00 C 6779 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 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: 4.10 Conformation dependent library (CDL) restraints added in 1.7 seconds 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2556 Finding SS restraints... Secondary structure from input PDB file: 66 helices and 8 sheets defined 61.3% alpha, 6.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.63 Creating SS restraints... Processing helix chain 'B' and resid 6 through 32 Processing helix chain 'B' and resid 33 through 35 No H-bonds generated for 'chain 'B' and resid 33 through 35' Processing helix chain 'B' and resid 36 through 42 Processing helix chain 'B' and resid 49 through 54 removed outlier: 3.550A pdb=" N PHE B 54 " --> pdb=" O GLU B 50 " (cutoff:3.500A) Processing helix chain 'B' and resid 54 through 83 Proline residue: B 60 - end of helix removed outlier: 3.551A pdb=" N LEU B 83 " --> pdb=" O PHE B 79 " (cutoff:3.500A) Processing helix chain 'B' and resid 94 through 108 removed outlier: 3.813A pdb=" N LEU B 108 " --> pdb=" O THR B 104 " (cutoff:3.500A) Processing helix chain 'B' and resid 114 through 134 removed outlier: 3.533A pdb=" N GLN B 118 " --> pdb=" O ALA B 114 " (cutoff:3.500A) Processing helix chain 'B' and resid 142 through 168 Processing helix chain 'B' and resid 172 through 179 Processing helix chain 'B' and resid 188 through 207 Proline residue: B 196 - end of helix Processing helix chain 'B' and resid 208 through 217 removed outlier: 3.555A pdb=" N LEU B 212 " --> pdb=" O LEU B 208 " (cutoff:3.500A) Processing helix chain 'B' and resid 218 through 225 Processing helix chain 'B' and resid 228 through 251 Processing helix chain 'B' and resid 256 through 259 Processing helix chain 'B' and resid 260 through 269 Processing helix chain 'B' and resid 272 through 302 removed outlier: 3.697A pdb=" N ARG B 276 " --> pdb=" O ARG B 272 " (cutoff:3.500A) Processing helix chain 'B' and resid 307 through 323 Proline residue: B 317 - end of helix Processing helix chain 'B' and resid 349 through 370 Processing helix chain 'B' and resid 380 through 387 removed outlier: 3.753A pdb=" N ASP B 384 " --> pdb=" O GLY B 380 " (cutoff:3.500A) removed outlier: 3.992A pdb=" N ASP B 385 " --> pdb=" O ALA B 381 " (cutoff:3.500A) Processing helix chain 'B' and resid 390 through 416 Processing helix chain 'B' and resid 421 through 425 Processing helix chain 'B' and resid 426 through 442 Processing helix chain 'B' and resid 454 through 473 removed outlier: 3.511A pdb=" N GLY B 473 " --> pdb=" O ALA B 469 " (cutoff:3.500A) Processing helix chain 'B' and resid 475 through 500 removed outlier: 3.528A pdb=" N SER B 500 " --> pdb=" O MET B 496 " (cutoff:3.500A) Processing helix chain 'B' and resid 504 through 513 Processing helix chain 'B' and resid 520 through 538 Processing helix chain 'B' and resid 541 through 548 Processing helix chain 'B' and resid 552 through 559 removed outlier: 3.593A pdb=" N ARG B 556 " --> pdb=" O GLY B 552 " (cutoff:3.500A) removed outlier: 3.609A pdb=" N GLY B 559 " --> pdb=" O ALA B 555 " (cutoff:3.500A) Processing helix chain 'B' and resid 561 through 584 removed outlier: 3.760A pdb=" N ALA B 568 " --> pdb=" O PRO B 564 " (cutoff:3.500A) Processing helix chain 'B' and resid 588 through 602 Proline residue: B 595 - end of helix removed outlier: 3.759A pdb=" N GLY B 602 " --> pdb=" O ALA B 598 " (cutoff:3.500A) Processing helix chain 'B' and resid 605 through 631 Processing helix chain 'B' and resid 638 through 658 removed outlier: 3.697A pdb=" N LEU B 642 " --> pdb=" O PRO B 638 " (cutoff:3.500A) removed outlier: 3.577A pdb=" N THR B 644 " --> pdb=" O GLY B 640 " (cutoff:3.500A) Proline residue: B 649 - end of helix Processing helix chain 'D' and resid 34 through 36 No H-bonds generated for 'chain 'D' and resid 34 through 36' Processing helix chain 'D' and resid 42 through 53 Processing helix chain 'D' and resid 61 through 69 removed outlier: 3.623A pdb=" N TYR D 65 " --> pdb=" O ASP D 61 " (cutoff:3.500A) Processing helix chain 'D' and resid 88 through 95 Processing helix chain 'D' and resid 109 through 115 Processing helix chain 'D' and resid 129 through 146 removed outlier: 3.702A pdb=" N MET D 133 " --> pdb=" O GLN D 129 " (cutoff:3.500A) Processing helix chain 'D' and resid 147 through 165 removed outlier: 3.518A pdb=" N MET D 165 " --> pdb=" O PHE D 161 " (cutoff:3.500A) Processing helix chain 'D' and resid 166 through 170 removed outlier: 3.780A pdb=" N VAL D 170 " --> pdb=" O PRO D 167 " (cutoff:3.500A) Processing helix chain 'D' and resid 197 through 204 Processing helix chain 'D' and resid 225 through 230 Processing helix chain 'D' and resid 242 through 251 Processing helix chain 'D' and resid 252 through 257 Processing helix chain 'D' and resid 258 through 263 removed outlier: 3.513A pdb=" N ALA D 263 " --> pdb=" O PRO D 259 " (cutoff:3.500A) Processing helix chain 'D' and resid 277 through 293 Processing helix chain 'U' and resid 49 through 58 removed outlier: 3.947A pdb=" N LEU U 53 " --> pdb=" O GLY U 49 " (cutoff:3.500A) removed outlier: 3.968A pdb=" N LEU U 54 " --> pdb=" O LYS U 50 " (cutoff:3.500A) Processing helix chain 'U' and resid 79 through 86 Processing helix chain 'U' and resid 101 through 109 Processing helix chain 'U' and resid 110 through 115 removed outlier: 4.136A pdb=" N HIS U 114 " --> pdb=" O ARG U 110 " (cutoff:3.500A) Processing helix chain 'U' and resid 124 through 135 removed outlier: 3.725A pdb=" N VAL U 128 " --> pdb=" O ASP U 124 " (cutoff:3.500A) Processing helix chain 'U' and resid 149 through 164 Processing helix chain 'U' and resid 179 through 198 Processing helix chain 'U' and resid 207 through 214 Processing helix chain 'U' and resid 236 through 246 removed outlier: 4.196A pdb=" N GLU U 242 " --> pdb=" O GLY U 238 " (cutoff:3.500A) removed outlier: 4.409A pdb=" N MET U 243 " --> pdb=" O GLU U 239 " (cutoff:3.500A) Processing helix chain 'C' and resid 46 through 48 No H-bonds generated for 'chain 'C' and resid 46 through 48' Processing helix chain 'C' and resid 49 through 58 Processing helix chain 'C' and resid 79 through 84 Processing helix chain 'C' and resid 104 through 110 Processing helix chain 'C' and resid 110 through 115 Processing helix chain 'C' and resid 125 through 135 removed outlier: 3.508A pdb=" N GLU C 129 " --> pdb=" O ARG C 125 " (cutoff:3.500A) Processing helix chain 'C' and resid 149 through 165 removed outlier: 3.758A pdb=" N ARG C 155 " --> pdb=" O GLY C 151 " (cutoff:3.500A) removed outlier: 3.950A pdb=" N MET C 160 " --> pdb=" O ALA C 156 " (cutoff:3.500A) removed outlier: 3.747A pdb=" N LEU C 161 " --> pdb=" O TRP C 157 " (cutoff:3.500A) removed outlier: 3.613A pdb=" N ASP C 165 " --> pdb=" O LEU C 161 " (cutoff:3.500A) Processing helix chain 'C' and resid 179 through 197 Processing helix chain 'C' and resid 207 through 213 Processing helix chain 'C' and resid 231 through 236 Processing helix chain 'C' and resid 237 through 246 Processing sheet with id=AA1, first strand: chain 'D' and resid 79 through 80 removed outlier: 3.746A pdb=" N GLY D 58 " --> pdb=" O ILE D 38 " (cutoff:3.500A) removed outlier: 8.019A pdb=" N PHE D 122 " --> pdb=" O MET D 100 " (cutoff:3.500A) removed outlier: 6.899A pdb=" N TRP D 102 " --> pdb=" O PHE D 122 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'D' and resid 220 through 223 removed outlier: 6.710A pdb=" N HIS D 187 " --> pdb=" O ILE D 183 " (cutoff:3.500A) removed outlier: 3.760A pdb=" N ILE D 183 " --> pdb=" O HIS D 187 " (cutoff:3.500A) removed outlier: 6.531A pdb=" N LEU D 178 " --> pdb=" O LEU D 236 " (cutoff:3.500A) removed outlier: 7.951A pdb=" N PHE D 238 " --> pdb=" O LEU D 178 " (cutoff:3.500A) removed outlier: 7.806A pdb=" N THR D 180 " --> pdb=" O PHE D 238 " (cutoff:3.500A) removed outlier: 6.900A pdb=" N VAL D 235 " --> pdb=" O GLN D 267 " (cutoff:3.500A) removed outlier: 7.557A pdb=" N VAL D 269 " --> pdb=" O VAL D 235 " (cutoff:3.500A) removed outlier: 6.196A pdb=" N CYS D 237 " --> pdb=" O VAL D 269 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'U' and resid 30 through 33 removed outlier: 4.344A pdb=" N ALA U 13 " --> pdb=" O LEU U 69 " (cutoff:3.500A) removed outlier: 3.799A pdb=" N LEU U 69 " --> pdb=" O ALA U 13 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'U' and resid 87 through 89 removed outlier: 6.271A pdb=" N ALA U 88 " --> pdb=" O LEU U 170 " (cutoff:3.500A) removed outlier: 6.305A pdb=" N VAL U 39 " --> pdb=" O TYR U 217 " (cutoff:3.500A) removed outlier: 7.411A pdb=" N VAL U 219 " --> pdb=" O VAL U 39 " (cutoff:3.500A) removed outlier: 6.066A pdb=" N GLY U 41 " --> pdb=" O VAL U 219 " (cutoff:3.500A) removed outlier: 3.593A pdb=" N ALA U 228 " --> pdb=" O ALA U 220 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'U' and resid 250 through 252 Processing sheet with id=AA6, first strand: chain 'C' and resid 32 through 33 Processing sheet with id=AA7, first strand: chain 'C' and resid 87 through 90 removed outlier: 5.935A pdb=" N ALA C 88 " --> pdb=" O LEU C 170 " (cutoff:3.500A) removed outlier: 7.535A pdb=" N ASP C 172 " --> pdb=" O ALA C 88 " (cutoff:3.500A) removed outlier: 6.604A pdb=" N LEU C 90 " --> pdb=" O ASP C 172 " (cutoff:3.500A) removed outlier: 6.587A pdb=" N LEU C 169 " --> pdb=" O ILE C 202 " (cutoff:3.500A) removed outlier: 7.733A pdb=" N VAL C 204 " --> pdb=" O LEU C 169 " (cutoff:3.500A) removed outlier: 7.355A pdb=" N LEU C 171 " --> pdb=" O VAL C 204 " (cutoff:3.500A) removed outlier: 7.089A pdb=" N VAL C 219 " --> pdb=" O VAL C 39 " (cutoff:3.500A) removed outlier: 5.981A pdb=" N GLY C 41 " --> pdb=" O VAL C 219 " (cutoff:3.500A) removed outlier: 7.099A pdb=" N LEU C 221 " --> pdb=" O GLY C 41 " (cutoff:3.500A) removed outlier: 6.436A pdb=" N ILE C 43 " --> pdb=" O LEU C 221 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'C' and resid 249 through 252 598 hydrogen bonds defined for protein. 1749 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: 3.14 Time building geometry restraints manager: 4.42 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.34: 3367 1.34 - 1.46: 2164 1.46 - 1.58: 5023 1.58 - 1.69: 0 1.69 - 1.81: 117 Bond restraints: 10671 Sorted by residual: bond pdb=" CA LEU B 550 " pdb=" C LEU B 550 " ideal model delta sigma weight residual 1.525 1.571 -0.045 1.40e-02 5.10e+03 1.04e+01 bond pdb=" N PRO B 549 " pdb=" CD PRO B 549 " ideal model delta sigma weight residual 1.473 1.506 -0.033 1.40e-02 5.10e+03 5.47e+00 bond pdb=" N PRO B 549 " pdb=" CA PRO B 549 " ideal model delta sigma weight residual 1.467 1.483 -0.016 1.21e-02 6.83e+03 1.70e+00 bond pdb=" N LEU B 550 " pdb=" CA LEU B 550 " ideal model delta sigma weight residual 1.461 1.480 -0.018 1.49e-02 4.50e+03 1.49e+00 bond pdb=" CA ALA U 72 " pdb=" C ALA U 72 " ideal model delta sigma weight residual 1.530 1.515 0.015 1.30e-02 5.92e+03 1.32e+00 ... (remaining 10666 not shown) Histogram of bond angle deviations from ideal: 99.86 - 106.72: 421 106.72 - 113.58: 5989 113.58 - 120.44: 4119 120.44 - 127.30: 3873 127.30 - 134.16: 159 Bond angle restraints: 14561 Sorted by residual: angle pdb=" C PRO B 549 " pdb=" CA PRO B 549 " pdb=" CB PRO B 549 " ideal model delta sigma weight residual 111.21 100.99 10.22 1.33e+00 5.65e-01 5.91e+01 angle pdb=" N LEU B 550 " pdb=" CA LEU B 550 " pdb=" CB LEU B 550 " ideal model delta sigma weight residual 112.30 103.05 9.25 1.52e+00 4.33e-01 3.70e+01 angle pdb=" N LEU B 550 " pdb=" CA LEU B 550 " pdb=" C LEU B 550 " ideal model delta sigma weight residual 112.24 118.73 -6.49 1.28e+00 6.10e-01 2.57e+01 angle pdb=" N LEU B 548 " pdb=" CA LEU B 548 " pdb=" C LEU B 548 " ideal model delta sigma weight residual 113.16 106.36 6.80 1.42e+00 4.96e-01 2.29e+01 angle pdb=" N GLY B 551 " pdb=" CA GLY B 551 " pdb=" C GLY B 551 " ideal model delta sigma weight residual 113.18 121.39 -8.21 2.37e+00 1.78e-01 1.20e+01 ... (remaining 14556 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.97: 5621 17.97 - 35.95: 531 35.95 - 53.92: 124 53.92 - 71.90: 15 71.90 - 89.87: 8 Dihedral angle restraints: 6299 sinusoidal: 2319 harmonic: 3980 Sorted by residual: dihedral pdb=" C PRO B 549 " pdb=" N PRO B 549 " pdb=" CA PRO B 549 " pdb=" CB PRO B 549 " ideal model delta harmonic sigma weight residual -120.70 -107.88 -12.82 0 2.50e+00 1.60e-01 2.63e+01 dihedral pdb=" CB GLU U 65 " pdb=" CG GLU U 65 " pdb=" CD GLU U 65 " pdb=" OE1 GLU U 65 " ideal model delta sinusoidal sigma weight residual 0.00 89.87 -89.87 1 3.00e+01 1.11e-03 1.06e+01 dihedral pdb=" CA ASN U 16 " pdb=" CB ASN U 16 " pdb=" CG ASN U 16 " pdb=" OD1 ASN U 16 " ideal model delta sinusoidal sigma weight residual -90.00 -158.21 68.21 2 2.00e+01 2.50e-03 1.03e+01 ... (remaining 6296 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.043: 1327 0.043 - 0.086: 330 0.086 - 0.128: 81 0.128 - 0.171: 4 0.171 - 0.214: 1 Chirality restraints: 1743 Sorted by residual: chirality pdb=" CA PRO B 549 " pdb=" N PRO B 549 " pdb=" C PRO B 549 " pdb=" CB PRO B 549 " both_signs ideal model delta sigma weight residual False 2.72 2.93 -0.21 2.00e-01 2.50e+01 1.15e+00 chirality pdb=" CB ILE U 247 " pdb=" CA ILE U 247 " pdb=" CG1 ILE U 247 " pdb=" CG2 ILE U 247 " both_signs ideal model delta sigma weight residual False 2.64 2.48 0.16 2.00e-01 2.50e+01 6.67e-01 chirality pdb=" CA ILE C 251 " pdb=" N ILE C 251 " pdb=" C ILE C 251 " pdb=" CB ILE C 251 " both_signs ideal model delta sigma weight residual False 2.43 2.57 -0.14 2.00e-01 2.50e+01 4.94e-01 ... (remaining 1740 not shown) Planarity restraints: 1827 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C LEU B 548 " -0.043 5.00e-02 4.00e+02 6.30e-02 6.34e+00 pdb=" N PRO B 549 " 0.109 5.00e-02 4.00e+02 pdb=" CA PRO B 549 " -0.031 5.00e-02 4.00e+02 pdb=" CD PRO B 549 " -0.035 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA LEU B 550 " 0.010 2.00e-02 2.50e+03 1.99e-02 3.98e+00 pdb=" C LEU B 550 " -0.034 2.00e-02 2.50e+03 pdb=" O LEU B 550 " 0.014 2.00e-02 2.50e+03 pdb=" N GLY B 551 " 0.011 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C LEU D 44 " 0.028 5.00e-02 4.00e+02 4.20e-02 2.83e+00 pdb=" N PRO D 45 " -0.073 5.00e-02 4.00e+02 pdb=" CA PRO D 45 " 0.021 5.00e-02 4.00e+02 pdb=" CD PRO D 45 " 0.024 5.00e-02 4.00e+02 ... (remaining 1824 not shown) Histogram of nonbonded interaction distances: 2.09 - 2.65: 118 2.65 - 3.22: 9889 3.22 - 3.78: 16347 3.78 - 4.34: 21405 4.34 - 4.90: 35357 Nonbonded interactions: 83116 Sorted by model distance: nonbonded pdb=" O PRO B 549 " pdb=" CG LEU B 550 " model vdw 2.093 3.470 nonbonded pdb=" O ILE B 103 " pdb=" OG1 THR B 107 " model vdw 2.193 2.440 nonbonded pdb=" OG1 THR C 231 " pdb=" OE2 GLU C 234 " model vdw 2.207 2.440 nonbonded pdb=" OH TYR U 217 " pdb=" OE1 GLN U 229 " model vdw 2.215 2.440 nonbonded pdb=" O THR B 546 " pdb=" CG PRO B 549 " model vdw 2.235 3.440 ... (remaining 83111 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'C' and (resid 11 through 72 or (resid 73 through 75 and (name N or name \ CA or name C or name O or name CB )) or resid 76 through 137 or (resid 138 and ( \ name N or name CA or name C or name O or name CB )) or resid 139 or (resid 140 t \ hrough 141 and (name N or name CA or name C or name O or name CB )) or resid 142 \ through 265)) selection = (chain 'U' and (resid 11 through 19 or (resid 20 and (name N or name CA or name \ C or name O or name CB )) or resid 21 or (resid 22 and (name N or name CA or nam \ e C or name O or name CB )) or (resid 23 through 25 and (name N or name CA or na \ me C or name O or name CB )) or resid 26 or (resid 27 through 28 and (name N or \ name CA or name C or name O or name CB )) or resid 29 through 37 or (resid 38 an \ d (name N or name CA or name C or name O or name CB )) or resid 39 through 59 or \ (resid 60 and (name N or name CA or name C or name O or name CB )) or resid 61 \ through 93 or (resid 94 and (name N or name CA or name C or name O or name CB )) \ or resid 95 through 118 or (resid 119 through 120 and (name N or name CA or nam \ e C or name O or name CB )) or (resid 121 through 125 and (name N or name CA or \ name C or name O or name CB )) or resid 126 through 195 or (resid 196 and (name \ N or name CA or name C or name O or name CB )) or resid 197 through 247 or (resi \ d 248 and (name N or name CA or name C or name O or name CB )) or resid 249 thro \ ugh 259 or (resid 260 and (name N or name CA or name C or name O or name CB )) o \ r resid 261 through 265)) } 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=0.00 max=1.00 mean=0.98 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: 1.860 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: 8.100 Check model and map are aligned: 0.150 Set scattering table: 0.090 Process input model: 30.330 Find NCS groups from input model: 0.450 Set up NCS constraints: 0.050 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.800 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 43.840 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6978 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.002 0.045 10671 Z= 0.157 Angle : 0.620 10.221 14561 Z= 0.328 Chirality : 0.039 0.214 1743 Planarity : 0.004 0.063 1827 Dihedral : 14.950 89.874 3743 Min Nonbonded Distance : 2.093 Molprobity Statistics. All-atom Clashscore : 9.76 Ramachandran Plot: Outliers : 0.29 % Allowed : 4.73 % Favored : 94.98 % Rotamer Outliers : 2.66 % Cbeta Deviations : 0.08 % Peptide Plane: Cis-proline : 0.00 % 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: -1.98 (0.20), residues: 1395 helix: -0.21 (0.17), residues: 796 sheet: -2.06 (0.41), residues: 97 loop : -2.96 (0.22), residues: 502 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 249 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 221 time to evaluate : 1.228 Fit side-chains outliers start: 28 outliers final: 10 residues processed: 247 average time/residue: 0.2554 time to fit residues: 85.5316 Evaluate side-chains 184 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 174 time to evaluate : 1.181 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.1016 time to fit residues: 3.5947 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 117 optimal weight: 10.0000 chunk 105 optimal weight: 5.9990 chunk 58 optimal weight: 4.9990 chunk 35 optimal weight: 0.0370 chunk 71 optimal weight: 0.0570 chunk 56 optimal weight: 3.9990 chunk 108 optimal weight: 4.9990 chunk 42 optimal weight: 8.9990 chunk 66 optimal weight: 2.9990 chunk 81 optimal weight: 2.9990 chunk 126 optimal weight: 0.9980 overall best weight: 1.4180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 51 GLN B 55 HIS B 171 GLN U 16 ASN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7329 moved from start: 0.2441 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.052 10671 Z= 0.266 Angle : 0.658 12.894 14561 Z= 0.316 Chirality : 0.044 0.215 1743 Planarity : 0.005 0.048 1827 Dihedral : 4.260 35.700 1494 Min Nonbonded Distance : 2.134 Molprobity Statistics. All-atom Clashscore : 11.64 Ramachandran Plot: Outliers : 0.07 % Allowed : 4.09 % Favored : 95.84 % Rotamer Outliers : 3.70 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.08 % 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.91 (0.22), residues: 1395 helix: 0.67 (0.18), residues: 793 sheet: -1.76 (0.44), residues: 96 loop : -2.47 (0.23), residues: 506 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 233 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 39 poor density : 194 time to evaluate : 1.195 Fit side-chains outliers start: 39 outliers final: 23 residues processed: 224 average time/residue: 0.2211 time to fit residues: 70.9125 Evaluate side-chains 197 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 174 time to evaluate : 1.094 Switching outliers to nearest non-outliers outliers start: 23 outliers final: 0 residues processed: 23 average time/residue: 0.1016 time to fit residues: 6.0097 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 70 optimal weight: 0.9980 chunk 39 optimal weight: 0.0010 chunk 105 optimal weight: 9.9990 chunk 85 optimal weight: 0.9980 chunk 34 optimal weight: 5.9990 chunk 126 optimal weight: 4.9990 chunk 136 optimal weight: 0.9990 chunk 112 optimal weight: 0.9980 chunk 125 optimal weight: 2.9990 chunk 43 optimal weight: 0.6980 chunk 101 optimal weight: 3.9990 overall best weight: 0.7386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 51 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7316 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.003 0.039 10671 Z= 0.183 Angle : 0.607 15.069 14561 Z= 0.282 Chirality : 0.041 0.230 1743 Planarity : 0.005 0.046 1827 Dihedral : 3.992 31.546 1494 Min Nonbonded Distance : 2.141 Molprobity Statistics. All-atom Clashscore : 10.24 Ramachandran Plot: Outliers : 0.14 % Allowed : 4.44 % Favored : 95.41 % Rotamer Outliers : 2.56 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.08 % 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.37 (0.22), residues: 1395 helix: 1.14 (0.19), residues: 790 sheet: -1.71 (0.43), residues: 112 loop : -2.23 (0.25), residues: 493 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 213 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 186 time to evaluate : 1.224 Fit side-chains outliers start: 27 outliers final: 14 residues processed: 204 average time/residue: 0.1986 time to fit residues: 59.8948 Evaluate side-chains 191 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 177 time to evaluate : 1.197 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 0 residues processed: 14 average time/residue: 0.1377 time to fit residues: 4.7891 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 124 optimal weight: 0.9990 chunk 95 optimal weight: 0.3980 chunk 65 optimal weight: 1.9990 chunk 13 optimal weight: 10.0000 chunk 60 optimal weight: 0.6980 chunk 84 optimal weight: 0.4980 chunk 126 optimal weight: 4.9990 chunk 134 optimal weight: 7.9990 chunk 66 optimal weight: 0.9980 chunk 120 optimal weight: 0.0980 chunk 36 optimal weight: 0.9980 overall best weight: 0.5380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 51 GLN B 81 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7310 moved from start: 0.2972 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.051 10671 Z= 0.164 Angle : 0.586 10.575 14561 Z= 0.273 Chirality : 0.041 0.206 1743 Planarity : 0.004 0.048 1827 Dihedral : 3.935 30.119 1494 Min Nonbonded Distance : 2.130 Molprobity Statistics. All-atom Clashscore : 10.58 Ramachandran Plot: Outliers : 0.14 % Allowed : 4.01 % Favored : 95.84 % Rotamer Outliers : 2.09 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.08 % 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.23), residues: 1395 helix: 1.40 (0.19), residues: 790 sheet: -1.55 (0.45), residues: 112 loop : -2.09 (0.25), residues: 493 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 202 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 22 poor density : 180 time to evaluate : 1.157 Fit side-chains outliers start: 22 outliers final: 10 residues processed: 197 average time/residue: 0.1962 time to fit residues: 57.6305 Evaluate side-chains 187 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 177 time to evaluate : 1.152 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.1025 time to fit residues: 3.4810 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 111 optimal weight: 0.9980 chunk 76 optimal weight: 3.9990 chunk 1 optimal weight: 5.9990 chunk 100 optimal weight: 7.9990 chunk 55 optimal weight: 0.8980 chunk 114 optimal weight: 10.0000 chunk 92 optimal weight: 8.9990 chunk 0 optimal weight: 10.0000 chunk 68 optimal weight: 1.9990 chunk 120 optimal weight: 0.4980 chunk 33 optimal weight: 0.6980 overall best weight: 1.0182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: D 187 HIS C 195 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7394 moved from start: 0.3435 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.045 10671 Z= 0.213 Angle : 0.610 11.491 14561 Z= 0.286 Chirality : 0.041 0.148 1743 Planarity : 0.005 0.047 1827 Dihedral : 4.049 28.970 1494 Min Nonbonded Distance : 2.150 Molprobity Statistics. All-atom Clashscore : 11.21 Ramachandran Plot: Outliers : 0.14 % Allowed : 4.01 % Favored : 95.84 % Rotamer Outliers : 2.47 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.23), residues: 1395 helix: 1.38 (0.19), residues: 791 sheet: -1.34 (0.46), residues: 114 loop : -2.05 (0.25), residues: 490 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 209 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 26 poor density : 183 time to evaluate : 1.170 Fit side-chains revert: symmetry clash outliers start: 26 outliers final: 15 residues processed: 197 average time/residue: 0.2047 time to fit residues: 59.2060 Evaluate side-chains 191 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 176 time to evaluate : 1.262 Switching outliers to nearest non-outliers outliers start: 15 outliers final: 0 residues processed: 15 average time/residue: 0.1213 time to fit residues: 4.8394 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 45 optimal weight: 7.9990 chunk 120 optimal weight: 5.9990 chunk 26 optimal weight: 8.9990 chunk 78 optimal weight: 3.9990 chunk 33 optimal weight: 0.7980 chunk 134 optimal weight: 3.9990 chunk 111 optimal weight: 0.0000 chunk 62 optimal weight: 1.9990 chunk 11 optimal weight: 9.9990 chunk 44 optimal weight: 0.8980 chunk 70 optimal weight: 0.0370 overall best weight: 0.7464 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7373 moved from start: 0.3627 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.040 10671 Z= 0.181 Angle : 0.612 11.978 14561 Z= 0.283 Chirality : 0.042 0.295 1743 Planarity : 0.005 0.048 1827 Dihedral : 4.026 28.183 1494 Min Nonbonded Distance : 2.150 Molprobity Statistics. All-atom Clashscore : 10.72 Ramachandran Plot: Outliers : 0.14 % Allowed : 3.73 % Favored : 96.13 % Rotamer Outliers : 1.80 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.08 (0.23), residues: 1395 helix: 1.43 (0.19), residues: 793 sheet: -1.25 (0.45), residues: 118 loop : -1.98 (0.26), residues: 484 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 201 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 19 poor density : 182 time to evaluate : 1.252 Fit side-chains outliers start: 19 outliers final: 4 residues processed: 197 average time/residue: 0.2072 time to fit residues: 59.9960 Evaluate side-chains 179 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 175 time to evaluate : 1.052 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.1053 time to fit residues: 2.3078 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 129 optimal weight: 9.9990 chunk 15 optimal weight: 2.9990 chunk 76 optimal weight: 3.9990 chunk 98 optimal weight: 1.9990 chunk 113 optimal weight: 30.0000 chunk 75 optimal weight: 0.7980 chunk 133 optimal weight: 9.9990 chunk 83 optimal weight: 0.8980 chunk 81 optimal weight: 5.9990 chunk 61 optimal weight: 0.9990 chunk 82 optimal weight: 0.9980 overall best weight: 1.1384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7429 moved from start: 0.3939 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.043 10671 Z= 0.224 Angle : 0.639 11.931 14561 Z= 0.298 Chirality : 0.042 0.169 1743 Planarity : 0.005 0.048 1827 Dihedral : 4.164 27.427 1494 Min Nonbonded Distance : 2.156 Molprobity Statistics. All-atom Clashscore : 10.92 Ramachandran Plot: Outliers : 0.14 % Allowed : 4.80 % Favored : 95.05 % Rotamer Outliers : 1.71 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.23), residues: 1395 helix: 1.38 (0.19), residues: 793 sheet: -1.10 (0.44), residues: 120 loop : -2.06 (0.26), residues: 482 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 195 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 177 time to evaluate : 1.114 Fit side-chains revert: symmetry clash outliers start: 18 outliers final: 9 residues processed: 188 average time/residue: 0.2044 time to fit residues: 56.4444 Evaluate side-chains 180 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 171 time to evaluate : 1.165 Switching outliers to nearest non-outliers outliers start: 9 outliers final: 0 residues processed: 9 average time/residue: 0.1128 time to fit residues: 3.4065 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 53 optimal weight: 1.9990 chunk 80 optimal weight: 4.9990 chunk 40 optimal weight: 0.8980 chunk 26 optimal weight: 4.9990 chunk 25 optimal weight: 2.9990 chunk 85 optimal weight: 0.5980 chunk 91 optimal weight: 0.9990 chunk 66 optimal weight: 1.9990 chunk 12 optimal weight: 5.9990 chunk 105 optimal weight: 2.9990 chunk 121 optimal weight: 5.9990 overall best weight: 1.2986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7448 moved from start: 0.4203 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.046 10671 Z= 0.242 Angle : 0.651 13.638 14561 Z= 0.306 Chirality : 0.043 0.194 1743 Planarity : 0.005 0.048 1827 Dihedral : 4.284 27.266 1494 Min Nonbonded Distance : 2.188 Molprobity Statistics. All-atom Clashscore : 11.16 Ramachandran Plot: Outliers : 0.14 % Allowed : 4.52 % Favored : 95.34 % Rotamer Outliers : 1.33 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.22), residues: 1395 helix: 1.27 (0.19), residues: 791 sheet: -1.04 (0.44), residues: 119 loop : -2.08 (0.26), residues: 485 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 195 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 181 time to evaluate : 1.239 Fit side-chains outliers start: 14 outliers final: 6 residues processed: 189 average time/residue: 0.1993 time to fit residues: 56.1956 Evaluate side-chains 179 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 173 time to evaluate : 1.208 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.1009 time to fit residues: 2.8262 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 128 optimal weight: 0.9980 chunk 117 optimal weight: 10.0000 chunk 124 optimal weight: 2.9990 chunk 75 optimal weight: 3.9990 chunk 54 optimal weight: 0.7980 chunk 98 optimal weight: 1.9990 chunk 38 optimal weight: 0.9990 chunk 112 optimal weight: 0.8980 chunk 118 optimal weight: 1.9990 chunk 81 optimal weight: 4.9990 chunk 132 optimal weight: 9.9990 overall best weight: 1.1384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 164 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7446 moved from start: 0.4328 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.044 10671 Z= 0.227 Angle : 0.666 14.979 14561 Z= 0.308 Chirality : 0.043 0.204 1743 Planarity : 0.005 0.049 1827 Dihedral : 4.263 26.702 1494 Min Nonbonded Distance : 2.188 Molprobity Statistics. All-atom Clashscore : 11.40 Ramachandran Plot: Outliers : 0.14 % Allowed : 5.16 % Favored : 94.70 % Rotamer Outliers : 0.38 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.23), residues: 1395 helix: 1.30 (0.19), residues: 795 sheet: -1.07 (0.44), residues: 119 loop : -2.08 (0.26), residues: 481 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 175 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 171 time to evaluate : 1.202 Fit side-chains outliers start: 4 outliers final: 0 residues processed: 175 average time/residue: 0.2027 time to fit residues: 52.9156 Evaluate side-chains 165 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 165 time to evaluate : 1.245 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.7396 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 80 optimal weight: 0.5980 chunk 62 optimal weight: 0.2980 chunk 91 optimal weight: 2.9990 chunk 138 optimal weight: 0.9990 chunk 127 optimal weight: 0.0770 chunk 110 optimal weight: 7.9990 chunk 11 optimal weight: 9.9990 chunk 85 optimal weight: 0.8980 chunk 67 optimal weight: 0.7980 chunk 87 optimal weight: 0.5980 chunk 117 optimal weight: 10.0000 overall best weight: 0.4738 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 183 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7383 moved from start: 0.4313 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.038 10671 Z= 0.165 Angle : 0.632 13.727 14561 Z= 0.291 Chirality : 0.041 0.195 1743 Planarity : 0.005 0.049 1827 Dihedral : 4.087 25.002 1494 Min Nonbonded Distance : 2.173 Molprobity Statistics. All-atom Clashscore : 10.24 Ramachandran Plot: Outliers : 0.14 % Allowed : 4.44 % Favored : 95.41 % Rotamer Outliers : 0.28 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.21 (0.23), residues: 1395 helix: 1.52 (0.19), residues: 794 sheet: -0.90 (0.45), residues: 119 loop : -2.01 (0.26), residues: 482 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2790 Ramachandran restraints generated. 1395 Oldfield, 0 Emsley, 1395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 185 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 182 time to evaluate : 1.287 Fit side-chains outliers start: 3 outliers final: 2 residues processed: 183 average time/residue: 0.2229 time to fit residues: 59.4340 Evaluate side-chains 177 residues out of total 1114 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 175 time to evaluate : 1.239 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0989 time to fit residues: 2.0017 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 139 random chunks: chunk 33 optimal weight: 1.9990 chunk 101 optimal weight: 6.9990 chunk 16 optimal weight: 1.9990 chunk 30 optimal weight: 4.9990 chunk 110 optimal weight: 7.9990 chunk 46 optimal weight: 0.8980 chunk 113 optimal weight: 7.9990 chunk 13 optimal weight: 6.9990 chunk 20 optimal weight: 3.9990 chunk 96 optimal weight: 5.9990 chunk 6 optimal weight: 0.9980 overall best weight: 1.9786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4021 r_free = 0.4021 target = 0.187459 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 54)----------------| | r_work = 0.3803 r_free = 0.3803 target = 0.165804 restraints weight = 11817.915| |-----------------------------------------------------------------------------| r_work (start): 0.3807 rms_B_bonded: 0.92 r_work: 0.3622 rms_B_bonded: 1.78 restraints_weight: 0.5000 r_work: 0.3446 rms_B_bonded: 4.59 restraints_weight: 0.2500 r_work (final): 0.3446 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7875 moved from start: 0.4628 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.055 10671 Z= 0.329 Angle : 0.735 14.304 14561 Z= 0.344 Chirality : 0.045 0.189 1743 Planarity : 0.005 0.048 1827 Dihedral : 4.533 29.124 1494 Min Nonbonded Distance : 2.176 Molprobity Statistics. All-atom Clashscore : 12.46 Ramachandran Plot: Outliers : 0.14 % Allowed : 5.30 % Favored : 94.55 % Rotamer Outliers : 0.28 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % 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.21 (0.22), residues: 1395 helix: 1.14 (0.19), residues: 790 sheet: -1.07 (0.44), residues: 119 loop : -2.08 (0.26), residues: 486 =============================================================================== Job complete usr+sys time: 2258.40 seconds wall clock time: 41 minutes 30.70 seconds (2490.70 seconds total)