Non-Equilibrium Switching [MDRun] [FECalc]

  • Purpose:

    • This Floe performs Relative Binding Free Energy (RBFE) calculations using the Non-Equilibrium Switching (NES) method refined by the de Groot lab (Gapsys et al., Chem. Sci., 2020, 11, 1140-1152).

  • Method Recommendations/Requirements:

    • Four inputs are required:

      • An Orion dataset containing an equilibrium run for each bound ligand participating in the NES run.

      • An Orion dataset containing an equilibrium run for each unbound ligand participating in the NES run.

      • A Orion Dataset containing the ligand edges to run the RBFE calculation generated by using the Mapper floe

      • [Optional] a text file containing experimental binding free energies for at least one ligand, one experimental datapoint per line, of form “ligA_name {deltaG(exptl)} {error_deltaG(exptl)} {units}” for example, “gn1c -8.56 0.17 kcal/mol”.

  • Limitations

    • If no experimental binding free energies (the fourth input above) are given, the estimation of ligand binding free energies has no reference value so the relative values will be centered around the mean.

    • Currently there is no mitigation for the effects of changes in buried waters, protein sidechain flips, or large protein movements between ligA and ligB.

  • Expertise Level:

    • Regular

  • Compute Resource:

    • High

  • Keywords:

    • MD, FECalc

  • Related Floes:

    • Ligand Bound and Unbound Equilibration for NES [MDPrep] [MD]

    • Equilibration and Non-Equilibrium Switching [MDPrep] [MD] [FECalc]

    • Compare Experimental Affinity with NES Results [Utility] [FECalc]

    • Non-Equilibrium Switching Recovery [Utility] [FECalc]

The floe will draw a number of starting snapshots from the bound and unbound trajectories of the ligands. Then for each edge in the edge file, it will generate an RBFE alchemical transformation from ligA into ligB, and carry out the NES fast transformation of ligA into ligB, and vice versa, for each of the snapshots. The resulting relative free energy change, or DeltaDeltaG, for each edge is the primary output of this method. A maximum likelihood estimator is then used to derive a predicted binding affinity (free energy, or DeltaG) for each ligand. The mean value of the input experimental binding free energies is used as the reference value for the computed ones.

The speed of the NES transformation and the number of snapshots transformed can be adjusted from default values by the user at runtime. The floe outputs two floe report/dataset pairs, one for the calculated RBFE edges (DeltaDeltaGs), and one for the derived affinity predictions (DeltaGs) of ligand.

Promoted Parameters

  • fail (dataset_out) : NES Dataset Failures out
  • recovery (dataset_out) : Recovery Out
  • exp (file_in) : The ligand affinity experimental file with affinities in units of kcal/mol or kJ/mol
  • nes_time (decimal) : NPT simulation time in nanoseconds
    Default: 0.05
  • out (dataset_out) : Output dataset to write to
  • DG (dataset_out) : Affinity Out
  • max_mapper_edges (integer) : The max number of mapper edges allowed
    Default: 100
  • cpu_count_md (integer) : The number of CPUs to run this cube with
    Default: 16 Min: 1 Max: 128
  • gpu_count_md (integer) : The number of GPUs to run this cube with
    Default: 1 Max: 16
  • spot_policy_md (string) : Control cube placement on spot market instances
    Default: Required
    Choices: Allowed, Preferred, NotPreferred, Prohibited, Required
  • bound (data_source) : Bound Input Dataset
  • map (data_source) : Mapper Input Dataset
  • trajectory_frames (integer) : The total number of trajectory frames to be used along the NE switching
    Default: 80
  • unbound (data_source) : Unbound Input Dataset

Extra Required Parameters

  • Log Field (Field Type: String) : The field to store messages to floe report
    Default: Log Field
  • Log Field (Field Type: String) : The field to store messages to floe report
    Default: Log Field
  • Log Field (Field Type: String) : The field to store messages to floe report
    Default: Log Field
  • Log Field (Field Type: String) : The field to store messages to floe report
    Default: Log Field
  • Log Field (Field Type: String) : The field to store messages to floe report
    Default: Log Field