Starting phenix.real_space_refine (version: dev) on Mon Dec 12 07:52:10 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/12_2022/6lr0_31583.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/12_2022/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/12_2022/6lr0_31583.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/12_2022/6lr0_31583.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/12_2022/6lr0_31583.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6lr0_31583/12_2022/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 Convert atoms to be neutral 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-4807/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.78, per 1000 atoms: 0.60 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.44 Conformation dependent library (CDL) restraints added in 1.4 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.60 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.62 Time building geometry restraints manager: 3.74 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.740 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.960 Check model and map are aligned: 0.150 Convert atoms to be neutral: 0.090 Process input model: 27.260 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.640 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 34.170 ------------------------------------------------------------------------------- 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.214 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 9 outliers final: 2 residues processed: 214 average time/residue: 0.2499 time to fit residues: 72.0565 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.101 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0907 time to fit residues: 1.8854 ------------------------------------------------------------------------------- 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.020 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 36 outliers final: 12 residues processed: 171 average time/residue: 0.1917 time to fit residues: 47.7430 Evaluate side-chains 126 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 114 time to evaluate : 1.160 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 0 residues processed: 12 average time/residue: 0.0940 time to fit residues: 3.7038 ------------------------------------------------------------------------------- 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: 4.9990 chunk 34 optimal weight: 3.9990 chunk 92 optimal weight: 0.9980 chunk 76 optimal weight: 0.8980 chunk 30 optimal weight: 0.9990 chunk 111 optimal weight: 0.5980 chunk 120 optimal weight: 30.0000 chunk 99 optimal weight: 0.1980 chunk 110 optimal weight: 2.9990 chunk 38 optimal weight: 5.9990 chunk 89 optimal weight: 2.9990 overall best weight: 0.7382 ------------------------------------------------------------------------------- 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 938 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U1124 GLN U1258 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7717 moved from start: 0.3168 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.036 9761 Z= 0.185 Angle : 0.674 9.505 13199 Z= 0.340 Chirality : 0.043 0.156 1502 Planarity : 0.004 0.042 1689 Dihedral : 6.193 35.420 1327 Min Nonbonded Distance : 2.180 Molprobity Statistics. All-atom Clashscore : 11.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.76 % Favored : 92.24 % Rotamer Outliers : 2.63 % 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.13 (0.19), residues: 1224 helix: -3.91 (0.14), residues: 659 sheet: -2.53 (0.67), residues: 43 loop : -2.56 (0.25), residues: 522 *********************** 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 156 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 129 time to evaluate : 1.163 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 27 outliers final: 8 residues processed: 149 average time/residue: 0.1953 time to fit residues: 43.2346 Evaluate side-chains 111 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 103 time to evaluate : 1.108 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0962 time to fit residues: 2.9512 ------------------------------------------------------------------------------- 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: 20.0000 chunk 84 optimal weight: 6.9990 chunk 58 optimal weight: 20.0000 chunk 12 optimal weight: 0.7980 chunk 53 optimal weight: 3.9990 chunk 75 optimal weight: 3.9990 chunk 112 optimal weight: 10.0000 chunk 118 optimal weight: 9.9990 chunk 106 optimal weight: 2.9990 chunk 32 optimal weight: 2.9990 chunk 99 optimal weight: 9.9990 overall best weight: 2.9588 ------------------------------------------------------------------------------- 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 361 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U 546 GLN ** 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: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7868 moved from start: 0.3449 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.007 0.077 9761 Z= 0.452 Angle : 0.773 13.080 13199 Z= 0.385 Chirality : 0.048 0.188 1502 Planarity : 0.004 0.048 1689 Dihedral : 6.651 39.462 1327 Min Nonbonded Distance : 2.175 Molprobity Statistics. All-atom Clashscore : 12.75 Ramachandran Plot: Outliers : 0.08 % Allowed : 9.72 % Favored : 90.20 % Rotamer Outliers : 3.90 % 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.91 (0.21), residues: 1224 helix: -3.83 (0.15), residues: 654 sheet: -2.15 (0.72), residues: 43 loop : -2.36 (0.27), residues: 527 *********************** 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 150 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 40 poor density : 110 time to evaluate : 1.170 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 40 outliers final: 15 residues processed: 142 average time/residue: 0.2262 time to fit residues: 46.7444 Evaluate side-chains 121 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 106 time to evaluate : 1.154 Switching outliers to nearest non-outliers outliers start: 15 outliers final: 0 residues processed: 15 average time/residue: 0.0964 time to fit residues: 4.2325 ------------------------------------------------------------------------------- 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: 0.9980 chunk 1 optimal weight: 0.0980 chunk 88 optimal weight: 0.7980 chunk 49 optimal weight: 7.9990 chunk 101 optimal weight: 10.0000 chunk 82 optimal weight: 0.5980 chunk 0 optimal weight: 0.9990 chunk 60 optimal weight: 20.0000 chunk 106 optimal weight: 0.6980 chunk 30 optimal weight: 0.6980 chunk 40 optimal weight: 0.2980 overall best weight: 0.4780 ------------------------------------------------------------------------------- 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 U1258 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7736 moved from start: 0.3834 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.033 9761 Z= 0.157 Angle : 0.656 9.494 13199 Z= 0.324 Chirality : 0.042 0.166 1502 Planarity : 0.003 0.037 1689 Dihedral : 5.802 35.615 1327 Min Nonbonded Distance : 2.169 Molprobity Statistics. All-atom Clashscore : 11.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.43 % Favored : 92.57 % Rotamer Outliers : 1.85 % 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.62 (0.21), residues: 1224 helix: -3.59 (0.15), residues: 665 sheet: -1.75 (0.75), residues: 43 loop : -2.20 (0.27), residues: 516 *********************** 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 141 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 19 poor density : 122 time to evaluate : 1.116 Fit side-chains revert: symmetry clash outliers start: 19 outliers final: 7 residues processed: 136 average time/residue: 0.1984 time to fit residues: 39.3594 Evaluate side-chains 115 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 108 time to evaluate : 1.125 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0961 time to fit residues: 2.7591 ------------------------------------------------------------------------------- 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: 20.0000 chunk 23 optimal weight: 1.9990 chunk 69 optimal weight: 0.8980 chunk 29 optimal weight: 2.9990 chunk 119 optimal weight: 0.0070 chunk 98 optimal weight: 40.0000 chunk 55 optimal weight: 4.9990 chunk 9 optimal weight: 5.9990 chunk 39 optimal weight: 0.0020 chunk 62 optimal weight: 0.0970 chunk 114 optimal weight: 0.7980 overall best weight: 0.3604 ------------------------------------------------------------------------------- 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 765 ASN U 813 GLN U1258 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7694 moved from start: 0.4149 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.033 9761 Z= 0.140 Angle : 0.626 9.578 13199 Z= 0.306 Chirality : 0.042 0.165 1502 Planarity : 0.003 0.036 1689 Dihedral : 5.366 31.972 1327 Min Nonbonded Distance : 2.178 Molprobity Statistics. All-atom Clashscore : 10.35 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.27 % Favored : 92.73 % Rotamer Outliers : 1.17 % 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.45 (0.21), residues: 1224 helix: -3.36 (0.16), residues: 681 sheet: -1.76 (0.75), residues: 43 loop : -2.18 (0.27), residues: 500 *********************** 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 139 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 127 time to evaluate : 1.089 Fit side-chains revert: symmetry clash outliers start: 12 outliers final: 6 residues processed: 135 average time/residue: 0.2103 time to fit residues: 40.6294 Evaluate side-chains 115 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 109 time to evaluate : 1.050 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.0936 time to fit residues: 2.5590 ------------------------------------------------------------------------------- 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: 4.9990 chunk 86 optimal weight: 2.9990 chunk 100 optimal weight: 7.9990 chunk 66 optimal weight: 1.9990 chunk 118 optimal weight: 0.9990 chunk 74 optimal weight: 0.3980 chunk 72 optimal weight: 3.9990 chunk 54 optimal weight: 9.9990 chunk 73 optimal weight: 0.9980 chunk 47 optimal weight: 8.9990 overall best weight: 1.2786 ------------------------------------------------------------------------------- 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 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7760 moved from start: 0.4243 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.041 9761 Z= 0.235 Angle : 0.657 8.621 13199 Z= 0.319 Chirality : 0.044 0.157 1502 Planarity : 0.003 0.036 1689 Dihedral : 5.497 34.133 1327 Min Nonbonded Distance : 2.203 Molprobity Statistics. All-atom Clashscore : 10.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.60 % Favored : 92.40 % Rotamer Outliers : 1.07 % 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.25 (0.22), residues: 1224 helix: -3.26 (0.16), residues: 681 sheet: -1.30 (0.84), residues: 38 loop : -2.03 (0.28), residues: 505 *********************** 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 126 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 115 time to evaluate : 1.168 Fit side-chains revert: symmetry clash outliers start: 11 outliers final: 5 residues processed: 123 average time/residue: 0.2125 time to fit residues: 37.7724 Evaluate side-chains 116 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 111 time to evaluate : 0.996 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0838 time to fit residues: 2.2044 ------------------------------------------------------------------------------- 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.8980 chunk 22 optimal weight: 2.9990 chunk 75 optimal weight: 3.9990 chunk 80 optimal weight: 2.9990 chunk 58 optimal weight: 10.0000 chunk 11 optimal weight: 0.6980 chunk 93 optimal weight: 0.6980 chunk 107 optimal weight: 0.0170 chunk 113 optimal weight: 0.0670 overall best weight: 0.4756 ------------------------------------------------------------------------------- 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 Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7691 moved from start: 0.4549 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.153 Angle : 0.639 11.659 13199 Z= 0.306 Chirality : 0.042 0.198 1502 Planarity : 0.003 0.035 1689 Dihedral : 5.136 31.535 1327 Min Nonbonded Distance : 2.207 Molprobity Statistics. All-atom Clashscore : 9.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.78 % Favored : 93.22 % Rotamer Outliers : 1.07 % 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.06 (0.22), residues: 1224 helix: -3.11 (0.16), residues: 683 sheet: -0.85 (0.87), residues: 38 loop : -1.97 (0.28), residues: 503 *********************** 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 128 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 117 time to evaluate : 1.192 Fit side-chains revert: symmetry clash outliers start: 11 outliers final: 5 residues processed: 126 average time/residue: 0.1964 time to fit residues: 36.6310 Evaluate side-chains 112 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 107 time to evaluate : 1.104 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0947 time to fit residues: 2.4272 ------------------------------------------------------------------------------- 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.0470 chunk 110 optimal weight: 8.9990 chunk 113 optimal weight: 0.0000 chunk 66 optimal weight: 4.9990 chunk 48 optimal weight: 4.9990 chunk 86 optimal weight: 3.9990 chunk 33 optimal weight: 0.6980 chunk 99 optimal weight: 0.0870 chunk 104 optimal weight: 3.9990 chunk 72 optimal weight: 2.9990 chunk 116 optimal weight: 0.0770 overall best weight: 0.1818 ------------------------------------------------------------------------------- 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 424 ASN U 442 ASN U 466 GLN ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U1258 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7625 moved from start: 0.4867 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.038 9761 Z= 0.134 Angle : 0.631 11.272 13199 Z= 0.305 Chirality : 0.041 0.182 1502 Planarity : 0.003 0.041 1689 Dihedral : 4.815 30.001 1327 Min Nonbonded Distance : 2.219 Molprobity Statistics. All-atom Clashscore : 9.52 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.29 % Favored : 93.71 % Rotamer Outliers : 0.49 % 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.22), residues: 1224 helix: -2.77 (0.17), residues: 673 sheet: -0.76 (0.87), residues: 38 loop : -1.80 (0.27), residues: 513 *********************** 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 130 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 125 time to evaluate : 1.146 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 5 outliers final: 1 residues processed: 127 average time/residue: 0.1948 time to fit residues: 36.0892 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.053 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1049 time to fit residues: 1.7033 ------------------------------------------------------------------------------- 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.6980 chunk 55 optimal weight: 1.9990 chunk 81 optimal weight: 0.5980 chunk 122 optimal weight: 20.0000 chunk 112 optimal weight: 8.9990 chunk 97 optimal weight: 0.0980 chunk 10 optimal weight: 2.9990 chunk 75 optimal weight: 0.6980 chunk 59 optimal weight: 5.9990 chunk 77 optimal weight: 2.9990 chunk 103 optimal weight: 0.0030 overall best weight: 0.4190 ------------------------------------------------------------------------------- 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 492 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U 697 GLN U1258 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7634 moved from start: 0.5005 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.035 9761 Z= 0.152 Angle : 0.651 12.192 13199 Z= 0.313 Chirality : 0.042 0.174 1502 Planarity : 0.003 0.035 1689 Dihedral : 4.745 29.599 1327 Min Nonbonded Distance : 2.169 Molprobity Statistics. All-atom Clashscore : 9.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.54 % Favored : 93.46 % Rotamer Outliers : 0.10 % 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.44 (0.22), residues: 1224 helix: -2.64 (0.18), residues: 662 sheet: -0.58 (0.90), residues: 38 loop : -1.71 (0.27), residues: 524 ********************** 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 118 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 117 time to evaluate : 1.049 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 117 average time/residue: 0.2206 time to fit residues: 37.5533 Evaluate side-chains 109 residues out of total 1026 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 109 time to evaluate : 1.145 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.6645 ------------------------------------------------------------------------------- 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: 3.9990 chunk 27 optimal weight: 0.8980 chunk 97 optimal weight: 0.5980 chunk 40 optimal weight: 7.9990 chunk 100 optimal weight: 3.9990 chunk 12 optimal weight: 0.9980 chunk 18 optimal weight: 0.4980 chunk 85 optimal weight: 0.4980 chunk 5 optimal weight: 0.6980 overall best weight: 0.6380 ------------------------------------------------------------------------------- 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 492 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** U 563 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** U 867 GLN U1258 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4468 r_free = 0.4468 target = 0.223178 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 68)----------------| | r_work = 0.4067 r_free = 0.4067 target = 0.179329 restraints weight = 10756.395| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 35)----------------| | r_work = 0.4023 r_free = 0.4023 target = 0.175495 restraints weight = 18920.440| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 36)----------------| | r_work = 0.4011 r_free = 0.4011 target = 0.174349 restraints weight = 18171.238| |-----------------------------------------------------------------------------| r_work (final): 0.3729 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7413 moved from start: 0.5069 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.037 9761 Z= 0.167 Angle : 0.642 12.686 13199 Z= 0.310 Chirality : 0.042 0.165 1502 Planarity : 0.003 0.038 1689 Dihedral : 4.747 29.487 1327 Min Nonbonded Distance : 2.140 Molprobity Statistics. All-atom Clashscore : 9.89 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.29 % Favored : 93.71 % Rotamer Outliers : 0.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: -3.34 (0.23), residues: 1224 helix: -2.51 (0.18), residues: 670 sheet: -0.18 (0.90), residues: 36 loop : -1.73 (0.27), residues: 518 =============================================================================== Job complete usr+sys time: 1884.75 seconds wall clock time: 35 minutes 2.74 seconds (2102.74 seconds total)