- phenixbb - phenix-online.org

refinement problem
by Sneha Rangarajan 01 Oct '14

01 Oct '14

Hello, I have a 1.9A dataset for a single chain Fv-peptide complex that shows psedomerohedral twin operators in Xtriage (1 in p212121, 3 in p21 and 7 in p1). Although most programs (in CCP4) identify the best SG to be p212121, MR cannot find a solution with p212121 due to packing clashes. I tried reprocessing it in p21 and MR can find only one molecule in the asu. The map looks quite good with good density for the sidechains. However the Rfactors wont come down. After a few rounds of refinement (rigidbody, Xray, TLS and Ind B fac) it is still 38/48. I tried processing it in p1 and have managed to find 3 molecules in the asu (missing one molecule). The map looks good. While the Rfactors are lower, the gap between Rwork and Rfree is still high (27/40). I don't know what space group I should choose to ultimately refine it and get the Rfactors low. If you need any more information, I would be happy to provide it. Does anyone have any suggestions? Thanks a lot, Sneha

2 2

Target function for refinement in case of Pseudo-translational symmetry
by Varun Bhaskar 30 Sep '14

30 Sep '14

Hi, I collected a dataset for a binary complex (72kDa size) to 3.62A resolution which could be processed in P32 2 1 space group. Data quality analysis with phenix.xtriage showed a presence of Pseudotranslational symmetry (Non-origin peak in the patterson with 60% intensity at 0.33, -0.33, -0.013). I determined the structure with MR. I have 6 copies in ASU. I tried Zanuda for SG validation and this looks like the right space group. For refinement when I use Maximum likelihood (ML) as the target function, my Rfactor/Rfree is around 31.3/32 but when I change it to Least square (LS) my Rfactor/Rfree drops to 27.3/29.8. The map quality after ML or LS refinement looks very similar for most of the molecule, except for few places where the map from ML has slightly better features. I read people use LS instead of ML only in the case of twinning but I am not sure if it could also be used in the case of PST. Would the use of LS as the target function for refinement be right in this case? Thanks in advance Best Regards Varun

2 1

Histidine as a hydrogen bond acceptor in Phenix refine
by Phan, Jason 27 Sep '14

27 Sep '14

Hi All, I've been refining a 1.1 Å structure with explicit hydrogen atoms and noticed that the serine HG atom (hydrogen atom attached to oxygen gamma) is rotated more than 90 degrees from the OG - NE2 (histidine) line (5 cases observed). In this position, the HG is staggered with the two carbon beta hydrogens and should be the preferred configuration. However, there is a histidine nearby with atom NE2 at 2.7 Å from the serine OG and in good geometry for Hbond formation. The other nitrogen, ND1, is protonated and interacting with an aspartic acid (similar to the catalytic triad configuration). In this protomer, the NE2 atom has a lone pair and should act as an acceptor. Thus, the OG - H - NE2 angle should be as close to 180 as possible for an ideal hydrogen bond angle. The OG - H - NE2 angle of more than 90 degrees looks a little off to me. I checked the other program, e.g. Schrodinger Prime to see how it treats this situation and found that all the OG - H - NE2 angles are within 10-60 degrees from the ideal value of 180, certainly not 90 as observed in Phenix refine. My question is is NE2 or ND1 atom of the HID or HIE protomer treated as a potential hydrogen bond acceptor in Phenix? For other acceptors like the backbone carbonyl oxygen, carboxylic acid oxygen, and water oxygen, the hydrogen bond angle OG -H -O looks good. Not when it involves histidine atom NE2 or ND1 though. The pH of the crystallization is 7 by the way. Also, there's no water nearby that would draw the HG away from the histidine. Jason

4 5

Fo-Fc omit map for ligand in the protein-ligand structure
by Wei Shi 21 Sep '14

21 Sep '14

Hi all, I am working with a protein dimer-ligand structure (resolution: 2.8Å) and am going to make a figure showing the bound ligand superimposed with the electron density (Fo-Fc) map obtained after omit refinement. What I did is to make the ligand-free structure model and then use Phenix.refine using the following parameters plus simulated annealing (torsion angle) or (Cartesian): Default: XYZ coordinates + Real-space+ Individual B factors + Occupancies. Also click: Optimize X-ray/stereochemistry weight + Optimize X-ray/ADP weight + Secondary structure restraints + NCS. Contoured at 2σ, the Fo-Fc map generated show green density for most part of the ligand, but for several atoms in the middle of the ligand, the green density is missing... It's like this: some atoms have green density and the next few no green density and then a couple few has green density and then a couple more no green denisty.... The missing green density usually means that part is disordered in the structure, right? I am not sure why the density is missing for some atoms in the middle of the ligand. It would makes more sense if the missing green density is at the end of the ligand... Does anyone happen to have any clue about what it means for the green density missing in my case? Also, I am not sure whether I did the omit refinement right: (1). For simulated annealing, it has two options, one is Cartesian and the other is Torsion angle, and it's said that torsion angle is more suitable of low resolution data. Is there any difference about which one to choose in my case? (2). Is it Okay to run the omit refinement in the presence of the following parameters as I did: Optimize X-ray/stereochemistry weight + Optimize X-ray/ADP weight + Secondary structure restraints + NCS. Thank you so much! Best, Wei

5 24

Increase in R/R(free) after phenix.real_space_refine
by George Devaniranjan 20 Sep '14

20 Sep '14

Hi, I refined (phenix.refine) a 1.0Å structure but it had a few CB outliers and some rotamer outliers. I decided to try real_space_refine, to check if they could be fixed (given that it has a higher radius of convergence). (using run=local_grid_search+minimization_global in real_space_refine) But, the R/R(free) increased from 11.39/13.24( for the refine model) to 16.51/17.55 (real_space_refine model). It seemed a big increase, when I looked at the map against the structure the protein itself seemed fine. The rotamers outliers had decreased, the CB outliers were fixed..etc The solvent fit to the density had worsened in COOT. What am I doing wrong? Any suggestions, advice would be much appreciated. Thank you, George

3 2

RMS angles and comprehensive validation in phenix
by C 20 Sep '14

20 Sep '14

5 6

Problems getting rigid and tls to run
by David Shin 19 Sep '14

19 Sep '14

Hi - Running versions: Version: 1.9.1692 mac-intel-osx-x86_64 osx-10.6.8 Version: 1.9.1692 intel-linux-2.6-x86_64 redhat-e5.10 In both cases, I can run an input model I used previously and get rigid and tls to work. With a newer model, I cannot get either to work (and yes, there was a reason I tested rigid on a model that was refined once). Weirder is I have gone from gui on mac to produce an .eff file to transfer to linux - in the .eff file I get: strategy = *individual_sites *individual_sites_real_space rigid_body \ *individual_adp group_adp tls *occupancies group_anomalous despite having checked the rigid and tls boxes Weirder yet, is I hand edit the file to add *rigid *tls and run on the linux machine, and both still get turned off. I diff the input refine.eff files between using the older model and newer, and nothing appears very different. The only thing I can think of is that there could be some legacy .hidden files in the working directory that change something, but it makes no sense once I port everything to a different computer. Anyone have this trouble? Thanks, Dave

1 0

Refining structure with overlapping conformations of different residues
by Igor Petrik 19 Sep '14

19 Sep '14

Hi, I am trying to refine a structure in which there are two sets of linked conformations of two residues. Around half of the occupancy (altloc A) corresponds to the sidechain of residue 29 pointing into the pocket and residue 65 pointing out and the other conformation (altloc B) has the side chain of residue 65 in the pocket (where residue 29 is in altloc A) and residue 29 pointing out of the pocket. I am refining them by using occupancy refinement groups in, but phenix refine seems to think that altloc A of residue 29 and altloc B of residue 65 may exist at the same time and pushes them apart to avoid clashes. Is there any way around this? - Igor

2 1

calculate TLS contribution to total B-factor
by Michael Thompson 19 Sep '14

19 Sep '14

Hello All, I have a structure that was refined with phenix.refine using TLS, so my PDB file contains B-total in the isotropic B-factor column as well as ANISOU lines. For the sake of analysis, I am interested in determining the TLS contribution to these total B-factors. How can I split B-total into isotropic and TLS contributions? I'm fairly certain this is possible, as it is the inverse of what one would do after refining a structure with TLS in Refmac, which requires that the isotropic and TLS components be added together prior to deposition. I feel like perhaps there has been some discussion of this in the past, but I'm not able to dig anything up from the bb archives. Any help would be much appreciated. Thanks, Mike

2 1

Problem finding and running iotbx.explore_metric_symmetry
by Arka Chakraborty 18 Sep '14

18 Sep '14

Dear BB members, I am having a problem finding and running the iotbx.explore_metric_symmetry program. I have tried accessing it from the phenix cctbx directory and also from a separate installation of cctbx but could not locate the program. Could anyone point out how exactly to access it?. Thanks a lot in advance, Best regards, Arko -- *Arka Chakraborty* *ibmb (Institut de Biologia Molecular de Barcelona)* *BARCELONA, SPAIN*

2 1