Comments (5)
Hi, thank you for your use and report.
I checked the options for electrostatics.
But, the program correctly imported the changed Mg concentration value.
If the effect is not noticeable, there are three possible reasons.
- The insensitive structure to electrostatic interaction: highly dense connectivity in the structure could prevent the structural deformation from the electrostatic force.
- The low cutoff value: the default cutoff value was set to 2.5 nm, but if the interhelical distance (or distance between substructures) is larger than the cutoff value, the electrostatic interaction would not affect the deformation much.
- The accuracy of the present model: in the present model (based on the Debye-Huckel model), the effective charge was expressed in terms of Mg concentration and based on the data for Mg 20 mM. So, the prediction for other Mg concentrations could not accurate because we do not have much structural data for different concentrations.
If you need my help more, please send me the files. I will check it.
Thank you,
Jae Young Lee
from snupi.
Yes, you are correct I missed setting the cutoff value.
Apologies for erroneously reporting it as a bug.
I performed a small screen for Mg concentrations from 5 to 25 with a cutoff of 5.0 nm.
Comparison with the (2.5nm cutoff - 20Mg) run showed that the resulting structure more strongly resemble the (5nm cutoff - 5Mg) run.
Qualitative comparison with experimental cryo-EM also shows significantly better agreement with the 5nm runs, and correct Mg concentration.
A note on compute time:
I tested a simple, 1033bp scaffold structure. unsurprisingly doubling the cutoff significantly affect performance.
CPU TIME (ANALYSIS) jumped from 16 [s] to 32.17 [s].
which in my opinion it still incredibly fast :)
Based on this limited analysis I would propose a higher default value for electrostatic interaction.
what do you think?
from snupi.
Sorry for the late reply. I did not recognize the update.
The proper cutoff distance was introduced for efficient analysis, so increasing the cutoff distance could provide a correct result.
But, increasing the cutoff distance also brings much calculation of searching interacting nodes and generating electrostatic elements, making the system heavy.
This could especially be severe for a highly packed system.
The current value (2.5 nm) was selected in consideration of accuracy and efficiency after analyzing many structures, so changing the default value should be careful..
But of course, you can choose the value depending on your structure/environment.
Thank you,
Jae Young Lee
from snupi.
For most of our single structure use cases a factor 2-3 in runtime hardly makes much difference.
But of course, you are right that highly packed systems might require significantly longer computation.
- I now use 3.75nm electrostatics cutoff as my default.
- If faster computation is required I reduce the cutoff to the default value of 2.5nm
best, elija
from snupi.
Thank you for your trials and opinion. The effects of the electrostatic parameters including cutoff are of importance and interest. Your report is surely helpful for future improvement.
Thank you,
Jae Young
from snupi.
Related Issues (12)
- Design error HOT 2
- Import Sequences HOT 1
- SNUPI on Linux? HOT 1
- single-stranded loops HOT 3
- Mac - Protocol.pdf missing HOT 1
- Mac - SNUPI2 - gen_PDB error HOT 2
- Mac - SNUPI2.1 - RMSF HOT 2
- SNUPI2.1 - Stacking HOT 3
- MAC v3_00 Error: Cannot find CTF archive HOT 10
- Inconsistency between the design files and the oxDNA output. HOT 6
- SNUPI on M1 chip HOT 1
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 snupi.