Ligand Bound and Unbound Equilibration for NES [MDPrep] [MDRun]¶
Purpose:
This Floe performs MD simulations given a prepared protein and a set of posed and prepared ligands, running both bound and unbound simulations of each ligand in preparation for Relative Binding Free Energy Calculations using Non-Equilibrium Switching.
Method Recommendations/Requirements:
The ligands need to have reasonable 3D coordinates, all atoms, and correct chemistry (in particular bond orders and formal charges).
Each ligand can have multiple conformers but each conformer will be run separately as a different ligand.
The starting poses should not have very high gradients, in particular no bad clashes with the protein.
The protein needs to be prepared to MD standards: protein chains must be capped, all atoms in protein residues (including hydrogens) must be present, and missing protein loops resolved or capped.
Crystallographic internal waters should be retained where possible.
Limitations
Currently this floe cannot handle covalent bonds between different components such as ligand, protein, and cofactors.
Glycosylation on proteins is truncated and the amino acid is capped with H.
Expertise Level:
Regular/Intermediate/Advanced
Compute Resource:
Depends on simulation length; Minimal resources for default 6 ns.
Keywords:
MD, MDPrep
Related Floes:
Bound Protein-Ligand MD [MDPrep] [MD]
Short Trajectory MD with Analysis [MDPrep] [MD]
Given the inputs of the protein and posed ligands, the complex is formed with each ligand/conformer separately, and the bound and unbound simulations are then carried out. Currently only one of the conformers will be used in the NES calculations. A minimization stage is performed on the system followed by a warm up (NVT ensemble) and several equilibration stages (NPT ensemble). In the minimization, warm up, and equilibration stages, positional harmonic restraints are applied on the ligand and protein. At the end of the equilibration stages a production run (by default 6 ns) is performed on the unrestrained system. Two datasets are written, one for the bound and one for the unbound ligands.
Promoted Parameters
cpu_count_md (integer) : The number of CPUs to run this cube withDefault: 16 Min: 1 Max: 128 gpu_count_md (integer) : The number of GPUs to run this cube withDefault: 1 Max: 16 spot_policy_md (string) : Control cube placement on spot market instancesDefault: PreferredChoices: Allowed, Preferred, NotPreferred, Prohibited, Required md_engine (string) : Select the MD available engineDefault: OpenMMChoices: OpenMM, Gromacs prod_us_ns (decimal) : The Equilibration production running time in nsDefault: 6.0 hmr_us (boolean) : On enables Hydrogen Mass Repartitioning. Not currently implemented in GromacsDefault: True fail (dataset_out) : Dataset Failures out ligands (data_source) : Ligand Dataset charge_ligands (boolean) : Assign Ligand Partial Charges or notDefault: True flask_title (string) : Prefix name used to identity the Protein. If not specified, it will use the title of the input protein.Default: “” out_bound (dataset_out) : Bound Dataset Out protein_ff (string) : Force field parameters to be applied to the proteinDefault: Amber14SBChoices: Amber14SB, Amber99SB, Amber99SBildn, AmberFB15 ligand_ff (string) : Force field to be applied to the ligandDefault: OpenFF_2.0.0Choices: Gaff_1.81, Gaff_2.11, OpenFF_1.1.1, OpenFF_1.2.1, OpenFF_1.3.1, OpenFF_2.0.0, Smirnoff99Frosst protein (data_source) : Protein Dataset out_unbound (dataset_out) : Unbound Dataset Out
Extra Required Parameters
Log Field (Field Type: String) : The field to store messages to floe reportDefault: Log Field Log Field (Field Type: String) : The field to store messages to floe reportDefault: Log Field Log Field (Field Type: String) : The field to store messages to floe reportDefault: Log Field Log Field (Field Type: String) : The field to store messages to floe reportDefault: Log Field