vacuum minimization

Author: ndonyapour

cwl_adapters/combine_structure.cwl

@@ -12,7 +12,7 @@ baseCommand: ['python', '/combine_structure.py']
 
 hints:
   DockerRequirement:
-    dockerPull: ndonyapour/combine_structure
+    dockerPull: ndonyapour/combine_structure_zip
 
 inputs:
   input_structure1:
@@ -40,6 +40,20 @@ inputs:
     - edam:format_3877 
     inputBinding:
       prefix: --input_structure2
+  
+  input_top_zip_path:
+    label: Input zip file
+    doc: |-
+      Input zip file
+      Type: string
+      File type: input
+      Accepted formats: zip
+      Example file: https://github.com/bioexcel/biobb_md/blob/master/biobb_md/test/data/gromacs/genion.zip
+    type: File
+    format:
+    - edam:format_3987
+    inputBinding:
+      prefix: --input_top_zip_path
 
   output_structure_path:
     label: Output combined PDB file path

examples/docking/docking.yml

@@ -90,6 +90,7 @@ steps:
       in:
         input_structure1: '*receptor.xyz' # '*receptor_hydrogens.pdb'
         input_structure2: '*pose_ligand.xyz'
+        input_top_zip_path: '*complex_vac.zip'
         output_structure_path: '&complex_vac.pdb'
 
 wic:

examples/docking/dsb.yml

@@ -9,23 +9,17 @@ steps:
       sdf_path: ~sdf_path
 - stability.yml:
     in:
-      # NOTE: crd_path and top_zip_path are currently only necessary
-      # due to a limitation of the inlineing implementation.
-      # Moreover, crd_path and top_zip_path are overridden below
       crd_path: '*complex_vac.pdb'
       top_zip_path: '*complex_vac.zip'
-      nsteps: 100000
+      nsteps: 10000
       dt: 0.002
       temperature: 298.0
       pressure: 1.0
 - stability.yml:
     in:
-      # NOTE: crd_path and top_zip_path are currently only necessary
-      # due to a limitation of the inlineing implementation.
-      # Moreover, crd_path and top_zip_path are overridden below
-      crd_path: '*complex_vac.pdb'
-      top_zip_path: '*complex_vac.zip'
-      nsteps: 100000
+      crd_path: '*ligand_GMX.gro'
+      top_zip_path: '*ligand_GMX.zip'
+      nsteps: 10000
       dt: 0.002
       temperature: 298.0
       pressure: 1.0
@@ -41,13 +35,18 @@ wic:
         environment:
           action: checkpoint
         steps:
-          (1, setup.yml):
+          (1, setup_vac_min.yml):
             wic:
               steps:
-                (4, genion):
-                  in:
-                    output_top_zip_path: '&genion_complex.zip'
+                (3, solv_ion.yml):
+                  wic:
+                    steps:
+                      (4, genion):
+                        in:
+                          output_top_zip_path: '&genion_complex.zip'
           (2, basic.yml):
+            in:
+              top_zip_path: '*genion_complex.zip'
             wic:
               steps:
                 (3, prod.yml):
@@ -82,13 +81,24 @@ wic:
           action: restore
           save_defs: ['genion_complex.zip', 'prod_complex.gro']
         steps:
-          (1, setup.yml):
+          (1, setup_vac_min.yml):
             wic:
               steps:
-                (4, genion):
-                  in:
-                    output_top_zip_path: '&genion_ligand.zip'
+                (1, topology.yml):
+                  wic:
+                    steps:
+                      (1, editconf):
+                        in: 
+                          output_crd_path: '&ligand_box.g96'
+                (3, solv_ion.yml):
+                  wic:
+                    steps:
+                      (4, genion):
+                        in:
+                          output_top_zip_path: '&genion_ligand.zip'
           (2, basic.yml):
+            in:
+              top_zip_path: '*genion_ligand.zip'
             wic:
               steps:
                 (3, prod.yml):

examples/docking/gen_topol_params.yml

@@ -70,6 +70,7 @@ steps:
     in:
       input_structure1: ~input_receptor_xyz_path
       input_structure2: '*pose_ligand.xyz'
+      input_top_zip_path: '*complex_vac.zip'
       output_structure_path: '&complex_vac.pdb'
 
 wic:

examples/docking/vs_demo_2.yml

@@ -146,13 +146,18 @@ wic:
       wic:
         inlineable: False
         steps:
-          (1, setup.yml):
+          (1, setup_vac_min.yml):
             wic:
               steps:
-                (4, genion):
-                  in:
-                    output_top_zip_path: '&genion_complex.zip'
+                (3, solv_ion.yml):
+                  wic:
+                    steps:
+                      (4, genion):
+                        in:
+                          output_top_zip_path: '&genion_complex.zip'
           (2, basic.yml):
+            in:
+              top_zip_path: '*genion_complex.zip'
             wic:
               steps:
                 (3, prod.yml):

examples/docking/vs_demo_3.yml

@@ -124,13 +124,18 @@ wic:
       wic:
         inlineable: False
         steps:
-          (1, setup.yml):
+          (1, setup_vac_min.yml):
             wic:
               steps:
-                (4, genion):
-                  in:
-                    output_top_zip_path: '&genion_complex.zip'
+                (3, solv_ion.yml):
+                  wic:
+                    steps:
+                      (4, genion):
+                        in:
+                          output_top_zip_path: '&genion_complex.zip'
           (2, basic.yml):
+            in:
+              top_zip_path: '*genion_complex.zip'
             wic:
               steps:
                 (3, prod.yml):

examples/docking/vs_demo_4.yml

@@ -191,13 +191,18 @@ wic:
       wic:
         inlineable: False
         steps:
-          (1, setup.yml):
+          (1, setup_vac_min.yml):
             wic:
               steps:
-                (4, genion):
-                  in:
-                    output_top_zip_path: '&genion_complex.zip'
+                (3, solv_ion.yml):
+                  wic:
+                    steps:
+                      (4, genion):
+                        in:
+                          output_top_zip_path: '&genion_complex.zip'
           (2, basic.yml):
+            in:
+              top_zip_path: '*genion_complex.zip'
             wic:
               steps:
                 (3, prod.yml):

examples/gromacs/basic.yml

@@ -1,4 +1,8 @@
 inputs:
+  top_zip_path:
+    type: File
+    format:
+    - edam:format_3987
   nsteps:
     type: int
   dt:
@@ -16,6 +20,8 @@ steps:
 # to this level and use the following syntax to pass it to only one call site.
 # (Moreover, the former will create multiple definitions of &min.tpr, which is not allowed.)
   - min.yml:
+      in:
+        top_zip_path: ~top_zip_path
   - equil.yml:
   - prod.yml:
       in:

examples/gromacs/cg.yml

@@ -1,6 +1,13 @@
+inputs:
+  top_zip_path:
+    type: File
+    format:
+    - edam:format_3987
+
 steps:
   - grompp:
       in:
+        input_top_zip_path: ~top_zip_path
         config:
           mdp:
             integrator: cg

examples/gromacs/min.yml

@@ -1,9 +1,19 @@
+inputs:
+  top_zip_path:
+    type: File
+    format:
+    - edam:format_3987
+
 steps:
   - steep.yml:
+      in:
+        top_zip_path: ~top_zip_path
 # Fatal error:
 # The coordinates could not be constrained. Minimizer 'cg' can not handle
 # constraint failures, use minimizer 'steep' before using 'cg'.
   - cg.yml:
+      in:
+        top_zip_path: ~top_zip_path
 # Fatal error:
 # The combination of constraints and L-BFGS minimization is not implemented.
 # Either do not use constraints, or use another minimizer (e.g. steepest