============================== Collecting inputs ============================== ----------Processing X-ray data---------- ----------Processing PDB file(s)---------- Symmetric amino acids flipped Residue "B PHE 225": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Time to flip residues: 0.00s Monomer Library directory: "/dev/shm/nigel/phenix-1.20.1-4487/modules/chem_data/mon_lib" Total number of atoms: 9333 Number of models: 1 Model: "" Number of chains: 3 Chain: "A" Number of atoms: 3111 Number of conformers: 1 Conformer: "" Number of residues, atoms: 215, 3111 Classifications: {'peptide': 215} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 1, 'TRANS': 213} Chain: "B" Number of atoms: 3111 Number of conformers: 1 Conformer: "" Number of residues, atoms: 215, 3111 Classifications: {'peptide': 215} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 1, 'TRANS': 213} Chain: "C" Number of atoms: 3111 Number of conformers: 1 Conformer: "" Number of residues, atoms: 215, 3111 Classifications: {'peptide': 215} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 1, 'TRANS': 213} Time building chain proxies: 1.64, per 1000 atoms: 0.18 Number of scatterers: 9333 At special positions: 0 Unit cell: (110.36, 115.96, 164.12, 90, 90, 90) Space group: I 2 2 2 (No. 23) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 9 16.00 O 966 8.00 N 909 7.00 C 2865 6.00 H 4584 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Custom bonds: Warning: Ignoring bond with distance_ideal = None: atom_selection_1 = None atom_selection_2 = None Total number of added/changed bonds: 0 Custom angles: Warning: Ignoring angle with angle_ideal = None: atom_selection_1 = None atom_selection_2 = None atom_selection_3 = None Total number of new custom angles: 0 Total number of changed angles: 0 Custom dihedrals: Warning: Ignoring dihedral with angle_ideal = None: atom_selection_1 = None atom_selection_2 = None atom_selection_3 = None atom_selection_4 = None Total number of custom dihedrals: 0 Custom planarities: Warning: Ignoring planarity with with sigma <= 0: atom_selection = None None Total number of custom planarities: 0 Custom parallelities: Warning: Ignoring parallelity with empty atom selection. Total number of custom parallelities: 0 Custom bonds: Warning: Ignoring bond with distance_ideal = None: atom_selection_1 = None atom_selection_2 = None Total number of added/changed bonds: 0 Custom angles: Warning: Ignoring angle with angle_ideal = None: atom_selection_1 = None atom_selection_2 = None atom_selection_3 = None Total number of new custom angles: 0 Total number of changed angles: 0 Custom dihedrals: Warning: Ignoring dihedral with angle_ideal = None: atom_selection_1 = None atom_selection_2 = None atom_selection_3 = None atom_selection_4 = None Total number of custom dihedrals: 0 Custom planarities: Warning: Ignoring planarity with with sigma <= 0: atom_selection = None None Total number of custom planarities: 0 Custom parallelities: Warning: Ignoring parallelity with empty atom selection. Total number of custom parallelities: 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: 2.33 Conformation dependent library (CDL) restraints added in 225.5 milliseconds Adding C-beta torsion restraints... Number of C-beta restraints generated: 1200 Time building geometry restraints manager: 2.55 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: 4581 1.15 - 1.32: 1124 1.32 - 1.49: 1726 1.49 - 1.66: 1935 1.66 - 1.84: 15 Bond restraints: 9381 Sorted by residual: bond pdb=" N GLU A 44 " pdb=" H GLU A 44 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" N ARG A 186 " pdb=" H ARG A 186 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" N ALA A 210 " pdb=" H ALA A 210 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" N LEU B 177 " pdb=" H LEU B 177 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 bond pdb=" N ALA A 228 " pdb=" H ALA A 228 " ideal model delta sigma weight residual 0.860 1.021 -0.161 2.00e-02 2.50e+03 6.45e+01 ... (remaining 9376 not shown) Histogram of bond angle deviations from ideal: 98.29 - 105.44: 223 105.44 - 112.60: 10348 112.60 - 119.75: 2632 119.75 - 126.90: 3592 126.90 - 134.06: 50 Bond angle restraints: 16845 Sorted by residual: angle pdb=" N LEU A 233 " pdb=" CA LEU A 233 " pdb=" C LEU A 233 " ideal model delta sigma weight residual 110.80 119.87 -9.07 2.13e+00 2.20e-01 1.81e+01 angle pdb=" N THR B 209 " pdb=" CA THR B 209 " pdb=" C THR B 209 " ideal model delta sigma weight residual 113.88 119.07 -5.19 1.23e+00 6.61e-01 1.78e+01 angle pdb=" N SER A 211 " pdb=" CA SER A 211 " pdb=" CB SER A 211 " ideal model delta sigma weight residual 110.06 116.02 -5.96 1.51e+00 4.39e-01 1.56e+01 angle pdb=" CG ARG A 191 " pdb=" CD ARG A 191 " pdb=" NE ARG A 191 " ideal model delta sigma weight residual 112.00 103.57 8.43 2.20e+00 2.07e-01 1.47e+01 angle pdb=" CA ALA C 210 " pdb=" C ALA C 210 " pdb=" O ALA C 210 " ideal model delta sigma weight residual 121.47 117.12 4.35 1.15e+00 7.56e-01 1.43e+01 ... (remaining 16840 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 15.79: 3372 15.79 - 31.57: 240 31.57 - 47.36: 85 47.36 - 63.15: 18 63.15 - 78.93: 8 Dihedral angle restraints: 3723 sinusoidal: 1713 harmonic: 2010 Sorted by residual: dihedral pdb=" C ASP A 218 " pdb=" N ASP A 218 " pdb=" CA ASP A 218 " pdb=" CB ASP A 218 " ideal model delta harmonic sigma weight residual -122.60 -113.38 -9.22 0 2.50e+00 1.60e-01 1.36e+01 dihedral pdb=" CA LEU A 233 " pdb=" C LEU A 233 " pdb=" N ALA A 234 " pdb=" CA ALA A 234 " ideal model delta harmonic sigma weight residual -180.00 -162.17 -17.83 0 5.00e+00 4.00e-02 1.27e+01 dihedral pdb=" C VAL C 205 " pdb=" N VAL C 205 " pdb=" CA VAL C 205 " pdb=" CB VAL C 205 " ideal model delta harmonic sigma weight residual -122.00 -130.88 8.88 0 2.50e+00 1.60e-01 1.26e+01 ... (remaining 3720 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.093: 559 0.093 - 0.186: 157 0.186 - 0.278: 32 0.278 - 0.371: 16 0.371 - 0.464: 1 Chirality restraints: 765 Sorted by residual: chirality pdb=" CA ARG A 208 " pdb=" N ARG A 208 " pdb=" C ARG A 208 " pdb=" CB ARG A 208 " both_signs ideal model delta sigma weight residual False 2.51 2.05 0.46 2.00e-01 2.50e+01 5.38e+00 chirality pdb=" CA MET B 171 " pdb=" N MET B 171 " pdb=" C MET B 171 " pdb=" CB MET B 171 " both_signs ideal model delta sigma weight residual False 2.51 2.17 0.34 2.00e-01 2.50e+01 2.90e+00 chirality pdb=" CA GLU A 109 " pdb=" N GLU A 109 " pdb=" C GLU A 109 " pdb=" CB GLU A 109 " both_signs ideal model delta sigma weight residual False 2.51 2.18 0.33 2.00e-01 2.50e+01 2.67e+00 ... (remaining 762 not shown) Planarity restraints: 1539 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CD ARG A 208 " 0.046 9.50e-02 1.11e+02 3.78e-02 2.72e+01 pdb=" NE ARG A 208 " -0.032 2.00e-02 2.50e+03 pdb=" CZ ARG A 208 " -0.052 2.00e-02 2.50e+03 pdb=" NH1 ARG A 208 " -0.001 2.00e-02 2.50e+03 pdb=" NH2 ARG A 208 " 0.001 2.00e-02 2.50e+03 pdb="HH11 ARG A 208 " 0.058 2.00e-02 2.50e+03 pdb="HH12 ARG A 208 " -0.020 2.00e-02 2.50e+03 pdb="HH21 ARG A 208 " 0.056 2.00e-02 2.50e+03 pdb="HH22 ARG A 208 " -0.013 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR A 59 " -0.022 2.00e-02 2.50e+03 2.61e-02 2.05e+01 pdb=" CG TYR A 59 " -0.013 2.00e-02 2.50e+03 pdb=" CD1 TYR A 59 " -0.006 2.00e-02 2.50e+03 pdb=" CD2 TYR A 59 " 0.016 2.00e-02 2.50e+03 pdb=" CE1 TYR A 59 " 0.006 2.00e-02 2.50e+03 pdb=" CE2 TYR A 59 " -0.016 2.00e-02 2.50e+03 pdb=" CZ TYR A 59 " -0.006 2.00e-02 2.50e+03 pdb=" OH TYR A 59 " 0.012 2.00e-02 2.50e+03 pdb=" HD1 TYR A 59 " -0.010 2.00e-02 2.50e+03 pdb=" HD2 TYR A 59 " 0.062 2.00e-02 2.50e+03 pdb=" HE1 TYR A 59 " 0.024 2.00e-02 2.50e+03 pdb=" HE2 TYR A 59 " -0.048 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR C 59 " -0.005 2.00e-02 2.50e+03 2.61e-02 2.05e+01 pdb=" CG TYR C 59 " -0.016 2.00e-02 2.50e+03 pdb=" CD1 TYR C 59 " -0.011 2.00e-02 2.50e+03 pdb=" CD2 TYR C 59 " 0.013 2.00e-02 2.50e+03 pdb=" CE1 TYR C 59 " 0.008 2.00e-02 2.50e+03 pdb=" CE2 TYR C 59 " -0.016 2.00e-02 2.50e+03 pdb=" CZ TYR C 59 " -0.002 2.00e-02 2.50e+03 pdb=" OH TYR C 59 " 0.014 2.00e-02 2.50e+03 pdb=" HD1 TYR C 59 " -0.023 2.00e-02 2.50e+03 pdb=" HD2 TYR C 59 " 0.057 2.00e-02 2.50e+03 pdb=" HE1 TYR C 59 " 0.030 2.00e-02 2.50e+03 pdb=" HE2 TYR C 59 " -0.049 2.00e-02 2.50e+03 ... (remaining 1536 not shown) Histogram of nonbonded interaction distances: 1.27 - 1.94: 322 1.94 - 2.60: 9707 2.60 - 3.27: 30939 3.27 - 3.93: 40839 3.93 - 4.60: 59795 Nonbonded interactions: 141602 Sorted by model distance: nonbonded pdb=" HE1 HIS A 198 " pdb=" HB3 ALA A 234 " model vdw sym.op. 1.270 2.270 x,-y+2,-z+2 nonbonded pdb=" HD3 ARG A 208 " pdb=" O ARG A 208 " model vdw sym.op. 1.352 2.620 -x+2,y,-z+2 nonbonded pdb=" O ARG A 208 " pdb=" HD3 ARG A 208 " model vdw sym.op. 1.352 2.620 -x+2,y,-z+2 nonbonded pdb=" HB3 ALA B 234 " pdb=" HE ARG C 216 " model vdw sym.op. 1.398 2.270 x,-y+2,-z+2 nonbonded pdb=" OXT ALA B 234 " pdb="HH21 ARG C 216 " model vdw sym.op. 1.425 1.850 x,-y+2,-z+2 ... (remaining 141597 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 9 15.91 1 O 966 7.97 1 N 909 6.97 1 C 2865 5.97 1 H 4584 1.00 1 sf(0) = scattering factor at diffraction angle 0. Number of scatterers: 9333 At special positions: 0 Unit cell: (110.36, 115.96, 164.12, 90, 90, 90) Space group: I 2 2 2 (No. 23) ----------F(model) initialization---------- Twinning will be detected automatically. start: r(all,work,free)=0.4182 0.4213 0.3451 n_refl.: 12716 re-set all scales: r(all,work,free)=0.4182 0.4213 0.3451 n_refl.: 12716 remove outliers: r(all,work,free)=0.4086 0.4113 0.3451 n_refl.: 12706 bulk-solvent and scaling: r(all,work,free)=0.2233 0.2222 0.2445 n_refl.: 12706 add H (0.00, 0.00): r(all,work,free)=0.2283 0.2270 0.2526 n_refl.: 12706 remove outliers: r(all,work,free)=0.2282 0.2269 0.2526 n_refl.: 12705 |--(resolution: 3.55 - 32.24 A, n_refl.=12705 (all), 4.90 % free)------------| | | | r_work= 0.2269 r_free= 0.2526 coordinate error (max.-lik. estimate): 0.25 A | | | | normalized target function (ml) (work): 6.896248 | | target function (ml) not normalized (work): 83320.463398 | | target function (ml) not normalized (free): 4336.834498 | |-----------------------------------------------------------------------------| End of input processing ======================= ROSETTA/PHENIX X-ray refinement ======================= Output directory: /net/cci-filer3/home/nigel/work/ANI/result/xray/pairs/rosetta_12 ----------Setting up input files for Rosetta---------- |-starting model--------------------------------------------------------------| | target_work(ml) = 6.89625 r_work = 0.2269 r_free = 0.2526 | |-----------------------------------------------------------------------------| Validation statistics (starting model): Ramachandran outliers = 0.00 % favored = 99.69 % Rotamer outliers = 1.49 % C-beta deviations = 0 Clashscore = 1.82 RMS(bonds) = 0.0167 RMS(angles) = 1.56 MolProbity score = 1.07 Rosetta command-line arguments: -parser:protocol /dev/shm/bkpoon/software/rosetta_src_2021.16.61629_bundle/main/source/src/apps/public/crystal_refinement/low_resolution_refine.xml -s 4yei_rosetta_in.pdb -mtzfile 4yei_rosetta_data.mtz -run:preserve_header -crystal_refine -parser:script_vars symmdef=4yei_rosetta.symm -parser:script_vars bfactstrat=individual -parser:script_vars map_type=Auto -nstruct 1 Generating 5 models on 1 processors... 1: r_work = 0.2908 r_free = 0.3438 energy = -4887.74 rmsd = 0.897 *** 2: r_work = 0.2944 r_free = 0.3513 energy = -4956.01 rmsd = 0.908 3: r_work = 0.2931 r_free = 0.3503 energy = -4933.73 rmsd = 0.903 4: r_work = 0.2931 r_free = 0.3533 energy = -4894.63 rmsd = 0.886 5: r_work = 0.2943 r_free = 0.3480 energy = -4914.38 rmsd = 0.910 |-after ROSETTA---------------------------------------------------------------| | target_work(ml) = 7.16633 r_work = 0.2908 r_free = 0.3438 | |-----------------------------------------------------------------------------| Validation statistics after ROSETTA: Ramachandran outliers = 0.00 % favored = 92.96 % Rotamer outliers = 0.00 % C-beta deviations = 1 Clashscore = 1.29 RMS(bonds) = 0.0248 RMS(angles) = 1.88 MolProbity score = 1.30 ----------Running phenix.refine (with null strategy)---------- start r_work=0.4690 r_free=0.4523 1_bss r_work=0.2202 r_free=0.2780 end r_work=0.2202 r_free=0.2780 Ramachandran outliers = 0.00 % favored = 92.96 % Rotamer outliers = 0.00 % C-beta deviations = 1 Clashscore = 4.82 RMS(bonds) = 0.0906 RMS(angles) = 2.07 MolProbity score = 1.70 Resolution = 3.55 R-work = 0.2202 R-free = 0.2780 ----------Final results---------- Refined model: 4yei_rosetta_phenix_001.pdb Final maps: 4yei_rosetta_phenix_001.mtz Elapsed time: 18059.0s Citation: DiMaio F., Echols N., Headd J.J., Terwilliger T.C., Adams P.D., and Baker D. (2013). Improved low-resolution crystallographic refinement with Phenix and Rosetta. Nat Methods 10, 1102-4.