Giter VIP home page Giter VIP logo

gmpavanlab / swarm-cg Goto Github PK

View Code? Open in Web Editor NEW
41.0 41.0 9.0 32.86 MB

Swarm-CG: Automatic Parametrization of Bonded Terms in MARTINI-based Coarse-Grained Models of Simple to Complex Molecules via Fuzzy Self-Tuning Particle Swarm Optimization

Home Page: https://pubs.acs.org/doi/10.1021/acsomega.0c05469

License: MIT License

Python 99.81% Shell 0.19%
bonded-parameters bonded-terms cg-model coarse-grained coarse-graining gromacs molecular-dynamics molecular-modeling optimization optimization-tools swarm-cg

swarm-cg's Issues

Swarm-CG & Gromacs 2021

Hi, I'd like to use your model to optimize bonded parameters for a drug molecule. Just wanted to check if your tool is compatible with the latest version of gromacs (2021) and does it currently handle virtual sites?

scg_optimize fails on WSL due to permission denied

Dear Swarm-CG creators,

Running the scg_optimize executable on WSL crashes with the following error:

shutil.Error: [('.internal/input_CG_simulation_files', 'CG_sim_files_eval_step_1', "[Errno 13] Permission denied: 
'CG_sim_files_eval_step_1'")]

We traced the error to line 2690 in swarmCG.py:

# create new directory for new parameters evaluation
current_eval_dir = f'{config.iteration_sim_files_dirname}_eval_step_{ns.nb_eval}'
shutil.copytree(config.input_sim_files_dirname, current_eval_dir)

We found that there is a bug in the copytree function on WSL (https://bugs.python.org/issue38633). A temporary fix could be to add the following lines at the beginning of the file:

import errno
orig_copyxattr = shutil._copyxattr
def patched_copyxattr(src, dst, *, follow_symlinks=True):
  try:
    orig_copyxattr(src, dst, follow_symlinks=follow_symlinks)
  except OSError as ex:
    if ex.errno != errno.EACCES: raise
shutil._copyxattr = patched_copyxattr

Thanks to @franciscoadasme for his help in the debugging.

I hope this helps!

Constraints optimization

Dear Swarm-CG creators,
First of all, great project!
I have issue with constraints optimization. The addition of constraints into the itp file always crushes the optimization with the following error:

-- ! ERROR ! -- In the provided CG ITP file constraints have been grouped, but constraints group 1 holds lines that have different parameters. Parameters should be identical within a group, only CG beads IDs should differ. Please correct the CG ITP file and separate groups using a blank or commented line.

But all parameters inside the group are identical in my itp file, moreover even addition only one bond constraint results in the same error. Without any constraint all works perfectly fine.

itp_file.zip

VirtualSites are silently dropped from itp

Hi,

First of all I want to say that this is a very interesting project. I'm very much looking forward to using the program in some high throughout parametrization of models. However, it seems that the program fails to include virtual sites. I assume they are dropped because they cannot be optimized. However, the program also does not issue a warning or error message at the pre process stage.

From a martini centric perspective it would be very useful to have the code simply keep virtual sites, because a lot of martini3 models are using virtual sites. Can this be implemented in the current code?

EDIT: Actually after hacking in VS treating them like exclusions, I figured out that this will not be possible to have because MDAnalysis make whole does not support making whole molecules that have VS (i.e. not contiguous bonds).

Extended Usage Request: Adding GROMACS angle function 10 - Restriced Bending Angle Potential

Hello again,

I'm using your program and it is awesome!

I do have one suggestion, that I think may help many CG users. That is adding the restricted bending angle potential (gromacs function 10) to the options for angular potential optimization. The restricted bending angle potential stabilizes systems and prevents them from blowing up in certain cases, especially when the CG model includes a dihedral angle. The reasoning behind the stabilization is illustrated in this exchange on the martini forum, where adding the restricted bending angle was needed to prevent my system from blowing up from lincs errors related co-linear atoms within the dihedral that caused the potential to diverge.

http://www.cgmartini.nl/index.php/component/kunena/15-water/5894-yet-another-thread-about-lincs-warnings-in-polarized-martini-water#8632

Cheers!
Mike

Newer gromacs version than 2019 have a different tpr and are not supported?

Hoi,

I run into an MDAnalysis error due to a non-supported tpr format. I understand the issue with the ever changing format. Therefore I was wondering if you have some ideas to work around this limitation?

I would really like to use your software, but 2019 is 5 years ago and I guess it would be cool if we are not locked into outdated software?

Love to hear your idea on this issue.

Cheers,

Bart

ValueError: Failed to construct topology from file 4-prod.tpr with parser <class 'MDAnalysis.topology.TPRParser.TPRParser'>.
Error: Your tpx version is 133, which this parser does not support, yet

md run failing

Starting iteration 4 at 14:14:40 on 12-10-2022
MD run failed (simulation process terminated with error)
Iteration time: 0.5 min

Even the distribution plots are not appropriate.
Kindly tell me how to resolve this issue.

AttributeError: 'Namespace' object has no attribute 'gmx_cmd'

Hello,

I'm running into a pesky AttributeError: 'Namespace' object has no attribute 'gmx_cmd' when running scg_optimize on a polymer I am trying to map using your program.

It seems everything is running smoothly until the minimization step, but then an error is triggered when it tries run the equilibration step of the first CG iteration. I attached my input and output files here. I didn't include the atomistic trajectory as the file is large, and I think this might be a bug on the cg-swarm side given that the minimization runs fine and this is a python based error. If you need the trajectory to reproduce the error, let me know and I can drop a google drive link here.

I use slurm to run the bash script 'cg-swarm-run' in the zip file below. I think if you look there first, the rest of the files should be arranged according to your example.

cg_swarm_testing(2).zip

Please let me know if I made a mistake anywhere!

Cheers,
Mike Boyle

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo 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.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.