pymol-users — Collaborative support and tips for PyMOL users

[PyMOL] Help - how mutation affects protein structure

[susanna p. <pag...@gm...>]

Hello,

I am quite new to PyMOL. I would like to assess how an amino acid change might affect the protein structure. I used the Wizard Mutagenesis to change the amino acid of interest but don’t know how I can see how the protein structure is affected by the mutation.

Can you please help me with this?

Many thanks in advance.

Best wishes,
Susanna

Re: [PyMOL] Virtual Annual CCP4/APS Crystallographic School (2022) in the US

[Xu, Q. <qx...@an...>]

Dear Colleagues

This is a final reminder that the registration for the upcoming CCP4/APS school will be closed on Apr 15, 2022. Please submit your application soon if you plan to attend.


Charles, Andrey, Garib and Qingping



________________________________
From: Xu, Qingping
Sent: Monday, January 24, 2022 10:07 AM
To: CC...@JI... <CC...@JI...>; pym...@li... <pym...@li...>; phe...@ph... <phe...@ph...>; Ahrens, Kristin <ka...@an...>
Subject: Virtual Annual CCP4/APS Crystallographic School (2022) in the US


Dear Colleagues,

We are pleased to announce the virtual annual CCP4/APS crystallographic school “From data collection to structure refinement and beyond” will be held on June 13-24, 2022 at the Advanced Photon Source (APS), Argonne National Laboratory (ANL), near Chicago, Illinois, USA. All details can be found at the school website: http://www.ccp4.ac.uk/schools/APS-2022/index.php.

Dates: June 6 (Bioinformatics & setup) and June 13 through 24, 2022 (with a weekend break on June 18 and 19th)

Location: Virtual/Remote

The school comprises two parts: remote data collection workshop and virtual crystallographic computing workshop. Data collection workshop includes beamline training, data collection on GM/CA at APS beamlines 23ID-D and 23ID-B equipped with Pilatus3 6M and Eiger 16M detectors respectively, and data processing. For data collection, only the participants’ crystals will be used. Crystallographic computation workshop will feature many modern crystallographic software packages taught by authors and other experts. The daily schedule will be organized in three sections – lectures, tutorials, and hands-on (interactive troubleshooting of the technical difficulties the participants face in their projects). We have had considerable success resolving these problems in past years, attested by resulting publications (see http://www.ccp4.ac.uk/schools/APS-school/publications.php). A draft of the program, contact info and other details can be found at the School website. Note that for the virtual school, an additional day (June 6th, 2022) was scheduled to prepare for the main event, with lectures on generating AlphaFold models to use in later part of the school.

Applicants: Graduate students, postdoctoral researchers and early-career faculty, along with commercial/industrial researchers are encouraged to apply. Only about 20 applicants will be selected for participation. Participants of the workshop are strongly encouraged to bring their own problem data sets or crystals so the problems can be addressed during data collection and/or computation workshops.

Application: Application deadline is April 15th, 2022. To apply, visit https://www.ccp4.ac.uk/schools/APS-2022/application.php.

Fees: There is no fee for attending the virtual school for the chosen applicants.

We hope to see you at the school.

Charles, Andrey, Garib and Qingping

[PyMOL] Problem with the legacy plugin CASTp pocket loader

[Lionel M. P. <lio...@ip...>]

Dear all,

I encounter problems using the legacy plugin CASTp pocket loader on 
PyMOL Version 2.4.1. It fails to open poc files and terminates with 
"Error: unsupported file type poc". Any idea on how to solve this problem?

Thanks a lot!

With kind regards,

Lionel
-- 
logo-ipbs


  Lionel Mourey


    CNRS research Director | Team leader

+33 5 61 17 54 36 | +33 6 73 50 95 51 | Lio...@ip... 
<mailto:Lio...@ip...>

UMR5089 | CNRS - UT3 | 205 Route de Narbonne BP 64182 - 31077 Toulouse 
Cedex 4

ipbs.fr <https://www.ipbs.fr> | cribligand.ipbs.fr 
<http://cribligand.ipbs.fr/>

@IpbsToulouse <https://twitter.com/IpbsToulouse>

ipbs.cnrs <https://www.facebook.com/IPBS.CNRS/>

[PyMOL] Ambient setting for a selection

[Daniel B. <db...@ia...>]

Hello,

Is there a way to get  “set ambient” to work on a selection?  

I am trying to find a way to have a single atom sphere representation appear to glow.

Thank you,

Dan

Re: [PyMOL] hydrogen-bonds detected in the protein-ligand complexes

[EDUARDO J. A. C. <ed...@ep...>]

Good Question
I like to know Too.
And about Hydrophobic interaction , Is it possible to show the
hydrophobic interactions on pymol too?

Em qui., 31 de mar. de 2022 às 06:49, Enrico Martinez <
jms...@gm...> escreveu:

> Dear Pymol users!
> Another question related to the analysis of the multi-model pdb
> obtained from protein-ligand docking calculations.
>
> I am using pymol to visualise protein-ligand interactions between each
> pose predicted by docking and protein in my script:
>
> preset.ligand_cartoon('my_pdb*',_self=cmd)
>
> then I may use the following command to switch to the state in
> multi_model pdb, which could be important:
> cmd.set('state','5','my_pdb*')
>
> Is it possible to print the number of the state which has
> 1) the maximal number of the hydrogen bonds?
> 2) has a particular hydrogen bond of my interest?
> Cheers,
> Enrico
>
>
> _______________________________________________
> PyMOL-users mailing list
> Archives: http://www.mail-archive.com/pym...@li...
> Unsubscribe:
> https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe
>


-- 

*Eduardo José Azevedo Corrêa*

Biólogo| Professor e Pesquisador


(37)9997-5771

(37)3271-4673

*www.epamig.br <http://www.epamig.br/>*

*EPAMIG Pitangui*
Rodovia BR - MG352 Km35 Zona Rural

Caixa Postal 43

CEP 35650-000 - Pitangui - MG

[PyMOL] hydrogen-bonds detected in the protein-ligand complexes

[Enrico M. <jms...@gm...>]

Dear Pymol users!
Another question related to the analysis of the multi-model pdb
obtained from protein-ligand docking calculations.

I am using pymol to visualise protein-ligand interactions between each
pose predicted by docking and protein in my script:

preset.ligand_cartoon('my_pdb*',_self=cmd)

then I may use the following command to switch to the state in
multi_model pdb, which could be important:
cmd.set('state','5','my_pdb*')

Is it possible to print the number of the state which has
1) the maximal number of the hydrogen bonds?
2) has a particular hydrogen bond of my interest?
Cheers,
Enrico

[PyMOL] Yale morph server results. All scenes displayed.

[Bruce O. <bco...@gm...>]

Hi!

I just got output from the Yale morph server
*http://morph2.molmovdb.org/submit.html
<http://morph2.molmovdb.org/submit.html>* that allows construction of
movies given structures of protein conformers.

I opened the result in PyMol and see all the scenes displayed at the same
time.

Does anyone have experience creating a movie from their output?

Thx!

Bruce

Re: [PyMOL] Analysis of docking results in multi-model format

[Enrico M. <jms...@gm...>]

just as an update of my first topic.
Is there any possibility to associate any plugin or external tool
(rdkit) to calculate pharmacophore model for each docking pose loaded
in multi-model format (i) as well as for the ligand in X-ray structure
(ii) and then compare i with ii in order to find instantly the docking
pose which is more similar to X-ray structure according to some part
of the pharmacophore shared between the both molecules ?

чт, 24 мар. 2022 г. в 14:11, Enrico Martinez <jms...@gm...>:
>
> Dear Pymol Users!
> I am dealing with the analysis of the results of protein-ligand
> docking poses representing the multi-model pdb.  I need to find a
> possibility (e.g. via some script that could be executed in the pymol)
> to compare each docking pose with the X-ray structure (loaded as the
> separate model in pymol) in order to find automatically the model (=
> docking solution) which may fit better to it (e.g. via estimating RMSD
> of some part of the ligand in each docking solution compared to the
> X-ray structure).
>
> Assuming that the both pdbs ((docking poses, and X-ray structure))
> have been superimposed (based on the protein atoms) how could I
> automatically switch to the model (in the ensemble) with the identical
> position of the ligand as in the X-ray structure?  I would be grateful
> for any suggestions
> With kind regards,
> Enrico

Re: [PyMOL] set_view based on origin and a vector

[Jared S. <jar...@co...>]

Hi Wei -

Fun problem!  Not a "smart" way to do it, but I think functional, using
pseudoatoms and orient to set the view angle:

```
# Create an origin pseudoatom
ori_coords = [0.4017522, 23.15367207, 6.68992345]
cmd.pseudoatom("ori", pos=ori_coords)
color red, ori
show nb_spheres, ori

# Second pseudoatom as a stand-in for the camera
axis = [0.74, 0.18, -0.65]
cam_coords = [sum(x) for x in zip(ori_coords, axis)]
cmd.pseudoatom("cam", pos=cam_coords)

# Adjust the view so the sphere is centered and crosshairs are aligned on
camera Z axis
orient
center ori
turn y, -90

# Print the current view matrix
get_view
```
yields
```
set_view (\
    -0.673201025,    0.023786053,    0.739076912,\
     0.231151611,    0.956164062,    0.179775774,\
    -0.702402592,    0.291863978,   -0.649188757,\
     0.000000000,    0.000000000,  -14.178204536,\
     0.401752204,   23.153671265,    6.689923286,\
    11.178205490,   17.178203583,  -20.000000000 )
```
[image: image.png]

Hope that helps.

Cheers,
Jared


On Fri, Mar 25, 2022 at 12:44 PM Wei Wang <ww...@co...> wrote:

> Dear PyMOL community,
>
> I have a need to adjust my view point in python script. I want to achieve
> the following:
>
>    1. I have a origin with coordinates [ 0.4017522 23.15367207 6.68992345]
>    2. I have a vector/rotation axis [0.74, 0.18, -0.65]
>    3. I want to look directly down that vector using above origin.
>
> Theoretically I should be able to use those information to determine the
> values needed to put in set_view. However, I’m not smart enough to do that.
> Would anyone please show me the way?
>
> Many thanks,
>
> Wei
> _______________________________________________
> PyMOL-users mailing list
> Archives: http://www.mail-archive.com/pym...@li...
> Unsubscribe:
> https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe

[PyMOL] set_view based on origin and a vector

[Wei W. <ww...@co...>]

Dear PyMOL community,

I have a need to adjust my view point in python script. I want to achieve the following:

1. I have a origin with coordinates [ 0.4017522 23.15367207 6.68992345]
2. I have a vector/rotation axis [0.74, 0.18, -0.65]
3. I want to look directly down that vector using above origin.

Theoretically I should be able to use those information to determine the values needed to put in set_view. However, I’m not smart enough to do that. Would anyone please show me the way?

Many thanks,

Wei

Re: [PyMOL] Analysis of docking results in multi-model format

[Enrico M. <jms...@gm...>]

Hello Joel,
Actually I am using VINA that indeed provides us with these RMSD
values. However the goal is to compare the docking poses with the
X-ray structure containing another (but very similar in terms of the
functional groups) ligand.
I've already tried some pymol built-in commands like pair_fit but I
could not adapt it to the multi-model pdb in order to select
"on-the-fly" the pose which has better similarity to the pose observed
in the X-ray structure.
Cheers,
Enrico

пт, 25 мар. 2022 г. в 01:28, Joel Tyndall <joe...@ot...>:
>
> Hi Enrico,
>
> Which docking protocol are you using as this may already have RMSD to a reference ligand and presentation of best poses (based on scoring)
>
> Joel
>
> -----Original Message-----
> From: Enrico Martinez <jms...@gm...>
> Sent: Friday, 25 March 2022 2:11 AM
> To: pymol-users <pym...@li...>
> Subject: [PyMOL] Analysis of docking results in multi-model format
>
> Dear Pymol Users!
> I am dealing with the analysis of the results of protein-ligand docking poses representing the multi-model pdb.  I need to find a possibility (e.g. via some script that could be executed in the pymol) to compare each docking pose with the X-ray structure (loaded as the separate model in pymol) in order to find automatically the model (= docking solution) which may fit better to it (e.g. via estimating RMSD of some part of the ligand in each docking solution compared to the X-ray structure).
>
> Assuming that the both pdbs ((docking poses, and X-ray structure)) have been superimposed (based on the protein atoms) how could I automatically switch to the model (in the ensemble) with the identical position of the ligand as in the X-ray structure?  I would be grateful for any suggestions With kind regards, Enrico
>
>
> _______________________________________________
> PyMOL-users mailing list
> Archives: https://apc01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.mail-archive.com%2Fpymol-users%40lists.sourceforge.net&amp;data=04%7C01%7Cjoel.tyndall%40otago.ac.nz%7C5801f6d4300c4429df4008da0d996aa8%7C0225efc578fe4928b1579ef24809e9ba%7C0%7C0%7C637837249971556989%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&amp;sdata=%2FH4rWyrcccM8UDMCjpXhzfhvbWIRK7iVxT3WTHhd11Q%3D&amp;reserved=0
> Unsubscribe: https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsourceforge.net%2Fprojects%2Fpymol%2Flists%2Fpymol-users%2Funsubscribe&amp;data=04%7C01%7Cjoel.tyndall%40otago.ac.nz%7C5801f6d4300c4429df4008da0d996aa8%7C0225efc578fe4928b1579ef24809e9ba%7C0%7C0%7C637837249971556989%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&amp;sdata=W1bYxrtg1zExUmoFtA3srAUsOc8PAODG2QimqwEY%2FQw%3D&amp;reserved=0

Re: [PyMOL] Analysis of docking results in multi-model format

[Joel T. <joe...@ot...>]

Hi Enrico,

Which docking protocol are you using as this may already have RMSD to a reference ligand and presentation of best poses (based on scoring)

Joel

-----Original Message-----
From: Enrico Martinez <jms...@gm...> 
Sent: Friday, 25 March 2022 2:11 AM
To: pymol-users <pym...@li...>
Subject: [PyMOL] Analysis of docking results in multi-model format

Dear Pymol Users!
I am dealing with the analysis of the results of protein-ligand docking poses representing the multi-model pdb.  I need to find a possibility (e.g. via some script that could be executed in the pymol) to compare each docking pose with the X-ray structure (loaded as the separate model in pymol) in order to find automatically the model (= docking solution) which may fit better to it (e.g. via estimating RMSD of some part of the ligand in each docking solution compared to the X-ray structure).

Assuming that the both pdbs ((docking poses, and X-ray structure)) have been superimposed (based on the protein atoms) how could I automatically switch to the model (in the ensemble) with the identical position of the ligand as in the X-ray structure?  I would be grateful for any suggestions With kind regards, Enrico


_______________________________________________
PyMOL-users mailing list
Archives: https://apc01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.mail-archive.com%2Fpymol-users%40lists.sourceforge.net&amp;data=04%7C01%7Cjoel.tyndall%40otago.ac.nz%7C5801f6d4300c4429df4008da0d996aa8%7C0225efc578fe4928b1579ef24809e9ba%7C0%7C0%7C637837249971556989%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&amp;sdata=%2FH4rWyrcccM8UDMCjpXhzfhvbWIRK7iVxT3WTHhd11Q%3D&amp;reserved=0
Unsubscribe: https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsourceforge.net%2Fprojects%2Fpymol%2Flists%2Fpymol-users%2Funsubscribe&amp;data=04%7C01%7Cjoel.tyndall%40otago.ac.nz%7C5801f6d4300c4429df4008da0d996aa8%7C0225efc578fe4928b1579ef24809e9ba%7C0%7C0%7C637837249971556989%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&amp;sdata=W1bYxrtg1zExUmoFtA3srAUsOc8PAODG2QimqwEY%2FQw%3D&amp;reserved=0

[PyMOL] Analysis of docking results in multi-model format

[Enrico M. <jms...@gm...>]

Dear Pymol Users!
I am dealing with the analysis of the results of protein-ligand
docking poses representing the multi-model pdb.  I need to find a
possibility (e.g. via some script that could be executed in the pymol)
to compare each docking pose with the X-ray structure (loaded as the
separate model in pymol) in order to find automatically the model (=
docking solution) which may fit better to it (e.g. via estimating RMSD
of some part of the ligand in each docking solution compared to the
X-ray structure).

Assuming that the both pdbs ((docking poses, and X-ray structure))
have been superimposed (based on the protein atoms) how could I
automatically switch to the model (in the ensemble) with the identical
position of the ligand as in the X-ray structure?  I would be grateful
for any suggestions
With kind regards,
Enrico

[PyMOL] fixing atom order in multi-model pdb

[Enrico M. <jms...@gm...>]

Dear Pymol Users!
I am dealing with the multi-model pdb containing the docking poses
that were merged with the protein atoms. Initially the part containing
the ligand poses also included some flexible side-chains of the
receptor (obtained from VINA) that were subsequently merged ( in each
model) with the static part of the protein using some script.

The complex may be correctly visualised in pymol, however (due to the
reorganizations of the atoms in the protein part) there are some
problems in the visualisations of the interatomic interactions between
the static part of the protein and the ligands.

May we use some pymol in no-gui mode to "fix" the atomic order problem
in the loaded pdb (in each models!) and then save it as a new pdb?
Here is the model of the workflow:

pymol -c -d "
from pymol import cmddb')
cmd.load('complex.pdb')
# execute some command to fix atom order in the multi-model pdb
cmd.save('complex_fixed_pdb.pdb')
"

Re: [PyMOL] pymol cartoon for hexamer

[<tna...@os...>]

Hi,

You have to predict how six monomers assemble into a hexamer.
This can be done by AlphaFold multimer or ColabFold notebook.

Best regards,

Takanori Nakane

________________________________________
差出人: vivek ambastha <vi...@mi...>
送信日時: 2022年3月22日 12:56
宛先: pym...@li...
件名: [PyMOL] pymol cartoon for hexamer

Hello all,
I am naïve to pymol or even protein structure. I need help in creating the hexameric model of a protein whose structure is in predicted through alpha folds.
I can download the monomer. My objective is to develop a cartoon for the hexamer using these six individual monomers. Is this possible using Pymol? If yes then how to do that

Best,

Vivek.A

[PyMOL] pymol cartoon for hexamer

[vivek a. <vi...@mi...>]

Hello all,
I am naïve to pymol or even protein structure. I need help in creating the hexameric model of a protein whose structure is in predicted through alpha folds.
I can download the monomer. My objective is to develop a cartoon for the hexamer using these six individual monomers. Is this possible using Pymol? If yes then how to do that

Best,

Vivek.A

[PyMOL] RMSD from DCD

[Mcguire, K. <klm...@UC...>]

I have been using VMD for most of my MD analysis, but I am trying to learn PyMol now.  I know how to load a pdb and a dcd onto that pdb.  But how do I measure the RMSD of each frame in the DCD aligned to the starting structure, and then output a graph like of RMSD vs frame? Doing this through scripting would be great. Thanks!

Re: [PyMOL] Combine two states into one state?

[Thomas H. <th...@th...>]

Hi Nicholas,

You can use the "split_states" command.
https://pymolwiki.org/index.php/Split_states

Cheers,
  Thomas

On Wed, Mar 16, 2022 at 6:57 AM Nicholas Gao <dom...@gm...> wrote:
>
> Dear PyMOL users,
>
> I have an antibody-antigen structure where the authors have put the antigen atoms in state 1 and put the antibody atoms in state 2.  Is there a good way via command line to collapse state 1 and 2 into a single state, where both antibody and antigen atoms are viewable together?
>
> Sincerely,
> Nicholas Gao
> _______________________________________________
> PyMOL-users mailing list
> Archives: http://www.mail-archive.com/pym...@li...
> Unsubscribe: https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe

[PyMOL] Combine two states into one state?

[Nicholas G. <dom...@gm...>]

Dear PyMOL users,

I have an antibody-antigen structure where the authors have put the antigen
atoms in state 1 and put the antibody atoms in state 2.  Is there a good
way via command line to collapse state 1 and 2 into a single state, where
both antibody and antigen atoms are viewable together?

Sincerely,
Nicholas Gao

Re: [PyMOL] Extract H-bond angle

[Neena S. E. <nee...@gm...>]

Hello PyMOL users,

With "list_hbonds analyte" command, PyMOL gives the list of all H-bonds
with distances in a protein.
I was wondering how to also extract H-bond angle information for those
h-bonds. Is there a command line for this?

Thank you so much,
Neena

On Thu, 3 Mar 2022 at 19:10, Neena Susan Eappen <nee...@gm...>
wrote:

> Hello PyMOL users,
>
> With "list_hbonds analyte" command, PyMOL gives the list of all H-bonds
> with distances in a protein.
> I was wondering how to also extract H-bond angle information for those
> h-bonds. Is there a command line for this?
>
> Thank you so much,
> Neena
>

Re: [PyMOL] Virtual Annual CCP4/APS Crystallographic School (2022) in the US

[Qingping Xu <qx...@an...>]

Dear Colleagues,

Just a friendly reminder that the registration of the upcoming CCP4/APS 
virtual school is open. If you plan to attend the school, please submit 
your application before the Apr 15th, 2022 deadline, preferably sooner.

The school schedule is available on the CCP4 school website 
(https://www.ccp4.ac.uk/schools/APS-2022/index.php).

Charles, Andrey, Garib and Qingping


On 1/24/22 10:07, Qingping Xu wrote:
>
> Dear Colleagues,
>
> We are pleased to announce the virtual annual CCP4/APS 
> crystallographic school “From data collection to structure refinement 
> and beyond” will be held on June 13-24, 2022 at the Advanced Photon 
> Source (APS), Argonne National Laboratory (ANL), near Chicago, 
> Illinois, USA. All details can be found at the school website: 
> http://www.ccp4.ac.uk/schools/APS-2022/index.php.
>
> Dates: June 6 (Bioinformatics & setup) and June 13 through 24, 2022 
> (with a weekend break on June 18 and 19th)
>
> Location: Virtual/Remote
>
> The school comprises two parts: remote data collection workshop and 
> virtual crystallographic computing workshop. Data collection workshop 
> includes beamline training, data collection on GM/CA at APS beamlines 
> 23ID-D and 23ID-B equipped with Pilatus3 6M and Eiger 16M detectors 
> respectively, and data processing. For data collection, only the 
> participants’ crystals will be used. Crystallographic computation 
> workshop will feature many modern crystallographic software packages 
> taught by authors and other experts. The daily schedule will be 
> organized in three sections – lectures, tutorials, and hands-on 
> (interactive troubleshooting of the technical difficulties the 
> participants face in their projects). We have had considerable success 
> resolving these problems in past years, attested by resulting 
> publications (see 
> http://www.ccp4.ac.uk/schools/APS-school/publications.php). A draft of 
> the program, contact info and other details can be found at the School 
> website. Note that for the virtual school, an additional day (June 
> 6th, 2022) was scheduled to prepare for the main event, with lectures 
> on generating AlphaFold models to use in later part of the school.
>
> Applicants: Graduate students, postdoctoral researchers and 
> early-career faculty, along with commercial/industrial researchers are 
> encouraged to apply. Only about 20 applicants will be selected for 
> participation. Participants of the workshop are strongly encouraged to 
> bring their own problem data sets or crystals so the problems can be 
> addressed during data collection and/or computation workshops.
>
> Application: Application deadline is April 15th, 2022. To apply, visit 
> https://www.ccp4.ac.uk/schools/APS-2022/application.php.
>
> Fees: There is no fee for attending the virtual school for the chosen 
> applicants.
>
> We hope to see you at the school.
>
> Charles, Andrey, Garib and Qingping
>

Re: [PyMOL] Pymol-Alignment Module

[Jared S. <jar...@co...>]

Hi Ecem -

If I understand correctly, you first want to align on one set of residues
and calculate their RMSD, then to calculate the RMSD for the remaining
residues without moving the structures.  This is straightforward using
PyMOL's internal command line.

```
# Create two selections called "sel1" and "sel2"
select sel1, resi 1-20
select sel2, resi 30-100

# Note: if the numbering is different between the structures, you may need
something like this instead:
# select sel1, (struc1 and resi 1-20) or (struc2 and resi 2-21)

# Superimpose the structures (this prints the RMSD to the console)
super struc1 and sel1, struc2 and sel1

# Calculate the "in-place" RMSD for the remaining residues (also printed to
the console)
rms_cur struc1 and sel2, struc2 and sel2
```

Note that if there are any residues with alternate conformations, or if the
numbering is different between the two structures, the `rms_cur` step may
need to be modified as indicated in the examples on its wiki page
<https://pymolwiki.org/index.php/Rms_cur#Examples>.

Hope that helps.

Cheers,
Jared


On Fri, Mar 4, 2022 at 5:37 AM Ecem Güngör <ece...@gm...> wrote:

> Dear Pymol community,
>
> I have a question about alignment module of Pymol. I have two different
> pdb ids and each have about 100 amino acid long. I would like to align the
> specific parts of these proteins and see both how the rest of the proteins
> and also these specific parts superimpose with each other. In other words,
> if I align the proteins based on their 20-30 amino acids, I would like to
> first see what the RMSD is for these specific regions (20-30 aa) and then
> the RMSD of the rest (1-20 and 30-100) is. I truly appreciate your
> suggestions and guidance on this. Thank you!
>
> Have a great week/weekend,
> Ecem
> _______________________________________________
> PyMOL-users mailing list
> Archives: http://www.mail-archive.com/pym...@li...
> Unsubscribe:
> https://sourceforge.net/projects/pymol/lists/pymol-users/unsubscribe

[PyMOL] saving info regarding protein-ligand interactions in the log

[Enrico M. <jms...@gm...>]

Dear Pymol users!
I am using a multi-model pdb containing the 80 predicted docking
solutions as well as the receptor (the both presented in each frame).
I need to calculate all hydrogen bonds (and possible other
interactions) observed in each frame (=docking solution)? I've already
used pymol for the visualizations of the interactions, and find it
very useful for the visual inspection. May I use some pymol script in
order that I could save the information regarding the observed
interactions in the log file?
Many thanks in advance
Cheers,
Enrico

[PyMOL] Pymol-Alignment Module

[Ecem G. <ece...@gm...>]

Dear Pymol community,

I have a question about alignment module of Pymol. I have two different pdb
ids and each have about 100 amino acid long. I would like to align the
specific parts of these proteins and see both how the rest of the proteins
and also these specific parts superimpose with each other. In other words,
if I align the proteins based on their 20-30 amino acids, I would like to
first see what the RMSD is for these specific regions (20-30 aa) and then
the RMSD of the rest (1-20 and 30-100) is. I truly appreciate your
suggestions and guidance on this. Thank you!

Have a great week/weekend,
Ecem

[PyMOL] Extract H-bond angle

[Neena S. E. <nee...@gm...>]

Hello PyMOL users,

With "list_hbonds analyte" command, PyMOL gives the list of all H-bonds
with distances in a protein.
I was wondering how to also extract H-bond angle information for those
h-bonds. Is there a command line for this?

Thank you so much,
Neena