When preparing molecular systems for simulations, one common challenge faced by molecular modelers is deciding whether to retain or remove crystal waters. These water molecules can be critical in active sites or pose unnecessary computational overhead if irrelevant. The GROMACS Wizard in SAMSON provides a methodical and intuitive way to deal with this issue.
In this blog post, we’ll guide you through the careful selection and removal of crystal waters outside of functional active sites. The focus is on ensuring that your molecular model is both accurate and efficient for simulation tasks.
Why Remove Crystal Waters?
In protein structures derived from the Protein Data Bank (PDB), crystal waters are often included. While tightly bound waters in active sites might be essential for the protein’s function, many other waters are non-essential for simulations. Removing irrelevant waters helps reduce complexity, improving performance and focusing the simulation on meaningful interactions.
Step-by-Step Guide to Filtering Crystal Waters
Here’s how you can carefully filter out unnecessary waters while retaining important active-site molecules:
- Select the active site: Highlight the relevant structures (e.g., ligand), residues, or atoms that are part of the active site. For example, you might focus selection using the Document view or Viewport tools in SAMSON.
- Expand the selection: Right-click on your current selection and choose Expand selection > Advanced. This helps in isolating unnecessary waters based on their distance from the active site.
- Configure parameters: In the dialogue box, use Water as the node type and specify a distance threshold, such as waters beyond a certain distance from the active site. You can enable the auto-update option to preview the selection in real time.
- Verify and erase: Confirm your selection of unnecessary waters. Once satisfied, right-click the selection and select Erase selection. By doing so, you’ll ensure only meaningful waters remain.
The interface provides a user-friendly visualization, helping to minimize errors while performing selections. Here’s a snapshot of the process for removing waters beyond a specified distance:
Tips for Robust Modeling
While the above steps streamline water molecule selection, here are some tips to maximize accuracy:
- Preserve tightly bound waters: Functional water molecules near active sites (e.g., catalytic waters) should always be retained.
- Use distance-based filters: Employ the distance-filtering feature to refine water selections carefully and ensure no critical water molecules are accidentally removed.
- Validate decisions: Double-check removed molecules by revisiting the visualization to avoid accidental deletions.
Next Steps
Once you’ve optimized your molecular system by removing unnecessary waters, you’re ready to proceed to the Preparation step. This phase involves defining the simulation box, solvent, ions, and topology files, paving the way for your GROMACS workflow.
For further details and step-by-step guidance, explore the full documentation of this process at GROMACS Wizard – Pre-processing of the system.
Note: SAMSON and all SAMSON Extensions are free for non-commercial use. You can download SAMSON at https://www.samson-connect.net.