Re: [ccp4bb] [ccpem] Making thioester bond in Chimera/ChimeraX

2023-11-09 Thread Eric Pettersen
Hello Firdous,
Probably the problem is that mol A and mol B are two different models.  You 
can't add a bond in Chimera or ChimeraX until after combining them into one 
model first.

ChimeraX combine command:
<https://rbvi.ucsf.edu/chimerax/docs/user/commands/combine.html>

ChimeraX bond command or Build Structure tool, Adjust Bonds section:
<https://rbvi.ucsf.edu/chimerax/docs/user/commands/bond.html>
<https://rbvi.ucsf.edu/chimerax/docs/user/tools/buildstructure.html#bonds>

For the best chance of getting an answer to a question about ChimeraX, please 
use the chimerax-us...@cgl.ucsf.edu <mailto:chimerax-us...@cgl.ucsf.edu>address:
<https://rbvi.ucsf.edu/chimerax/docs/contact.html>

I hope this helps,

--Eric

Eric Pettersen
UCSF Computer Graphics Lab

> From: Collaborative Computational Project in Electron cryo-Microscopy 
> mailto:cc...@jiscmail.ac.uk>> on behalf of Firdous 
> Tarique mailto:kahkashantari...@gmail.com>>
> Date: Thursday, November 9, 2023 at 12:27 PM
> To: cc...@jiscmail.ac.uk <mailto:cc...@jiscmail.ac.uk>  <mailto:cc...@jiscmail.ac.uk>>
> Subject: [ccpem] Making thioester bond in Chimera/ChimeraX
> 
> Hi
>  
> Does anybody know how to make thioester bonds in Chimera or ChimeraX ?
>  
> I am unable to make a bond between Cys321 in mol A with Gly1 in mol B 
> (S-glycyl-L-cysteine) using Chimera or ChimeraX.
>  
> I tried Isolde in ChimeraX and Build Structure in Chimera to make a bond 
> between selected atoms but so far have been unsuccessful. I don't know if I 
> am missing something.
>  
> Any help would be appreciated.
>  
> Thanks
>  
> Firdous
>  



To unsubscribe from the CCP4BB list, click the following link:
https://www.jiscmail.ac.uk/cgi-bin/WA-JISC.exe?SUBED1=CCP4BB=1

This message was issued to members of www.jiscmail.ac.uk/CCP4BB, a mailing list 
hosted by www.jiscmail.ac.uk, terms & conditions are available at 
https://www.jiscmail.ac.uk/policyandsecurity/


Re: [ccp4bb] CCP4BB Digest - 13 Feb 2021 to 14 Feb 2021 (#2021-51)

2021-02-14 Thread Eric Pettersen
> On Feb 14, 2021, at 4:00 PM, CCP4BB automatic digest system 
>  wrote:
> 
> Date:Sat, 13 Feb 2021 20:24:07 -0500
> From:Nicholas Larsen  <mailto:nicholas_lar...@h3biomedicine.com>>
> Subject: Re: Bug in mmCIF handling of UNK residues?
> 
> I hope this doesn't confuse the discussion, but my understanding was "UNK"
> stood for "unknown" residue and this will cause errors.  UNK naming
> convention is the default output of Schrodinger when generating ligand PDB
> files.  Coot will display the PDB containing "UNK" as a residue, but if you
> try to use the CIF file to real-space refine, the ligand will blow up.   I
> found that renaming the residue in the output PDB and regenerating the CIF
> file with the corrected RESID name solved the problem.  So in
> my experience, the problem is the name "UNK" and this just needs to be
> switched to something else.  Has anyone else seen this?
> Nick

The PDB defines UNK as unknown amino acid.  UNL is unknown ligand.  N is 
unknown nucleic acid.

--Eric

Eric Pettersen
UCSF Computer Graphics Lab





To unsubscribe from the CCP4BB list, click the following link:
https://www.jiscmail.ac.uk/cgi-bin/WA-JISC.exe?SUBED1=CCP4BB=1

This message was issued to members of www.jiscmail.ac.uk/CCP4BB, a mailing list 
hosted by www.jiscmail.ac.uk, terms & conditions are available at 
https://www.jiscmail.ac.uk/policyandsecurity/


Re: [ccp4bb] strict structure based alignment

2012-07-15 Thread Eric Pettersen
On Jul 13, 2012, at 4:00 PM, Christian Roth wrote:

 I want align a couple or protein structures by secondary structure matching 
 to 
 one target and want get a kind of aminoacid alignment file e.g. what residue 
 fit 
 the other, without adjustments due to sequence based alignments. 
 I tried Strap, but as far as I understood it, it takes also the sequence into 
 account. I tried also Rapido, but this does only a pairwise comparison. 
 Superpose does align it nicely (ccp4 based or Coot based) but there seems to 
 be no option to print the sequence alignment in a file and it is again  just 
 a 
 pairwise comparison .
 Is there an other program which does something similar?

If you use UCSF Chimera (www.cgl.ucsf.edu/chimera), you can use the MatchMaker 
tool to superimpose the structures.  MatchMaker allows you to adjust the weight 
of sequence similarity vs. secondary structure matching, so you can just make 
the sequence similarity 0% and the secondary structure 100%.  With the 
structures superimposed, you can use the Match-Align tool to generate a 
sequence alignment based solely on the proximity of residues to one another in 
space.  Be warned that Match-Align will be very slow for 10+ structures, but 
is fine for half a dozen or so.  The generated alignment will be displayed in a 
window that will have a File menu where you can save the alignment to a 
variety of common formats.

--Eric

Eric Pettersen
UCSF Computer Graphics Lab




Re: [ccp4bb] Tool for calculating RMSD

2012-06-22 Thread Eric Pettersen
If you open the structures in Chimera, the measure rotation command will 
report the transformation matrix as well as the screw axis and shift/rotation 
along/around the screw axis needed.  There are also options to that command for 
actually depicting the screw axis and/or rotation.

www.cgl.ucsf.edu/chimera
www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/measure.html

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu

On Jun 22, 2012, at 4:00 PM, Jeremy Tame wrote:

 Yes, LSQKAB does, but not the equivalent screw axis (which exists by Chasles 
 theorem). I think a screw
 axis may well be the best way to describe the rigid body movement which began 
 the thread. A simple
 distance and angle (2 numbers) are easier to understand, and often of more 
 biological relevance, than
 RT matrices, for example in the case of relative domain motions.
 
 On Jun 20, 2012, at 5:20 PM, Mark J van Raaij wrote:
 
 LSQKAB (superpose in CCP4i GUI) also outputs in its log-file the 
 translation parameters and rotation matrix it used to superpose the 
 structure.
 




Re: [ccp4bb] residuewise rmsd for multiple chain superposition

2012-02-21 Thread Eric Pettersen

On Feb 21, 2012, at 4:00 PM, sreetama das wrote:

   Is there any software which will print the RMSDs (residue  
wise, and not chain wise) for more-than-2 chains (having similar/ 
same sequences) superposed together?


In UCSF Chimera you can use the Match-Align tool to create a  
structure-based sequence alignment (or you can open an alignment file  
if you already have one).  The the alignment will show per-residue  
RMSD as a bar chart across the top of the alignment (if you opened  
your own alignment file you'd have to use the Headers-RMSD menu entry  
to show the bar chart).  You can save the numeric RMSD values to a  
file with Headers-Save.


Also, Match-Align will report a couple of additional normalized RMSD  
scores, namely Structural Distance Measure and Q-score.


Chimera home page:  www.cgl.ucsf.edu/chimera
Match-Align page: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/matchalign/matchalign.html
sequence aignment viewer: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/multalignviewer/framemav.html

--Eric


Eric Pettersen
http://www.cgl.ucsf.edu/home/pett

And isn't sanity really just a one trick pony anyway? I mean
all you get is one trick, rational thinking, but when
you're good and crazy, oooh oooh oooh, the sky is the
limit!-The Tick




Re: [ccp4bb] Putting Text into Movies

2012-02-21 Thread Eric Pettersen

On Feb 21, 2012, at 4:00 PM, Jacob Keller wrote:


is there a good way to put text labels into movies from pymol or
otherwise? It would be helpful for conveying some ideas...


As Jason pointed out, PyMOL has some useful facilities for this.   
Chimera has a tool for placing standalone text and arrows called 2D  
Labels.  The arrows and text can be any color and the text can also be  
various sizes, styles (italics, etc.) and typefaces (serif etc.), all  
on a per-character basis.  You can place text and arrows interactively  
with the mouse or via commands.  Text can also be faded in or out as a  
movie plays/records.


A pretty nice example use of 2D Labels in the DNA basics movie on  
this page: http://crc.nd.edu/index.php/cyberinfrastructure/visualization/80


Chimera home page:  www.cgl.ucsf.edu/chimera
2D Labels page: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/2dlabels/2dlabels.html

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] secondary structure output

2012-01-24 Thread Eric Pettersen

On Jan 24, 2012, at 4:00 PM, Ed Pozharski wrote:


I am looking for a program/server that would determine secondary
structure from a pdb file and then output a new pdb file with
HELIX/SHEET records.  I have a model for which pymol fails to produce
correct secondary structure.  DSSP and STRIDE identify the secondary
structure correctly but I'd need to convert their output myself to
either pdb header or pymol script.  So I wonder maybe something like
that already exists.


If you open the PDB in Chimera and the PDB lacks HELIX/SHEET records,  
Chimera will use DSSP to compute the secondary structure and if you  
then save a new PDB file (File-Save PDB) the saved file will have the  
appropriate HELIX/SHEET records.  You can get Chimera from www.cgl.ucsf.edu/chimera 
 .


--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu


Re: [ccp4bb] Software for showing crystal packing

2011-12-05 Thread Eric Pettersen


On Dec 3, 2011, at 4:00 PM, Yuri Pompeu wrote:


Hello everyone,
Whats a good software for showing crystal packing and unit cell,  
axes , etc...
I know pymol and coot will do it but would love to hear other  
possibilities/ideas.


Chimera's Unit Cell tool can do this.

http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/unitcell/unitcell.html

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] Superpositions: Deviation by Residue

2011-11-28 Thread Eric Pettersen

On Nov 28, 2011, at 4:03 PM, Jacob Keller wrote:


Let me refine my question (sorry for my lack of clarity): is there a
program that will output the distances between the corresponding ca's
of a superposition on a residue-by-residue basis, and not just a
global RMSD value (doubtless these numbers are part of the
superposition algorithm itself)? I want to plot these values as a
function of residue number to show which parts of the structures
deviate more or less from each other.


If you use the MatchMaker tool in UCSF Chimera to make the  
superposition, it has an option to show the corresponding sequence  
alignment.  The sequence alignment will have an RMSD header running  
across the top, which is a bar graph of the RMSD values.  You can the  
alignment's Headers-Save... menu item to save the numerical values to  
a file if you want.


If you already have the structures superimposed on your own, you can  
use Chimera's Match-Align tool to create a superposition-based  
sequence alignment, and do the same thing with its RMSD header.


Some links:

Chimera home page: http://www.cgl.ucsf.edu/chimera
alignment tool's RMSD header: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/multalignviewer/multalignviewer.html#assessment

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] [CCP4] identify a rotation centre: domain rotation

2011-11-28 Thread Eric Pettersen

Perhaps this video would be helpful: 
http://www.cgl.ucsf.edu/chimera/videodoc/AlignDomains/index.html

The axis point shown in the reply log at the very end is what I  
think you want.


BTW, the video is a little out of date, in that if you get the very  
latest daily build the MatchMaker tool can match the smaller domains  
with having to delete part of the structure.


--Eric

On Nov 23, 2011, at 4:00 PM, WENHE ZHONG wrote:


Dear members,

I would like to have your ideas if there is any way to identify a  
rotation

centre of domain in two different states using CCP4 or other program.

The situation is: the domain of the protein will rotate between two
different states (depending on substrate binding) around 8 degree,  
and it
is (nearly) clearly that the domain is rotated around a rotation  
centrel.
So the question is how to identify this rotation centre in this 3D  
model?
The ideal is to identify a region of residues in the domain which  
are most

closed to the rotation centre.

The tool I am using right now is the superpose tool in CCP4  
package. The

output which I think mightbe uesful is:

CENTROID OF WORKING MOLECULE  :  157.812 152.396 -70.778
CENTROID OF WORKING MOLECULE  :(fractional)  157.812 152.396  
-70.778
CENTROID OF REFERENCE MOLECULE:  157.251 151.877  
-70.874
CENTROID OF REFERENCE MOLECULE:(fractional)  157.251 151.877  
-70.874

Distance between CENTROIDS  :0.770
Direction cosines of vector between CENTROIDS:   0.729   0.674
0.124


I would say the “CENTROID it mentioned above, such as (157.251  
151.877
-70.874), is possibly near to the rotation centre. I would like to  
have

your opinion though. Thank you.






Re: [ccp4bb] Temperature Factor statistics

2011-08-31 Thread Eric Pettersen

On Aug 31, 2011, at 4:00 PM, Yuri Pompeu wrote:

After i get my output file from baverage containing the average b- 
factor and rms by residues,

How can I calculate and display the average (and or mean) B-factors?
Is there a way of calculating it by protein, ligands and solvent  
separately?


Chimera's Render by Attribute tool can depict atoms by their B-factor  
or residues by their average B-factor (assuming you have a PDB file  
with a B-factor column).  It also shows a histogram of the B-factor  
values.  You can use the Attribute Calculator tool to compute the  
average B-factor of any part of your structure and write it to a file.


Some links:

Chimera:  www.cgl.ucsf.edu/chimera
Render by Attribute: 
www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/render/render.html
Attribute Calculator: 
www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/calculator/calculator.html
a tutorial focused on coloring by B-factor: 
www.cgl.ucsf.edu/chimera/docs/UsersGuide/tutorials/bfactor.html

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] Defining an ellipse around a protein structure

2011-05-31 Thread Eric Pettersen

On May 28, 2011, at 4:00 PM, Vandu Murugan wrote:


Dear all,
Is there a program or server that would define a ellipsoid  
around a
given protein molecule?  I would also like to calculate the axis  
components

of the defined ellipse. Thanks in advance..


Chimera can do this, as per:

http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-May/006407.html

--Eric




Re: [ccp4bb] mutation and minimization

2011-05-13 Thread Eric Pettersen

On May 12, 2011, at 4:00 PM, CCP4BB automatic digest system wrote:


Hey all,

I would like to introduce point mutations in a structure and quickly
(and dirtily) minimize the new residue.  (Best rotamer dependent on
local environment, or the like.)  What are simple approaches that  
don't

involve VMD/NAMD or some such overkill.


Chimera is pretty good for this.  It has a Rotamer tool for making the  
substitution based on Dunbrack or Richardson libraries, and can screen  
based on probability / H-bonds formed / steric clashes:


http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/rotamers/rotamers.html

You can then use Chimera's Minimize Structure tool to minimize the  
side chain and/or the local environment or, if you're feeling frisky,  
the whole protein:


http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/minimize/minimize.html

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu

Chimera home page:  www.cgl.ucsf.edu/chimera

Re: [ccp4bb] mutation and minimization

2011-05-13 Thread Eric Pettersen

On May 13, 2011, at 4:00 PM, Tom J. Brett wrote:

On a similar extension to this topic, is there a good software out  
there for doing these kinds of modifications and minimizations for  
protein structures with chemicals entities (i.e., protein/inhibitor  
complexes). What I am looking to do is take a protein/inhbitor  
complex and add chemical modifications onto the inhibitor and see  
how they fit. so a little minimization/relaxation of the surrounding  
protein would be nice. It there a freely available software that  
will allow something like this, if not what commercial packages are  
cheapest/best at this? Thanks and sorry for hijacking this thread.


Chimera also has a tool (Build Structure) that can be used for simple  
chemical modifications:


http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/editing/editing.html

It's not as nice as some builders, but it can generally get the job  
done.  Chimera can compute the charges on your modified structure by  
calling out to AmberTools/antechamber (included with the Chimera  
download) and then perform minimization.  Charge computation might  
take awhile if the inhibitor is large-ish.


This isn't to say that you shouldn't follow the protocol outlined by  
Saugata Hazra in another posting.  It depends on your goals.


--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] Thermal Ellipsoids pic

2011-05-09 Thread Eric Pettersen

On May 9, 2011, at 4:00 PM, Kenneth A. Satyshur wrote:


YO!

I am making thermal ellipsoid plots of a highly flexibile region of a
chromophore using rastep and render from the raster3D package of  
Merritt el al.
But I cannot control the orientation. It obscures the rest of the  
atoms in the
plane of the molecule. Does anyone know a way to orientate molecules  
AND

get a publication quality png file? ccp4mg only produces vague traces.


Hey, you can use Chimera to orient the molecule as you like and then  
use Chimera's Thermal Ellipsoids tool to draw the ellipsoids.  That  
tool lacks rastep's cut-away quadrant depiction but has all the  
other depictions.  More info about Thermal Ellipsoids here:


http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/thermal/thermal.html

about Chimera itself here:

www.cgl.ucsf.edu/chimera

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] Comparing two proteins

2011-04-13 Thread Eric Pettersen

On Apr 13, 2011, at 4:00 PM, Rex Palmer wrote:

What is the best program to use for comparing two protein structures  
which are very similar both structurally and wrt aa sequence? ie to  
get the rms deviations both generally and in selected regions.


Best is kind of subjective, but you can use the MatchMaker tool in  
Chimera to superimpose the structures and show the superposition.   
MatchMaker also has an option to show the corresponding sequence  
alignment.  The alignment will have a histogram of the column RMSD  
values running across the top.  You can drag the mouse on the  
alignment to select/highlight the corresponding parts of the structure  
and to show the RMSD of the dragged region in the alignment's region- 
browser window.


MatchMaker described here: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/matchmaker/matchmaker.html

Chimera obtainable here:  http://www.cgl.ucsf.edu/chimera

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] program to calculate electron density at x,y,z

2011-04-02 Thread Eric Pettersen

On Apr 1, 2011, at 4:00 PM, Ed Pozharski wrote:


I need to calculate the electron density values for a list of spatial
locations (e.g. atom positions in a model) using an mtz-file that
already contains map coefficients.  To write my own code may be easier
than I think (if one can manipulate mtz columns, isn't the only  
problem
left how to incorporate symmetry-related reflections?), but I would  
need

an alternative at least for troubleshooting purposes. So,

Does anyone know of a software tool that can calculate point electron
density for every atom in a structure?

If I would have to bring a dependency into this, the best choice for  
me

would be clipper libs.


Hi Ed,
	You'd have to convert the map to ccp4 or XPLOR format, but once  
you've done that it's pretty straightforward to get the values at atom  
positions using Chimera, pretty much as described here:


http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-March/006202.html

except skipping the part where you average them across a residue.

--Eric

Eric Pettersen
UCSF Computer Graphics Lab


Re: [ccp4bb] Model Building: continuous update of distances as fragment moved

2011-02-15 Thread Eric Pettersen

On Feb 14, 2011, at 4:01 PM, Paul McLaughlin wrote:


In doing some model building I want to move a domain of a protein
manually, as a rigid body, and see in real time continuous updates of
some distances from points on the domain to the rest of the protein  
(we

have cross linking data).
I dimly remember being able to do with this O,at least a decade ago
(probably more),  but it doesn't seem to work in current versions (I
have even dimmer recollections that this might have been a special
feature of a particular graphics card in an SGI)
In any event, does anyone know of anything that might allow us to do
this? { I am not asking about computational ways of exploring  
alternate

packing of the domain to satisfy distances ( I know about these) -
rather we want to get a feel for what is possible by seeing it for
ourselves).


DIstance monitors in Chimera (www.cgl.ucsf.edu/chimera) automatically  
update as structures are moved.  Possibly also of interest, Chimera  
can monitor steric clashes/contacts as structures are moved:


http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/findclash/findclash.html

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] rmsd calculation for all atoms

2010-12-29 Thread Eric Pettersen


From:Michael Swan
Subject: rmsd calculation for all atoms.

Dear all,
I am having a bit of trouble finding a program to do an rmsd  
calculation and

give me the differences between all the atoms in the structures.
I have two structures which are identical in the sense that one is  
the apo
protein and one is the bound structure. I would like to superimpose  
them
using only one domain then output the rmsd values for all the atoms  
or at
least individual residues. The structures are not completely  
identical as
some extra residues were visible in the bound structure so I hope  
there is a

program that can ignore those differences.
If anyone knows of a program that will do this calculation I would  
much
appreciate hearing about it. There seem to be many programs that  
will give

an average rmsd over the whole structure or just the region used for
alignment but I haven't found anything that will give me distances  
between

atoms or residues that were not used in the alignment.
Thanks,
Mike.


If the domain is an entire chain, you could use Chimera's MatchMaker  
tool to do the alignment:


http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/matchmaker/matchmaker.html

It will report that RMSD for the atoms used for matching.  To get the  
RMSD for all atoms in common you would have to use the 'rmsd' command  
and specify the parts in common:


http://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/rmsd.html

So if all of chain A and residues 12-20 of chain B were in common, then:

rmsd #0:.a:12-20.b #1:.a:12-20.b

If there were differences in the side chains (say some were  
incomplete) then you might have to restrict the RMSD measurement to  
the backbone (which you might want to do anyway):


rmsd #0:@c,ca,n:12-2...@c,ca,n #1:@c,ca,n:12-2...@c,ca,n

--Eric


Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] CIF format frustrations

2010-07-29 Thread Eric Pettersen


On Jul 29, 2010, at 4:00 PM, CCP4BB automatic digest system wrote:


Date:Thu, 29 Jul 2010 13:49:12 +0100
From:Simon Kolstoe s.kols...@ucl.ac.uk
Subject: CIF format frustrations

Dear all,

What's the best way to convert a coordinate cif file to a pdb (without
compiling a program, installing new utilities, writing script  
files...)?


You could use Chimera (www.cgl.ucsf.edu/chimera).  File..Open to open  
the CIF file, and then File...Save PDB to write it as PDB (probably  
after seeing if it looks right in the graphics window).


--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] How to align a sequence to a know profile

2010-06-01 Thread Eric Pettersen
   I've a protein sequence of known domain. Based on structure  
alignment,
I've got a alignment of those with known structures. Then how to add  
my

sequence to the alignment?Any suggestions?


Hi Yuan,
	You can use UCSF Chimera's Multalign Viewer tool to do this.  Just  
open your alignment with that tool and use the Edit-Add Sequence  
menu item to add your sequence to the alignment.


UCSF Chimera: http://www.cgl.ucsf.edu/chimera/
Multalign Viewer tool: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/multalignviewer/framemav.html
specifically the Add Sequence part: 
http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/multalignviewer/multalignviewer.html#mavmenu-edit

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu




Re: [ccp4bb] programmatic symmetry mate generation

2010-04-26 Thread Eric Pettersen

Hi James,
	The symmetry copies of a molecule asymmetric unit can be saved to a  
new PDB file using the saveunitcell.py Python script with Chimera.


http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts

It uses the CRYST1 record in the PDB file to determine the symmetry  
group.  Unfortunately as mentioned in the script you need to use  
Chimera as the Python interpreter, which you may not be willing to do  
-- but you can just call out to the command-line version of Chimera  
(also documented in the script) via os.system() or whatnot.
	There is also a sym command in Chimera which can generate symmetry  
mates from a wide spectrum of possible criteria.  A Python script to  
call sym is just:


from chimera import runCommand
runCommand(sym options you want)

--Eric

Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu