Starting phenix.real_space_refine (version: 1.20rc4) on Wed Dec 1 22:27:53 2021 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838.map Found restraint, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838_ligands.cif Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.5 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { restraint_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838_ligands.cif" default_restraint = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838_ligands.cif" real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mi3_23838/12_2021/7mi3_23838_updated.pdb" } resolution = 3.5 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.001 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Construct map_model_manager Extract box with map and model Check model and map are aligned Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "A PHE 2368": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 2650": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 3364": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 3389": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 3457": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 3470": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.07s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.20rc4-4434/modules/chem_data/mon_lib" Total number of atoms: 9682 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 9682 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1196, 9682 Unusual residues: {' MG': 1, 'ATP': 1, 'ZG7': 2} Classifications: {'undetermined': 4, 'peptide': 1192} Link IDs: {'PTRANS': 39, 'TRANS': 1152, None: 4} Not linked: pdbres="GLU A3523 " pdbres="ZG7 A4201 " Not linked: pdbres="ZG7 A4201 " pdbres="ZG7 A4202 " Not linked: pdbres="ZG7 A4202 " pdbres="ATP A4203 " Not linked: pdbres="ATP A4203 " pdbres=" MG A4204 " Chain breaks: 1 Time building chain proxies: 5.78, per 1000 atoms: 0.60 Number of scatterers: 9682 At special positions: 0 Unit cell: (89.096, 137.27, 103.082, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 8 Type Number sf(0) Br 2 34.99 Cl 2 17.00 S 49 16.00 P 3 15.00 Mg 1 11.99 O 1804 8.00 N 1624 7.00 C 6197 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 Amimo 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.22 Conformation dependent library (CDL) restraints added in 1.5 seconds 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2272 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 50 helices and 8 sheets defined 51.5% alpha, 10.7% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.80 Creating SS restraints... Processing helix chain 'A' and resid 2032 through 2045 Processing helix chain 'A' and resid 2051 through 2065 removed outlier: 3.928A pdb=" N LYS A2065 " --> pdb=" O TYR A2061 " (cutoff:3.500A) Processing helix chain 'A' and resid 2080 through 2095 Processing helix chain 'A' and resid 2106 through 2108 No H-bonds generated for 'chain 'A' and resid 2106 through 2108' Processing helix chain 'A' and resid 2111 through 2114 No H-bonds generated for 'chain 'A' and resid 2111 through 2114' Processing helix chain 'A' and resid 2129 through 2139 Processing helix chain 'A' and resid 2160 through 2170 removed outlier: 3.564A pdb=" N MET A2166 " --> pdb=" O TYR A2162 " (cutoff:3.500A) removed outlier: 5.130A pdb=" N SER A2168 " --> pdb=" O GLU A2164 " (cutoff:3.500A) removed outlier: 6.918A pdb=" N VAL A2169 " --> pdb=" O ALA A2165 " (cutoff:3.500A) removed outlier: 3.881A pdb=" N LEU A2170 " --> pdb=" O MET A2166 " (cutoff:3.500A) Processing helix chain 'A' and resid 2204 through 2209 removed outlier: 3.689A pdb=" N ARG A2209 " --> pdb=" O ALA A2205 " (cutoff:3.500A) Processing helix chain 'A' and resid 2222 through 2240 removed outlier: 3.803A pdb=" N ALA A2236 " --> pdb=" O LYS A2232 " (cutoff:3.500A) removed outlier: 4.504A pdb=" N ASN A2239 " --> pdb=" O GLU A2235 " (cutoff:3.500A) removed outlier: 4.380A pdb=" N LYS A2240 " --> pdb=" O ALA A2236 " (cutoff:3.500A) Processing helix chain 'A' and resid 2243 through 2257 removed outlier: 4.074A pdb=" N LYS A2248 " --> pdb=" O PHE A2244 " (cutoff:3.500A) removed outlier: 3.545A pdb=" N ASP A2255 " --> pdb=" O ASP A2251 " (cutoff:3.500A) removed outlier: 4.648A pdb=" N SER A2256 " --> pdb=" O LEU A2252 " (cutoff:3.500A) Processing helix chain 'A' and resid 2260 through 2267 Processing helix chain 'A' and resid 2280 through 2282 No H-bonds generated for 'chain 'A' and resid 2280 through 2282' Processing helix chain 'A' and resid 2284 through 2304 removed outlier: 3.805A pdb=" N TYR A2298 " --> pdb=" O LEU A2294 " (cutoff:3.500A) removed outlier: 3.997A pdb=" N ARG A2299 " --> pdb=" O ILE A2295 " (cutoff:3.500A) removed outlier: 5.811A pdb=" N GLN A2300 " --> pdb=" O SER A2296 " (cutoff:3.500A) removed outlier: 6.377A pdb=" N TRP A2301 " --> pdb=" O SER A2297 " (cutoff:3.500A) removed outlier: 4.231A pdb=" N GLN A2303 " --> pdb=" O ARG A2299 " (cutoff:3.500A) removed outlier: 6.186A pdb=" N ASN A2304 " --> pdb=" O GLN A2300 " (cutoff:3.500A) Processing helix chain 'A' and resid 2307 through 2327 Processing helix chain 'A' and resid 2333 through 2345 removed outlier: 3.857A pdb=" N GLN A2340 " --> pdb=" O ARG A2336 " (cutoff:3.500A) removed outlier: 3.702A pdb=" N THR A2341 " --> pdb=" O ALA A2337 " (cutoff:3.500A) Processing helix chain 'A' and resid 2382 through 2386 Processing helix chain 'A' and resid 2395 through 2409 Processing helix chain 'A' and resid 2424 through 2434 removed outlier: 3.842A pdb=" N MET A2428 " --> pdb=" O LYS A2424 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N ASN A2434 " --> pdb=" O ASN A2430 " (cutoff:3.500A) Processing helix chain 'A' and resid 2451 through 2461 removed outlier: 4.234A pdb=" N ARG A2460 " --> pdb=" O SER A2456 " (cutoff:3.500A) removed outlier: 4.552A pdb=" N HIS A2461 " --> pdb=" O ALA A2457 " (cutoff:3.500A) Processing helix chain 'A' and resid 2501 through 2512 Processing helix chain 'A' and resid 2548 through 2553 removed outlier: 4.521A pdb=" N HIS A2553 " --> pdb=" O ARG A2549 " (cutoff:3.500A) Processing helix chain 'A' and resid 2565 through 2581 removed outlier: 4.569A pdb=" N LEU A2581 " --> pdb=" O ALA A2577 " (cutoff:3.500A) Processing helix chain 'A' and resid 2583 through 2606 removed outlier: 4.799A pdb=" N SER A2587 " --> pdb=" O GLU A2584 " (cutoff:3.500A) removed outlier: 3.837A pdb=" N GLU A2590 " --> pdb=" O SER A2587 " (cutoff:3.500A) Proline residue: A2591 - end of helix Processing helix chain 'A' and resid 2619 through 2634 removed outlier: 3.525A pdb=" N ARG A2624 " --> pdb=" O GLU A2621 " (cutoff:3.500A) removed outlier: 3.505A pdb=" N VAL A2629 " --> pdb=" O VAL A2626 " (cutoff:3.500A) Processing helix chain 'A' and resid 2641 through 2659 removed outlier: 4.906A pdb=" N ASP A2658 " --> pdb=" O ARG A2654 " (cutoff:3.500A) removed outlier: 5.215A pdb=" N ARG A2659 " --> pdb=" O ILE A2655 " (cutoff:3.500A) Processing helix chain 'A' and resid 2663 through 2680 Processing helix chain 'A' and resid 2709 through 2726 Processing helix chain 'A' and resid 2736 through 2750 Processing helix chain 'A' and resid 2766 through 2777 removed outlier: 3.526A pdb=" N ARG A2771 " --> pdb=" O THR A2767 " (cutoff:3.500A) Processing helix chain 'A' and resid 2792 through 2808 Processing helix chain 'A' and resid 2819 through 2821 No H-bonds generated for 'chain 'A' and resid 2819 through 2821' Processing helix chain 'A' and resid 2825 through 2836 Processing helix chain 'A' and resid 2841 through 2844 Processing helix chain 'A' and resid 2846 through 2863 removed outlier: 3.557A pdb=" N LYS A2851 " --> pdb=" O GLU A2848 " (cutoff:3.500A) removed outlier: 3.520A pdb=" N THR A2860 " --> pdb=" O ARG A2857 " (cutoff:3.500A) Processing helix chain 'A' and resid 2870 through 2884 Processing helix chain 'A' and resid 2897 through 2904 Processing helix chain 'A' and resid 2906 through 2909 No H-bonds generated for 'chain 'A' and resid 2906 through 2909' Processing helix chain 'A' and resid 2922 through 2935 removed outlier: 3.994A pdb=" N VAL A2935 " --> pdb=" O ASN A2931 " (cutoff:3.500A) Processing helix chain 'A' and resid 2942 through 2944 No H-bonds generated for 'chain 'A' and resid 2942 through 2944' Processing helix chain 'A' and resid 2961 through 2979 Processing helix chain 'A' and resid 2990 through 3010 removed outlier: 3.600A pdb=" N VAL A3000 " --> pdb=" O LEU A2996 " (cutoff:3.500A) removed outlier: 3.637A pdb=" N LYS A3001 " --> pdb=" O ARG A2997 " (cutoff:3.500A) Processing helix chain 'A' and resid 3317 through 3336 removed outlier: 3.789A pdb=" N ILE A3321 " --> pdb=" O GLN A3318 " (cutoff:3.500A) removed outlier: 3.991A pdb=" N TYR A3333 " --> pdb=" O TYR A3330 " (cutoff:3.500A) removed outlier: 4.988A pdb=" N PHE A3334 " --> pdb=" O GLU A3331 " (cutoff:3.500A) removed outlier: 5.221A pdb=" N GLY A3335 " --> pdb=" O THR A3332 " (cutoff:3.500A) removed outlier: 4.313A pdb=" N HIS A3336 " --> pdb=" O TYR A3333 " (cutoff:3.500A) Processing helix chain 'A' and resid 3339 through 3355 removed outlier: 3.746A pdb=" N LEU A3352 " --> pdb=" O ILE A3348 " (cutoff:3.500A) Processing helix chain 'A' and resid 3366 through 3370 Processing helix chain 'A' and resid 3373 through 3381 removed outlier: 4.105A pdb=" N LEU A3380 " --> pdb=" O LYS A3376 " (cutoff:3.500A) Processing helix chain 'A' and resid 3388 through 3399 Processing helix chain 'A' and resid 3414 through 3425 removed outlier: 4.237A pdb=" N SER A3419 " --> pdb=" O ILE A3415 " (cutoff:3.500A) removed outlier: 7.223A pdb=" N ASN A3424 " --> pdb=" O ASN A3420 " (cutoff:3.500A) removed outlier: 7.969A pdb=" N LYS A3425 " --> pdb=" O TYR A3421 " (cutoff:3.500A) Processing helix chain 'A' and resid 3436 through 3446 Processing helix chain 'A' and resid 3460 through 3466 removed outlier: 4.394A pdb=" N ARG A3464 " --> pdb=" O ILE A3461 " (cutoff:3.500A) Processing helix chain 'A' and resid 3507 through 3510 No H-bonds generated for 'chain 'A' and resid 3507 through 3510' Processing sheet with id= A, first strand: chain 'A' and resid 2211 through 2215 removed outlier: 8.827A pdb=" N LEU A2070 " --> pdb=" O PHE A2190 " (cutoff:3.500A) removed outlier: 6.469A pdb=" N ILE A2192 " --> pdb=" O LEU A2070 " (cutoff:3.500A) removed outlier: 7.784A pdb=" N LEU A2072 " --> pdb=" O ILE A2192 " (cutoff:3.500A) removed outlier: 6.532A pdb=" N PHE A2194 " --> pdb=" O LEU A2072 " (cutoff:3.500A) removed outlier: 7.863A pdb=" N GLY A2074 " --> pdb=" O PHE A2194 " (cutoff:3.500A) removed outlier: 6.558A pdb=" N THR A2196 " --> pdb=" O GLY A2074 " (cutoff:3.500A) removed outlier: 8.204A pdb=" N ARG A2191 " --> pdb=" O SER A2148 " (cutoff:3.500A) removed outlier: 6.318A pdb=" N ILE A2150 " --> pdb=" O ARG A2191 " (cutoff:3.500A) removed outlier: 7.497A pdb=" N LEU A2193 " --> pdb=" O ILE A2150 " (cutoff:3.500A) removed outlier: 6.671A pdb=" N VAL A2152 " --> pdb=" O LEU A2193 " (cutoff:3.500A) removed outlier: 7.357A pdb=" N GLU A2195 " --> pdb=" O VAL A2152 " (cutoff:3.500A) removed outlier: 6.547A pdb=" N ALA A2098 " --> pdb=" O ARG A2149 " (cutoff:3.500A) removed outlier: 7.748A pdb=" N TRP A2151 " --> pdb=" O ALA A2098 " (cutoff:3.500A) removed outlier: 6.696A pdb=" N VAL A2100 " --> pdb=" O TRP A2151 " (cutoff:3.500A) removed outlier: 7.698A pdb=" N VAL A2153 " --> pdb=" O VAL A2100 " (cutoff:3.500A) removed outlier: 6.652A pdb=" N TYR A2102 " --> pdb=" O VAL A2153 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 2175 through 2177 Processing sheet with id= C, first strand: chain 'A' and resid 2358 through 2361 removed outlier: 3.589A pdb=" N THR A2358 " --> pdb=" O PHE A2371 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'A' and resid 2554 through 2559 removed outlier: 6.762A pdb=" N GLY A2413 " --> pdb=" O ALA A2555 " (cutoff:3.500A) removed outlier: 7.703A pdb=" N LEU A2557 " --> pdb=" O GLY A2413 " (cutoff:3.500A) removed outlier: 6.260A pdb=" N ILE A2415 " --> pdb=" O LEU A2557 " (cutoff:3.500A) removed outlier: 7.293A pdb=" N LEU A2559 " --> pdb=" O ILE A2415 " (cutoff:3.500A) removed outlier: 6.901A pdb=" N CYS A2417 " --> pdb=" O LEU A2559 " (cutoff:3.500A) removed outlier: 3.637A pdb=" N ILE A2414 " --> pdb=" O ILE A2531 " (cutoff:3.500A) removed outlier: 6.640A pdb=" N LEU A2482 " --> pdb=" O HIS A2530 " (cutoff:3.500A) removed outlier: 7.496A pdb=" N VAL A2532 " --> pdb=" O LEU A2482 " (cutoff:3.500A) removed outlier: 5.944A pdb=" N LEU A2484 " --> pdb=" O VAL A2532 " (cutoff:3.500A) removed outlier: 7.224A pdb=" N ALA A2534 " --> pdb=" O LEU A2484 " (cutoff:3.500A) removed outlier: 6.778A pdb=" N CYS A2486 " --> pdb=" O ALA A2534 " (cutoff:3.500A) Processing sheet with id= E, first strand: chain 'A' and resid 2465 through 2467 Processing sheet with id= F, first strand: chain 'A' and resid 2780 through 2783 removed outlier: 6.863A pdb=" N MET A2756 " --> pdb=" O PHE A2889 " (cutoff:3.500A) removed outlier: 8.070A pdb=" N ILE A2891 " --> pdb=" O MET A2756 " (cutoff:3.500A) removed outlier: 6.622A pdb=" N LEU A2758 " --> pdb=" O ILE A2891 " (cutoff:3.500A) Processing sheet with id= G, first strand: chain 'A' and resid 3514 through 3517 removed outlier: 7.538A pdb=" N LEU A3494 " --> pdb=" O PRO A3405 " (cutoff:3.500A) removed outlier: 8.458A pdb=" N LEU A3407 " --> pdb=" O LEU A3494 " (cutoff:3.500A) removed outlier: 6.943A pdb=" N ILE A3496 " --> pdb=" O LEU A3407 " (cutoff:3.500A) Processing sheet with id= H, first strand: chain 'A' and resid 3478 through 3481 removed outlier: 3.591A pdb=" N HIS A3484 " --> pdb=" O ILE A3481 " (cutoff:3.500A) 432 hydrogen bonds defined for protein. 1110 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.55 Time building geometry restraints manager: 4.52 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.13 - 1.29: 1587 1.29 - 1.46: 3126 1.46 - 1.62: 5093 1.62 - 1.78: 18 1.78 - 1.94: 65 Bond restraints: 9889 Sorted by residual: bond pdb=" C26 ZG7 A4202 " pdb=" O27 ZG7 A4202 " ideal model delta sigma weight residual 1.407 1.166 0.241 2.00e-02 2.50e+03 1.45e+02 bond pdb=" C26 ZG7 A4201 " pdb=" O27 ZG7 A4201 " ideal model delta sigma weight residual 1.407 1.187 0.220 2.00e-02 2.50e+03 1.21e+02 bond pdb=" C14 ZG7 A4202 " pdb=" N15 ZG7 A4202 " ideal model delta sigma weight residual 1.376 1.531 -0.155 2.00e-02 2.50e+03 6.02e+01 bond pdb=" C14 ZG7 A4201 " pdb=" N15 ZG7 A4201 " ideal model delta sigma weight residual 1.376 1.520 -0.144 2.00e-02 2.50e+03 5.18e+01 bond pdb=" C5 ATP A4203 " pdb=" C6 ATP A4203 " ideal model delta sigma weight residual 1.409 1.479 -0.070 1.00e-02 1.00e+04 4.87e+01 ... (remaining 9884 not shown) Histogram of bond angle deviations from ideal: 59.63 - 83.58: 6 83.58 - 107.53: 456 107.53 - 131.48: 12908 131.48 - 155.43: 31 155.43 - 179.38: 2 Bond angle restraints: 13403 Sorted by residual: angle pdb=" PA ATP A4203 " pdb=" O3A ATP A4203 " pdb=" PB ATP A4203 " ideal model delta sigma weight residual 136.83 121.50 15.33 1.00e+00 1.00e+00 2.35e+02 angle pdb=" PB ATP A4203 " pdb=" O3B ATP A4203 " pdb=" PG ATP A4203 " ideal model delta sigma weight residual 139.87 128.63 11.24 1.00e+00 1.00e+00 1.26e+02 angle pdb=" C05 ZG7 A4201 " pdb=" C08 ZG7 A4201 " pdb=" C11 ZG7 A4201 " ideal model delta sigma weight residual 121.66 95.55 26.11 3.00e+00 1.11e-01 7.57e+01 angle pdb=" C5 ATP A4203 " pdb=" C4 ATP A4203 " pdb=" N3 ATP A4203 " ideal model delta sigma weight residual 126.80 118.41 8.39 1.00e+00 1.00e+00 7.04e+01 angle pdb=" C05 ZG7 A4202 " pdb=" C08 ZG7 A4202 " pdb=" C11 ZG7 A4202 " ideal model delta sigma weight residual 121.66 99.96 21.70 3.00e+00 1.11e-01 5.23e+01 ... (remaining 13398 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.88: 5530 17.88 - 35.76: 302 35.76 - 53.64: 46 53.64 - 71.53: 24 71.53 - 89.41: 11 Dihedral angle restraints: 5913 sinusoidal: 2411 harmonic: 3502 Sorted by residual: dihedral pdb=" CA ASP A2105 " pdb=" C ASP A2105 " pdb=" N THR A2106 " pdb=" CA THR A2106 " ideal model delta harmonic sigma weight residual -180.00 -153.51 -26.49 0 5.00e+00 4.00e-02 2.81e+01 dihedral pdb=" CA GLN A2639 " pdb=" C GLN A2639 " pdb=" N THR A2640 " pdb=" CA THR A2640 " ideal model delta harmonic sigma weight residual 180.00 154.29 25.71 0 5.00e+00 4.00e-02 2.64e+01 dihedral pdb=" CA GLU A2810 " pdb=" C GLU A2810 " pdb=" N SER A2811 " pdb=" CA SER A2811 " ideal model delta harmonic sigma weight residual -180.00 -156.77 -23.23 0 5.00e+00 4.00e-02 2.16e+01 ... (remaining 5910 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.085: 1157 0.085 - 0.169: 315 0.169 - 0.254: 37 0.254 - 0.338: 7 0.338 - 0.423: 4 Chirality restraints: 1520 Sorted by residual: chirality pdb=" CB ILE A2226 " pdb=" CA ILE A2226 " pdb=" CG1 ILE A2226 " pdb=" CG2 ILE A2226 " both_signs ideal model delta sigma weight residual False 2.64 2.22 0.42 2.00e-01 2.50e+01 4.47e+00 chirality pdb=" CB ILE A2141 " pdb=" CA ILE A2141 " pdb=" CG1 ILE A2141 " pdb=" CG2 ILE A2141 " both_signs ideal model delta sigma weight residual False 2.64 3.00 -0.35 2.00e-01 2.50e+01 3.09e+00 chirality pdb=" CA LEU A2366 " pdb=" N LEU A2366 " pdb=" C LEU A2366 " pdb=" CB LEU A2366 " both_signs ideal model delta sigma weight residual False 2.51 2.16 0.35 2.00e-01 2.50e+01 3.09e+00 ... (remaining 1517 not shown) Planarity restraints: 1686 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C08 ZG7 A4202 " -0.014 2.00e-02 2.50e+03 2.73e-02 2.61e+01 pdb=" C11 ZG7 A4202 " 0.067 2.00e-02 2.50e+03 pdb=" C14 ZG7 A4202 " -0.009 2.00e-02 2.50e+03 pdb=" C16 ZG7 A4202 " -0.015 2.00e-02 2.50e+03 pdb=" C17 ZG7 A4202 " -0.012 2.00e-02 2.50e+03 pdb=" C26 ZG7 A4202 " 0.002 2.00e-02 2.50e+03 pdb=" C28 ZG7 A4202 " -0.018 2.00e-02 2.50e+03 pdb=" C29 ZG7 A4202 " 0.003 2.00e-02 2.50e+03 pdb=" N12 ZG7 A4202 " -0.005 2.00e-02 2.50e+03 pdb=" N13 ZG7 A4202 " -0.029 2.00e-02 2.50e+03 pdb=" N15 ZG7 A4202 " -0.009 2.00e-02 2.50e+03 pdb=" N30 ZG7 A4202 " 0.010 2.00e-02 2.50e+03 pdb=" O18 ZG7 A4202 " 0.056 2.00e-02 2.50e+03 pdb=" O27 ZG7 A4202 " -0.029 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA LEU A2728 " -0.018 2.00e-02 2.50e+03 3.57e-02 1.28e+01 pdb=" C LEU A2728 " 0.062 2.00e-02 2.50e+03 pdb=" O LEU A2728 " -0.023 2.00e-02 2.50e+03 pdb=" N GLU A2729 " -0.021 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C02 ZG7 A4202 " -0.003 2.00e-02 2.50e+03 2.42e-02 1.17e+01 pdb=" C03 ZG7 A4202 " -0.003 2.00e-02 2.50e+03 pdb=" C04 ZG7 A4202 " -0.006 2.00e-02 2.50e+03 pdb=" C05 ZG7 A4202 " 0.059 2.00e-02 2.50e+03 pdb=" C06 ZG7 A4202 " -0.002 2.00e-02 2.50e+03 pdb=" C07 ZG7 A4202 " -0.005 2.00e-02 2.50e+03 pdb=" C08 ZG7 A4202 " -0.034 2.00e-02 2.50e+03 pdb="CL01 ZG7 A4202 " -0.004 2.00e-02 2.50e+03 ... (remaining 1683 not shown) Histogram of nonbonded interaction distances: 1.98 - 2.56: 68 2.56 - 3.15: 8279 3.15 - 3.73: 14789 3.73 - 4.32: 20691 4.32 - 4.90: 33976 Nonbonded interactions: 77803 Sorted by model distance: nonbonded pdb=" O3G ATP A4203 " pdb="MG MG A4204 " model vdw 1.978 2.170 nonbonded pdb=" OG1 THR A2081 " pdb="MG MG A4204 " model vdw 2.003 2.170 nonbonded pdb=" O1B ATP A4203 " pdb="MG MG A4204 " model vdw 2.038 2.170 nonbonded pdb=" OD1 ASN A2832 " pdb=" OG SER A2905 " model vdw 2.263 2.440 nonbonded pdb=" NH1 ARG A2763 " pdb=" O SER A3511 " model vdw 2.298 2.520 ... (remaining 77798 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Set scattering table Set to: electron Number of scattering types: 8 Type Number sf(0) Gaussians Br 2 7.06 5 P 3 5.49 5 Mg 1 5.21 5 S 49 5.16 5 Cl 2 4.86 5 C 6197 2.51 5 N 1624 2.21 5 O 1804 1.98 5 sf(0) = scattering factor at diffraction angle 0. Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=0.50 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Write .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.760 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Construct map_model_manager: 0.000 Extract box with map and model: 10.800 Check model and map are aligned: 0.170 Convert atoms to be neutral: 0.100 Process input model: 29.510 Find NCS groups from input model: 0.270 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Set scattering table: 0.010 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.270 Set ADP refinement strategy: 0.000 Write .geo file: 0.000 Internal consistency checks: 0.000 Total: 44.930 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8637 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.012 0.241 9889 Z= 0.807 Angle : 1.516 26.107 13403 Z= 0.833 Chirality : 0.078 0.423 1520 Planarity : 0.008 0.054 1686 Dihedral : 12.695 89.407 3641 Min Nonbonded Distance : 1.978 Molprobity Statistics. All-atom Clashscore : 5.17 Ramachandran Plot: Outliers : 0.34 % Allowed : 6.82 % Favored : 92.85 % Rotamer Outliers : 1.94 % 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: -1.80 (0.22), residues: 1188 helix: -0.54 (0.19), residues: 584 sheet: -1.71 (0.43), residues: 125 loop : -1.75 (0.26), residues: 479 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 242 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 21 poor density : 221 time to evaluate : 1.324 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 21 outliers final: 7 residues processed: 238 average time/residue: 0.2687 time to fit residues: 84.7760 Evaluate side-chains 109 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 102 time to evaluate : 1.215 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 7 outliers final: 1 residues processed: 7 average time/residue: 0.1381 time to fit residues: 3.2521 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 100 optimal weight: 0.9990 chunk 90 optimal weight: 6.9990 chunk 50 optimal weight: 8.9990 chunk 30 optimal weight: 6.9990 chunk 60 optimal weight: 20.0000 chunk 48 optimal weight: 1.9990 chunk 93 optimal weight: 0.8980 chunk 36 optimal weight: 0.9990 chunk 56 optimal weight: 0.8980 chunk 69 optimal weight: 7.9990 chunk 108 optimal weight: 1.9990 overall best weight: 1.1586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A2274 HIS ** A2300 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A2363 ASN A2481 ASN A2536 ASN A2569 GLN A2854 ASN ** A3399 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A3471 ASN Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8624 moved from start: 0.2520 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.033 9889 Z= 0.190 Angle : 0.676 8.505 13403 Z= 0.346 Chirality : 0.044 0.191 1520 Planarity : 0.005 0.052 1686 Dihedral : 6.159 60.651 1300 Min Nonbonded Distance : 2.007 Molprobity Statistics. All-atom Clashscore : 9.93 Ramachandran Plot: Outliers : 0.17 % Allowed : 2.78 % Favored : 97.05 % Rotamer Outliers : 2.21 % 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.25 (0.24), residues: 1188 helix: 0.79 (0.21), residues: 601 sheet: -1.07 (0.43), residues: 138 loop : -1.09 (0.29), residues: 449 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 164 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 140 time to evaluate : 1.217 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 24 outliers final: 15 residues processed: 157 average time/residue: 0.2311 time to fit residues: 50.9670 Evaluate side-chains 120 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 105 time to evaluate : 1.263 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 15 outliers final: 1 residues processed: 15 average time/residue: 0.1312 time to fit residues: 4.6552 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 60 optimal weight: 5.9990 chunk 33 optimal weight: 0.4980 chunk 89 optimal weight: 1.9990 chunk 73 optimal weight: 1.9990 chunk 29 optimal weight: 10.0000 chunk 108 optimal weight: 3.9990 chunk 116 optimal weight: 8.9990 chunk 96 optimal weight: 3.9990 chunk 107 optimal weight: 4.9990 chunk 36 optimal weight: 3.9990 chunk 86 optimal weight: 7.9990 overall best weight: 2.4988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A2300 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A2500 GLN A2598 HIS A3471 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8679 moved from start: 0.3013 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.076 9889 Z= 0.252 Angle : 0.621 8.292 13403 Z= 0.318 Chirality : 0.043 0.145 1520 Planarity : 0.004 0.043 1686 Dihedral : 5.894 59.228 1300 Min Nonbonded Distance : 1.928 Molprobity Statistics. All-atom Clashscore : 10.71 Ramachandran Plot: Outliers : 0.08 % Allowed : 3.45 % Favored : 96.46 % Rotamer Outliers : 2.40 % 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.07 (0.25), residues: 1188 helix: 1.04 (0.22), residues: 601 sheet: -0.90 (0.43), residues: 143 loop : -0.92 (0.29), residues: 444 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 136 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 26 poor density : 110 time to evaluate : 1.366 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 26 outliers final: 14 residues processed: 130 average time/residue: 0.2334 time to fit residues: 43.0729 Evaluate side-chains 111 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 97 time to evaluate : 1.297 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 14 outliers final: 1 residues processed: 14 average time/residue: 0.1191 time to fit residues: 4.5126 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 106 optimal weight: 0.7980 chunk 81 optimal weight: 5.9990 chunk 56 optimal weight: 0.9980 chunk 11 optimal weight: 0.7980 chunk 51 optimal weight: 9.9990 chunk 72 optimal weight: 4.9990 chunk 108 optimal weight: 1.9990 chunk 115 optimal weight: 0.0040 chunk 102 optimal weight: 0.5980 chunk 30 optimal weight: 10.0000 chunk 95 optimal weight: 0.5980 overall best weight: 0.5592 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A3399 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A3471 ASN ** A3497 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8630 moved from start: 0.3484 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.051 9889 Z= 0.151 Angle : 0.590 9.206 13403 Z= 0.293 Chirality : 0.041 0.161 1520 Planarity : 0.004 0.073 1686 Dihedral : 5.656 59.201 1300 Min Nonbonded Distance : 2.003 Molprobity Statistics. All-atom Clashscore : 8.74 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.11 % Favored : 96.89 % Rotamer Outliers : 1.01 % 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.33 (0.25), residues: 1188 helix: 1.16 (0.22), residues: 597 sheet: -0.56 (0.43), residues: 143 loop : -0.72 (0.30), residues: 448 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 138 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 127 time to evaluate : 1.294 Fit side-chains revert: symmetry clash outliers start: 11 outliers final: 4 residues processed: 134 average time/residue: 0.2227 time to fit residues: 42.7990 Evaluate side-chains 96 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 92 time to evaluate : 1.325 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 4 outliers final: 1 residues processed: 4 average time/residue: 0.1718 time to fit residues: 2.7715 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 65 optimal weight: 8.9990 chunk 1 optimal weight: 8.9990 chunk 85 optimal weight: 9.9990 chunk 47 optimal weight: 4.9990 chunk 98 optimal weight: 8.9990 chunk 79 optimal weight: 8.9990 chunk 0 optimal weight: 5.9990 chunk 58 optimal weight: 10.0000 chunk 103 optimal weight: 4.9990 chunk 29 optimal weight: 7.9990 chunk 38 optimal weight: 4.9990 overall best weight: 5.7990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A2634 ASN A3471 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8770 moved from start: 0.3855 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.007 0.106 9889 Z= 0.481 Angle : 0.719 9.952 13403 Z= 0.370 Chirality : 0.046 0.176 1520 Planarity : 0.005 0.065 1686 Dihedral : 5.978 58.348 1300 Min Nonbonded Distance : 1.857 Molprobity Statistics. All-atom Clashscore : 12.26 Ramachandran Plot: Outliers : 0.08 % Allowed : 4.29 % Favored : 95.62 % Rotamer Outliers : 2.49 % 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.16 (0.24), residues: 1188 helix: 1.01 (0.21), residues: 592 sheet: -0.46 (0.44), residues: 140 loop : -0.84 (0.29), residues: 456 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 121 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 94 time to evaluate : 1.255 Fit side-chains outliers start: 27 outliers final: 13 residues processed: 114 average time/residue: 0.2132 time to fit residues: 35.7389 Evaluate side-chains 103 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 90 time to evaluate : 1.217 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 13 outliers final: 1 residues processed: 13 average time/residue: 0.1153 time to fit residues: 4.1202 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 103 optimal weight: 3.9990 chunk 22 optimal weight: 6.9990 chunk 67 optimal weight: 5.9990 chunk 28 optimal weight: 3.9990 chunk 115 optimal weight: 2.9990 chunk 95 optimal weight: 8.9990 chunk 53 optimal weight: 4.9990 chunk 9 optimal weight: 4.9990 chunk 38 optimal weight: 0.7980 chunk 60 optimal weight: 0.8980 chunk 111 optimal weight: 0.9990 overall best weight: 1.9386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A3471 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8713 moved from start: 0.4060 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.043 9889 Z= 0.203 Angle : 0.602 9.107 13403 Z= 0.299 Chirality : 0.042 0.169 1520 Planarity : 0.004 0.055 1686 Dihedral : 5.692 58.781 1300 Min Nonbonded Distance : 1.974 Molprobity Statistics. All-atom Clashscore : 9.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.79 % Favored : 96.21 % Rotamer Outliers : 1.66 % 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.41 (0.25), residues: 1188 helix: 1.20 (0.22), residues: 589 sheet: -0.33 (0.43), residues: 144 loop : -0.69 (0.29), residues: 455 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 119 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 101 time to evaluate : 1.276 Fit side-chains revert: symmetry clash outliers start: 18 outliers final: 4 residues processed: 114 average time/residue: 0.2333 time to fit residues: 38.2419 Evaluate side-chains 94 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 90 time to evaluate : 1.364 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 4 outliers final: 1 residues processed: 4 average time/residue: 0.2048 time to fit residues: 2.9583 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 12 optimal weight: 0.0970 chunk 65 optimal weight: 9.9990 chunk 84 optimal weight: 5.9990 chunk 96 optimal weight: 8.9990 chunk 64 optimal weight: 9.9990 chunk 114 optimal weight: 10.0000 chunk 71 optimal weight: 0.7980 chunk 69 optimal weight: 5.9990 chunk 52 optimal weight: 4.9990 chunk 70 optimal weight: 0.9990 chunk 45 optimal weight: 7.9990 overall best weight: 2.5784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A3453 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A3497 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8726 moved from start: 0.4215 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.044 9889 Z= 0.245 Angle : 0.604 8.805 13403 Z= 0.301 Chirality : 0.042 0.170 1520 Planarity : 0.004 0.058 1686 Dihedral : 5.648 58.352 1300 Min Nonbonded Distance : 1.933 Molprobity Statistics. All-atom Clashscore : 9.93 Ramachandran Plot: Outliers : 0.08 % Allowed : 4.12 % Favored : 95.79 % Rotamer Outliers : 1.20 % 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.56 (0.25), residues: 1188 helix: 1.28 (0.22), residues: 584 sheet: -0.12 (0.41), residues: 155 loop : -0.56 (0.29), residues: 449 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 110 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 97 time to evaluate : 1.318 Fit side-chains revert: symmetry clash outliers start: 13 outliers final: 6 residues processed: 108 average time/residue: 0.2405 time to fit residues: 37.6639 Evaluate side-chains 95 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 89 time to evaluate : 1.321 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 6 outliers final: 1 residues processed: 6 average time/residue: 0.1514 time to fit residues: 3.1478 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 68 optimal weight: 0.9990 chunk 34 optimal weight: 0.2980 chunk 22 optimal weight: 4.9990 chunk 72 optimal weight: 10.0000 chunk 78 optimal weight: 10.0000 chunk 56 optimal weight: 0.9990 chunk 10 optimal weight: 2.9990 chunk 90 optimal weight: 0.6980 chunk 104 optimal weight: 4.9990 chunk 109 optimal weight: 8.9990 chunk 100 optimal weight: 0.9990 overall best weight: 0.7986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A3453 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A3471 ASN A3497 HIS ** A3521 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8682 moved from start: 0.4404 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.030 9889 Z= 0.143 Angle : 0.579 10.452 13403 Z= 0.285 Chirality : 0.041 0.156 1520 Planarity : 0.004 0.054 1686 Dihedral : 5.461 59.022 1300 Min Nonbonded Distance : 2.019 Molprobity Statistics. All-atom Clashscore : 8.69 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.62 % Favored : 96.38 % Rotamer Outliers : 0.83 % 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.67 (0.25), residues: 1188 helix: 1.41 (0.22), residues: 578 sheet: -0.26 (0.43), residues: 147 loop : -0.47 (0.29), residues: 463 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 116 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 107 time to evaluate : 1.287 Fit side-chains revert: symmetry clash outliers start: 9 outliers final: 1 residues processed: 113 average time/residue: 0.2311 time to fit residues: 37.5876 Evaluate side-chains 96 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 95 time to evaluate : 1.252 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 1 outliers final: 1 residues processed: 1 average time/residue: 0.4589 time to fit residues: 2.4511 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 106 optimal weight: 5.9990 chunk 109 optimal weight: 6.9990 chunk 64 optimal weight: 0.9990 chunk 46 optimal weight: 5.9990 chunk 83 optimal weight: 10.0000 chunk 32 optimal weight: 0.7980 chunk 96 optimal weight: 5.9990 chunk 101 optimal weight: 8.9990 chunk 70 optimal weight: 0.3980 chunk 113 optimal weight: 0.7980 chunk 68 optimal weight: 2.9990 overall best weight: 1.1984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A3453 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A3471 ASN ** A3521 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8682 moved from start: 0.4496 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.060 9889 Z= 0.168 Angle : 0.581 10.613 13403 Z= 0.288 Chirality : 0.041 0.172 1520 Planarity : 0.004 0.057 1686 Dihedral : 5.396 58.799 1300 Min Nonbonded Distance : 1.987 Molprobity Statistics. All-atom Clashscore : 9.57 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.87 % Favored : 96.13 % Rotamer Outliers : 1.11 % 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.72 (0.25), residues: 1188 helix: 1.48 (0.22), residues: 575 sheet: -0.16 (0.43), residues: 147 loop : -0.50 (0.29), residues: 466 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 111 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 99 time to evaluate : 1.234 Fit side-chains revert: symmetry clash outliers start: 12 outliers final: 5 residues processed: 106 average time/residue: 0.2157 time to fit residues: 34.0588 Evaluate side-chains 100 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 95 time to evaluate : 1.303 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 5 outliers final: 1 residues processed: 5 average time/residue: 0.2143 time to fit residues: 3.1601 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 53 optimal weight: 0.0000 chunk 78 optimal weight: 9.9990 chunk 118 optimal weight: 0.9990 chunk 109 optimal weight: 2.9990 chunk 94 optimal weight: 7.9990 chunk 9 optimal weight: 0.8980 chunk 72 optimal weight: 7.9990 chunk 57 optimal weight: 1.9990 chunk 74 optimal weight: 4.9990 chunk 100 optimal weight: 0.9990 chunk 28 optimal weight: 7.9990 overall best weight: 0.9790 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A2741 HIS ** A3453 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A3471 ASN ** A3521 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8669 moved from start: 0.4648 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.029 9889 Z= 0.148 Angle : 0.591 11.448 13403 Z= 0.289 Chirality : 0.041 0.255 1520 Planarity : 0.004 0.055 1686 Dihedral : 5.329 58.685 1300 Min Nonbonded Distance : 1.999 Molprobity Statistics. All-atom Clashscore : 9.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.70 % Favored : 96.30 % Rotamer Outliers : 0.55 % 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.76 (0.25), residues: 1188 helix: 1.44 (0.22), residues: 583 sheet: 0.01 (0.44), residues: 147 loop : -0.47 (0.29), residues: 458 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 112 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 106 time to evaluate : 1.556 Fit side-chains revert: symmetry clash outliers start: 6 outliers final: 5 residues processed: 109 average time/residue: 0.2108 time to fit residues: 34.5192 Evaluate side-chains 98 residues out of total 1084 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 93 time to evaluate : 1.333 Switching outliers to nearest non-outliers revert: symmetry clash outliers start: 5 outliers final: 1 residues processed: 5 average time/residue: 0.1696 time to fit residues: 2.9260 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 87 optimal weight: 0.9990 chunk 13 optimal weight: 0.9990 chunk 26 optimal weight: 7.9990 chunk 94 optimal weight: 7.9990 chunk 39 optimal weight: 4.9990 chunk 97 optimal weight: 9.9990 chunk 11 optimal weight: 6.9990 chunk 17 optimal weight: 0.9990 chunk 83 optimal weight: 10.0000 chunk 5 optimal weight: 6.9990 chunk 68 optimal weight: 0.7980 overall best weight: 1.7588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A3453 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A3521 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3272 r_free = 0.3272 target = 0.071933 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 38)----------------| | r_work = 0.2919 r_free = 0.2919 target = 0.056647 restraints weight = 31104.209| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 34)----------------| | r_work = 0.2970 r_free = 0.2970 target = 0.058817 restraints weight = 15592.346| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 34)----------------| | r_work = 0.3004 r_free = 0.3004 target = 0.060328 restraints weight = 9912.130| |-----------------------------------------------------------------------------| r_work (final): 0.3001 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3272 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3272 r_free = 0.3272 target_work(ls_wunit_k1) = 0.072 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3272 r_free = 0.3272 target_work(ls_wunit_k1) = 0.072 | | occupancies: max = 1.00 min = 0.00 number of occupancies < 0.1: 2 | |-----------------------------------------------------------------------------| r_final: 0.3272 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8686 moved from start: 0.4673 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.093 9889 Z= 0.216 Angle : 0.598 11.793 13403 Z= 0.299 Chirality : 0.042 0.247 1520 Planarity : 0.005 0.112 1686 Dihedral : 5.335 57.881 1300 Min Nonbonded Distance : 1.964 Molprobity Statistics. All-atom Clashscore : 10.40 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.96 % Favored : 96.04 % Rotamer Outliers : 0.37 % 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.81 (0.25), residues: 1188 helix: 1.54 (0.22), residues: 576 sheet: 0.05 (0.44), residues: 147 loop : -0.51 (0.29), residues: 465 =============================================================================== Job complete usr+sys time: 2355.24 seconds wall clock time: 46 minutes 19.24 seconds (2779.24 seconds total)