Closed blakemertz closed 6 months ago
Hi Blake,
If I understand your question correctly, the simplest way to fix it is to use the coordinates you want and paste the connectivity and atom type from final.xyz
As for the functionalities besides parametrize a small molecule
in poltype, they are not extensively tested and will be removed from poltype.
If the atom order is the same, you just need copy the atom type column.
Thanks, Pengyu
From: Chengwen Liu @.***> Sent: Friday, April 26, 2024 3:19 PM To: TinkerTools/poltype2 Cc: Subscribed Subject: Re: [TinkerTools/poltype2] retaining original coordinates of input molecules (Issue #44)
Hi Blake, If I understand your question correctly, the simplest way to fix it is to use the coordinates you want and paste the connectivity and atom type from final.xyz As for the functionalities besides parametrize a small molecule in poltype, they are not extensively tested and will be removed from poltype.
— Reply to this email directly, view it on GitHubhttps://github.com/TinkerTools/poltype2/issues/44#issuecomment-2080064236, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABNC6XTGXCJESWSVAPQ45KTY7KZERAVCNFSM6AAAAABG3LHEJCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDAOBQGA3DIMRTGY. You are receiving this because you are subscribed to this thread.Message ID: @.***>
@leucinw @pren thank you for the clarification -- copying the coordinates from my original sdf into the final.xyz worked.
I am trying to set up some MD runs to simulate a series of small molecules against a reference protein target using AMOEBA, and was able to generate parameter files for the small molecules in poltype2. However, when I compare my final.xyz and my input.sdf against each other, final.xyz has been translated away from the original position of input.sdf. My input.sdf was docked into the binding site of the target protein, so I was really hoping to avoid having to align my final.xyz to my input.sdf in order to set up my MD simulations.
The description of aligning a ligand to the coordinates of an unoccupied protein in the documentation here, but the impression I get from that documentation is that defining
still assumes that your ligand is docked in the correct spatial coordinates with respect to the protein. If I am missing something in my understanding of the documentation, I apologize.