- Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] carbon-chlorine bond
by Elaine Meng 17 Nov '08

17 Nov '08

Hi Francesco, I was going to write back that your PDB file from AMBER has the atom names mis-aligned so that Chimera probably thinks the chloride is a carbon, which would have a shorter automatic bond cutoff. However, when I put those 4 lines from your message below in a file and open it in Chimera, the elements are correctly identified and there is already a C-Cl bond; only the H is unconnected. Perhaps the characters on the far right were sufficient to identify element Cl. The H is far away and must be bonded to something else in the complete structure. Elaine On Nov 17, 2008, at 10:19 AM, Francesco Pietra wrote: > > Hi Elaine: > Yes, it is a pdb file written with program ptraj (Amber) for the > minimum energy. The chlorine atom is defined as (partial block): > > ATOM 14544 H1 AA1 341 76.488 57.341 > 49.169 hn > ATOM 14545 C3 AA1 341 75.036 56.270 > 51.132 ca > ATOM 14546 Cl AA1 341 73.973 56.486 > 49.780 cl > ATOM 14547 N1 AA1 341 74.619 55.535 > 52.192 nb > > and there are no CONECT. As you said above, prmtop/mdcrd with movie do > show a line connecting C-Cl, but I don't know if there is another way > to extract the minimum-energy structure than that pdb file. > > Why are all other bonds shown correctly by Chimera in the above > situation? Why just C-Cl not? > > francesco

1 0

Re: [Chimera-users] carbon-chlorine bond
by Elaine Meng 17 Nov '08

17 Nov '08

Hi Francesco, It only means that whatever file(s) you are reading into Chimera do not specify the bond explicitly. I can't tell from the information given whether your MD program knows the bond is there; it may know, but not write this information into the file that you input to Chimera. If the file(s) describe trajectory in the format of the MD program (AMBER? GROMOS?) and you are using the MD Movie tool in Chimera for display, all the bonds should be there. <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/movie/framemovie.h… > However, if you are doing something else like writing a PDB file from those programs, maybe the PDB file does not specify the bond explicitly, and the other viewers might have a more relaxed criterion for automatically identifying a bond. I don't think there is a way to change the automatic criteria, but you could just add a CONECT record specifying the bond to the PDB file, or just add the bond in Chimera with the command "bond" or the "Add/ Delete Bonds" section of Build Structure. <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/bond.html> <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.ht… > Best, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html On Nov 16, 2008, at 11:00 AM, Francesco Pietra wrote: > For a monochlorinated small alkaloid docked onto a protein model, > after molecular dynamics simulations Chimera shows the chlorine atom > as a dot, without a line to represent the bond. Before dynamics, the > C-Cl line connecting the two atoms was shown. Other viewers show the > bond both before and after MD. For the purpose of presentation, how > to set a longer C-Cl bond length as a temporary default for Chimera? > > Very likely the parameterization for chlorine in the MD suite is not > as it should. > > Thanks > francesco pietra

1 0

clustering of docking solutions
by Sebastian Kruggel 17 Nov '08

17 Nov '08

Hi all, I am looking for a possibility to cluster docking results (multimol2 file) in chimera. In DockView I can't find such tool and in the EnsembleCluster method I don't see a way to set rmsd limits. I found a hint to the MD Analysis tool in the mailing list, ok, I could convert my docking solutions into pdbs and use this, but probably there is another simple and convenient way that I am just too blind to see ;-) Maybe somebody can help? Thanks in advance, Sebastian -- Sebastian Kruggel Institut für Pharmazie Bundesstr. 45 | Raum 112 D 20146 Hamburg Tel +49(0) 40 42838 3626

1 0

Re: [Chimera-users] Running a job in chimera
by Elaine Meng 12 Nov '08

12 Nov '08

Dear Tuhin, The script to perform the actions (your #1 or #2) could be written in python or in Chimera commands. Just opening the python file (*.py) or Chimera command file (*.com) in Chimera will execute it. The trickier part would be how to automatically loop through your multiple structures. You could write a long Chimera command script that includes the names of the structures, or you could use a python script that loops through a list of filenames. The methods can be combined; for example, the python script could merely loop through the names and open the structures, then open a separate Chimera command script. You could also use a shell script to loop through filenames, but I guess that would be slower because it would involve starting a new Chimera for each input (#1) or set of inputs (#2). If that is easier for you, however, it is reasonable. Please see this previous posting for examples/discussion of Chimera scripts for processing multiple structures, calling a script from another script, and using aliases: <http://www.cgl.ucsf.edu/pipermail/chimera-users/2008-October/003218.html > See also Chimera nogui mode: <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/options.html#nogui> I believe all the actions you want to perform can be done with Chimera commands, for example open /location/location/frame20.pdb delete solvent & protein z>5 write format pdb relative 0 0 /location/location/frame20zone5.pdb close 0 or open /location/mystruct1.pdb open /location/mystruct2.pdb open /location/mystruct3.pdb matchmaker #0 #1 matchmaker #0 #2 These are just possibilities out of several ways to perform the tasks. For example, the top example assumes your structure is a protein, and in the bottom example, you could use "match" instead of "matchmaker" to superimpose structures. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html On Nov 11, 2008, at 1:41 AM, tuhin(a)iitk.ac.in wrote: > Dear All, > I have run several simulations. But, I need to analyse a selective > number > of structures from each simulation saved at a given time step "t". > Thus, > for a 20ns simulation, starting t=0 and incrementing it by 50ps, > I'll have > ~400 individual structures. Each structure is within a waterbox, which > contains approx. 20,000 water molecules. Once I get the structures, > I need > to do the following using CHIMERA: > > 1. Open each structure and save water molecules within a cutoff > distance of > 5.0 Angstrom from the protein. Thus, I'll end-up having a shell > (<=5.0A) > of water molecules around the protein. The saved structure will be > in > .pdb format and include the protein with 5.0A water shell. > I know how to do it in CHIMERA using the commandline option/ > pulldown > menu, but, how to do it for 400 structures. Is there a way to run it > in a "Batch Job" ? > > 2. I also need to superimpose Ca-Ca atoms of 400 structures. But, I > feel > there is a limit to the no. of structures that CHIMERA could handle, > thus, at max. I would be superimposing 20-30 structures at a time. > How > could I automatize the same through a script sothat I could run it > in > batches? > > Thanks in advance. Any suggestion(s) is welcomed. > Warm regards, > Tuhin

2 1

zone command
by Francesco Pietra 12 Nov '08

12 Nov '08

A polypeptide was docked onto a protein that already carried ligands. Then, molecular dynamics of the complex was carried out. To do this, a single PDB file was created with the program ptraj of Amber. Although ptraj removes chain labels, TER lines, etc., renumbering all the residues in a single sequence, I know, from the ligand positions, the corresponding numbering for the chains and for the polypeptide. Which is the best procedure to map around certain groups of amino acids of the protein, i.e., those that constitute the experimentally-evaluated binding site? I want to take notice of which residues of the polypeptide are most involved in the binding. Thanks francesco pietra

2 1

Re: [Chimera-users] Powerpoint
by CHIU,YEUNG 12 Nov '08

12 Nov '08

Hi John, According to my experience, you can not directly insert the video file from Chimera into your powerpoint slides, but only a superlink to the video. Usually what I do is I set up a buttom with the superlink, and everytime I click the buttom, the video can be played by Quicktime~ If you get some other way to do it, please let me know~ thank you~ On Wed Nov 12 07:28:07 EST 2008, "Beale, John" <jbeale(a)stlcop.edu> wrote: > Can movies created in Chimera be played in a Powerpoint > presentation? If > so, how do I set it up? > > John Beale > > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > Sincerely yours CHIU,YEUNG

1 0

Powerpoint
by Beale, John 12 Nov '08

12 Nov '08

Can movies created in Chimera be played in a Powerpoint presentation? If so, how do I set it up? John Beale

1 0

Running a job in chimera
by tuhin＠iitk.ac.in 11 Nov '08

11 Nov '08

Dear All, I have run several simulations. But, I need to analyse a selective number of structures from each simulation saved at a given time step "t". Thus, for a 20ns simulation, starting t=0 and incrementing it by 50ps, I'll have ~400 individual structures. Each structure is within a waterbox, which contains approx. 20,000 water molecules. Once I get the structures, I need to do the following using CHIMERA: 1. Open each structure and save water molecules within a cutoff distance of 5.0 Angstrom from the protein. Thus, I'll end-up having a shell (<=5.0A) of water molecules around the protein. The saved structure will be in .pdb format and include the protein with 5.0A water shell. I know how to do it in CHIMERA using the commandline option/pulldown menu, but, how to do it for 400 structures. Is there a way to run it in a "Batch Job" ? 2. I also need to superimpose Ca-Ca atoms of 400 structures. But, I feel there is a limit to the no. of structures that CHIMERA could handle, thus, at max. I would be superimposing 20-30 structures at a time. How could I automatize the same through a script sothat I could run it in batches? Thanks in advance. Any suggestion(s) is welcomed. Warm regards, Tuhin

1 0

How to make bonds between two pdb files
by Ann Spevacek 06 Nov '08

06 Nov '08

Hello - I am trying to make a bond between a protein and a polysaccharide chain in Chimera and am not quite sure how to do it. I am also curious as to whether the program will take into account interactions such as sterics when it generates the bond. Any guidance will be greatly appreciated! Thanks, Ann Spevacek Millhauser Lab UCSC

2 1

successive minimizations: can we use a unique protein setup?
by Jean Didier Pie Marechal 06 Nov '08

06 Nov '08

Hi everyone, I want to do the following: 1. minimize a region of a structure 2. change some part of the region by hand (like swap a residue) 3. minimize once more. and this 2 or three times. However, each time I want to minimize I have to pass through the all setup of the protein (add H, hist protonation, and so one). Couldn't we keep the information of the first setup and directly fo through the minimization? All the best, JD Dr. Jean-Didier Maréchal Lecturer Computational Bioorganic and Bioionorganic Chemistry @ Transmet Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dimecres, Novembre 5, 2008 2:00 am Assumpte: Chimera-users Digest, Vol 67, Issue 4 > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > or, via email, send a message with subject or body 'help' to > chimera-users-request(a)cgl.ucsf.edu > > You can reach the person managing the list at > chimera-users-owner(a)cgl.ucsf.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Chimera-users digest..." > > > Today's Topics: > > 1. Using match (Omar Davulcu) > 2. Re: Using match (Eric Pettersen) > 3. Re: Models or chain (Eric Pettersen) > 4. removing side chains (Matthew Dougherty) > 5. Re: removing side chains (Eric Pettersen) > 6. Re: Using match (Elaine Meng) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Mon, 3 Nov 2008 14:53:19 -0800 > From: "Omar Davulcu" <davulcuo(a)ohsu.edu> > Subject: [Chimera-users] Using match > To: chimera-users(a)cgl.ucsf.edu > Message-ID: <001301c93e06$f7394020$e5abc060$@edu> > Content-Type: text/plain; charset="us-ascii" > > Hi, > > > > I'm using the match command with a cutoff to align specific regions of > related enzymes. It appears to work fine and remove residue pairs > above the > cutoff, but is there a way to list those residue pairs which were > (or were > not) removed? > > > > Thanks for any help! > > Omar Davulcu > >

2 1