Comments (10)
Thanks for reporting. I can't see the attachment containing the ligands. Were you meant to attach, or are they available from the paper link?
from biosimspace.
Hi Lester, Sorry I forgot to attach the zip file...
issue.tar.gz
from biosimspace.
hi @jasmin-guven , can you also share the BSS script you used to generate a merged molecule so we can check the syntax used to generate atom mappings ?
It would be good to use https://github.com/OpenBioSim/biosimspace/issues for reporting future issues
from biosimspace.
I can't even parameterise the ligands using both GAFF or OpenFF, so can't perform the merge from the SDF files alone. That said, I've checked the mappings manually and all involving ligand11 include the P
element in the mapping, e.g. mapping atom index 15 in ligand11 to atom index 16 in ligand5 and ligand4.
mapping_11_5
Out[23]:
{0: 0,
1: 1,
2: 2,
3: 3,
4: 4,
5: 5,
6: 6,
7: 7,
8: 8,
9: 9,
10: 11,
11: 12,
12: 13,
13: 14,
14: 15,
15: 16,
16: 18,
17: 19,
18: 22,
19: 23,
20: 24,
21: 25,
22: 26,
23: 10}
mapping_11_4
Out[24]:
{0: 1,
1: 2,
2: 3,
3: 4,
4: 5,
5: 6,
6: 7,
7: 8,
8: 9,
9: 10,
10: 11,
11: 12,
12: 13,
13: 14,
14: 15,
15: 16,
16: 20,
17: 21,
18: 22,
19: 23,
20: 24,
21: 25,
22: 0,
23: 26}
from biosimspace.
hi @jasmin-guven , can you also share the BSS script you used to generate a merged molecule so we can check the syntax used to generate atom mappings ? It would be good to use https://github.com/OpenBioSim/biosimspace/issues for reporting future issues
Hi @jmichel80,
Here's my code:
https://github.com/meyresearch/metalloenzymes/blob/main/scripts/prepafe.py
And sorry, will report these in the above repo in the future!
from biosimspace.
No problem. Could you share the equlibrated prm/rst files for the ligands as well so that I can try to merge. It may be that the mapping is fine with the original SDF file, but is failing for the topology loaded from AMBER format. We now have a direct RDKit conversion, which should make things like this more robust.
from biosimspace.
Hi, @lohedges
Here are equilibrated prm7 and rst7 files for ligand 11 and ligand 3, 4, 12 and 14.
ligand_3.zip
ligand_11.zip
ligand_4.zip
ligand_14.zip
ligand_12.zip
ligand_5.zip
Thanks!
from biosimspace.
Hi again. I've tested locally and with the latest BioSimSpace from OpenBioSim I get the same mapping using either the SDF input or the PR7/RST7 AMBER files, which doesn't have dummies for the phosphonate group. I would expect this to be the case since we now use the native Sire to RDKit conversion behind the scenes, which has some tools to recover stereochemistry.
Prior to this addition, conversion to RDKit format went via an intermediate file format, so information could be lost depending on the format used. If molecules were loaded from SDF, then we would try use that format as the intermediate in order to preserve stereochemistry. (Which assumes that the original SDF file was well formatted.) If it wasn't loaded from SDF, then we go via PDB format and hope that RDKit can recover missing information during sanitisation. In your script you are performing the merge using the ligands loaded from AMBER format files, not the original SDFs, so are going via PDB format. I can only assume that some of the ligands aren't being loaded into RDKit in an optimal way, hence the strange MCS mappings. (Perhaps the equilibrated coordinates are sub-optimal, or some stereochemistry is lost.)
I'd advise the following:
- Try with the current BioSimSpace from OpenBioSim:
mamba create -n openbiosim -c conda-forge -c openbiosim/label/dev biosimspace ambertools gromacs
mamba activate openbiosim
(You can replace mamba
with conda
if you don't use it.)
- Try changing your script to load the SDF ligands to perform the MCS mappings.
from biosimspace.
For reference, this is what I see using your AMBER files with the latest code:
The two phosphonate groups map onto each other. The only things that don't map are:
- A hydrogen attached to atom index 10 in ligand11 becomes a dummy in the lambda = 1 end state.
- Three hydrogens attached to atom index 10 in ligand5 become dummies in the lambda = 0 end state.
from biosimspace.
Hi Lester,
Sorry for my delayed response.
I tried equilibrating the systems again and now it seems I don't get weird mappings. I also installed openbiosimspace, but was met with other errors along my workflow (which I will try to raise issues for on the correct site, if they persist).
If I get this error again, I will probably try using the SDF files to set up the mappings, unless I get my openbiosimspace to work.
Thanks again for helping, I will close this issue now.
from biosimspace.
Related Issues (20)
- Storage quota exceeded on the Anaconda cloud HOT 1
- HILLS file not updating during restarting simulation HOT 4
- Is it possible to add ions to an already solvated system, without adding any more water molecules? HOT 4
- Problem when converting Gro/Top to Rst7 HOT 12
- Making HMR work with dummies HOT 11
- Make _removeDummies() recognise dummy bonds HOT 10
- Would it be possible to include HiMap as an alternative to LoMap in the generateNetwork method? HOT 2
- Unable to import BioSimSpace after mamba installation HOT 12
- BioSimSpace.Align.viewMapping: py3Dmol.view not rendering in Jupyter Notebook HOT 11
- `_toRegularMolecule` failing with `convert_amber_dummies=True` HOT 15
- Absolute Binding Free Energy Calculations with GROMACS HOT 7
- Pickling a parametrised molecule from an SDF file with GAFF2 HOT 4
- Are atomtypes meant to be updated when using repartitionHydrogenMass? HOT 16
- [BUG] OSError "It looks like you failed to include a topology file." when reading SDF which doesn't contain redundant bond order information HOT 10
- [BUG] Cannot import BioSimSpace after installation of biosimspace 2023.1.2 using mamba HOT 3
- AnalysisError: SOMD free-energy analysis failed! HOT 5
- Reporting a vulnerability
- [BUG] Water molecule coordinates are setting to (0,0,0) HOT 4
- Issue in reproducing the funnel methadynamics tutorial HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from biosimspace.