Comments (9)
Do you have the unmodified FESetup output files (rather than the ones written by BioSImSpace). If those don't run, then it's an issue with the configurations, rather than BioSImSpace (although we might be able to modify the protocol to help). I usually see crashes like this when the system hasn't been minimised correctly.
from biosimspace.
For reference, sander can't even perform a miminisation using these input files (at least with SHAKE on, which is what I was told we want to use).
from biosimspace.
It looks like the original configuration is a bunch of junk (I don't know what FESetup uses for minimisation/equilibration). If I disable SHAKE and re-minimise the system using BioSimSpace, then production MD runs just fine.
from biosimspace.
If FESetup is doing minimization/equilibration its default is to use amber16/sander although you can use other programs if you specify the mdengine in the input file. However, for each of the three categories of molecules (ligand, protein, and complex) the minimization/equilibration must be set separately, i.e. if you minimize/equilibrate the ligand or protein it does not automatically minimize/equilibrate the complex. fesetup/_complexes/*/solvated.rst is the input for minimization not the result.
from biosimspace.
from biosimspace.
Thanks for the clear explanation of the protocol, that's interesting.
Sorry for any confusion, I wasn't trying to imply that you can just minimise with BioSimSpace and be good to go, only that minimisation fixes whatever steric clashes are left after the FESetup step. I just wanted to do the simplest thing to see if I could then run any protocol using the re-minimised system (since the FESetup configurations crashed within the first few iterations).
from biosimspace.
If I am reading the above correctly, it sounds like we can't rely on SOMD, i.e. OpenMM, to correctly minimise or equilibrate the system. At the moment I'm making using of the minimisation and equilibration options in SOMD during my free energy simulations, since it is currently the only MD engine that can run dynamics with a system containing perturbed molecules.
For the purposes of the D3R challenge I could hack the other engines to write atom properties at lambda = 0, then copy the updated coordinates back into the original system following minimisation/equilibration, but this would take a bit of work. (I discussed something like this with Christopher last week.)
from biosimspace.
from biosimspace.
Great, this is pretty much what I'm doing anyway. It sounds like I'll need to add another keyword option to my protocols to restrain all non-solvent atoms. I'll also need to deal with writing input files for the lambda=0 state (using the property map) then just re-adding the updated solvent coordinates back into the original system.
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
- [BUG] HOT 10
- 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.