BioExcel Building Blocks Workflows list
BioExcel Building Blocks Workflows
In this page, users can find a list of all available BioExcel Building Blocks Workflows. A list of main steps and a direct launch button are provided for each workflow.
They can also be found in a GitHub repository in all their versions with all the instructions and sample files needed for executing them at home in the users' own premises.
Please help us to improve the BioExcel Building Blocks Workflows by filling in our 2023 survey.
This workflow provides a pipeline to setup DNA structures following the recommended guidelines by the Ascona B-DNA Consortium (ABC) members. It follows the work started with the NAFlex tool to offer a single, reproducible pipeline for structure preparation, ensuring reproducibility and coherence between all the members of the consortium.
- Modeling Nucleic Acids 3D Structure
- Generating Topology
- Adding Water Box (with given input parameters. Default: 15Å, truncated octahedron box, opcbox water model (SPC/E + Joung-Chetham monovalent ions + Li/Merz highly charged ions)
- Neutralizing the system (with K+ ions) and Adding additional Ionic Concentration (with given input parameters. Default: 100 Mol/L)
- Randomizing Ions
- Generating Topology with Hydrogen Mass Partitioning (allowing 4fs timestep MD)
- Equilibrating the System (10-steps protocol Daniel R. Roe and Bernard R. Brooks):
- Eq Step 1 -- System Energetic Minimization, 5 Kcal/mol heavy atoms restraints (1000 cycles)
- Eq Step 2 -- NVT Equilibration, 5 Kcal/mol heavy atoms restraints, timestep 1fs (15ps)
- Eq Step 3 -- System Energetic Minimization, 2 Kcal/mol heavy atoms restraints (1000 cycles)
- Eq Step 4 -- System Energetic Minimization, 0.1 Kcal/mol heavy atoms restraints (1000 cycles)
- Eq Step 5 -- System Energetic Minimization (1000 cycles)
- Eq Step 6 -- NPT Equilibration, 1 Kcal/mol heavy atoms restraints, timestep 1fs (5ps)
- Eq Step 7 -- NPT Equilibration, 0.5 Kcal/mol heavy atoms restraints, timestep 1fs (5ps)
- Eq Step 8 -- NPT Equilibration, 0.5 Kcal/mol backbone atoms restraints, timestep 1fs (10ps)
- Eq Step 9 -- NPT Equilibration, timestep 2fs, unrestrained MD simulation, writing every 1ps (20ps, 20 frames)
- Eq Step 10 -- NPT Equilibration, timestep 2fs, unrestrained MD simulation, writing every 1ps (20ps, 20 frames)
- Running Final Free (unrestrained) Short MD Simulation, timestep 4fs, unrestrained MD simulation, writing every 1ps (500ps, 500 frames)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, the nucleic acid 3D structure is previously computed by BioBB Workflows.
This workflow performs an automatic ligand parameterization for a small molecule using GAFF force field, generating parameters compatible with the AMBER MD package.
- Energetically Minimize Hydrogen Atoms (with given input parameters. Default: Steepest Descent algorithm)
- Generating Ligand Parameters (with given input parameters. Default: GAFF ff, pH 7.4, charge depending on the protonation state)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, small molecule protonation states are previously computed by BioBB Workflows, including the major (de)protonated form of the molecule at a specified pH 7.4 (major microspecies) and a list of (de)protonated forms of the molecule at a specified pH in range from 0 to 14.
Launch AMBER Ligand Parameterization workflow
This workflow performs a simulation setup of a protein-ligand(s) complex system, compatible with the AMBER MD package.
- Extracting Protein Structure
- Preparing Protein PDB structure for AMBER
- Creating Protein-Ligand Complex Topology (with given input parameters. Default: protein.ff14SB + gaff ffs)
- Energetically Minimizing the Structure
- Step 1: Hydrogen minimization, 500 cycles, applying position restraints (50 Kcal/mol.Å2) to the protein heavy atoms.
- Step 2: System minimization, 500 cycles, applying position restraints (500 Kcal/mol.Å2) to the small molecule.
- Creating Solvent Box and Solvating the System (with given input parameters. Default: truncated octahedron, 9Å, TIP3P water model)
- Neutralizing the system and Adding an additional Ionic Concentration (Sodium (Na+) and Chloride (Cl-) counterions with a given input additional ionic concentration. Default: 150mM.)
- Energetically Minimizing the System (Steepest descent algorithm, 500 cycles, applying position restraints (15 Kcal/mol.Å2) to the solute atoms)
- Heating the System (From 0 to 300K. Solute atoms restrained (force constant of 10 Kcal/mol.Å2). Length 100ps.)
- Equilibrating the System (NVT. Heavy atoms restrained (force constant of 5 Kcal/mol.Å2). 25000 steps, 2fs timestep. Length 50ps)
- Equilibrating the System (NPT. Heavy atoms restrained (force constant of 2.5 Kcal/mol.Å2). 25000 steps, 2fs timestep. Length 50ps)
- Free (unrestrained) Molecular Dynamics Simulation (NPT. Unrestrained. 250000 steps, 2fs timestep. Length 500ps)
- Post-processing and Visualizing Resulting 3D Trajectory (Remove PBC, RMSd, Rgyr)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, an exhaustive protein structure checking is previously run by BioBB Workflows.
Launch AMBER Complex Set Up workflow
This workflow performs a simulation setup of a protein system, compatible with the AMBER MD package.
- Preparing Protein PDB structure for AMBER
- Creating Protein Topology (with given input parameters. Default: protein.ff14SB)
- Energetically Minimizing the Structure
- Step 1: Hydrogen minimization, 500 cycles, applying position restraints (50 Kcal/mol.Å2) to the protein heavy atoms.
- Step 2: System minimization, 500 cycles, applying position restraints (50 Kcal/mol.Å2) to the protein heavy atoms.
- Creating Solvent Box and Solvating the System (with given input parameters. Default: truncated octahedron, 9Å, TIP3P water model)
- Neutralizing the system and Adding an additional Ionic Concentration (Sodium (Na+) and Chloride (Cl-) counterions with a given input additional ionic concentration. Default: 150mM.)
- Energetically Minimizing the System (Steepest descent algorithm, 500 cycles, applying position restraints (50 Kcal/mol.Å2) to the solute atoms)
- Heating the System (From 0 to 300K. Solute atoms restrained (force constant of 10 Kcal/mol.Å2). Length 100ps.)
- Equilibrating the System (NVT. Heavy atoms restrained (force constant of 5 Kcal/mol.Å2). 25000 steps, 2fs timestep. Length 50ps)
- Equilibrating the System (NPT. Heavy atoms restrained (force constant of 2.5 Kcal/mol.Å2). 25000 steps, 2fs timestep. Length 50ps)
- Free (unrestrained) Molecular Dynamics Simulation (NPT. Unrestrained. 250000 steps, 2fs timestep. Length 500ps)
- Post-processing and Visualizing Resulting 3D Trajectory (Remove PBC, RMSd, Rgyr)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, an exhaustive protein structure checking is previously run by BioBB Workflows.
This workflow performs an automatic ligand parameterization for a small molecule using GAFF force field, generating parameters compatible with the CNS/XPLOR MD package.
- Energetically Minimize Hydrogen Atoms (with given input parameters. Default: Steepest Descent algorithm)
- Generating Ligand Parameters (with given input parameters. Default: GAFF ff, pH 7.4, charge depending on the protonation state)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, small molecule protonation states are previously computed by BioBB Workflows, including the major (de)protonated form of the molecule at a specified pH 7.4 (major microspecies) and a list of (de)protonated forms of the molecule at a specified pH in range from 0 to 14.
Launch CNS/XPLOR Ligand Parameterization workflow
This workflow performs a simulation setup of a protein-ligand complex system, compatible with the GROMACS MD package.
- Fixing Protein Structure
- Creating Protein Topology (with given input parameters. Default: amber99sb-ildn)
- Generating Protein-Ligand complex structure and topology
- Creating Solvent Box and Solvating the System (with given input parameters. Default: cubic, 1nm, SPC water model)
- Neutralizing the system and Adding an additional Ionic Concentration (Sodium (Na+) and Chloride (Cl-) counterions with a given input additional ionic concentration. Default: 50mM.)
- Energetically Minimizing the System (Steepest descent algorithm, 5000 steps, maximum force placed at 500 KJ/mol.nm2)
- Equilibrating the System (NVT. Heavy atoms restrained (force constant of 1000 KJ/mol.nm2). 50000 steps, 2fs timestep. Length 100ps)
- Equilibrating the System (NPT. Heavy atoms restrained (force constant of 1000 KJ/mol.nm2). 50000 steps, 2fs timestep. Length 100ps)
- Free (unrestrained) Molecular Dynamics Simulation (NPT. Unrestrained. 250000 steps, 2fs timestep. Length 500ps)
- Post-processing and Visualizing Resulting 3D Trajectory (Remove PBC, RMSd, Rgyr)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, an exhaustive protein structure checking is previously run by BioBB Workflows.
Launch GROMACS Complex Set Up workflow
This workflow extracts structural and dynamical properties from a DNA or RNA MD trajectory, including helical parameters, stiffness, bimodalities and correlations.
- Extracting Helical Parameters from input trajectory (with Curves+ and Canal tools)
- Computing Average Helical Parameters
- Base Pair Step (Inter Base Pair) Average Parameters
- Base Pair (Intra Base Pair) Average Parameters
- Axis Average Parameters
- Grooves Average Parameters
- Backbone Torsions Average Parameters
- Computing Time Series Helical Parameters
- Base Pair Step (Inter Base Pair) Time Series Parameters
- Base Pair (Intra Base Pair) Time Series Parameters
- Axis Time Series Parameters
- Grooves Time Series Parameters
- Backbone Torsions Time Series Parameters
- Computing Stiffness
- Base Pair Step (Inter Base Pair) Stiffness
- Base Pair (Intra Base Pair) Stiffness
- Computing Bimodality
- Base Pair Step (Inter Base Pair) Bimodality
- Computing Correlations
- Base Pair Step (Inter Base Pair) Sequence Correlations
- Base Pair (Intra Base Pair) Sequence Correlations
- Base Pair Step (Inter Base Pair) Helical Parameter Correlations
- Base Pair (Intra Base Pair) Helical Parameter Correlations
- Base Pair Step (Inter Base Pair) Neighboring Steps Correlations
- Base Pair (Intra Base Pair) Neighboring Steps Correlations
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, the nucleic acid sequence (needed for the Helical parameter analysis) is previously extracted from the structure by BioBB Workflow. Besides, the trajectory is shown in the NGL visualizer to quickly check for inconsistencies.
Launch DNA helical parameters workflow
This workflow performs a protein-ligand docking with AutoDock Vina.
- Generating Cavity Box from Pocket (with given input parameters. Default: 12Å offset)
- Preparing Small Molecule (ligand) for Docking
- Preparing Target Protein for Docking
- Running the Docking
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, a pocket search on the protein structure is previously run by BioBB Workflows, which helps finding the proper box to use in the Docking process.
Launch Docking Autodock Vina workflow
This workflow performs an automatic ligand parameterization for a small molecule using GAFF force field, generating parameters compatible with the GROMACS MD package.
- Energetically Minimize Hydrogen Atoms (with given input parameters. Default: Steepest Descent algorithm)
- Generating Ligand Parameters (with given input parameters. Default: GAFF ff, pH 7.4, charge depending on the protonation state)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, small molecule protonation states are previously computed by BioBB Workflows, including the major (de)protonated form of the molecule at a specified pH 7.4 (major microspecies) and a list of (de)protonated forms of the molecule at a specified pH in range from 0 to 14.
Launch GROMACS Ligand Parameterization workflow
This workflow performs an automatic ligand parameterization for a small molecule using GAFF force field, generating parameters compatible with the GROMACS MD package and the OPLS/AA force field.
- Energetically Minimize Hydrogen Atoms (with given input parameters. Default: Steepest Descent algorithm)
- Generating Ligand Parameters (with given input parameters. Default: GAFF ff, pH 7.4, charge depending on the protonation state)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, small molecule protonation states are previously computed by BioBB Workflows, including the major (de)protonated form of the molecule at a specified pH 7.4 (major microspecies) and a list of (de)protonated forms of the molecule at a specified pH in range from 0 to 14.
Launch GROMACS OPLS/AA Ligand Parameterization workflow
This workflow computes a set of Quality Control (QC) analyses on top of an uploaded MD trajectory. QC analyses include positional divergence (RMSd), change of shape (Radius of Gyration), identification of flexible regions (atomic/residue fluctuations), and identification of different molecular conformations (trajectory clustering).
- Calculating Average structure
- Computing Root Mean Square deviations (RMSd) against the first structure from the trajectory (equilibrated system) and against the average structure.
- Computing residue fluctuation (Bfactor)
- Computing Radius of Gyration
- Computing cluster families using RMSd between snapshots as metrics
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, the molecule sequence is previously extracted from the structure by BioBB Workflow. Besides, the trajectory is shown in the NGL visualizer to quickly check for inconsistencies.
Launch Protein MD Analysis workflow
This workflow performs a simulation setup of a protein system after modeling one (or many) protein residue mutations. Results are compatible with the GROMACS MD package.
These steps are executed for each mutation.
- Fixing Protein Structure
- Creating Protein Topology (with given input parameters. Default: amber99sb-ildn)
- Creating Solvent Box and Solvating the System (with given input parameters. Default: cubic, 1nm, SPC water model)
- Neutralizing the system and Adding an additional Ionic Concentration (Sodium (Na+) and Chloride (Cl-) counterions with a given input additional ionic concentration. Default: 50mM.)
- Energetically Minimizing the System (Steepest descent algorithm, 5000 steps, maximum force placed at 500 KJ/mol.nm2)
- Equilibrating the System (NVT. Heavy atoms restrained (force constant of 1000 KJ/mol.nm2). 50000 steps, 2fs timestep. Length 100ps)
- Equilibrating the System (NPT. Heavy atoms restrained (force constant of 1000 KJ/mol.nm2). 50000 steps, 2fs timestep. Length 100ps)
- Free (unrestrained) Molecular Dynamics Simulation (NPT. Unrestrained. 250000 steps, 2fs timestep. Length 500ps)
- Post-processing and Visualizing Resulting 3D Trajectory (Remove PBC, RMSd, Rgyr)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, an exhaustive protein structure checking is previously run by BioBB Workflows.
Launch GROMACS Set Up Mutations workflow
This workflow performs a simulation setup of a protein system, compatible with the GROMACS MD package.
- Fixing Protein Structure
- Creating Protein Topology (with given input parameters. Default: amber99sb-ildn)
- Creating Solvent Box and Solvating the System (with given input parameters. Default: cubic, 1nm, SPC water model)
- Neutralizing the system and Adding an additional Ionic Concentration (Sodium (Na+) and Chloride (Cl-) counterions with a given input additional ionic concentration. Default: 50mM.)
- Energetically Minimizing the System (Steepest descent algorithm, 5000 steps, maximum force placed at 500 KJ/mol.nm2)
- Equilibrating the System (NVT. Heavy atoms restrained (force constant of 1000 KJ/mol.nm2). 50000 steps, 2fs timestep. Length 100ps)
- Equilibrating the System (NPT. Heavy atoms restrained (force constant of 1000 KJ/mol.nm2). 50000 steps, 2fs timestep. Length 100ps)
- Free (unrestrained) Molecular Dynamics Simulation (NPT. Unrestrained. 250000 steps, 2fs timestep. Length 500ps)
- Post-processing and Visualizing Resulting 3D Trajectory (Remove PBC, RMSd, Rgyr)
Note: Please check steps done previously to launch this workflow in the corresponding help section. In this case, an exhaustive protein structure checking is previously run by BioBB Workflows.
Launch GROMACS Set Up workflow