Conformational Analysis for Small Molecules Using MacroModel and ConfGen
Tutorial Created with Software Release: 2025-3
Topics: Hit Discovery , Hit-to-Lead & Lead Optimization , Small Molecule Drug Discovery
Products Used: ConfGen , MacroModel
|
21.9 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:
In this tutorial, you will investigate tools that are useful for exploring conformational space for small molecules. MacroModel will be used to set up and analyze simple torsions of a single compound to gain a better understanding of the energetic barriers of a small molecule. ConfGen will then be used to screen molecules with diverse linker atoms as it is a fast method geared towards retrieving low energy bioactive conformations.
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 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 if you would like to see an expanded view of your project data.
-
Double click the Maestro icon to start Maestro.
- (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/macromodel_confgen.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.
-
Select the file
MM_ConfGen.prjzip. -
Click Open.
- Structures are shown in the Entriesa 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.
- In the Save scratch project dialog box, click OK.
Note: Please see the Glossary of Terms for the distinction 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 .
- Go to File > Save Project As.
- Change the File name to MM_Confgen_tutorial.
-
Click Save.
-
The project is named
MM_Confgen_tutorial.prj.
-
The project is named
2. Preparing to Explore Conformational Space
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, or whole loop regions. In order to make these structures suitable for modeling tasks, we use the Protein Preparation Workflow to resolve these 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 proteins and ligands have already been prepared in order to save time. However, these preparation steps are a necessary part of conformational analysis and must be done before using MacroModel or ConfGen. Please see the Introduction to Structure Preparation and Visualization tutorial for instructions on using the Protein Preparation Workflow and LigPrep.
Conformer generation is useful in many aspects of both molecular modeling in general and drug discovery in particular. Moreover, the ability to generate a bioactive conformer is a vital prerequisite to any successful computer-aided drug design project. Efficient conformer sets have wide-ranging ramifications in downstream applications. For example, with fewer irrelevant conformations to process, virtual database screens and shape-based similarity searches run to completion in a fraction of the time without sacrificing accuracy.
There are three tools for exploring the conformational space of a ligand within the Schrodinger Small Molecule Drug Discovery suite known as MacroModel, ConfGen, and Prime. Macromodel provides a wealth of conformational search options including those for small molecules, proteins, and protein-ligand complexes. ConfGen is embedded into a number of applications as it is suited to explore large scale conformer generations and is tailored to predict active conformations with successful ligand binding to a receptor. The third tool available is Prime, a protein centric tool used for large structures such as macrocycles. This tutorial will address MacroModel and ConfGen tools to explore the conformational space of small molecules. It is important to choose the most appropriate application for your research project. See the MacroModel Product Overview, ConfGen Product Overview,and Prime Product Overview for more information. For information about generating conformations for Glide Docking, see this KB Article .
3. Scanning Torsions Using MacroModel
In this section, we will use the MacroModel Coordinate Scan to set up and analyze simple torsions of a single compound to gain a better understanding of the energetic barriers of a small molecule. We will superimpose the output conformations to a common core, for easier visualization. Next, we will use the potential energy surface to find low energy conformations of a small molecule. These steps are a prelude to a full conformational search, which can explore many torsions of a small molecule simultaneously or systematically. We will use the potential energy surface to find low energy conformations of a small molecule. Investigating conformational space is important for understanding structural-property dependence, and for constructing initial models for Molecular Dynamics simulations.
3.1 Set up a torsional scan
-
Confirm the Ddrive entry is 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 in the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
- The structure is shown in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
Note: Atom number labels were added from the Style toolbox by clicking the Apply Labels arrow and choosing Element + Atom Number.
-
Go to Tasks > Browse > MacroModel > Coordinate Scan.
- The Coordinate Scan panel opens.
- For Use structures from, choose Workspace (included entry).
- Go to the Scan tab.
-
For Coordinate Type, choose Dihedral.
- A banner appears prompting you to pick four atoms to add dihedral coordinate.
-
In the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed, 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 numbers 1, 14, 15, and 16.
- Atoms are highlighted in cubes.
- Torsion 1 is defined.
- 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 numbers 1, 6, 25, and 26.
- Atoms are highlighted in cubes.
- Torsion 2 is defined.
- In the Coordinate Scan panel, go to the Potential tab.
-
For Solvent, choose None.
- A gas phase minimization will be run.
- For Job name, type mmod_Ddrive.
-
Click Run.
- This job takes ~30 seconds.
- A banner appears and a new group mmod_Ddrive-out is added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion once the job is completed.
Optional: Run this job in water and then chloroform solvents to compare intramolecular interactions in the conformations.
Note: MacroModel also has substantial options for controlling force fields, constraints, and settings associated with navigating a potential energy surface.
- Close the Coordinate Scan panel.
3.2 Analyze torsional scan results
-
In the top right corner of 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.
- In the dialog box, click Choose.
- Search and select Potential Energy-OPLS3 and Relative Energy-OPLS3 from the list.
-
Click OK.
- The Potential Energy and Relative Energy properties are added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion table.
- Both properties are given in kJ/mol.
- You may need to click and drag the edge of the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion table to see the added properties.
Note: OPLS4 calculations will still use this OPLS3 property category for the time being while we update the remainder of the panel. Notice that once imported, the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion property headings are more general.
The relative energies of small molecule conformations play a crucial role in determining the shape, function, and activity of a ligand.
Note: The last entry in the mmod_Ddrive-out
group is the lowest energy structure found from scanning the two dihedrals. The lowest energy structure gives us an indication of how adaptable the molecule is and how much strain is needed for it to adopt a particular conformation.
- Double-click the Ddrive entry title and rename it to Ddrive_Starting_Structure.
- Right-click the renamed entry and choose Move to > Other row.
- The Move to Row dialog box opens.
- For the Row number, type 10.
-
Click Move.
- The Ddrive_Staring_Structure is moved to the mmod_Ddrive-out group as its top entry.
- It is included and selected.
- Select the mmod_Ddrive-out group by clicking on the group heading.
-
Go to Tasks > Browse > Structure Alignment > Superposition.
- The Superposition panel opens.
- For Superimpose entries from, choose Project Table (selected).
- For Reference structure, choose 9: Ddrive_Starting_Structure.
- Check the box for Add property to Project Table option.
- For Choose method, choose Substructures.
- For Define structures for superposition using, choose SMARTS.
- Return to the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
-
Type A on the keyboard.
- Atom selection mode is activated.
- In the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed, use shift-click to 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 all atoms in the pyrimidine core.
- In the Superposition panel, click Get from Selection.
-
Click Superimpose Structures.
- The selected structures are aligned to the pyrimidine core.
- RSMDs for the structures are shown in the Results table at the bottom of the panel.
- Close the Superposition panel.
Figure 3-14. Fixing the Ddrive_Starting_Structure in the Workspace and steeping through the conformations.
-
In the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, double-click the In circle next to Ddrive_Starting_Structure.
- The structure is fixed 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 second entry.
- Use the left and right arrow keys to step through the conformations.
Note: Use shift-click to include contiguous entries or ctrl-click (cmd-click) to include non-contiguous entries in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed. The 300 lowest energy conformations are shown in the adjacent figure.
3.3 View the potential energy surface
A coordinate scan produces a contour diagram describing the molecular mechanics potential energy of a structure, relative to the value of either one or two coordinates (distances, angles, or dihedrals), can be generated with MacroModel. These exercises demonstrate how to produce a contour diagram describing the variation in energy of a molecule with respect to rotation of two dihedral angles.
-
Go to Tasks > Browse > MacroModel > Plot Coordinate Scan.
- The Plot Coordinate Scan Results panel opens.
Note: The potential energy surface of a compound can help with understanding key aspects of a molecule's behavior.
-
At the top of the panel, click Load Results.
- Select a grid file dialog box opens.
-
Double-click mmod_Ddrive folder and choose mmod_
Ddrive-out.grd. -
Click Open.
- The potential energy landscape is shown, colored rainbow.
- Low energies are shown in blue, high energies are shown in red.
- Set the Energy scale to Absolute.
Note: The energy scale is absolute because we are only studying a single molecule, rather than comparing several molecules.
Note: This plot shows there are large patches of lower energy with relatively sparse contour spacing between them. This suggests a reasonable freedom of movement of the two torsions relative to each other.
- Close the Plot Coordinate Scan Results panel.
-
Scroll to the bottom of the mmod_Ddrive-out group and includethe entry is represented in the Workspace, the circle in the In column is blue the Minimum E. structure during drive entry.
- The Ddrive_Starting_Structure should still be fixed in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Go to the Style Toolbox and choose Apply CPK representation.
![]()
Figure 3-20. Tiled view of a low and high energy conformation, showing a CH2 steric clash in yellow.
- Tile the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed by using the Workspace Layout Toggles (bottom right corner).
- The Workspacethe 3D display area in the center of the main window, where molecular structures are displayed is tiled.
- The steric clash between the CH2 group and phenyl ring can be visualized.
Note: Looking into the energetic components of the system can help explain a potential energy surface.
- Un-tile the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed by deselecting Tile by.
- Go to Workspace and choose Clear Workspace.
-
Type ctrl-T (cmd-T).
- The 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 dataopens.
-
Click Tree to open the Property Tree.
- Different calculated properties can be toggled on and off.
- Click the arrow next to each application to view more properties.
-
Double-click All.
- All the properties are now hidden.
- Open the tree to MacroModel > Secondary.
-
Under Secondary, check Bend Energy-S-OPLS, RMS Derivative-S-OPLS and Van der Waal Energy- S-OPLS.
- 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 is updated with the chosen properties.
There is a high VdW energy between the phenyl group and the nearby CH2 group of Torsion 1. This causes the pyrimidine core to distort and result in a high energy bend.
4. Screening for Active Conformations using ConfGen
As part of a medicinal chemistry exploration project, it is often helpful to change a single-atom linker in a molecule to evaluate the impact it has. We will change the linker atom in a series of molecules and visualize the variation of active conformations. These torsion differences often have implications for the shape of the ligand. We will be using ConfGen instead of MacroModel due to its speed and ability to generate bioactive conformations. You may choose to also perform the steps in this section using MacroModel to compare your results to Section 3 and to the ConfGen conformations generated here in Section 4. While it's impossible for a conformer search algorithm to determine a flexible ligand's bioactive conformer with absolute confidence, carefully considered search criteria do allow an algorithm to reject conformers likely to be high energy or inactive. Beyond merely expediting the conformer search process, this approach creates efficiently sized conformer sets that nevertheless contain a reasonable approximation of the bioactive geometry.
4.1 Set up the screening
- In the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, expand 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 Phenyl-linker-Pyridyl group by clicking on the group heading.
- Right-click the selection and choose Include.
- All the entries in the group are included in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
-
Go to Tasks > Browse > Structure Analysis > Bioactive Search.
- The ConfGen panel opens.
- For Use structures from, choose Workspace (4 entries).
- Click Load.
- Set the Target number of conformations to 50.
- Change the Job name to Phenyl-linker-Pyridyl_ConfGen
-
Click Run.
- This job takes ~30 seconds.
- A banner appears and a new group Phenyl-linker-Pyridyl_ConfGen-in-out 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 confGen panel.
-
In the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion, Ctrl+Click (Cmd+Click) to includethe entry is represented in the Workspace, the circle in the In column is blue a generated conformation of Ph-S-Py group and Ddrive_Starting_Structure.
- C-S has a greater bond length than the original C-S bond and the original Ph-S-Py compound.
Optional: Select Measure to measure and compare bond lengths in the structures.
4.2 Analyze the results of the linker series
-
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 table, click Change table settings (three vertical dots) and choose Show Property.
- The Show Properties in Table dialog box opens.
- In the dialog box, click Choose.
- Search and select energy and relative energy properties from the list.
-
Click OK.
- The Energy and Relative Energy properties are added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion table.
- Both properties are given in kJ/mol.
- You may need to click and drag the edge of the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion table to see the added properties.
- Right-click energy and choose Sort Selected (Ascending).
- Go to Workspace and choose Clear Workspace.
- 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 Phenyl-linker-Pyridyl_ConfGen group by clicking on the group heading.
- Right-click the selection and choose Include.
-
Click Continue to the warning box that appears.
- The entire group is included in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
-
Click the Workspace Configuration Toggle (bottom right corner) and tile your Workspacethe 3D display area in the center of the main window, where molecular structures are displayed by group.
- Four groups with overlaid structures are present in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed.
- Rotate the molecules in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed
- Right-click and select Clear Workspace
Note: These structures are not aligned, but already you can see differences in the conformation due to changes in the linker atoms.
Optional: Align each group in turn to a common phenyl core to more easily see changes in one conformation to the next.
- Using cmd/ctrl-click, includethe entry is represented in the Workspace, the circle in the In column is blue the lowest energy conformer in Ph-NH2-Py, Ph-CH2-Py, and Ph-S-Py groups using the energy property.
-
Double click the Presets button to apply Maestro presets to these compounds
- Notice differences in these active conformations.
- Un-tile the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed by expanding the Workspace Configuration (plus icon in the bottom right corner) toggle and deselecting Tile by.
-
With the lowest energy conformations still included in the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed, go to Tasks > Browse > Structure Alignment > Superimposition.
- The Superposition panel opens.
- For Superimpose entries from, choose Workspace.
- For Reference Structure, choose Ph-NH2-Py.
- For Choose method, select Substructures.
- For Define structures for superposition using, choose SMARTS.
- In the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed, use shift-click to 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 all pyridine carbon atoms.
- In the Superposition panel, click Get from Selection.
- Choose Superimpose Structure.
Note: Changing the linker from an oxygen (in the structure with grey carbons) to a sulfur (in the structure with green carbons), alters the phenyl ring orientation with respect to the pyridyl nitrogen. The CH2 linker induces the least planar conformation, due to sterics, while the NH2 linker is most coplanar with the pyridyl ring.
Note: This was done via ConfGen so that only active conformations were chosen. Different results will be returned if the same comparison is done using MacroModel.
Optional: A Phenyl-O-Group has also been included in the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion. If you would like more practice with MacroModel or ConfGen please use these structures to investigate the CH3,CFH2,CF2H, and CF3 groups and how they affect torsion parameters and active conformations of these ligands.
5. Conclusions and References
MacroModel is a powerful tool to explore the full range of conformations for small molecules, and can also be used to analyze protein-ligand complexes. You learned how specific settings can make a difference to the conformational output, so understanding the choices you make are important to the results. Through the exploration of small molecule active conformations using ConfGen, you observed how conformations can vary due to a single atom change and how this information can be used in lead optimization.
For further learning:
- Introduction to Structure Preparation and Visualization
- Structure-Based Virtual Screening Using Glide
- Understanding and Visualizing Target Flexibility
- Approximating Protein Flexibility without Molecular Dynamics
- Directions on how to use the Embrace.mae and Embrace.sbc files included in the tutorial zip file can be found in the Protein-Ligand Complex Calculations with MacroModel Embrace section in the documentation.
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