For molecular modelers diving into advanced Molecular Dynamics (MD) simulations, center-of-mass (COM) pulling is often essential. However, the setup process—ensuring the system is ready, the pulling parameters are correctly defined, and the model is prepared—can feel like a daunting task. With SAMSON and the GROMACS Wizard, the process becomes significantly streamlined. In this post, we’ll focus on preparing a system for COM pulling efficiently.
Why is COM Pulling Important?
COM pulling allows molecular dynamics simulations to explore binding, unbinding, or mechanical stretching events by applying directed forces. These insights are pivotal for studying protein-ligand interactions, DNA mechanics, or inter-molecular dynamics observed under specific forces. The challenge, though, lies in setting up a system that runs efficiently, provides plausible initial geometry, and is physiologically relevant. Here’s where SAMSON eases the pain.
A Streamlined Workflow Using SAMSON GROMACS Wizard
The key is breaking the setup into manageable steps. Let’s focus on chain-specific pulling: separating chain A from chain B in the PDB model 2BEG, as detailed in the documentation.
Step 1: Define the Simulation Box With Pulling in Mind
Box dimensions are critical for simulations involving periodic boundary conditions. If parts of a molecule overlap due to periodicity during the pull, your simulation results will be distorted. In this example:
- The pull direction (z-axis) requires ample space, at least double the pulling distance, accounting for buffer room to prevent boundary-related artifacts.
- Set the box size to
6.5 nm × 4.5 nm × 12 nm, offering enough length for a 5 nm pull in the z-direction.
Here’s how the box looks in the SAMSON viewport:
Step 2: Adding Solvent and Neutralizing the System
Your molecular system functions within a simulated solution. Add solvent (e.g., water) before balancing charges using specific counterions:
- Apply a solvent model like
TIP3P. - Add Na and Cl ions for neutralization, alongside a salt concentration of
0.1 mol/liter.
By incorporating these conditions, the simulation mimics biologically relevant environments more effectively.
Step 3: Specify Index Groups for Pulling
To perform COM pulling, you need to define molecular groups that interact during the simulation. For chain-specific pulling:
- Create an index group for
CHAIN A(the group being pulled). - Create another index group for
CHAIN B(the anchor chain).
SAMSON’s “Generate based on current selection in document” feature makes it intuitive to assign regions in the model:
Step 4: Configure Pulling Parameters
Once the system and index groups are ready, the next step is defining COM pulling parameters. Here’s what was applied in our setup:
- Pull
CHAIN AfromCHAIN Bin the z-direction for5 nm. - Rate:
0.01 nm/psover0.5 ns. - Force constant:
1000 kJ mol^-1 nm^-2.
The diagram below shows how SAMSON simplifies grouping and parameter editing.
What's Next?
With these steps, the system is ready for equilibration and production runs emphasizing COM pulling. After running the simulation, you’ll observe force evolution and displacement generated by the pulling force:
COM pulling results like these can also be refined further using advanced analysis, such as umbrella sampling, to compute potentials of mean force along reaction coordinates. This opens doors to deeper insight into molecular interactions under force.
To explore the full workflow in more detail, consult the original SAMSON COM Pulling Documentation.
Note: SAMSON and all SAMSON Extensions are free for non-commercial use. You can get SAMSON at https://www.samson-connect.net.