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
|
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
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
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.
-
Double-click the Maestro icon
- (No icon? See Starting Maestro)
- Go to File > Change Working Directory
- Find your directory, and click Choose
- 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
- After downloading the zip file, unzip the contents in your Working Directorythe location that files are saved for ease of access throughout the tutorial
-
Go to File > Open Project and choose
Ligand_Designer_WM.prjzip -
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
- Go to File > Save Project As
-
Change the File name to Ligand_Designer_WM, click Save
-
The project is now named
Ligand_Designer_WM.prj
-
The project is now named
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).
- 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
-
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
-
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
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.
-
Under Display, select Displaceable Waters
- These are high energy waters, are previously determined by WaterMap
- 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
-
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
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.
-
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
-
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.
Note: You may need to expand the table to see all of the columns
- Change the Threshold for PSA to 120 by double clicking
-
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
-
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.
- Continue going through the poses with the right arrow in the Ligand Designer
- 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
- Under Post-Processing, click the dropdown and select Sort by MPO
Note: You can also sort by Clashes and Good Interactions
-
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
- 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
- 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
- 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.
-
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
- 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
- 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
- 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
- Click View similar purchasable compounds
Note: If you do not have an FPsim-GPU license you will not be able to access this panel
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
-
Click Predict Poses
- The poses for the 10 most similar compounds are predicted and added to the Entry List
3.5 Export favorite compounds
- Click the star icon to add the top purchasable compound to favorites
- Under Post-Processing, click the dropdown and choose Save to File
- Click the dropdown that says All, and choose Favorites
-
Click Save
- The Save File panel opens
- For File name, type ND-022_036_enum_favorite
-
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.
For further learning:
- Enumeration Tools for Library Design
- Identifying Binding Site Requirements and Lead Optimization with WaterMap
- Forming Protein-Ligand Interactions with the Ligand Designer
- Introduction to Molecular Modeling in Drug Discovery Online Course
- High-Throughput Virtual Screening for Hit Finding and Evaluation Online Course
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