Starting phenix.real_space_refine on Thu Mar 13 02:28:19 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/8wfs_37498/03_2025/8wfs_37498.cif Found real_map, /net/cci-nas-00/data/ceres_data/8wfs_37498/03_2025/8wfs_37498.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.36 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/8wfs_37498/03_2025/8wfs_37498.map" default_real_map = "/net/cci-nas-00/data/ceres_data/8wfs_37498/03_2025/8wfs_37498.map" model { file = "/net/cci-nas-00/data/ceres_data/8wfs_37498/03_2025/8wfs_37498.cif" } default_model = "/net/cci-nas-00/data/ceres_data/8wfs_37498/03_2025/8wfs_37498.cif" } resolution = 3.36 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.000 sd= 0.015 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 2864 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 3973 2.51 5 N 1123 2.21 5 O 1124 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 22 residue(s): 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/modules/chem_data/mon_lib" Total number of atoms: 6274 Number of models: 1 Model: "" Number of chains: 13 Chain: "c" Number of atoms: 884 Number of conformers: 1 Conformer: "" Number of residues, atoms: 118, 884 Classifications: {'peptide': 118} Link IDs: {'PCIS': 1, 'PTRANS': 14, 'TRANS': 102} Chain: "C" Number of atoms: 882 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 117, 877 Classifications: {'peptide': 117} Link IDs: {'PCIS': 1, 'PTRANS': 13, 'TRANS': 102} Conformer: "B" Number of residues, atoms: 117, 877 Classifications: {'peptide': 117} Link IDs: {'PCIS': 1, 'PTRANS': 13, 'TRANS': 102} bond proxies already assigned to first conformer: 894 Chain: "E" Number of atoms: 957 Number of conformers: 1 Conformer: "" Number of residues, atoms: 125, 957 Classifications: {'peptide': 125} Link IDs: {'PCIS': 1, 'PTRANS': 10, 'TRANS': 113} Chain: "B" Number of atoms: 1035 Number of conformers: 1 Conformer: "" Number of residues, atoms: 140, 1035 Classifications: {'peptide': 140} Link IDs: {'PCIS': 1, 'PTRANS': 15, 'TRANS': 123} Chain: "b" Number of atoms: 1035 Number of conformers: 1 Conformer: "" Number of residues, atoms: 140, 1035 Classifications: {'peptide': 140} Link IDs: {'PCIS': 1, 'PTRANS': 15, 'TRANS': 123} Chain: "X" Number of atoms: 1112 Number of conformers: 1 Conformer: "" Number of residues, atoms: 146, 1112 Classifications: {'peptide': 146} Link IDs: {'PCIS': 1, 'PTRANS': 11, 'TRANS': 133} Chain: "a" Number of atoms: 179 Number of conformers: 1 Conformer: "" Number of residues, atoms: 22, 179 Classifications: {'peptide': 22} Link IDs: {'PTRANS': 1, 'TRANS': 20} Chain: "A" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "D" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "F" Number of atoms: 39 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 39 Unusual residues: {'BMA': 1, 'NAG': 2} Classifications: {'undetermined': 3} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 3 Unresolved non-hydrogen angles: 6 Unresolved non-hydrogen dihedrals: 9 Unresolved non-hydrogen chiralities: 3 Chain: "G" Number of atoms: 39 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 39 Unusual residues: {'BMA': 1, 'NAG': 2} Classifications: {'undetermined': 3} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 3 Unresolved non-hydrogen angles: 6 Unresolved non-hydrogen dihedrals: 9 Unresolved non-hydrogen chiralities: 3 Chain: "H" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "I" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Residues with excluded nonbonded symmetry interactions: 1 residue: pdb=" N AALA C 114 " occ=0.50 ... (8 atoms not shown) pdb=" CB BALA C 114 " occ=0.50 Time building chain proxies: 5.23, per 1000 atoms: 0.83 Number of scatterers: 6274 At special positions: 0 Unit cell: (83.7876, 97.7522, 134.275, 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 1124 8.00 N 1123 7.00 C 3973 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=26, symmetry=0 Simple disulfide: pdb=" SG CYS c 1 " - pdb=" SG CYS c 7 " distance=2.03 Simple disulfide: pdb=" SG CYS c 5 " - pdb=" SG CYS c 14 " distance=2.03 Simple disulfide: pdb=" SG CYS c 68 " - pdb=" SG CYS c 93 " distance=2.04 Simple disulfide: pdb=" SG CYS c 70 " - pdb=" SG CYS c 116 " distance=2.03 Simple disulfide: pdb=" SG CYS C 1 " - pdb=" SG CYS C 7 " distance=2.03 Simple disulfide: pdb=" SG CYS C 5 " - pdb=" SG CYS C 14 " distance=2.03 Simple disulfide: pdb=" SG CYS C 68 " - pdb=" SG CYS C 93 " distance=2.03 Simple disulfide: pdb=" SG CYS C 70 " - pdb=" SG CYS C 116 " distance=2.03 Simple disulfide: pdb=" SG CYS E 4 " - pdb=" SG CYS E 10 " distance=2.03 Simple disulfide: pdb=" SG CYS E 8 " - pdb=" SG CYS E 22 " distance=2.03 Simple disulfide: pdb=" SG CYS E 73 " - pdb=" SG CYS E 97 " distance=2.03 Simple disulfide: pdb=" SG CYS E 75 " - pdb=" SG CYS E 119 " distance=2.03 Simple disulfide: pdb=" SG CYS B 1 " - pdb=" SG CYS B 7 " distance=2.03 Simple disulfide: pdb=" SG CYS B 5 " - pdb=" SG CYS B 14 " distance=2.03 Simple disulfide: pdb=" SG CYS B 68 " - pdb=" SG CYS B 93 " distance=2.03 Simple disulfide: pdb=" SG CYS B 70 " - pdb=" SG CYS B 116 " distance=2.03 Simple disulfide: pdb=" SG CYS B 122 " - pdb=" SG CYS a 484 " distance=2.03 Simple disulfide: pdb=" SG CYS b 1 " - pdb=" SG CYS b 7 " distance=2.03 Simple disulfide: pdb=" SG CYS b 5 " - pdb=" SG CYS b 14 " distance=2.03 Simple disulfide: pdb=" SG CYS b 68 " - pdb=" SG CYS b 93 " distance=2.03 Simple disulfide: pdb=" SG CYS b 70 " - pdb=" SG CYS b 116 " distance=2.03 Simple disulfide: pdb=" SG CYS b 122 " - pdb=" SG CYS a 485 " distance=2.03 Simple disulfide: pdb=" SG CYS X 4 " - pdb=" SG CYS X 10 " distance=2.03 Simple disulfide: pdb=" SG CYS X 8 " - pdb=" SG CYS X 22 " distance=2.03 Simple disulfide: pdb=" SG CYS X 73 " - pdb=" SG CYS X 97 " distance=2.03 Simple disulfide: pdb=" SG CYS X 75 " - pdb=" SG CYS X 119 " 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 Links applied BETA1-4 " NAG A 1 " - " NAG A 2 " " NAG D 1 " - " NAG D 2 " " NAG F 1 " - " NAG F 2 " " NAG F 2 " - " BMA F 3 " " NAG G 1 " - " NAG G 2 " " NAG G 2 " - " BMA G 3 " " NAG H 1 " - " NAG H 2 " " NAG I 1 " - " NAG I 2 " NAG-ASN " NAG A 1 " - " ASN c 41 " " NAG D 1 " - " ASN C 41 " " NAG F 1 " - " ASN E 44 " " NAG G 1 " - " ASN B 41 " " NAG H 1 " - " ASN b 41 " " NAG I 1 " - " ASN X 44 " Time building additional restraints: 2.58 Conformation dependent library (CDL) restraints added in 2.4 seconds 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1494 Finding SS restraints... Secondary structure from input PDB file: 33 helices and 11 sheets defined 31.2% alpha, 8.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.73 Creating SS restraints... Processing helix chain 'c' and resid 49 through 53 removed outlier: 3.545A pdb=" N ALA c 52 " --> pdb=" O LEU c 49 " (cutoff:3.500A) Processing helix chain 'c' and resid 69 through 71 No H-bonds generated for 'chain 'c' and resid 69 through 71' Processing helix chain 'c' and resid 72 through 81 removed outlier: 4.190A pdb=" N GLY c 81 " --> pdb=" O ALA c 77 " (cutoff:3.500A) Processing helix chain 'c' and resid 85 through 89 removed outlier: 3.541A pdb=" N TYR c 88 " --> pdb=" O ARG c 85 " (cutoff:3.500A) removed outlier: 3.913A pdb=" N ARG c 89 " --> pdb=" O ALA c 86 " (cutoff:3.500A) No H-bonds generated for 'chain 'c' and resid 85 through 89' Processing helix chain 'c' and resid 96 through 100 removed outlier: 3.662A pdb=" N ARG c 100 " --> pdb=" O PRO c 97 " (cutoff:3.500A) Processing helix chain 'c' and resid 108 through 115 Processing helix chain 'C' and resid 69 through 71 No H-bonds generated for 'chain 'C' and resid 69 through 71' Processing helix chain 'C' and resid 72 through 81 Processing helix chain 'C' and resid 96 through 100 removed outlier: 3.651A pdb=" N ARG C 100 " --> pdb=" O PRO C 97 " (cutoff:3.500A) Processing helix chain 'C' and resid 103 through 107 Processing helix chain 'C' and resid 108 through 115 removed outlier: 3.754A pdb=" N LEU C 112 " --> pdb=" O ALA C 108 " (cutoff:3.500A) removed outlier: 3.861A pdb=" N ALA C 115 " --> pdb=" O GLU C 111 " (cutoff:3.500A) Processing helix chain 'E' and resid 74 through 76 No H-bonds generated for 'chain 'E' and resid 74 through 76' Processing helix chain 'E' and resid 77 through 85 Processing helix chain 'E' and resid 88 through 95 removed outlier: 4.089A pdb=" N VAL E 95 " --> pdb=" O LEU E 92 " (cutoff:3.500A) Processing helix chain 'B' and resid 69 through 71 No H-bonds generated for 'chain 'B' and resid 69 through 71' Processing helix chain 'B' and resid 72 through 80 removed outlier: 3.594A pdb=" N ALA B 80 " --> pdb=" O ARG B 76 " (cutoff:3.500A) Processing helix chain 'B' and resid 85 through 89 removed outlier: 3.524A pdb=" N TYR B 88 " --> pdb=" O ARG B 85 " (cutoff:3.500A) Processing helix chain 'B' and resid 96 through 100 removed outlier: 3.589A pdb=" N ARG B 100 " --> pdb=" O PRO B 97 " (cutoff:3.500A) Processing helix chain 'B' and resid 103 through 107 Processing helix chain 'B' and resid 108 through 117 removed outlier: 4.135A pdb=" N ALA B 115 " --> pdb=" O GLU B 111 " (cutoff:3.500A) Processing helix chain 'B' and resid 125 through 140 removed outlier: 3.829A pdb=" N GLN B 129 " --> pdb=" O ALA B 125 " (cutoff:3.500A) removed outlier: 3.993A pdb=" N ALA B 131 " --> pdb=" O ALA B 127 " (cutoff:3.500A) Processing helix chain 'b' and resid 49 through 53 removed outlier: 3.908A pdb=" N LEU b 53 " --> pdb=" O LEU b 50 " (cutoff:3.500A) Processing helix chain 'b' and resid 69 through 71 No H-bonds generated for 'chain 'b' and resid 69 through 71' Processing helix chain 'b' and resid 72 through 81 removed outlier: 3.505A pdb=" N GLY b 81 " --> pdb=" O ALA b 77 " (cutoff:3.500A) Processing helix chain 'b' and resid 96 through 100 removed outlier: 3.551A pdb=" N ARG b 100 " --> pdb=" O PRO b 97 " (cutoff:3.500A) Processing helix chain 'b' and resid 103 through 107 Processing helix chain 'b' and resid 108 through 116 removed outlier: 3.798A pdb=" N LEU b 112 " --> pdb=" O ALA b 108 " (cutoff:3.500A) removed outlier: 3.927A pdb=" N CYS b 116 " --> pdb=" O LEU b 112 " (cutoff:3.500A) Processing helix chain 'b' and resid 125 through 140 Processing helix chain 'X' and resid 74 through 76 No H-bonds generated for 'chain 'X' and resid 74 through 76' Processing helix chain 'X' and resid 77 through 88 Processing helix chain 'X' and resid 89 through 94 Processing helix chain 'X' and resid 130 through 146 removed outlier: 3.662A pdb=" N ALA X 137 " --> pdb=" O LEU X 133 " (cutoff:3.500A) removed outlier: 3.839A pdb=" N GLY X 145 " --> pdb=" O VAL X 141 " (cutoff:3.500A) removed outlier: 3.919A pdb=" N LEU X 146 " --> pdb=" O ALA X 142 " (cutoff:3.500A) Processing helix chain 'a' and resid 483 through 503 removed outlier: 3.688A pdb=" N LEU a 487 " --> pdb=" O PHE a 483 " (cutoff:3.500A) Proline residue: a 488 - end of helix removed outlier: 4.359A pdb=" N PHE a 491 " --> pdb=" O LEU a 487 " (cutoff:3.500A) removed outlier: 3.556A pdb=" N VAL a 493 " --> pdb=" O LEU a 489 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'c' and resid 6 through 8 Processing sheet with id=AA2, first strand: chain 'c' and resid 66 through 67 removed outlier: 7.480A pdb=" N TRP c 66 " --> pdb=" O VAL c 94 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA2 Processing sheet with id=AA3, first strand: chain 'C' and resid 6 through 8 removed outlier: 3.658A pdb=" N VAL C 36 " --> pdb=" O VAL C 12 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'C' and resid 66 through 67 removed outlier: 7.132A pdb=" N TRP C 66 " --> pdb=" O VAL C 94 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA4 Processing sheet with id=AA5, first strand: chain 'E' and resid 9 through 13 removed outlier: 3.660A pdb=" N GLY E 17 " --> pdb=" O LEU E 13 " (cutoff:3.500A) removed outlier: 5.594A pdb=" N LEU E 18 " --> pdb=" O HIS E 39 " (cutoff:3.500A) removed outlier: 7.273A pdb=" N LEU E 41 " --> pdb=" O LEU E 18 " (cutoff:3.500A) removed outlier: 6.476A pdb=" N VAL E 20 " --> pdb=" O LEU E 41 " (cutoff:3.500A) removed outlier: 7.731A pdb=" N LEU E 40 " --> pdb=" O ASP E 65 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'E' and resid 71 through 72 removed outlier: 7.077A pdb=" N TRP E 71 " --> pdb=" O ALA E 98 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA6 Processing sheet with id=AA7, first strand: chain 'B' and resid 6 through 8 removed outlier: 6.827A pdb=" N LEU B 35 " --> pdb=" O HIS B 60 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'b' and resid 6 through 8 removed outlier: 3.522A pdb=" N VAL b 36 " --> pdb=" O VAL b 12 " (cutoff:3.500A) removed outlier: 7.071A pdb=" N LEU b 35 " --> pdb=" O HIS b 60 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'b' and resid 66 through 67 removed outlier: 7.041A pdb=" N TRP b 66 " --> pdb=" O VAL b 94 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA9 Processing sheet with id=AB1, first strand: chain 'X' and resid 9 through 13 removed outlier: 6.947A pdb=" N LEU X 18 " --> pdb=" O ARG X 38 " (cutoff:3.500A) removed outlier: 5.638A pdb=" N LEU X 40 " --> pdb=" O LEU X 18 " (cutoff:3.500A) removed outlier: 6.900A pdb=" N VAL X 20 " --> pdb=" O LEU X 40 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'X' and resid 71 through 72 removed outlier: 6.640A pdb=" N TRP X 71 " --> pdb=" O ALA X 98 " (cutoff:3.500A) No H-bonds generated for sheet with id=AB2 149 hydrogen bonds defined for protein. 351 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: 1.49 Time building geometry restraints manager: 2.11 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 2060 1.34 - 1.47: 1427 1.47 - 1.59: 2898 1.59 - 1.71: 1 1.71 - 1.83: 56 Bond restraints: 6442 Sorted by residual: bond pdb=" CB PRO c 29 " pdb=" CG PRO c 29 " ideal model delta sigma weight residual 1.492 1.651 -0.159 5.00e-02 4.00e+02 1.01e+01 bond pdb=" CB PRO C 2 " pdb=" CG PRO C 2 " ideal model delta sigma weight residual 1.492 1.572 -0.080 5.00e-02 4.00e+02 2.54e+00 bond pdb=" C ALA B 86 " pdb=" N PRO B 87 " ideal model delta sigma weight residual 1.334 1.347 -0.013 8.40e-03 1.42e+04 2.38e+00 bond pdb=" C1 BMA G 3 " pdb=" C2 BMA G 3 " ideal model delta sigma weight residual 1.519 1.548 -0.029 2.00e-02 2.50e+03 2.17e+00 bond pdb=" CG PRO C 2 " pdb=" CD PRO C 2 " ideal model delta sigma weight residual 1.503 1.551 -0.048 3.40e-02 8.65e+02 1.95e+00 ... (remaining 6437 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.61: 8639 2.61 - 5.22: 181 5.22 - 7.83: 18 7.83 - 10.45: 9 10.45 - 13.06: 1 Bond angle restraints: 8848 Sorted by residual: angle pdb=" CA PRO c 29 " pdb=" N PRO c 29 " pdb=" CD PRO c 29 " ideal model delta sigma weight residual 112.00 98.94 13.06 1.40e+00 5.10e-01 8.70e+01 angle pdb=" N GLY b 124 " pdb=" CA GLY b 124 " pdb=" C GLY b 124 " ideal model delta sigma weight residual 115.66 105.40 10.26 1.56e+00 4.11e-01 4.33e+01 angle pdb=" CA PRO C 2 " pdb=" N PRO C 2 " pdb=" CD PRO C 2 " ideal model delta sigma weight residual 112.00 106.47 5.53 1.40e+00 5.10e-01 1.56e+01 angle pdb=" N PRO c 29 " pdb=" CD PRO c 29 " pdb=" CG PRO c 29 " ideal model delta sigma weight residual 103.20 98.14 5.06 1.50e+00 4.44e-01 1.14e+01 angle pdb=" N PRO c 29 " pdb=" CA PRO c 29 " pdb=" C PRO c 29 " ideal model delta sigma weight residual 111.68 116.99 -5.31 1.67e+00 3.59e-01 1.01e+01 ... (remaining 8843 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 22.26: 3602 22.26 - 44.52: 419 44.52 - 66.78: 71 66.78 - 89.03: 20 89.03 - 111.29: 16 Dihedral angle restraints: 4128 sinusoidal: 1818 harmonic: 2310 Sorted by residual: dihedral pdb=" CB CYS X 73 " pdb=" SG CYS X 73 " pdb=" SG CYS X 97 " pdb=" CB CYS X 97 " ideal model delta sinusoidal sigma weight residual 93.00 14.03 78.97 1 1.00e+01 1.00e-02 7.76e+01 dihedral pdb=" CB CYS b 122 " pdb=" SG CYS b 122 " pdb=" SG CYS a 485 " pdb=" CB CYS a 485 " ideal model delta sinusoidal sigma weight residual 93.00 171.92 -78.92 1 1.00e+01 1.00e-02 7.76e+01 dihedral pdb=" CB CYS B 122 " pdb=" SG CYS B 122 " pdb=" SG CYS a 484 " pdb=" CB CYS a 484 " ideal model delta sinusoidal sigma weight residual 93.00 168.07 -75.07 1 1.00e+01 1.00e-02 7.13e+01 ... (remaining 4125 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.056: 870 0.056 - 0.111: 160 0.111 - 0.167: 22 0.167 - 0.223: 3 0.223 - 0.278: 5 Chirality restraints: 1060 Sorted by residual: chirality pdb=" C5 BMA F 3 " pdb=" C4 BMA F 3 " pdb=" C6 BMA F 3 " pdb=" O5 BMA F 3 " both_signs ideal model delta sigma weight residual False -2.27 -2.55 0.28 2.00e-01 2.50e+01 1.94e+00 chirality pdb=" C5 BMA G 3 " pdb=" C4 BMA G 3 " pdb=" C6 BMA G 3 " pdb=" O5 BMA G 3 " both_signs ideal model delta sigma weight residual False -2.27 -2.54 0.27 2.00e-01 2.50e+01 1.76e+00 chirality pdb=" CA PRO c 29 " pdb=" N PRO c 29 " pdb=" C PRO c 29 " pdb=" CB PRO c 29 " both_signs ideal model delta sigma weight residual False 2.72 2.47 0.25 2.00e-01 2.50e+01 1.51e+00 ... (remaining 1057 not shown) Planarity restraints: 1135 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C PHE c 28 " 0.089 5.00e-02 4.00e+02 1.24e-01 2.47e+01 pdb=" N PRO c 29 " -0.214 5.00e-02 4.00e+02 pdb=" CA PRO c 29 " 0.059 5.00e-02 4.00e+02 pdb=" CD PRO c 29 " 0.066 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C CYS C 1 " -0.059 5.00e-02 4.00e+02 8.78e-02 1.23e+01 pdb=" N PRO C 2 " 0.152 5.00e-02 4.00e+02 pdb=" CA PRO C 2 " -0.045 5.00e-02 4.00e+02 pdb=" CD PRO C 2 " -0.047 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C PRO X 51 " -0.053 5.00e-02 4.00e+02 8.04e-02 1.04e+01 pdb=" N PRO X 52 " 0.139 5.00e-02 4.00e+02 pdb=" CA PRO X 52 " -0.042 5.00e-02 4.00e+02 pdb=" CD PRO X 52 " -0.044 5.00e-02 4.00e+02 ... (remaining 1132 not shown) Histogram of nonbonded interaction distances: 2.26 - 2.79: 1511 2.79 - 3.32: 5736 3.32 - 3.84: 10220 3.84 - 4.37: 11281 4.37 - 4.90: 19207 Nonbonded interactions: 47955 Sorted by model distance: nonbonded pdb=" NH2 ARG c 102 " pdb=" OE2 GLU c 111 " model vdw 2.261 3.120 nonbonded pdb=" NH1 ARG c 67 " pdb=" OD2 ASP c 69 " model vdw 2.273 3.120 nonbonded pdb=" O SER X 126 " pdb=" OG SER X 126 " model vdw 2.289 3.040 nonbonded pdb=" NH1 ARG C 89 " pdb=" OD1 ASP E 86 " model vdw 2.319 3.120 nonbonded pdb=" O CYS b 122 " pdb=" N GLY b 124 " model vdw 2.330 3.120 ... (remaining 47950 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.00 Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'D' selection = chain 'H' selection = chain 'I' } ncs_group { reference = chain 'B' selection = chain 'b' } ncs_group { reference = (chain 'C' and (resid 1 through 113 or resid 115 through 117)) selection = (chain 'c' and (resid 1 through 113 or resid 115 through 117)) } ncs_group { reference = chain 'E' selection = (chain 'X' and resid 1 through 125) } ncs_group { reference = chain 'F' selection = chain 'G' } 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=0.50 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: 2.830 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.060 Construct map_model_manager: 0.010 Extract box with map and model: 0.270 Check model and map are aligned: 0.040 Set scattering table: 0.060 Process input model: 19.800 Find NCS groups from input model: 0.240 Set up NCS constraints: 0.030 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.000 Load rotamer database and sin/cos tables:2.640 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 25.980 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7303 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.005 0.159 6442 Z= 0.293 Angle : 0.866 13.057 8848 Z= 0.414 Chirality : 0.048 0.278 1060 Planarity : 0.008 0.124 1129 Dihedral : 19.138 111.292 2556 Min Nonbonded Distance : 2.261 Molprobity Statistics. All-atom Clashscore : 13.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.79 % Favored : 94.21 % Rotamer: Outliers : 0.94 % Allowed : 32.97 % Favored : 66.09 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.77 (0.29), residues: 797 helix: -1.26 (0.39), residues: 161 sheet: -1.45 (0.69), residues: 68 loop : -1.15 (0.26), residues: 568 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.037 0.001 TRP b 123 HIS 0.002 0.001 HIS E 72 PHE 0.008 0.001 PHE c 28 TYR 0.032 0.002 TYR a 492 ARG 0.008 0.001 ARG b 89 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 137 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 131 time to evaluate : 0.652 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: b 113 ARG cc_start: 0.7548 (mtp85) cc_final: 0.7257 (mtp85) outliers start: 6 outliers final: 2 residues processed: 134 average time/residue: 0.3884 time to fit residues: 59.7489 Evaluate side-chains 132 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 130 time to evaluate : 0.726 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain C residue 12 VAL Chi-restraints excluded: chain a residue 485 CYS Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 65 optimal weight: 0.5980 chunk 59 optimal weight: 0.9990 chunk 32 optimal weight: 1.9990 chunk 20 optimal weight: 0.8980 chunk 39 optimal weight: 0.8980 chunk 31 optimal weight: 0.7980 chunk 61 optimal weight: 0.8980 chunk 23 optimal weight: 0.2980 chunk 37 optimal weight: 0.5980 chunk 45 optimal weight: 0.0020 chunk 70 optimal weight: 0.8980 overall best weight: 0.4588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: c 64 ASN ** C 40 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** E 115 GLN B 64 ASN b 40 ASN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3778 r_free = 0.3778 target = 0.137105 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3557 r_free = 0.3557 target = 0.120868 restraints weight = 10031.314| |-----------------------------------------------------------------------------| r_work (start): 0.3561 rms_B_bonded: 2.50 r_work: 0.3449 rms_B_bonded: 2.75 restraints_weight: 0.5000 r_work: 0.3298 rms_B_bonded: 4.90 restraints_weight: 0.2500 r_work (final): 0.3298 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3299 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3299 r_free = 0.3299 target_work(ls_wunit_k1) = 0.104 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3299 r_free = 0.3299 target_work(ls_wunit_k1) = 0.104 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3299 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7632 moved from start: 0.0786 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.046 6442 Z= 0.278 Angle : 0.802 9.761 8848 Z= 0.367 Chirality : 0.046 0.179 1060 Planarity : 0.006 0.064 1129 Dihedral : 12.657 87.122 1177 Min Nonbonded Distance : 2.431 Molprobity Statistics. All-atom Clashscore : 10.68 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.93 % Favored : 93.07 % Rotamer: Outliers : 3.61 % Allowed : 30.14 % Favored : 66.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.39 (0.30), residues: 797 helix: -0.25 (0.41), residues: 165 sheet: -1.48 (0.68), residues: 68 loop : -1.10 (0.26), residues: 564 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.029 0.002 TRP b 123 HIS 0.003 0.001 HIS c 60 PHE 0.005 0.001 PHE X 55 TYR 0.035 0.002 TYR a 492 ARG 0.010 0.001 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 23 poor density : 130 time to evaluate : 0.545 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: E 21 ASP cc_start: 0.7440 (OUTLIER) cc_final: 0.7079 (t70) REVERT: E 36 ARG cc_start: 0.8485 (OUTLIER) cc_final: 0.7812 (ptt90) REVERT: b 76 ARG cc_start: 0.8334 (tpt90) cc_final: 0.8112 (tmt170) REVERT: b 113 ARG cc_start: 0.7859 (mtp85) cc_final: 0.7620 (mtp85) REVERT: b 138 LEU cc_start: -0.2571 (OUTLIER) cc_final: -0.3131 (tp) REVERT: X 62 GLN cc_start: 0.8557 (OUTLIER) cc_final: 0.8162 (pt0) outliers start: 23 outliers final: 7 residues processed: 147 average time/residue: 0.4027 time to fit residues: 68.0450 Evaluate side-chains 138 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 127 time to evaluate : 0.612 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain C residue 26 THR Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 36 ARG Chi-restraints excluded: chain b residue 138 LEU Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 62 GLN Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 116 LEU Chi-restraints excluded: chain X residue 146 LEU Chi-restraints excluded: chain a residue 485 CYS Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 46 optimal weight: 0.8980 chunk 31 optimal weight: 0.0010 chunk 4 optimal weight: 0.4980 chunk 64 optimal weight: 0.4980 chunk 15 optimal weight: 1.9990 chunk 44 optimal weight: 0.9990 chunk 16 optimal weight: 0.8980 chunk 0 optimal weight: 1.9990 chunk 42 optimal weight: 0.8980 chunk 20 optimal weight: 0.8980 chunk 3 optimal weight: 0.5980 overall best weight: 0.4986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** c 40 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 40 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 115 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3771 r_free = 0.3771 target = 0.136534 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.3545 r_free = 0.3545 target = 0.120169 restraints weight = 10175.741| |-----------------------------------------------------------------------------| r_work (start): 0.3548 rms_B_bonded: 2.52 r_work: 0.3437 rms_B_bonded: 2.75 restraints_weight: 0.5000 r_work: 0.3285 rms_B_bonded: 4.92 restraints_weight: 0.2500 r_work (final): 0.3285 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3288 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3288 r_free = 0.3288 target_work(ls_wunit_k1) = 0.103 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3288 r_free = 0.3288 target_work(ls_wunit_k1) = 0.103 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3288 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7643 moved from start: 0.1017 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.044 6442 Z= 0.283 Angle : 0.776 9.620 8848 Z= 0.354 Chirality : 0.046 0.194 1060 Planarity : 0.006 0.057 1129 Dihedral : 10.395 76.162 1176 Min Nonbonded Distance : 2.428 Molprobity Statistics. All-atom Clashscore : 9.80 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.05 % Favored : 92.95 % Rotamer: Outliers : 6.28 % Allowed : 28.10 % Favored : 65.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.29 (0.30), residues: 797 helix: -0.02 (0.41), residues: 165 sheet: -1.57 (0.67), residues: 68 loop : -1.05 (0.26), residues: 564 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.029 0.001 TRP b 123 HIS 0.002 0.001 HIS B 60 PHE 0.005 0.001 PHE a 491 TYR 0.028 0.002 TYR a 492 ARG 0.011 0.001 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 170 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 40 poor density : 130 time to evaluate : 0.663 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: C 72 LEU cc_start: 0.8835 (tp) cc_final: 0.8544 (tp) REVERT: E 21 ASP cc_start: 0.7472 (OUTLIER) cc_final: 0.7072 (t70) REVERT: E 36 ARG cc_start: 0.8525 (OUTLIER) cc_final: 0.7852 (ptt90) REVERT: E 95 VAL cc_start: 0.8169 (OUTLIER) cc_final: 0.7858 (t) REVERT: b 76 ARG cc_start: 0.8337 (tpt90) cc_final: 0.8131 (tmt170) REVERT: b 138 LEU cc_start: -0.2893 (OUTLIER) cc_final: -0.3434 (tp) REVERT: X 124 GLN cc_start: 0.7827 (OUTLIER) cc_final: 0.6955 (pp30) outliers start: 40 outliers final: 20 residues processed: 158 average time/residue: 0.3590 time to fit residues: 66.3251 Evaluate side-chains 148 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 25 poor density : 123 time to evaluate : 0.657 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 6 SER Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain c residue 91 LEU Chi-restraints excluded: chain C residue 7 CYS Chi-restraints excluded: chain C residue 26 THR Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 36 ARG Chi-restraints excluded: chain E residue 41 LEU Chi-restraints excluded: chain E residue 61 LEU Chi-restraints excluded: chain E residue 62 GLN Chi-restraints excluded: chain E residue 95 VAL Chi-restraints excluded: chain E residue 116 LEU Chi-restraints excluded: chain E residue 123 LEU Chi-restraints excluded: chain B residue 11 LEU Chi-restraints excluded: chain B residue 71 ARG Chi-restraints excluded: chain B residue 103 LEU Chi-restraints excluded: chain b residue 138 LEU Chi-restraints excluded: chain X residue 4 CYS Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 16 MET Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 116 LEU Chi-restraints excluded: chain X residue 124 GLN Chi-restraints excluded: chain X residue 146 LEU Chi-restraints excluded: chain a residue 485 CYS Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 61 optimal weight: 0.3980 chunk 64 optimal weight: 0.9990 chunk 7 optimal weight: 0.9980 chunk 72 optimal weight: 0.2980 chunk 25 optimal weight: 2.9990 chunk 52 optimal weight: 1.9990 chunk 34 optimal weight: 3.9990 chunk 11 optimal weight: 1.9990 chunk 77 optimal weight: 0.9980 chunk 46 optimal weight: 0.0870 chunk 23 optimal weight: 0.4980 overall best weight: 0.4558 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 40 ASN ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 115 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3770 r_free = 0.3770 target = 0.136583 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 37)----------------| | r_work = 0.3545 r_free = 0.3545 target = 0.120249 restraints weight = 10187.059| |-----------------------------------------------------------------------------| r_work (start): 0.3551 rms_B_bonded: 2.54 r_work: 0.3438 rms_B_bonded: 2.77 restraints_weight: 0.5000 r_work: 0.3286 rms_B_bonded: 4.93 restraints_weight: 0.2500 r_work (final): 0.3286 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3289 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3289 r_free = 0.3289 target_work(ls_wunit_k1) = 0.103 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3289 r_free = 0.3289 target_work(ls_wunit_k1) = 0.103 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3289 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7624 moved from start: 0.1146 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.043 6442 Z= 0.266 Angle : 0.767 9.494 8848 Z= 0.350 Chirality : 0.045 0.214 1060 Planarity : 0.006 0.053 1129 Dihedral : 9.301 67.696 1176 Min Nonbonded Distance : 2.429 Molprobity Statistics. All-atom Clashscore : 10.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.80 % Favored : 93.20 % Rotamer: Outliers : 5.18 % Allowed : 29.67 % Favored : 65.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.30 (0.30), residues: 797 helix: -0.07 (0.41), residues: 165 sheet: -1.58 (0.67), residues: 68 loop : -1.05 (0.26), residues: 564 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.021 0.001 TRP b 123 HIS 0.002 0.001 HIS E 72 PHE 0.004 0.001 PHE c 28 TYR 0.027 0.002 TYR a 492 ARG 0.008 0.001 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 163 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 33 poor density : 130 time to evaluate : 0.719 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: E 21 ASP cc_start: 0.7473 (OUTLIER) cc_final: 0.7060 (t70) REVERT: E 36 ARG cc_start: 0.8527 (OUTLIER) cc_final: 0.7851 (ptt90) REVERT: E 95 VAL cc_start: 0.8135 (OUTLIER) cc_final: 0.7821 (t) REVERT: b 21 TRP cc_start: 0.8396 (OUTLIER) cc_final: 0.7934 (t-100) REVERT: b 113 ARG cc_start: 0.7745 (mtp85) cc_final: 0.7476 (mtp85) REVERT: b 138 LEU cc_start: -0.2940 (OUTLIER) cc_final: -0.3483 (tp) REVERT: X 124 GLN cc_start: 0.7780 (OUTLIER) cc_final: 0.6868 (pp30) outliers start: 33 outliers final: 21 residues processed: 154 average time/residue: 0.3621 time to fit residues: 64.4921 Evaluate side-chains 151 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 27 poor density : 124 time to evaluate : 0.711 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 6 SER Chi-restraints excluded: chain c residue 26 THR Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain c residue 91 LEU Chi-restraints excluded: chain C residue 5 CYS Chi-restraints excluded: chain C residue 26 THR Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 36 ARG Chi-restraints excluded: chain E residue 41 LEU Chi-restraints excluded: chain E residue 61 LEU Chi-restraints excluded: chain E residue 62 GLN Chi-restraints excluded: chain E residue 95 VAL Chi-restraints excluded: chain E residue 102 LEU Chi-restraints excluded: chain E residue 116 LEU Chi-restraints excluded: chain E residue 123 LEU Chi-restraints excluded: chain B residue 71 ARG Chi-restraints excluded: chain B residue 103 LEU Chi-restraints excluded: chain b residue 21 TRP Chi-restraints excluded: chain b residue 138 LEU Chi-restraints excluded: chain X residue 4 CYS Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 16 MET Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 116 LEU Chi-restraints excluded: chain X residue 124 GLN Chi-restraints excluded: chain X residue 146 LEU Chi-restraints excluded: chain a residue 485 CYS Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 40 optimal weight: 0.7980 chunk 23 optimal weight: 0.7980 chunk 20 optimal weight: 0.5980 chunk 75 optimal weight: 0.5980 chunk 72 optimal weight: 0.9990 chunk 54 optimal weight: 0.0970 chunk 67 optimal weight: 0.6980 chunk 27 optimal weight: 0.0470 chunk 45 optimal weight: 2.9990 chunk 37 optimal weight: 0.9980 chunk 57 optimal weight: 0.3980 overall best weight: 0.3476 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 115 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3785 r_free = 0.3785 target = 0.137638 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.3559 r_free = 0.3559 target = 0.121172 restraints weight = 10168.137| |-----------------------------------------------------------------------------| r_work (start): 0.3561 rms_B_bonded: 2.55 r_work: 0.3450 rms_B_bonded: 2.78 restraints_weight: 0.5000 r_work: 0.3299 rms_B_bonded: 4.95 restraints_weight: 0.2500 r_work (final): 0.3299 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3303 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3303 r_free = 0.3303 target_work(ls_wunit_k1) = 0.104 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3303 r_free = 0.3303 target_work(ls_wunit_k1) = 0.104 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3303 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7610 moved from start: 0.1297 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 6442 Z= 0.226 Angle : 0.756 9.605 8848 Z= 0.343 Chirality : 0.044 0.221 1060 Planarity : 0.006 0.050 1129 Dihedral : 8.038 59.127 1176 Min Nonbonded Distance : 2.392 Molprobity Statistics. All-atom Clashscore : 9.48 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.17 % Favored : 93.83 % Rotamer: Outliers : 5.97 % Allowed : 29.04 % Favored : 64.99 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.28 (0.29), residues: 797 helix: 0.14 (0.42), residues: 159 sheet: -1.55 (0.67), residues: 68 loop : -1.10 (0.26), residues: 570 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.022 0.001 TRP b 123 HIS 0.002 0.000 HIS B 60 PHE 0.006 0.001 PHE c 28 TYR 0.025 0.002 TYR a 492 ARG 0.011 0.001 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 164 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 38 poor density : 126 time to evaluate : 0.789 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: E 21 ASP cc_start: 0.7462 (OUTLIER) cc_final: 0.7075 (t70) REVERT: E 95 VAL cc_start: 0.8091 (OUTLIER) cc_final: 0.7808 (t) REVERT: B 17 ARG cc_start: 0.7931 (mtm-85) cc_final: 0.7618 (mtm180) REVERT: B 72 LEU cc_start: 0.8748 (tt) cc_final: 0.8461 (tp) REVERT: b 21 TRP cc_start: 0.8394 (OUTLIER) cc_final: 0.7929 (t-100) REVERT: b 113 ARG cc_start: 0.7670 (mtp85) cc_final: 0.7335 (mtp85) REVERT: X 11 ARG cc_start: 0.7988 (OUTLIER) cc_final: 0.7759 (mtp180) REVERT: X 62 GLN cc_start: 0.8558 (OUTLIER) cc_final: 0.8266 (pt0) REVERT: X 124 GLN cc_start: 0.7768 (OUTLIER) cc_final: 0.6903 (pp30) outliers start: 38 outliers final: 23 residues processed: 154 average time/residue: 0.3909 time to fit residues: 69.9235 Evaluate side-chains 151 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 29 poor density : 122 time to evaluate : 0.686 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 6 SER Chi-restraints excluded: chain c residue 26 THR Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain c residue 91 LEU Chi-restraints excluded: chain C residue 5 CYS Chi-restraints excluded: chain C residue 7 CYS Chi-restraints excluded: chain C residue 26 THR Chi-restraints excluded: chain C residue 36 VAL Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 41 LEU Chi-restraints excluded: chain E residue 61 LEU Chi-restraints excluded: chain E residue 62 GLN Chi-restraints excluded: chain E residue 95 VAL Chi-restraints excluded: chain E residue 102 LEU Chi-restraints excluded: chain E residue 116 LEU Chi-restraints excluded: chain E residue 123 LEU Chi-restraints excluded: chain B residue 71 ARG Chi-restraints excluded: chain B residue 103 LEU Chi-restraints excluded: chain b residue 21 TRP Chi-restraints excluded: chain X residue 4 CYS Chi-restraints excluded: chain X residue 11 ARG Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 16 MET Chi-restraints excluded: chain X residue 20 VAL Chi-restraints excluded: chain X residue 62 GLN Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 116 LEU Chi-restraints excluded: chain X residue 124 GLN Chi-restraints excluded: chain X residue 146 LEU Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 18 optimal weight: 0.2980 chunk 40 optimal weight: 1.9990 chunk 30 optimal weight: 0.9990 chunk 58 optimal weight: 0.9980 chunk 52 optimal weight: 0.4980 chunk 73 optimal weight: 0.8980 chunk 15 optimal weight: 1.9990 chunk 17 optimal weight: 0.0370 chunk 47 optimal weight: 0.9980 chunk 75 optimal weight: 0.5980 chunk 4 optimal weight: 0.0470 overall best weight: 0.2956 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 115 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3791 r_free = 0.3791 target = 0.138250 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.3569 r_free = 0.3569 target = 0.121959 restraints weight = 10130.497| |-----------------------------------------------------------------------------| r_work (start): 0.3573 rms_B_bonded: 2.52 r_work: 0.3463 rms_B_bonded: 2.74 restraints_weight: 0.5000 r_work: 0.3314 rms_B_bonded: 4.90 restraints_weight: 0.2500 r_work (final): 0.3314 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3315 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3315 r_free = 0.3315 target_work(ls_wunit_k1) = 0.105 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3315 r_free = 0.3315 target_work(ls_wunit_k1) = 0.105 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3315 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7615 moved from start: 0.1404 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 6442 Z= 0.211 Angle : 0.740 9.160 8848 Z= 0.338 Chirality : 0.043 0.244 1060 Planarity : 0.006 0.048 1129 Dihedral : 7.089 58.450 1175 Min Nonbonded Distance : 2.406 Molprobity Statistics. All-atom Clashscore : 9.48 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.42 % Favored : 93.58 % Rotamer: Outliers : 5.97 % Allowed : 28.41 % Favored : 65.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.23 (0.30), residues: 797 helix: 0.23 (0.42), residues: 159 sheet: -1.59 (0.67), residues: 68 loop : -1.07 (0.26), residues: 570 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.018 0.001 TRP b 123 HIS 0.002 0.000 HIS B 139 PHE 0.004 0.001 PHE c 28 TYR 0.024 0.001 TYR a 492 ARG 0.010 0.000 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 162 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 38 poor density : 124 time to evaluate : 0.639 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: E 21 ASP cc_start: 0.7380 (OUTLIER) cc_final: 0.6961 (t70) REVERT: E 95 VAL cc_start: 0.8128 (OUTLIER) cc_final: 0.7867 (t) REVERT: B 17 ARG cc_start: 0.7854 (mtm-85) cc_final: 0.7645 (mtm-85) REVERT: B 72 LEU cc_start: 0.8746 (tt) cc_final: 0.8467 (tp) REVERT: b 21 TRP cc_start: 0.8376 (OUTLIER) cc_final: 0.7938 (t-100) REVERT: X 11 ARG cc_start: 0.7965 (OUTLIER) cc_final: 0.7739 (mtp180) REVERT: X 124 GLN cc_start: 0.7721 (OUTLIER) cc_final: 0.6831 (pp30) outliers start: 38 outliers final: 23 residues processed: 152 average time/residue: 0.3554 time to fit residues: 62.6952 Evaluate side-chains 149 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 28 poor density : 121 time to evaluate : 0.637 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 6 SER Chi-restraints excluded: chain c residue 26 THR Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain c residue 91 LEU Chi-restraints excluded: chain C residue 5 CYS Chi-restraints excluded: chain C residue 7 CYS Chi-restraints excluded: chain C residue 26 THR Chi-restraints excluded: chain C residue 36 VAL Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 61 LEU Chi-restraints excluded: chain E residue 62 GLN Chi-restraints excluded: chain E residue 95 VAL Chi-restraints excluded: chain E residue 102 LEU Chi-restraints excluded: chain E residue 116 LEU Chi-restraints excluded: chain E residue 123 LEU Chi-restraints excluded: chain B residue 71 ARG Chi-restraints excluded: chain B residue 103 LEU Chi-restraints excluded: chain b residue 21 TRP Chi-restraints excluded: chain X residue 4 CYS Chi-restraints excluded: chain X residue 11 ARG Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 18 LEU Chi-restraints excluded: chain X residue 20 VAL Chi-restraints excluded: chain X residue 110 ARG Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 124 GLN Chi-restraints excluded: chain X residue 146 LEU Chi-restraints excluded: chain a residue 485 CYS Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 25 optimal weight: 0.9980 chunk 55 optimal weight: 0.4980 chunk 60 optimal weight: 0.2980 chunk 58 optimal weight: 0.9990 chunk 38 optimal weight: 0.0270 chunk 40 optimal weight: 0.5980 chunk 67 optimal weight: 0.9980 chunk 37 optimal weight: 0.9990 chunk 4 optimal weight: 0.0050 chunk 65 optimal weight: 0.5980 chunk 72 optimal weight: 0.9990 overall best weight: 0.2852 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 115 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3796 r_free = 0.3796 target = 0.138734 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 41)----------------| | r_work = 0.3584 r_free = 0.3584 target = 0.122806 restraints weight = 10081.111| |-----------------------------------------------------------------------------| r_work (start): 0.3579 rms_B_bonded: 2.52 r_work: 0.3468 rms_B_bonded: 2.75 restraints_weight: 0.5000 r_work: 0.3317 rms_B_bonded: 4.93 restraints_weight: 0.2500 r_work (final): 0.3317 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3321 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3321 r_free = 0.3321 target_work(ls_wunit_k1) = 0.105 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3321 r_free = 0.3321 target_work(ls_wunit_k1) = 0.105 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3321 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7593 moved from start: 0.1480 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 6442 Z= 0.206 Angle : 0.742 9.082 8848 Z= 0.336 Chirality : 0.044 0.294 1060 Planarity : 0.006 0.048 1129 Dihedral : 6.620 57.475 1175 Min Nonbonded Distance : 2.390 Molprobity Statistics. All-atom Clashscore : 9.96 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.17 % Favored : 93.83 % Rotamer: Outliers : 5.49 % Allowed : 29.36 % Favored : 65.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.24 (0.30), residues: 797 helix: 0.18 (0.42), residues: 159 sheet: -1.51 (0.67), residues: 68 loop : -1.08 (0.26), residues: 570 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.043 0.002 TRP a 498 HIS 0.002 0.001 HIS B 139 PHE 0.004 0.001 PHE c 28 TYR 0.023 0.001 TYR a 492 ARG 0.010 0.000 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 161 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 35 poor density : 126 time to evaluate : 1.165 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: C 72 LEU cc_start: 0.8813 (tp) cc_final: 0.8402 (tp) REVERT: E 21 ASP cc_start: 0.7331 (OUTLIER) cc_final: 0.6910 (t70) REVERT: E 95 VAL cc_start: 0.8081 (OUTLIER) cc_final: 0.7841 (t) REVERT: B 17 ARG cc_start: 0.7911 (mtm-85) cc_final: 0.7563 (mtm180) REVERT: B 72 LEU cc_start: 0.8738 (OUTLIER) cc_final: 0.8454 (tp) REVERT: b 21 TRP cc_start: 0.8371 (OUTLIER) cc_final: 0.7909 (t-100) REVERT: X 62 GLN cc_start: 0.8528 (OUTLIER) cc_final: 0.8231 (pt0) REVERT: X 124 GLN cc_start: 0.7762 (OUTLIER) cc_final: 0.6959 (pp30) outliers start: 35 outliers final: 24 residues processed: 153 average time/residue: 0.3948 time to fit residues: 70.4075 Evaluate side-chains 153 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 30 poor density : 123 time to evaluate : 0.681 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 6 SER Chi-restraints excluded: chain c residue 26 THR Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain c residue 91 LEU Chi-restraints excluded: chain C residue 5 CYS Chi-restraints excluded: chain C residue 7 CYS Chi-restraints excluded: chain C residue 36 VAL Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 61 LEU Chi-restraints excluded: chain E residue 62 GLN Chi-restraints excluded: chain E residue 95 VAL Chi-restraints excluded: chain E residue 102 LEU Chi-restraints excluded: chain E residue 116 LEU Chi-restraints excluded: chain E residue 123 LEU Chi-restraints excluded: chain B residue 71 ARG Chi-restraints excluded: chain B residue 72 LEU Chi-restraints excluded: chain B residue 103 LEU Chi-restraints excluded: chain b residue 21 TRP Chi-restraints excluded: chain X residue 4 CYS Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 16 MET Chi-restraints excluded: chain X residue 18 LEU Chi-restraints excluded: chain X residue 20 VAL Chi-restraints excluded: chain X residue 47 LEU Chi-restraints excluded: chain X residue 62 GLN Chi-restraints excluded: chain X residue 110 ARG Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 124 GLN Chi-restraints excluded: chain X residue 146 LEU Chi-restraints excluded: chain a residue 485 CYS Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 68 optimal weight: 1.9990 chunk 21 optimal weight: 0.7980 chunk 44 optimal weight: 0.2980 chunk 57 optimal weight: 0.6980 chunk 29 optimal weight: 0.4980 chunk 35 optimal weight: 0.5980 chunk 39 optimal weight: 1.9990 chunk 8 optimal weight: 0.5980 chunk 11 optimal weight: 1.9990 chunk 76 optimal weight: 0.0670 chunk 15 optimal weight: 0.9980 overall best weight: 0.4118 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 64 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 115 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** X 115 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3785 r_free = 0.3785 target = 0.137875 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.3578 r_free = 0.3578 target = 0.121986 restraints weight = 10179.808| |-----------------------------------------------------------------------------| r_work (start): 0.3573 rms_B_bonded: 2.52 r_work: 0.3459 rms_B_bonded: 2.76 restraints_weight: 0.5000 r_work: 0.3307 rms_B_bonded: 4.95 restraints_weight: 0.2500 r_work (final): 0.3307 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3309 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3309 r_free = 0.3309 target_work(ls_wunit_k1) = 0.104 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3309 r_free = 0.3309 target_work(ls_wunit_k1) = 0.104 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3309 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7616 moved from start: 0.1542 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.040 6442 Z= 0.248 Angle : 0.752 9.152 8848 Z= 0.341 Chirality : 0.045 0.295 1060 Planarity : 0.006 0.048 1129 Dihedral : 6.340 54.811 1175 Min Nonbonded Distance : 2.408 Molprobity Statistics. All-atom Clashscore : 10.04 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.93 % Favored : 93.07 % Rotamer: Outliers : 6.28 % Allowed : 28.26 % Favored : 65.46 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 7.06 % Cis-general : 0.00 % Twisted Proline : 1.18 % Twisted General : 0.14 % 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: -1.29 (0.30), residues: 797 helix: 0.13 (0.42), residues: 159 sheet: -1.47 (0.67), residues: 68 loop : -1.12 (0.26), residues: 570 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.034 0.002 TRP a 498 HIS 0.002 0.001 HIS B 60 PHE 0.005 0.001 PHE a 491 TYR 0.023 0.002 TYR a 492 ARG 0.011 0.001 ARG C 89 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1594 Ramachandran restraints generated. 797 Oldfield, 0 Emsley, 797 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 168 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 40 poor density : 128 time to evaluate : 0.687 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: c 113 ARG cc_start: 0.8405 (tpp-160) cc_final: 0.8153 (tpp80) REVERT: C 72 LEU cc_start: 0.8828 (tp) cc_final: 0.8413 (tp) REVERT: E 21 ASP cc_start: 0.7306 (OUTLIER) cc_final: 0.6888 (t70) REVERT: E 95 VAL cc_start: 0.8110 (OUTLIER) cc_final: 0.7826 (t) REVERT: B 72 LEU cc_start: 0.8742 (OUTLIER) cc_final: 0.8464 (tp) REVERT: b 21 TRP cc_start: 0.8397 (OUTLIER) cc_final: 0.7872 (t-100) REVERT: X 11 ARG cc_start: 0.7924 (OUTLIER) cc_final: 0.7717 (mtp180) REVERT: X 62 GLN cc_start: 0.8547 (OUTLIER) cc_final: 0.8249 (pt0) REVERT: X 124 GLN cc_start: 0.7790 (OUTLIER) cc_final: 0.7018 (pp30) outliers start: 40 outliers final: 29 residues processed: 155 average time/residue: 0.3311 time to fit residues: 59.8701 Evaluate side-chains 163 residues out of total 637 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 36 poor density : 127 time to evaluate : 0.674 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain c residue 6 SER Chi-restraints excluded: chain c residue 26 THR Chi-restraints excluded: chain c residue 35 LEU Chi-restraints excluded: chain c residue 91 LEU Chi-restraints excluded: chain C residue 5 CYS Chi-restraints excluded: chain C residue 7 CYS Chi-restraints excluded: chain C residue 26 THR Chi-restraints excluded: chain C residue 36 VAL Chi-restraints excluded: chain E residue 21 ASP Chi-restraints excluded: chain E residue 61 LEU Chi-restraints excluded: chain E residue 62 GLN Chi-restraints excluded: chain E residue 95 VAL Chi-restraints excluded: chain E residue 102 LEU Chi-restraints excluded: chain E residue 116 LEU Chi-restraints excluded: chain E residue 123 LEU Chi-restraints excluded: chain B residue 35 LEU Chi-restraints excluded: chain B residue 71 ARG Chi-restraints excluded: chain B residue 72 LEU Chi-restraints excluded: chain B residue 103 LEU Chi-restraints excluded: chain b residue 21 TRP Chi-restraints excluded: chain X residue 4 CYS Chi-restraints excluded: chain X residue 11 ARG Chi-restraints excluded: chain X residue 14 GLU Chi-restraints excluded: chain X residue 15 THR Chi-restraints excluded: chain X residue 16 MET Chi-restraints excluded: chain X residue 18 LEU Chi-restraints excluded: chain X residue 20 VAL Chi-restraints excluded: chain X residue 47 LEU Chi-restraints excluded: chain X residue 62 GLN Chi-restraints excluded: chain X residue 110 ARG Chi-restraints excluded: chain X residue 112 THR Chi-restraints excluded: chain X residue 115 GLN Chi-restraints excluded: chain X residue 124 GLN Chi-restraints excluded: chain X residue 146 LEU Chi-restraints excluded: chain a residue 485 CYS Chi-restraints excluded: chain a residue 501 PHE Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/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-5624/modules/cctbx_project/iotbx/cli_parser.py", line 980, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/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-5624/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/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-5624/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-5624/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-5624/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-5624/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-5624/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-5624/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-5624/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-5624/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/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-5624/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-5624/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-5624/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-5624/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: 60.5756 > 50: distance: 36 - 86: 9.199 distance: 50 - 54: 8.156 distance: 54 - 55: 12.924 distance: 55 - 56: 4.063 distance: 55 - 58: 8.053 distance: 56 - 57: 8.289 distance: 56 - 62: 7.051 distance: 58 - 59: 15.810 distance: 59 - 60: 7.451 distance: 59 - 61: 8.479 distance: 62 - 63: 3.301 distance: 63 - 64: 9.245 distance: 63 - 66: 4.176 distance: 64 - 65: 9.031 distance: 64 - 70: 11.504 distance: 66 - 67: 3.148 distance: 67 - 68: 5.174 distance: 67 - 69: 12.066 distance: 70 - 71: 7.279 distance: 71 - 72: 7.801 distance: 71 - 74: 11.695 distance: 72 - 73: 17.085 distance: 72 - 81: 5.216 distance: 74 - 75: 7.713 distance: 75 - 76: 12.230 distance: 76 - 77: 12.733 distance: 77 - 78: 7.912 distance: 78 - 79: 12.848 distance: 78 - 80: 12.007 distance: 81 - 82: 8.107 distance: 82 - 83: 7.694 distance: 83 - 84: 15.606 distance: 83 - 87: 16.779 distance: 85 - 86: 11.653 distance: 87 - 88: 15.212 distance: 88 - 89: 10.247 distance: 88 - 91: 12.688 distance: 89 - 90: 5.091 distance: 89 - 94: 6.082 distance: 91 - 92: 20.249 distance: 91 - 93: 24.142 distance: 94 - 95: 32.261 distance: 95 - 96: 41.593 distance: 95 - 98: 21.538 distance: 96 - 99: 21.975 distance: 99 - 100: 40.628 distance: 99 - 105: 39.083 distance: 100 - 101: 19.526 distance: 100 - 103: 18.029 distance: 101 - 102: 31.624 distance: 101 - 106: 14.663 distance: 102 - 118: 30.812 distance: 104 - 105: 4.389 distance: 106 - 107: 8.875 distance: 106 - 112: 15.189 distance: 107 - 108: 6.877 distance: 107 - 110: 3.395 distance: 108 - 109: 8.662 distance: 108 - 113: 15.157 distance: 110 - 111: 7.924 distance: 111 - 112: 11.424 distance: 113 - 114: 31.079 distance: 114 - 115: 4.727 distance: 114 - 117: 23.016 distance: 115 - 116: 19.261 distance: 115 - 118: 15.256 distance: 118 - 119: 10.614 distance: 119 - 120: 9.982 distance: 120 - 121: 4.248 distance: 120 - 126: 21.593 distance: 122 - 123: 13.711 distance: 123 - 124: 23.826 distance: 123 - 125: 28.726