WaterMap-Guided Lead Optimization with the Ligand Designer

Tutorial Created with Software Release: 2023-3
Topics: Hit-to-Lead & Lead Optimization, Medicinal Chemistry Design, Small Molecule Drug Discovery
Products Used: Ligand Designer

Tutorial files

0.4 MB

This tutorial is written for use with a 3-button mouse with a scroll wheel.
Words found in the Glossary of Terms are shown like this: Workspacethe 3D display area in the center of the main window, where molecular structures are displayedthe 3D display area in the center of the main window, where molecular structures are displayed

 

Tip: You can hover over a glossary term to display its definition. You can click on an image to expand it in the page.
Abstract:

 

In this tutorial, you will learn how to use the Ligand Designer to optimize a hit molecule.

 

Tutorial Content
  1. Creating Projects and Importing Structures

  1. Displacing Unstable Waters with the Ligand Designer

  1. Analyzing Enumeration Outputs with the Ligand Designer

  1. Conclusion and References

  1. Glossary of Terms

1. Creating Projects and Importing Structures

At the start of the session, change the file path to your chosen Working Directorythe location that files are saved in Maestro to make file navigation easier. Each session in Maestro begins with a default Scratch Projecta temporary project in which work is not saved, closing a scratch project removes all current work and begins a new scratch project, which is not saved. A Maestro project stores all your data and has a .prj extension. A project may contain numerous entries corresponding to imported structures, as well as the output of modeling-related tasks. Once a project is created, the project is automatically saved each time a change is made.

Structures can be imported from the PDB directly, or from your Working Directorythe location that files are saved using File > Import Structures, and are added to the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and Project Tabledisplays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data. The Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion is located to the left of the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed. The Project Tabledisplays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data can be accessed by Ctrl+T (Cmd+T) or Window > Project Table if you would like to see an expanded view of your project data.

  1. Double-click the Maestro icon

Figure 1-1. Change Working Directory option.

  1. Go to File > Change Working Directory
  2. Find your directory, and click Choose
  3. Pre-generated input and results files are included for running jobs or examining output. Download the zip file here: https://www.schrodinger.com/sites/default/files/s3/release/current/Tutorials/zip/ligand_designer_wm.zip
  4. After downloading the zip file, unzip the contents in your Working Directorythe location that files are saved for ease of access throughout the tutorial  

Figure 1-2. Open Project panel.

  1. Go to File > Open Project and choose Ligand_Designer_WM.prjzip
  2. Click Open
    • Structures are shown in the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
    • A structure is includedthe entry is represented in the Workspace, the circle in the In column is blue in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Figure 1-3. Save Project panel.

  1. Go to File > Save Project As
  2. Change the File name to Ligand_Designer_WM, click Save
    • The project is now named Ligand_Designer_WM.prj

2. Displacing Unstable Waters with the Ligand Designer

In this section, we are starting with the hit compound ND-022, and using the results from the WaterMap analysis to further optimize the hit by displacing a high energy water in the binding site (see https://www.pnas.org/content/113/13/E1796 for more information).

Figure 2-1. Include entries in the Workspace.

  1. Use Ctrl+click to include ACC Receptor, ND-022, and watermap_ND-022 in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Note: Recent advances in the Ligand Designer now allow two receptors to be used at once

Figure 2-2. Ligand Designer in Lead Optimization.

  1. Go to Tasks > Browse > Lead Optimization > Ligand Designer
    • The Ligand Designer panel opens

The Ligand Designer panel does not currently support covalently-bound ligands. For covalent-bound ligand enumerations and docking, see Covalent Docking for Virtual Screening and Pose Prediction

Figure 2-3. Analyze Workspace.

  1. Click Analyze Workspace
    • New entries are added to the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion in the Ligand Designer group
    • The Workspacethe 3D display area in the center of the main window, where molecular structures are displayed now has just a single ligand and a cloud surrounding it to represent the growth space

Note: Uncheck Adjust view when analyzing to retain the previous view in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Note: The darker blue growth space indicates the region is solvent exposes

 

Figure 2-4. Select a Primary Receptor.

Note: You can also design ligands with two receptors. One is the primary receptor, for which you want to optimize binding. The other can be a receptor for which you want to minimize off-target effects, for example; or a conformer of the primary receptor, for which you want to optimize binding for both receptors.

 

The example shown in Figure 2-4 shows two identical receptors. In this case, the choice of which is primary is irrelevant.

Figure 2-5. Visualize Displaceable waters in the Workspace.

  1. Under Display, select Displaceable Waters
    • These are high energy waters, are previously determined by WaterMap
  2. Click the Workflows dropdown, and select Displace Unstable Water

Note: The protein structure has been undisplayed using the protein quick-select and the style panel

Figure 2-6. Select high energy water in the Workspace.

  1. Select the marked high energy water by the phenyl-group of ND-022
    • Vectors appear from different possible attachment points
    • A pop-up explaining your options appears

Figure 2-7. Select attachment point.

  1. Select the marked attachment point
    • A pop-up up appears

 

By clicking the filter icon next to Enumerate, you can filter out your outputs by various properties.

 

By clicking the cog icon you can customize which library is used in the enumeration. By default, the Ligand Designer chooses the library that was curated for the particular type of enumeration you have selected. If you change the library and click save in the panel, it will keep that library as the default (and the cog wheel will be highlighted blue to indicate that a custom library is being used). See the final section of the Enumeration Tools for Library Design tutorial for more information on how to create custom R-group libraries.

Figure 2-8. Launch the enumeration.

  1. Click Enumerate
    • The enumeration job is launched
    • New entries are added to the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion after they are successfully docked
    • You can follow the status of the job with the status bar at the bottom of the Ligand Designer panel

3. Analyzing Enumeration Outputs with the Ligand Designer

All 71 outputs from the enumeration have been added to a new group in the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion. We will now use various features of the Ligand Designer to evaluate and inspect the poses of our compounds and decide which we would like to share with colleagues.

3.1 Multi-parameter optimization

Figure 3-1. Radar plot.

  1. Next to Multi-Parameter Optimization, click the cog
    • A pop-up appears

Note: You can customize the relative importance of various parameters, the expected value distribution, and the thresholds across that distribution.

Figure 3-2. Customize MPO.

Note: You may need to expand the table to see all of the columns

 

  1. Change the Threshold for PSA to 120 by double clicking
  2. Click OK
    • The Scoring function is changed to “Custom”
    • The radar plot in the Ligand Designer is now populated based on the updated MPO
    • The MPO value is a geometric mean of the score assigned to each of the parameters (weighted by the specified weights)

Click the dropdown that says Lipinski to choose from other pre-generated MPOs. Also, if you click Edit Properties you can add additional properties to your MPO.

3.2 Pose inspection

Figure 3-3. Analyze ND-022_1 in the Workspace.

  1. Click the right arrow in the Ligand Designer
    • ND-022_1 is now included in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed
    • The radar plot has been updated with the properties of the new compound

If desired, toggle on the protein structure in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed by checking ACC Receptor in the Structure Hierarchy. Toggle the Growth Space back on if it disappears.

Figure 3-4. Pose inspection and favoriting.

  1. Continue going through the poses with the right arrow in the Ligand Designer
  2. When you see a pose you like, click the star icon to add it to your favorites

You can review favorite compounds by clicking the dropdown next to the icon and selecting Group Favorites in Project, which will create a new group on the top of the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion with only the favorite compounds

Figure 3-5. Sort by MPO.

  1. Under Post-Processing, click the dropdown and select Sort by MPO

Note: You can also sort by Clashes and Good Interactions

  1. Click Sort
    • Compounds in the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion are now sorted by MPO score

3.3 Modify the compound with the 2D sketcher

Figure 3-6. Include ND-022_1.

  1. Include ND-022_36 (the enumerated compound with the best MPO score), in the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion

Figure 3-7. 2D Sketcher in the Ligand Designer.

  1. For Display, choose 2D/3D Editing

Note: You can make edits either directly in the sketcher or in the Workspace using the 3D Builder. Any change you make in the Workspace with the 3D builder will be mirrored in the 2D sketcher

  1. Make a change to ND-022_36. Here we have converted an alkane into an alkene. Use predicted properties while creating your analog. If you choose to use the 2DSketcher, make sure you are in Draw Mode indicated by the pencil icon.

Note: The properties under the sketcher are synced with the sketcher and will update automatically before your compound is docked.

Figure 3-8. Predict Pose.

  1. Click Predict Pose
    • The Pose is predicted and added to the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
  2. Includethe entry is represented in the Workspace, the circle in the In column is blue the new Entry in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed
  3. Click the X icon to close the 2D Sketcher

 

If you click the arrow on the Predict Pose icon you can switch from using the recommended MCS docking protocol to a minimization protocol.

3.4 View similar purchasable compounds

Figure 3-9. View similar purchasable compounds.

  1. Include the user-designed_ligand in the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
  2. Click View similar purchasable compounds

Note: If you do not have an FPsim-GPU license you will not be able to access this panel

 

Figure 3-10. Predict poses of similar purchasable compounds

Note: The Purchasable Compounds panel performs a fingerprint-based search against over a billion purchasable compounds and returns the 10 most similar to the compound that is included in the Workspace

  1. Click Predict Poses
    • The poses for the 10 most similar compounds are predicted and added to the Entry List

3.5 Export favorite compounds

Figure 3-11. Favorite the first purchasable compound.

  1. Click the star icon to add the top purchasable compound to favorites

Figure 3-12. Save the favorite poses to File.

  1. Under Post-Processing, click the dropdown and choose Save to File
  2. Click the dropdown that says All, and choose Favorites
  3. Click Save
    • The Save File panel opens

Figure 3-13. Save to your Working Directory.

  1. For File name, type ND-022_036_enum_favorite
  2. Click Save
    • The favorite compounds are saved to your Working Directorythe location that files are saved

4. Conclusion and References

In this tutorial, we first enumerated and docked 64 analogs of the hit molecule that would displace a key, high energy water. We then modified the MPO to thresholds that would work best for our system, inspected and selected favorite poses, and modified one of the compounds using the 2D sketcher. Finally, we saved the favorite subset of compounds to file so we could share them with our colleagues.

5. Glossary of Terms

Entry List - a simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion

included - the entry is represented in the Workspace, the circle in the In column is blue

Project Table - displays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data

Scratch Project - a temporary project in which work is not saved, closing a scratch project removes all current work and begins a new scratch project

selected - (1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries

Working Directory - the location that files are saved

Workspace - the 3D display area in the center of the main window, where molecular structures are displayed