Search results for "look through"

Re: [phenixbb] Autosol: MIRAS

by ash.k＠aol.com

@Gabor: Thanks Gabor! I'll try MR-SAD. Hope it works. @Andre: I have tried Shelx and Phenix. Se sites from both of them were not the same. I tried giving Shelxd output to Phenix also. That also did not help. I have not tried Autosharp yet. I'll see that too.. However, what now I am primarily looking at is the way to get Se-substructure from the MIRAS run. Is there a way to get that? autosol gives only one heavy atom in the output! Since I have a map which looks promising, so before going to the output of another run I want to explore all the possibilities from this. Se locations in this map should help. AK -----Original Message----- From: Gabor Bunkoczi <gb360(a)cam.ac.uk> To: PHENIX user mailing list <phenixbb(a)phenix-online.org> Sent: Thu, Sep 6, 2012 5:52 am Subject: Re: [phenixbb] Autosol: MIRAS Dear AK, if you just need the Se positions, you can run Phaser MR-SAD (available from the GUI) using your partial model as a starting point, together with the Se-Met data. Dependending on the model quality (and the model does not have to be very good), it may be able to locate the Se-sites. BW, Gabor On Sep 6 2012, ash.k(a)aol.com wrote: > Satisfactory map means: I am expecting a coiled coil or helix bundle type > of assembly and I could see some densities appropriate for helices. There > are proper solvent channels and continuous stretches of densities. > > More to add about data on this: the data is anisotropic and the longer > 'c' axis and alignment of helical density along c axis support this. This > also makes me think that perhaps map is sensible. > > Few rounds of model building, refinement and DM has been successful to > assign around 10 polyala helices and their distribution looks sensible > from packing point of view. Now the problem is to assign the side chain. > I was hoping to make use of Se locations for this. R and Rfree are still > random, in the range of 50%. I am not sure, but it could be perhaps > because most of the scattering material is still unassigned. > > > > > > >-----Original Message----- >From: Francis E Reyes <Francis.Reyes(a)colorado.edu> >To: PHENIX user mailing list <phenixbb(a)phenix-online.org> >Cc: phenixbb <phenixbb(a)phenix-online.org> >Sent: Thu, Sep 6, 2012 8:07 am >Subject: Re: [phenixbb] Autosol: MIRAS > > > This one is solvable, but with extreme difficulty. I recently completed a > structure solution with experimental phases starting at 5.0 A using phase > information from multiple derivatives. > > >How would you describe a somewhat satisfactory map? > > >F > >On Sep 5, 2012, at 7:08 PM, ash.k(a)aol.com wrote: > > > > > > > >Hi Shya, > > I did wavelength scan, got a good signal for Se and used appropriate > wavelengths for data collection and also used experimental f' and f'' > values for phasing. I think the reasons SAD or MAD for SeMet data is not > working are (i) low resolution: 3.7 for SeMet (Anomalous data is up to > 4.8A) (ii) I should have mentioned this earlier: 3 out of 6 Se are very > close to N-terminus, possible they are disordered. Unit cell is also some > what big..100, 120 and 320A; F222 space group. > >AK > > > >-----Original Message----- >From: Shya Biswas <shyabiswas(a)gmail.com> >To: PHENIX user mailing list <phenixbb(a)phenix-online.org> >Sent: Thu, Sep 6, 2012 7:29 am >Subject: Re: [phenixbb] Autosol: MIRAS > > >Hi AK, >Did you do a wavelength scan when you collected the SE dataset you >need to put the values of f' and f'' from your wavelength scan in >order to locate the heavy atom sites, 6 methionine should be enough to >phase your molecule. >Shya > >On Wed, Sep 5, 2012 at 9:25 PM, <ash.k(a)aol.com> wrote: >> Hi all, >> >> I am trying to solve a structure through experimental phasing using >> AUTOSOL. I have a couple of heavy atom derivative datasets (Hg, La, Eu, >> Cd) and also a SeMet data. Unfortunately all the datasets are of low >> resolution (3.7-4.2A) and there are possibly 4-8 molecules in the asu. >> MIR, SAD and MAD alone did not give any convincing solution. >> >> However, MIRAS, with a combination of few heavy atom datasets and the >> anomalous data from SeMet crystals, gave a somewhat satisfactory map. >> But the heavy atom site picked by AUTOSOL list only one of the heavy >> atoms i.e. Lanthanum. In another set of run, the solution of which was >> not convincing, the heavy atom substructure had only Hg. There are 6 Met >> out of 200 residues in one molecule and mass spec results show that Se >> incorporation is 100%. >> >> Now, my doubt is that why does the heavy atom substructure contain only >> La and how can I get the substructure involving Se from this solution >> (or the datasets used)? Se location is going to help me a lot for >> finding a starting point to assign side chains. >> >> Any suggestion would be greatly appreciated. >> >> Thanks >> AK >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >> >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb > > > > >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb > > > >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb > > > -- ################################################## Dr Gabor Bunkoczi Cambridge Institute for Medical Research Wellcome Trust/MRC Building Addenbrooke's Hospital Hills Road Cambridge CB2 0XY ################################################## _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb

12 years, 9 months

Re: [phenixbb] temp files

by James Holton

What if people are using the GUI ? On 4/24/2021 2:51 PM, Tom Terwilliger wrote: > Hi James, > > Not the code (I hope), just any scripts that use these methods. Like this: > > mkdir /var/tmp/autosol > phenix.autosol p9.sca 2 se temp_dir=/var/tmp/autosol > > Now the temp files go in /var/tmp/autosol and the output files go in > AutoSol_run_xxx/ as usual > > All the best, > Tom > > All the best, > Tom T > > On Sat, Apr 24, 2021 at 11:59 AM James Holton <jmholton(a)lbl.gov > <mailto:[email protected]>> wrote: > > Thank you Tom! > > Ok. So, in order to change the default I need to go through the > code looking for "temp_dir" and change things? > > > On 4/24/2021 10:16 AM, Tom Terwilliger wrote: >> Hi James, >> >> There is no overall Phenix temp directory specification, but most >> of the temp_dir usage is from >> autosol/autobuild/ligandfit/map_to_model. Each of these has the >> keyword "temp_dir=xxxx" which you should be >> able to set to any directory you want (and local is better as you >> note). Most programs using a temp_dir also have a keyword >> clean_up=True as well. >> >> All the best, >> Tom T >> >> On Sat, Apr 24, 2021 at 10:51 AM James Holton <jmholton(a)lbl.gov >> <mailto:[email protected]>> wrote: >> >> Thank you Li-Wei >> >> Definitely not placing blame on one program. Phenix.autobuild >> is another >> big temp file producer. So is XDS. Clearly this ligand run >> was a case >> of a misconfigured, runaway task that never finished. >> However, the files >> lingered on disk, eating up inodes for 3 years! >> >> The reason I'm asking is I think there are significant >> performance >> increases to be gained by using fast, local storage for >> scratch files. >> This is not just in speed but storage and overall system/cluster >> performance. Very few things are more expensive than an NFS >> write! >> >> Does anyone know how to change the default temp file location >> across >> phenix ? Is this a cctbx thing? >> >> Thanks >> >> -James >> >> >> On 4/23/2021 9:38 PM, Li-Wei Hung wrote: >> > Hi James, >> > >> > I'll leave the global Phenix temp aspect to Billy. >> > For ligand identification specifically, the working >> directory is where >> > all the files are located. The program will purge most of the >> > intermediate files upon completion. If the user interrupted >> the runs >> > or if the program crashed at certain spots, the purge >> mechanism might >> > not kick in. Even so, it'd take many runs to accumulate >> 20e6 (2e7?) >> > files. In any case, you've got a point and I'll look into >> salvaging >> > intermediate files of ligand identification as soon as they >> are not >> > needed in the process. >> > >> > Thanks, >> > >> > Li-Wei >> > >> > On 4/23/2021 7:03 PM, James Holton wrote: >> >> Hello all, >> >> >> >> Is there a way to configure phenix at install time (or >> perhaps >> >> post-install) to put temporary files under /tmp ? I just >> had to >> >> delete 20e6 temp files over NFS from a single user's >> phenix ligand >> >> identification run. The delete took almost a month. >> >> >> >> Apologies if I am neglecting to look somewhere obvious in the >> >> documentation, >> >> >> >> Happy Weekend! >> >> >> >> -James Holton >> >> MAD Scientist >> >> >> >> _______________________________________________ >> >> phenixbb mailing list >> >> phenixbb(a)phenix-online.org <mailto:[email protected]> >> >> http://phenix-online.org/mailman/listinfo/phenixbb >> <http://phenix-online.org/mailman/listinfo/phenixbb> >> >> Unsubscribe: phenixbb-leave(a)phenix-online.org >> <mailto:[email protected]> >> > >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org <mailto:[email protected]> >> http://phenix-online.org/mailman/listinfo/phenixbb >> <http://phenix-online.org/mailman/listinfo/phenixbb> >> Unsubscribe: phenixbb-leave(a)phenix-online.org >> <mailto:[email protected]> >> >> >> >> -- >> Thomas C Terwilliger >> Laboratory Fellow, Los Alamos National Laboratory >> Senior Scientist, New Mexico Consortium >> 100 Entrada Dr, Los Alamos, NM 87544 >> Email: tterwilliger(a)newmexicoconsortium.org >> <mailto:[email protected]> >> Tel: 505-431-0010 >> > > > > -- > Thomas C Terwilliger > Laboratory Fellow, Los Alamos National Laboratory > Senior Scientist, New Mexico Consortium > 100 Entrada Dr, Los Alamos, NM 87544 > Email: tterwilliger(a)newmexicoconsortium.org > <mailto:[email protected]> > Tel: 505-431-0010 >

4 years, 2 months

Re: [cctbxbb] Niggli-reduced cell C++ implementation

by Martin Uhrin

Dear cctbxers, I've finally found the time to play around with a C++ version of the KG algorithm and I've come across a result I don't understand. I've tried both David's C++ and the cctbx python niggli_cell() implementations and they both give the roughly the same answer. I'm reducing the following cell with two, equivalent, representations (a, b, c, alpha, beta, gamma): Before: 1: 4.630811 4.630811 4.630811 90 90 90 2: 3.27448 5.67156 5.67156 99.5941 106.779 90 After: 1: 4.63081 4.63081 4.63081 90 90 90 2: 3.27448 5.67154 5.67156 99.5941 90 106.778 Looking at the trace, cell 1 undergoes step 3 and finishes while cell 2 undergoes steps 2, 3, 7 and 4. Does anyone know why these haven't converged to the same cell? Many thanks, -Martin On 23 March 2012 17:12, Ralf Grosse-Kunstleve <rwgrosse-kunstleve(a)lbl.gov>wrote: > Hi Martin, > Let me know if you need svn write access to check in your changes. All I > need is your sourceforge user id. > Ralf > > > On Fri, Mar 23, 2012 at 3:35 AM, Martin Uhrin <martin.uhrin.10(a)ucl.ac.uk>wrote: > >> Dear David and Rolf, >> >> thank you for your encouragement. >> >> David: I'm more than happy to port your implementation to cctbx if you're >> happy with this. Of course I don't want to step on your toes so if you'd >> rather do it yourself (or not at all) that's cool. >> >> There may be some licensing issues to sort out as it looks like cctbx has >> a custom (non viral) license but the BSD license is likely compatible. >> >> On first impression I think a new class would be the way to go but I'd >> have to look at the two algorithms in greater detail to be sure. >> >> All the best, >> -Martin >> >> >> On 22 March 2012 22:00, Ralf Grosse-Kunstleve <rwgrosse-kunstleve(a)lbl.gov >> > wrote: >> >>> Hi Martin, >>> You're very welcome to add a C++ version of the Krivy-Gruber algorithm >>> to cctbx if that's what you had in mind. >>> I'm not sure what's better, generalizing the fast-minimum-reduction >>> code, or just having an independent implementation. >>> Ralf >>> >>> On Thu, Mar 22, 2012 at 2:24 PM, Martin Uhrin <martin.uhrin.10(a)ucl.ac.uk >>> > wrote: >>> >>>> Dear Cctbx community, >>>> >>>> Firstly I'd like to say thank you to Rolf, Nicholas and Paul for their >>>> expertly thought through implementation of the reduced cell algorithm. >>>> I've found it to be extremely useful for my work. >>>> >>>> My code is all in C++ and I'd like to be able to use the Krivy-Gruber >>>> algorithm. My understanding is that only the reduced (Buerger) unit cell >>>> algorithm is implemented in C++ [1] which guarantees shortest lengths but >>>> not unique angles. From my understanding the Krivy-Gruber would also >>>> guarantee me uniqueness of unit cell angles, however this is only >>>> implemented in Python [2]. Sorry to be so verbose, I just wanted to check >>>> that I was on the right page. >>>> >>>> Would it be possible for me to implement the Krivy-Gruber in C++ by >>>> adding in the epsilon_relative to the parameter and following the procedure >>>> found in the python version? >>>> >>>> Many thanks, >>>> -Martin >>>> >>>> [1] >>>> http://cctbx.sourceforge.net/current/c_plus_plus/classcctbx_1_1uctbx_1_1fas… >>>> [2] >>>> http://cctbx.sourceforge.net/current/python/cctbx.uctbx.krivy_gruber_1976.h… >>>> >>>> >>>> >>>> -- >>>> Martin Uhrin Tel: +44 >>>> 207 679 3466 >>>> Department of Physics & Astronomy Fax:+44 207 679 0595 >>>> University College London >>>> martin.uhrin.10(a)ucl.ac.uk >>>> Gower St, London, WC1E 6BT, U.K. http://www.cmmp.ucl.ac.uk >>>> >>>> _______________________________________________ >>>> cctbxbb mailing list >>>> cctbxbb(a)phenix-online.org >>>> http://phenix-online.org/mailman/listinfo/cctbxbb >>>> >>>> >>> >>> _______________________________________________ >>> cctbxbb mailing list >>> cctbxbb(a)phenix-online.org >>> http://phenix-online.org/mailman/listinfo/cctbxbb >>> >>> >> >> >> -- >> Martin Uhrin Tel: +44 >> 207 679 3466 >> Department of Physics & Astronomy Fax:+44 207 679 0595 >> University College London >> martin.uhrin.10(a)ucl.ac.uk >> Gower St, London, WC1E 6BT, U.K. http://www.cmmp.ucl.ac.uk >> >> _______________________________________________ >> cctbxbb mailing list >> cctbxbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/cctbxbb >> >> > > _______________________________________________ > cctbxbb mailing list > cctbxbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/cctbxbb > > -- Martin Uhrin Tel: +44 207 679 3466 Department of Physics & Astronomy Fax:+44 207 679 0595 University College London martin.uhrin.10(a)ucl.ac.uk Gower St, London, WC1E 6BT, U.K. http://www.cmmp.ucl.ac.uk

13 years, 2 months

Re: [cctbxbb] some thoughts on cctbx and pip

by Luc Bourhis

Hi, I did look into that many years ago, and even toyed with building a pip installer. What stopped me is the exact conclusion you reached too: the user would not have the pip experience he expects. You are right that it is a lot of effort but is it worth it? Considering that remark, I don’t think so. Now, Conda was created specifically to go beyond pip pure-python-only support. Since cctbx has garnered support for Conda, the best avenue imho is to go the extra length to have a package on Anaconda.org <http://anaconda.org/>, and then to advertise it hard to every potential user out there. Best wishes, Luc > On 16 Aug 2019, at 21:45, Aaron Brewster <asbrewster(a)lbl.gov> wrote: > > Hi, to avoid clouding Dorothee's documentation email thread, which I think is a highly useful enterprise, here's some thoughts about putting cctbx into pip. Pip doesn't install non-python dependencies well. I don't think boost is available as a package on pip (at least the package version we use). wxPython4 isn't portable through pip (https://wiki.wxpython.org/How%20to%20install%20wxPython#Installing_wxPython… <https://wiki.wxpython.org/How%20to%20install%20wxPython#Installing_wxPython…>). MPI libraries are system dependent. If cctbx were a pure python package, pip would be fine, but cctbx is not. > > All that said, we could build a manylinux1 version of cctbx and upload it to PyPi (I'm just learning about this). For a pip package to be portable (which is a requirement for cctbx), it needs to conform to PEP513, the manylinux1 standard (https://www.python.org/dev/peps/pep-0513/ <https://www.python.org/dev/peps/pep-0513/>). For example, numpy is built according to this standard (see https://pypi.org/project/numpy/#files <https://pypi.org/project/numpy/#files>, where you'll see the manylinux1 wheel). Note, the manylinux1 standard is built with Centos 5.11 which we no longer support. > > There is also a manylinux2010 standard, which is based on Centos 6 (https://www.python.org/dev/peps/pep-0571/ <https://www.python.org/dev/peps/pep-0571/>). This is likely a more attainable target (note though by default C++11 is not supported on Centos 6). > > If we built a manylinuxX version of cctbx and uploaded it to PyPi, the user would need all the non-python dependencies. There's no way to specify these in pip. For example, cctbx requires boost 1.63 or better. The user will need to have it in a place their python can find it, or we could package it ourselves and supply it, similar to how the pip h5py package now comes with an hd5f library, or how the pip numpy package includes an openblas library. We'd have to do the same for any packages we depend on that aren't on pip using the manylinux standards, such as wxPython4. > > Further, we need to think about how dials and other cctbx-based packages interact. If pip install cctbx is set up, how does pip install dials work, such that any dials shared libraries can find the cctbx libraries? Can shared libraries from one pip package link against libraries in another pip package? Would each package need to supply its own boost? Possibly this is well understood in the pip field, but not by me :) > > Finally, there's the option of providing a source pip package. This would require the full compiler toolchain for any given platform (macOS, linux, windows). These are likely available for developers, but not for general users. > > Anyway, these are some of the obstacles. Not saying it isn't possible, it's just a lot of effort. > > Thanks, > -Aaron > > _______________________________________________ > cctbxbb mailing list > cctbxbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/cctbxbb

5 years, 10 months

Re: [phenixbb] phaser MR

by Randy Read

Hi, If you're unlucky, a combination of pseudo-symmetry and twinning gives only weak indications of twinning, because the twin-related reflections are correlated but the twinning tests assume that they are not. It might be worth considering that possibility. You could try solving the structure in lower symmetry space groups to see whether you get a solution that does not exactly obey the higher symmetry. Are you really certain about the space group? If one symmetry operator is wrong (e.g. a 2-fold instead of a 2-fold screw), some of the symmetry operators are still right, and you can get an apparent solution that partly explains the data. It might be worth running tests in all the variants of P422 as well. Best wishes Randy Read On 23 Jul 2009, at 22:54, sbiswas2(a)ncsu.edu wrote: > Hi Francis, > Thanks for your response. The matthews coefficient suggests two > molecules > in the AU. Phaser also finds two molecules. I ran the dataset through > phenix xtriage it did not indicate twinning though. The molecule also > exists in nature as a monomer. > Shya > > >> Twinning? What's your matthews coefficient say? Do you know if your >> structure is a multimer (biochemistry, etc)? Does it agree with the >> matthews coefficient? >> >> If the unit cell is not big enough to hold all of the contents,then >> this is an indicator for twinning . >> >> FR >> >> On Jul 23, 2009, at 3:09 PM, sbiswas2(a)ncsu.edu wrote: >> >>> Hi all, >>> >>> I was trying to solve a structure by molecular replacement. I scaled >>> the >>> data in P4222 space group (resolution 2.7A) with two molecules in >>> the >>> assymmetric unit (molecule A and B) I ran phaser with my model and >>> got a >>> Zscore of 5.1. When I look at the map that I got from phaser I could >>> easily see good electron density for both molecules, However upon >>> inspection of the electron density map there were considerable >>> interaction >>> or clashes with molecule B and a symmetry atom. Molecule A had no >>> clashes >>> however with the symmetry atoms. I was wondering if anyone knows how >>> to >>> resolve this. Could it be a problem of space group. The statistics >>> are >>> good for space group P4222 and the I/sigI was good till 2.7A. >>> Any advice is appreciated, >>> Shya >>> >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> --------------------------------------------- >> Francis Reyes M.Sc. >> 215 UCB >> University of Colorado at Boulder >> >> gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D >> >> 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb ------ Randy J. Read Department of Haematology, University of Cambridge Cambridge Institute for Medical Research Tel: + 44 1223 336500 Wellcome Trust/MRC Building Fax: + 44 1223 336827 Hills Road E-mail: rjr27(a)cam.ac.uk Cambridge CB2 0XY, U.K. www- structmed.cimr.cam.ac.uk

15 years, 11 months

Re: [phenixbb] Selecting ellipsoid of data

by Kay Diederichs

Hi Frank, "such a tool" is at http://strucbio.biologie.uni-konstanz.de/xdswiki/index.php/Aniso_cutoff where it's meant to be applied to INTEGRATE.HKL which comes out of XDS. Doing it this way has the benefit that the statistics printed out by XDS' CORRECT (or SCALA/TRUNCATE; there are people who prefer that route) match the data you refine agains. HTH, Kay > Message: 1 > Date: Thu, 15 Jul 2010 06:32:40 +0100 > From: Frank von Delft<frank.vondelft(a)sgc.ox.ac.uk> > To: PHENIX user mailing list<phenixbb(a)phenix-online.org> > Subject: [phenixbb] Selecting ellipsoid of data > Message-ID:<4C3E9D78.5030203(a)sgc.ox.ac.uk> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi, is there a tool in phenix that allows me to select an ellipsoid of > data -- specified e.g. by the highest resolutions in three reciprocal > lattice directions. (Yes, I'm playing with anisotropy, "playing" being > the operative word.) > > phx > > > ------------------------------ > > Message: 2 > Date: Wed, 14 Jul 2010 22:47:30 -0700 > From: "Ralf W. Grosse-Kunstleve"<rwgk(a)cci.lbl.gov> > To: phenixbb(a)phenix-online.org > Subject: Re: [phenixbb] Selecting ellipsoid of data > Message-ID:<201007150547.o6F5lUoL018490(a)cci.lbl.gov> > Content-Type: text/plain; charset=us-ascii > > Hi Frank, > >> Hi, is there a tool in phenix that allows me to select an ellipsoid of >> data -- specified e.g. by the highest resolutions in three reciprocal >> lattice directions. (Yes, I'm playing with anisotropy, "playing" being >> the operative word.) > > I'm not aware of such a tool. > > Ralf > > > ------------------------------ > > Message: 3 > Date: Thu, 15 Jul 2010 07:15:58 +0100 > From: Frank von Delft<frank.vondelft(a)sgc.ox.ac.uk> > To: "Ralf W. Grosse-Kunstleve"<rwgk(a)cci.lbl.gov>, PHENIX user mailing > list<phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] Selecting ellipsoid of data > Message-ID:<4C3EA79E.2030503(a)sgc.ox.ac.uk> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi Ralf > > Yeah, I figured. So if I want to use cctbx, where do I start? Just a > pointer to a) package and b) function where I'll see the syntax. > > So equation for ellipsoid is x^2/a^2 + y^2/b^2 + z^2/c^2 = 1; so I > imagine I take each reflection, convert each of h,k,l to 1/reso, and > with a = 1/res(a*), I just check whether the above is< 1. > > The main thing I still need is to convert h,k,l into orthogonal > coordinates.... or do I? I suppose I don't, as what I care for is not > whether it's "really" an ellipsoid, only whether it cuts through miller > index space anisotropically. > > Hmmmm... I may be able to do it in sftools; but if you can in<1minute > give me a link to where to look to get started in cctbx, that would be > awesome. > > (Thanks for listening :) > > > > >> Hi Frank, >> >> >>> Hi, is there a tool in phenix that allows me to select an ellipsoid of >>> data -- specified e.g. by the highest resolutions in three reciprocal >>> lattice directions. (Yes, I'm playing with anisotropy, "playing" being >>> the operative word.) >>> >> I'm not aware of such a tool. >> >> Ralf

14 years, 11 months

Re: [phenixbb] Third issue (...) - P.S.

by Frank von Delft

Thanks for the elaboration. All I meant was that figure 2 did not show improvements (points pretty much on the y=x line), while figure 1 showed huge improvements (points below the y=x line). I should have added that in the past I've had puzzled questions around here about phenix results at lower resolutions, and I always wondered whether weight optimization may be an issue. So maybe what I was looking for was a figure like figure 2, but of the old optimization. Cheers Frank On 29/07/2011 17:41, Pavel Afonine wrote: > Hi Frank, > > thanks for prompting me to have another look at this! Yes, a quick > section reiterating the results shown in both figures would be > definitely helpful. Ok, I'm doing it now (better late than never!): > > Figure 1 (old version vs new version), going from left to right, top to > bottom: > > - the new procedure always produces smaller Rfree-Rwork gap (except two > outliers out of ~110 cases). In a few dozens cases this means going from > Rfree-Rwork ~10-14% down to 8% and less. Which is I believe an improvement. > > - Molprobity clashscore gets much lower in all cases, and in many cases > it get down from clashscore>150-200 to less than 50. Again, I take this > as an improvement. > > - Ramachandran plot: always less outliers, and in a few dozens cases it > get down from 10-25% to less than 5%. It's good too! > > - Always less Rotamer outliers. Not bad. > > -<Bij> for bonded atoms gets always smaller. In many cases it goes down > from ridiculously high ~30-50A**2 to less than 20 or so. There are > outliers which, as pointed out in the text, need investigation. > > - Rfree difference wise (the last graphic at the bottom), almost all the > differences stay within "equally good range of values", which is defined > in the text. Some bad outliers, are 1) inoptimal NCS groups selection > (which was done automatically using the old procedure, not the new one > that does it in torsion angle space), 2) not using TLS even though this > was used originally, and other issues that needs an investigation. There > are improvements too. > > Figure 2 (with weight optimization vs without weight optimization): > going from left to right, top to bottom: > > - Rfree-Rfree gets much better in a few dozens of cases, or stays in a > "equally good range of values". > > - clashscore either stays the same or gets smaller, which is good. Same > for Ramachandran and Rotamer outliers. > > - Except two outliers (out of total ~110 cases), Rfree(opt)-Rfree(no > opt) is consistently smaller, with a number of cases having significant > improvement (over ~2% or so). > > Thanks again for your comments, and prompting me to spell out a proper > analysis of the presented pictures. > > Pavel. > > > On 7/29/11 1:21 AM, Frank von Delft wrote: >> Hi, just scanned through the article on automatic weight adjustment. >> >> If I had to summarise figures 1 and 2, I'd have to conclude that the >> latest weight optimization only occasionally produces better results >> than not doing it (fig 2); and that therefore the old optimization >> was considerably worse than doing nothing at all. >> >> Or have I misinterpreted it? >> Cheers >> phx >> >> >> >> On 28/07/2011 19:59, Nigel Moriarty wrote: >>> Dear Colleagues, >>> >>> I am pleased to announce the publication of the third issue of the >>> Computational Crystallography Newsletter: >>> >>> http://www.phenix-online.org/newsletter/ >>> >>> A listing of the articles and short communications is given below. >>> Please note that the newsletter accepts articles of a general nature >>> of interest to all crystallographers. Please send any articles to me >>> at >>> NWMoriarty(a)lbl.gov noting that there is a Word Template on the website >>> to streamline production. >>> >>> Articles >>> -------- >>> >>> Improved target weight optimization in phenix.refine >>> Mite-y lysozyme crystal and structure >>> >>> >>> Short communications >>> -------------------- >>> >>> A lightweight, versatile framework for visualizing reciprocal-space data >>> An extremely fast spotfinder for real-time beamline applications >>> Hints for running phenix.mr_rosetta >>> >>> >>> Cheers >>> Nigel >>> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/phenixbb

13 years, 11 months

Re: [phenixbb] phenix.automr high zscore but too many clashes

by Georg Zocher

Dear Zach, are you sure about your spacegroup? And what about the matthews parameter (http://www.ruppweb.org/Mattprob/)? Best Regards, Georg Andreas Förster schrieb: > Hey Zach, > > remove the loops from your search model and rerun the MR. > Alternatively, increase the maximum number of clashes allowed until you > get a solution written out. See where your solutions clash. Again, if > it's in a loop region, cut the loop and rerun. Your z-scores look > great. (Notice that you have zed-scores when you run phaser through > ccp4i and zee-scores when you run phaser through phenix.) > > > Andreas > > > > zach powers wrote: > >> Hi, >> >> >> I have a problem and i wonder if someone has had a similar experience. I >> have been using phenix.automr program to find an MR solution for a >> protein. When I ask Phaser to use copies=2 I find a number of solutions >> with poor z-scores (3-4). When I ask it to use copies=3, I get great >> Z-scores (see below) but no solution due to the high number of clashes >> (>100!). >> >> As an x-ray newbie I am a bit perplexed of what to make of this: the >> solution is not good because there is not enough space in the unit cell >> to comfortably accommodate all three molecules yet the Z-score indicates >> the structure is good. >> >> My structure does have several loops and these may be contributing to >> the clashing residues. As a newbie I have a newbie question - what does >> this mean (great z-score but no solutions)? Is this non-solution a >> possible solution if I play with it, or is the packing simply too tight >> and the Z-scores are not valid? >> >> In the meantime I am chopping my protein into sub-domains as recommended >> in the Phaser documentation, but if anyone has seen something like this >> before, any suggestions are welcome. >> >> thanks >> zach charlop-powers >> >> >> >> >> Packing Table: Space Group P 3 >> ------------------------------ >> Solutions accepted if number of clashes = 135 (lowest number of >> clashes in list) >> provided this number of clashes <= 10 (maximum number of allowed clashes) >> # #Clashes # Accepted Annotation >> 1 188 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=10.4 >> 2 156 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=10.2 >> 3 151 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=9.8 >> 4 146 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=9.7 >> 5 187 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=16.3 >> TFZ=10.5 >> 6 202 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=9.2 >> 7 170 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=9.1 >> 8 144 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=16.3 >> TFZ=10.2 >> 9 177 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=9.1 >> 10 153 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.9 >> 11 150 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.8 >> 12 217 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.6 >> 13 182 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.6 >> 14 135 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.5 >> 15 175 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.3 >> 16 159 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.3 >> 17 192 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.3 >> 18 165 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=16.3 >> TFZ=9.2 >> 19 169 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=16.3 >> TFZ=9.1 >> 20 158 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.0 >> 21 173 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=8.0 >> 22 154 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=16.3 >> TFZ=8.9 >> 23 206 NO RFZ=20.3 TFZ=18.2 PAK=10 LLG=431 RFZ=18.5 >> TFZ=7.8 >> >> 0 accepted of 23 solutions >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> > > -- Universität Tübingen Interfakultäres Institut für Biochemie Dr. Georg Zocher Hoppe-Seyler-Str. 4 72076 Tuebingen Germany Fon: +49(0)-7071-2973374 Mail: georg.zocher(a)uni-tuebingen.de http://www.ifib.uni-tuebingen.de

16 years, 2 months

Re: [cctbxbb] Git

by Graeme.Winter＠diamond.ac.uk

Nick I offer an opinion here - git is like C++ - in the same way as Fortran 77 is like SVN. You can code fortran in any language syntax you like just fine. If you write it using C++ syntax you need to use g++ to compile it but that’s it. You can do that without using any C++ idioms. Or you can go the full monty and code with STL metaprogramming & such - much more complex and more powerful but not mandatory. You can use git like svn which works a bit better, that’s OK. It does not stop other people from using branches and such - you just run off the trunk like you do with SVN and let them get on with it. Every so often you will see a massive commit come in, from someone who did use git more powerfully, but that is no different to someone running over svn. What it does mean though is if you have a bunch of changes and one is a bug fix, you can easily cherry pick just the bug fix to push back for example, leaving the rest of the work in progress … in progress. I was that old fashioned developer Luis alluded to, I did not like it much when I started and now I find it *awesome* - you can commit as you go and only push when you are ready, which makes it really easy to keep your head clearer on big things Cheers Graeme On 12 Jan 2016, at 17:52, Nicholas Sauter <nksauter(a)lbl.gov<mailto:[email protected]>> wrote: >From the DIALS-West perspective, the switch to git has been a stumbling block to participation. I would have to agree with Ralf that git seems to be a tool for very smart people but not for folks who just want a simple tool for managing code. I can't agree with Markus that a linear history is dispensable. In fact, this is one feature in svn that I've found very helpful over the years. If a feature is broken today, but I know for a fact that it worked sometime in the past, I simply do a svn update -r "[datestamp]" to narrow down the exact date when the feature became broken, then I isolate the exact commit and look at the code. Can this be done with git or does it even make sense if there is no concept of linear change? The git environment seems a bit chaotic, and not conducive to close cooperation, from what I can see now. More experience with it may convince me otherwise. Nick Nicholas K. Sauter, Ph. D. Computer Staff Scientist, Molecular Biophysics and Integrated Bioimaging Division Lawrence Berkeley National Laboratory 1 Cyclotron Rd., Bldg. 33R0345 Berkeley, CA 94720 (510) 486-5713 On Tue, Jan 12, 2016 at 2:22 AM, <luis.fuentes-montero(a)diamond.ac.uk<mailto:[email protected]>> wrote: Hi fellows, I believe it is not a bad idea to move to https://github.com<https://github.com/> (the servers) and later on decide if move to use git (the CLI tool). This can an incremental way to do transition. And will allow really old fashion developers (those who refuse to learn git) or conservative policies to survive and still move to a new more reliable code repository. Just my thoughts, Greetings, Luiso From: cctbxbb-bounces(a)phenix-online.org<mailto:[email protected]> [mailto:[email protected]<mailto:[email protected]>] On Behalf Of markus.gerstel(a)diamond.ac.uk<mailto:[email protected]> Sent: 12 January 2016 09:56 To: cctbxbb(a)phenix-online.org<mailto:[email protected]> Subject: Re: [cctbxbb] Git Dear Luc, We actually had a look at the quality of the Github svn interface when we prepared for the DIALS move. My verdict was that it is surprisingly useful. All the possible git operations are mapped rather well onto a linear SVN history. All the branches and tags are accessible either by checking out the root of the repository, ie. svn co https://github.com/dials/dials.git (not recommended, as, say, any new branch/tag will probably result in a huge update operation for you) or by checking them out directly, ie. svn co https://github.com/dials/dials.git/branches/fix_export or svn co https://github.com/dials/dials.git/tags/1.0 Forks by other people are, as is the case when using git, simply not visible. I agree that looking at the project history through SVN may not be as clear as the git history. I wonder how relevant this is though, as you can explore the history, without running any commands, on the Github website, eg. https://github.com/dials/dials/network, https://github.com/dials/dials/commits/master, etc. Creating new branches and tags and merging stuff via SVN would be a major operation. However, that has always been the case with SVN, and for that reason one generally just does not do these things in SVN. But for an SVN user group those operations would not be important – you every only need to merge stuff if you create branches – so I don’t really see the problem. If you want to tag releases you can do that on the Github website as well as with a git command. In summary, I do recognize that SVN users will have difficulties participating in branched development, and in particular will not be able to quickly switch between branches. But I don’t think that this will be a problem, or that there is a need for a policy to keep the history linear. Best wishes -Markus From: cctbxbb-bounces(a)phenix-online.org<mailto:[email protected]> [mailto:[email protected]] On Behalf Of Luc Bourhis Sent: 11 January 2016 16:45 To: cctbx mailing list Subject: Re: [cctbxbb] Gi, Hi Graeme, Can we revisit the idea of moving to git for cctbx? This brings to mind a question I have been asking myself since the subject has been brought forth. The idea Paul wrote about on this list was a move to Github. I guess some, perhaps many, developers will keep interacting with the repository using subversion. I am worried this would clash with the workflow of those of us who would go the native git way. By that I mean creating many branches and merge points, which one would merge into the official repository when ready. I am worried the history would look very opaque for the subversion users. I would even probably create a fork on Github, making it even more opaque for a tool like subversion. Has anybody thought about such issues? My preferred solution would be for everybody to move to git but I don’t think that’s realistic. At the other end of the spectrum, there is putting in place a policy to keep the history linear. Best wishes, Luc -- This e-mail and any attachments may contain confidential, copyright and or privileged material, and are for the use of the intended addressee only. If you are not the intended addressee or an authorised recipient of the addressee please notify us of receipt by returning the e-mail and do not use, copy, retain, distribute or disclose the information in or attached to the e-mail. Any opinions expressed within this e-mail are those of the individual and not necessarily of Diamond Light Source Ltd. Diamond Light Source Ltd. cannot guarantee that this e-mail or any attachments are free from viruses and we cannot accept liability for any damage which you may sustain as a result of software viruses which may be transmitted in or with the message. Diamond Light Source Limited (company no. 4375679). Registered in England and Wales with its registered office at Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire, OX11 0DE, United Kingdom -- This e-mail and any attachments may contain confidential, copyright and or privileged material, and are for the use of the intended addressee only. If you are not the intended addressee or an authorised recipient of the addressee please notify us of receipt by returning the e-mail and do not use, copy, retain, distribute or disclose the information in or attached to the e-mail. Any opinions expressed within this e-mail are those of the individual and not necessarily of Diamond Light Source Ltd. Diamond Light Source Ltd. cannot guarantee that this e-mail or any attachments are free from viruses and we cannot accept liability for any damage which you may sustain as a result of software viruses which may be transmitted in or with the message. Diamond Light Source Limited (company no. 4375679). Registered in England and Wales with its registered office at Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire, OX11 0DE, United Kingdom _______________________________________________ cctbxbb mailing list cctbxbb(a)phenix-online.org<mailto:[email protected]> http://phenix-online.org/mailman/listinfo/cctbxbb _______________________________________________ cctbxbb mailing list cctbxbb(a)phenix-online.org<mailto:[email protected]> http://phenix-online.org/mailman/listinfo/cctbxbb

9 years, 5 months

Re: [cctbxbb] bootstrap.py build on Ubuntu

by David Waterman

Hi Billy, I'm replying on this old thread because I have finally got round to trying a bootstrap build for DIALS out again on Ubuntu, having waited for updates to the dependencies and updating the OS to 16.04. The good news is, the build ran through fine. This is the first time I've had a bootstrap build complete without error on Ubuntu, so thanks to you and the others who have worked on improving the build in the last few months! The bad news is I'm getting two failures in the DIALS tests: dials/test/command_line/tst_export_bitmaps.py dials_regression/test.py Both are from PIL File "/home/fcx32934/dials_test_build/base/lib/python2.7/site-packages/PIL/Image.py", line 401, in _getencoder raise IOError("encoder %s not available" % encoder_name) IOError: encoder zip not available Indeed, from base_tmp/imaging_install_log it looks like PIL is not configured properly -------------------------------------------------------------------- PIL 1.1.7 SETUP SUMMARY -------------------------------------------------------------------- version 1.1.7 platform linux2 2.7.8 (default_cci, Jun 10 2016, 16:04:32) [GCC 5.3.1 20160413] -------------------------------------------------------------------- *** TKINTER support not available *** JPEG support not available *** ZLIB (PNG/ZIP) support not available *** FREETYPE2 support not available *** LITTLECMS support not available -------------------------------------------------------------------- Any ideas? I have zlib headers but perhaps PIL can't find them. On a related note, the free version of PIL has not been updated for years. The replacement Pillow has started to diverge. I first noticed this when Ubuntu 16.04 gave me Pillow 3.1.2 and my cctbx build with the system python produced failures because it no longer supports certain deprecated methods from PIL. I worked around that in r24587, but these things are a losing battle. Is it time to switch cctbx over to Pillow instead of PIL? Cheers -- David On 7 January 2016 at 18:12, Billy Poon <bkpoon(a)lbl.gov> wrote: > Hi all, > > Since wxPython was updated to 3.0.2, I have been thinking about updating > the other GUI-related packages to more recent versions. I would probably > update to the latest, stable version that does not involve major changes to > the API so that backwards compatibility is preserved. Let me know if that > would be helpful and I can prioritize the migration and testing. > > -- > Billy K. Poon > Research Scientist, Molecular Biophysics and Integrated Bioimaging > Lawrence Berkeley National Laboratory > 1 Cyclotron Road, M/S 33R0345 > Berkeley, CA 94720 > Tel: (510) 486-5709 > Fax: (510) 486-5909 > Web: https://phenix-online.org > > On Thu, Jan 7, 2016 at 9:30 AM, Nicholas Sauter <nksauter(a)lbl.gov> wrote: > >> David, >> >> I notice that the Pango version, 1.16.1, was released in 2007, so perhaps >> it is no surprise that the latest Ubuntu does not support it. Maybe this >> calls for stepping forward the Pango version until you find one that works. >> I see that the latest stable release is 1.39. >> >> This would be valuable information for us..Billy Poon in the Phenix group >> is supporting the Phenix GUI, so it might be advisable for him to update >> the Pango version in the base installer. >> >> Nick >> >> Nicholas K. Sauter, Ph. D. >> Computer Staff Scientist, Molecular Biophysics and Integrated Bioimaging >> Division >> Lawrence Berkeley National Laboratory >> 1 Cyclotron Rd., Bldg. 33R0345 >> Berkeley, CA 94720 >> (510) 486-5713 >> >> On Thu, Jan 7, 2016 at 8:54 AM, David Waterman <dgwaterman(a)gmail.com> >> wrote: >> >>> Hi again >>> >>> Another data point: I just tried this on a different Ubuntu machine, >>> this time running 14.04. In this case pango installed just fine. In fact >>> all other packages installed too and the machine is now compiling cctbx. >>> >>> I might have enough for comparison between the potentially working 14.04 >>> and failed 15.04 builds to figure out what is wrong in the second case. >>> >>> Cheers >>> >>> -- David >>> >>> On 7 January 2016 at 09:56, David Waterman <dgwaterman(a)gmail.com> wrote: >>> >>>> Hi folks >>>> >>>> I recently tried building cctbx+dials on Ubuntu 15.04 following the >>>> instructions here: >>>> http://dials.github.io/documentation/installation_developer.html >>>> >>>> This failed during installation of pango-1.16.1. Looking >>>> at pango_install_log, I see the command that failed was as follows: >>>> >>>> gcc -DHAVE_CONFIG_H -I. -I. -I../.. >>>> -DSYSCONFDIR=\"/home/fcx32934/sw/dials_bootstrap_test/base/etc\" >>>> -DLIBDIR=\"/home/fcx32934/sw/dials_bootstrap_test/base/lib\" >>>> -DG_DISABLE_CAST_CHECKS -I../.. -DG_DISABLE_DEPRECATED >>>> -I/home/fcx32934/sw/dials_bootstrap_test/base/include >>>> -I/home/fcx32934/sw/dials_bootstrap_test/base/include/freetype2 -g -O2 >>>> -Wall -MT fribidi.lo -MD -MP -MF .deps/fribidi.Tpo -c fribidi.c -fPIC >>>> -DPIC -o .libs/fribidi.o >>>> In file included from fribidi.h:31:0, >>>> from fribidi.c:28: >>>> fribidi_config.h:1:18: fatal error: glib.h: No such file or directory >>>> >>>> The file glib.h appears to be in base/include/glib-2.0/, however this >>>> directory was not explicitly included in the command above, only its >>>> parent. This suggests a configuration failure in pango to me. Taking a look >>>> at base_tmp/pango-1.16.1/config.log, I see what look like the relevant >>>> lines: >>>> >>>> configure:22227: checking for GLIB >>>> configure:22235: $PKG_CONFIG --exists --print-errors "$GLIB_MODULES" >>>> configure:22238: $? = 0 >>>> configure:22253: $PKG_CONFIG --exists --print-errors "$GLIB_MODULES" >>>> configure:22256: $? = 0 >>>> configure:22304: result: yes >>>> >>>> but this doesn't tell me very much. Does anyone have any suggestions as >>>> to how I might proceed? >>>> >>>> Many thanks >>>> >>>> -- David >>>> >>> >>> >>> _______________________________________________ >>> cctbxbb mailing list >>> cctbxbb(a)phenix-online.org >>> http://phenix-online.org/mailman/listinfo/cctbxbb >>> >>> >> >> _______________________________________________ >> cctbxbb mailing list >> cctbxbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/cctbxbb >> >> > > _______________________________________________ > cctbxbb mailing list > cctbxbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/cctbxbb > >

9 years

Search results for query "look through"