MM-ISMSA

MM-ISMSA provides an ultrafast scoring function to evaluate protein-protein docking binding energies and per-residue energy decompositions using molecular mechanics, implicit solvation, and surface-area contributions.

Key Features:

• Molecular Mechanics (MM) Component: Employs a 12-6 Lennard-Jones potential to compute non-bonded interactions in protein complexes.
• Electrostatics and Solvation Model: Incorporates an implicit solvent model (ISM) with individual desolvation penalties per partner and a hydrogen-bonding term for electrostatic assessment.
• Surface Area (SA) Contribution: Accounts for loss of water contacts upon complex formation via a surface-area term.
• Energy Decomposition: Computes per-residue and all-vs-all pairwise decompositions of total binding energies and interaction terms.
• Performance Evaluation: Benchmarked against APBSA and MM-PB(GB)SA on 1,242 decoy complexes comparing total binding energies, electrostatic terms, charge-charge interactions, desolvation energies, and per-residue contributions.
• Docking Model Selection Accuracy: On a test set of 15 protein-protein complexes with decoys, achieved a correlation of r² ∼ 0.9 with established methods and identified the best docking solution in 80% of cases.
• Computational Efficiency: Requires approximately 0.2–5.0 seconds per calculation depending on residue count, enabling extensive energy decompositions.
• Molecular Dynamics Trajectory Processing: Processes MD trajectories to produce averaged results with standard deviations.

Scientific Applications:

• Protein-Protein Docking Scoring: Scores and ranks docking models by total binding energy.
• Binding Energy and Interaction Analysis: Quantifies Lennard-Jones, electrostatic, charge-charge, desolvation, hydrogen-bonding, and surface-area contributions to binding.
• Per-Residue Energetic Decomposition: Identifies residue-level contributions to complex stability.
• Docking Model Validation: Selects docking solutions closest to experimental structures based on binding energy.
• MD-Based Interaction Analysis: Analyzes averaged energetic properties and their variability over molecular dynamics trajectories.

Methodology:

Scoring combines a Molecular Mechanics 12-6 Lennard-Jones potential, electrostatics modeled with an implicit solvent model (ISM) including individual desolvation penalties and a hydrogen-bonding term, a surface-area loss-of-water contacts term, and computes per-residue and all-vs-all pairwise energy decompositions and averages over molecular dynamics trajectories.