Edge Mapper for RBFE calculations [MDRun] [FECalc]

  • Purpose:

    • This Floe produces ligand pairs or edges to run RBFE calculations

  • Method Recommendations/Requirements:

    • One inputs is required:

      • A ligand dataset. The provided ligands must have reasonable 3D coordinates, all atoms, and correct chemistry such as bond orders and formal charges should be correctly assigned

      • Optionally, the user can also import into the mapper floe a map generated externally, to generate a mapper output dataset of only those edges defined by the external map. The syntax of this external map file is described below

  • Limitations

    • The provided ligand set should contain similar ligands where a common scaffold is used and functional groups are added to it. Large scaffold mutations should be avoided to improve RBFE and ABFE result errors

  • Expertise Level:

    • Regular

  • Compute Resource:

    • low

  • Keywords:

    • MD, FECalc

  • Related Floes: * Non-Equilibrium Switching Recovery [Utility] [FECalc]

In general given N ligands N(N-1)/2 edges can be used to run RBFE calculations. The mapper floe selects a reasonable subset of these edges with the goal to keep the computational demand low. The edges are selected based on set of heuristics like ligand mcs similarities, equal charges and OpenEye rocs score. The produced floe report shows the similarity scores as numbers in the range [0, 1] where zero means very unlikely and 1 means very likely to succeed. The floe mapper output can be used with the NES floe to run RBFE calculations. In order to have successful RBFE results the ligands should be similar and do not attempt large mutations. The Mapper floe can accept a edge text input file as well. In this case the produced map is defined in the provide text file and no attempt to use internal heuristics is done. In addition, in this case the edge scores are not calculated and on each edge the string “N/A” will be displayed in the floe report.The text file of edges, describes the map of desired alchemical transformations of one ligand into another; each transformation forms an edges of a connected graph of ligands. The file must have one line per transformation, of format

ligA_name >> ligB_name

where “ligA_name” and “ligB_name” are the ligand title names for the ligands to be transformed. These ligand names must correspond exactly to those in the provided ligand set

Promoted Parameters

  • mapper (dataset_out) : Mapper Output text file
  • edges (file_in) : The ligand edge mapping file
  • ligands (data_source) : Ligand or Bound Equilibrium Reader

Extra Required Parameters

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