diff --git a/_build/.doctrees/cli_tutorials.doctree b/_build/.doctrees/cli_tutorials.doctree new file mode 100644 index 0000000..dd1e94c Binary files /dev/null and b/_build/.doctrees/cli_tutorials.doctree differ diff --git a/_build/.doctrees/configuration.doctree b/_build/.doctrees/configuration.doctree new file mode 100644 index 0000000..ad3bf38 Binary files /dev/null and b/_build/.doctrees/configuration.doctree differ diff --git a/_build/.doctrees/environment.pickle b/_build/.doctrees/environment.pickle index 36c8459..cd7d1e1 100644 Binary files a/_build/.doctrees/environment.pickle and b/_build/.doctrees/environment.pickle differ diff --git a/_build/.doctrees/index.doctree b/_build/.doctrees/index.doctree index 15b3e1f..883eff4 100644 Binary files a/_build/.doctrees/index.doctree and b/_build/.doctrees/index.doctree differ diff --git a/_build/.doctrees/installation.doctree b/_build/.doctrees/installation.doctree new file mode 100644 index 0000000..b13b7c8 Binary files /dev/null and b/_build/.doctrees/installation.doctree differ diff --git a/_build/.doctrees/python_tutorials.doctree b/_build/.doctrees/python_tutorials.doctree new file mode 100644 index 0000000..8f3127f Binary files /dev/null and b/_build/.doctrees/python_tutorials.doctree differ diff --git a/_build/.doctrees/tutorial_cli_abfe.doctree b/_build/.doctrees/tutorial_cli_abfe.doctree new file mode 100644 index 0000000..eef2502 Binary files /dev/null and b/_build/.doctrees/tutorial_cli_abfe.doctree differ diff --git a/_build/.doctrees/tutorial_cli_production_md.doctree b/_build/.doctrees/tutorial_cli_production_md.doctree new file mode 100644 index 0000000..3d5cdc3 Binary files /dev/null and b/_build/.doctrees/tutorial_cli_production_md.doctree differ diff --git a/_build/.doctrees/tutorial_cli_system_preparation.doctree b/_build/.doctrees/tutorial_cli_system_preparation.doctree new file mode 100644 index 0000000..df784c5 Binary files /dev/null and b/_build/.doctrees/tutorial_cli_system_preparation.doctree differ diff --git a/_build/.doctrees/tutorial_python_production.doctree b/_build/.doctrees/tutorial_python_production.doctree new file mode 100644 index 0000000..4a01c51 Binary files /dev/null and b/_build/.doctrees/tutorial_python_production.doctree differ diff --git a/_build/_sources/cli_tutorials.rst.txt b/_build/_sources/cli_tutorials.rst.txt new file mode 100644 index 0000000..f0c8539 --- /dev/null +++ b/_build/_sources/cli_tutorials.rst.txt @@ -0,0 +1,22 @@ +CLI Tutorials +============== + +maize-biosimspace provides a large number of nodes for performing common operations such as parameterisation, minimisation, equilibration, production molecular +dynamics, and alchemical free energy calculations. Nodes are generally engine-specific as maize requires that nodes have a list of ``required_callables`` whose presence +in the environment is checked before the node is run. To see a list of all available nodes, type ``bss_`` and hit tab to list them. To see the options for each, pass +pass the ``-h`` flag e.g. + +.. code-block:: bash + + bss_parameterise -h + +These tutorials give specific examples of using BioSimSpace maize nodes to run production molecular dynamics, to create equilibrated systems starting from unparameterised +input structures, and to run absolute binding free energy calculations starting from a protein pdb and an sdf file containing multiple ligands. + +.. toctree:: + :maxdepth: 1 + + tutorial_cli_production_md + tutorial_cli_system_preparation + tutorial_cli_abfe + \ No newline at end of file diff --git a/_build/_sources/configuration.rst.txt b/_build/_sources/configuration.rst.txt new file mode 100644 index 0000000..e0277a5 --- /dev/null +++ b/_build/_sources/configuration.rst.txt @@ -0,0 +1,22 @@ +Configuration +------------- +For a detailed guide to configuring Maize, please see the `maize documentation `_ +and the `maize-contrib documentation `_. + +By default, Maize looks for `$XDG_CONFIG_HOME/maize.toml` (`~/.config/maize.toml`) for configuration information. To set up Maize to run through slurm, run + +.. code-block:: bash + + export XDG_CONFIG_HOME=~/.config + +Then create `~/.config/maize.toml` containing the following + +.. code-block:: toml + + system = "slurm" # Can be one of {'cobalt', 'flux', 'local', 'lsf', 'pbspro', 'rp', 'slurm'} + max_jobs = 100 # The maximum number of jobs that can be submitted by a node at once + queue = "gpu" #CHANGEME to your desired GPU queue + launcher = "srun" # The launcher to use for the command, usually one of {'srun', 'mpirun', 'mpiexec'} + walltime = "24:00:00" # Job walltime limit, shorter times may improve queueing times + +Now when you execute a Maize workflow, it will submit jobs to the slurm scheduler. \ No newline at end of file diff --git a/_build/_sources/index.rst.txt b/_build/_sources/index.rst.txt index 585dc11..0e0acee 100644 --- a/_build/_sources/index.rst.txt +++ b/_build/_sources/index.rst.txt @@ -10,9 +10,18 @@ maize-biosimspace .. toctree:: :hidden: :maxdepth: 1 - :caption: Examples + :caption: Setup - production-md + installation + configuration + +.. toctree:: + :hidden: + :maxdepth: 2 + :caption: Tutorials + + cli_tutorials + python_tutorials .. toctree:: :hidden: @@ -29,37 +38,37 @@ maize-biosimspace Steps Maize -Installation ------------- -To install, simply clone this repository and run: +Quick Start +----------- +Clone the repository and install: .. code-block:: bash + git clone https://github.com/fjclark/maize-biosimspace.git + cd maize-biosimspace mamba env create -f env-users.yml mamba activate maize-biosimspace pip install --no-deps ./ -If you want to keep up-to-date with the latest changes to the core, clone `maize `_, switch to the directory, and run (in the same conda environment): +Ensure that the required AMBER, GROMACS, and NAMD executables are available in your environment if you plan to use the relevant nodes. -.. code-block:: bash +Many maize-biosimspace workflows will now be available through CLIs. Try typing `bss_` and hitting tab to list them. To see the options for each, type e.g. - pip install --no-deps ./ - -If you plan on developing, you should use ``env-dev.yml`` instead and use the ``-e`` flag for ``pip``. +.. code-block:: bash -Configuration -------------- -Each step documentation will contain information on how to setup and run the node, as well as install the required dependencies. Dependencies can be managed in several ways, depending on the node and workflow you are running: + bss_parameterise -h -* Through a ``module`` system: +To run production 0.1 ns of production MD on a protein-ligand complex using GROMACS, dumping all of the intermediate files to the current directory: - Specify a module providing an executable in the ``config.toml`` (see `Configuring workflows `_). This module will then be loaded in the process running the node. +.. code-block:: bash -* With a separate python environment: + mkdir gmx_md_example + cd gmx_md_example + cp ../tests/data/complex.* . + bss_production_gromacs --inp complex.prm7 complex.rst7 --runtime 0.1 --save_name gmx_md_out --dump_to . - Some nodes will require custom python environments that are likely to be incompatible with the other environments. In those cases, the node process can be spawned in a custom environment. Note that this environment must still contain *maize*. +This will run locally. To configure maize nodes to run through a scheduler such as slurm, see :doc:`configuration`. -* By specifying the executable location and possibly script interpreter. This can also be accomplished using ``config.toml`` (see `Configuring workflows `_). Indices and tables diff --git a/_build/_sources/installation.rst.txt b/_build/_sources/installation.rst.txt new file mode 100644 index 0000000..b46d2e0 --- /dev/null +++ b/_build/_sources/installation.rst.txt @@ -0,0 +1,29 @@ +Installation +------------ +For a basic installation, simply clone this repository and run: + +.. code-block:: bash + + git clone https://github.com/fjclark/maize-biosimspace.git + cd maize-biosimspace + mamba env create -f env-users.yml + mamba activate maize-biosimspace + pip install --no-deps ./ + +Ensure that the required AMBER, GROMACS, and NAMD executables are available in your environment if you plan to use the relevant nodes. + +If you want to keep up-to-date with the latest changes to the core, clone `maize `_, switch to the directory, and run (in the same conda environment): + +.. code-block:: bash + + pip install --no-deps ./ + +If you plan on developing, you should use ``env-dev.yml`` instead and use the ``-e`` flag for ``pip``. + +.. code-block:: bash + + git clone https://github.com/fjclark/maize-biosimspace.git + cd maize-biosimspace + mamba env create -f env-dev.yml + mamba activate maize-biosimspace + pip install -e --no-deps ./ \ No newline at end of file diff --git a/_build/_sources/python_tutorials.rst.txt b/_build/_sources/python_tutorials.rst.txt new file mode 100644 index 0000000..cf2d4c0 --- /dev/null +++ b/_build/_sources/python_tutorials.rst.txt @@ -0,0 +1,13 @@ +Python Tutorials +================ + +Nodes and subgraphs from maize-biosimspace can be combined into more complex +workflows through python, and can benefit from other nodes available through +maize-contrib. You can also use pre-made maize-biosimspace workflows through python +to gain fuller control over all the options (compared to the CLI) and to write +more reusable scripts. + +.. toctree:: + :maxdepth: 1 + + tutorial_python_production \ No newline at end of file diff --git a/_build/_sources/tutorial_cli_abfe.rst.txt b/_build/_sources/tutorial_cli_abfe.rst.txt new file mode 100644 index 0000000..a084df1 --- /dev/null +++ b/_build/_sources/tutorial_cli_abfe.rst.txt @@ -0,0 +1,55 @@ +Absolute Binding Free Energy Calculations +========================================= + +Here, we'll run a quick absolute binding free energy calculation for benzene bound to T4 Lysozyme. For this, we'll +use the ``bss_abfe_multi_isomer`` workflow through its CLI, which requires only an SDF file containing all required +ligands, and the pdb of the protein. Check the options with + +.. code-block:: bash + + bss_abfe_multi_isomer -h + +Copy over the required input: + +.. code-block:: bash + + mkdir bss_abfe_example + cd bss_abfe_example + cp ../tests/data/benzene.sdf . + cp ../tests/data/t4l.pdb . + +Now, let's run a relatively very short (but still fairly expensive) +ABFE calculation with 2 replicates: + +.. code-block:: bash + + bss_abfe_multi_isomer --lig_sdfs_file benzene.sdf \ + --protein_pdb benzene.t4l \ + --ligand_force_field gaff2 \ + --protein_force_field ff14SB \ + --abfe_timestep 4 \ + --abfe_n_replicates 2 \ + --abfe_runtime 0.1 \ + --abfe_runtime_generate_boresch_restraint 0.1 \ + --prep_runtime_restrained_npt 0.05 \ + --prep_runtime_unrestrained_npt 0.05 \ + --abfe_estimator TI \ + --results_file_name abfe_out \ + +The ``abfe_out`` file should show results around -4 kcal / mol. + +Running through the command line with this many arguments is unweildy, +and some options aren't available through the CLI (for example, the lambda +spacing). It's likely a better option to write a quick script - using the +pre-made workflow directly in python - simply import the workflow factory +, customise the options, and run (all in a python script): + +.. code-block:: python + + from maize.graphs.exs.biosimspace.afe import getabfe_multi_isomer_workflow + workflow = getabfe_multi_isomer_workflow() + + # Set workflow options... + + # Run + workflow.execute() diff --git a/_build/_sources/tutorial_cli_production_md.rst.txt b/_build/_sources/tutorial_cli_production_md.rst.txt new file mode 100644 index 0000000..c2e8d5d --- /dev/null +++ b/_build/_sources/tutorial_cli_production_md.rst.txt @@ -0,0 +1,30 @@ +Production Molecular Dynamics +============================= + +Here, we'll run production molecular dynamics on a protein-ligand complex. To check the available production CLIs, type ``bss_production`` and hit tab: + +.. code-block:: bash + + bss_production_gromacs bss_production_pmemd bss_production_pmemd_cuda bss_production_sander bss_production_somd + +We'll pick gromacs. To check the available arguments and defaults, run + +.. code-block:: bash + + bss_production_gromacs -h + +We'll run a quick 0.1 ns of production molecular dynamics on the protein-ligand complex included with `maize-biosimspace` for testing. We'll specify the output +name to be ``gmx_md_out`` and we'll save all of the intermediate files (including input scripts, logs, and trajectory files) to a subdirectory in the current +working directory by specifying ``--dump_to .`` + +.. code-block:: bash + + mkdir gmx_md_example + cd gmx_md_example + cp ../tests/data/complex.* . + bss_production_gromacs --inp complex.prm7 complex.rst7 --runtime 0.1 --save_name gmx_md_out --dump_to . + +You should now have the final coordinate file, ``gmx_md_out.rst7``, a copy of the input topology file ``gmx_md_out.prm7``, and a sub-directory containing all of the +intermediate files. Note that despite running through GROMACS, we were able to pass in AMBER files as input. This is because BioSimSpace automatically converts +between file formats (using Sire under the hood). + diff --git a/_build/_sources/tutorial_cli_system_preparation.rst.txt b/_build/_sources/tutorial_cli_system_preparation.rst.txt new file mode 100644 index 0000000..b9879d6 --- /dev/null +++ b/_build/_sources/tutorial_cli_system_preparation.rst.txt @@ -0,0 +1,53 @@ +System Preparation +================== + +Often, we have structure files for a protein and / or ligand and we would like to parameterise, solvate, minimise, heat, and equilibrate them to obtain systems +suitable for production molecular dynamics simulations or free energy calculations. While individual nodes are provided for all of these steps, +`maize-biosimspace` also provides CLIs for two complete system preparation workflows: ``bss_system_prep_free``, which is designed to prepare a ligand in a box +of water (but can also be used for an apo protein) and ``bss_system_prep_bound``, which is designed to set up protein-ligand complexes. + +Here, we'll prepare the complex of T4 lysozyme L99A bound to benzene, a very common test system for absolute binding free energy calculations. First, we'll copy +over the required input: + +.. code-block:: bash + + mkdir sysprep_bound_example + cd sysprep_bound_example + cp ../tests/data/benzene.sdf . + cp ../tests/data/t4l.pdb . + +.. tip:: + + This pdb has been sanitised and will work with BioSimSpace first time, but often pdbs will require some tweaking before they are accepted by ``tleap`` (which + BioSimSpace uses behind the scenes). The recommended workflow is + + * Clean your unsanitised pdb using pdb4amber, e.g. ``pdb4amber -i protein.pdb -o protein_sanitised.pdb`` + * Attempt to run the workflow below + * If the workflow raises an error, attempt parameterisation directly with ``tleap`` to get more detailed error messages. E.g., type ``tleap``, then + + .. code-block:: bash + + source leaprc.protein.ff14SB + source leaprc.water.tip3p + # Loading an unsanitised pdb will likely raise an error + prot = loadpdb protein_sanitised.pdb + saveamberparm prot protein.parm7 protein.rst7 + savepdb prot protein_fully_sanitised.pdb + + * If the above fails, this is often due to residue/ atom names which do not match the templates. Read the errors to find out which residues / atoms are causing the issues, then check the expected names in library which was loaded after typing ``source leaprc.protein.ff14SB`` e.g. ``cat $AMBERHOME/dat/leap/lib/amino12.lib``. Rename the offending atoms/ residues and repeat the above step. + + BioSimSpace is very fussy about parameterisation and will fail if tleap raises any warnings. To get round this, run the tleap script above and use the + ``protein_full_sanitised.pdb`` file as your input, which will not raise any errors. + +To run system preparation for our protein-ligand complex, we'll use the ``bss_system_prep_bound`` CLI, saving the output system to "t4l_benzene_complex_equilibrated" +and using the gaff2 ff14SB force fields. There are a large number of other parameters which can be modified (see ``bss_system_prep_bound -``) but we'll run the +defaults for now. + +.. note:: + + Make sure that you have access to a GPU locally, or have configured `Maize` to submit to a queue with gpu access (see :doc:`configuration`) + +.. code-block:: bash + + bss_system_prep_bound --inp benzene.sdf --protein_pdb t4l.pdb --ligand_force_field gaff2 --protein_force_field ff14SB --save_name t4l_benezene_complex_equilibrated + diff --git a/_build/_sources/tutorial_python_production.rst.txt b/_build/_sources/tutorial_python_production.rst.txt new file mode 100644 index 0000000..46e1f42 --- /dev/null +++ b/_build/_sources/tutorial_python_production.rst.txt @@ -0,0 +1,63 @@ +Creating a Workflow with Production MD +======================================= + +Here, we'll look at a basic example of a custom workflow which +uses ``pmemd.cuda`` to run some production MD. In reality, you +would want to add some extra steps before or afterwards: + +.. code-block:: python + + """Run production Molecular Dynamics using PMEMD.CUDA through BioSimSpace.""" + + from pathlib import Path + + from maize.core.workflow import Workflow + from maize.steps.exs.biosimspace import ProductionPmemdCuda + from maize.steps.io import LoadData, Return + from maize.utilities.execution import JobResourceConfig + + # Build the graph + flow = Workflow(name="Prod_BSS_AMBER_Test", cleanup_temp=False, level="debug") + + # Add the nodes + load_sys = flow.add(LoadData[list[Path]]) + prod_pmemd = flow.add( + ProductionPmemdCuda, + name="Production_Amber", + parameters={ + "runtime": 1.0, # ns + }, + ) + retu = flow.add(Return[list[Path]]) + + # Set parameters + load_sys.data.set( + [ + Path( + "< path to complex.prm7>" # CHANGEME + ), + Path( + "< path to complex.rst7>" # CHANGEME + ), + ] + ) + + # Connect the nodes + flow.connect(load_sys.out, prod_pmemd.inp) + flow.connect(prod_pmemd.out, retu.inp) + + # Check and run! + flow.check() + flow.visualize() + flow.execute() + + mols = retu.get() + + # Load a BioSimSpace system from the returned paths + import BioSimSpace as BSS + + sys = BSS.IO.readMolecules([str(mols[0]), str(mols[1])]) + print(40 * "#") + print(sys) + # In reality, you would do something here... + diff --git a/_build/cli_tutorials.html b/_build/cli_tutorials.html new file mode 100644 index 0000000..bf0eda8 --- /dev/null +++ b/_build/cli_tutorials.html @@ -0,0 +1,301 @@ + + + + + + + + + CLI Tutorials - maize biosimspace + + + + + + + + + + + + + + + + Contents + + + + + + Menu + + + + + + + + Expand + + + + + + Light mode + + + + + + + + + + + + + + Dark mode + + + + + + + Auto light/dark mode + + + + + + + + + + + + + + + + + + + +
+
+
+ +
+
+
maize biosimspace
 +
+
+
+ +
+ +
+
+ +
+
+
+ + + + + Back to top + +
+ +
+ +
+ +
+
+
+

CLI Tutorials#

+

maize-biosimspace provides a large number of nodes for performing common operations such as parameterisation, minimisation, equilibration, production molecular +dynamics, and alchemical free energy calculations. Nodes are generally engine-specific as maize requires that nodes have a list of required_callables whose presence +in the environment is checked before the node is run. To see a list of all available nodes, type bss_ and hit tab to list them. To see the options for each, pass +pass the -h flag e.g.

+
bss_parameterise -h
+
+
+

These tutorials give specific examples of using BioSimSpace maize nodes to run production molecular dynamics, to create equilibrated systems starting from unparameterised +input structures, and to run absolute binding free energy calculations starting from a protein pdb and an sdf file containing multiple ligands.

+
+ +
+
+ +
+
+ +
+ +
+
+ + + + + + \ No newline at end of file diff --git a/_build/configuration.html b/_build/configuration.html new file mode 100644 index 0000000..06c7ce3 --- /dev/null +++ b/_build/configuration.html @@ -0,0 +1,300 @@ + + + + + + + + + Configuration - maize biosimspace + + + + + + + + + + + + + + + + Contents + + + + + + Menu + + + + + + + + Expand + + + + + + Light mode + + + + + + + + + + + + + + Dark mode + + + + + + + Auto light/dark mode + + + + + + + + + + + + + + + + + + + +
+
+
+ +
+
+
maize biosimspace
 +
+
+
+ +
+ +
+
+ +
+
+
+ + + + + Back to top + +
+ +
+ +
+ +
+
+
+

Configuration#

+

For a detailed guide to configuring Maize, please see the maize documentation +and the maize-contrib documentation.

+

By default, Maize looks for $XDG_CONFIG_HOME/maize.toml (~/.config/maize.toml) for configuration information. To set up Maize to run through slurm, run

+
export XDG_CONFIG_HOME=~/.config
+
+
+

Then create ~/.config/maize.toml containing the following

+
system = "slurm"  # Can be one of {'cobalt', 'flux', 'local', 'lsf', 'pbspro', 'rp', 'slurm'}
+max_jobs = 100  # The maximum number of jobs that can be submitted by a node at once
+queue = "gpu" #CHANGEME to your desired GPU queue
+launcher = "srun"  # The launcher to use for the command, usually one of {'srun', 'mpirun', 'mpiexec'}
+walltime = "24:00:00"  # Job walltime limit, shorter times may improve queueing times
+
+
+

Now when you execute a Maize workflow, it will submit jobs to the slurm scheduler.

+
+ +
+
+ +
+ +
+
+ + + + + + \ No newline at end of file diff --git a/_build/genindex.html b/_build/genindex.html index 8708618..466f125 100644 --- a/_build/genindex.html +++ b/_build/genindex.html @@ -158,13 +158,27 @@
-

Examples

+

Setup

+

Tutorials

+

Reference

    -
  • Steps
      +
    • Steps
      • diff --git a/_build/index.html b/_build/index.html index 9d1561d..0d4d14c 100644 --- a/_build/index.html +++ b/_build/index.html @@ -3,7 +3,7 @@ - + maize biosimspace @@ -160,13 +160,27 @@
      -

      Examples

      +

      Setup

      +

      Tutorials

      +

      Reference

        -
      • Steps
          +
        • Steps
          • @@ -218,32 +232,31 @@

          maize-biosimspace

      -
      -

      Installation#

      -

      To install, simply clone this repository and run:

      -
      mamba env create -f env-users.yml
+
      
+
      
+
      
+
      Quick Start#
      
+
      Clone the repository and install:
      
+
      git clone https://github.com/fjclark/maize-biosimspace.git
+cd maize-biosimspace
+mamba env create -f env-users.yml
 mamba activate maize-biosimspace
 pip install --no-deps ./
 
      
 
      
-
      If you want to keep up-to-date with the latest changes to the core, clone maize, switch to the directory, and run (in the same conda environment):
      
-
      pip install --no-deps ./
+
      Ensure that the required AMBER, GROMACS, and NAMD executables are available in your environment if you plan to use the relevant nodes.
      
+
      Many maize-biosimspace workflows will now be available through CLIs. Try typing bss_ and hitting tab to list them. To see the options for each, type e.g.
      
+
      bss_parameterise -h
 
      
 
      
-
      If you plan on developing, you should use env-dev.yml instead and use the -e flag for pip.
      
-
      
-
      
-
      Configuration#
      
-
      Each step documentation will contain information on how to setup and run the node, as well as install the required dependencies. Dependencies can be managed in several ways, depending on the node and workflow you are running:
      
-
        
-
      • Through a module system:
        
-
        Specify a module providing an executable in the config.toml (see Configuring workflows). This module will then be loaded in the process running the node.
        
-
        • 
-
      • With a separate python environment:
        
-
        Some nodes will require custom python environments that are likely to be incompatible with the other environments. In those cases, the node process can be spawned in a custom environment. Note that this environment must still contain maize.
        
-
        • 
-
      • By specifying the executable location and possibly script interpreter. This can also be accomplished using config.toml (see Configuring workflows).
        • 
-
      
+
      To run production 0.1 ns of production MD on a protein-ligand complex using GROMACS, dumping all of the intermediate files to the current directory:
      
+
      mkdir gmx_md_example
+cd gmx_md_example
+cp ../tests/data/complex.* .
+bss_production_gromacs --inp complex.prm7 complex.rst7 --runtime 0.1 --save_name gmx_md_out --dump_to .
+
      
+
      
+
      This will run locally. To configure maize nodes to run through a scheduler such as slurm, see Configuration.
      
 
      
 
      
 
      Indices and tables#
      
@@ -260,12 +273,12 @@ 
      Indices and tables
-          
+          
               
      
                 
      
                   Next
                 
      
-                
      Production MD
      
+                
      Installation
      
               
      
               
                   
@@ -301,8 +314,7 @@ 
      Indices and tables
             
        
 
      • maize-biosimspace
 
        • 
diff --git a/_build/installation.html b/_build/installation.html
new file mode 100644
index 0000000..60cc9c6
--- /dev/null
+++ b/_build/installation.html
@@ -0,0 +1,306 @@
+
+
+  
+    
+    
+
+
+    
+        Installation - maize biosimspace
+      
+    
+    
+    
+    
+
+
+
+  
+    
+    
+    
+
+
+  
+    Contents
+    
+      
+    
+  
+  
+    Menu
+    
+      
+      
+      
+    
+  
+  
+    Expand
+    
+      
+    
+  
+  
+    Light mode
+    
+      
+      
+      
+      
+      
+      
+      
+      
+      
+    
+  
+  
+    Dark mode
+    
+      
+      
+    
+  
+  
+    Auto light/dark mode
+    
+      
+      
+      
+      
+      
+      
+      
+    
+  
+
+
+
+
+
+
+
+
+
+
        
+  
        
+    
        
+      
+    
        
+    
        
+      
        maize biosimspace
        
+    
        
+    
        
+      
        
+        
+      
        
+      
+    
        
+  
        
+  
+  
        
+    
        
+      
        
+        
+          
+            
+          
+          Back to top
+        
+        
        
+          
+
        
+            
+          
        
+          
+        
        
+        
        
+          
        
+
        Installation#
        
+
        For a basic installation, simply clone this repository and run:
        
+
        git clone https://github.com/fjclark/maize-biosimspace.git
+cd maize-biosimspace
+mamba env create -f env-users.yml
+mamba activate maize-biosimspace
+pip install --no-deps ./
+
        
+
        
+
        Ensure that the required AMBER, GROMACS, and NAMD executables are available in your environment if you plan to use the relevant nodes.
        
+
        If you want to keep up-to-date with the latest changes to the core, clone maize, switch to the directory, and run (in the same conda environment):
        
+
        pip install --no-deps ./
+
        
+
        
+
        If you plan on developing, you should use env-dev.yml instead and use the -e flag for pip.
        
+
        git clone https://github.com/fjclark/maize-biosimspace.git
+cd maize-biosimspace
+mamba env create -f env-dev.yml
+mamba activate maize-biosimspace
+pip install -e --no-deps ./
+
        
+
        
+
        
+
+        
        
+      
        
+      
+    
        
+    
+  
        
+
        
+    
+    
+    
+    
+    
+
\ No newline at end of file
diff --git a/_build/objects.inv b/_build/objects.inv
index 68de4f7..404704e 100644
Binary files a/_build/objects.inv and b/_build/objects.inv differ
diff --git a/_build/production-md.html b/_build/production-md.html
index fae6836..5238174 100644
--- a/_build/production-md.html
+++ b/_build/production-md.html
@@ -3,7 +3,7 @@
   
     
     
-
+
 
     
         Production MD - maize biosimspace
@@ -160,13 +160,27 @@
   
 
 
        
-  
        Examples
        
-
        
         
        
           
        
diff --git a/_build/python_tutorials.html b/_build/python_tutorials.html
new file mode 100644
index 0000000..2c56391
--- /dev/null
+++ b/_build/python_tutorials.html
@@ -0,0 +1,295 @@
+
+
+  
+    
+    
+
+
+    
+        Python Tutorials - maize biosimspace
+      
+    
+    
+    
+    
+
+
+
+  
+    
+    
+    
+
+
+  
+    Contents
+    
+      
+    
+  
+  
+    Menu
+    
+      
+      
+      
+    
+  
+  
+    Expand
+    
+      
+    
+  
+  
+    Light mode
+    
+      
+      
+      
+      
+      
+      
+      
+      
+      
+    
+  
+  
+    Dark mode
+    
+      
+      
+    
+  
+  
+    Auto light/dark mode
+    
+      
+      
+      
+      
+      
+      
+      
+    
+  
+
+
+
+
+
+
+
+
+
+
        
+  
        
+    
        
+      
+    
        
+    
        
+      
        maize biosimspace
        
+    
        
+    
        
+      
        
+        
+      
        
+      
+    
        
+  
        
+  
+  
        
+    
        
+      
        
+        
+          
+            
+          
+          Back to top
+        
+        
        
+          
+
        
+            
+          
        
+          
+        
        
+        
        
+          
        
+
        Python Tutorials#
        
+
        Nodes and subgraphs from maize-biosimspace can be combined into more complex
+workflows through python, and can benefit from other nodes available through
+maize-contrib. You can also use pre-made maize-biosimspace workflows through python
+to gain fuller control over all the options (compared to the CLI) and to write
+more reusable scripts.
        
+
        
+
+
        
+
        
+
+        
        
+      
        
+      
+    
        
+    
+  
        
+
        
+    
+    
+    
+    
+    
+
\ No newline at end of file
diff --git a/_build/search.html b/_build/search.html
index 4264d49..1067658 100644
--- a/_build/search.html
+++ b/_build/search.html
@@ -157,13 +157,27 @@
   
 
 
        
-  
        Examples
        
+  
        Setup
        
 
+
        Tutorials
        
+
 
        Reference
        
 
          
-
        • Steps
            
+
          • Steps
 
            • 
diff --git a/_build/searchindex.js b/_build/searchindex.js
index 89ea088..6b5ab02 100644
--- a/_build/searchindex.js
+++ b/_build/searchindex.js
@@ -1 +1 @@
-Search.setIndex({"docnames": ["index", "production-md", "steps/biosimspace", "steps/index"], "filenames": ["index.rst", "production-md.rst", "steps/biosimspace.rst", "steps/index.rst"], "titles": ["maize-biosimspace", "Production MD", "BioSimSpace", "Steps"], "terms": {"i": [0, 2, 3], "graph": [0, 1], "base": 0, "workflow": [0, 1, 2], "manag": 0, "comput": 0, "chemistri": 0, "pipelin": 0, "thi": [0, 2, 3], "repositori": 0, "contain": [0, 3], "namespac": 0, "packag": [0, 2], "provid": [0, 2], "extens": 0, "you": [0, 1], "can": 0, "find": 0, "core": [0, 1], "document": [0, 3], "here": [0, 1], "To": 0, "simpli": 0, "clone": 0, "run": [0, 1], "mamba": 0, "env": 0, "creat": 0, "f": 0, "user": 0, "yml": 0, "activ": 0, "pip": 0, "dep": 0, "If": 0, "want": 0, "keep": 0, "up": 0, "date": 0, "latest": 0, "chang": 0, "switch": 0, "directori": 0, "same": 0, "conda": 0, "environ": 0, "plan": 0, "develop": 0, "should": 0, "us": [0, 1, 2], "dev": 0, "instead": 0, "e": 0, "flag": 0, "each": 0, "step": [0, 1], "inform": 0, "how": 0, "setup": 0, "node": [0, 1, 2], "well": 0, "requir": 0, "depend": 0, "sever": 0, "wai": 0, "ar": [0, 2], "through": [0, 1], "modul": 0, "system": [0, 1], "specifi": 0, "an": [0, 1], "execut": [0, 1], "config": 0, "toml": 0, "see": 0, "load": [0, 1], "process": 0, "With": [0, 2], "separ": 0, "python": [0, 2], "some": 0, "custom": 0, "like": 0, "incompat": 0, "other": 0, "In": [0, 1], "those": 0, "case": 0, "spawn": 0, "note": 0, "must": 0, "still": 0, "By": 0, "locat": 0, "possibli": 0, "script": 0, "interpret": 0, "also": 0, "accomplish": 0, "index": 0, "search": 0, "page": 0, "exampl": 1, "molecular": 1, "dynam": 1, "biosimspac": [1, 3], "maiz": 1, "pmemd": 1, "cuda": 1, "from": 1, "pathlib": 1, "import": 1, "path": 1, "ex": 1, "productionpmemdcuda": 1, "io": 1, "loaddata": 1, "return": 1, "util": 1, "jobresourceconfig": 1, "build": 1, "flow": 1, "name": 1, "prod_bss_amber_test": 1, "cleanup_temp": 1, "fals": 1, "level": 1, "debug": 1, "add": 1, "load_si": 1, "list": 1, "prod_pmemd": 1, "production_amb": 1, "paramet": 1, "runtim": 1, "1": 1, "0": 1, "n": 1, "retu": 1, "set": 1, "data": 1, "complex": 1, "prm7": 1, "changem": 1, "rst7": 1, "connect": 1, "out": 1, "inp": 1, "check": 1, "visual": 1, "mol": 1, "get": 1, "bss": 1, "sy": 1, "readmolecul": 1, "str": 1, "print": 1, "40": 1, "realiti": 1, "would": 1, "do": 1, "someth": 1, "http": 2, "openbiosim": 2, "org": 2, "framework": 2, "engin": 2, "agnost": 2, "simul": 2, "implement": 2, "allow": 2, "parameteris": 2, "solvat": 2, "minimis": 2, "equilibr": 2, "product": 2, "rel": 2, "absolut": 2, "bind": 2, "free": 2, "energi": 2, "calcul": 2, "all": 3}, "objects": {}, "objtypes": {}, "objnames": {}, "titleterms": {"maiz": 0, "biosimspac": [0, 2], "instal": 0, "configur": 0, "indic": 0, "tabl": 0, "product": 1, "md": 1, "step": 3}, "envversion": {"sphinx.domains.c": 3, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 9, "sphinx.domains.index": 1, "sphinx.domains.javascript": 3, "sphinx.domains.math": 2, "sphinx.domains.python": 4, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "sphinx.ext.intersphinx": 1, "sphinx.ext.viewcode": 1, "nbsphinx": 4, "sphinx": 60}, "alltitles": {"maize-biosimspace": [[0, "maize-biosimspace"]], "Installation": [[0, "installation"]], "Configuration": [[0, "configuration"]], "Indices and tables": [[0, "indices-and-tables"]], "Production MD": [[1, "production-md"]], "BioSimSpace": [[2, "biosimspace"]], "biosimspace": [[2, "id1"]], "Steps": [[3, "steps"], [3, null]]}, "indexentries": {}})
\ No newline at end of file
+Search.setIndex({"docnames": ["cli_tutorials", "configuration", "index", "installation", "production-md", "python_tutorials", "steps/biosimspace", "steps/index", "tutorial_cli_abfe", "tutorial_cli_production_md", "tutorial_cli_system_preparation", "tutorial_python_production"], "filenames": ["cli_tutorials.rst", "configuration.rst", "index.rst", "installation.rst", "production-md.rst", "python_tutorials.rst", "steps/biosimspace.rst", "steps/index.rst", "tutorial_cli_abfe.rst", "tutorial_cli_production_md.rst", "tutorial_cli_system_preparation.rst", "tutorial_python_production.rst"], "titles": ["CLI Tutorials", "Configuration", "maize-biosimspace", "Installation", "Production MD", "Python Tutorials", "BioSimSpace", "Steps", "Absolute Binding Free Energy Calculations", "Production Molecular Dynamics", "System Preparation", "Creating a Workflow with Production MD"], "terms": {"maiz": [0, 1, 3, 4, 5, 8, 9, 10, 11], "biosimspac": [0, 3, 4, 5, 7, 8, 9, 10, 11], "provid": [0, 2, 6, 10], "larg": [0, 10], "number": [0, 1, 10], "node": [0, 1, 2, 3, 4, 5, 6, 10, 11], "perform": 0, "common": [0, 10], "oper": 0, "parameteris": [0, 6, 10], "minimis": [0, 6, 10], "equilibr": [0, 6, 10], "product": [0, 2, 5, 6, 10], "molecular": [0, 4, 10, 11], "dynam": [0, 4, 10, 11], "alchem": 0, "free": [0, 6, 10], "energi": [0, 6, 10], "calcul": [0, 6, 10], "ar": [0, 2, 3, 6, 10], "gener": 0, "engin": [0, 6], "specif": 0, "requir": [0, 2, 3, 8, 10], "have": [0, 9, 10], "list": [0, 2, 4, 11], "required_cal": 0, "whose": 0, "presenc": 0, "environ": [0, 2, 3], "i": [0, 2, 6, 7, 8, 9, 10], "check": [0, 4, 8, 9, 10, 11], "befor": [0, 10, 11], "run": [0, 1, 2, 3, 4, 8, 9, 10, 11], "To": [0, 1, 2, 9, 10], "see": [0, 1, 2, 10], "all": [0, 2, 5, 7, 8, 9, 10], "avail": [0, 2, 3, 5, 8, 9], "type": [0, 2, 9, 10], "bss_": [0, 2], "hit": [0, 2, 9], "tab": [0, 2, 9], "them": [0, 2, 10], "option": [0, 2, 5, 8], "each": [0, 2], "pass": [0, 9], "h": [0, 2, 8, 9], "flag": [0, 3], "e": [0, 2, 3, 10], "g": [0, 2, 10], "bss_parameteris": [0, 2], "These": 0, "give": 0, "exampl": [0, 4, 8, 11], "us": [0, 1, 2, 3, 4, 5, 6, 8, 9, 10, 11], "creat": [0, 1, 2, 3, 5], "system": [0, 1, 4, 11], "start": 0, "from": [0, 4, 5, 8, 11], "unparameteris": 0, "input": [0, 8, 9, 10], "structur": [0, 10], "absolut": [0, 6, 10], "bind": [0, 6, 10], "protein": [0, 2, 8, 9, 10], "pdb": [0, 8, 10], "an": [0, 4, 8, 10], "sdf": [0, 8, 10], "file": [0, 2, 8, 9, 10], "contain": [0, 1, 2, 7, 8, 9], "multipl": 0, "ligand": [0, 2, 8, 9, 10], "prepar": 0, "For": [1, 3, 8], "detail": [1, 10], "guid": 1, "pleas": 1, "document": [1, 2, 7], "contrib": [1, 5], "By": 1, "default": [1, 9, 10], "look": [1, 11], "xdg_config_hom": 1, "toml": 1, "config": 1, "inform": 1, "set": [1, 4, 8, 10, 11], "up": [1, 3, 10], "through": [1, 2, 4, 5, 8, 9, 11], "slurm": [1, 2], "export": 1, "Then": 1, "follow": 1, "can": [1, 2, 5, 10], "one": 1, "cobalt": 1, "flux": 1, "local": [1, 2, 10], "lsf": 1, "pbspro": 1, "rp": 1, "max_job": 1, "100": 1, "The": [1, 8, 10], "maximum": 1, "job": 1, "submit": [1, 10], "onc": 1, "queue": [1, 10], "gpu": [1, 10], "changem": [1, 4, 11], "your": [1, 2, 3, 10], "desir": 1, "launcher": 1, "srun": 1, "command": [1, 8], "usual": 1, "mpirun": 1, "mpiexec": 1, "walltim": 1, "24": 1, "00": 1, "limit": 1, "shorter": 1, "time": [1, 10], "mai": 1, "improv": 1, "now": [1, 2, 8, 9, 10], "when": 1, "you": [1, 2, 3, 4, 5, 9, 10, 11], "execut": [1, 2, 3, 4, 8, 11], "workflow": [1, 2, 4, 5, 6, 8, 10], "schedul": [1, 2], "graph": [2, 4, 8, 11], "base": 2, "manag": 2, "comput": 2, "chemistri": 2, "pipelin": 2, "thi": [2, 3, 6, 7, 8, 9, 10], "repositori": [2, 3], "namespac": 2, "packag": [2, 6], "extens": 2, "find": [2, 10], "core": [2, 3, 4, 11], "here": [2, 4, 8, 9, 10, 11], "clone": [2, 3], "instal": 2, "git": [2, 3], "http": [2, 3, 6], "github": [2, 3], "com": [2, 3], "fjclark": [2, 3], "cd": [2, 3, 8, 9, 10], "mamba": [2, 3], "env": [2, 3], "f": [2, 3], "user": [2, 3], "yml": [2, 3], "activ": [2, 3], "pip": [2, 3], "dep": [2, 3], "ensur": [2, 3], "amber": [2, 3, 9], "gromac": [2, 3, 9], "namd": [2, 3], "plan": [2, 3], "relev": [2, 3], "mani": [2, 8], "cli": [2, 5, 8, 9, 10], "try": 2, "0": [2, 4, 8, 9, 11], "1": [2, 4, 8, 9, 11], "n": [2, 4, 9, 11], "md": [2, 5], "complex": [2, 4, 5, 9, 10, 11], "dump": 2, "intermedi": [2, 9], "current": [2, 9], "directori": [2, 3, 9], "mkdir": [2, 8, 9, 10], "gmx_md_exampl": [2, 9], "cp": [2, 8, 9, 10], "test": [2, 8, 9, 10], "data": [2, 4, 8, 9, 10, 11], "bss_production_gromac": [2, 9], "inp": [2, 4, 9, 10, 11], "prm7": [2, 4, 9, 11], "rst7": [2, 4, 9, 10, 11], "runtim": [2, 4, 9, 11], "save_nam": [2, 9, 10], "gmx_md_out": [2, 9], "dump_to": [2, 9], "configur": [2, 10], "index": 2, "modul": 2, "search": 2, "page": 2, "basic": [3, 11], "simpli": [3, 8], "If": [3, 10], "want": [3, 11], "keep": 3, "date": 3, "latest": 3, "chang": 3, "switch": 3, "same": 3, "conda": 3, "develop": 3, "should": [3, 8, 9], "dev": 3, "instead": 3, "pmemd": [4, 11], "cuda": [4, 11], "pathlib": [4, 11], "import": [4, 8, 11], "path": [4, 11], "step": [4, 10, 11], "ex": [4, 8, 11], "productionpmemdcuda": [4, 11], "io": [4, 11], "loaddata": [4, 11], "return": [4, 11], "util": [4, 11], "jobresourceconfig": [4, 11], "build": [4, 11], "flow": [4, 11], "name": [4, 9, 10, 11], "prod_bss_amber_test": [4, 11], "cleanup_temp": [4, 11], "fals": [4, 11], "level": [4, 11], "debug": [4, 11], "add": [4, 11], "load_si": [4, 11], "prod_pmemd": [4, 11], "production_amb": [4, 11], "paramet": [4, 10, 11], "retu": [4, 11], "connect": [4, 11], "out": [4, 10, 11], "visual": [4, 11], "mol": [4, 8, 11], "get": [4, 10, 11], "load": [4, 10, 11], "bss": [4, 11], "sy": [4, 11], "readmolecul": [4, 11], "str": [4, 11], "print": [4, 11], "40": [4, 11], "In": [4, 11], "realiti": [4, 11], "would": [4, 10, 11], "do": [4, 10, 11], "someth": [4, 11], "subgraph": 5, "combin": 5, "more": [5, 10], "benefit": 5, "other": [5, 10], "also": [5, 10], "pre": [5, 8], "made": [5, 8], "gain": 5, "fuller": 5, "control": 5, "over": [5, 8, 10], "compar": 5, "write": [5, 8], "reusabl": 5, "script": [5, 8, 9, 10], "python": [6, 8], "openbiosim": 6, "org": 6, "With": 6, "framework": 6, "agnost": 6, "simul": [6, 10], "implement": 6, "allow": 6, "solvat": [6, 10], "rel": [6, 8], "we": [8, 9, 10, 11], "ll": [8, 9, 10, 11], "quick": [8, 9], "benzen": [8, 10], "bound": [8, 10], "t4": [8, 10], "lysozym": [8, 10], "bss_abfe_multi_isom": 8, "its": 8, "which": [8, 10, 11], "onli": 8, "copi": [8, 9, 10], "bss_abfe_exampl": 8, "t4l": [8, 10], "let": 8, "": 8, "veri": [8, 10], "short": 8, "still": 8, "fairli": 8, "expens": 8, "abf": 8, "2": 8, "replic": 8, "lig_sdfs_fil": 8, "protein_pdb": [8, 10], "ligand_force_field": [8, 10], "gaff2": [8, 10], "protein_force_field": [8, 10], "ff14sb": [8, 10], "abfe_timestep": 8, "4": 8, "abfe_n_repl": 8, "abfe_runtim": 8, "abfe_runtime_generate_boresch_restraint": 8, "prep_runtime_restrained_npt": 8, "05": 8, "prep_runtime_unrestrained_npt": 8, "abfe_estim": 8, "ti": 8, "results_file_nam": 8, "abfe_out": 8, "show": 8, "result": 8, "around": 8, "kcal": 8, "line": 8, "argument": [8, 9], "unweildi": 8, "some": [8, 10, 11], "aren": 8, "t": 8, "lambda": 8, "space": 8, "It": 8, "like": [8, 10], "better": 8, "directli": [8, 10], "factori": 8, "customis": 8, "af": 8, "getabfe_multi_isomer_workflow": 8, "bss_product": 9, "bss_production_pmemd": 9, "bss_production_pmemd_cuda": 9, "bss_production_sand": 9, "bss_production_somd": 9, "pick": 9, "includ": 9, "specifi": 9, "output": [9, 10], "save": [9, 10], "log": 9, "trajectori": 9, "subdirectori": 9, "work": [9, 10], "final": 9, "coordin": 9, "topologi": 9, "sub": 9, "note": 9, "despit": 9, "were": 9, "abl": 9, "becaus": 9, "automat": 9, "convert": 9, "between": 9, "format": 9, "sire": 9, "under": 9, "hood": 9, "often": 10, "heat": 10, "obtain": 10, "suitabl": 10, "while": 10, "individu": 10, "two": 10, "complet": 10, "bss_system_prep_fre": 10, "design": 10, "box": 10, "water": 10, "apo": 10, "bss_system_prep_bound": 10, "l99a": 10, "first": 10, "sysprep_bound_exampl": 10, "ha": 10, "been": 10, "sanitis": 10, "tweak": 10, "thei": 10, "accept": 10, "tleap": 10, "behind": 10, "scene": 10, "recommend": 10, "clean": 10, "unsanitis": 10, "pdb4amber": 10, "o": 10, "protein_sanitis": 10, "attempt": 10, "below": 10, "rais": 10, "error": 10, "messag": 10, "sourc": 10, "leaprc": 10, "tip3p": 10, "prot": 10, "loadpdb": 10, "saveamberparm": 10, "parm7": 10, "savepdb": 10, "protein_fully_sanitis": 10, "abov": 10, "fail": 10, "due": 10, "residu": 10, "atom": 10, "match": 10, "templat": 10, "read": 10, "caus": 10, "issu": 10, "expect": 10, "librari": 10, "wa": 10, "after": 10, "cat": 10, "amberhom": 10, "dat": 10, "leap": 10, "lib": 10, "amino12": 10, "renam": 10, "offend": 10, "repeat": 10, "fussi": 10, "about": 10, "ani": 10, "warn": 10, "round": 10, "protein_full_sanitis": 10, "our": 10, "t4l_benzene_complex_equilibr": 10, "forc": 10, "field": 10, "There": 10, "modifi": 10, "make": 10, "sure": 10, "access": 10, "t4l_benezene_complex_equilibr": 10, "custom": 11, "extra": 11, "afterward": 11}, "objects": {}, "objtypes": {}, "objnames": {}, "titleterms": {"cli": 0, "tutori": [0, 5], "configur": 1, "maiz": 2, "biosimspac": [2, 6], "quick": 2, "start": 2, "indic": 2, "tabl": 2, "instal": 3, "product": [4, 9, 11], "md": [4, 11], "python": 5, "step": 7, "absolut": 8, "bind": 8, "free": 8, "energi": 8, "calcul": 8, "molecular": 9, "dynam": 9, "system": 10, "prepar": 10, "creat": 11, "workflow": 11}, "envversion": {"sphinx.domains.c": 3, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 9, "sphinx.domains.index": 1, "sphinx.domains.javascript": 3, "sphinx.domains.math": 2, "sphinx.domains.python": 4, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "sphinx.ext.intersphinx": 1, "sphinx.ext.viewcode": 1, "nbsphinx": 4, "sphinx": 60}, "alltitles": {"CLI Tutorials": [[0, "cli-tutorials"]], "Configuration": [[1, "configuration"]], "maize-biosimspace": [[2, "maize-biosimspace"]], "Quick Start": [[2, "quick-start"]], "Indices and tables": [[2, "indices-and-tables"]], "Installation": [[3, "installation"]], "Production MD": [[4, "production-md"]], "Python Tutorials": [[5, "python-tutorials"]], "BioSimSpace": [[6, "biosimspace"]], "biosimspace": [[6, "id1"]], "Steps": [[7, "steps"], [7, null]], "Absolute Binding Free Energy Calculations": [[8, "absolute-binding-free-energy-calculations"]], "Production Molecular Dynamics": [[9, "production-molecular-dynamics"]], "Creating a Workflow with Production MD": [[11, "creating-a-workflow-with-production-md"]], "System Preparation": [[10, "system-preparation"]]}, "indexentries": {}})
\ No newline at end of file
diff --git a/_build/steps/biosimspace.html b/_build/steps/biosimspace.html
index f665d4c..dbf551a 100644
--- a/_build/steps/biosimspace.html
+++ b/_build/steps/biosimspace.html
@@ -160,13 +160,27 @@
   
 
 
            
-  
            Examples
            
+  
            Setup
            
 
+
            Tutorials
            
+
 
            Reference
            
 
              
-
            • Steps
                
+
              • Steps
 
                • 
diff --git a/_build/steps/index.html b/_build/steps/index.html
index 98e76e3..06505a3 100644
--- a/_build/steps/index.html
+++ b/_build/steps/index.html
@@ -3,7 +3,7 @@
   
     
     
-
+
 
     
         Steps - maize biosimspace
@@ -160,13 +160,27 @@
   
 
 
                
-  
                Examples
                
+  
                Setup
                
 
+
                Tutorials
                
+
 
                Reference
                
 
                
               
             
-          
+          
               
               
                
                 
                
                   Previous
                 
                
                 
-                
                Production MD
                
+                
                Creating a Workflow with Production MD
                
                 
               
                
                             
diff --git a/_build/tutorial_cli_abfe.html b/_build/tutorial_cli_abfe.html
new file mode 100644
index 0000000..013a024
--- /dev/null
+++ b/_build/tutorial_cli_abfe.html
@@ -0,0 +1,329 @@
+
+
+  
+    
+    
+
+
+    
+        Absolute Binding Free Energy Calculations - maize biosimspace
+      
+    
+    
+    
+    
+
+
+
+  
+    
+    
+    
+
+
+  
+    Contents
+    
+      
+    
+  
+  
+    Menu
+    
+      
+      
+      
+    
+  
+  
+    Expand
+    
+      
+    
+  
+  
+    Light mode
+    
+      
+      
+      
+      
+      
+      
+      
+      
+      
+    
+  
+  
+    Dark mode
+    
+      
+      
+    
+  
+  
+    Auto light/dark mode
+    
+      
+      
+      
+      
+      
+      
+      
+    
+  
+
+
+
+
+
+
+
+
+
+
                
+  
                
+    
                
+      
+    
                
+    
                
+      
                maize biosimspace
                
+    
                
+    
                
+      
                
+        
+      
                
+      
+    
                
+  
                
+  
+  
                
+    
                
+      
                
+        
+          
+            
+          
+          Back to top
+        
+        
                
+          
+
                
+            
+          
                
+          
+        
                
+        
                
+          
                
+
                Absolute Binding Free Energy Calculations#
                
+
                Here, we’ll run a quick absolute binding free energy calculation for benzene bound to T4 Lysozyme. For this, we’ll
+use the bss_abfe_multi_isomer workflow through its CLI, which requires only an SDF file containing all required
+ligands, and the pdb of the protein. Check the options with
                
+
                bss_abfe_multi_isomer -h
+
                
+
                
+
                Copy over the required input:
                
+
                mkdir bss_abfe_example
+cd bss_abfe_example
+cp ../tests/data/benzene.sdf .
+cp ../tests/data/t4l.pdb .
+
                
+
                
+
                Now, let’s run a relatively very short (but still fairly expensive)
+ABFE calculation with 2 replicates:
                
+
                bss_abfe_multi_isomer --lig_sdfs_file benzene.sdf \
+                      --protein_pdb benzene.t4l \
+                      --ligand_force_field gaff2 \
+                      --protein_force_field ff14SB \
+                      --abfe_timestep 4 \
+                      --abfe_n_replicates 2 \
+                      --abfe_runtime 0.1 \
+                      --abfe_runtime_generate_boresch_restraint 0.1 \
+                      --prep_runtime_restrained_npt 0.05 \
+                      --prep_runtime_unrestrained_npt 0.05 \
+                      --abfe_estimator TI \
+                      --results_file_name abfe_out \
+
                
+
                
+
                The abfe_out file should show results around -4 kcal / mol.
                
+
                Running through the command line with this many arguments is unweildy,
+and some options aren’t available through the CLI (for example, the lambda
+spacing). It’s likely a better option to write a quick script - using the
+pre-made workflow directly in python - simply import the workflow factory
+, customise the options, and run (all in a python script):
                
+
                from maize.graphs.exs.biosimspace.afe import getabfe_multi_isomer_workflow
+workflow = getabfe_multi_isomer_workflow()
+
+# Set workflow options...
+
+# Run
+workflow.execute()
+
                
+
                
+
                
+
+        
                
+      
                
+      
+    
                
+    
+  
                
+
                
+    
+    
+    
+    
+    
+
\ No newline at end of file
diff --git a/_build/tutorial_cli_production_md.html b/_build/tutorial_cli_production_md.html
new file mode 100644
index 0000000..a71bebf
--- /dev/null
+++ b/_build/tutorial_cli_production_md.html
@@ -0,0 +1,305 @@
+
+
+  
+    
+    
+
+
+    
+        Production Molecular Dynamics - maize biosimspace
+      
+    
+    
+    
+    
+
+
+
+  
+    
+    
+    
+
+
+  
+    Contents
+    
+      
+    
+  
+  
+    Menu
+    
+      
+      
+      
+    
+  
+  
+    Expand
+    
+      
+    
+  
+  
+    Light mode
+    
+      
+      
+      
+      
+      
+      
+      
+      
+      
+    
+  
+  
+    Dark mode
+    
+      
+      
+    
+  
+  
+    Auto light/dark mode
+    
+      
+      
+      
+      
+      
+      
+      
+    
+  
+
+
+
+
+
+
+
+
+
+
                
+  
                
+    
                
+      
+    
                
+    
                
+      
                maize biosimspace
                
+    
                
+    
                
+      
                
+        
+      
                
+      
+    
                
+  
                
+  
+  
                
+    
                
+      
                
+        
+          
+            
+          
+          Back to top
+        
+        
                
+          
+
                
+            
+          
                
+          
+        
                
+        
                
+          
                
+
                Production Molecular Dynamics#
                
+
                Here, we’ll run production molecular dynamics on a protein-ligand complex. To check the available production CLIs, type bss_production and hit tab:
                
+
                bss_production_gromacs     bss_production_pmemd       bss_production_pmemd_cuda  bss_production_sander      bss_production_somd
+
                
+
                
+
                We’ll pick gromacs. To check the available arguments and defaults, run
                
+
                bss_production_gromacs -h
+
                
+
                
+
                We’ll run a quick 0.1 ns of production molecular dynamics on the protein-ligand complex included with maize-biosimspace for testing. We’ll specify the output
+name to be gmx_md_out and we’ll save all of the intermediate files (including input scripts, logs, and trajectory files) to a subdirectory in the current
+working directory by specifying --dump_to .
                
+
                mkdir gmx_md_example
+cd gmx_md_example
+cp ../tests/data/complex.* .
+bss_production_gromacs --inp complex.prm7 complex.rst7 --runtime 0.1 --save_name gmx_md_out --dump_to .
+
                
+
                
+
                You should now have the final coordinate file, gmx_md_out.rst7, a copy of the input topology file gmx_md_out.prm7, and a sub-directory containing all of the
+intermediate files. Note that despite running through GROMACS, we were able to pass in AMBER files as input. This is because BioSimSpace automatically converts
+between file formats (using Sire under the hood).
                
+
                
+
+        
                
+      
                
+      
+    
                
+    
+  
                
+
                
+    
+    
+    
+    
+    
+
\ No newline at end of file
diff --git a/_build/tutorial_cli_system_preparation.html b/_build/tutorial_cli_system_preparation.html
new file mode 100644
index 0000000..081450e
--- /dev/null
+++ b/_build/tutorial_cli_system_preparation.html
@@ -0,0 +1,330 @@
+
+
+  
+    
+    
+
+
+    
+        System Preparation - maize biosimspace
+      
+    
+    
+    
+    
+
+
+
+  
+    
+    
+    
+
+
+  
+    Contents
+    
+      
+    
+  
+  
+    Menu
+    
+      
+      
+      
+    
+  
+  
+    Expand
+    
+      
+    
+  
+  
+    Light mode
+    
+      
+      
+      
+      
+      
+      
+      
+      
+      
+    
+  
+  
+    Dark mode
+    
+      
+      
+    
+  
+  
+    Auto light/dark mode
+    
+      
+      
+      
+      
+      
+      
+      
+    
+  
+
+
+
+
+
+
+
+
+
+
                
+  
                
+    
                
+      
+    
                
+    
                
+      
                maize biosimspace
                
+    
                
+    
                
+      
                
+        
+      
                
+      
+    
                
+  
                
+  
+  
                
+    
                
+      
                
+        
+          
+            
+          
+          Back to top
+        
+        
                
+          
+
                
+            
+          
                
+          
+        
                
+        
                
+          
                
+
                System Preparation#
                
+
                Often, we have structure files for a protein and / or ligand and we would like to parameterise, solvate, minimise, heat, and equilibrate them to obtain systems
+suitable for production molecular dynamics simulations or free energy calculations. While individual nodes are provided for all of these steps,
+maize-biosimspace also provides CLIs for two complete system preparation workflows: bss_system_prep_free, which is designed to prepare a ligand in a box
+of water (but can also be used for an apo protein) and bss_system_prep_bound, which is designed to set up protein-ligand complexes.
                
+
                Here, we’ll prepare the complex of T4 lysozyme L99A bound to benzene, a very common test system for absolute binding free energy calculations. First, we’ll copy
+over the required input:
                
+
                mkdir sysprep_bound_example
+cd sysprep_bound_example
+cp ../tests/data/benzene.sdf .
+cp ../tests/data/t4l.pdb .
+
                
+
                
+
                
+
                Tip
                
+
                This pdb has been sanitised and will work with BioSimSpace first time, but often pdbs will require some tweaking before they are accepted by tleap (which
+BioSimSpace uses behind the scenes). The recommended workflow is
                
+
                  
+
                • Clean your unsanitised pdb using pdb4amber, e.g. pdb4amber -i protein.pdb -o protein_sanitised.pdb
                  • 
+
                • Attempt to run the workflow below
                  • 
+
                • If the workflow raises an error, attempt parameterisation directly with tleap to get more detailed error messages. E.g., type tleap, then
                  • 
+
                
+
                source leaprc.protein.ff14SB
+source leaprc.water.tip3p
+# Loading an unsanitised pdb will likely raise an error
+prot = loadpdb protein_sanitised.pdb
+saveamberparm prot protein.parm7 protein.rst7
+savepdb prot protein_fully_sanitised.pdb
+
                
+
                
+
                  
+
                • If the above fails, this is often due to residue/ atom names which do not match the templates. Read the errors to find out which residues / atoms are causing the issues, then check the expected names in library which was loaded after typing source leaprc.protein.ff14SB e.g. cat $AMBERHOME/dat/leap/lib/amino12.lib.  Rename the offending atoms/ residues and repeat the above step.
                  • 
+
                
+
                BioSimSpace is very fussy about parameterisation and will fail if tleap raises any warnings. To get round this, run the tleap script above and use the
+protein_full_sanitised.pdb file as your input, which will not raise any errors.
                
+
                
+
                To run system preparation for our protein-ligand complex, we’ll use the bss_system_prep_bound CLI, saving the output system to “t4l_benzene_complex_equilibrated”
+and using the gaff2 ff14SB force fields. There are a large number of other parameters which can be modified (see bss_system_prep_bound -) but we’ll run the
+defaults for now.
                
+
                
+
                Note
                
+
                Make sure that you have access to a GPU locally, or have configured Maize to submit to a queue with gpu access (see Configuration)
                
+
                
+
                bss_system_prep_bound --inp benzene.sdf --protein_pdb t4l.pdb --ligand_force_field gaff2 --protein_force_field ff14SB --save_name t4l_benezene_complex_equilibrated
+
                
+
                
+
                
+
+        
                
+      
                
+      
+    
                
+    
+  
                
+
                
+    
+    
+    
+    
+    
+
\ No newline at end of file
diff --git a/_build/tutorial_python_production.html b/_build/tutorial_python_production.html
new file mode 100644
index 0000000..7954106
--- /dev/null
+++ b/_build/tutorial_python_production.html
@@ -0,0 +1,343 @@
+
+
+  
+    
+    
+
+
+    
+        Creating a Workflow with Production MD - maize biosimspace
+      
+    
+    
+    
+    
+
+
+
+  
+    
+    
+    
+
+
+  
+    Contents
+    
+      
+    
+  
+  
+    Menu
+    
+      
+      
+      
+    
+  
+  
+    Expand
+    
+      
+    
+  
+  
+    Light mode
+    
+      
+      
+      
+      
+      
+      
+      
+      
+      
+    
+  
+  
+    Dark mode
+    
+      
+      
+    
+  
+  
+    Auto light/dark mode
+    
+      
+      
+      
+      
+      
+      
+      
+    
+  
+
+
+
+
+
+
+
+
+
+
                
+  
                
+    
                
+      
+    
                
+    
                
+      
                maize biosimspace
                
+    
                
+    
                
+      
                
+        
+      
                
+      
+    
                
+  
                
+  
+  
                
+    
                
+      
                
+        
+          
+            
+          
+          Back to top
+        
+        
                
+          
+
                
+            
+          
                
+          
+        
                
+        
                
+          
                
+
                Creating a Workflow with Production MD#
                
+
                Here, we’ll look at a basic example of a custom workflow which
+uses pmemd.cuda to run some production MD. In reality, you
+would want to add some extra steps before or afterwards:
                
+
                """Run production Molecular Dynamics using PMEMD.CUDA through BioSimSpace."""
+
+from pathlib import Path
+
+from maize.core.workflow import Workflow
+from maize.steps.exs.biosimspace import ProductionPmemdCuda
+from maize.steps.io import LoadData, Return
+from maize.utilities.execution import JobResourceConfig
+
+# Build the graph
+flow = Workflow(name="Prod_BSS_AMBER_Test", cleanup_temp=False, level="debug")
+
+# Add the nodes
+load_sys = flow.add(LoadData[list[Path]])
+prod_pmemd = flow.add(
+    ProductionPmemdCuda,
+    name="Production_Amber",
+    parameters={
+        "runtime": 1.0, # ns
+    },
+)
+retu = flow.add(Return[list[Path]])
+
+# Set parameters
+load_sys.data.set(
+    [
+        Path(
+            "< path to complex.prm7>" # CHANGEME
+        ),
+        Path(
+            "< path to complex.rst7>" # CHANGEME
+        ),
+    ]
+)
+
+# Connect the nodes
+flow.connect(load_sys.out, prod_pmemd.inp)
+flow.connect(prod_pmemd.out, retu.inp)
+
+# Check and run!
+flow.check()
+flow.visualize()
+flow.execute()
+
+mols = retu.get()
+
+# Load a BioSimSpace system from the returned paths
+import BioSimSpace as BSS
+
+sys = BSS.IO.readMolecules([str(mols[0]), str(mols[1])])
+print(40 * "#")
+print(sys)
+# In reality, you would do something here...
+
                
+
                
+
                
+
+        
                
+      
                
+      
+    
                
+    
+  
                
+
                
+    
+    
+    
+    
+    
+
\ No newline at end of file