Re: [Rdkit-discuss] Alignment using LIgpargen file
Dear Phuong,
you might just update the coordinates of the LigParGen PDB with the
aligned coordinates.
Cheers,
p.
On 08/23/18 18:23, Phuong Chau wrote:
Dear Paolo,
The new chemical chemical after alignment (UNK.pdb) file is in format
of RDKIT pdb file. Is it possible for the output in format of
LIgParGen? I am trying to run stimulation with the aligned chemical
but it seems like it can not read pdb file generated from RDKit.
Best,
Phuong Chau
On Wed, Aug 22, 2018 at 12:32 PM, Paolo Tosco
mailto:paolo.tosco.m...@gmail.com>> wrote:
Dear Phuong,
it works for me with the latest RDKit release. I am wondering if
you are using an older RDKit version where maybe the PDB parser
was not assigning formal charges based on connectivity? Then you
would have to manually set the formal charge of the protonated
nitrogen atom (# 2 N) to +1.
Cheers,
p.
On 08/22/18 19:06, Phuong Chau wrote:
Hi Paolo,
I got this new error when I tried with different chemicals:
python 3DAlignmentwith3OAwithGRO.py CCNC\(N\)N_afterMD.pdb
CC\[NH2\+\]CCN_Ligpargen.pdb
[12:02:45] Explicit valence for atom # 2 N, 4, is greater than
permitted
Traceback (most recent call last):
File "3DAlignmentwith3OAwithGRO.py", line 12, in
rdDistGeom.EmbedMolecule(prbMolwithH)
Boost.Python.ArgumentError: Python argument types in
rdkit.Chem.rdDistGeom.EmbedMolecule(NoneType)
did not match C++ signature:
EmbedMolecule(RDKit::ROMol {lvalue} mol, unsigned int
maxAttempts=0, int randomSeed=-1, bool clearConfs=True, bool
useRandomCoords=False, double boxSizeMult=2.0, bool
randNegEig=True, unsigned int numZeroFail=1, boost::python::dict
{lvalue} coordMap={}, double forceTol=0.001, bool
ignoreSmoothingFailures=False)
The second argument I also downloaded it from Ligpargen but it
gave me this error. I am not so sure why it happened.
I also attached the two files :
On Tue, Aug 21, 2018 at 3:19 PM, Paolo Tosco
mailto:paolo.tosco.m...@gmail.com>>
wrote:
Hi Phuong,
it does have hydrogens after the alignment:
$ ls
UNK_A1C198.pdb align.py drug.pdb
$ python align.py drug.pdb UNK_A1C198.pdb
$ ls
UNK.pdb UNK_A1C198.pdb align.py drug.pdb
$ cat UNK.pdb
ATOM 1 C00 UNK 1 37.884 56.016 51.678 1.00
0.00 C
ATOM 2 C01 UNK 1 37.501 54.554 51.968 1.00
0.00 C
ATOM 3 C02 UNK 1 36.085 54.316 52.531 1.00
0.00 C
ATOM 4 C03 UNK 1 35.697 55.221 53.715 1.00
0.00 C
ATOM 5 C04 UNK 1 34.907 56.485 53.304 1.00
0.00 C
ATOM 6 C05 UNK 1 35.788 57.733 53.108 1.00
0.00 C
ATOM 7 C06 UNK 1 36.227 57.998 51.654 1.00
0.00 C
ATOM 8 C07 UNK 1 36.813 56.799 50.889 1.00
0.00 C
ATOM 9 H08 UNK 1 38.821 56.001 51.080 1.00
0.00 H
ATOM 10 H09 UNK 1 38.160 56.522 52.625 1.00
0.00 H
ATOM 11 H0A UNK 1 37.598 53.969 51.027 1.00
0.00 H
ATOM 12 H0B UNK 1 38.243 54.138 52.684 1.00
0.00 H
ATOM 13 H0C UNK 1 35.325 54.368 51.723 1.00
0.00 H
ATOM 14 H0D UNK 1 36.057 53.262 52.884 1.00
0.00 H
ATOM 15 H0E UNK 1 35.031 54.622 54.375 1.00
0.00 H
ATOM 16 H0F UNK 1 36.586 55.471 54.333 1.00
0.00 H
ATOM 17 H0G UNK 1 34.203 56.712 54.134 1.00
0.00 H
ATOM 18 H0H UNK 1 34.270 56.304 52.411 1.00
0.00 H
ATOM 19 H0I UNK 1 35.198 58.620 53.428 1.00
0.00 H
ATOM 20 H0J UNK 1 36.665 57.714 53.788 1.00
0.00 H
ATOM 21 H0K UNK 1 36.982 58.814 51.667 1.00
0.00 H
ATOM 22 H0M UNK 1 35.352 58.378 51.082 1.00
0.00 H
ATOM 23 H0N UNK 1 35.991 56.144 50.537 1.00
0.00 H
ATOM 24 H0O UNK 1 37.280 57.194 49.960 1.00
0.00 H
CONECT 1 2 8 9 10
CONECT 2 3 11 12
CONECT 3 4 13 14
CONECT 4 5 15 16
CONECT 5 6 17 18
CONECT 6 7 19 20
CONECT 7 8 21 22
CONECT 8 23 24
END
p.
On 08/21/18 20:46, Phuong Chau wrote:
Hi Paolo,
I tried the code, it worked with the alignment but the new
chemical after alignment still does not have any H atoms.
Re: [Rdkit-discuss] Alignment using LIgpargen file
Dear Paolo,
The new chemical chemical after alignment (UNK.pdb) file is in format of
RDKIT pdb file. Is it possible for the output in format of LIgParGen? I am
trying to run stimulation with the aligned chemical but it seems like it
can not read pdb file generated from RDKit.
Best,
Phuong Chau
On Wed, Aug 22, 2018 at 12:32 PM, Paolo Tosco
wrote:
> Dear Phuong,
>
> it works for me with the latest RDKit release. I am wondering if you are
> using an older RDKit version where maybe the PDB parser was not assigning
> formal charges based on connectivity? Then you would have to manually set
> the formal charge of the protonated nitrogen atom (# 2 N) to +1.
>
> Cheers,
> p.
>
> On 08/22/18 19:06, Phuong Chau wrote:
>
> Hi Paolo,
>
> I got this new error when I tried with different chemicals:
>
> python 3DAlignmentwith3OAwithGRO.py CCNC\(N\)N_afterMD.pdb
> CC\[NH2\+\]CCN_Ligpargen.pdb
> [12:02:45] Explicit valence for atom # 2 N, 4, is greater than permitted
> Traceback (most recent call last):
> File "3DAlignmentwith3OAwithGRO.py", line 12, in
> rdDistGeom.EmbedMolecule(prbMolwithH)
> Boost.Python.ArgumentError: Python argument types in
> rdkit.Chem.rdDistGeom.EmbedMolecule(NoneType)
> did not match C++ signature:
> EmbedMolecule(RDKit::ROMol {lvalue} mol, unsigned int maxAttempts=0,
> int randomSeed=-1, bool clearConfs=True, bool useRandomCoords=False, double
> boxSizeMult=2.0, bool randNegEig=True, unsigned int numZeroFail=1,
> boost::python::dict {lvalue} coordMap={}, double forceTol=0.001, bool
> ignoreSmoothingFailures=False)
>
> The second argument I also downloaded it from Ligpargen but it gave me
> this error. I am not so sure why it happened.
>
> I also attached the two files :
>
> On Tue, Aug 21, 2018 at 3:19 PM, Paolo Tosco
> wrote:
>
>> Hi Phuong,
>>
>> it does have hydrogens after the alignment:
>>
>> $ ls
>> UNK_A1C198.pdb align.py drug.pdb
>>
>> $ python align.py drug.pdb UNK_A1C198.pdb
>>
>> $ ls
>> UNK.pdb UNK_A1C198.pdb align.py drug.pdb
>> $ cat UNK.pdb
>> ATOM 1 C00 UNK 1 37.884 56.016 51.678 1.00
>> 0.00 C
>> ATOM 2 C01 UNK 1 37.501 54.554 51.968 1.00
>> 0.00 C
>> ATOM 3 C02 UNK 1 36.085 54.316 52.531 1.00
>> 0.00 C
>> ATOM 4 C03 UNK 1 35.697 55.221 53.715 1.00
>> 0.00 C
>> ATOM 5 C04 UNK 1 34.907 56.485 53.304 1.00
>> 0.00 C
>> ATOM 6 C05 UNK 1 35.788 57.733 53.108 1.00
>> 0.00 C
>> ATOM 7 C06 UNK 1 36.227 57.998 51.654 1.00
>> 0.00 C
>> ATOM 8 C07 UNK 1 36.813 56.799 50.889 1.00
>> 0.00 C
>> ATOM 9 H08 UNK 1 38.821 56.001 51.080 1.00
>> 0.00 H
>> ATOM 10 H09 UNK 1 38.160 56.522 52.625 1.00
>> 0.00 H
>> ATOM 11 H0A UNK 1 37.598 53.969 51.027 1.00
>> 0.00 H
>> ATOM 12 H0B UNK 1 38.243 54.138 52.684 1.00
>> 0.00 H
>> ATOM 13 H0C UNK 1 35.325 54.368 51.723 1.00
>> 0.00 H
>> ATOM 14 H0D UNK 1 36.057 53.262 52.884 1.00
>> 0.00 H
>> ATOM 15 H0E UNK 1 35.031 54.622 54.375 1.00
>> 0.00 H
>> ATOM 16 H0F UNK 1 36.586 55.471 54.333 1.00
>> 0.00 H
>> ATOM 17 H0G UNK 1 34.203 56.712 54.134 1.00
>> 0.00 H
>> ATOM 18 H0H UNK 1 34.270 56.304 52.411 1.00
>> 0.00 H
>> ATOM 19 H0I UNK 1 35.198 58.620 53.428 1.00
>> 0.00 H
>> ATOM 20 H0J UNK 1 36.665 57.714 53.788 1.00
>> 0.00 H
>> ATOM 21 H0K UNK 1 36.982 58.814 51.667 1.00
>> 0.00 H
>> ATOM 22 H0M UNK 1 35.352 58.378 51.082 1.00
>> 0.00 H
>> ATOM 23 H0N UNK 1 35.991 56.144 50.537 1.00
>> 0.00 H
>> ATOM 24 H0O UNK 1 37.280 57.194 49.960 1.00
>> 0.00 H
>> CONECT1289 10
>> CONECT23 11 12
>> CONECT34 13 14
>> CONECT45 15 16
>> CONECT56 17 18
>> CONECT67 19 20
>> CONECT78 21 22
>> CONECT8 23 24
>> END
>>
>> p.
>>
>>
>> On 08/21/18 20:46, Phuong Chau wrote:
>>
>> Hi Paolo,
>>
>> I tried the code, it worked with the alignment but the new chemical after
>> alignment still does not have any H atoms. If I add H after alignment, it
>> has weird structure with H atoms. Would you please show me how to fix this
>> problem?
>>
>> Best,
>> Phuong Chau
>>
>> On Tue, Aug 21, 2018 at 8:01 AM, Paolo Tosco
>> wrote:
>>
>>> Hi Phuong,
>>>
>>> If you wish to retain Hs you just need to set removeHs = False when you
>>> call MolFromPDBFile():
>>>
>>> # align.py
>>>
>>> import sys
>>> from rdkit import Chem
>>> from rdkit.Chem import rdDistGeom, rdMolAlign, rdForceFieldHelpers
>>> from rdkit import Chem
>>>
>>>
Re: [Rdkit-discuss] Alignment using LIgpargen file
Dear Phuong,
it works for me with the latest RDKit release. I am wondering if you are
using an older RDKit version where maybe the PDB parser was not
assigning formal charges based on connectivity? Then you would have to
manually set the formal charge of the protonated nitrogen atom (# 2 N)
to +1.
Cheers,
p.
On 08/22/18 19:06, Phuong Chau wrote:
Hi Paolo,
I got this new error when I tried with different chemicals:
python 3DAlignmentwith3OAwithGRO.py CCNC\(N\)N_afterMD.pdb
CC\[NH2\+\]CCN_Ligpargen.pdb
[12:02:45] Explicit valence for atom # 2 N, 4, is greater than permitted
Traceback (most recent call last):
File "3DAlignmentwith3OAwithGRO.py", line 12, in
rdDistGeom.EmbedMolecule(prbMolwithH)
Boost.Python.ArgumentError: Python argument types in
rdkit.Chem.rdDistGeom.EmbedMolecule(NoneType)
did not match C++ signature:
EmbedMolecule(RDKit::ROMol {lvalue} mol, unsigned int
maxAttempts=0, int randomSeed=-1, bool clearConfs=True, bool
useRandomCoords=False, double boxSizeMult=2.0, bool randNegEig=True,
unsigned int numZeroFail=1, boost::python::dict {lvalue} coordMap={},
double forceTol=0.001, bool ignoreSmoothingFailures=False)
The second argument I also downloaded it from Ligpargen but it gave me
this error. I am not so sure why it happened.
I also attached the two files :
On Tue, Aug 21, 2018 at 3:19 PM, Paolo Tosco
mailto:paolo.tosco.m...@gmail.com>> wrote:
Hi Phuong,
it does have hydrogens after the alignment:
$ ls
UNK_A1C198.pdb align.py drug.pdb
$ python align.py drug.pdb UNK_A1C198.pdb
$ ls
UNK.pdb UNK_A1C198.pdb align.py drug.pdb
$ cat UNK.pdb
ATOM 1 C00 UNK 1 37.884 56.016 51.678 1.00
0.00 C
ATOM 2 C01 UNK 1 37.501 54.554 51.968 1.00
0.00 C
ATOM 3 C02 UNK 1 36.085 54.316 52.531 1.00
0.00 C
ATOM 4 C03 UNK 1 35.697 55.221 53.715 1.00
0.00 C
ATOM 5 C04 UNK 1 34.907 56.485 53.304 1.00
0.00 C
ATOM 6 C05 UNK 1 35.788 57.733 53.108 1.00
0.00 C
ATOM 7 C06 UNK 1 36.227 57.998 51.654 1.00
0.00 C
ATOM 8 C07 UNK 1 36.813 56.799 50.889 1.00
0.00 C
ATOM 9 H08 UNK 1 38.821 56.001 51.080 1.00
0.00 H
ATOM 10 H09 UNK 1 38.160 56.522 52.625 1.00
0.00 H
ATOM 11 H0A UNK 1 37.598 53.969 51.027 1.00
0.00 H
ATOM 12 H0B UNK 1 38.243 54.138 52.684 1.00
0.00 H
ATOM 13 H0C UNK 1 35.325 54.368 51.723 1.00
0.00 H
ATOM 14 H0D UNK 1 36.057 53.262 52.884 1.00
0.00 H
ATOM 15 H0E UNK 1 35.031 54.622 54.375 1.00
0.00 H
ATOM 16 H0F UNK 1 36.586 55.471 54.333 1.00
0.00 H
ATOM 17 H0G UNK 1 34.203 56.712 54.134 1.00
0.00 H
ATOM 18 H0H UNK 1 34.270 56.304 52.411 1.00
0.00 H
ATOM 19 H0I UNK 1 35.198 58.620 53.428 1.00
0.00 H
ATOM 20 H0J UNK 1 36.665 57.714 53.788 1.00
0.00 H
ATOM 21 H0K UNK 1 36.982 58.814 51.667 1.00
0.00 H
ATOM 22 H0M UNK 1 35.352 58.378 51.082 1.00
0.00 H
ATOM 23 H0N UNK 1 35.991 56.144 50.537 1.00
0.00 H
ATOM 24 H0O UNK 1 37.280 57.194 49.960 1.00
0.00 H
CONECT 1 2 8 9 10
CONECT 2 3 11 12
CONECT 3 4 13 14
CONECT 4 5 15 16
CONECT 5 6 17 18
CONECT 6 7 19 20
CONECT 7 8 21 22
CONECT 8 23 24
END
p.
On 08/21/18 20:46, Phuong Chau wrote:
Hi Paolo,
I tried the code, it worked with the alignment but the new
chemical after alignment still does not have any H atoms. If I
add H after alignment, it has weird structure with H atoms. Would
you please show me how to fix this problem?
Best,
Phuong Chau
On Tue, Aug 21, 2018 at 8:01 AM, Paolo Tosco
mailto:paolo.tosco.m...@gmail.com>>
wrote:
Hi Phuong,
If you wish to retain Hs you just need to set removeHs =
False when you call MolFromPDBFile():
# align.py
import sys
from rdkit import Chem
from rdkit.Chem import rdDistGeom, rdMolAlign,
rdForceFieldHelpers
from rdkit import Chem
refMolwithH = Chem.MolFromPDBFile(sys.argv[1], removeHs = False)
s = sys.argv[2]
prbMolwithH = Chem.MolFromPDBFile(s, removeHs = False)
idx=s.find('_')
chemB= s[:idx]
rdDistGeom.EmbedMolecule(prbMolwithH)
Re: [Rdkit-discuss] Alignment using LIgpargen file
Hi Paolo,
I got this new error when I tried with different chemicals:
python 3DAlignmentwith3OAwithGRO.py CCNC\(N\)N_afterMD.pdb
CC\[NH2\+\]CCN_Ligpargen.pdb
[12:02:45] Explicit valence for atom # 2 N, 4, is greater than permitted
Traceback (most recent call last):
File "3DAlignmentwith3OAwithGRO.py", line 12, in
rdDistGeom.EmbedMolecule(prbMolwithH)
Boost.Python.ArgumentError: Python argument types in
rdkit.Chem.rdDistGeom.EmbedMolecule(NoneType)
did not match C++ signature:
EmbedMolecule(RDKit::ROMol {lvalue} mol, unsigned int maxAttempts=0,
int randomSeed=-1, bool clearConfs=True, bool useRandomCoords=False, double
boxSizeMult=2.0, bool randNegEig=True, unsigned int numZeroFail=1,
boost::python::dict {lvalue} coordMap={}, double forceTol=0.001, bool
ignoreSmoothingFailures=False)
The second argument I also downloaded it from Ligpargen but it gave me this
error. I am not so sure why it happened.
I also attached the two files :
On Tue, Aug 21, 2018 at 3:19 PM, Paolo Tosco
wrote:
> Hi Phuong,
>
> it does have hydrogens after the alignment:
>
> $ ls
> UNK_A1C198.pdb align.py drug.pdb
>
> $ python align.py drug.pdb UNK_A1C198.pdb
>
> $ ls
> UNK.pdb UNK_A1C198.pdb align.py drug.pdb
> $ cat UNK.pdb
> ATOM 1 C00 UNK 1 37.884 56.016 51.678 1.00
> 0.00 C
> ATOM 2 C01 UNK 1 37.501 54.554 51.968 1.00
> 0.00 C
> ATOM 3 C02 UNK 1 36.085 54.316 52.531 1.00
> 0.00 C
> ATOM 4 C03 UNK 1 35.697 55.221 53.715 1.00
> 0.00 C
> ATOM 5 C04 UNK 1 34.907 56.485 53.304 1.00
> 0.00 C
> ATOM 6 C05 UNK 1 35.788 57.733 53.108 1.00
> 0.00 C
> ATOM 7 C06 UNK 1 36.227 57.998 51.654 1.00
> 0.00 C
> ATOM 8 C07 UNK 1 36.813 56.799 50.889 1.00
> 0.00 C
> ATOM 9 H08 UNK 1 38.821 56.001 51.080 1.00
> 0.00 H
> ATOM 10 H09 UNK 1 38.160 56.522 52.625 1.00
> 0.00 H
> ATOM 11 H0A UNK 1 37.598 53.969 51.027 1.00
> 0.00 H
> ATOM 12 H0B UNK 1 38.243 54.138 52.684 1.00
> 0.00 H
> ATOM 13 H0C UNK 1 35.325 54.368 51.723 1.00
> 0.00 H
> ATOM 14 H0D UNK 1 36.057 53.262 52.884 1.00
> 0.00 H
> ATOM 15 H0E UNK 1 35.031 54.622 54.375 1.00
> 0.00 H
> ATOM 16 H0F UNK 1 36.586 55.471 54.333 1.00
> 0.00 H
> ATOM 17 H0G UNK 1 34.203 56.712 54.134 1.00
> 0.00 H
> ATOM 18 H0H UNK 1 34.270 56.304 52.411 1.00
> 0.00 H
> ATOM 19 H0I UNK 1 35.198 58.620 53.428 1.00
> 0.00 H
> ATOM 20 H0J UNK 1 36.665 57.714 53.788 1.00
> 0.00 H
> ATOM 21 H0K UNK 1 36.982 58.814 51.667 1.00
> 0.00 H
> ATOM 22 H0M UNK 1 35.352 58.378 51.082 1.00
> 0.00 H
> ATOM 23 H0N UNK 1 35.991 56.144 50.537 1.00
> 0.00 H
> ATOM 24 H0O UNK 1 37.280 57.194 49.960 1.00
> 0.00 H
> CONECT1289 10
> CONECT23 11 12
> CONECT34 13 14
> CONECT45 15 16
> CONECT56 17 18
> CONECT67 19 20
> CONECT78 21 22
> CONECT8 23 24
> END
>
> p.
>
>
> On 08/21/18 20:46, Phuong Chau wrote:
>
> Hi Paolo,
>
> I tried the code, it worked with the alignment but the new chemical after
> alignment still does not have any H atoms. If I add H after alignment, it
> has weird structure with H atoms. Would you please show me how to fix this
> problem?
>
> Best,
> Phuong Chau
>
> On Tue, Aug 21, 2018 at 8:01 AM, Paolo Tosco
> wrote:
>
>> Hi Phuong,
>>
>> If you wish to retain Hs you just need to set removeHs = False when you
>> call MolFromPDBFile():
>>
>> # align.py
>>
>> import sys
>> from rdkit import Chem
>> from rdkit.Chem import rdDistGeom, rdMolAlign, rdForceFieldHelpers
>> from rdkit import Chem
>>
>> refMolwithH = Chem.MolFromPDBFile(sys.argv[1], removeHs = False)
>> s = sys.argv[2]
>> prbMolwithH = Chem.MolFromPDBFile(s, removeHs = False)
>> idx=s.find('_')
>> chemB= s[:idx]
>>
>> rdDistGeom.EmbedMolecule(prbMolwithH)
>> rdForceFieldHelpers.UFFOptimizeMolecule(prbMolwithH)
>>
>> ##Alignment
>> pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)
>> score = pyO3A.Align()
>>
>> ##3D coords of Chem B after alignmnet
>> Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))
>>
>> and then
>>
>> ./align.py drug.pdb UNK_A1C198.pdb
>> Cheers,
>> p.
>>
>> On 08/20/18 19:02, Phuong Chau wrote:
>>
>> Hello everyone,
>>
>> I am trying to align two chemicals using their pdb files with the
>> following script:
>>
>> *refMolwithH = Chem.MolFromPDBFile(sys.argv[1])*
>> *s = sys.argv[2]*
>> *prbMolwithH = Chem.MolFromPDBFile(s)*
>> *idx=s.find('_')*
>> *chemB= s[:idx]*
>>
>>
Re: [Rdkit-discuss] Alignment using LIgpargen file
Hi Phuong,
it does have hydrogens after the alignment:
$ ls
UNK_A1C198.pdb align.py drug.pdb
$ python align.py drug.pdb UNK_A1C198.pdb
$ ls
UNK.pdb UNK_A1C198.pdb align.py drug.pdb
$ cat UNK.pdb
ATOM 1 C00 UNK 1 37.884 56.016 51.678 1.00
0.00 C
ATOM 2 C01 UNK 1 37.501 54.554 51.968 1.00
0.00 C
ATOM 3 C02 UNK 1 36.085 54.316 52.531 1.00
0.00 C
ATOM 4 C03 UNK 1 35.697 55.221 53.715 1.00
0.00 C
ATOM 5 C04 UNK 1 34.907 56.485 53.304 1.00
0.00 C
ATOM 6 C05 UNK 1 35.788 57.733 53.108 1.00
0.00 C
ATOM 7 C06 UNK 1 36.227 57.998 51.654 1.00
0.00 C
ATOM 8 C07 UNK 1 36.813 56.799 50.889 1.00
0.00 C
ATOM 9 H08 UNK 1 38.821 56.001 51.080 1.00
0.00 H
ATOM 10 H09 UNK 1 38.160 56.522 52.625 1.00
0.00 H
ATOM 11 H0A UNK 1 37.598 53.969 51.027 1.00
0.00 H
ATOM 12 H0B UNK 1 38.243 54.138 52.684 1.00
0.00 H
ATOM 13 H0C UNK 1 35.325 54.368 51.723 1.00
0.00 H
ATOM 14 H0D UNK 1 36.057 53.262 52.884 1.00
0.00 H
ATOM 15 H0E UNK 1 35.031 54.622 54.375 1.00
0.00 H
ATOM 16 H0F UNK 1 36.586 55.471 54.333 1.00
0.00 H
ATOM 17 H0G UNK 1 34.203 56.712 54.134 1.00
0.00 H
ATOM 18 H0H UNK 1 34.270 56.304 52.411 1.00
0.00 H
ATOM 19 H0I UNK 1 35.198 58.620 53.428 1.00
0.00 H
ATOM 20 H0J UNK 1 36.665 57.714 53.788 1.00
0.00 H
ATOM 21 H0K UNK 1 36.982 58.814 51.667 1.00
0.00 H
ATOM 22 H0M UNK 1 35.352 58.378 51.082 1.00
0.00 H
ATOM 23 H0N UNK 1 35.991 56.144 50.537 1.00
0.00 H
ATOM 24 H0O UNK 1 37.280 57.194 49.960 1.00
0.00 H
CONECT 1 2 8 9 10
CONECT 2 3 11 12
CONECT 3 4 13 14
CONECT 4 5 15 16
CONECT 5 6 17 18
CONECT 6 7 19 20
CONECT 7 8 21 22
CONECT 8 23 24
END
p.
On 08/21/18 20:46, Phuong Chau wrote:
Hi Paolo,
I tried the code, it worked with the alignment but the new chemical
after alignment still does not have any H atoms. If I add H after
alignment, it has weird structure with H atoms. Would you please show
me how to fix this problem?
Best,
Phuong Chau
On Tue, Aug 21, 2018 at 8:01 AM, Paolo Tosco
mailto:paolo.tosco.m...@gmail.com>> wrote:
Hi Phuong,
If you wish to retain Hs you just need to set removeHs = False
when you call MolFromPDBFile():
# align.py
import sys
from rdkit import Chem
from rdkit.Chem import rdDistGeom, rdMolAlign, rdForceFieldHelpers
from rdkit import Chem
refMolwithH = Chem.MolFromPDBFile(sys.argv[1], removeHs = False)
s = sys.argv[2]
prbMolwithH = Chem.MolFromPDBFile(s, removeHs = False)
idx=s.find('_')
chemB= s[:idx]
rdDistGeom.EmbedMolecule(prbMolwithH)
rdForceFieldHelpers.UFFOptimizeMolecule(prbMolwithH)
##Alignment
pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)
score = pyO3A.Align()
##3D coords of Chem B after alignmnet
Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))
and then
./align.py drug.pdb UNK_A1C198.pdb
Cheers,
p.
On 08/20/18 19:02, Phuong Chau wrote:
Hello everyone,
I am trying to align two chemicals using their pdb files with the
following script:
*refMolwithH = Chem.MolFromPDBFile(sys.argv[1])*
*s = sys.argv[2]*
*prbMolwithH = Chem.MolFromPDBFile(s)*
*idx=s.find('_')*
*chemB= s[:idx]*
*
*
*rdDistGeom.EmbedMolecule(prbMolwithH)*
*AllChem.UFFOptimizeMolecule(prbMolwithH)*
*
*
*##Alignment*
*pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)*
*score = pyO3A.Align()*
*
*
*##3D coords of Chem B after alignmnet*
*Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))*
The probe chemical pdb is generated from GRO file generated from
LigParGen browser by this script:
*gmx editconf -f input.gro -o output.pdb*
The problem is the new aligned chemical is not aligned with the
refMol but with probMol and the new chemical does not have H
atoms on it.
Would you please help me with this problem? Are there other ways
to align the pdb files generated from LIgParGen?
Thank you so much for your help
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org!http://sdm.link/slashdot
___
Rdkit-discuss mailing list
Re: [Rdkit-discuss] Alignment using LIgpargen file
Hi Paolo,
I tried the code, it worked with the alignment but the new chemical after
alignment still does not have any H atoms. If I add H after alignment, it
has weird structure with H atoms. Would you please show me how to fix this
problem?
Best,
Phuong Chau
On Tue, Aug 21, 2018 at 8:01 AM, Paolo Tosco
wrote:
> Hi Phuong,
>
> If you wish to retain Hs you just need to set removeHs = False when you
> call MolFromPDBFile():
>
> # align.py
>
> import sys
> from rdkit import Chem
> from rdkit.Chem import rdDistGeom, rdMolAlign, rdForceFieldHelpers
> from rdkit import Chem
>
> refMolwithH = Chem.MolFromPDBFile(sys.argv[1], removeHs = False)
> s = sys.argv[2]
> prbMolwithH = Chem.MolFromPDBFile(s, removeHs = False)
> idx=s.find('_')
> chemB= s[:idx]
>
> rdDistGeom.EmbedMolecule(prbMolwithH)
> rdForceFieldHelpers.UFFOptimizeMolecule(prbMolwithH)
>
> ##Alignment
> pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)
> score = pyO3A.Align()
>
> ##3D coords of Chem B after alignmnet
> Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))
>
> and then
>
> ./align.py drug.pdb UNK_A1C198.pdb
> Cheers,
> p.
>
> On 08/20/18 19:02, Phuong Chau wrote:
>
> Hello everyone,
>
> I am trying to align two chemicals using their pdb files with the
> following script:
>
> *refMolwithH = Chem.MolFromPDBFile(sys.argv[1])*
> *s = sys.argv[2]*
> *prbMolwithH = Chem.MolFromPDBFile(s)*
> *idx=s.find('_')*
> *chemB= s[:idx]*
>
> *rdDistGeom.EmbedMolecule(prbMolwithH)*
> *AllChem.UFFOptimizeMolecule(prbMolwithH)*
>
> *##Alignment*
> *pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)*
> *score = pyO3A.Align()*
>
> *##3D coords of Chem B after alignmnet*
> *Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))*
>
> The probe chemical pdb is generated from GRO file generated from LigParGen
> browser by this script:
> *gmx editconf -f input.gro -o output.pdb*
>
> The problem is the new aligned chemical is not aligned with the refMol but
> with probMol and the new chemical does not have H atoms on it.
>
> Would you please help me with this problem? Are there other ways to align
> the pdb files generated from LIgParGen?
>
> Thank you so much for your help
>
>
>
> --
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
>
>
>
> ___
> Rdkit-discuss mailing
> listRdkit-discuss@lists.sourceforge.nethttps://lists.sourceforge.net/lists/listinfo/rdkit-discuss
>
>
>
--
Phuong Chau
Smith College '20
Engineering Major
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot___
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
Re: [Rdkit-discuss] Alignment using LIgpargen file
Hi Phuong,
If you wish to retain Hs you just need to set removeHs = False when you
call MolFromPDBFile():
# align.py
import sys
from rdkit import Chem
from rdkit.Chem import rdDistGeom, rdMolAlign, rdForceFieldHelpers
from rdkit import Chem
refMolwithH = Chem.MolFromPDBFile(sys.argv[1], removeHs = False)
s = sys.argv[2]
prbMolwithH = Chem.MolFromPDBFile(s, removeHs = False)
idx=s.find('_')
chemB= s[:idx]
rdDistGeom.EmbedMolecule(prbMolwithH)
rdForceFieldHelpers.UFFOptimizeMolecule(prbMolwithH)
##Alignment
pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)
score = pyO3A.Align()
##3D coords of Chem B after alignmnet
Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))
and then
./align.py drug.pdb UNK_A1C198.pdb
Cheers,
p.
On 08/20/18 19:02, Phuong Chau wrote:
Hello everyone,
I am trying to align two chemicals using their pdb files with the
following script:
*refMolwithH = Chem.MolFromPDBFile(sys.argv[1])*
*s = sys.argv[2]*
*prbMolwithH = Chem.MolFromPDBFile(s)*
*idx=s.find('_')*
*chemB= s[:idx]*
*
*
*rdDistGeom.EmbedMolecule(prbMolwithH)*
*AllChem.UFFOptimizeMolecule(prbMolwithH)*
*
*
*##Alignment*
*pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)*
*score = pyO3A.Align()*
*
*
*##3D coords of Chem B after alignmnet*
*Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))*
The probe chemical pdb is generated from GRO file generated from
LigParGen browser by this script:
*gmx editconf -f input.gro -o output.pdb*
The problem is the new aligned chemical is not aligned with the refMol
but with probMol and the new chemical does not have H atoms on it.
Would you please help me with this problem? Are there other ways to
align the pdb files generated from LIgParGen?
Thank you so much for your help
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
___
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot___
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
Re: [Rdkit-discuss] Alignment using LIgpargen file
Hi Phuong,
could you please send me the PDB files you are trying to align? You may
reply to me directly.
Cheers,
p.
On 08/20/18 19:02, Phuong Chau wrote:
Hello everyone,
I am trying to align two chemicals using their pdb files with the
following script:
*refMolwithH = Chem.MolFromPDBFile(sys.argv[1])*
*s = sys.argv[2]*
*prbMolwithH = Chem.MolFromPDBFile(s)*
*idx=s.find('_')*
*chemB= s[:idx]*
*
*
*rdDistGeom.EmbedMolecule(prbMolwithH)*
*AllChem.UFFOptimizeMolecule(prbMolwithH)*
*
*
*##Alignment*
*pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)*
*score = pyO3A.Align()*
*
*
*##3D coords of Chem B after alignmnet*
*Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))*
The probe chemical pdb is generated from GRO file generated from
LigParGen browser by this script:
*gmx editconf -f input.gro -o output.pdb*
The problem is the new aligned chemical is not aligned with the refMol
but with probMol and the new chemical does not have H atoms on it.
Would you please help me with this problem? Are there other ways to
align the pdb files generated from LIgParGen?
Thank you so much for your help
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
___
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot___
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
[Rdkit-discuss] Alignment using LIgpargen file
Hello everyone,
I am trying to align two chemicals using their pdb files with the following
script:
*refMolwithH = Chem.MolFromPDBFile(sys.argv[1])*
*s = sys.argv[2]*
*prbMolwithH = Chem.MolFromPDBFile(s)*
*idx=s.find('_')*
*chemB= s[:idx]*
*rdDistGeom.EmbedMolecule(prbMolwithH)*
*AllChem.UFFOptimizeMolecule(prbMolwithH)*
*##Alignment*
*pyO3A = rdMolAlign.GetO3A(prbMolwithH, refMolwithH)*
*score = pyO3A.Align()*
*##3D coords of Chem B after alignmnet*
*Chem.MolToPDBFile(prbMolwithH,'{}.pdb'.format(chemB))*
The probe chemical pdb is generated from GRO file generated from LigParGen
browser by this script:
*gmx editconf -f input.gro -o output.pdb*
The problem is the new aligned chemical is not aligned with the refMol but
with probMol and the new chemical does not have H atoms on it.
Would you please help me with this problem? Are there other ways to align
the pdb files generated from LIgParGen?
Thank you so much for your help
--
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot___
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
