MartiniSurf: Automated Simulations of Surface-Immobilized Biomolecular Systems with Martini

MartiniSurf is an open-source command-line framework that automates the construction of GROMACS-ready, coarse-grained simulation systems for biomolecules immobilized on solid supports, enabling high-throughput exploration of surface-tethered protein and DNA dynamics with controlled orientations and attachment geometries.

The Gist

Imagine you are a master chef trying to bake the perfect cake. You have a delicious recipe (the biomolecule, like an enzyme or DNA), but you need to stick it to a specific spot on your baking tray (the solid surface) so it doesn't slide around.

The problem is, if you stick the cake in the wrong spot, or if you glue it down too tightly, it might not bake properly. It might lose its shape, or the heat might not reach the center. In the real world, scientists do this with tiny biological machines to make medicines or sensors, but figuring out the perfect way to stick them down is incredibly difficult and time-consuming.

Enter MartiniSurf: The "Smart Glue" Robot

This paper introduces a new tool called MartiniSurf. Think of it as a highly sophisticated, automated robot chef that helps scientists design exactly how to stick these tiny biological machines to surfaces.

Here is how it works, broken down into simple concepts:

1. The Problem: The "Jenga Tower" Dilemma

Scientists want to study how these biological machines behave when stuck to a surface. To do this, they use computer simulations.

• The Old Way: It was like building a Jenga tower by hand, one block at a time, using different tools for every single step. You had to manually grab the molecule, manually glue it, manually check the angles, and manually build the surface. If you made one tiny mistake, the whole tower fell, and you had to start over. It was slow, messy, and hard to repeat.
• The New Way (MartiniSurf): This tool is like a 3D printer that knows the blueprint perfectly. You tell it, "I want to stick this protein to a graphene surface using these three specific points," and it builds the entire digital world for you in seconds.

2. The "Martini" Scale: Zooming Out

To simulate these systems, scientists can't look at every single atom (that's like trying to count every grain of sand on a beach). Instead, they use the Martini Force Field, which is like looking at the beach from a helicopter.

• Instead of seeing individual sand grains, you see "beads" of sand.
• MartiniSurf takes the detailed blueprint of a protein or DNA and automatically converts it into these "beads." This makes the computer simulation run thousands of times faster, allowing scientists to watch the molecule move and dance over time.

3. The Three Magic Modes

MartiniSurf offers three different ways to "glue" the molecule to the surface, giving scientists total control:

• The "Velcro" Mode (Anchor-based): Imagine sticking the molecule down with invisible elastic bands attached to specific points (like specific amino acids). The molecule is held in place but can still wiggle a bit. This is great for testing how different attachment points affect the molecule's stability.
• The "Rope" Mode (Linker-based): Sometimes, you don't want to stick the molecule directly to the surface. You want a little rope (a chemical linker) in between. MartiniSurf can build this rope for you. This is crucial because the length and flexibility of the rope can change how well the enzyme works.
• The "Magnet" Mode (Adsorption): Sometimes you just want to see what happens if you drop the molecule onto the surface and let it stick naturally, without forcing it. This mode simulates that natural attraction.

4. Building the Stage

Before the molecule can be stuck down, you need a stage. MartiniSurf can build different types of stages:

• Flat Surfaces: Like a smooth table (graphene) or a bumpy carpet (agarose).
• Curved Surfaces: Like wrapping the molecule around a cylinder (carbon nanotubes).
• Custom Decor: You can even tell the robot to paint the surface with specific chemical "stickers" to see how they interact with the molecule.

5. Why This Matters

Before MartiniSurf, designing these experiments was like trying to assemble a complex Lego set without instructions, using parts from five different boxes.

• Reproducibility: Now, if one scientist in Spain builds a simulation, a scientist in Japan can type the exact same command and get the exact same result. No more "it worked on my computer" issues.
• Speed: What used to take days of manual setup now takes minutes.
• Discovery: Because it's so fast and easy, scientists can now test hundreds of different "gluing" strategies to find the one that makes the enzyme work best. This helps in designing better medical treatments, more efficient industrial enzymes, and advanced biosensors.

The Bottom Line

MartiniSurf is the ultimate "set-it-and-forget-it" tool for scientists. It takes the messy, manual work of building digital models of sticky biological molecules and turns it into a simple, automated command. It allows researchers to focus less on building the simulation and more on understanding how these tiny machines work, ultimately helping us create better technologies for health and industry.

It's like giving scientists a remote control for the microscopic world, letting them instantly test "What if we stick it here?" or "What if we use a longer rope?" without ever having to leave their chairs.

Technical Summary

1. Problem Statement

The rational design of biomolecule immobilization strategies (e.g., enzymes or DNA on solid supports) requires a molecular-level understanding of how surface properties, tethering geometry, and structural dynamics influence stability and function. While coarse-grained (CG) molecular dynamics (MD) simulations using the Martini force field offer an efficient framework for studying these interactions, constructing reproducible, orientation-controlled immobilized systems remains technically challenging.

• Current Limitations: Existing workflows often rely on manual preprocessing, custom scripts, and non-standardized orientation procedures.
• The Gap: There is a lack of a unified, automated framework capable of systematically generating surface-immobilized systems with controlled multivalent anchoring, explicit linkers, and diverse surface chemistries (e.g., graphene, agarose) for high-throughput exploration.

2. Methodology: The MartiniSurf Framework

MartiniSurf is an open-source, command-line interface (CLI) tool written in Python that automates the entire workflow for preparing surface-immobilized biomolecular systems for GROMACS simulations. It integrates tools from the Martini ecosystem into a single pipeline.

Core Workflow Components:

1. Input & Preprocessing:

   • Accepts local PDB files, RCSB PDB IDs, or UniProt accessions.
   • Automatically retrieves (e.g., via AlphaFold models) and cleans structures.
   • Supports Proteins (via martinize2), DNA (via martinize-dna), and pre-coarse-grained complexes.
   • Allows for structure-based potentials (GōMartini) or standard Martini 2/3 force fields.

2. Surface Generation:

   • Generates model solid supports using 2D hexagonal lattices.
   • Built-in Modes: Supports planar surfaces (e.g., graphene, graphite, agarose-like) with tunable lattice spacing, bead identity, charge distribution, and multilayer configurations.
   • Functionalization: Allows the addition of user-defined ligands or linker-like moieties to the surface via the --surface-linkers option.
   • Flexibility: Accepts user-provided coordinate/topology files for heterogeneous or experimentally inspired architectures.

3. Biomolecule Orientation & Immobilization:
   The framework aligns the biomolecule relative to the surface using three distinct, parameter-driven regimes:

   • Anchor-based (Implicit): Selected residues are restrained to the surface at a defined distance using harmonic restraints (simulating multivalent attachment).
   • Linker-based (Explicit): Introduces a coarse-grained linker molecule between the surface and the biomolecule, allowing for the modeling of linker flexibility and geometry.
   • Adsorption Mode: Orients the biomolecule toward the surface without explicit anchoring constraints.

4. System Assembly:

   • Automates solvation, ionization (with adjustable salt concentration), and optional substrate insertion.
   • Supports "frozen water" layers for DNA simulations.
   • Outputs a complete, GROMACS-ready directory containing topology files (.top), coordinate files (.gro), index groups, and molecular dynamics parameter files (.mdp).

3. Key Contributions

• Unified Automation: Consolidates structure retrieval, cleaning, coarse-graining, surface generation, and orientation into a single command-line workflow, minimizing user-dependent variability.
• Reproducibility: Explicit input parameters ensure that identical configurations always generate the same coordinate and topology files, facilitating systematic comparison across different immobilization geometries.
• Versatility:
  • Supports both Proteins (including GōMartini native-contact potentials) and DNA.
  • Handles diverse surface chemistries (carbon nanomaterials, polysaccharides) and tethering strategies (implicit restraints vs. explicit linkers).
• Accessibility: Provides cloud-based Google Colab notebooks for users without local HPC access, enabling system preparation and short validation simulations directly in the browser.

4. Results & Validation

The authors validated MartiniSurf through two representative case studies:

• Protein Immobilization (Alcohol Dehydrogenase - ADH):

  • Successfully reconstructed a previously studied system where ADH was immobilized on an agarose-like surface via multivalent harmonic restraints.
  • Extended the model to include chemically functionalized surfaces with explicit surface-bound groups and a complete catalytic assembly (Enzyme + NADH cofactor + Ethanol substrate).
  • Demonstrated the ability to model complex, multi-component systems consistently.

• DNA Immobilization:

  • Modeled a DNA decamer tethered to graphene surfaces (neutral and positively charged) using an explicit aliphatic linker.
  • Findings: The simulations reproduced known all-atom trends:
    • On neutral graphene, the DNA remained upright (small tilt angle).
    • On positively charged graphene, the DNA bent toward the surface due to electrostatic attraction (increased tilt angle, reduced distance).
  • This confirmed that MartiniSurf captures essential charge-dependent orientational responses while enabling simulations at extended temporal and spatial scales.

5. Significance

MartiniSurf addresses a critical bottleneck in computational biotechnology by lowering the technical barrier to constructing complex, surface-immobilized simulation systems.

• Rational Design: It enables the systematic, high-throughput in silico exploration of how surface chemistry and tethering geometry modulate biomolecular stability, dynamics, and function.
• Scalability: By leveraging coarse-grained models, it allows for the simulation of large assemblies and long timescales that are computationally prohibitive with all-atom methods.
• Future Impact: The modular architecture allows for easy integration of new surface materials (e.g., lipid bilayers, heterogeneous supports) and advanced force-field developments, positioning it as a robust platform for the rational design of immobilized biocatalysts, biosensors, and nanobiotechnological devices.

Availability: The source code, documentation, and usage examples are open-source and available at https://github.com/BioKT/MartiniSurf.