Endonuclease PvuII (1PVI) DNA - GATTACAGATTACA
CAP - Catabolite gene Activating Protein (1BER)
DNA - GATTACAGATTACAGATTACA Endonuclease PvuII bound to palindromic DNA recognition site CAGCTG (1PVI) DNA - GATTACAGATTACAGATTACA TBP - TATA box  Binding Protein (1C9B)
CAP - Catabolite gene Activating Protein (1BER)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
TBP - TATA box  Binding Protein (1C9B)

° 

Macros

Macros are written using the Yanaconda language. If you already know a programming language, you can learn Yanaconda in 15 minutes, otherwise it may take 30.

To install, just save the macro in the directory yasara/mcr or wherever you want it. The syntax of Yanaconda is described in the YASARA documentation.

° 

Molecular dynamics

° 

Analyzing the ligand binding energy during a molecular dynamics simulation

This macro calculates the binding energy of a ligand during a molecular dynamics simulation, including the time average. If you didn't use AMBER14 for the MD, adapt the ForceField setting below for consistency. By default, the receptor must be object 1, the ligand must be object 2 (or change the settings). More positive energies indicate better binding, negative energies DO NOT indicate no binding, see the 'BindEnergy' command in the user manual for details, also for words of caution, a much more sophisticated approach is in preparation.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/10/04

° 

Analyzing a molecular dynamics trajectory

This macro analyzes a simulation and creates a detailed report with a large number of plots, e.g. energies, RMSDs, hydrogen bonds. It also tries to identify the main ligand and provides ligand-specific data. All results are additionally written to a simple text table, which can be imported into your favorite spreadsheet program. Your own analysis can often be added with just one line of code, search for 'Example:'.

Written by: Elmar Krieger and Kornel Ozvoldik

License: GNU GPL

Last modified: 2024/10/04

Download: md_analyze.mcr

° 

Running a molecular dynamics simulation of a protein on top of a metal surface with maximum speed

This macro sets up and runs a simulation of a protein on top of a monoelemental face-centerd cubic crystal metal surface as quickly as possible

Written by: Elmar Krieger & Kornel Ozvoldik

License: GNU GPL

Last modified: 2024/10/04

° 

Running a molecular dynamics simulation of a protein on top of a metal surface with normal or fast speed

This macro sets up and runs a simulation of a protein on top of a monoelemental face-centerd cubic crystal metal surface

Written by: Elmar Krieger & Kornel Ozvoldik

License: GNU GPL

Last modified: 2024/10/04

° 

Running the DHFR benchmark to measure MD performance

This macro runs the DHFR benchmark for about 10 seconds and displays the results

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/09/10

° 

Running a molecular dynamics simulation of a membrane protein with normal or fast speed

This macro sets up and runs a simulation of a membrane protein. It scans the protein for secondary structure elements with hydrophobic surface residues, orients it accordingly and embeds it in a membrane of adjustable lipid composition. Finally a 250 ps restrained equilibration simulation is run, which ensures that the membrane can adapt to the newly embedded protein. Then the real simulation starts.

Written by: Elmar Krieger & Kornel Ozvoldik

License: GNU GPL

Last modified: 2024/05/14

° 

Running a steered molecular dynamics simulation in water

This macro sets up and runs a steered simulation to pull a ligand from its binding site. As soon as they have been pulled apart completely and atoms start crossing periodic boundaries, the macro works no longer correctly. You can set the pulling 'acceleration' below, which is not changed during the simulation. So this macro can currently not tell you how strong you need to pull, it can only tell you whether the chosen acceleration is strong enough or not to pull the ligand off.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/10/16

° 

Play back a molecular dynamics trajectory

This macro provides an interactive molecular dynamics trajectory player

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/10/16

Download: md_play.mcr

° 

Analyzing a molecular dynamics trajectory in blocks

This macro analyzes a simulation in blocks. YASARA runs for each block the md_analyze macro and creates in the end a single report containing detailed information about the system, e.g. energies, RMSDs, hydrogen bonds per block. By default the simulation snapshots are split into 3 blocks, so the second and third block can be compared to see if convergence has been reached. The number of blocks can be freely choosen as well as setting the number of snapshots to be compiled into one block instead. For more information about the methods used and how to add your own analysis see the md_analyze macro.

Written by: Elmar Krieger and Kornel Ozvoldik

License: GNU GPL

Last modified: 2023/10/16

° 

Refining a homology model by molecular dynamics simulation in water

This macro runs a 500 ps simulation of a homology model using the protocol described in Proteins 57,678-683. It saves one PDB file every 25 picoseconds, including a table with force field energies to identify the best snapshot. In YASARA Structure, this table additionally contains Dihedrals, Packing1D and Packing3D checks (see the Check command). And in the Twinset, three more checks are included: PhiPsi, Backbone and Packing3. Note that MD simulations tend to make models worse if you simulate too long.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

Download: md_refine.mcr

° 

Convert between Sim, XTC, MDCrd and PDB simulation trajectories

This macro converts an existing MD trajectory between various formats. Supported are conversions between YASARA Sim, GROMACS XTC and AMBER MDCrd trajectories, as well as conversion to PDB files

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

Download: md_convert.mcr

° 

Analyzing the change of secondary structure in a molecular dynamics trajectory

This macro analyzes a simulation and creates a table with per-residue secondary structure as well as a graphical plot (save a ray-traced screenshot with 'Ball zoom' set to 2).

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

° 

Analyzing a molecular dynamics trajectory per residue

This macro analyzes a simulation and creates a table with average RMSFs and RMSDs from the starting structure per residue

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

° 

Analyzing a molecular dynamics trajectory of multiple objects

This macro analyzes a simulation and creates a table with energies and RMSDs from the starting structures for multiple objects

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

° 

Running a molecular dynamics simulation of a membrane protein with maximum speed

This macro sets up and runs a simulation of a membrane protein as quickly as possible. It scans the protein for secondary structure elements with hydrophobic surface residues, orients it accordingly and embeds it in a membrane of adjustable lipid composition. Finally a 250 ps restrained equilibration simulation is run, which ensures that the membrane can adapt to the newly embedded protein. Then the real simulation starts.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2020/11/10

° 

Running an accurate molecular dynamics simulation in water with maximum speed

This macro runs a simulation as quickly as possible, using the standard macro md_run.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2015/03/24

Download: md_runfast.mcr

° 

Molecular modeling

° 

Packing library

This macro provides a function library for assembling giant structures with millions to billions of atoms, to be 'included' by other macros

Written by: Elmar Krieger & Kornel Ozvoldik

License: GNU GPL

Last modified: 2023/10/16

° 

Running an energy minimization

This macro runs an energy minimization of the soup without cleaning it first. Compared to the normal minimization experiment, it temporarily adds a water shell so that all force fields can be used directly, also those optimized for use with explicit solvent

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2019/11/22

Download: em_run.mcr

° 

Running an energy minimization

This macro cleans the soup and then runs an energy minimization. Compared to the normal minimization experiment, it temporarily adds a water shell so that all force fields can be used directly, also those optimized for use with explicit solvent

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2019/07/06

Download: em_runclean.mcr

° 

Structure determination

° 

Analyze and merge an NMR ensemble

This macro analyzes an NMR ensemble and merges the members into a single PDB file

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2022/10/17

Download: nmr_analyze.mcr

° 

Set default parameters for NMR structure determination

This macro is included by others to set default parameters

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2022/09/20

° 

Refine an NMR ensemble in explicit solvent

This macro refines an NMR ensemble in explicit solvent (HOH or DMSO) to increase the accuracy

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2022/01/22

° 

Solve an NMR structure with fixed quality

This macro solves an NMR structure based on a protein sequence and a restraint file in XPLOR format. Contrary to the normal nmr_solve macro, it refines only those structures in water, that match a quality criterion, and performs an unlimited number of folding runs until the ensemble has the requested size

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2019/09/15

° 

Solve an NMR structure

This macro solves an NMR structure based on a protein sequence and a restraint file in XPLOR format

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2019/09/15

Download: nmr_solve.mcr

° 

Refine an NMR ensemble in vacuo

This macro refines a roughly folded structure in vacuo to create a realistic fold

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2018/02/05

° 

Fold an ensemble of NMR structures

This macro folds a structure from the stretched-out conformation using NMR restraints

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2018/02/05

Download: nmr_fold.mcr

° 

Solve an NMR structure

This macro refines a single structure or structure ensemble (provided in PDB format) in vacuo & explicit solvent and finally analyzes the result

Written by: Elmar Krieger, Sander Nabuurs, Chris Spronk

License: GNU GPL

Last modified: 2016/02/27

Download: nmr_refine.mcr

° 

Structure prediction

° 

Docking a ligand to a receptor locally

This macro runs VINA or AutoDock to sample local ligand conformations in a user-supplied complex. An analysis log file is written at the end.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/09/18

° 

Docking a ligand to a receptor ensemble with flexible side-chains

This macro runs VINA or AutoDock to predict the structure of a ligand-receptor complex. Receptor flexibility is considered by creating a receptor ensemble with alternative high-scoring solutions of the side-chain rotamer network. An analysis log file is written at the end. The macro can also continue a docking run that got interrupted, especially if the scene has not been moved or rotated manually during the docking.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/09/18

° 

Docking a ligand to a receptor

This macro runs VINA or AutoDock to predict the structure of a ligand-receptor complex. It can also continue a docking run that got interrupted. An analysis log file is written at the end.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/09/18

Download: dock_run.mcr

° 

Docking a covalently bound ligand to a receptor

This macro runs VINA or AutoDock to predict the structure of a covalent ligand-receptor complex. It can also continue a docking run that got interrupted. An analysis log file is written at the end.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/09/18

° 

Docking many ligands to a receptor to perform virtual screening

This macro runs VINA or AutoDock to dock multiple ligands against a receptor and saves a sorted table of their binding energies, as well as the corresponding complexes

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/07/30

° 

Building a homology model quickly

This macro builds a homology model using a FASTA sequence of the target, and optionally template structures and alignments. It takes some shortcuts (alignment quality, number of templates, refinement) to finish as quickly as possible

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/04/10

° 

Building a protein model

This macro builds a protein model using a FASTA sequence of the target, and optionally template structures and alignments. AI-based ab initio structure predictions are made locally using ESMFold and OmegaFold (7 GB of disk space needed). If 'foldmethod' is set to 'All', AlphaFold predictions will also be used (75 GB of disk space needed).

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/04/10

Download: pm_build.mcr

° 

Building a homology model

This macro builds a homology model using a FASTA sequence of the target, and optionally template structures and alignments. To include AI-based ab initio structure prediction, use the 'pm_build' macro.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2024/04/10

Download: hm_build.mcr

° 

Rescore ligand poses with an energy minimization and local docking

This macro loads a receptor-ligand complex obtained with YASARA's own docking macros or other software, energy-minimizes each ligand pose and rescores it using AutoDock VINA's local search, confined closely to the original ligand pose.

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

° 

Showing the docking results interactively

This macro provides an interactive player to cycle through the docking poses and clusters

Written by: Elmar Krieger

License: GNU GPL

Last modified: 2023/05/08

Download: dock_play.mcr