Structure-Based Virtual Screening Using Phase
Tutorial Created with Software Release: 2025-3
Topics: Hit Discovery , Small Molecule Drug Discovery , Virtual Screening
Products Used: Phase
|
69.2 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 displayed
Abstract:
This tutorial demonstrates the creation, validation, and application of pharmacophore hypotheses to recognize common protein-ligand interactions and use them in virtual screening. You will learn how to create a pharmacophore hypothesis using a protein-ligand complex, how to modify a pharmacophore hypothesis to bias by experimental observables and to screen against the hypothesis to identify Leukotriene-A4 hydrolase inhibitors.
Tutorial Content
1. Pharmacophore Modeling Prerequisites
Structure files obtained from the PDB, vendors, and other sources often lack necessary information for performing modeling-related tasks. Typically, these files are missing hydrogens, partial charges, side chains, and/or whole loop regions. In order to make these structures suitable for modeling tasks, we use the Protein Preparation Workflow to resolve issues. Similarly, ligand files can be sourced from numerous places, such as vendors or databases, often in the form of 1D or 2D structures with unstandardized chemistry. LigPrep can convert ligand files to 3D structures, with the chemistry properly standardized and extrapolated, ready for use in virtual screening.
In this tutorial, the protein and ligands have already been prepared in order to save time. Please see the Introduction to Structure Preparation and Visualization tutorial for instructions on using the Protein Preparation Workflow and LigPrep.
2. 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 automatically added to the Entriesa 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 Entriesa 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 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/sbvs_phase.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.
- Choose lkha4_pharmacophore.prjzip.
- Click Open.
- In the Save scratch project dialog box, click OK.
- Go to File > Save Project As.
- Change the File name to Phase_tutorial.
- Click Save.
- The project is now named
Phase_tutorial.prj. - Structures are shown in the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
- The project is now named
Note: Please see the Glossary of Terms for the difference between includedthe entry is represented in the Workspace, the circle in the In column is blue and 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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries.
3. Creating, Viewing, and Manipulating Hypotheses
In this section, we will create a pharmacophore hypothesis from a protein-ligand complex. We will learn how to merge hypotheses, add a site to a hypothesis that wasn’t recognized by the feature definitions, as well as generate and manipulate excluded volumes.
3.1 Create a pharmacophore hypothesis from a protein-ligand complex
- Confirm that the 3CHP_prepared 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.
- Double-click Presets.
- The structure is rendered in the default custom preset and the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed zooms to the ligand.
- Go to Tasks > Browse > Ligand-Based Virtual Screening > Develop Pharmacophore Hypothesis.
- The Develop Pharmacophore Model panel opens.
- Under Create pharmacophore model using, choose Receptor-ligand complex (Workspace).
- For Method, choose Auto (E-Pharmacophore).
- Under Choose ligand, click A:4BO (901).
- Click the Hypothesis Settings button.
- The Hypothesis Settings dialog box opens.
- In the Features tab, choose Donors as vectors.
- Go to the Excluded Volumes tab.
- Check Create receptor-based excluded volumes shell.
- Click Save.
- In the Develop Pharmacophore Model panel, change the Job name to lkha4_epharm_exvol.
- Click Run.
- This job takes a few minutes.
- A banner appears in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed when the job has been incorporatedonce a job is finished, output files from the working directory are added to the project and shown in the Entry List and Project Table.
- Ctrl+Click (Cmd+Click) to includethe entry is represented in the Workspace, the circle in the In column is blue lkha4_epharm_exvol and 3CHP_prepared in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- The generated E-Pharmacophore is now overlaid in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed with the input protein-ligand complex.
3.2 Create a pharmacophore hypothesis manually
Several studies have noted the presence of hydrogen-bonding between the TRP 311 residue backbone and the bound ligand. Since that feature is not present in the previously generated pharmacophore hypothesis, you can add it through a manual hypothesis.
- In the Develop Pharmacophore Model panel, choose Manual.
Note: There are some features that were not included in the E-Pharmacophore that can be chosen when generating a Manual Hypothesis. This is due to a difference in cutoff between the two tools.
- Under Show features, drag the slider right to More.
- Features D3, N4, P5, and R6 are shown.
- Click Hypothesis Settings.
- The Hypothesis Settings dialog box opens.
- Check Create receptor-based excluded volumes shell.
- Click Save.
Note: Selecting Show inverse features from receptor generates additional features based on idealized geometries of specific interactions and hydrophobic features of the binding site.
- In the Workspace Configuration Toolbar (bottom right), right-click Interactions Toggle.
- The Interactions panel opens.
- Next to Non-covalent bonds, choose Ligand-receptor.
- Toggle on Aromatic H-bonds.
- Aromatic H-bonds are shown in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Under Contacts/Clashes, toggle off Bad.
- Bad clashes are hidden in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- In the Hierarchy, search for TRP 311 using the search toggle.
- Select(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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries TRP 311.
- The TRP 311 residue is selected in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Click Fit view to selected atoms.
- The Workspacethe 3D display area in the center of the main window, where molecular structures are displayed zooms to TRP 311.
Note: You may also double-click TRP 311 in the Hierarchy to zoom to the selected residue.
- In the Develop Pharmacophore Model panel, under Choose features, check Add and choose Donor in Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- A banner appears prompting you to pick an atom.
- Select(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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries atom H20 on the ligand, which makes an aromatic H-bond with the backbone of TRP 311.
- A donor vector D7 appears in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- D7 is checked in the Choose features box in the Develop Pharmacophore Model panel.
- In the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed, right-click the D7 vector.
- The Edit Feature D7 dialog box opens.
- Check Manually reposition feature.
- Click and drag the scroll wheel (Shift-click and drag on a Trackpad) to reposition the vector along the hydrogen bond.
- Click OK.
Note: See the Custom Mouse Actions Panel help topic for further information on mouse actions.
- In the Develop Pharmacophore Model panel, under Choose features, check all the features (D3, D7, N4, P5, R6).
- Change the Job name to lkha4_manual.
- Click Create.
- This job takes a few seconds.
- A banner appears in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed confirming the job has been successfully imported.
- A new entry titled lkha4_manual is added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
- Close the Develop Pharmacophore Model panel.
- Ctrl+Click (Cmd+Click) to includethe entry is represented in the Workspace, the circle in the In column is blue lkha4_manual and 3CHP_prepared in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- The generated E-Pharmacophore is now overlaid in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed with the input protein-ligand complex.
3.3 Manipulate a hypothesis
Note: For ligand screening, it is at times desirable to broaden or narrow the types of moieties a particular pharmacophore feature is allowed to match.
- Includethe entry is represented in the Workspace, the circle in the In column is blue lkha4_manual in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Next to lkha4_manual, click the H.
- Choose Show > Feature Properties.
- Properties of each feature are displayed in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
Note: Excluded volumes can be edited manually. As a pedagogical practice, we will demonstrate adding and removing excluded volumes in the next steps.
- Ctrl-click (Cmd-click) to includethe entry is represented in the Workspace, the circle in the In column is blue lkha4_manual and 3CHP_prepared in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- In the Hierarchy, search for HIP 139 using the Search toggle.
- Select(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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries HIP 139.
- Click Fit view to selected atoms.
- The Workspacethe 3D display area in the center of the main window, where molecular structures are displayed zooms to HIP 139 residue.
- Next to lkha4_manual, click the H.
- Choose Show > Manage excluded Volumes.
- Manage Excluded Volumes Dialog box opens.
- Click Selected Atoms.
- Excluded volumes are added for atoms in HIP 139.
- Information is added to the bottom of the table.
- Scroll to the bottom of the table and select(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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries the last excluded volume.
- The 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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries excluded volume is highlighted in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Click Delete.
- The highlighted excluded volume is deleted from the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Click Close.
3.4 Merge hypotheses
Merging hypotheses allows you to develop multiple models and choose common/desired features manually.
- Ctrl+Click (Cmd+Click) to includethe entry is represented in the Workspace, the circle in the In column is blue lkha4_epharm_exvol and lkha4_manual in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- For each includedthe entry is represented in the Workspace, the circle in the In column is blue entry, click the H and choose Show then Reference ligand
- Reference ligands are added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and are includedthe entry is represented in the Workspace, the circle in the In column is blue and fixed in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Go to Tasks > Browse > Ligand-Based Virtual Screening > Develop Pharmacophore Hypothesis.
- The Develop Pharmacophore Model panel opens.
- For Create pharmacophore model using, choose Merged hypotheses (Workspace).
- For Choose features, from lkha4_epharm_exvol, choose R6 and R7.
- For Choose features, from lkha4_manual, choose D3, N5, P6, and XVols.
Note: The D3 feature is the donor that we added manually. The R6 feature in the E-Pharmacophore aligns well with the R7 from the manual pharmacophore, so only one of them needs to be added to the merged hypothesis.
- Change the Job name to lkha4_merge.
- Click Create.
- This job takes a few seconds.
- A banner appears in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed confirming the job has been successfully imported.
- A new entry titled lkha4_merge is added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
- Close the Develop Pharmacophore Model panel.
4. Screening Ligands Against a Hypothesis
In this section, you will cover screening ligands against a hypothesis. You will then visualize the results and examine the PhaseScreenScore to evaluate how well the ligands align to the pharmacophore features of the hypothesis it was screened against.
4.1 Screen ligands from a file
- In the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, select(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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries the DUD-E_Actives and DUDE_Inactives groups.
- Includethe entry is represented in the Workspace, the circle in the In column is blue lkha4_merge in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
DUD-E is a ‘Database of Useful Decoys: Enhanced.’ For each active in the database, there are 50 decoys with similar 1-D physio-chemical properties, but different 2-D topology. The DUD-E ligand set can thus be used in pharmacophore screening to avoid artificial enrichment.
- Go to Tasks > Browse > Ligand-Based Virtual Screening > Ligand Screening.
- The Phase Ligand Screening panel opens.
- For Ligands to screen, choose Project Table (selected entries).
- 215 ligands will be screened.
- Click Add Hypothesis.
- Choose Workspace.
- lhka4_merge is added to the table.
- For the lkha_merge hypothesis, click the Matches column and choose 3 of 5.
- Ligands must satisfy at least 3 features of each hypothesis.
- Click the Screening Settings button.
- The Screening Settings dialog box opens.
- Go to the Groups tab.
- Check Pick features.
- A banner appears in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed prompting you to pick features to be grouped.
- Select features D1 and R5 in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Click Create Group.
- The table is populated.
Note: Any ligand that does not match at least one of the grouped features will not be returned in the screen, even if it matches the desired 3 of 5 total features.
- Click Save.
- Change the Job name to lkha4_screen.
- Click Run.
- This job takes ~1 minute.
- A banner appears in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed confirming the job has been incorporated.
- A new group lkha4_screen is added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
- Close the Phase Ligand Screening panel.
4.2 Visualize pharmacophore screening results
- Select(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 Entries (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries the lkha4_screen group by clicking on the group heading.
- Double-click the In circle of lhka4_merge to fix the hypothesis in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Includethe entry is represented in the Workspace, the circle in the In column is blue the top ligand.
- Click Style selected entries (paint brush icon) and choose ball-and-stick representation.
- Use the right and left arrow keys to step through the ligands.
- At the top right corner of the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, click the Change table settings (three vertical dots) and choose Show Property.
- The Show Properties in Table dialog box opens.
- Click Choose.
- Search for and select PhaseScreenScore and Matched Ligand Sites from the list.
- Click OK.
- The properties are added to the table.
- Hover over the Matched Ligand Sites circles to see what features are present.
Note: Only 46 of the 215 screened ligands were returned following the Phase Ligand Screen due to the matching criteria. Of those returned, only 6 were inactives, and those had low Phase Screen Scores relative to the rest of the dataset.
Please note that your results may vary slightly if you are using a different version of the software.
Both the quantity and quality of feature matching is taken into account in the PhaseScreenScore. A ligand can thus both match fewer sites and have a higher PhaseScreenScore than another ligand in the screening set.
5. Conclusion and References
In this tutorial, you learned how to create a pharmacophore hypothesis from a receptor-ligand complex, manipulate a hypothesis by adding features in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed and merge hypotheses. Finally, you screened a ligand set against the merged hypothesis, and examined their PhaseScreenScores and Matched Ligand Sites.
For further learning:
- Introduction to Structure Preparation and Visualization
- Forming Protein-Ligand Interactions with the Ligand Designer
- Ligand-Based Virtual Screening Using Phase
- Introduction to Molecular Modeling in Drug Discovery Online Course
- Designing quality ligand libraries online course
- Target enablement, preparation, & validation online course
- Virtual screening with integrated physics & machine learning online course
For further reading:
- Phase User Manual
- A Database of Useful Decoys : Enhanced
- Pharmacophore Screening with Phase - Schrӧdinger webpage with an overview of our pharmacophore screening technologies.
- Schrödinger has partnered with Enamine, MilliporeSigma, and MolPort to provide a Phase database of fragments, lead-like, near drug-like, and drug-like compounds. Options include: 1. Aldrich’s Market Direct Screening Collection 2. Molport’s Screening Compound Database 3. Enamine’s Stock Screening Compounds Collection.
6. Glossary of Terms
Entries - 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
incorporated - once a job is finished, output files from the working directory are added to the project and shown in the Entry List and Project Table
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 Entries (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