Filtering and Validating Protein-Protein Docked Poses with Macromolecular Pose Filtering and MM-GBSA Residue Scanning
Tutorial Created with Software Release: 2026-2
Topics: Antibody Design , Biologics Drug Discovery , Structure Prediction & Target Enablement
Products Used: PIPER
|
51 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 how to use alanine scanning data to filter and validate poses from a protein-protein docking job for prospective mutagenesis.
You will learn how to combine Macromolecular Pose Filtering and MM-GBSA Residue Scanning to triage and validate an ensemble of docked poses and identify those which most closely match with the mutational affinity patterns.
Tutorial Content
- Introduction to Validating Predicted Protein Complexes
- Creating Projects and Importing Structures
- Filtering the Docked Poses with the Macromolecular Pose Filtering Panel
- Preparing the Filtered Poses
- Validating the Filtered Poses with MM-GBSA Residue Scanning
- Conclusion and References
- Glossary of Terms
1. Introduction to Validating Predicted Protein Complexes
Accurate structural information is necessary for understanding biomolecular recognition and guiding rational design. While experimentally determined structures from X-ray crystallography, NMR, or cryo-EM provide the most reliable starting point for structural analysis, such data are not always available for every protein or protein-protein complex of interest. In many cases, computational structure prediction methods can provide useful alternatives, particularly when reliable homologous templates with similar sequences and folds exist. Advances in homology modeling, co-folding approaches, and machine learning-based predictors have significantly improved the reliability of predicted structures. However, predicting protein-protein complexes remains inherently more challenging than modeling individual proteins.
Given that experimental structures are not typically sought until much later in the lead optimization process, generating better complex models earlier can be highly valuable in driving forward a rational structure-based design approach. The most effective way to build confidence in a predicted model is to compare it against available experimental evidence or prior knowledge about the system. Experimental observations such as mutational binding data, cross-reactivity patterns, or known epitope and paratope residues can help determine whether the predicted interface is consistent with biological expectations. Similarly, experimental techniques such as hydrogen–deuterium exchange (HDX) can provide valuable constraints that support or challenge the predicted binding mode. Incorporating such information during model assessment helps identify biologically plausible conformations and ensures that subsequent computational analyses are based on structures that are consistent with known experimental behavior.
This tutorial uses the 2VXT structure, which is an antigen-bound murine antibody. In this tutorial, you will learn how to filter and validate an ensemble of antibody-antigen docked poses by leveraging data from alanine scanning. Here, 51 CDR residues of the 2VXT antibody were mutated site-by-site to alanine and the effect on target binding was measured. The assumption underlying this approach is that the residues that showed a significant loss in binding affinity are likely to play an important role in the interaction interface. Subsequently, you will map the binding poses with these interfacial residues and perform an in silico mutagenesis to alanine on these poses to validate whether the proposed binding modes are consistent with the observed mutational effects.
2. Creating Projects and Importing Structures
-
Open BioLuminate and create a new project named
pose_filtering_tutorial.prjfor this tutorial.- Don’t know how? See First steps in Maestro.
- Download the tutorial zip file including input files and reference outputs here: Tutorial Files
- After downloading the zip file, unzip the contents in your Working Directorythe location where files are saved for ease of access throughout the tutorial.
We have provided 30 cluster centroid structures, which are the default output of the protein-protein docking job in BioLuminate. To proceed with this tutorial, you will import these structures from the tutorial archive.
- Go to File > Import Structures.
-
Navigate to
2VXT_Docking_Example-out.maegzin your Working Directorythe location where files are saved and click Open.- The structures are added to the Entriesa simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and the top entry 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.
3. Filtering the Docked Poses with the Macromolecular Pose Filtering Panel
In the alanine scanning analysis performed on 51 CDR residues of the 2VXT antibody, a total of 11 residues were identified as critical for antigen recognition because their mutation to alanine resulted in significant loss of binding affinity, while the remaining 40 mutations were neutral with respect to binding affinity. The critical residues are distributed between the heavy chain and the light chain:
- Heavy chain residues: H: TYR 32, H: PHE 33, H: ASP 52, H: TYR 53, H: ASP 56, H: ARG 94, H: GLY 95, H: ARG 101
- Light chain residues: L: LYS 32, L: TYR 49, L: TYR 96
In this section, you will use the Macromolecular Pose Filtering panel to screen the 30 docked poses, ensuring that the filtered poses maintain contacts with the critical residues. Note that the wide geometric distribution of the mutations across both the heavy and the light chain is quite advantageous for filtering the docked poses. It spreads the ‘binding points’ across a larger area, thereby significantly limiting the number of possible orientations that can be satisfied.
The following interactive tutorial guides you through setting up and running the pose filtering job using specific residue constraints. You will first define the critical heavy and light chain CDR residues at the interface, the mutations to which were identified as deleterious. You will then specify the contact cutoff and the minimum number of residues that must match the filter. Finally, you will run the filtering job and view the filtered poses.
Defining the cutoff and the minimum number of residues is the most critical part in the above workflow. If the parameters are set too stringently (for instance, a 2.0 Å cutoff and requiring all 11 residues to match the filter), the filtering may fail to return any poses that meet the criteria. Conversely, setting the parameters too loosely could result in too many poses being returned, negating the purpose of the filtering procedure.
As a rule of thumb, the cutoff values should be kept within 3.0 – 5.0 Å. This range is recommended to ensure residues are interfacial while providing leniency for slightly longer distances that may result from the imperfect geometry of rigid-body docking. Regarding the minimum number of residues, a suitable range is typically 50-75 % of the desired residues. This specific setting is less impactful on the final results because the analysis will display how many residues meet the criteria for all poses that satisfy the minimum requirement (e.g., if the minimum is set to 7, the results will include poses with 7, 8, 9,10, or 11 interfacial residues).
In instances where the constrained residues are all adjacent to each other, or if only a few residues are specified, the filtering may return a large number of poses that meet the criteria. In such cases, we recommend gradually decreasing the cutoff distance to select only those poses that exhibit the closest binding interactions. A cutoff of 2.0 Å is usually sufficient, but 1.5 Å is the recommended limit, as this distance approaches the typical equilibrium distances of many non-covalent interactions.
4. Preparing the Filtered Poses
Warning
Before validating the filtered poses through mutagenesis, the structures must be reprepared. This is necessary because the docking calculation converts the structures into .pdb format, which strips away some force field parameter assignments that need to be reinstated for subsequent MM-GBSA residue scanning calculations. Furthermore, since rigid-body docking does not reorient side chains or relax the structures, a small minimization step is required to eliminate clashes or misaligned non-covalent bond orientations.
In this section, you will reprepare the filtered poses using the Protein Preparation Workflow panel. It is important to note that the pKa of a particular titratable residue at the interface may differ between poses due to a change in its local environment. Therefore, to avoid calculating binding affinity changes from inconsistent wild-type states, you will modify the default settings to ensure that the wild-type state used for mutation is identical across all poses.
The following interactive tutorial guides you to set up and run the preparation job for the filtered poses with modified settings. To learn about the protein preparation workflow in detail, please consult the Introduction to Structure Preparation and Visualization tutorial.
5. Validating the Filtered Poses with MM-GBSA Residue Scanning
In this section, you will set up and run MM-GBSA residue scanning calculations on the filtered poses. The goal of running the MM-GBSA residue scanning calculations is to confirm consistency with experimental results: mutations predicted as deleterious by MM-GBSA should match those found to be deleterious in experiments, and likewise, neutral MM-GBSA mutations should align with experimentally neutral mutations that do not impact binding affinity.
Recall that of 51 CDR residues mutated to alanine, 11 mutations were found to be deleterious, causing a loss of binding, while the remaining 40 mutations were neutral with respect to binding affinity. For validation purposes, it is beneficial to incorporate as much experimental information as possible. Therefore, in addition to the deleterious mutations used for pose filtering, you will also test 6 neutral mutations, which are expected to have minimal impact on binding affinity when mutated to alanine. Neutral mutations can have no impact on binding affinity for one of two reasons: they have no impact on binding or that they are not at the interfaces. Since we do not know a priori whether these sites are interfacial or not, we do not include them in the initial filtering criteria; however, since we know their effect on binding we can use them as controls in the validation step.
To maximize the effectiveness of the mutational validation using MM-GBSA, you will select one residue per CDR loop to ensure broader geometric coverage of the binding interface:
- Heavy chain CDRs: H1 - SER 28, H2 - GLY 55, H3 - TYR 90
- Light chain CDRs: L1 - GLN 27, L2 - SER 52, L3 - SER 93
The following interactive tutorial guides you through running the MM-GBSA residue scanning calculations and analyzing the results for the first of the poses which passed the filtering.
Now that you have learned how to run an MM-GBSA Residue Scanning Calculation for one of the docked poses, follow the same steps to set up and run the calculations for pose 9 and pose 13. If you do not wish to run the job yourself, the pre-generated results for both the calculations are available in the tutorial archive. The results of the MM-GBSA ∆ Affinity predictions are summarised below:
|
Mutation |
Pose 1 |
Pose 9 |
Pose 13 |
|
H: TYR 32 – ALA |
13.67 |
2.85 |
6.59 |
|
H: PHE 33 – ALA |
11.06 |
6.75 |
6.25 |
|
H: ASP 52 – ALA |
6.70 |
0.67 |
-0.63 |
|
H: TYR 53 – ALA |
6.52 |
12.26 |
3.13 |
|
H: ASP 56 – ALA |
3.02 |
10.67 |
-2.67 |
|
H: ARG 94 – ALA |
19.99 |
2.58 |
10.33 |
|
H: GLY 95 – ALA |
19.53 |
-1.78 |
8.99 |
|
H: ARG 101 – ALA |
23.25 |
1.58 |
17.60 |
|
L: LYS 32 – ALA |
12.91 |
8.68 |
1.53 |
|
L: TYR 49 – ALA |
10.89 |
0.71 |
11.98 |
|
L: TYR 96 – ALA |
2.01 |
3.19 |
-1.02 |
|
H: SER 28 – ALA |
0.17 |
0.14 |
0.23 |
|
H: GLY 95 – ALA |
-0.06 |
-0.23 |
-0.05 |
|
H: TYR 90 – ALA |
-0.22 |
-0.14 |
-0.18 |
|
L: GLN 27 – ALA |
-0.16 |
-0.17 |
-0.17 |
|
L: SER 52 – ALA |
-0.53 |
-0.10 |
-0.22 |
|
L: SER 93 – ALA |
-0.39 |
-0.20 |
-0.31 |
|
Deleterious (11) |
10 |
5 |
7 |
|
Neutral Controls (6) |
6 |
6 |
6 |
Pose 1 showed the best match to the experimental data, with 9 out of 11 deleterious mutations and all 6 neutral mutations matching the desired profile. Poses 9 and 13 both matched the experimental profile for all 6 neutral mutations. However, the match for deleterious mutations was lower, with pose 9 matching 5 out of 11 and pose 13 matching 7 out of 11. This indicates that pose 1 is the most representative of the experimental structure out of the three models and can likely be used in prospective calculations with high confidence.
Given that the true bound state for this system is already known, we can directly compare the filtered and validated docked pose (Pose 1) against the crystal structure. This comparison reveals an interface RMSD of 2.18 Å between the two structures.
Figure: Pose 1 (green) aligned with the crystal structure pose (gray).
6. Conclusion and References
In this tutorial, you learned how to use prior experimental knowledge to effectively filter and validate predicted protein-protein complex structures. In many practical protein-protein modeling projects, the challenge extends beyond predicting the correct binding orientation of two partners and includes identifying the biologically relevant conformations of the individual proteins themselves. In such cases, each binding partner may first need to be modeled independently using approaches such as AlphaFold, homology modeling, or related structure prediction methods. Because proteins can adopt multiple plausible conformations, several structures variants of one or both partners are often generated. Each possible combination of these structures can then be subjected to protein-protein docking, producing an ensemble of complexes. This substantially increases the size and diversity of the pose set and makes robust filtering strategies essential.
For further learning:
- Learning Path: Antibody Modeling
- Introduction to Computational Antibody Engineering online course (Course Page | Preview)
For further reading:
7. 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
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
Recent actions - This is a list of your recent actions, which you can use to reopen a panel, displayed below the Browse row. (Right-click to delete.)
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 table (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries
Working Directory - the location where files are saved
Workspace - the 3D display area in the center of the main window, where molecular structures are displayed


