About Restrain Files for ligands
Hello all, I hope everybody is doing well. Happy New Year 2021!!! I am stuck with some unusual behavior of restrain files generated by Phenix.ReadySet. I crystallized the protein with ligands which have double bonds between the atoms. As I see clear electron density for the ligand, I drew a ligand molecule in Coot (Ligand builder) and fitted in electron density very nicely. Until here, I can see the double bonds in the modelled ligand. I generated the restrain files using Phenix.ReadySet to further proceed for refinement. Unfortunately, after refinement the double bonds were converted to single bonds. I tried to generate new restrain files and refined the model, resulting in similar output. But, fortunately in one case, the double bonds were retained after refinement also. I tried to repeat with other ligands, where I was unsuccessful to retain the double bonds. please help me. Thank you so much. Regards Sharan -- Sharan Karade Postdoc-fellow IBBR-UMD, 9600 Gudelsky Dr, Rockville Maryland 20850
Hi Sharan, I'm guessing your are over-relying on graphic representations of your molecules which is discussed here to some extent: http://phenix-online.org/newsletter/CCN_2016_01.pdf#page=10 What really matters is that the ligand fits the map and its geometry matches what's expected. Pavel On 1/7/21 10:20, Sharan Karade wrote:
Hello all,
I hope everybody is doing well. Happy New Year 2021!!!
I am stuck with some unusual behavior of restrain files generated by Phenix.ReadySet. I crystallized the protein with ligands which have double bonds between the atoms. As I see clear electron density for the ligand, I drew a ligand molecule in Coot (Ligand builder) and fitted in electron density very nicely. Until here, I can see the double bonds in the modelled ligand. I generated the restrain files using Phenix.ReadySet to further proceed for refinement. Unfortunately, after refinement the double bonds were converted to single bonds. I tried to generate new restrain files and refined the model, resulting in similar output. But, fortunately in one case, the double bonds were retained after refinement also. I tried to repeat with other ligands, where I was unsuccessful to retain the double bonds.
please help me. Thank you so much.
Regards Sharan
-- Sharan Karade Postdoc-fellow IBBR-UMD, 9600 Gudelsky Dr, Rockville Maryland 20850
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
Hello,
Naively, Perhaps the issue is that the molecule used to generate the .cif
restraints does not have the hydrogens in the right place/ right
hybridization?
My suggestion would be to check it i.e. with Avogadro or any other chemical
structure editor.
Cheers
Maxime
Le jeu. 7 janv. 2021 à 14:04, Pavel Afonine
Hi Sharan,
I'm guessing your are over-relying on graphic representations of your molecules which is discussed here to some extent:
http://phenix-online.org/newsletter/CCN_2016_01.pdf#page=10
What really matters is that the ligand fits the map and its geometry matches what's expected.
Pavel
On 1/7/21 10:20, Sharan Karade wrote:
Hello all,
I hope everybody is doing well. Happy New Year 2021!!!
I am stuck with some unusual behavior of restrain files generated by Phenix.ReadySet. I crystallized the protein with ligands which have double bonds between the atoms. As I see clear electron density for the ligand, I drew a ligand molecule in Coot (Ligand builder) and fitted in electron density very nicely. Until here, I can see the double bonds in the modelled ligand. I generated the restrain files using Phenix.ReadySet to further proceed for refinement. Unfortunately, after refinement the double bonds were converted to single bonds. I tried to generate new restrain files and refined the model, resulting in similar output. But, fortunately in one case, the double bonds were retained after refinement also. I tried to repeat with other ligands, where I was unsuccessful to retain the double bonds.
please help me. Thank you so much.
Regards Sharan
-- Sharan Karade Postdoc-fellow IBBR-UMD, 9600 Gudelsky Dr, Rockville Maryland 20850
_______________________________________________ phenixbb mailing [email protected]http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
To add to this: best way to generate restraints file within Phenix framework is to use eLBOW with a SMILES string (using GUI or CL), then quickly inspect the restraints file (aka 'ligand CIF file') to make sure it reflects what you know about the ligand. Pavel On 1/8/21 08:02, Maxime Cuypers wrote:
Hello, Naively, Perhaps the issue is that the molecule used to generate the .cif restraints does not have the hydrogens in the right place/ right hybridization? My suggestion would be to check it i.e. with Avogadro or any other chemical structure editor.
Cheers Maxime
Le jeu. 7 janv. 2021 à 14:04, Pavel Afonine
mailto:[email protected]> a écrit : Hi Sharan,
I'm guessing your are over-relying on graphic representations of your molecules which is discussed here to some extent:
http://phenix-online.org/newsletter/CCN_2016_01.pdf#page=10 http://phenix-online.org/newsletter/CCN_2016_01.pdf#page=10
What really matters is that the ligand fits the map and its geometry matches what's expected.
Pavel
On 1/7/21 10:20, Sharan Karade wrote:
Hello all,
I hope everybody is doing well. Happy New Year 2021!!!
I am stuck with some unusual behavior of restrain files generated by Phenix.ReadySet. I crystallized the protein with ligands which have double bonds between the atoms. As I see clear electron density for the ligand, I drew a ligand molecule in Coot (Ligand builder) and fitted in electron density very nicely. Until here, I can see the double bonds in the modelled ligand. I generated the restrain files using Phenix.ReadySet to further proceed for refinement. Unfortunately, after refinement the double bonds were converted to single bonds. I tried to generate new restrain files and refined the model, resulting in similar output. But, fortunately in one case, the double bonds were retained after refinement also. I tried to repeat with other ligands, where I was unsuccessful to retain the double bonds.
please help me. Thank you so much.
Regards Sharan
-- Sharan Karade Postdoc-fellow IBBR-UMD, 9600 Gudelsky Dr, Rockville Maryland 20850
_______________________________________________ phenixbb mailing list [email protected] mailto:[email protected] http://phenix-online.org/mailman/listinfo/phenixbb http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe:[email protected] mailto:[email protected]
_______________________________________________ phenixbb mailing list [email protected] mailto:[email protected] http://phenix-online.org/mailman/listinfo/phenixbb http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected] mailto:[email protected]
Hi Pavel, and Maxime
Thank you so much for your valuable suggestions. I will try to generate the
restrain file as you suggested and re-run the refinement.
Regards
Sharan
On Fri, Jan 8, 2021 at 1:16 PM Pavel Afonine
To add to this: best way to generate restraints file within Phenix framework is to use eLBOW with a SMILES string (using GUI or CL), then quickly inspect the restraints file (aka 'ligand CIF file') to make sure it reflects what you know about the ligand.
Pavel
On 1/8/21 08:02, Maxime Cuypers wrote:
Hello, Naively, Perhaps the issue is that the molecule used to generate the .cif restraints does not have the hydrogens in the right place/ right hybridization? My suggestion would be to check it i.e. with Avogadro or any other chemical structure editor.
Cheers Maxime
Le jeu. 7 janv. 2021 à 14:04, Pavel Afonine
a écrit : Hi Sharan,
I'm guessing your are over-relying on graphic representations of your molecules which is discussed here to some extent:
http://phenix-online.org/newsletter/CCN_2016_01.pdf#page=10
What really matters is that the ligand fits the map and its geometry matches what's expected.
Pavel
On 1/7/21 10:20, Sharan Karade wrote:
Hello all,
I hope everybody is doing well. Happy New Year 2021!!!
I am stuck with some unusual behavior of restrain files generated by Phenix.ReadySet. I crystallized the protein with ligands which have double bonds between the atoms. As I see clear electron density for the ligand, I drew a ligand molecule in Coot (Ligand builder) and fitted in electron density very nicely. Until here, I can see the double bonds in the modelled ligand. I generated the restrain files using Phenix.ReadySet to further proceed for refinement. Unfortunately, after refinement the double bonds were converted to single bonds. I tried to generate new restrain files and refined the model, resulting in similar output. But, fortunately in one case, the double bonds were retained after refinement also. I tried to repeat with other ligands, where I was unsuccessful to retain the double bonds.
please help me. Thank you so much.
Regards Sharan
-- Sharan Karade Postdoc-fellow IBBR-UMD, 9600 Gudelsky Dr, Rockville Maryland 20850
_______________________________________________ phenixbb mailing [email protected]http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
-- Sharan Karade Postdoc-fellow IBBR-UMD, 9600 Gudelsky Dr, Rockville Maryland 20850
On 07/01/2021 18:20, Sharan Karade wrote:
I am stuck with some unusual behavior of restrain files generated by Phenix.ReadySet. I crystallized the protein with ligands which have double bonds between the atoms. As I see clear electron density for the ligand, I drew a ligand molecule in Coot (Ligand builder) and fitted in electron density very nicely. Until here, I can see the double bonds in the modelled ligand.
It is highly likely that you have restraints at this point.
I generated the restrain files using Phenix.ReadySet to further proceed for refinement.
How?
Unfortunately, after refinement the double bonds were converted to single bonds.
refinement won't change the bond orders.
I tried to generate new restrain files and refined the model,
How?
resulting in similar output.
Similar but different?
But, fortunately in one case, the double bonds were retained after refinement also. I tried to repeat with other ligands, where I was unsuccessful to retain the double bonds.
Hmm... Probably the best way of interacting between Coot and ELBOW at the moment is, after you have drawn your ligand in Lidia, click on the red-lipped smiley icon to capture the SMILES string - and then feed that to ELBOW. On 07/01/2021 20:03, Pavel Afonine wrote:
I'm guessing your are over-relying on graphic representations of your molecules which is discussed here to some extent:
http://phenix-online.org/newsletter/CCN_2016_01.pdf#page=10
What really matters is that the ligand fits the map and its geometry matches what's expected.
Re Figure 3b: Oh yikes, that doesn't look right - the elements have been lost. Can you send me that dictionary so that I can try to reproduce the problem? Paul.
participants (4)
-
Maxime Cuypers
-
Paul Emsley
-
Pavel Afonine
-
Sharan Karade