Starting phenix.real_space_refine (version: dev) on Sun Feb 19 06:41:22 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/02_2023/6lr0_31583.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/02_2023/6lr0_31583.map 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 { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/02_2023/6lr0_31583.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/02_2023/6lr0_31583.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/02_2023/6lr0_31583.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/02_2023/6lr0_31583.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.007 sd= 0.098 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Process input model Symmetric amino acids flipped Residue "U ARG 44": "NH1" <-> "NH2" Residue "U PHE 51": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "U ARG 184": "NH1" <-> "NH2" Residue "U ARG 237": "NH1" <-> "NH2" Residue "U ARG 289": "NH1" <-> "NH2" Residue "U ARG 299": "NH1" <-> "NH2" Residue "U ARG 317": "NH1" <-> "NH2" Residue "U ARG 415": "NH1" <-> "NH2" Residue "U ARG 432": "NH1" <-> "NH2" Residue "U ARG 487": "NH1" <-> "NH2" Residue "U ARG 571": "NH1" <-> "NH2" Residue "U ARG 575": "NH1" <-> "NH2" Residue "U ARG 637": "NH1" <-> "NH2" Residue "U ARG 646": "NH1" <-> "NH2" Residue "U ARG 692": "NH1" <-> "NH2" Residue "U ARG 832": "NH1" <-> "NH2" Residue "U ARG 860": "NH1" <-> "NH2" Residue "U ARG 920": "NH1" <-> "NH2" Residue "U ARG 958": "NH1" <-> "NH2" Residue "U ARG 1019": "NH1" <-> "NH2" Residue "U ARG 1187": "NH1" <-> "NH2" Residue "U ARG 1221": "NH1" <-> "NH2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/chem_data/mon_lib" Total number of atoms: 9576 Number of models: 1 Model: "" Number of chains: 1 Chain: "U" Number of atoms: 9576 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1230, 9576 Classifications: {'peptide': 1230} Link IDs: {'PTRANS': 36, 'TRANS': 1193} Chain breaks: 2 Time building chain proxies: 5.25, per 1000 atoms: 0.55 Number of scatterers: 9576 At special positions: 0 Unit cell: (90.1, 160.06, 112.36, 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 54 16.00 O 1766 8.00 N 1627 7.00 C 6129 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: 3.91 Conformation dependent library (CDL) restraints added in 1.6 seconds 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2272 Finding SS restraints... Secondary structure from input PDB file: 61 helices and 4 sheets defined 36.3% alpha, 2.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.59 Creating SS restraints... Processing helix chain 'U' and resid 46 through 50 Processing helix chain 'U' and resid 79 through 84 Processing helix chain 'U' and resid 88 through 100 removed outlier: 3.716A pdb=" N GLU U 96 " --> pdb=" O ASP U 92 " (cutoff:3.500A) removed outlier: 3.636A pdb=" N GLU U 99 " --> pdb=" O VAL U 95 " (cutoff:3.500A) Processing helix chain 'U' and resid 147 through 155 removed outlier: 4.004A pdb=" N LEU U 153 " --> pdb=" O ALA U 149 " (cutoff:3.500A) Processing helix chain 'U' and resid 173 through 183 removed outlier: 3.655A pdb=" N PHE U 179 " --> pdb=" O ARG U 175 " (cutoff:3.500A) Processing helix chain 'U' and resid 187 through 192 removed outlier: 3.734A pdb=" N PHE U 190 " --> pdb=" O ILE U 187 " (cutoff:3.500A) removed outlier: 3.560A pdb=" N CYS U 192 " --> pdb=" O TRP U 189 " (cutoff:3.500A) Processing helix chain 'U' and resid 230 through 238 removed outlier: 3.769A pdb=" N PHE U 236 " --> pdb=" O LEU U 232 " (cutoff:3.500A) removed outlier: 3.943A pdb=" N ARG U 237 " --> pdb=" O LEU U 233 " (cutoff:3.500A) removed outlier: 3.653A pdb=" N GLY U 238 " --> pdb=" O GLY U 234 " (cutoff:3.500A) Processing helix chain 'U' and resid 240 through 245 removed outlier: 3.740A pdb=" N VAL U 244 " --> pdb=" O LYS U 240 " (cutoff:3.500A) Processing helix chain 'U' and resid 248 through 257 removed outlier: 4.183A pdb=" N ILE U 254 " --> pdb=" O PRO U 250 " (cutoff:3.500A) removed outlier: 3.830A pdb=" N GLY U 255 " --> pdb=" O LEU U 251 " (cutoff:3.500A) Processing helix chain 'U' and resid 258 through 263 Processing helix chain 'U' and resid 266 through 275 removed outlier: 3.567A pdb=" N LEU U 271 " --> pdb=" O THR U 267 " (cutoff:3.500A) removed outlier: 4.413A pdb=" N LYS U 272 " --> pdb=" O ASP U 268 " (cutoff:3.500A) removed outlier: 4.021A pdb=" N ALA U 275 " --> pdb=" O LEU U 271 " (cutoff:3.500A) Processing helix chain 'U' and resid 277 through 282 removed outlier: 3.706A pdb=" N ASP U 282 " --> pdb=" O GLY U 278 " (cutoff:3.500A) Processing helix chain 'U' and resid 287 through 294 removed outlier: 3.678A pdb=" N VAL U 291 " --> pdb=" O SER U 287 " (cutoff:3.500A) removed outlier: 3.718A pdb=" N ALA U 293 " --> pdb=" O ARG U 289 " (cutoff:3.500A) removed outlier: 3.741A pdb=" N PHE U 294 " --> pdb=" O THR U 290 " (cutoff:3.500A) Processing helix chain 'U' and resid 308 through 312 removed outlier: 3.684A pdb=" N ALA U 311 " --> pdb=" O LEU U 308 " (cutoff:3.500A) Processing helix chain 'U' and resid 320 through 325 Processing helix chain 'U' and resid 329 through 334 removed outlier: 3.919A pdb=" N PHE U 334 " --> pdb=" O TRP U 330 " (cutoff:3.500A) Processing helix chain 'U' and resid 335 through 344 removed outlier: 3.978A pdb=" N LEU U 339 " --> pdb=" O LEU U 335 " (cutoff:3.500A) Processing helix chain 'U' and resid 363 through 376 removed outlier: 3.621A pdb=" N ILE U 367 " --> pdb=" O PHE U 363 " (cutoff:3.500A) removed outlier: 3.992A pdb=" N GLY U 369 " --> pdb=" O SER U 365 " (cutoff:3.500A) removed outlier: 4.639A pdb=" N LEU U 373 " --> pdb=" O GLY U 369 " (cutoff:3.500A) removed outlier: 3.849A pdb=" N GLY U 374 " --> pdb=" O ALA U 370 " (cutoff:3.500A) Processing helix chain 'U' and resid 377 through 386 removed outlier: 4.929A pdb=" N PHE U 383 " --> pdb=" O CYS U 379 " (cutoff:3.500A) removed outlier: 4.219A pdb=" N ALA U 384 " --> pdb=" O LEU U 380 " (cutoff:3.500A) Processing helix chain 'U' and resid 392 through 397 removed outlier: 3.716A pdb=" N THR U 396 " --> pdb=" O SER U 392 " (cutoff:3.500A) Processing helix chain 'U' and resid 462 through 467 Processing helix chain 'U' and resid 490 through 495 Processing helix chain 'U' and resid 525 through 535 removed outlier: 3.749A pdb=" N GLN U 530 " --> pdb=" O GLU U 526 " (cutoff:3.500A) Processing helix chain 'U' and resid 540 through 544 removed outlier: 3.508A pdb=" N LEU U 544 " --> pdb=" O ILE U 541 " (cutoff:3.500A) Processing helix chain 'U' and resid 563 through 575 removed outlier: 3.615A pdb=" N ALA U 570 " --> pdb=" O ARG U 566 " (cutoff:3.500A) removed outlier: 4.398A pdb=" N ARG U 575 " --> pdb=" O ARG U 571 " (cutoff:3.500A) Processing helix chain 'U' and resid 593 through 603 removed outlier: 3.766A pdb=" N VAL U 597 " --> pdb=" O SER U 593 " (cutoff:3.500A) removed outlier: 3.886A pdb=" N VAL U 600 " --> pdb=" O MET U 596 " (cutoff:3.500A) Processing helix chain 'U' and resid 652 through 656 Processing helix chain 'U' and resid 700 through 704 removed outlier: 3.660A pdb=" N SER U 704 " --> pdb=" O SER U 701 " (cutoff:3.500A) Processing helix chain 'U' and resid 742 through 747 Processing helix chain 'U' and resid 748 through 752 Processing helix chain 'U' and resid 775 through 780 removed outlier: 3.556A pdb=" N LEU U 780 " --> pdb=" O PHE U 776 " (cutoff:3.500A) Processing helix chain 'U' and resid 789 through 802 removed outlier: 3.621A pdb=" N ILE U 795 " --> pdb=" O GLN U 791 " (cutoff:3.500A) removed outlier: 3.542A pdb=" N ASN U 796 " --> pdb=" O ARG U 792 " (cutoff:3.500A) removed outlier: 3.583A pdb=" N CYS U 799 " --> pdb=" O ILE U 795 " (cutoff:3.500A) Processing helix chain 'U' and resid 803 through 806 Processing helix chain 'U' and resid 828 through 833 removed outlier: 3.685A pdb=" N LYS U 833 " --> pdb=" O LYS U 829 " (cutoff:3.500A) Processing helix chain 'U' and resid 834 through 836 No H-bonds generated for 'chain 'U' and resid 834 through 836' Processing helix chain 'U' and resid 838 through 842 removed outlier: 3.983A pdb=" N GLN U 842 " --> pdb=" O MET U 839 " (cutoff:3.500A) Processing helix chain 'U' and resid 843 through 847 removed outlier: 4.072A pdb=" N TRP U 846 " --> pdb=" O ASP U 843 " (cutoff:3.500A) Processing helix chain 'U' and resid 875 through 881 Processing helix chain 'U' and resid 885 through 897 removed outlier: 3.825A pdb=" N ALA U 889 " --> pdb=" O ASN U 885 " (cutoff:3.500A) removed outlier: 3.884A pdb=" N ILE U 891 " --> pdb=" O THR U 887 " (cutoff:3.500A) removed outlier: 3.644A pdb=" N PHE U 894 " --> pdb=" O MET U 890 " (cutoff:3.500A) removed outlier: 3.967A pdb=" N PHE U 896 " --> pdb=" O ILE U 892 " (cutoff:3.500A) removed outlier: 3.529A pdb=" N SER U 897 " --> pdb=" O ALA U 893 " (cutoff:3.500A) Processing helix chain 'U' and resid 904 through 907 Processing helix chain 'U' and resid 908 through 913 Processing helix chain 'U' and resid 925 through 930 removed outlier: 3.774A pdb=" N LYS U 930 " --> pdb=" O ALA U 926 " (cutoff:3.500A) Processing helix chain 'U' and resid 935 through 940 Processing helix chain 'U' and resid 940 through 945 removed outlier: 3.555A pdb=" N SER U 945 " --> pdb=" O ASN U 941 " (cutoff:3.500A) Processing helix chain 'U' and resid 946 through 951 Processing helix chain 'U' and resid 956 through 962 Processing helix chain 'U' and resid 994 through 1007 removed outlier: 3.639A pdb=" N SER U 999 " --> pdb=" O ALA U 995 " (cutoff:3.500A) removed outlier: 3.766A pdb=" N TYR U1000 " --> pdb=" O ASN U 996 " (cutoff:3.500A) removed outlier: 3.759A pdb=" N ARG U1001 " --> pdb=" O SER U 997 " (cutoff:3.500A) Processing helix chain 'U' and resid 1013 through 1015 No H-bonds generated for 'chain 'U' and resid 1013 through 1015' Processing helix chain 'U' and resid 1016 through 1021 removed outlier: 3.614A pdb=" N ILE U1021 " --> pdb=" O VAL U1017 " (cutoff:3.500A) Processing helix chain 'U' and resid 1031 through 1035 removed outlier: 3.917A pdb=" N PHE U1035 " --> pdb=" O GLY U1032 " (cutoff:3.500A) Processing helix chain 'U' and resid 1037 through 1043 removed outlier: 3.558A pdb=" N TYR U1041 " --> pdb=" O TYR U1037 " (cutoff:3.500A) removed outlier: 3.808A pdb=" N ALA U1042 " --> pdb=" O THR U1038 " (cutoff:3.500A) removed outlier: 3.690A pdb=" N LYS U1043 " --> pdb=" O PRO U1039 " (cutoff:3.500A) No H-bonds generated for 'chain 'U' and resid 1037 through 1043' Processing helix chain 'U' and resid 1048 through 1053 removed outlier: 3.843A pdb=" N PHE U1052 " --> pdb=" O ALA U1048 " (cutoff:3.500A) Processing helix chain 'U' and resid 1120 through 1125 removed outlier: 3.553A pdb=" N GLN U1124 " --> pdb=" O SER U1120 " (cutoff:3.500A) Processing helix chain 'U' and resid 1148 through 1153 removed outlier: 3.509A pdb=" N LEU U1152 " --> pdb=" O ASN U1148 " (cutoff:3.500A) Processing helix chain 'U' and resid 1171 through 1176 removed outlier: 3.714A pdb=" N LYS U1175 " --> pdb=" O MET U1171 " (cutoff:3.500A) Processing helix chain 'U' and resid 1185 through 1191 removed outlier: 4.306A pdb=" N ILE U1189 " --> pdb=" O MET U1185 " (cutoff:3.500A) removed outlier: 3.746A pdb=" N ALA U1190 " --> pdb=" O GLU U1186 " (cutoff:3.500A) Processing helix chain 'U' and resid 1200 through 1204 removed outlier: 3.822A pdb=" N LEU U1204 " --> pdb=" O VAL U1201 " (cutoff:3.500A) Processing helix chain 'U' and resid 1220 through 1223 Processing helix chain 'U' and resid 1224 through 1233 removed outlier: 3.830A pdb=" N ALA U1228 " --> pdb=" O LYS U1224 " (cutoff:3.500A) Processing helix chain 'U' and resid 1250 through 1264 removed outlier: 3.537A pdb=" N VAL U1257 " --> pdb=" O SER U1253 " (cutoff:3.500A) removed outlier: 3.651A pdb=" N VAL U1259 " --> pdb=" O LYS U1255 " (cutoff:3.500A) removed outlier: 3.685A pdb=" N ALA U1260 " --> pdb=" O THR U1256 " (cutoff:3.500A) Processing helix chain 'U' and resid 1308 through 1314 removed outlier: 3.519A pdb=" N LYS U1312 " --> pdb=" O GLY U1308 " (cutoff:3.500A) removed outlier: 3.577A pdb=" N LEU U1313 " --> pdb=" O ALA U1309 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'U' and resid 438 through 445 removed outlier: 6.806A pdb=" N ASN U 439 " --> pdb=" O THR U 426 " (cutoff:3.500A) removed outlier: 6.771A pdb=" N THR U 426 " --> pdb=" O ASN U 439 " (cutoff:3.500A) removed outlier: 5.710A pdb=" N LEU U 441 " --> pdb=" O ASN U 424 " (cutoff:3.500A) removed outlier: 3.724A pdb=" N ASN U 424 " --> pdb=" O LEU U 441 " (cutoff:3.500A) removed outlier: 3.554A pdb=" N ILE U 420 " --> pdb=" O ILE U 445 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'U' and resid 450 through 451 removed outlier: 6.175A pdb=" N THR U 451 " --> pdb=" O SER U 612 " (cutoff:3.500A) removed outlier: 6.082A pdb=" N LEU U 580 " --> pdb=" O ILE U 611 " (cutoff:3.500A) removed outlier: 7.556A pdb=" N VAL U 613 " --> pdb=" O LEU U 580 " (cutoff:3.500A) removed outlier: 6.191A pdb=" N LEU U 582 " --> pdb=" O VAL U 613 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'U' and resid 1079 through 1080 Processing sheet with id=AA4, first strand: chain 'U' and resid 1156 through 1159 removed outlier: 6.714A pdb=" N GLY U1157 " --> pdb=" O LEU U1241 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA4 158 hydrogen bonds defined for protein. 420 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 2.91 Time building geometry restraints manager: 3.79 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.20 - 1.33: 2048 1.33 - 1.45: 2223 1.45 - 1.57: 5402 1.57 - 1.69: 0 1.69 - 1.81: 88 Bond restraints: 9761 Sorted by residual: bond pdb=" C GLN U 76 " pdb=" N PRO U 77 " ideal model delta sigma weight residual 1.333 1.378 -0.045 7.80e-03 1.64e+04 3.37e+01 bond pdb=" CG1 ILE U 171 " pdb=" CD1 ILE U 171 " ideal model delta sigma weight residual 1.513 1.363 0.150 3.90e-02 6.57e+02 1.49e+01 bond pdb=" CB VAL U1024 " pdb=" CG2 VAL U1024 " ideal model delta sigma weight residual 1.521 1.421 0.100 3.30e-02 9.18e+02 9.18e+00 bond pdb=" CG ARG U 175 " pdb=" CD ARG U 175 " ideal model delta sigma weight residual 1.520 1.445 0.075 3.00e-02 1.11e+03 6.22e+00 bond pdb=" CB ILE U 160 " pdb=" CG2 ILE U 160 " ideal model delta sigma weight residual 1.521 1.442 0.079 3.30e-02 9.18e+02 5.66e+00 ... (remaining 9756 not shown) Histogram of bond angle deviations from ideal: 97.33 - 104.71: 145 104.71 - 112.08: 4357 112.08 - 119.46: 3782 119.46 - 126.83: 4806 126.83 - 134.20: 109 Bond angle restraints: 13199 Sorted by residual: angle pdb=" N ILE U1170 " pdb=" CA ILE U1170 " pdb=" C ILE U1170 " ideal model delta sigma weight residual 112.83 104.69 8.14 9.90e-01 1.02e+00 6.75e+01 angle pdb=" N ILE U 213 " pdb=" CA ILE U 213 " pdb=" C ILE U 213 " ideal model delta sigma weight residual 113.71 107.32 6.39 9.50e-01 1.11e+00 4.53e+01 angle pdb=" N LEU U 67 " pdb=" CA LEU U 67 " pdb=" C LEU U 67 " ideal model delta sigma weight residual 113.38 105.94 7.44 1.23e+00 6.61e-01 3.66e+01 angle pdb=" C ASN U 946 " pdb=" N ILE U 947 " pdb=" CA ILE U 947 " ideal model delta sigma weight residual 120.53 127.26 -6.73 1.41e+00 5.03e-01 2.28e+01 angle pdb=" N GLY U 924 " pdb=" CA GLY U 924 " pdb=" C GLY U 924 " ideal model delta sigma weight residual 115.46 107.75 7.71 1.63e+00 3.76e-01 2.24e+01 ... (remaining 13194 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 16.28: 5279 16.28 - 32.56: 428 32.56 - 48.84: 99 48.84 - 65.12: 10 65.12 - 81.40: 7 Dihedral angle restraints: 5823 sinusoidal: 2288 harmonic: 3535 Sorted by residual: dihedral pdb=" CA LEU U 712 " pdb=" C LEU U 712 " pdb=" N ALA U 713 " pdb=" CA ALA U 713 " ideal model delta harmonic sigma weight residual 180.00 135.33 44.67 0 5.00e+00 4.00e-02 7.98e+01 dihedral pdb=" CA TYR U1037 " pdb=" C TYR U1037 " pdb=" N THR U1038 " pdb=" CA THR U1038 " ideal model delta harmonic sigma weight residual -180.00 -142.24 -37.76 0 5.00e+00 4.00e-02 5.70e+01 dihedral pdb=" CA ILE U1271 " pdb=" C ILE U1271 " pdb=" N VAL U1272 " pdb=" CA VAL U1272 " ideal model delta harmonic sigma weight residual -180.00 -146.00 -34.00 0 5.00e+00 4.00e-02 4.62e+01 ... (remaining 5820 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.093: 1305 0.093 - 0.187: 177 0.187 - 0.280: 16 0.280 - 0.374: 3 0.374 - 0.467: 1 Chirality restraints: 1502 Sorted by residual: chirality pdb=" CB ILE U 785 " pdb=" CA ILE U 785 " pdb=" CG1 ILE U 785 " pdb=" CG2 ILE U 785 " both_signs ideal model delta sigma weight residual False 2.64 2.18 0.47 2.00e-01 2.50e+01 5.45e+00 chirality pdb=" CB VAL U 808 " pdb=" CA VAL U 808 " pdb=" CG1 VAL U 808 " pdb=" CG2 VAL U 808 " both_signs ideal model delta sigma weight residual False -2.63 -2.29 -0.33 2.00e-01 2.50e+01 2.80e+00 chirality pdb=" CB VAL U 360 " pdb=" CA VAL U 360 " pdb=" CG1 VAL U 360 " pdb=" CG2 VAL U 360 " both_signs ideal model delta sigma weight residual False -2.63 -2.32 -0.31 2.00e-01 2.50e+01 2.35e+00 ... (remaining 1499 not shown) Planarity restraints: 1689 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C LYS U 969 " 0.060 5.00e-02 4.00e+02 8.99e-02 1.29e+01 pdb=" N PRO U 970 " -0.155 5.00e-02 4.00e+02 pdb=" CA PRO U 970 " 0.045 5.00e-02 4.00e+02 pdb=" CD PRO U 970 " 0.050 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CB PHE U 383 " -0.023 2.00e-02 2.50e+03 2.06e-02 7.44e+00 pdb=" CG PHE U 383 " 0.046 2.00e-02 2.50e+03 pdb=" CD1 PHE U 383 " -0.012 2.00e-02 2.50e+03 pdb=" CD2 PHE U 383 " 0.002 2.00e-02 2.50e+03 pdb=" CE1 PHE U 383 " 0.002 2.00e-02 2.50e+03 pdb=" CE2 PHE U 383 " -0.013 2.00e-02 2.50e+03 pdb=" CZ PHE U 383 " -0.002 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C SER U 249 " -0.044 5.00e-02 4.00e+02 6.62e-02 7.01e+00 pdb=" N PRO U 250 " 0.114 5.00e-02 4.00e+02 pdb=" CA PRO U 250 " -0.033 5.00e-02 4.00e+02 pdb=" CD PRO U 250 " -0.037 5.00e-02 4.00e+02 ... (remaining 1686 not shown) Histogram of nonbonded interaction distances: 2.29 - 2.81: 2970 2.81 - 3.33: 8225 3.33 - 3.86: 14900 3.86 - 4.38: 17263 4.38 - 4.90: 29290 Nonbonded interactions: 72648 Sorted by model distance: nonbonded pdb=" O GLN U 546 " pdb=" OG1 THR U 550 " model vdw 2.291 2.440 nonbonded pdb=" OG1 THR U 227 " pdb=" O GLY U 369 " model vdw 2.295 2.440 nonbonded pdb=" O SER U1027 " pdb=" OG SER U1027 " model vdw 2.308 2.440 nonbonded pdb=" O GLN U1196 " pdb=" NH2 ARG U1226 " model vdw 2.316 2.520 nonbonded pdb=" O THR U 155 " pdb=" OG1 THR U 155 " model vdw 2.330 2.440 ... (remaining 72643 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: 4 Type Number sf(0) Gaussians S 54 5.16 5 C 6129 2.51 5 N 1627 2.21 5 O 1766 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=1.00 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 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: 0.660 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: 2.910 Check model and map are aligned: 0.150 Process input model: 26.830 Find NCS groups from input model: 0.260 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Set scattering table: 0.100 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.560 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 33.520 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7773 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.011 0.150 9761 Z= 0.698 Angle : 1.380 14.085 13199 Z= 0.753 Chirality : 0.066 0.467 1502 Planarity : 0.007 0.090 1689 Dihedral : 13.072 81.405 3551 Min Nonbonded Distance : 2.291 Molprobity Statistics. All-atom Clashscore : 5.77 Ramachandran Plot: Outliers : 0.08 % Allowed : 11.03 % Favored : 88.89 % Rotamer Outliers : 0.88 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.50 % 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: -6.35 (0.16), residues: 1224 helix: -4.96 (0.08), residues: 637 sheet: -3.98 (0.66), residues: 36 loop : -3.12 (0.23), residues: 551 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 217 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 208 time to evaluate : 1.076 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 9 outliers final: 2 residues processed: 214 average time/residue: 0.2646 time to fit residues: 76.0917 Evaluate side-chains 124 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 122 time to evaluate : 1.071 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0920 time to fit residues: 1.9510 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 103 optimal weight: 9.9990 chunk 93 optimal weight: 0.8980 chunk 51 optimal weight: 10.0000 chunk 31 optimal weight: 0.7980 chunk 62 optimal weight: 9.9990 chunk 49 optimal weight: 7.9990 chunk 96 optimal weight: 0.7980 chunk 37 optimal weight: 7.9990 chunk 58 optimal weight: 1.9990 chunk 71 optimal weight: 0.6980 chunk 111 optimal weight: 9.9990 overall best weight: 1.0382 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 170 GLN U 207 ASN U 466 GLN U 492 GLN ** U 515 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 536 ASN U 546 GLN U 576 ASN U 605 GLN U 640 HIS U 697 GLN U 708 HIS U 778 GLN U 791 GLN U 842 GLN U 875 GLN ** U 938 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U1097 ASN ** U1124 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 16 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7731 moved from start: 0.2605 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.041 9761 Z= 0.220 Angle : 0.726 8.690 13199 Z= 0.373 Chirality : 0.044 0.243 1502 Planarity : 0.004 0.052 1689 Dihedral : 7.003 38.284 1327 Min Nonbonded Distance : 2.179 Molprobity Statistics. All-atom Clashscore : 11.19 Ramachandran Plot: Outliers : 0.08 % Allowed : 8.74 % Favored : 91.18 % Rotamer Outliers : 3.51 % 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: -5.50 (0.19), residues: 1224 helix: -4.36 (0.12), residues: 637 sheet: -3.25 (0.70), residues: 36 loop : -2.59 (0.25), residues: 551 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 180 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 36 poor density : 144 time to evaluate : 1.207 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 36 outliers final: 12 residues processed: 171 average time/residue: 0.1881 time to fit residues: 46.6257 Evaluate side-chains 126 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 114 time to evaluate : 0.962 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 0 residues processed: 12 average time/residue: 0.0848 time to fit residues: 3.6134 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 62 optimal weight: 2.9990 chunk 34 optimal weight: 3.9990 chunk 92 optimal weight: 0.9980 chunk 76 optimal weight: 0.5980 chunk 30 optimal weight: 0.9980 chunk 111 optimal weight: 7.9990 chunk 120 optimal weight: 30.0000 chunk 99 optimal weight: 0.0070 chunk 110 optimal weight: 8.9990 chunk 38 optimal weight: 4.9990 chunk 89 optimal weight: 2.9990 overall best weight: 1.1200 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 515 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 546 GLN ** U 938 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U1124 GLN U1258 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7753 moved from start: 0.3104 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.039 9761 Z= 0.224 Angle : 0.693 9.458 13199 Z= 0.349 Chirality : 0.043 0.156 1502 Planarity : 0.004 0.044 1689 Dihedral : 6.340 36.555 1327 Min Nonbonded Distance : 2.188 Molprobity Statistics. All-atom Clashscore : 11.45 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.42 % Favored : 91.58 % Rotamer Outliers : 3.22 % 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: -5.17 (0.19), residues: 1224 helix: -3.96 (0.14), residues: 664 sheet: -2.51 (0.68), residues: 43 loop : -2.53 (0.26), residues: 517 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 162 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 33 poor density : 129 time to evaluate : 1.173 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 33 outliers final: 11 residues processed: 155 average time/residue: 0.1939 time to fit residues: 44.6585 Evaluate side-chains 117 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 106 time to evaluate : 0.982 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.0973 time to fit residues: 3.4854 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 110 optimal weight: 6.9990 chunk 84 optimal weight: 4.9990 chunk 58 optimal weight: 20.0000 chunk 12 optimal weight: 0.0770 chunk 53 optimal weight: 3.9990 chunk 75 optimal weight: 3.9990 chunk 112 optimal weight: 0.3980 chunk 118 optimal weight: 10.0000 chunk 106 optimal weight: 10.0000 chunk 32 optimal weight: 0.0470 chunk 99 optimal weight: 7.9990 overall best weight: 1.7040 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 765 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U1258 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7799 moved from start: 0.3481 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.053 9761 Z= 0.290 Angle : 0.695 10.723 13199 Z= 0.348 Chirality : 0.045 0.153 1502 Planarity : 0.004 0.041 1689 Dihedral : 6.227 37.246 1327 Min Nonbonded Distance : 2.180 Molprobity Statistics. All-atom Clashscore : 11.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.66 % Favored : 91.34 % Rotamer Outliers : 2.83 % 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: -4.81 (0.21), residues: 1224 helix: -3.75 (0.15), residues: 657 sheet: -2.07 (0.73), residues: 43 loop : -2.29 (0.27), residues: 524 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 143 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 29 poor density : 114 time to evaluate : 1.125 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 29 outliers final: 10 residues processed: 136 average time/residue: 0.1989 time to fit residues: 40.2006 Evaluate side-chains 116 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 106 time to evaluate : 1.067 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.1008 time to fit residues: 3.4581 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 67 optimal weight: 3.9990 chunk 1 optimal weight: 0.6980 chunk 88 optimal weight: 0.9980 chunk 49 optimal weight: 8.9990 chunk 101 optimal weight: 4.9990 chunk 82 optimal weight: 1.9990 chunk 0 optimal weight: 2.9990 chunk 60 optimal weight: 20.0000 chunk 106 optimal weight: 5.9990 chunk 30 optimal weight: 0.9990 chunk 40 optimal weight: 7.9990 overall best weight: 1.5386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 765 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U1258 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7797 moved from start: 0.3769 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.051 9761 Z= 0.265 Angle : 0.681 10.240 13199 Z= 0.338 Chirality : 0.044 0.159 1502 Planarity : 0.004 0.037 1689 Dihedral : 6.034 37.008 1327 Min Nonbonded Distance : 2.195 Molprobity Statistics. All-atom Clashscore : 11.29 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.25 % Favored : 91.75 % Rotamer Outliers : 2.53 % 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: -4.65 (0.21), residues: 1224 helix: -3.63 (0.15), residues: 663 sheet: -1.83 (0.75), residues: 43 loop : -2.20 (0.27), residues: 518 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 146 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 26 poor density : 120 time to evaluate : 1.130 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 26 outliers final: 6 residues processed: 138 average time/residue: 0.1887 time to fit residues: 39.0637 Evaluate side-chains 116 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 110 time to evaluate : 1.053 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.0970 time to fit residues: 2.3981 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 107 optimal weight: 0.3980 chunk 23 optimal weight: 1.9990 chunk 69 optimal weight: 0.7980 chunk 29 optimal weight: 0.7980 chunk 119 optimal weight: 0.9990 chunk 98 optimal weight: 30.0000 chunk 55 optimal weight: 5.9990 chunk 9 optimal weight: 4.9990 chunk 39 optimal weight: 8.9990 chunk 62 optimal weight: 9.9990 chunk 114 optimal weight: 0.9980 overall best weight: 0.7982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 466 GLN ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U 765 ASN U1196 GLN U1258 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7751 moved from start: 0.4005 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.035 9761 Z= 0.180 Angle : 0.634 8.871 13199 Z= 0.313 Chirality : 0.043 0.162 1502 Planarity : 0.003 0.036 1689 Dihedral : 5.644 34.879 1327 Min Nonbonded Distance : 2.183 Molprobity Statistics. All-atom Clashscore : 11.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.35 % Favored : 92.65 % Rotamer Outliers : 1.75 % 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: -4.46 (0.21), residues: 1224 helix: -3.45 (0.16), residues: 673 sheet: -1.54 (0.77), residues: 43 loop : -2.13 (0.27), residues: 508 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 137 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 119 time to evaluate : 1.046 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 18 outliers final: 7 residues processed: 133 average time/residue: 0.2098 time to fit residues: 41.0805 Evaluate side-chains 113 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 106 time to evaluate : 1.032 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0963 time to fit residues: 2.8091 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 13 optimal weight: 1.9990 chunk 67 optimal weight: 3.9990 chunk 86 optimal weight: 3.9990 chunk 100 optimal weight: 10.0000 chunk 66 optimal weight: 2.9990 chunk 118 optimal weight: 0.7980 chunk 74 optimal weight: 0.5980 chunk 72 optimal weight: 0.0570 chunk 54 optimal weight: 20.0000 chunk 73 optimal weight: 0.6980 chunk 47 optimal weight: 9.9990 overall best weight: 0.8300 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 765 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 813 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7737 moved from start: 0.4204 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.038 9761 Z= 0.186 Angle : 0.646 11.723 13199 Z= 0.314 Chirality : 0.043 0.168 1502 Planarity : 0.003 0.035 1689 Dihedral : 5.448 33.444 1327 Min Nonbonded Distance : 2.193 Molprobity Statistics. All-atom Clashscore : 10.67 Ramachandran Plot: Outliers : 0.08 % Allowed : 7.19 % Favored : 92.73 % Rotamer Outliers : 0.88 % 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: -4.28 (0.22), residues: 1224 helix: -3.28 (0.16), residues: 681 sheet: -1.25 (0.79), residues: 43 loop : -2.05 (0.28), residues: 500 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 129 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 120 time to evaluate : 1.046 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 9 outliers final: 2 residues processed: 128 average time/residue: 0.2171 time to fit residues: 40.6256 Evaluate side-chains 116 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 114 time to evaluate : 1.085 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0881 time to fit residues: 1.8656 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 70 optimal weight: 2.9990 chunk 35 optimal weight: 0.9980 chunk 23 optimal weight: 0.6980 chunk 22 optimal weight: 1.9990 chunk 75 optimal weight: 4.9990 chunk 80 optimal weight: 2.9990 chunk 58 optimal weight: 30.0000 chunk 11 optimal weight: 0.9980 chunk 93 optimal weight: 0.7980 chunk 107 optimal weight: 9.9990 chunk 113 optimal weight: 7.9990 overall best weight: 1.0982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 375 ASN U 442 ASN ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7753 moved from start: 0.4371 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.041 9761 Z= 0.213 Angle : 0.665 11.633 13199 Z= 0.319 Chirality : 0.043 0.192 1502 Planarity : 0.003 0.035 1689 Dihedral : 5.415 34.023 1327 Min Nonbonded Distance : 2.196 Molprobity Statistics. All-atom Clashscore : 10.30 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.01 % Favored : 91.99 % Rotamer Outliers : 1.46 % 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: -4.12 (0.22), residues: 1224 helix: -3.17 (0.16), residues: 682 sheet: -0.97 (0.82), residues: 43 loop : -1.97 (0.28), residues: 499 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 126 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 111 time to evaluate : 1.195 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 15 outliers final: 11 residues processed: 120 average time/residue: 0.2060 time to fit residues: 36.7530 Evaluate side-chains 118 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 107 time to evaluate : 1.108 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.0937 time to fit residues: 3.4898 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 103 optimal weight: 0.0170 chunk 110 optimal weight: 9.9990 chunk 113 optimal weight: 0.0060 chunk 66 optimal weight: 0.9990 chunk 48 optimal weight: 9.9990 chunk 86 optimal weight: 4.9990 chunk 33 optimal weight: 0.7980 chunk 99 optimal weight: 30.0000 chunk 104 optimal weight: 0.0980 chunk 72 optimal weight: 2.9990 chunk 116 optimal weight: 3.9990 overall best weight: 0.3836 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 172 GLN U 424 ASN U 442 ASN ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U 576 ASN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7682 moved from start: 0.4626 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.041 9761 Z= 0.147 Angle : 0.643 11.686 13199 Z= 0.308 Chirality : 0.042 0.177 1502 Planarity : 0.003 0.038 1689 Dihedral : 5.062 31.511 1327 Min Nonbonded Distance : 2.206 Molprobity Statistics. All-atom Clashscore : 9.78 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.70 % Favored : 93.30 % Rotamer Outliers : 0.39 % 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: -3.91 (0.22), residues: 1224 helix: -2.93 (0.17), residues: 688 sheet: -0.42 (0.93), residues: 38 loop : -2.01 (0.27), residues: 498 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 122 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 118 time to evaluate : 1.201 Fit side-chains revert: symmetry clash outliers start: 4 outliers final: 1 residues processed: 122 average time/residue: 0.2145 time to fit residues: 37.9251 Evaluate side-chains 114 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 113 time to evaluate : 1.085 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0889 time to fit residues: 1.6449 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 71 optimal weight: 0.9980 chunk 55 optimal weight: 9.9990 chunk 81 optimal weight: 0.9980 chunk 122 optimal weight: 20.0000 chunk 112 optimal weight: 0.5980 chunk 97 optimal weight: 0.0870 chunk 10 optimal weight: 2.9990 chunk 75 optimal weight: 0.5980 chunk 59 optimal weight: 5.9990 chunk 77 optimal weight: 3.9990 chunk 103 optimal weight: 0.0670 overall best weight: 0.4696 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 442 ASN ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U 576 ASN U1258 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7673 moved from start: 0.4816 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.040 9761 Z= 0.155 Angle : 0.649 12.566 13199 Z= 0.311 Chirality : 0.042 0.171 1502 Planarity : 0.003 0.035 1689 Dihedral : 4.924 30.927 1327 Min Nonbonded Distance : 2.192 Molprobity Statistics. All-atom Clashscore : 9.52 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.54 % Favored : 93.46 % Rotamer Outliers : 0.29 % 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: -3.64 (0.23), residues: 1224 helix: -2.75 (0.18), residues: 679 sheet: -0.70 (0.84), residues: 43 loop : -1.83 (0.28), residues: 502 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2448 Ramachandran restraints generated. 1224 Oldfield, 0 Emsley, 1224 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 115 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 112 time to evaluate : 1.004 Fit side-chains revert: symmetry clash outliers start: 3 outliers final: 0 residues processed: 114 average time/residue: 0.2228 time to fit residues: 36.9361 Evaluate side-chains 107 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 107 time to evaluate : 1.078 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4814 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 123 random chunks: chunk 29 optimal weight: 0.9990 chunk 89 optimal weight: 1.9990 chunk 14 optimal weight: 2.9990 chunk 27 optimal weight: 2.9990 chunk 97 optimal weight: 0.1980 chunk 40 optimal weight: 2.9990 chunk 100 optimal weight: 6.9990 chunk 12 optimal weight: 0.4980 chunk 18 optimal weight: 0.0570 chunk 85 optimal weight: 0.9990 chunk 5 optimal weight: 0.0050 overall best weight: 0.3514 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** U 116 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** U 466 GLN ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U1258 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4476 r_free = 0.4476 target = 0.224029 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 64)----------------| | r_work = 0.4080 r_free = 0.4080 target = 0.180354 restraints weight = 10795.958| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 44)----------------| | r_work = 0.4036 r_free = 0.4036 target = 0.176789 restraints weight = 18099.613| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 52)----------------| | r_work = 0.4001 r_free = 0.4001 target = 0.173241 restraints weight = 16366.774| |-----------------------------------------------------------------------------| r_work (final): 0.4000 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7474 moved from start: 0.5016 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.042 9761 Z= 0.150 Angle : 0.646 12.203 13199 Z= 0.310 Chirality : 0.042 0.167 1502 Planarity : 0.003 0.040 1689 Dihedral : 4.774 29.472 1327 Min Nonbonded Distance : 2.162 Molprobity Statistics. All-atom Clashscore : 9.47 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.54 % Favored : 93.46 % Rotamer Outliers : 0.39 % 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: -3.41 (0.23), residues: 1224 helix: -2.57 (0.18), residues: 671 sheet: -0.74 (0.84), residues: 43 loop : -1.74 (0.27), residues: 510 =============================================================================== Job complete usr+sys time: 1939.45 seconds wall clock time: 36 minutes 15.00 seconds (2175.00 seconds total)