Comments (3)
This is what will be done by BioSimSpace when writing the pert file and AMBER prm/rst files as part of the SOMD free energy protocol. I plan on making a copy of the user input molecule, then applying my pert file writer (which has to rename the atoms in order to make them unique). You can look at the _toPertFile
method in BioSimSpace/python/_SireWrappers/_molecue.py
to see how the renaming works.
I'll leave this issue as a reminder and close when the pert file / free energy protocol is fully implemented.
from biosimspace.
On a related note, check also def make_pert_file in
https://github.com/CCPBioSim/fesetup/blob/dev/mutate/topol/pertfile.py
I don't know how BSS merge molecules handle this, but there was some complicated logic in FESetup to make sure that dihedrals and impropers were correctly written in the pert file. The issue is that the same dihedral or improper may be defined with different atom indices in the respective end-states, and care must be taken to avoid accidentally duplicating forcefield terms.
There are also some rules around the values of bond/angles or dihedral/impropers for terms involving dummy atoms that could be controlled with keyword arguments.
from biosimspace.
I believe I handle this correctly in both the Gromacs and pert file writer. Since all of the terms are stored as Sire objects, with corresponding BondIDs, AngleIDs, and DihedralDs, it's easy to check for matching terms in the end-states, e.g. by checking the the ID and its mirror.
I'll take a look through the FESetup logic to see if I'm missing anything. For the Gromacs writer, I believe that my procedure is consistent with everything from the latest manual.
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.