Coupling Monte Carlo, Variational Implicit Solvation, and Binary Level-Set for Simulations of Biomolecular Binding

Area of Science:

• Computational chemistry and biophysics
• Molecular modeling and simulation
• Biomolecular interactions

Background:

• Accurate simulation of biomolecular binding in aqueous environments is crucial for understanding biological processes.
• Existing methods often face challenges in efficiently calculating solvation free energy and exploring binding configurations.
• Developing robust computational tools is essential for advancing drug discovery and molecular design.

Purpose of the Study:

• To introduce a novel hybrid computational approach for simulating biomolecular binding in aqueous solvents.
• To efficiently estimate the solvation free energy of biomolecular complexes using a variational implicit-solvent model (VISM) and a binary level-set method.
• To provide a foundation for more realistic simulations of ligand-protein interactions and binding pathways.

Main Methods:

• A hybrid approach combining the Monte Carlo (MC) method, a variational implicit-solvent model (VISM), and a binary level-set method.
• VISM functional minimization to estimate solvation free energy, considering solute volumetric, interfacial, van der Waals, and electrostatic contributions.
• A fast binary level-set method for efficient minimization of the VISM functional during MC moves, approximating surface area via convolution.

Main Results:

• The hybrid approach was applied to the p53-MDM2 system, approximating molecules as rigid bodies.
• The method successfully captured preliminary binding poses before the final bound state.
• Subsequent all-atom molecular dynamics simulations confirmed that these poses rapidly converged to the final bound state.

Conclusions:

• The developed hybrid approach represents a significant advancement toward realistic simulations of biomolecular interactions.
• The efficient estimation of solvation free energy and capture of binding poses pave the way for studying complex binding phenomena.
• Further refinements in coarse-graining, MC sampling, and integration with other models will enable detailed analysis of free-energy landscapes and kinetic pathways for ligand binding.