MM-ISMSA: An Ultrafast and Accurate Scoring Function for Protein–Protein Docking

Abstract

An ultrafast and accurate scoring function for protein-protein docking is presented. It includes (1) a molecular mechanics (MM) part based on a 12-6 Lennard-Jones potential; (2) an electrostatic component based on an implicit solvent model (ISM) with individual desolvation penalties for each partner in the protein-protein complex plus a hydrogen bonding term; and (3) a surface area (SA) contribution to account for the loss of water contacts upon protein-protein complex formation. The accuracy and performance of the scoring function, termed MM-ISMSA, have been assessed by (1) comparing the total binding energies, the electrostatic term, and its components (charge-charge and individual desolvation energies), as well as the per residue contributions, to results obtained with well-established methods such as APBSA or MM-PB(GB)SA for a set of 1242 decoy protein-protein complexes and (2) testing its ability to recognize the docking solution closest to the experimental structure as that providing the most favorable total binding energy. For this purpose, a test set consisting of 15 protein-protein complexes with known 3D structure mixed with 10 decoys for each complex was used. The correlation between the values afforded by MM-ISMSA and those from the other methods is quite remarkable (r(2) ∼ 0.9), and only 0.2-5.0 s (depending on the number of residues) are spent on a single calculation including an all vs all pairwise energy decomposition. On the other hand, MM-ISMSA correctly identifies the best docking solution as that closest to the experimental structure in 80% of the cases. Finally, MM-ISMSA can process molecular dynamics trajectories and reports the results as averaged values with their standard deviations. MM-ISMSA has been implemented as a plugin to the widely used molecular graphics program PyMOL, although it can also be executed in command-line mode. MM-ISMSA is distributed free of charge to nonprofit organizations.

Efficient Synthesis of Muramic and Glucuronic Acid Glycodendrimers as Dengue Virus Antagonists

“…[23] Finally the 3D geometry of the complexf ormed by the DENV2 with glycodendrimer 25 was deduced by using AutoDock Vina. [25] DENV2-glycodendrimer 25 interaction energy was measured during the simulation and it remained nearly constant around À110kcal mol À1 ,i ndicating the pose of the dendrimer in the protein was energeticallyf avorable (Table 1). Structure stabilityo ft he complex was evaluated by calculating the root-mean-square-deviation (RMSD) of the Ca atoms along the trajectory.…”

Activation of Human Toll-like Receptor 4 (TLR4)·Myeloid Differentiation Factor 2 (MD-2) by Hypoacylated Lipopolysaccharide from a Clinical Isolate of Burkholderia cenocepacia

“…The resulting structure was submitted to MDS with implicit water and AMBER* force field. Energy analysis was performed by means of the MM-ISMSA method (50).…”

Design and synthesis of potent hydroxamate inhibitors with increased selectivity within the gelatinase family

“…44,45 Chemistry General procedures. Periodic boundary conditions were applied and electrostatic interactions were treated using the smooth particle mesh Ewald method 42 with a grid spacing of 1 Å.…”