Re: [PyMOL] coordination bonds

2010-05-23 Thread Nathaniel Echols
On Sun, May 23, 2010 at 6:33 PM, Maia Cherney  wrote:

> Thank you for your help. However, I have 4  Mg ions  in the au, so all Mg
> ions (and all coordinated ligands) got connected. I need only coordination
> for one Mg.
>

Then you need to qualify the selection with the residue number, and repeat
up to 4 times.  Same principle still applies, though:

dist (elem mg and resid 341), (all within 2.2 of (elem mg and resid 341))

this particular example may be unnecessarily verbose if the residue number
is unique.

-Nat
--

___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] coordination bonds

2010-05-23 Thread Nathaniel Echols
On Sun, May 23, 2010 at 3:55 PM, Maia Cherney  wrote:

> I would like to show coordination bonds between Mg ion and its ligands
> (oxygens) as broken lines, the distance of Mg-O bond is between 1.9-2.2 A.
>

dist elem mg, (all within 2.2 of elem mg)

Might need to try a slightly longer length to catch everything.

-Nat
--

___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] for evaluation only

2010-04-29 Thread Nathaniel Echols
On Thu, Apr 29, 2010 at 3:03 PM, victor kenyon  wrote:

> i am an academic using pymol and the new build places a "for
> evaluation only" graphic on top of the rendered structure. not very
> helpful for presentations, publications and the like. how do i
> circumvent this without using chimera to generate my images for pubs.
>

Use the open-source version, I would be very surprised if that overlays
anything.   Or download one of the builds of version 0.99, which is missing
some nifty features but still more than adequate for publication graphics as
long as you don't need to show anisotropic Bs.

http://www.pymol.org/rel/099/
--
___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] Continously check the file of a molecule, and do action when it gets modified

2010-03-28 Thread Nathaniel Echols
On Fri, Mar 26, 2010 at 5:14 PM, Werner Schroedinger <
werner.schroedin...@googlemail.com> wrote:

> For making this procedure faster and more automaticm I wonder if somehow
> can I tell pymol to pay atention to some external file, or for the
> ligand being visualized now, to its copy on the hard disk, and
> continuosly monitor it, so when it detects a change (for instance,
> beacuse the quantum program did something) reload the molecule and
> display it.
>

There may already be a built-in function to do this in PyMOL, but I'm not
aware of it.  It's easy to script in Python, however.  The code below
definitely works on Mac.  Depending on how fast your calculation is, you may
want to increase or decrease the interval, but 2s is a reasonable place to
start.

-Nat


from pymol import cmd
import threading
import os.path
import time

class monitor_file (object) :
  def __init__ (self, file_name, object_name, interval=2, verbose=True) :
self.file_name = file_name
self.object_name = object_name
self.interval = interval
self.abort = False
self.last_mtime = 0
if os.path.isfile(file_name) :
  self.reload_file()
  self.last_mtime = os.path.getmtime(file_name)

  def reload_file (self) :
if self.verbose :
  print "Reloading %s as '%s'" % (self.file_name, self.object_name)
cmd.load(self.file_name, self.object_name, state=1)

  def start (self) :
if self.verbose :
  print "starting thread"
t = threading.Thread(target=self.run)
t.start()

  def run (self) :
while not self.abort :
  if os.path.isfile(self.file_name) :
if self.verbose :
  print "checking file mtime. . ."
current_mtime = os.path.getmtime(self.file_name)
if current_mtime > self.last_mtime :
  self.reload_file()
  self.last_mtime = current_mtime
  time.sleep(self.interval)
return True

  def stop (self) :
self.abort = True

# Usage
m = monitor_file("/var/tmp/caffeine.pdb", "caffeine")
m.start()
--
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

[PyMOL] recentering and zooming on xyz coordinates

2010-03-25 Thread Nathaniel Echols
Is there a way to do this in such a way that the rest of the view matrix
changes appropriately, including clipping planes, etc.?  There is a handy
pair of functions in Coot that do this:

set_rotation_centre(x, y, z)
set_zoom(30)

which has the advantage of adjusting the perspective and focusing on exactly
the residue/molecule of interest.  I tried to emulate the first part in
PyMOL:

  def recenter_and_zoom (self, x, y, z) :
view = list(cmd.get_view())
view[12] = x
view[13] = y
view[14] = z
cmd.set_view(view)

which sets the center of rotation, but the perspective is screwy because the
rest of the view matrix is unchanged.  I can figure out the equivalent
'zoom' command for the residue of interest if I really need to, but I
already have the coordinates computed, so I'd rather try to figure out how
to adjust the zoom level through the matrix instead.  (Sorry if this is a
little incoherent - I still don't understand matrix arithmetic in OpenGL.)

thanks,
Nat
--
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] getting Residue environment distances and residue-residue distances programatically

2010-03-23 Thread Nathaniel Echols
On Tue, Mar 23, 2010 at 12:48 PM, hari jayaram  wrote:

> Hi
> I want to write a python/pymol  script that will give me residue -residue
> distances within a pdb file programmatically
>
> I know that within pymol I can get a distance object using
>
> select r55 , chain A and resi 55
> select r 22 , chain A and resi 22
>
> distance (r55) ,(r22)
>
> This creates the dist01 object . If I wanted to do this using a python
> script and get a pretty print of all the distances in the distance object ,
> how do I do that?
> Also pointers to other python toolkits that can allow me to make such
> measurements outside pymol will be greatly appreciated.
>

CCTBX will do this for sure - it supports a selection syntax very similar to
PyMOL's.  It isn't particularly intuitive, however; if you're interested in
using it, contact me offline and I can send you some examples.

There is also BioPython, which is written for bioinformatics rather than
crystallography and appears to have a much simpler API (but will definitely
read PDB files).

-Nat
--
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] PyMOL-users Digest, Vol 45, Issue 5

2010-02-15 Thread Nathaniel Echols
On Sun, Feb 14, 2010 at 4:39 PM, David Hall  wrote:

> I don't think DynoPlot can't do this.  Are there even callbacks to
> watch for things like changing state in the Pymol GUI?  Maybe that
> should be a feature request as it would allow for a lot more powerful
> plugins.  Callbacks that we can bind to through TKinter to watch for
> modifications in the OpenGL window of Pymol like changing state,
> moving atoms, etc, then we can easily build boxes that reflect these
> changes.  If this is somehow already possible, I would love someone to
> point me to an example.


Check out pymol.importing.load_callback - is that what you were thinking of?
 (I think I've tried this before, but I can't remember exactly what the
behavior was.)

-Nat
--
SOLARIS 10 is the OS for Data Centers - provides features such as DTrace,
Predictive Self Healing and Award Winning ZFS. Get Solaris 10 NOW
http://p.sf.net/sfu/solaris-dev2dev___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] color by element

2010-02-09 Thread Nathaniel Echols
On Tue, Feb 9, 2010 at 2:46 PM, Ariel Talavera  wrote:

> How can I use the color by element (HNOS) function present in the GUI,
> in a script?
>

I don't think there's a single command that does this; I've always done
something like this:

color grey80, elem c
color white, elem h
color red, elem o
color blue, elem n
color orange, elem s

and so on (in python, these become 'cmd.color("grey80", "elem c")', etc.).
 I just discovered that there is also a separate color name for each
element, e.g. "color carbon, elem c" will work, but the default coloring
scheme seems a bit wacky sometimes (carbon is green, etc.).  If there was a
convenient way to switch between available color sets (e.g. all of the
options available in the menus), that would be awesome.  Currently the
atom-name colors are hardcoded in C, unfortunately.

Nat
--
SOLARIS 10 is the OS for Data Centers - provides features such as DTrace,
Predictive Self Healing and Award Winning ZFS. Get Solaris 10 NOW
http://p.sf.net/sfu/solaris-dev2dev___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] New Feature Ideas

2010-01-21 Thread Nathaniel Echols
On Thu, Jan 21, 2010 at 6:10 AM, Schubert, Carsten [PRDUS] <
cschu...@its.jnj.com> wrote:

> 4) A better density wizard, let's just copy coot and be done with it.
> Ability to dynamically bind density levels or some other properties to
> the scroll-wheel for that matter.
>

+1


> 6) Integration/bundling of wxPython as an alternative to TCL/Tk??
>

Also +1.  Warren mentioned this last year as a possibility for integrating
PyMOL with the Phenix GUI, and in theory support for embedding a PyMOL view
in wxPython is already present.  However, that code appears to be very old
and not maintained - I couldn't get it to work even after updating it to use
the newer wx modules.  I suspect there is not much work to do in order to
resurrect it, but I don't know enough about PyMOL internals, or have time to
learn.  Note that this feature (if it were open-source) would allow anyone
to make a native Mac PyMOL GUI and distribute that, so it isn't something
I'd ask Schrodinger to subsidize.  (If anyone else is interested in helping
figure this out, however, let me know.)

7) Internal FFT routines to be able to read map-coefficients
>

+2.  I believe the CCP4 libraries for this are open-source, so those are one
option.  CCTBX is another (I have tried this and it works), but this will be
a lot more overhead.
--
Throughout its 18-year history, RSA Conference consistently attracts the
world's best and brightest in the field, creating opportunities for Conference
attendees to learn about information security's most important issues through
interactions with peers, luminaries and emerging and established companies.
http://p.sf.net/sfu/rsaconf-dev2dev___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] Preserving sec. structure in states for animation

2010-01-12 Thread Nathaniel Echols
On Tue, Jan 12, 2010 at 11:13 AM, Schubert, Carsten [PRDUS] <
cschu...@its.jnj.com> wrote:

>  I am trying to make a movie form a set of morphed pdb files. Currently
> the individual files are loaded into the same object with increasing states,
> all pdb files contain the SHEET and HELIX record created by dssp.  The
> object is then displayed as a cartoon. When I loop over the states I
> noticed that Pymol only seems to take secondary structure information from
> the first pdb file into account. For instance I have a helix which is split
> into 2 during the morph, but only displayed as one helix during the
> animation, albeit a bit stretched. Loading the individual files produces the
> right sec. structure. Any idea how I can convince Pymol to update the
> cartoon representation with each frame?
>

I don't think there's any way to do this that doesn't involve some
scripting.  My experience may be out of date, since I last tried this around
2002, but I don't think the movie API has changed much since then.  I
ultimately rendered each frame separately (i.e. different object for each
PDB file) using the same view matrix.  However, depending on how finely you
want to show the changes in secondary structure, something like the
following should work:

load frame1.pdb, morph
. . .
load frame20.pdb, morph
mset 1 -10
mpng movie, first=1
alter morph and resi 100:120, ss='H'
mset 11 -20
mpng movie, first=11

This is just off the top of my head - I can't guarantee that it will work
exactly as written, but it's pretty close to what you need.  In theory, it
will produce 20 frames, movie0001.png through movie0020.png, which should
appear continuous.

-Nat
--
This SF.Net email is sponsored by the Verizon Developer Community
Take advantage of Verizon's best-in-class app development support
A streamlined, 14 day to market process makes app distribution fast and easy
Join now and get one step closer to millions of Verizon customers
http://p.sf.net/sfu/verizon-dev2dev ___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] Electron density grid

2009-12-18 Thread Nathaniel Echols
On Thu, Dec 17, 2009 at 10:54 AM, Jason C Porta  wrote:

> I have been trying to modify the grid spacing on 2Fo-Fc map that was loaded
> into PyMol directly from phenix.refine. The mesh object is generated from
> the isomesh command. Basically, I would like there to be fewer grid points
> so that the density isn't so thick. If anyone has any advice, it would be
> greatly appreciated.


The grid spacing is built in to XPLOR maps.  The only workaround I can
suggest is to start from the map coefficients in MTZ format, and run CCP4
FFT to output a CCP4 format map with the grid spacing you prefer.
 (0.33*resolution is the default in most programs, but phenix.refine uses
0.25*resolution.)  I'll add a program to the Phenix GUI to do this at some
point.  For now, I'll see if I can make it easier to change the resolution
factor from the GUI.
--
This SF.Net email is sponsored by the Verizon Developer Community
Take advantage of Verizon's best-in-class app development support
A streamlined, 14 day to market process makes app distribution fast and easy
Join now and get one step closer to millions of Verizon customers
http://p.sf.net/sfu/verizon-dev2dev ___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] Alignment Objects

2009-09-11 Thread Nathaniel Echols
On Fri, Sep 11, 2009 at 1:59 PM, Maia Cherney  wrote:

> Is it posible in pymol to select all residues of a certain type
> automatically, like all arginines or all lysines + arginines in a chain?


select resnam arg
select chain A and (resnam lys or resnam arg)
--
Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day 
trial. Simplify your report design, integration and deployment - and focus on 
what you do best, core application coding. Discover what's new with 
Crystal Reports now.  http://p.sf.net/sfu/bobj-july___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] about secondary structure assignement: help

2009-07-21 Thread Nathaniel Echols
On Tue, Jul 21, 2009 at 12:23 PM, Paola Beassoni wrote:

> Dear Pymol Users,
> i am working with some pdb files with no secondary structures assignments
> embedded on it. Then when i try to see my molecules on cartoon diagrams, i
> see the secondary structure different to another programs, VMD for example.
> I used the command alter (for example alter 10-34/, ss='S'), and i am
> capable of assign the secondary structure that i want but this is just
> temporary. If i save the molecule, when i oppen it again it has the original
> secondary structure and no the altered one.
> Can anyone help me about this??? I really prefer pymol instead other
> programs for visualization!


The most compatible solution is probably to generate the HELIX and SHEET
records for your PDB file, which PyMOL will read (not sure about
VMD).  There is an open-source program called KSDSSP that will do this:

http://www.cgl.ucsf.edu/Overview/software.html#ksdssp

I haven't used it very much, but it's supposedly a re-implementation of the
original DSSP algorithm, so theoretically very robust.  If you're
continually generating new PDB files (e.g. refining the structure, etc.)
you'll probably need to keep these records in a separate file and paste them
into every new PDB file, since most programs won't retain that information
when writing files.  If you don't like the results, it's fairly simple to
modify the records to suit your taste.

Alternately, you can save a PyMOL script containing the alter commands and
just invoke that every time you load your structure.

-Nat
--
___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] Homology Modeling in Pymol

2009-06-30 Thread Nathaniel Echols
On Tue, Jun 30, 2009 at 2:24 PM, Buz Barstow  wrote:

> I have a very large library (almost 1000) of similar genetic
> sequences, for which I would like to generate homology models. Does
> anyone know of a way to automate requests to a homology modeling
> server, such as SWISS-Model through pymol, or alternatively through
> python?


I believe Modeller uses Python - never tried it myself, though.  There is a
server for Modeller but academic users can download the program itself.

There is in fact an API for using SWISS-Model from a Perl script (which
could probably be converted to Python, or wrapped by a Python script), which
I used years ago for a similar project.  As far as I know this is
unpublished, probably because they don't want people to trash their server,
but they were happy to supply it when I contacted the site maintainer.  They
asked that I limit the frequency of my requests, however.  I would recommend
emailing them.  (I have no idea how Modeller compares to SWISS-Model, but I
suspect it's more rigorous, especially when sequence homology is remote.)

There are other ways to go about this - bioinformaticists have been abusing
online databases and servers like this for years - but it's considered bad
etiquette to do this without permission (and some groups explicitly forbid
this use).

-Nat
--
___
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Re: [PyMOL] Newbie attempting to python script: ImportError: No module named _cmd

2009-05-01 Thread Nathaniel Echols
On Fri, May 1, 2009 at 3:57 AM, Luke Goodsell wrote:

> I'm trying to set up my machine to be able to run python scripts using
> the PyMol API, but whenever I try to import the pymol module, I get the
> following output:
> . . .
>  I am running Mac OS X 10.5.6. I downloaded MacPyMol 1.1r1 from
> http://delsci.com/macpymol/ (the educational edition, as I am a student).
> My PYMOL_PATH: /Applications/PyMOLX11Hybrid.app/pymol
> My PYMOL_EXE: /Applications/PyMOLX11Hybrid.app/Contents/MacOS/MacPyMOL
> . . .
>  I copied the contents of /Applications/PyMOLX11Hybrid.app/pymol/modules
> to /System/Library/Frameworks/Python.framework/Versions/2.5/lib/pymol,
> as per my understanding of the top of __init__.py
>

It's not entirely clear to me what you're trying to do, but I'm pretty
certain it's guaranteed not to work.  You can't just take compiled objects
from inside app bundles and move them around - especially on Mac.  You
definitely shouldn't be installing PyMOL modules that way.  If you need them
available for importing by /usr/bin/python, download the open-source version
from SourceForge and install it using python distutils (i.e. "python
setup.py build" and so on).  I just tried this and it seems to work fine -
only takes about 10 minutes.

-Nat
--
Register Now & Save for Velocity, the Web Performance & Operations 
Conference from O'Reilly Media. Velocity features a full day of 
expert-led, hands-on workshops and two days of sessions from industry 
leaders in dedicated Performance & Operations tracks. Use code vel09scf 
and Save an extra 15% before 5/3. http://p.sf.net/sfu/velocityconf___
PyMOL-users mailing list
PyMOL-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pymol-users


Re: [PyMOL] Roving detail and Roving density

2009-04-23 Thread Nathaniel Echols
On Wed, Apr 22, 2009 at 6:37 PM, Craig Smith  wrote:

> I was looking at the demos in pymol (wizard --> demo) and saw
> interesting features called roving detail and roving density.  I
> looked in the manual on how to do this feature but couldn't find it.
> Could anyone point me in the right direction or show me how to
> interactively display electron density as I move through my structure.
>

This is what I do (sorry, it's in Python, but easily converted to PyMOL
commands).  It assumes that you have two maps loaded, named "2fofc" and
"fofc".(Some of it is cribbed from the C code for the roving density
demo.)

cmd.set("roving_detail", 20)
cmd.set("roving_isomesh", 20)
cmd.set("roving_origin", 1)
cmd.set("roving_sticks", 0)
cmd.set("roving_ribbon", 0)
cmd.set("roving_lines", 0)
cmd.set("roving_spheres", 0)
cmd.set("roving_nb_spheres", 0)
cmd.set("roving_polar_contacts", 0)
cmd.set("roving_polar_cutoff", 0)
cmd.set("roving_map1_name", "2fofc")
cmd.set("roving_map1_level", 1)
cmd.set("roving_map2_name", "fofc")
cmd.set("roving_map3_name", "fofc")
cmd.set("roving_map2_level", 3)
cmd.set("roving_map3_level", -3)

cmd.isomesh("rov_m1", "2fofc", 1.0, "center", 20)
cmd.isomesh("rov_m2", "fofc", 3.0, "center", 20)
cmd.isomesh("rov_m3", "fofc", -3.0, "center", 20)
cmd.color("skyblue", "rov_m1")
cmd.color("green", "rov_m2")
cmd.color("red", "rov_m3")

This might not be ideal - some of those 'set' commands may be unnecessary or
incorrect.  However, the end result works more or less like Coot, which is
what I was aiming for.

-Nat
--
Stay on top of everything new and different, both inside and 
around Java (TM) technology - register by April 22, and save
$200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
300 plus technical and hands-on sessions. Register today. 
Use priority code J9JMT32. http://p.sf.net/sfu/p___
PyMOL-users mailing list
PyMOL-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pymol-users


[PyMOL] Roving density

2008-12-31 Thread Nathaniel Echols

Hi--

Can the roving mesh be used without changing the model appearance?   
I'd prefer to just show the entire structure as lines only, but it  
looks like the roving_detail setting controls both model and mesh.


thanks,
Nat



Re: [PyMOL] Ray Tracing A Protein Gallery

2008-12-15 Thread Nathaniel Echols
On Mon, Dec 15, 2008 at 11:41 AM, Buz Barstow  wrote:

> I'm making a gallery of protein molecules for my PhD thesis. I'd like
> to find an automatic way to ensure that all of the ray traced images
> have the same scale. Is there an easy way to do this?
>

Translate every model so that the center is at the origin, then use the
set_view command to apply the same viewing matrix for each model.  (You can
still rotate the models, just don't zoom.)  I'm not sure if there's already
a simple way to do the translation from within PyMOL, but it should be quite
straightforward using cctbx or CNS or something similar.

-Nat


Re: [PyMOL] remove cartesian coordinate from w/in transparent metal ions

2008-10-29 Thread Nathaniel Echols
On Wed, Oct 29, 2008 at 11:39 AM, Thomas S. Leyh, Ph. D.
wrote:

>  Every once in a while, I try to show a metal ion as a tranparent
> sphere and inevetiably get stuck trying to make the tiny cartesian
> coordinate system that appears at the center of the sphere disappear so it
> doesn't "show through".  Suggestions?
>

I think you want "hide nonbonded".

I have a related question: is there a way (or could there be) to control the
default appearance of molecules upon loading?  Personally, when I use PyMOL
I almost always want to look at cartoons, without nonbonded atoms.
 Something as simple as a few new settings that I could throw in ~/.pymolrc
would be fine.


Re: [PyMOL] displaying anisotropy ellipsoids in PyMOL v0.99rc6 running on Windows XP

2008-09-22 Thread Nathaniel Echols
Tim Fenn wrote something to do this a few years ago, but I don't think you
can ray-trace them:
http://www.stanford.edu/~fenn/pdb_aniso.py

On Mon, Sep 22, 2008 at 2:48 PM, "Javier M. González"
wrote:

> Hi all, could someone advise me on the simplest way to display anisotropy
> ellipsoids in PyMOL v0.99rc6 running on Windows XP?
> Is there any script to do this?
> I haven't moved to the latest versions of PyMOL yet.
>
> Any help would be appreciated. Thanks in advance.
>
> Javier.
>
> --
> Lic. Javier M. González
> Instituto de Biología Molecular y Celular de Rosario, IBR-CONICET
> Universidad Nacional Rosario, Facultad de Ciencias Bioquímicas y
> Farmacéuticas, Area Biofísica
> Suipacha 531 CP S2002LRK
> Tel: +54 341 4350661 / 4350596 / 4351235 Ext.109
> FAX: +54 341 4390465
> Rosario, Santa Fe, ARGENTINA
> gonza...@ibr.gov.ar
>
>
> -
> This SF.Net email is sponsored by the Moblin Your Move Developer's
> challenge
> Build the coolest Linux based applications with Moblin SDK & win great
> prizes
> Grand prize is a trip for two to an Open Source event anywhere in the world
> http://moblin-contest.org/redirect.php?banner_id=100&url=/
> ___
> PyMOL-users mailing list
> PyMOL-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/pymol-users
>
>


Re: [PyMOL] export 2 molecules to one pdb file?

2008-08-13 Thread Nathaniel Echols
Just use the (PyMOL) command line: "save complex.pdb" will save all objects.
 You'll need to figure out what directory you're in first - by default, this
is probably /Applications, but "cd ~" will get you home.  (Hey Warren, could
this be changed so "cd" alone works, like in a Unix shell?)
(Or if you have more than just those two molecules: "save complex.pdb,
object1 or object2")

On Wed, Aug 13, 2008 at 12:01 PM, Paul Shannon
wrote:

> Thanks to Warren for setting me straight on independent positioning of
> molecules.
>
> Now I wish to write out the two molecules, together, to one pdb file,
> as required by RosettaDock.  I must be blind.  They say:
>
>...from the main menu, File→Export Molecule can be used to write
> a PDB file containing
>the starting structure with both docking partners.
>
> The closest menu option I can find (using MacPyMol) is File->Save
> Molecule...
> The resulting dialog apparently only allows one object to be written.
>
> Any advice?  Thanks!
>
>  - Paul
> -
> This SF.Net email is sponsored by the Moblin Your Move Developer's
> challenge
> Build the coolest Linux based applications with Moblin SDK & win great
> prizes
> Grand prize is a trip for two to an Open Source event anywhere in the world
> http://moblin-contest.org/redirect.php?banner_id=100&url=/
> ___
> PyMOL-users mailing list
> PyMOL-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/pymol-users
>


Re: [PyMOL] Mac OS X install, general PyMOL function - badly need overview

2002-09-01 Thread Nathaniel Echols
>   I have downloaded the version (0.79.dmg.gz) of the PyMol package 
> suggested by Luca Jovine, but don't know how to get a working version of 
> Python installed.  The only links I can find are to Python .exe files 
> for versions 2.1 and 1.5.2, and I have no idea which of the other 
> downloads (xx-contrib, ext) are needed.  For that matter, is the v. 0.79 

It's been a while since I played with the Mac version (0.79), but I recall
it not needing any additional installations.  Python is pretty tightly
integrated into PyMOL in the "default" setup, from what I can tell.
Anyway, the installed application should run with a simple double-click,
though I ended up invoking it from a shell.  I imagine the non-native
version would be considerably more complicated to set up, but I haven't
tried it.

> tutorials back in May - did anything come of that?  A basic intro to 
> PyMOL is what I, a confirmed nonprogrammer (and damned proud of it) 
> need, no doubt like many others.  Maybe some academics have already put 

I've found the manual does an adequate job of introduction, even if it's
very outdated.  My only real complaint is that it barely touches on the
subject of programming with/for PyMOL (I'm a programmer, not a
crystallographer).

There's also this:

http://pymol.sourceforge.net/tutorial/index.html

It's geared towards Windows users, but this shouldn't make too much of a
difference.

-Nat






Re: [PyMOL] Secondary structure assignments

2002-08-20 Thread Nathaniel Echols
> The best way to avoid "derivative work" concern is to code in a "clean room"
> fashion, without reference to someone else's code.  In principle, this could
> be done working from the original paper(s) that describe an algorithm.

I assume this is what authors of other packages have done in the past.
If there's enough interest, and no one is working on this already, I'd be
willing to take a stab at it.

> Not that I'm volunteering ;-)  Just that I would hate to see someone who was
> willing and interested in implementing secondary structure code for PyMol to
> be frightened off the job if, indeed, it could be done in a way that respects
> others' rights.

The other thing I'd find useful is some molecular mechanics tools,
particularly energy minimization.  Blending PyMOL with MMTK seems like a
good way to do this, since MMTK's license looks fairly liberal.
Personally, I'd have to spend two weeks looking at the code for both
packages before I could even start thinking about how this would need to
be done.  I'm certainly not capable of implementing anything like that
from scratch- I slept through math class in college.

-Nat Echols




[PyMOL] clipping planes

2002-08-06 Thread Nathaniel Echols

I don't understand how these work at all.  Is there any straightforward 
way, given a structure's position and dimensions and the desired camera
position, to calculate the ideal position for the plane?   How does
PyMOL decide where to put the planes?  Right now I'm missing by a small
amount and atoms disappear randomly.

thanks,
Nat




Re: [PyMOL] Secondary structure assignments

2002-07-29 Thread Nathaniel Echols
> Ideally util.ss could become a shell that would know (from the environment?) 
> if dssp or similar was installed, then either run dssp and parse the output 
> to assign the ss, or else run the ss asignment code in the current util.ss if 
> no other options were available.

I've basically done this, though right now it uses STRIDE because the
output is easier to parse- it just uses popen2() to read the output from
the STRIDE executable (Does anyone actually know what the difference is
between these two programs?). But since I already wrote a dssp parser in
Perl, it shouldn't take long to modify.  My feeling is that it would be
better to have a robust secondary structure algorithm built in, since I've
already hit limitations in the method I'm using.  Though to be honest,
I've been getting away with using util.ss for some time now.

> The ability to write a pdb with the newly generated ss information would be 
> useful as well (or perhaps pymol already can do this if it knows the 
> protein's ss?).

I've done this too, though I think it breaks with more than one chain
right now.  If I get around to cleaning this up I'd be more than happy to
pass it on.

-Nat




Re: [PyMOL] Secondary structure assignments

2002-07-29 Thread Nathaniel Echols
> Yes, please. There are some structures that have to be generated by using
> symmetry operations defined in the "REMARK 350" lines in the pdb file 
> (example:
> 1jjy.pdb). The only way (I know of, if you know a better one, please tell me) 
> to
> do this is going to http://pqs.ebi.ac.uk where the complete structure is
> generated for you. However, the caveat is: This service removes the secondary
> structure information, and up to now I have manually copy-pasted it into the 
> new
> pdb file (Annoying work, really).

Unfortunately, the script I came up with won't deal with that; it
basically uses a third-party program to rescan the input file and uses the
alter command to correctly set the structure (or combines this with file
loading as a single command). If anyone would like this, I'm happy to
email it to them, but it's not really a good general solution.

How hard would it be to incorporate the DSSP algorithm?  I was looking at
the source for RasMol, which has a complete implementation.  My feeling is
that this could be ported without too much misery.  However, RasMol's
license terms are different, and I suspect even a complete translation
into Python with the ChemPy API  would still be considered a "derivative
work", meaning it wouldn't be able to be part of the main distribution,
and meaning someone would still have to code DSSP from scratch.  (I have
utterly no idea how it works, but RasMol's code seems clear enough despite
lack of comments)

-Nat




Re: [PyMOL] Secondary structure assignments

2002-07-28 Thread Nathaniel Echols
> programs.  So I am in search of a solution.  Can one manually create a
> secondary structure assignment matrix to insert into the PDB file?  I must
> admit that as a novice I do not understand the formatting of the secondary
> structure matrices that I have seen embedded  in some PDB files.  Is there a
> resource somewhere that explains the formatting of such a matrix?

I've been fiddling with this for my own work; I wrote a parser in Perl to
embed structural annotation in a PDB file based on output from DSSP,
STRIDE, or RasMol.  Unfortunately it doesn't handle multiple chains yet,
but I'm working on making it more general.  I suppose since I only use
this for loading files into PyMOL, I could just write it in Python as a
PyMOL module. . . would that be useful to anyone else?

The complete PDB specification for SS annotation is here:
http://www.rcsb.org/pdb/docs/format/pdbguide2.2/part_42.html
http://www.rcsb.org/pdb/docs/format/pdbguide2.2/part_44.html
http://www.rcsb.org/pdb/docs/format/pdbguide2.2/part_45.html

However, PyMOL's built-in commands end up being much easier to use... I
wish I'd realized that before I wrote my parser.

> Alternatively, is there a manual series of Pymol commands to directly assign
> secondary structure for a protein for insertion into the .pml script file?

alter atom_selection, ss='L' (or 'H' or 'S', I think)

> Even better, is there a separate program that can assign secondary structure
> and then allow you to save the revised PDB file with the secondary structure
> assignments inserted into the PDB output.  I would be willing to accept any

If there is, I'd like to know too!

-Nat Echols
Gerstein Lab
Yale University




Re: [PyMOL] Suppress some output

2002-07-25 Thread Nathaniel Echols
> Is there a way to keep my program 'quiet',
> or I can customize the output message myself?

Depending on your OS and shell, something like this may work:

denethor:~ > /usr/local/pymol/pymol.com >/dev/null 2>/dev/null

This works under Bash in any modern Linux distribution, from what I've
seen.  I suspect it would work on Irix too (still using Bash).  I've had
the exact same problem- mainly with util.rainbow().

-Nat




Re: [PyMOL] pymol without GL

2002-07-16 Thread Nathaniel Echols
> pymol.com -c scriptname
> But you have to get your view first somehow.

I can't do this until I compile it; if I try running a pre-existing build
I get errors because one of the GL libraries isn't found.  I've ended up
just installing GL on my system, but the crack technical wizards at SuSE
have decided I need Qt and KDE as well. . . on web server. *sigh*  It's
not really important that I be able to do this, and there's plenty of disk
space to spare, but I prefer limited installs on a server.

By the way, has anyone else experimented with hardware antialiasing on
NVidia cards under Linux?  I've had mixed results so far, and I'm not sure
what I'm doing wrong (this is with a GeForce 4 Ti 4600).

-Nat





[PyMOL] pymol without GL

2002-07-16 Thread Nathaniel Echols
quick question: I'm trying to do batch jobs on a server.  Obviously,
there's no X or OpenGL installed.  Is there any way to build PyMOL without
any graphical frontend, and just use it as a raytracer?  I don't want to
go through the trouble of adding in these libraries...

thanks,
Nat




[PyMOL] view area

2002-07-07 Thread Nathaniel Echols
Hi--

I'm trying to render a large number of images of a trajectory
automatically. I have two problems:

1. There doesn't seem to be any way to manually speciify the view area.
Thus, if I want all images to be in the same frame of reference, I need to
load them all as different states of the same object, e.g. "load
frame01.pdb, mov" many times over.  This works fine, but even with 4GB of
memory PyMOL crashes when I load everything in.  I've found a way around
this, but it would be much simpler if I could specifiy exactly what angle
and window to use for rendering.  I'm pretty sure this doesn't exist right
now, but will it in the future?

2. Using pretty much any method, PyMOL frequently lops off part of the
molecule.  This has become a large problem; I'm using PyMOL quite a bit to
make animations (replacing Molscript entirely), and frequently an
interesting part of the structure will be out of frame.  How can I avoid
this?  This is all being done on the command line (usually invoked by some
CGI script), so I have no way of adjusting the zoom.  'zoom' hasn't been
too much help.

Any advice would be much appreciated.  I'm trying to do stuff that would
be messy in pretty much any environment (e.g. I can't even load all the
structures with less than 1GB of memory), but the viewport problems apply
to pretty much anything.

thanks,
Nat