Comments (20)
@nhenriksen just a heads up that I'll work on the mol2
coordinates this afternoon. I'll fork the repo and open a PR to the nieldev
branch when I'm done.
from benchmarksets.
I have a brief comment in the README, although maybe it needs more explanation: " To account for the two possible orientations of the guest within the CD cavity, simulation files with the '-p' suffix indicate that the guest is bound with the polar functional group oriented out of the primary (narrow) face of the CD, whereas the '-s' suffix indicates the guest polar functional group is oriented out of the secondary (wider) face of the CD."
Basically, the host in asymmetric, so there are two obvious orientations for a polar guest to bind.
I'll look into the issue with the missing rst7s ....
I agree that we should update all the coordinates.
from benchmarksets.
From experimental literature, the default expectation is that most guests will prefer the "-s" orientation, although there is not a whole lot of data to support that. We've been doing some NMR in our lab that paints a complicated picture about guest orientation, and for simulations, the ammoniums definitely prefer the "-p" orientation for some force fields. So both orientations will want to be considered ultimately. But for the purposes of parameters, etc., the "-s" orientation is fine.
from benchmarksets.
I just want to thank you guys for the great work and dialogue here. This is exactly the sort of thing I hoped getting more of these types of systems up on GitHub would nucleate -- the ability to reuse them, with knowledge about what should happen/what's important, without having to read 20 papers on the topic. :)
from benchmarksets.
Sounds good!
from benchmarksets.
Yes, please share yaml!!
from benchmarksets.
@davidlmobley I didn't have time to finish this weekend, sorry, but the script should be ready (and uploaded) by the end of the day.
from benchmarksets.
Hi Andrea. Yes, the mol2 files have the expected protonation based on pH 6.9, which was used experimentally. I agree that it would be better to have the mol2 coordinates match the rst7 coordinates. I need to focus on finishing up edits I'm making to the benchmarksets paper first, but I could work on this later this week. If you get to it sooner than that, we can use your files.
from benchmarksets.
Ok, thanks! I'll let you know if/when I'll get to this so that we won't both work on the same thing.
from benchmarksets.
Question: what is the difference between the -p
and -s
version of the prmtop
and rst7
files (apologies if it's written somewhere and I've missed it)? And is there one I should prefer for the single-molecule files coordinates? I've noticed that the following files in cd-set1/
are missing a rst7
companion: acd-s15-s.prmtop
, acd-s11-p.prmtop
, acd-s19-p.prmtop
, acd-s21-p.prmtop
.
I've also noticed that the pdb and sd files have the same coordinates of the mol2. If you agree, @nhenriksen, I'd change the coordinates of those files too to make them consistent.
from benchmarksets.
I have a brief comment in the README
Ah! Sorry, I completely missed it, thanks! I'll arbitrarily pick all the s
as a reference then unless you have a reason to prefer the other (or to pick them 50-50).
from benchmarksets.
The setup simulations for those missing files had an error that evaded my detectors. I've added them now. daa8b2e
from benchmarksets.
I've added them now
Thanks!
Since the host's conformations of the rst7
files after equilibration are all slightly different, I was thinking about using structural alignment to get the new mol2 positions. Would this be ok? I'm not sure I can do much better without the docked poses.
from benchmarksets.
@andrrizzi - have you gotten around to generating "ligand in solution" files for these yet as well? One of my rotation students needs these and I'm trying to sort out whether she should do that (and probably deposit here) or whether that would duplicate work you are already doing for Yank.
from benchmarksets.
have you gotten around to generating "ligand in solution" files for these yet as well?
That's my next step for the Yank calculations I need to run. I plan to set them up during the weekend.
from benchmarksets.
OK, thanks. Then I won't have my rotation student do it. But can you share code?
from benchmarksets.
I was actually planning to use the Yank pipeline to set them up. With these mol2 files, it can handle the whole preparation, so I don't think there will be code to share.
from benchmarksets.
Ah, OK, thanks. So it's a "provide a mol2 file of the ligand and a PDB file of the complex" kind of thing, then? Sorry, I forget there are so many modes for using Yank.
from benchmarksets.
provide a mol2 file of the ligand and a PDB file of the complex
Almost. The input files of the pipeline are only the mol2
files of the single molecules that I've modified here (both host and guests). The pipeline then produces (unminimized/unequilibrated) amber files and pdb files of both the complex and solvent phases that YANK uses to run the calculations. But if your student plans to use YANK, having the YAML script I'll share + the mol2 files should suffice (of course they can use prepared system files if they prefer as you remember).
from benchmarksets.
@andrrizzi - any updates on the setup of your calculations? I will need to get my rotation student going soon, so I should probably start her working on setting some up for her purposes (cross comparing different forms of restraints) if yours aren't coming along.
from benchmarksets.
Related Issues (20)
- Add README.md for host-guest inputs; add additional info
- Update text to reflect host-guest sets added in #22 HOT 1
- Add additional supplementary data, perhaps in markdown files? HOT 5
- Get CD benchmarks in? HOT 3
- Provide suitable README files for CB7, GDCC benchmarks
- Provide bound-state starting structures for hosts HOT 5
- Deposit calculated values for provided files when available HOT 5
- Should we use table titles and not numbers in markdown files? HOT 6
- Provide compound ID for all files HOT 5
- Parmed bug affecting cyclodextrin files HOT 6
- Update README.md to reflect how to cite (once discussion is final) HOT 4
- Mixed 1-4 scaling? HOT 2
- Typo in reference for Lysozyme L99A/M102Q - catechol experimental value? HOT 3
- Errors using cyclic host `mol2` files with non-unique atom names HOT 11
- Add machine-readable tables for all the sets HOT 3
- The SDF files for charged compounds aren't listing the formal charges of the atoms HOT 2
- PDB files for BRD4 parameterized systems are missing CONECT records HOT 2
- Include additional data, references
- gdcc-set2 and CB experimental Delta H is missing HOT 4
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 benchmarksets.