Starting phenix.real_space_refine on Fri Jun 13 14:47:29 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/8v5a_42981/06_2025/8v5a_42981.cif Found real_map, /net/cci-nas-00/data/ceres_data/8v5a_42981/06_2025/8v5a_42981.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.0 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/8v5a_42981/06_2025/8v5a_42981.map" default_real_map = "/net/cci-nas-00/data/ceres_data/8v5a_42981/06_2025/8v5a_42981.map" model { file = "/net/cci-nas-00/data/ceres_data/8v5a_42981/06_2025/8v5a_42981.cif" } default_model = "/net/cci-nas-00/data/ceres_data/8v5a_42981/06_2025/8v5a_42981.cif" } resolution = 3.0 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } 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 Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- 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.001 sd= 0.019 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Anisotropic ADP refinement not supported. Converting 2628 atoms to isotropic. Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 54 5.16 5 C 8226 2.51 5 N 2232 2.21 5 O 2544 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 93 residue(s): 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/chem_data/mon_lib" Total number of atoms: 13056 Number of models: 1 Model: "" Number of chains: 2 Chain: "A" Number of atoms: 3476 Number of conformers: 1 Conformer: "" Number of residues, atoms: 450, 3476 Classifications: {'peptide': 450} Link IDs: {'PTRANS': 19, 'TRANS': 430} Chain breaks: 2 Chain: "U" Number of atoms: 876 Number of conformers: 1 Conformer: "" Number of residues, atoms: 117, 876 Classifications: {'peptide': 117} Link IDs: {'PTRANS': 3, 'TRANS': 113} Restraints were copied for chains: W, Y, C, B Time building chain proxies: 8.66, per 1000 atoms: 0.66 Number of scatterers: 13056 At special positions: 0 Unit cell: (125.55, 132.84, 132.84, 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 2544 8.00 N 2232 7.00 C 8226 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=16, symmetry=0 Simple disulfide: pdb=" SG CYS A 63 " - pdb=" SG CYS A 192 " distance=2.03 Simple disulfide: pdb=" SG CYS A 331 " - pdb=" SG CYS A 340 " distance=2.02 Simple disulfide: pdb=" SG CYS A 355 " - pdb=" SG CYS A 363 " distance=2.03 Simple disulfide: pdb=" SG CYS A 387 " - pdb=" SG CYS A 392 " distance=2.04 Simple disulfide: pdb=" SG CYS A 394 " - pdb=" SG CYS A 417 " distance=2.03 Simple disulfide: pdb=" SG CYS Y 22 " - pdb=" SG CYS Y 95 " distance=2.03 Simple disulfide: pdb=" SG CYS C 63 " - pdb=" SG CYS C 192 " distance=2.03 Simple disulfide: pdb=" SG CYS B 63 " - pdb=" SG CYS B 192 " distance=2.03 Simple disulfide: pdb=" SG CYS C 331 " - pdb=" SG CYS C 340 " distance=2.02 Simple disulfide: pdb=" SG CYS B 331 " - pdb=" SG CYS B 340 " distance=2.02 Simple disulfide: pdb=" SG CYS C 355 " - pdb=" SG CYS C 363 " distance=2.03 Simple disulfide: pdb=" SG CYS B 355 " - pdb=" SG CYS B 363 " distance=2.03 Simple disulfide: pdb=" SG CYS C 387 " - pdb=" SG CYS C 392 " distance=2.04 Simple disulfide: pdb=" SG CYS B 387 " - pdb=" SG CYS B 392 " distance=2.04 Simple disulfide: pdb=" SG CYS C 394 " - pdb=" SG CYS C 417 " distance=2.03 Simple disulfide: pdb=" SG CYS B 394 " - pdb=" SG CYS B 417 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 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.26 Conformation dependent library (CDL) restraints added in 1.9 seconds 3354 Ramachandran restraints generated. 1677 Oldfield, 0 Emsley, 1677 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 3138 Finding SS restraints... Secondary structure from input PDB file: 48 helices and 15 sheets defined 34.2% alpha, 21.3% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.62 Creating SS restraints... Processing helix chain 'A' and resid 21 through 26 Processing helix chain 'A' and resid 64 through 89 removed outlier: 4.545A pdb=" N ARG A 77 " --> pdb=" O ARG A 73 " (cutoff:3.500A) removed outlier: 4.116A pdb=" N LEU A 78 " --> pdb=" O LEU A 74 " (cutoff:3.500A) Proline residue: A 81 - end of helix removed outlier: 3.712A pdb=" N LEU A 88 " --> pdb=" O ASP A 84 " (cutoff:3.500A) removed outlier: 3.818A pdb=" N MET A 89 " --> pdb=" O GLY A 85 " (cutoff:3.500A) Processing helix chain 'A' and resid 125 through 138 Processing helix chain 'A' and resid 139 through 152 removed outlier: 4.568A pdb=" N ASP A 143 " --> pdb=" O GLN A 139 " (cutoff:3.500A) removed outlier: 3.829A pdb=" N GLU A 149 " --> pdb=" O GLU A 145 " (cutoff:3.500A) removed outlier: 4.031A pdb=" N ALA A 150 " --> pdb=" O LYS A 146 " (cutoff:3.500A) Processing helix chain 'A' and resid 175 through 182 Processing helix chain 'A' and resid 182 through 190 removed outlier: 4.141A pdb=" N ALA A 188 " --> pdb=" O VAL A 184 " (cutoff:3.500A) Processing helix chain 'A' and resid 190 through 215 removed outlier: 4.445A pdb=" N ALA A 194 " --> pdb=" O LEU A 190 " (cutoff:3.500A) Processing helix chain 'A' and resid 217 through 224 Processing helix chain 'A' and resid 228 through 236 Processing helix chain 'A' and resid 238 through 244 Processing helix chain 'A' and resid 253 through 262 Processing helix chain 'A' and resid 328 through 331 Processing helix chain 'A' and resid 351 through 357 Processing helix chain 'A' and resid 359 through 363 Processing helix chain 'A' and resid 452 through 487 Processing helix chain 'Y' and resid 86 through 90 Processing helix chain 'B' and resid 21 through 26 Processing helix chain 'B' and resid 64 through 89 removed outlier: 4.545A pdb=" N ARG B 77 " --> pdb=" O ARG B 73 " (cutoff:3.500A) removed outlier: 4.117A pdb=" N LEU B 78 " --> pdb=" O LEU B 74 " (cutoff:3.500A) Proline residue: B 81 - end of helix removed outlier: 3.712A pdb=" N LEU B 88 " --> pdb=" O ASP B 84 " (cutoff:3.500A) removed outlier: 3.818A pdb=" N MET B 89 " --> pdb=" O GLY B 85 " (cutoff:3.500A) Processing helix chain 'B' and resid 125 through 138 Processing helix chain 'B' and resid 139 through 152 removed outlier: 4.568A pdb=" N ASP B 143 " --> pdb=" O GLN B 139 " (cutoff:3.500A) removed outlier: 3.830A pdb=" N GLU B 149 " --> pdb=" O GLU B 145 " (cutoff:3.500A) removed outlier: 4.031A pdb=" N ALA B 150 " --> pdb=" O LYS B 146 " (cutoff:3.500A) Processing helix chain 'B' and resid 175 through 182 Processing helix chain 'B' and resid 182 through 190 removed outlier: 4.142A pdb=" N ALA B 188 " --> pdb=" O VAL B 184 " (cutoff:3.500A) Processing helix chain 'B' and resid 190 through 215 removed outlier: 4.445A pdb=" N ALA B 194 " --> pdb=" O LEU B 190 " (cutoff:3.500A) Processing helix chain 'B' and resid 217 through 224 Processing helix chain 'B' and resid 228 through 236 Processing helix chain 'B' and resid 238 through 244 Processing helix chain 'B' and resid 253 through 262 Processing helix chain 'B' and resid 328 through 331 Processing helix chain 'B' and resid 351 through 357 Processing helix chain 'B' and resid 359 through 363 Processing helix chain 'B' and resid 452 through 487 Processing helix chain 'W' and resid 86 through 90 Processing helix chain 'C' and resid 21 through 26 Processing helix chain 'C' and resid 64 through 89 removed outlier: 4.545A pdb=" N ARG C 77 " --> pdb=" O ARG C 73 " (cutoff:3.500A) removed outlier: 4.116A pdb=" N LEU C 78 " --> pdb=" O LEU C 74 " (cutoff:3.500A) Proline residue: C 81 - end of helix removed outlier: 3.712A pdb=" N LEU C 88 " --> pdb=" O ASP C 84 " (cutoff:3.500A) removed outlier: 3.820A pdb=" N MET C 89 " --> pdb=" O GLY C 85 " (cutoff:3.500A) Processing helix chain 'C' and resid 125 through 138 Processing helix chain 'C' and resid 139 through 152 removed outlier: 4.568A pdb=" N ASP C 143 " --> pdb=" O GLN C 139 " (cutoff:3.500A) removed outlier: 3.829A pdb=" N GLU C 149 " --> pdb=" O GLU C 145 " (cutoff:3.500A) removed outlier: 4.032A pdb=" N ALA C 150 " --> pdb=" O LYS C 146 " (cutoff:3.500A) Processing helix chain 'C' and resid 175 through 182 Processing helix chain 'C' and resid 182 through 190 removed outlier: 4.141A pdb=" N ALA C 188 " --> pdb=" O VAL C 184 " (cutoff:3.500A) Processing helix chain 'C' and resid 190 through 215 removed outlier: 4.446A pdb=" N ALA C 194 " --> pdb=" O LEU C 190 " (cutoff:3.500A) Processing helix chain 'C' and resid 217 through 224 Processing helix chain 'C' and resid 228 through 236 Processing helix chain 'C' and resid 238 through 244 Processing helix chain 'C' and resid 253 through 262 Processing helix chain 'C' and resid 328 through 331 Processing helix chain 'C' and resid 351 through 357 Processing helix chain 'C' and resid 359 through 363 Processing helix chain 'C' and resid 452 through 487 Processing helix chain 'U' and resid 86 through 90 Processing sheet with id=AA1, first strand: chain 'A' and resid 30 through 52 WARNING: can't find start of bonding for strands! previous: chain 'A' and resid 30 through 52 current: chain 'A' and resid 169 through 173 WARNING: can't find start of bonding for strands! previous: chain 'A' and resid 169 through 173 current: chain 'A' and resid 263 through 270 removed outlier: 7.304A pdb=" N VAL A 264 " --> pdb=" O ARG A 281 " (cutoff:3.500A) removed outlier: 7.793A pdb=" N ARG A 281 " --> pdb=" O VAL A 264 " (cutoff:3.500A) removed outlier: 5.969A pdb=" N VAL A 266 " --> pdb=" O GLN A 279 " (cutoff:3.500A) removed outlier: 6.328A pdb=" N GLN A 279 " --> pdb=" O VAL A 266 " (cutoff:3.500A) removed outlier: 4.492A pdb=" N HIS A 315 " --> pdb=" O ALA A 326 " (cutoff:3.500A) WARNING: can't find start of bonding for strands! previous: chain 'A' and resid 322 through 326 current: chain 'A' and resid 337 through 339 WARNING: can't find start of bonding for strands! previous: chain 'A' and resid 337 through 339 current: chain 'Y' and resid 46 through 52 WARNING: can't find start of bonding for strands! previous: chain 'Y' and resid 57 through 59 current: chain 'Y' and resid 103 through 106 Processing sheet with id=AA2, first strand: chain 'A' and resid 92 through 93 removed outlier: 3.986A pdb=" N GLY A 116 " --> pdb=" O LEU B 426 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 122 through 124 Processing sheet with id=AA4, first strand: chain 'A' and resid 301 through 303 Processing sheet with id=AA5, first strand: chain 'A' and resid 410 through 412 Processing sheet with id=AA6, first strand: chain 'A' and resid 392 through 393 removed outlier: 7.232A pdb=" N LEU A 426 " --> pdb=" O GLY C 116 " (cutoff:3.500A) removed outlier: 7.865A pdb=" N ILE C 118 " --> pdb=" O LEU A 426 " (cutoff:3.500A) removed outlier: 7.182A pdb=" N ASN A 428 " --> pdb=" O ILE C 118 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'Y' and resid 2 through 5 removed outlier: 3.537A pdb=" N THR Y 77 " --> pdb=" O ASP Y 72 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 30 through 52 WARNING: can't find start of bonding for strands! previous: chain 'B' and resid 30 through 52 current: chain 'B' and resid 169 through 173 WARNING: can't find start of bonding for strands! previous: chain 'B' and resid 169 through 173 current: chain 'B' and resid 263 through 270 removed outlier: 7.303A pdb=" N VAL B 264 " --> pdb=" O ARG B 281 " (cutoff:3.500A) removed outlier: 7.793A pdb=" N ARG B 281 " --> pdb=" O VAL B 264 " (cutoff:3.500A) removed outlier: 5.969A pdb=" N VAL B 266 " --> pdb=" O GLN B 279 " (cutoff:3.500A) removed outlier: 6.327A pdb=" N GLN B 279 " --> pdb=" O VAL B 266 " (cutoff:3.500A) removed outlier: 4.492A pdb=" N HIS B 315 " --> pdb=" O ALA B 326 " (cutoff:3.500A) WARNING: can't find start of bonding for strands! previous: chain 'B' and resid 322 through 326 current: chain 'B' and resid 337 through 339 WARNING: can't find start of bonding for strands! previous: chain 'B' and resid 337 through 339 current: chain 'W' and resid 46 through 52 WARNING: can't find start of bonding for strands! previous: chain 'W' and resid 57 through 59 current: chain 'W' and resid 103 through 106 Processing sheet with id=AA9, first strand: chain 'B' and resid 92 through 93 removed outlier: 4.035A pdb=" N GLY B 116 " --> pdb=" O LEU C 426 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'B' and resid 122 through 124 Processing sheet with id=AB2, first strand: chain 'B' and resid 301 through 303 Processing sheet with id=AB3, first strand: chain 'W' and resid 2 through 5 removed outlier: 3.537A pdb=" N THR W 77 " --> pdb=" O ASP W 72 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 30 through 52 WARNING: can't find start of bonding for strands! previous: chain 'C' and resid 30 through 52 current: chain 'C' and resid 169 through 173 WARNING: can't find start of bonding for strands! previous: chain 'C' and resid 169 through 173 current: chain 'C' and resid 263 through 270 removed outlier: 7.304A pdb=" N VAL C 264 " --> pdb=" O ARG C 281 " (cutoff:3.500A) removed outlier: 7.792A pdb=" N ARG C 281 " --> pdb=" O VAL C 264 " (cutoff:3.500A) removed outlier: 5.970A pdb=" N VAL C 266 " --> pdb=" O GLN C 279 " (cutoff:3.500A) removed outlier: 6.327A pdb=" N GLN C 279 " --> pdb=" O VAL C 266 " (cutoff:3.500A) removed outlier: 4.492A pdb=" N HIS C 315 " --> pdb=" O ALA C 326 " (cutoff:3.500A) WARNING: can't find start of bonding for strands! previous: chain 'C' and resid 322 through 326 current: chain 'C' and resid 337 through 339 WARNING: can't find start of bonding for strands! previous: chain 'C' and resid 337 through 339 current: chain 'U' and resid 46 through 52 WARNING: can't find start of bonding for strands! previous: chain 'U' and resid 57 through 59 current: chain 'U' and resid 103 through 106 Processing sheet with id=AB5, first strand: chain 'C' and resid 301 through 303 Processing sheet with id=AB6, first strand: chain 'U' and resid 2 through 5 removed outlier: 3.537A pdb=" N THR U 77 " --> pdb=" O ASP U 72 " (cutoff:3.500A) 523 hydrogen bonds defined for protein. 1425 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: 4.13 Time building geometry restraints manager: 3.64 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 3593 1.33 - 1.45: 2143 1.45 - 1.57: 7440 1.57 - 1.70: 0 1.70 - 1.82: 72 Bond restraints: 13248 Sorted by residual: bond pdb=" CG1 ILE B 262 " pdb=" CD1 ILE B 262 " ideal model delta sigma weight residual 1.513 1.431 0.082 3.90e-02 6.57e+02 4.43e+00 bond pdb=" CG1 ILE C 262 " pdb=" CD1 ILE C 262 " ideal model delta sigma weight residual 1.513 1.431 0.082 3.90e-02 6.57e+02 4.42e+00 bond pdb=" CG1 ILE A 262 " pdb=" CD1 ILE A 262 " ideal model delta sigma weight residual 1.513 1.431 0.082 3.90e-02 6.57e+02 4.37e+00 bond pdb=" CB ARG B 365 " pdb=" CG ARG B 365 " ideal model delta sigma weight residual 1.520 1.457 0.063 3.00e-02 1.11e+03 4.36e+00 bond pdb=" CB ARG C 365 " pdb=" CG ARG C 365 " ideal model delta sigma weight residual 1.520 1.458 0.062 3.00e-02 1.11e+03 4.32e+00 ... (remaining 13243 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.80: 16584 1.80 - 3.60: 1145 3.60 - 5.39: 199 5.39 - 7.19: 30 7.19 - 8.99: 9 Bond angle restraints: 17967 Sorted by residual: angle pdb=" N SER C 174 " pdb=" CA SER C 174 " pdb=" C SER C 174 " ideal model delta sigma weight residual 111.37 116.86 -5.49 1.64e+00 3.72e-01 1.12e+01 angle pdb=" N SER A 174 " pdb=" CA SER A 174 " pdb=" C SER A 174 " ideal model delta sigma weight residual 111.37 116.84 -5.47 1.64e+00 3.72e-01 1.11e+01 angle pdb=" N SER B 174 " pdb=" CA SER B 174 " pdb=" C SER B 174 " ideal model delta sigma weight residual 111.37 116.81 -5.44 1.64e+00 3.72e-01 1.10e+01 angle pdb=" CA PRO B 81 " pdb=" N PRO B 81 " pdb=" CD PRO B 81 " ideal model delta sigma weight residual 112.00 107.39 4.61 1.40e+00 5.10e-01 1.09e+01 angle pdb=" CA TYR W 59 " pdb=" CB TYR W 59 " pdb=" CG TYR W 59 " ideal model delta sigma weight residual 113.90 119.81 -5.91 1.80e+00 3.09e-01 1.08e+01 ... (remaining 17962 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.54: 7500 17.54 - 35.08: 483 35.08 - 52.62: 90 52.62 - 70.16: 15 70.16 - 87.70: 18 Dihedral angle restraints: 8106 sinusoidal: 3219 harmonic: 4887 Sorted by residual: dihedral pdb=" CA LYS A 173 " pdb=" C LYS A 173 " pdb=" N SER A 174 " pdb=" CA SER A 174 " ideal model delta harmonic sigma weight residual 180.00 -152.63 -27.37 0 5.00e+00 4.00e-02 3.00e+01 dihedral pdb=" CA LYS C 173 " pdb=" C LYS C 173 " pdb=" N SER C 174 " pdb=" CA SER C 174 " ideal model delta harmonic sigma weight residual -180.00 -152.66 -27.34 0 5.00e+00 4.00e-02 2.99e+01 dihedral pdb=" CA LYS B 173 " pdb=" C LYS B 173 " pdb=" N SER B 174 " pdb=" CA SER B 174 " ideal model delta harmonic sigma weight residual -180.00 -152.70 -27.30 0 5.00e+00 4.00e-02 2.98e+01 ... (remaining 8103 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.066: 1538 0.066 - 0.132: 541 0.132 - 0.198: 43 0.198 - 0.264: 5 0.264 - 0.330: 9 Chirality restraints: 2136 Sorted by residual: chirality pdb=" CG LEU U 4 " pdb=" CB LEU U 4 " pdb=" CD1 LEU U 4 " pdb=" CD2 LEU U 4 " both_signs ideal model delta sigma weight residual False -2.59 -2.26 -0.33 2.00e-01 2.50e+01 2.72e+00 chirality pdb=" CG LEU W 4 " pdb=" CB LEU W 4 " pdb=" CD1 LEU W 4 " pdb=" CD2 LEU W 4 " both_signs ideal model delta sigma weight residual False -2.59 -2.26 -0.33 2.00e-01 2.50e+01 2.72e+00 chirality pdb=" CG LEU Y 4 " pdb=" CB LEU Y 4 " pdb=" CD1 LEU Y 4 " pdb=" CD2 LEU Y 4 " both_signs ideal model delta sigma weight residual False -2.59 -2.26 -0.33 2.00e-01 2.50e+01 2.70e+00 ... (remaining 2133 not shown) Planarity restraints: 2307 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C ILE B 80 " 0.073 5.00e-02 4.00e+02 1.10e-01 1.92e+01 pdb=" N PRO B 81 " -0.189 5.00e-02 4.00e+02 pdb=" CA PRO B 81 " 0.056 5.00e-02 4.00e+02 pdb=" CD PRO B 81 " 0.060 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C ILE C 80 " -0.073 5.00e-02 4.00e+02 1.09e-01 1.92e+01 pdb=" N PRO C 81 " 0.189 5.00e-02 4.00e+02 pdb=" CA PRO C 81 " -0.056 5.00e-02 4.00e+02 pdb=" CD PRO C 81 " -0.060 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C ILE A 80 " 0.073 5.00e-02 4.00e+02 1.09e-01 1.92e+01 pdb=" N PRO A 81 " -0.189 5.00e-02 4.00e+02 pdb=" CA PRO A 81 " 0.056 5.00e-02 4.00e+02 pdb=" CD PRO A 81 " 0.060 5.00e-02 4.00e+02 ... (remaining 2304 not shown) Histogram of nonbonded interaction distances: 2.03 - 2.61: 102 2.61 - 3.18: 11732 3.18 - 3.75: 19833 3.75 - 4.33: 28650 4.33 - 4.90: 47515 Nonbonded interactions: 107832 Sorted by model distance: nonbonded pdb=" SG CYS W 22 " pdb=" SG CYS W 95 " model vdw 2.032 3.760 nonbonded pdb=" SG CYS U 22 " pdb=" SG CYS U 95 " model vdw 2.032 3.760 nonbonded pdb=" NE2 GLN A 139 " pdb=" O GLU Y 1 " model vdw 2.106 3.120 nonbonded pdb=" OE1 GLU C 136 " pdb=" OG SER U 99 " model vdw 2.219 3.040 nonbonded pdb=" OE1 GLU A 136 " pdb=" OG SER Y 99 " model vdw 2.242 3.040 ... (remaining 107827 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.09 Found NCS groups: ncs_group { reference = chain 'U' selection = chain 'W' selection = chain 'Y' } ncs_group { reference = chain 'A' selection = chain 'C' selection = chain 'B' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators 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 individual isotropic 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: 3.120 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.100 Construct map_model_manager: 0.010 Extract box with map and model: 0.570 Check model and map are aligned: 0.100 Set scattering table: 0.140 Process input model: 32.690 Find NCS groups from input model: 0.270 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:13.350 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 50.400 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7003 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.009 0.082 13264 Z= 0.338 Angle : 1.047 8.992 17999 Z= 0.580 Chirality : 0.064 0.330 2136 Planarity : 0.007 0.110 2307 Dihedral : 13.034 87.697 4920 Min Nonbonded Distance : 2.032 Molprobity Statistics. All-atom Clashscore : 6.89 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.70 % Favored : 96.30 % Rotamer: Outliers : 0.00 % Allowed : 0.41 % Favored : 99.59 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.64 (0.20), residues: 1677 helix: -0.35 (0.22), residues: 492 sheet: -0.23 (0.25), residues: 393 loop : -0.42 (0.23), residues: 792 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.032 0.005 TRP C 308 HIS 0.015 0.005 HIS C 350 PHE 0.024 0.004 PHE W 67 TYR 0.042 0.006 TYR A 438 ARG 0.016 0.003 ARG A 77 Details of bonding type rmsd hydrogen bonds : bond 0.16649 ( 523) hydrogen bonds : angle 7.62169 ( 1425) SS BOND : bond 0.00473 ( 16) SS BOND : angle 1.57536 ( 32) covalent geometry : bond 0.00855 (13248) covalent geometry : angle 1.04535 (17967) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3354 Ramachandran restraints generated. 1677 Oldfield, 0 Emsley, 1677 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3354 Ramachandran restraints generated. 1677 Oldfield, 0 Emsley, 1677 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 338 residues out of total 1467 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 338 time to evaluate : 1.380 Fit side-chains revert: symmetry clash REVERT: A 362 GLN cc_start: 0.7622 (mm-40) cc_final: 0.7382 (mm110) REVERT: A 413 THR cc_start: 0.7865 (p) cc_final: 0.7643 (m) REVERT: Y 31 ILE cc_start: 0.8385 (mm) cc_final: 0.8166 (mp) REVERT: B 84 ASP cc_start: 0.6974 (m-30) cc_final: 0.6669 (m-30) REVERT: W 99 SER cc_start: 0.7601 (t) cc_final: 0.7085 (t) REVERT: C 148 LYS cc_start: 0.7388 (mttt) cc_final: 0.7174 (mmmm) REVERT: C 362 GLN cc_start: 0.7581 (mm-40) cc_final: 0.7369 (mm-40) REVERT: C 413 THR cc_start: 0.7945 (p) cc_final: 0.7653 (m) outliers start: 0 outliers final: 5 residues processed: 338 average time/residue: 1.3199 time to fit residues: 483.7993 Evaluate side-chains 325 residues out of total 1467 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 320 time to evaluate : 1.507 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 218 ILE Chi-restraints excluded: chain A residue 401 ILE Chi-restraints excluded: chain B residue 218 ILE Chi-restraints excluded: chain B residue 474 ILE Chi-restraints excluded: chain C residue 218 ILE Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 168 random chunks: chunk 141 optimal weight: 0.0970 chunk 127 optimal weight: 0.9990 chunk 70 optimal weight: 0.5980 chunk 43 optimal weight: 0.6980 chunk 85 optimal weight: 0.0980 chunk 68 optimal weight: 0.6980 chunk 131 optimal weight: 0.6980 chunk 50 optimal weight: 0.7980 chunk 80 optimal weight: 0.7980 chunk 98 optimal weight: 0.0870 chunk 152 optimal weight: 0.6980 overall best weight: 0.3156 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 176 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 180 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 350 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 399 ASN ** B 176 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 180 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 198 GLN B 229 GLN B 292 GLN ** B 350 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** B 399 ASN B 407 GLN ** W 111 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 176 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 180 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 292 GLN U 103 ASN U 111 GLN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3859 r_free = 0.3859 target = 0.145173 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3545 r_free = 0.3545 target = 0.122640 restraints weight = 18402.589| |-----------------------------------------------------------------------------| r_work (start): 0.3549 rms_B_bonded: 2.11 r_work: 0.3422 rms_B_bonded: 2.71 restraints_weight: 0.5000 r_work: 0.3247 rms_B_bonded: 4.65 restraints_weight: 0.2500 r_work (final): 0.3247 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7596 moved from start: 0.1561 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 13264 Z= 0.140 Angle : 0.641 7.667 17999 Z= 0.338 Chirality : 0.046 0.147 2136 Planarity : 0.005 0.069 2307 Dihedral : 5.424 29.440 1834 Min Nonbonded Distance : 2.568 Molprobity Statistics. All-atom Clashscore : 6.24 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.43 % Favored : 98.57 % Rotamer: Outliers : 2.39 % Allowed : 9.41 % Favored : 88.21 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.05 (0.21), residues: 1677 helix: 0.59 (0.23), residues: 495 sheet: -0.43 (0.24), residues: 411 loop : 0.12 (0.25), residues: 771 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.002 TRP A 308 HIS 0.006 0.001 HIS B 350 PHE 0.012 0.001 PHE U 67 TYR 0.021 0.002 TYR C 48 ARG 0.005 0.001 ARG B 87 Details of bonding type rmsd hydrogen bonds : bond 0.04556 ( 523) hydrogen bonds : angle 5.65132 ( 1425) SS BOND : bond 0.00141 ( 16) SS BOND : angle 0.58356 ( 32) covalent geometry : bond 0.00323 (13248) covalent geometry : angle 0.64109 (17967) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3354 Ramachandran restraints generated. 1677 Oldfield, 0 Emsley, 1677 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3354 Ramachandran restraints generated. 1677 Oldfield, 0 Emsley, 1677 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 378 residues out of total 1467 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 35 poor density : 343 time to evaluate : 1.432 Fit side-chains REVERT: A 90 LYS cc_start: 0.6495 (mppt) cc_final: 0.6048 (mppt) REVERT: A 362 GLN cc_start: 0.8212 (mm-40) cc_final: 0.7967 (mm110) REVERT: A 422 ILE cc_start: 0.8322 (OUTLIER) cc_final: 0.7948 (mt) REVERT: A 472 GLU cc_start: 0.7643 (mm-30) cc_final: 0.7315 (mm-30) REVERT: Y 31 ILE cc_start: 0.8444 (mm) cc_final: 0.8224 (mp) REVERT: B 21 ASP cc_start: 0.8061 (t0) cc_final: 0.7822 (t0) REVERT: B 73 ARG cc_start: 0.8159 (ttt-90) cc_final: 0.7902 (ttt180) REVERT: B 84 ASP cc_start: 0.7984 (m-30) cc_final: 0.7774 (m-30) REVERT: B 319 LYS cc_start: 0.8217 (mtpp) cc_final: 0.7956 (ttpt) REVERT: B 323 LEU cc_start: 0.8248 (tp) cc_final: 0.7810 (tt) REVERT: B 353 GLU cc_start: 0.7104 (tp30) cc_final: 0.6618 (tp30) REVERT: B 375 ARG cc_start: 0.8609 (mpt180) cc_final: 0.8398 (mpt180) REVERT: B 422 ILE cc_start: 0.8369 (OUTLIER) cc_final: 0.8040 (mt) REVERT: W 29 SER cc_start: 0.8391 (p) cc_final: 0.8129 (m) REVERT: C 24 LYS cc_start: 0.8440 (mtpt) cc_final: 0.8158 (mtpt) REVERT: C 227 LYS cc_start: 0.8394 (mtpp) cc_final: 0.8155 (mtpp) REVERT: C 413 THR cc_start: 0.8663 (p) cc_final: 0.8340 (m) REVERT: C 464 LYS cc_start: 0.8428 (mtpt) cc_final: 0.8142 (mtmp) REVERT: C 468 GLU cc_start: 0.8051 (mm-30) cc_final: 0.7840 (tp30) outliers start: 35 outliers final: 15 residues processed: 349 average time/residue: 1.3029 time to fit residues: 492.1691 Evaluate side-chains 331 residues out of total 1467 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 17 poor density : 314 time to evaluate : 1.370 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 23 THR Chi-restraints excluded: chain A residue 328 VAL Chi-restraints excluded: chain A residue 360 ILE Chi-restraints excluded: chain A residue 420 ILE Chi-restraints excluded: chain A residue 422 ILE Chi-restraints excluded: chain A residue 474 ILE Chi-restraints excluded: chain Y residue 27 SER Chi-restraints excluded: chain B residue 23 THR Chi-restraints excluded: chain B residue 328 VAL Chi-restraints excluded: chain B residue 330 GLU Chi-restraints excluded: chain B residue 413 THR Chi-restraints excluded: chain B residue 422 ILE Chi-restraints excluded: chain B residue 474 ILE Chi-restraints excluded: chain W residue 53 SER Chi-restraints excluded: chain C residue 420 ILE Chi-restraints excluded: chain C residue 422 ILE Chi-restraints excluded: chain U residue 31 ILE Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/iotbx/cli_parser.py", line 994, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 767, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1525, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1427, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1308, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 65.1021 > 50: distance: 41 - 162: 6.554 distance: 44 - 159: 9.552 distance: 75 - 94: 4.834 distance: 87 - 140: 5.886 distance: 94 - 95: 6.443 distance: 95 - 96: 6.520 distance: 95 - 98: 10.941 distance: 96 - 97: 26.567 distance: 96 - 100: 35.002 distance: 98 - 99: 32.797 distance: 100 - 101: 13.295 distance: 101 - 102: 19.343 distance: 102 - 103: 14.873 distance: 102 - 104: 11.907 distance: 104 - 105: 7.021 distance: 105 - 106: 8.921 distance: 105 - 108: 6.363 distance: 106 - 107: 8.692 distance: 106 - 112: 5.349 distance: 107 - 126: 10.211 distance: 108 - 109: 5.624 distance: 109 - 110: 11.357 distance: 109 - 111: 15.713 distance: 112 - 113: 11.919 distance: 113 - 114: 3.213 distance: 113 - 116: 6.587 distance: 114 - 115: 8.520 distance: 114 - 120: 5.149 distance: 115 - 135: 13.347 distance: 116 - 117: 8.994 distance: 116 - 118: 7.471 distance: 117 - 119: 22.572 distance: 120 - 121: 9.550 distance: 121 - 122: 6.422 distance: 121 - 124: 13.995 distance: 122 - 123: 30.642 distance: 122 - 126: 5.766 distance: 124 - 125: 16.737 distance: 126 - 127: 4.385 distance: 127 - 128: 3.285 distance: 127 - 130: 15.234 distance: 128 - 129: 6.558 distance: 128 - 135: 9.678 distance: 130 - 131: 17.893 distance: 131 - 132: 15.748 distance: 132 - 133: 4.498 distance: 132 - 134: 4.442 distance: 135 - 136: 6.199 distance: 136 - 137: 14.497 distance: 136 - 139: 7.661 distance: 137 - 138: 9.628 distance: 137 - 141: 8.841 distance: 139 - 140: 9.313 distance: 141 - 142: 4.484 distance: 141 - 147: 10.474 distance: 142 - 143: 7.921 distance: 142 - 145: 9.729 distance: 143 - 144: 10.811 distance: 143 - 148: 12.015 distance: 145 - 146: 9.769 distance: 146 - 147: 4.715 distance: 148 - 149: 6.463 distance: 149 - 152: 5.746 distance: 150 - 151: 6.311 distance: 150 - 159: 3.531 distance: 152 - 153: 10.043 distance: 153 - 154: 3.447 distance: 154 - 155: 3.438 distance: 156 - 157: 4.514 distance: 156 - 158: 6.075 distance: 159 - 160: 8.083 distance: 160 - 161: 3.440 distance: 160 - 163: 5.452 distance: 161 - 162: 3.667 distance: 163 - 164: 9.652 distance: 163 - 165: 9.231