Starting phenix.refine on Mon May 20 17:17:37 2024 by afonine =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/finalised/1fp9.pdb Found miller_array, /net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/original/1fp9/1fp9.mtz Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- optimize_xyz=true optimize_adp=true main.number_of_mac=3 main.max_number_of_iterations=100 main.nproc=20 nonbonded_weight=100 modify_start_model.modify.selection=element H modify_start_model.modify.occupancies.set=1 strategy=individual_sites+individual_adp ramachandran_plot_restraints.enabled=False secondary_structure.enabled=True main.nqh_flips=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { model { file = "/net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/finalised/1fp9.pdb" } default_model = "/net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/finalised/1fp9.pdb" miller_array { file = "/net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/original/1fp9/1fp9.mtz" } default_miller_array = "/net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/original/1fp9/1fp9.mtz" } refinement { electron_density_maps { map_coefficients { map_type = 2mFo-DFc mtz_label_amplitudes = 2FOFCWT mtz_label_phases = PH2FOFCWT fill_missing_f_obs = True } map_coefficients { map_type = 2mFo-DFc mtz_label_amplitudes = 2FOFCWT_no_fill mtz_label_phases = PH2FOFCWT_no_fill } map_coefficients { map_type = mFo-DFc mtz_label_amplitudes = FOFCWT mtz_label_phases = PHFOFCWT } map_coefficients { map_type = anomalous mtz_label_amplitudes = ANOM mtz_label_phases = PANOM } map { map_type = 2mFo-DFc fill_missing_f_obs = True } map { map_type = 2mFo-DFc } map { map_type = mFo-DFc } } refine { strategy = *individual_sites individual_sites_real_space rigid_body \ *individual_adp group_adp tls occupancies group_anomalous \ den } main { nqh_flips = False max_number_of_iterations = 100 nproc = 20 } modify_start_model { modify { selection = element H occupancies { set = 1 } } } pdb_interpretation { secondary_structure { enabled = True } nonbonded_weight = 100 } target_weights { optimize_xyz_weight = true optimize_adp_weight = true } } Starting job =============================================================================== Symmetric amino acids flipped None Time to flip residues: 0.01s Monomer Library directory: "/net/anaconda/raid1/afonine/phenix3/modules/chem_data/mon_lib" Total number of atoms: 7982 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 7982 Number of conformers: 1 Conformer: "" Number of residues, atoms: 500, 7982 Classifications: {'peptide': 500} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 38, 'TRANS': 461} Time building chain proxies: 2.55, per 1000 atoms: 0.32 Number of scatterers: 7982 At special positions: 0 Unit cell: (104.913, 52.445, 104.908, 90, 96.42, 90) Space group: C 1 2 1 (No. 5) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 11 16.00 O 707 8.00 N 703 7.00 C 2641 6.00 H 3920 1.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: 7.23 Conformation dependent library (CDL) restraints added in 998.0 milliseconds Adding C-beta torsion restraints... Number of C-beta restraints generated: 910 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 26 helices and 3 sheets defined 45.2% alpha, 6.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.83 Creating SS restraints... Processing helix chain 'A' and resid 13 through 15 No H-bonds generated for 'chain 'A' and resid 13 through 15' Processing helix chain 'A' and resid 27 through 38 Processing helix chain 'A' and resid 69 through 71 No H-bonds generated for 'chain 'A' and resid 69 through 71' Processing helix chain 'A' and resid 75 through 80 removed outlier: 3.994A pdb=" N ARG A 80 " --> pdb=" O PRO A 76 " (cutoff:3.500A) Processing helix chain 'A' and resid 97 through 118 removed outlier: 3.803A pdb=" N TRP A 105 " --> pdb=" O TYR A 101 " (cutoff:3.500A) Proline residue: A 106 - end of helix removed outlier: 3.572A pdb=" N GLU A 117 " --> pdb=" O ARG A 113 " (cutoff:3.500A) removed outlier: 3.605A pdb=" N LYS A 118 " --> pdb=" O GLY A 114 " (cutoff:3.500A) Processing helix chain 'A' and resid 123 through 133 Processing helix chain 'A' and resid 135 through 150 removed outlier: 4.644A pdb=" N ASP A 140 " --> pdb=" O TRP A 136 " (cutoff:3.500A) Processing helix chain 'A' and resid 156 through 158 No H-bonds generated for 'chain 'A' and resid 156 through 158' Processing helix chain 'A' and resid 161 through 164 No H-bonds generated for 'chain 'A' and resid 161 through 164' Processing helix chain 'A' and resid 168 through 177 Processing helix chain 'A' and resid 179 through 205 Processing helix chain 'A' and resid 223 through 227 Processing helix chain 'A' and resid 229 through 231 No H-bonds generated for 'chain 'A' and resid 229 through 231' Processing helix chain 'A' and resid 265 through 270 Processing helix chain 'A' and resid 274 through 286 Processing helix chain 'A' and resid 295 through 298 Processing helix chain 'A' and resid 321 through 332 Processing helix chain 'A' and resid 347 through 356 Processing helix chain 'A' and resid 364 through 366 No H-bonds generated for 'chain 'A' and resid 364 through 366' Processing helix chain 'A' and resid 378 through 380 No H-bonds generated for 'chain 'A' and resid 378 through 380' Processing helix chain 'A' and resid 399 through 403 Processing helix chain 'A' and resid 408 through 421 removed outlier: 3.814A pdb=" N ASP A 420 " --> pdb=" O ARG A 416 " (cutoff:3.500A) Processing helix chain 'A' and resid 431 through 441 Processing helix chain 'A' and resid 451 through 454 No H-bonds generated for 'chain 'A' and resid 451 through 454' Processing helix chain 'A' and resid 459 through 461 No H-bonds generated for 'chain 'A' and resid 459 through 461' Processing helix chain 'A' and resid 484 through 496 removed outlier: 3.687A pdb=" N MET A 493 " --> pdb=" O ARG A 489 " (cutoff:3.500A) Processing sheet with id= A, first strand: chain 'A' and resid 446 through 450 removed outlier: 7.402A pdb=" N ALA A 6 " --> pdb=" O ALA A 447 " (cutoff:3.500A) removed outlier: 8.223A pdb=" N TYR A 449 " --> pdb=" O ALA A 6 " (cutoff:3.500A) removed outlier: 6.515A pdb=" N GLY A 8 " --> pdb=" O TYR A 449 " (cutoff:3.500A) removed outlier: 6.857A pdb=" N TYR A 43 " --> pdb=" O LEU A 9 " (cutoff:3.500A) removed outlier: 6.575A pdb=" N ARG A 209 " --> pdb=" O TRP A 44 " (cutoff:3.500A) removed outlier: 6.325A pdb=" N LEU A 289 " --> pdb=" O GLY A 212 " (cutoff:3.500A) removed outlier: 7.967A pdb=" N MET A 214 " --> pdb=" O LEU A 289 " (cutoff:3.500A) removed outlier: 6.170A pdb=" N ARG A 291 " --> pdb=" O MET A 214 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 300 through 303 Processing sheet with id= C, first strand: chain 'A' and resid 360 through 363 129 hydrogen bonds defined for protein. 381 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: 5.88 Time building geometry restraints manager: 4.30 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.98 - 1.15: 3919 1.15 - 1.32: 661 1.32 - 1.48: 1834 1.48 - 1.65: 1688 1.65 - 1.82: 21 Bond restraints: 8123 Sorted by residual: bond pdb=" NH1 ARG A 84 " pdb="HH11 ARG A 84 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.46e+01 bond pdb=" NH2 ARG A 476 " pdb="HH22 ARG A 476 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" N THR A 399 " pdb=" H THR A 399 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" NH2 ARG A 279 " pdb="HH22 ARG A 279 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" NE1 TRP A 159 " pdb=" HE1 TRP A 159 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 ... (remaining 8118 not shown) Histogram of bond angle deviations from ideal: 100.21 - 106.96: 414 106.96 - 113.72: 8838 113.72 - 120.48: 2791 120.48 - 127.23: 2538 127.23 - 133.99: 74 Bond angle restraints: 14655 Sorted by residual: angle pdb=" OE1 GLN A 27 " pdb=" CD GLN A 27 " pdb=" NE2 GLN A 27 " ideal model delta sigma weight residual 122.60 112.79 9.81 1.00e+00 1.00e+00 9.63e+01 angle pdb=" N THR A 391 " pdb=" CA THR A 391 " pdb=" C THR A 391 " ideal model delta sigma weight residual 111.36 102.64 8.72 1.09e+00 8.42e-01 6.40e+01 angle pdb=" CG GLN A 27 " pdb=" CD GLN A 27 " pdb=" NE2 GLN A 27 " ideal model delta sigma weight residual 116.40 126.40 -10.00 1.50e+00 4.44e-01 4.44e+01 angle pdb=" N GLY A 21 " pdb=" CA GLY A 21 " pdb=" C GLY A 21 " ideal model delta sigma weight residual 114.90 124.70 -9.80 1.54e+00 4.22e-01 4.05e+01 angle pdb=" N GLY A 114 " pdb=" CA GLY A 114 " pdb=" C GLY A 114 " ideal model delta sigma weight residual 114.67 108.55 6.12 1.10e+00 8.26e-01 3.10e+01 ... (remaining 14650 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 14.82: 3435 14.82 - 29.65: 193 29.65 - 44.47: 97 44.47 - 59.29: 91 59.29 - 74.12: 23 Dihedral angle restraints: 3839 sinusoidal: 2018 harmonic: 1821 Sorted by residual: dihedral pdb=" CG ARG A 84 " pdb=" CD ARG A 84 " pdb=" NE ARG A 84 " pdb=" CZ ARG A 84 " ideal model delta sinusoidal sigma weight residual 180.00 136.22 43.78 2 1.50e+01 4.44e-03 1.02e+01 dihedral pdb=" N GLU A 180 " pdb=" CA GLU A 180 " pdb=" CB GLU A 180 " pdb=" CG GLU A 180 " ideal model delta sinusoidal sigma weight residual 180.00 -129.53 -50.47 3 1.50e+01 4.44e-03 8.90e+00 dihedral pdb=" CA ARG A 296 " pdb=" CB ARG A 296 " pdb=" CG ARG A 296 " pdb=" CD ARG A 296 " ideal model delta sinusoidal sigma weight residual 60.00 109.93 -49.93 3 1.50e+01 4.44e-03 8.84e+00 ... (remaining 3836 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.056: 255 0.056 - 0.112: 233 0.112 - 0.168: 47 0.168 - 0.224: 23 0.224 - 0.280: 6 Chirality restraints: 564 Sorted by residual: chirality pdb=" CA PRO A 469 " pdb=" N PRO A 469 " pdb=" C PRO A 469 " pdb=" CB PRO A 469 " both_signs ideal model delta sigma weight residual False 2.72 2.44 0.28 2.00e-01 2.50e+01 1.97e+00 chirality pdb=" CA VAL A 349 " pdb=" N VAL A 349 " pdb=" C VAL A 349 " pdb=" CB VAL A 349 " both_signs ideal model delta sigma weight residual False 2.44 2.69 -0.25 2.00e-01 2.50e+01 1.50e+00 chirality pdb=" CA VAL A 182 " pdb=" N VAL A 182 " pdb=" C VAL A 182 " pdb=" CB VAL A 182 " both_signs ideal model delta sigma weight residual False 2.44 2.68 -0.24 2.00e-01 2.50e+01 1.40e+00 ... (remaining 561 not shown) Planarity restraints: 1211 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TYR A 141 " 0.051 2.00e-02 2.50e+03 2.06e-02 1.27e+01 pdb=" CG TYR A 141 " -0.010 2.00e-02 2.50e+03 pdb=" CD1 TYR A 141 " -0.011 2.00e-02 2.50e+03 pdb=" CD2 TYR A 141 " -0.012 2.00e-02 2.50e+03 pdb=" CE1 TYR A 141 " -0.006 2.00e-02 2.50e+03 pdb=" CE2 TYR A 141 " -0.005 2.00e-02 2.50e+03 pdb=" CZ TYR A 141 " 0.000 2.00e-02 2.50e+03 pdb=" OH TYR A 141 " 0.037 2.00e-02 2.50e+03 pdb=" HD1 TYR A 141 " -0.015 2.00e-02 2.50e+03 pdb=" HD2 TYR A 141 " -0.018 2.00e-02 2.50e+03 pdb=" HE1 TYR A 141 " -0.008 2.00e-02 2.50e+03 pdb=" HE2 TYR A 141 " -0.005 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR A 59 " -0.051 2.00e-02 2.50e+03 1.96e-02 1.15e+01 pdb=" CG TYR A 59 " 0.012 2.00e-02 2.50e+03 pdb=" CD1 TYR A 59 " 0.012 2.00e-02 2.50e+03 pdb=" CD2 TYR A 59 " 0.011 2.00e-02 2.50e+03 pdb=" CE1 TYR A 59 " 0.004 2.00e-02 2.50e+03 pdb=" CE2 TYR A 59 " 0.004 2.00e-02 2.50e+03 pdb=" CZ TYR A 59 " -0.002 2.00e-02 2.50e+03 pdb=" OH TYR A 59 " -0.032 2.00e-02 2.50e+03 pdb=" HD1 TYR A 59 " 0.017 2.00e-02 2.50e+03 pdb=" HD2 TYR A 59 " 0.016 2.00e-02 2.50e+03 pdb=" HE1 TYR A 59 " 0.004 2.00e-02 2.50e+03 pdb=" HE2 TYR A 59 " 0.004 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE A 192 " 0.052 2.00e-02 2.50e+03 1.93e-02 1.11e+01 pdb=" CG PHE A 192 " -0.016 2.00e-02 2.50e+03 pdb=" CD1 PHE A 192 " -0.012 2.00e-02 2.50e+03 pdb=" CD2 PHE A 192 " -0.013 2.00e-02 2.50e+03 pdb=" CE1 PHE A 192 " -0.001 2.00e-02 2.50e+03 pdb=" CE2 PHE A 192 " 0.000 2.00e-02 2.50e+03 pdb=" CZ PHE A 192 " 0.007 2.00e-02 2.50e+03 pdb=" HD1 PHE A 192 " -0.018 2.00e-02 2.50e+03 pdb=" HD2 PHE A 192 " -0.022 2.00e-02 2.50e+03 pdb=" HE1 PHE A 192 " 0.002 2.00e-02 2.50e+03 pdb=" HE2 PHE A 192 " 0.005 2.00e-02 2.50e+03 pdb=" HZ PHE A 192 " 0.017 2.00e-02 2.50e+03 ... (remaining 1208 not shown) Histogram of nonbonded interaction distances: 1.28 - 1.94: 353 1.94 - 2.61: 8938 2.61 - 3.27: 26654 3.27 - 3.94: 35141 3.94 - 4.60: 53231 Nonbonded interactions: 124317 Sorted by model distance: nonbonded pdb=" O PRO A 215 " pdb="HH22 ARG A 279 " model vdw 1.279 1.850 nonbonded pdb=" HG2 ARG A 75 " pdb=" HD3 PRO A 76 " model vdw 1.320 2.440 nonbonded pdb=" HB3 GLN A 27 " pdb="HH21 ARG A 30 " model vdw 1.321 2.270 nonbonded pdb=" HG3 GLU A 334 " pdb=" HD3 PRO A 336 " model vdw 1.355 2.440 nonbonded pdb="HH22 ARG A 291 " pdb=" HE1 HIS A 394 " model vdw 1.358 2.100 ... (remaining 124312 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. ============================== Scattering factors ============================= ----------X-ray scattering dictionary---------- Number of scattering types: 5 Type Number sf(0) Gaussians S 11 15.96 2 O 707 7.97 1 N 703 6.97 1 C 2641 5.97 2 H 3920 1.00 1 sf(0) = scattering factor at diffraction angle 0. =========================== Regularize ADP if needed ========================== ========================== Anomalous scatterer groups ========================= All atoms refined with f_prime=0 and f_double_prime=0. ====================== Modifying start model if requested ===================== Setting occupancies to: 1.000: selected atoms: 3920 of 7982 ******************************************************************************* Automatic adjustment: hydrogens.refine=riding ******************************************************************************* n_use = 7982 n_use_u_iso = 7982 n_use_u_aniso = 0 n_grad_site = 0 n_grad_u_iso = 0 n_grad_u_aniso = 0 n_grad_occupancy = 0 n_grad_fp = 0 n_grad_fdp = 0 total number of scatterers = 7982 ================== Extract refinement strategy and selections ================= Refinement flags and selection counts: individual_sites = True (7982 atoms) torsion_angles = False (0 atoms) rigid_body = False (0 atoms in 0 groups) individual_adp = True (iso = 7982 aniso = 0) group_adp = False (0 atoms in 0 groups) tls = False (0 atoms in 0 groups) occupancies = False (0 atoms) group_anomalous = False ============================ Non-default parameters =========================== A complete record of all parameters was written to the .eff file above. Below are only the non-defaults. #phil __ON__ refinement { electron_density_maps { map_coefficients { map_type = "2mFo-DFc" mtz_label_amplitudes = "2FOFCWT" mtz_label_phases = "PH2FOFCWT" fill_missing_f_obs = True } map_coefficients { map_type = "2mFo-DFc" mtz_label_amplitudes = "2FOFCWT_no_fill" mtz_label_phases = "PH2FOFCWT_no_fill" } map_coefficients { map_type = "mFo-DFc" mtz_label_amplitudes = "FOFCWT" mtz_label_phases = "PHFOFCWT" } map_coefficients { map_type = "anomalous" mtz_label_amplitudes = "ANOM" mtz_label_phases = "PANOM" } map { map_type = "2mFo-DFc" fill_missing_f_obs = True } map { map_type = "2mFo-DFc" } map { map_type = "mFo-DFc" } } refine { strategy = *individual_sites individual_sites_real_space rigid_body \ *individual_adp group_adp tls occupancies group_anomalous den } main { nqh_flips = False max_number_of_iterations = 100 nproc = 20 } modify_start_model { modify { selection = element H occupancies { set = 1 } } } hydrogens { refine = individual *riding Auto } pdb_interpretation { secondary_structure { protein { helix { serial_number = "1" helix_identifier = "1" selection = chain 'A' and resid 13 through 15 helix_type = alpha pi *3_10 unknown } helix { serial_number = "2" helix_identifier = "2" selection = chain 'A' and resid 27 through 38 } helix { serial_number = "3" helix_identifier = "3" selection = chain 'A' and resid 69 through 71 helix_type = alpha pi *3_10 unknown } helix { serial_number = "4" helix_identifier = "4" selection = chain 'A' and resid 75 through 80 } helix { serial_number = "5" helix_identifier = "5" selection = chain 'A' and resid 97 through 118 } helix { serial_number = "6" helix_identifier = "6" selection = chain 'A' and resid 123 through 133 } helix { serial_number = "7" helix_identifier = "7" selection = chain 'A' and resid 135 through 150 } helix { serial_number = "8" helix_identifier = "8" selection = chain 'A' and resid 156 through 158 helix_type = alpha pi *3_10 unknown } helix { serial_number = "9" helix_identifier = "9" selection = chain 'A' and resid 161 through 164 } helix { serial_number = "10" helix_identifier = "10" selection = chain 'A' and resid 168 through 177 } helix { serial_number = "11" helix_identifier = "11" selection = chain 'A' and resid 179 through 205 } helix { serial_number = "12" helix_identifier = "12" selection = chain 'A' and resid 223 through 227 } helix { serial_number = "13" helix_identifier = "13" selection = chain 'A' and resid 229 through 231 helix_type = alpha pi *3_10 unknown } helix { serial_number = "14" helix_identifier = "14" selection = chain 'A' and resid 265 through 270 } helix { serial_number = "15" helix_identifier = "15" selection = chain 'A' and resid 274 through 286 } helix { serial_number = "16" helix_identifier = "16" selection = chain 'A' and resid 295 through 298 helix_type = alpha pi *3_10 unknown } helix { serial_number = "17" helix_identifier = "17" selection = chain 'A' and resid 321 through 332 } helix { serial_number = "18" helix_identifier = "18" selection = chain 'A' and resid 347 through 356 } helix { serial_number = "19" helix_identifier = "19" selection = chain 'A' and resid 364 through 366 helix_type = alpha pi *3_10 unknown } helix { serial_number = "20" helix_identifier = "20" selection = chain 'A' and resid 378 through 380 helix_type = alpha pi *3_10 unknown } helix { serial_number = "21" helix_identifier = "21" selection = chain 'A' and resid 399 through 403 } helix { serial_number = "22" helix_identifier = "22" selection = chain 'A' and resid 408 through 421 } helix { serial_number = "23" helix_identifier = "23" selection = chain 'A' and resid 431 through 441 } helix { serial_number = "24" helix_identifier = "24" selection = chain 'A' and resid 451 through 454 } helix { serial_number = "25" helix_identifier = "25" selection = chain 'A' and resid 459 through 461 helix_type = alpha pi *3_10 unknown } helix { serial_number = "26" helix_identifier = "26" selection = chain 'A' and resid 484 through 496 } sheet { first_strand = chain 'A' and resid 446 through 450 sheet_id = " A" strand { selection = chain 'A' and resid 6 through 10 sense = *parallel antiparallel unknown bond_start_current = chain 'A' and resid 6 and name N bond_start_previous = chain 'A' and resid 447 and name O } strand { selection = chain 'A' and resid 43 through 45 sense = *parallel antiparallel unknown bond_start_current = chain 'A' and resid 43 and name N bond_start_previous = chain 'A' and resid 9 and name O } strand { selection = chain 'A' and resid 209 through 215 sense = *parallel antiparallel unknown bond_start_current = chain 'A' and resid 209 and name N bond_start_previous = chain 'A' and resid 44 and name O } strand { selection = chain 'A' and resid 289 through 293 sense = *parallel antiparallel unknown bond_start_current = chain 'A' and resid 289 and name N bond_start_previous = chain 'A' and resid 212 and name O } strand { selection = chain 'A' and resid 337 through 339 sense = *parallel antiparallel unknown bond_start_current = chain 'A' and resid 338 and name N bond_start_previous = chain 'A' and resid 290 and name O } } sheet { first_strand = chain 'A' and resid 300 through 303 sheet_id = " B" strand { selection = chain 'A' and resid 315 through 318 sense = parallel *antiparallel unknown bond_start_current = chain 'A' and resid 317 and name N bond_start_previous = chain 'A' and resid 301 and name O } } sheet { first_strand = chain 'A' and resid 360 through 363 sheet_id = " C" strand { selection = chain 'A' and resid 387 through 390 sense = *parallel antiparallel unknown bond_start_current = chain 'A' and resid 388 and name N bond_start_previous = chain 'A' and resid 360 and name O } } } enabled = True } nonbonded_weight = 100 } target_weights { optimize_xyz_weight = True optimize_adp_weight = True } } output { prefix = "1fp9_refine" serial = 1 target_output_format = None *pdb mmcif } #phil __OFF__ ================================== Hydrogens ================================== Total: count: 3920 occupancy sum: 3920.00 (% of total atoms 49.11) Rotatable: count: 903 occupancy sum: 903.00 (% of total atoms 11.31) ----------structure factors based statistics (before refinement)---------- ----------X-ray data---------- |--(resolution: 3.10 - 52.13 A, n_refl.=10132 (all), 10.01 % free)------------| | | | r_work= 0.4298 r_free= 0.5329 coordinate error (max.-lik. estimate): 0.35 A | | | | normalized target function (ml) (work): 5.559026 | | target function (ml) not normalized (work): 50687.195951 | | target function (ml) not normalized (free): 5601.206097 | |-----------------------------------------------------------------------------| |-----------------------------------------------------------------------------| | Bin Resolution Compl. No. Refl. R-factors Targets | |number range work test work test work test| | 1: 52.13 - 5.93 0.97 1355 151 0.4603 0.4878 6.2066 6.1554| | 2: 5.93 - 4.71 0.96 1323 147 0.3547 0.4794 5.7546 5.7028| | 3: 4.71 - 4.11 0.96 1288 144 0.4084 0.5271 5.5696 5.4451| | 4: 4.11 - 3.74 0.96 1288 143 0.4551 0.5640 5.5142 5.619| | 5: 3.74 - 3.47 0.96 1281 142 0.4754 0.5967 5.4726 5.2661| | 6: 3.47 - 3.26 0.97 1282 142 0.4418 0.5860 5.2538 5.3947| | 7: 3.26 - 3.10 0.97 1301 145 0.4200 0.5257 5.1055 5.0481| |-----------------------------------------------------------------------------| |-----------------------------------------------------------------------------| |R-free likelihood based estimates for figures of merit, absolute phase error,| |and distribution parameters alpha and beta (Acta Cryst. (1995). A51, 880-887)| | | | Bin Resolution No. Refl. FOM Phase Scale Alpha Beta | | # range work test error factor | | 1: 52.13 - 5.93 1355 151 0.16 77.62 0.60 0.03 59863.03| | 2: 5.93 - 4.71 1323 147 0.58 46.10 1.30 0.17 26644.58| | 3: 4.71 - 4.11 1288 144 0.83 24.31 1.53 0.26 9202.68| | 4: 4.11 - 3.74 1288 143 0.82 25.71 1.67 0.28 8178.16| | 5: 3.74 - 3.47 1281 142 0.80 27.27 1.71 0.29 7239.08| | 6: 3.47 - 3.26 1282 142 0.77 30.37 1.60 0.27 6252.41| | 7: 3.26 - 3.10 1301 145 0.73 33.46 1.45 0.23 4577.34| |alpha: min = 0.03 max = 0.29 mean = 0.22| |beta: min = 4577.34 max = 59863.03 mean = 17766.58| |figures of merit: min = 0.00 max = 1.00 mean = 0.67| |phase err.(work): min = 0.00 max = 90.00 mean = 38.16| |phase err.(test): min = 0.00 max = 89.99 mean = 38.55| |-----------------------------------------------------------------------------| ----------Initial model statistics (before refinement)---------- Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.046 4203 Z= 0.337 Angle : 1.249 9.998 5718 Z= 0.861 Chirality : 0.082 0.280 564 Planarity : 0.005 0.040 750 Dihedral : 14.688 74.118 1524 Min Nonbonded Distance : 2.289 Molprobity Statistics. All-atom Clashscore : 28.95 Ramachandran Plot: Outliers : 1.81 % Allowed : 12.25 % Favored : 85.94 % Rotamer: Outliers : 2.24 % Allowed : 2.24 % Favored : 95.51 % 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: -5.15 (0.30), residues: 498 helix: -4.21 (0.20), residues: 229 sheet: -1.62 (0.66), residues: 60 loop : -2.83 (0.38), residues: 209 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.006 0.001 ARG A 279 TYR 0.043 0.007 TYR A 141 PHE 0.032 0.005 PHE A 192 TRP 0.033 0.006 TRP A 159 HIS 0.002 0.000 HIS A 384 Individual atomic B min max mean iso aniso Overall: 12.70 77.70 33.40 2.74 4062 0 Protein: 12.70 77.70 33.40 2.74 4062 0 Chain A: 12.70 77.70 33.40 N/A 4062 0 Histogram: Values Number of atoms 12.70 - 19.20 560 19.20 - 25.70 851 25.70 - 32.20 919 32.20 - 38.70 595 38.70 - 45.20 378 45.20 - 51.70 244 51.70 - 58.20 183 58.20 - 64.70 129 64.70 - 71.20 85 71.20 - 77.70 118 |-Occupancies statistics------------------------------------------------------| | occupancies: max = 1.00 min = 1.00 number of occupancies < 0.1 = 0 | |-----------------------------------------------------------------------------| =================== Write initial parameters into .eff file =================== *********************** REFINEMENT MACRO_CYCLE 1 OF 3 ************************* |--(resolution: 3.10 - 52.13 A, n_refl.=10132 (all), 10.01 % free)------------| | | | r_work= 0.4298 r_free= 0.5329 coordinate error (max.-lik. estimate): 0.35 A | | | | normalized target function (ml) (work): 5.559026 | | target function (ml) not normalized (work): 50687.195951 | | target function (ml) not normalized (free): 5601.206097 | |-----------------------------------------------------------------------------| =========================== Bulk-solvent and scaling ========================== ============================= updating all scales ============================= start: r(all,work,free)=0.4407 0.4298 0.5329 n_refl.: 10132 re-set all scales: r(all,work,free)=0.4407 0.4298 0.5329 n_refl.: 10132 remove outliers: r(all,work,free)=0.4160 0.4027 0.5329 n_refl.: 10127 overall B=3.30 to atoms: r(all,work,free)=0.4309 0.4168 0.5513 n_refl.: 10127 bulk-solvent and scaling: r(all,work,free)=0.2399 0.2398 0.2416 n_refl.: 10127 remove outliers: r(all,work,free)=0.2397 0.2397 0.2400 n_refl.: 10125 ======================== Statistics in resolution bins ======================== Total model structure factor: F_model = k_total * (F_calc + k_mask * F_mask) k_total = k_isotropic * k_anisotropic Resolution Compl Nwork Nfree R_work kiso kani kmask 52.125-14.829 89.91 87 11 0.2450 1487.728 1446.815 1.051 1.275 0.404 14.763-11.705 97.09 91 9 0.1774 1550.646 1541.206 0.930 1.260 0.398 11.688-9.266 96.82 191 22 0.1307 1812.092 1781.536 0.927 1.236 0.395 9.245-7.329 97.77 355 40 0.1746 1351.329 1327.100 0.906 1.200 0.392 7.324-5.801 97.35 727 80 0.2319 959.455 924.737 0.843 1.146 0.389 5.796-4.593 96.03 1414 157 0.2081 1149.827 1113.607 0.910 1.055 0.352 4.593-3.637 95.90 2818 313 0.2376 1063.229 1016.166 1.121 0.912 0.295 3.636-3.100 96.41 3429 381 0.2953 729.417 669.032 1.293 0.739 0.111 Approximation of k_total with k_overall*exp(-b_overall*s**2/4) k_overall=1.0577 b_overall=4.8055 ============================ Set refinement target ============================ using refinement target ls_wunit_k1 ==================== Compute data/restraints target weights =================== |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 13.603384 wxc_scale = 0.500 wc = 1.000 | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxu = 0.301602 wxu_scale = 1.000 wu = 1.000 | |-----------------------------------------------------------------------------| ====================== XYZ individual (reciprocal space) ====================== Weight selection criteria: R_free deviation: 1.5000 R_free-R_work gap: 6.0000 Bond RMSD: 0.0150 Angle RMSD: 2.0000 R-FACTORS RMSD CLASH RAMA ROTA CBET WEIGHT work free delta bonds angl 0.1992 0.2730 0.0739 0.013 1.3 21.4 1.0 3.7 0 6.802 ====================== ADP refinement (reciprocal space) ====================== ----------Individual ADP refinement---------- R-FACTORS WEIGHT TARGETS work free delta data restr 19.92 27.30 7.39 1.968 36.488 0.302 0.047 20.31 27.17 6.87 3.989 37.021 0.009 0.049 19.86 27.20 7.34 3.925 36.970 0.038 0.047 18.80 27.35 8.54 4.382 36.883 0.151 0.042 18.22 27.51 9.29 5.311 36.950 0.302 0.040 17.88 27.59 9.71 6.291 37.162 0.452 0.038 17.71 27.66 9.95 6.935 37.296 0.603 0.037 17.59 27.71 10.12 7.468 37.415 0.754 0.037 17.50 27.77 10.26 7.906 37.517 0.905 0.036 17.44 27.83 10.39 8.289 37.612 1.056 0.036 17.39 27.86 10.47 8.583 37.684 1.206 0.036 17.30 27.92 10.62 9.153 37.848 1.357 0.036 17.25 27.97 10.72 9.535 37.954 1.508 0.035 max suggested for this run: 20.00 max allowed Rfree-Rwork gap: 6.0 range of equivalent Rfree: 1.5 Best ADP weight: 0.038 Accepted refinement result: 19.86 27.20 7.34 3.925 36.970 0.038 0.047 Individual atomic B min max mean iso aniso Overall: 18.64 82.30 37.65 2.51 4062 0 Protein: 18.64 82.30 37.65 2.51 4062 0 Chain A: 18.64 82.30 37.65 N/A 4062 0 Histogram: Values Number of atoms 18.64 - 25.01 350 25.01 - 31.38 1152 31.38 - 37.74 1035 37.74 - 44.11 575 44.11 - 50.47 375 50.47 - 56.84 195 56.84 - 63.20 153 63.20 - 69.57 89 69.57 - 75.93 78 75.93 - 82.30 60 =========================== Idealize ADP of riding H ========================== r_work=0.1986 r_free=0.2720 r_work=0.1985 r_free=0.2713 *********************** REFINEMENT MACRO_CYCLE 2 OF 3 ************************* |--(resolution: 3.10 - 52.13 A, n_refl.=10125 (all), 10.00 % free)------------| | | | r_work= 0.1985 r_free= 0.2713 coordinate error (max.-lik. estimate): 0.48 A | | | | normalized target function (ls_wunit_k1) (work): 0.047057 | | target function (ls_wunit_k1) not normalized (work): 428.786232 | | target function (ls_wunit_k1) not normalized (free): 93.353484 | |-----------------------------------------------------------------------------| =========================== Bulk-solvent and scaling ========================== ============================= updating all scales ============================= start: r(all,work,free)=0.2058 0.1985 0.2713 n_refl.: 10125 re-set all scales: r(all,work,free)=0.4216 0.4043 0.5550 n_refl.: 10125 remove outliers: r(all,work,free)=0.4216 0.4043 0.5550 n_refl.: 10125 overall B=0.00 to atoms: r(all,work,free)=0.4216 0.4043 0.5550 n_refl.: 10125 bulk-solvent and scaling: r(all,work,free)=0.2055 0.1981 0.2728 n_refl.: 10125 remove outliers: r(all,work,free)=0.2055 0.1981 0.2728 n_refl.: 10125 ======================== Statistics in resolution bins ======================== Total model structure factor: F_model = k_total * (F_calc + k_mask * F_mask) k_total = k_isotropic * k_anisotropic Resolution Compl Nwork Nfree R_work kiso kani kmask 52.125-14.829 89.91 87 11 0.2050 811.069 797.548 1.174 1.204 0.408 14.763-11.705 97.09 91 9 0.1310 845.370 848.086 1.033 1.192 0.407 11.688-9.266 96.82 191 22 0.1090 987.904 970.853 1.008 1.174 0.403 9.245-7.329 97.77 355 40 0.1461 736.708 726.743 0.974 1.146 0.400 7.324-5.801 97.35 727 80 0.1873 523.069 513.529 0.906 1.103 0.400 5.796-4.593 96.03 1414 157 0.1703 626.854 618.623 0.955 1.032 0.370 4.593-3.637 95.90 2818 313 0.1952 579.643 568.298 1.136 0.919 0.280 3.636-3.100 96.41 3429 381 0.2481 397.658 377.060 1.300 0.783 0.050 Approximation of k_total with k_overall*exp(-b_overall*s**2/4) k_overall=1.0610 b_overall=2.4036 ============================ Set refinement target ============================ using refinement target ml ============================ Update CDL restraints ============================ ==================== Compute data/restraints target weights =================== |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 2.971602 wxc_scale = 0.500 wc = 1.000 | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxu = 0.081802 wxu_scale = 1.000 wu = 1.000 | |-----------------------------------------------------------------------------| ====================== XYZ individual (reciprocal space) ====================== Weight selection criteria: R_free deviation: 1.5000 R_free-R_work gap: 6.0000 Bond RMSD: 0.0150 Angle RMSD: 2.0000 R-FACTORS RMSD CLASH RAMA ROTA CBET WEIGHT work free delta bonds angl 0.2254 0.2717 0.0463 0.003 0.6 4.6 0.0 2.7 0 0.125 0.2194 0.2713 0.0519 0.004 0.6 5.8 0.4 3.2 0 0.250 0.2111 0.2716 0.0605 0.006 0.8 9.4 0.8 4.2 0 0.500 0.2040 0.2788 0.0748 0.010 1.1 15.9 1.4 5.2 0 1.000 0.1949 0.2869 0.0921 0.018 1.6 26.7 3.2 7.5 0 2.000 0.1896 0.2897 0.1001 0.023 2.0 31.5 4.6 10.2 0 3.000 0.1870 0.2918 0.1048 0.029 2.3 38.1 5.8 13.5 0 4.000 0.1841 0.2928 0.1088 0.034 2.6 41.9 6.4 15.5 0 5.000 0.1819 0.2930 0.1111 0.038 2.9 44.9 6.2 17.0 0 6.000 0.1799 0.2955 0.1156 0.042 3.1 48.0 6.6 18.0 0 7.000 0.1778 0.2978 0.1200 0.046 3.4 53.0 7.6 18.0 1 8.000 0.1763 0.3010 0.1246 0.050 3.6 57.4 8.2 19.2 1 9.000 0.1745 0.3010 0.1266 0.053 3.7 59.4 8.8 19.2 2 10.000 0.1991 0.2836 0.0846 0.014 1.3 22.3 2.4 6.7 0 1.486 Best weight: R-FACTORS RMSD CLASH RAMA ROTA CBET WEIGHT work free delta bonds angl 0.2254 0.2717 0.0463 0.003 0.6 4.6 0.0 2.7 0 0.125 ====================== ADP refinement (reciprocal space) ====================== ----------Individual ADP refinement---------- R-FACTORS WEIGHT TARGETS work free delta data restr 22.54 27.17 4.63 5.769 41.297 0.082 6.382 23.61 27.76 4.15 6.377 43.332 0.002 6.416 22.39 27.17 4.78 6.006 41.212 0.010 6.374 21.84 27.21 5.37 6.370 41.046 0.041 6.359 20.81 27.44 6.63 8.389 41.263 0.082 6.330 20.52 27.53 7.01 9.250 41.639 0.123 6.323 20.33 27.61 7.29 10.054 41.941 0.164 6.319 20.19 27.70 7.51 10.769 42.188 0.205 6.317 20.09 27.76 7.68 11.392 42.391 0.245 6.315 20.01 27.84 7.83 11.937 42.564 0.286 6.314 19.95 27.90 7.95 12.394 42.706 0.327 6.313 19.91 27.96 8.05 12.808 42.833 0.368 6.313 19.87 28.00 8.13 13.152 42.939 0.409 6.313 max suggested for this run: 20.00 max allowed Rfree-Rwork gap: 6.0 range of equivalent Rfree: 1.5 Best ADP weight: 0.010 Accepted refinement result: 22.39 27.17 4.78 6.006 41.212 0.010 6.374 Individual atomic B min max mean iso aniso Overall: 15.64 83.27 37.49 2.54 4062 0 Protein: 15.64 83.27 37.49 2.54 4062 0 Chain A: 15.64 83.27 37.49 N/A 4062 0 Histogram: Values Number of atoms 15.64 - 22.40 99 22.40 - 29.17 940 29.17 - 35.93 1255 35.93 - 42.69 721 42.69 - 49.45 459 49.45 - 56.22 221 56.22 - 62.98 161 62.98 - 69.74 91 69.74 - 76.51 75 76.51 - 83.27 40 =========================== Idealize ADP of riding H ========================== r_work=0.2239 r_free=0.2717 r_work=0.2237 r_free=0.2717 *********************** REFINEMENT MACRO_CYCLE 3 OF 3 ************************* |--(resolution: 3.10 - 52.13 A, n_refl.=10125 (all), 10.00 % free)------------| | | | r_work= 0.2237 r_free= 0.2717 coordinate error (max.-lik. estimate): 0.50 A | | | | normalized target function (ml) (work): 6.373857 | | target function (ml) not normalized (work): 58078.580673 | | target function (ml) not normalized (free): 6576.740314 | |-----------------------------------------------------------------------------| =========================== Bulk-solvent and scaling ========================== ============================= updating all scales ============================= start: r(all,work,free)=0.2285 0.2237 0.2717 n_refl.: 10125 re-set all scales: r(all,work,free)=0.4283 0.4121 0.5569 n_refl.: 10125 remove outliers: r(all,work,free)=0.4283 0.4121 0.5569 n_refl.: 10125 overall B=0.00 to atoms: r(all,work,free)=0.4283 0.4121 0.5569 n_refl.: 10125 bulk-solvent and scaling: r(all,work,free)=0.2292 0.2244 0.2724 n_refl.: 10125 remove outliers: r(all,work,free)=0.2292 0.2244 0.2724 n_refl.: 10125 ======================== Statistics in resolution bins ======================== Total model structure factor: F_model = k_total * (F_calc + k_mask * F_mask) k_total = k_isotropic * k_anisotropic Resolution Compl Nwork Nfree R_work kiso kani kmask 52.125-14.829 89.91 87 11 0.2110 811.069 804.963 1.301 1.106 0.408 14.763-11.705 97.09 91 9 0.1385 845.370 851.052 1.132 1.100 0.407 11.688-9.266 96.82 191 22 0.1172 987.904 975.160 1.094 1.089 0.407 9.245-7.329 97.77 355 40 0.1536 736.708 731.968 1.050 1.072 0.403 7.324-5.801 97.35 727 80 0.2123 523.069 512.557 0.960 1.047 0.400 5.796-4.593 96.03 1414 157 0.1890 626.854 613.645 0.958 1.004 0.390 4.593-3.637 95.90 2818 313 0.2260 579.643 557.430 1.069 0.936 0.310 3.636-3.100 96.41 3429 381 0.2826 397.658 367.072 1.131 0.855 0.080 Approximation of k_total with k_overall*exp(-b_overall*s**2/4) k_overall=1.0577 b_overall=4.4071 ============================ Set refinement target ============================ using refinement target ml ============================ Update CDL restraints ============================ ============================= Setup riding H model ============================ Before: r_work=0.2244 r_free=0.2724 After: r_work=0.2245 r_free=0.2724 ==================== Compute data/restraints target weights =================== |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 3.162007 wxc_scale = 0.500 wc = 1.000 | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxu = 0.086166 wxu_scale = 1.000 wu = 1.000 | |-----------------------------------------------------------------------------| ====================== XYZ individual (reciprocal space) ====================== Weight selection criteria: R_free deviation: 1.5000 R_free-R_work gap: 6.0000 Bond RMSD: 0.0150 Angle RMSD: 2.0000 R-FACTORS RMSD CLASH RAMA ROTA CBET WEIGHT work free delta bonds angl 0.2248 0.2728 0.0479 0.003 0.6 5.0 0.0 2.7 0 0.125 0.2238 0.2728 0.0490 0.004 0.6 5.3 0.0 2.7 0 0.250 0.2197 0.2729 0.0532 0.007 0.7 5.5 0.0 2.7 0 0.500 0.2084 0.2748 0.0664 0.015 1.3 10.5 0.0 2.7 0 1.000 0.2000 0.2791 0.0791 0.017 1.9 23.7 0.8 4.2 0 2.000 0.1980 0.2787 0.0807 0.021 1.8 24.2 1.0 4.7 0 3.000 0.1951 0.2810 0.0859 0.023 2.2 28.2 1.4 6.2 0 4.000 0.1919 0.2845 0.0926 0.029 2.7 35.6 3.0 8.0 0 5.000 0.1869 0.2908 0.1039 0.038 3.6 48.7 3.8 9.5 0 6.000 0.1852 0.2920 0.1068 0.040 3.8 50.5 4.0 10.0 1 7.000 0.1844 0.2922 0.1078 0.042 3.9 50.4 4.8 12.0 1 8.000 0.1830 0.2938 0.1108 0.046 4.2 54.8 5.0 13.2 2 9.000 0.1820 0.2965 0.1145 0.049 4.5 61.0 6.2 16.5 3 10.000 0.2022 0.2773 0.0751 0.016 1.7 20.3 0.6 4.0 0 1.581 Best weight: R-FACTORS RMSD CLASH RAMA ROTA CBET WEIGHT work free delta bonds angl 0.2248 0.2728 0.0479 0.003 0.6 5.0 0.0 2.7 0 0.125 ====================== ADP refinement (reciprocal space) ====================== ----------Individual ADP refinement---------- R-FACTORS WEIGHT TARGETS work free delta data restr 22.48 27.28 4.79 5.742 41.129 0.086 6.380 23.87 27.94 4.07 6.245 42.747 0.003 6.424 22.48 27.30 4.82 5.787 41.101 0.011 6.380 21.99 27.30 5.31 6.197 40.924 0.043 6.366 20.84 27.65 6.81 8.508 40.935 0.086 6.335 20.57 27.77 7.20 9.370 41.332 0.129 6.330 20.39 27.93 7.53 10.170 41.651 0.172 6.327 20.27 28.05 7.78 10.882 41.909 0.215 6.325 20.18 28.17 7.99 11.546 42.133 0.258 6.324 20.12 28.26 8.14 12.048 42.300 0.302 6.323 20.08 28.33 8.26 12.516 42.451 0.345 6.323 20.04 28.40 8.36 12.920 42.579 0.388 6.323 20.01 28.46 8.45 13.277 42.691 0.431 6.323 max suggested for this run: 20.00 max allowed Rfree-Rwork gap: 6.0 range of equivalent Rfree: 1.5 Best ADP weight: 0.011 Accepted refinement result: 22.48 27.30 4.82 5.787 41.101 0.011 6.380 Individual atomic B min max mean iso aniso Overall: 15.99 83.46 37.43 2.52 4062 0 Protein: 15.99 83.46 37.43 2.52 4062 0 Chain A: 15.99 83.46 37.43 N/A 4062 0 Histogram: Values Number of atoms 15.99 - 22.73 121 22.73 - 29.48 989 29.48 - 36.23 1237 36.23 - 42.98 692 42.98 - 49.73 457 49.73 - 56.47 215 56.47 - 63.22 154 63.22 - 69.97 92 69.97 - 76.72 70 76.72 - 83.46 35 =========================== Idealize ADP of riding H ========================== r_work=0.2248 r_free=0.2730 r_work=0.2248 r_free=0.2730 ----------X-ray data---------- |--(resolution: 3.10 - 52.13 A, n_refl.=10125 (all), 10.00 % free)------------| | | | r_work= 0.2248 r_free= 0.2730 coordinate error (max.-lik. estimate): 0.49 A | | | | normalized target function (ml) (work): 6.380342 | | target function (ml) not normalized (work): 58137.673780 | | target function (ml) not normalized (free): 6582.498793 | |-----------------------------------------------------------------------------| |-----------------------------------------------------------------------------| | Bin Resolution Compl. No. Refl. R-factors Targets | |number range work test work test work test| | 1: 52.13 - 5.93 0.97 1352 151 0.1700 0.2023 6.3733 6.5428| | 2: 5.93 - 4.71 0.96 1323 147 0.1929 0.2451 6.4134 6.5236| | 3: 4.71 - 4.11 0.96 1288 144 0.2078 0.2309 6.5081 6.5525| | 4: 4.11 - 3.74 0.96 1286 142 0.2315 0.3088 6.4731 6.714| | 5: 3.74 - 3.47 0.96 1281 142 0.2686 0.3093 6.4536 6.4771| | 6: 3.47 - 3.26 0.97 1281 142 0.2696 0.3363 6.2972 6.4782| | 7: 3.26 - 3.10 0.97 1301 145 0.3033 0.3696 6.1455 6.1998| |-----------------------------------------------------------------------------| |-----------------------------------------------------------------------------| |R-free likelihood based estimates for figures of merit, absolute phase error,| |and distribution parameters alpha and beta (Acta Cryst. (1995). A51, 880-887)| | | | Bin Resolution No. Refl. FOM Phase Scale Alpha Beta | | # range work test error factor | | 1: 52.13 - 5.93 1352 151 0.88 18.46 0.99 0.97 58555.12| | 2: 5.93 - 4.71 1323 147 0.86 21.30 1.00 0.93 66792.63| | 3: 4.71 - 4.11 1288 144 0.85 21.96 1.01 0.91 72343.23| | 4: 4.11 - 3.74 1286 142 0.82 25.51 0.99 0.92 77837.32| | 5: 3.74 - 3.47 1281 142 0.78 29.00 1.06 0.93 74966.43| | 6: 3.47 - 3.26 1281 142 0.73 33.44 1.03 0.87 66123.59| | 7: 3.26 - 3.10 1301 145 0.69 37.02 0.91 0.73 48678.41| |alpha: min = 0.73 max = 0.97 mean = 0.89| |beta: min = 48678.41 max = 77837.32 mean = 66382.46| |figures of merit: min = 0.00 max = 1.00 mean = 0.80| |phase err.(work): min = 0.00 max = 89.70 mean = 26.60| |phase err.(test): min = 0.00 max = 88.40 mean = 27.36| |-----------------------------------------------------------------------------| ==================================== Final ==================================== ============================= updating all scales ============================= start: r(all,work,free)=0.2296 0.2248 0.2730 n_refl.: 10125 re-set all scales: r(all,work,free)=0.4283 0.4121 0.5569 n_refl.: 10125 remove outliers: r(all,work,free)=0.4283 0.4121 0.5569 n_refl.: 10125 overall B=0.00 to atoms: r(all,work,free)=0.4283 0.4121 0.5569 n_refl.: 10125 bulk-solvent and scaling: r(all,work,free)=0.2303 0.2254 0.2739 n_refl.: 10125 remove outliers: r(all,work,free)=0.2303 0.2254 0.2739 n_refl.: 10125 ======================== Statistics in resolution bins ======================== Total model structure factor: F_model = k_total * (F_calc + k_mask * F_mask) k_total = k_isotropic * k_anisotropic Resolution Compl Nwork Nfree R_work kiso kani kmask 52.125-14.829 89.91 87 11 0.2085 811.069 806.507 1.350 1.072 0.408 14.763-11.705 97.09 91 9 0.1401 845.370 852.103 1.170 1.068 0.407 11.688-9.266 96.82 191 22 0.1169 987.904 975.998 1.128 1.060 0.407 9.245-7.329 97.77 355 40 0.1556 736.708 732.514 1.079 1.047 0.403 7.324-5.801 97.35 727 80 0.2132 523.069 512.702 0.981 1.028 0.400 5.796-4.593 96.03 1414 157 0.1894 626.854 613.469 0.967 0.995 0.390 4.593-3.637 95.90 2818 313 0.2269 579.643 557.208 1.063 0.943 0.304 3.636-3.100 96.41 3429 381 0.2844 397.658 366.656 1.097 0.882 0.080 Approximation of k_total with k_overall*exp(-b_overall*s**2/4) k_overall=1.0577 b_overall=4.2874 ================= overall refinement statistics: step by step ================= ****************** REFINEMENT STATISTICS STEP BY STEP ****************** leading digit, like 1_, means number of macro-cycle 0 : statistics at the very beginning when nothing is done yet 1_bss: bulk solvent correction and/or (anisotropic) scaling 1_xyz: refinement of coordinates 1_adp: refinement of ADPs (Atomic Displacement Parameters) ------------------------------------------------------------------------ stage r-work r-free bonds angles b_min b_max b_ave n_water shift 0 : 0.4298 0.5329 0.005 1.249 12.7 77.7 33.4 0 0.000 1_bss: 0.2397 0.2400 0.005 1.249 16.0 81.0 36.7 0 0.000 1_settarget: 0.2397 0.2400 0.005 1.249 16.0 81.0 36.7 0 0.000 1_weight: 0.2397 0.2400 0.005 1.249 16.0 81.0 36.7 0 0.000 1_xyzrec: 0.1992 0.2730 0.013 1.275 16.0 81.0 36.7 0 0.272 1_adp: 0.1986 0.2720 0.013 1.275 18.6 82.3 37.7 0 0.272 1_regHadp: 0.1985 0.2713 0.013 1.275 18.6 82.3 37.7 0 0.272 2_bss: 0.1981 0.2728 0.013 1.275 18.6 82.3 37.7 0 0.272 2_settarget: 0.1981 0.2728 0.013 1.275 18.6 82.3 37.7 0 0.272 2_updatecdl: 0.1981 0.2728 0.013 1.316 18.6 82.3 37.7 0 0.272 2_weight: 0.1981 0.2728 0.013 1.316 18.6 82.3 37.7 0 0.272 2_xyzrec: 0.2254 0.2717 0.003 0.550 18.6 82.3 37.7 0 0.260 2_adp: 0.2239 0.2717 0.003 0.550 15.6 83.3 37.5 0 0.260 2_regHadp: 0.2237 0.2717 0.003 0.550 15.6 83.3 37.5 0 0.260 3_bss: 0.2244 0.2724 0.003 0.550 15.6 83.3 37.5 0 0.260 3_settarget: 0.2244 0.2724 0.003 0.550 15.6 83.3 37.5 0 0.260 3_updatecdl: 0.2244 0.2724 0.003 0.619 15.6 83.3 37.5 0 0.260 3_setrh: 0.2245 0.2724 0.003 0.619 15.6 83.3 37.5 0 0.260 3_weight: 0.2245 0.2724 0.003 0.619 15.6 83.3 37.5 0 0.260 3_xyzrec: 0.2248 0.2728 0.003 0.551 15.6 83.3 37.5 0 0.260 3_adp: 0.2248 0.2730 0.003 0.551 16.0 83.5 37.4 0 0.260 3_regHadp: 0.2248 0.2730 0.003 0.551 16.0 83.5 37.4 0 0.260 end: 0.2254 0.2739 0.003 0.551 16.0 83.5 37.4 0 0.260 ------------------------------------------------------------------------ ============================== Exporting results ============================== Writing refined structure to CIF file: /net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/phenix_V3/1fp9_refine_001.cif Writing default parameters for subsequent refinement: /net/anaconda/raid1/afonine/work/ANI/result/xray/pairs/phenix_V3/1fp9_refine_002.def Timing summary: Write initial files (write_before_run_outputs): 11.8300 Refinement macro-cycles (run) : 300.5000 Write final files (write_after_run_outputs) : 22.9500 Total : 335.2800 Total CPU time: 6.36 minutes =========================== phenix.refine: finished =========================== # Date 2024-05-20 Time 17:35:11 PDT -0700 (1716251711.58 s) Start R-work = 0.2397, R-free = 0.2400 Final R-work = 0.2254, R-free = 0.2739 =============================================================================== Job complete usr+sys time: 398.09 seconds wall clock time: 17 minutes 55.31 seconds (1075.31 seconds total)