You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I have a problem with the systemPrepare function since I installed the new version. When I switched to another environment which still had an older version of HTMD installed the error did not occur. Any help would be very much appreciated.:)
The script I am usig for building:
fromhtmd.uiimport*#Python based env to prepare, simulate systemsfrommoleculekit.utilimportmaxDistance#Python library for manipulation of structures. importglobimportosprotein=Molecule('../proteins/1jh9_test.pdb')
protein=systemPrepare(protein,pH=7.4)
protein=autoSegment(protein, sel='protein')
ligand=Molecule('../ligands_preparation/1jh9/HPA.mol2')
ligand.set('segid','L')
ligand.set('resname' , 'MOL')
mol=Molecule(name='combo')
mol.append(protein)
mol.append(ligand)
mol.reps.add(sel='protein', style='NewCartoon', color='Secondary Structure')
mol.reps.add(sel='resname MOL', style='Licorice')
mol.center()
D=maxDistance(mol, 'all') # let's find the maximum distance from 0,0,0print(f'max distance is: {D}')
DW=D+8# We add a buffer of 8 around the peptide smol=solvate(mol, minmax=[[-DW, -DW, -DW], [DW, DW, DW]])
topos_amber=amber.defaultTopo() + ['../ligands_preparation/1jh9/HPA.mol2']
frcmods_amber=amber.defaultParam() + ['../ligands_preparation/1jh9/HPA.frcmod']
bmol_amber=amber.build(smol,outdir='1jh9_DNA_build_test',topo=topos_amber, param=frcmods_amber, saltconc=0.15)
andtheerrormessage:
----Moleculechainreport----ChainA:
Firstresidue: ALA2Finalresidue: ARG340ChainB:
Firstresidue: DA699Finalresidue: DT715----Endofchainreport----Traceback (mostrecentcalllast):
File"/agh/scratch/Boje/MD/diverse_scaffold_generation/building/build_batches.py", line11, in<module>protein=systemPrepare(protein,pH=7.4)
File"/agh/projects/users/Merve/miniconda3/envs/acemd/lib/python3.10/site-packages/moleculekit/tools/preparation.py", line733, insystemPreparedefinition, forcefield=_get_custom_ff(molkit_ff=_molkit_ff)
File"/agh/projects/users/Merve/miniconda3/envs/acemd/lib/python3.10/site-packages/moleculekit/tools/preparation_customres.py", line415, in_get_custom_ffdefinition=get_definitions(aa_path=out_aa)
File"/sw/sci/app/pdb2pqr/3.1.0-focal/lib/python3.8/site-packages/pdb2pqr/io.py", line492, inget_definitionsdefinitions=defns.Definition(
File"/sw/sci/app/pdb2pqr/3.1.0-focal/lib/python3.8/site-packages/pdb2pqr/definitions.py", line157, in__init__regexp=re.compile(patch.applyto).match(name)
The text was updated successfully, but these errors were encountered:
Hey, if you look at the traceback of the error you can see that you have a messed up PYTHONPATH. It's using pdb2pqr from /sw/sci/app/pdb2pqr instead of using the one in /agh/projects/users/Merve/miniconda3/envs/acemd/. That's most definitely a different version from the one we are using in HTMD
Hi,
I have a problem with the systemPrepare function since I installed the new version. When I switched to another environment which still had an older version of HTMD installed the error did not occur. Any help would be very much appreciated.:)
The script I am usig for building:
The text was updated successfully, but these errors were encountered: