Binding free energy calculation for duocarmycin/DNA complex based on the QPLD-derived partial charge model.

Area of Science:

• Computational Chemistry
• Molecular Dynamics Simulations
• Drug Discovery

Background:

• Molecular dynamics (MD) simulations are crucial for understanding drug-target interactions.
• Accurate charge models are essential for reliable MD simulation results.
• The duocarmycin-DNA complex is a relevant system for cancer therapy research.

Purpose of the Study:

• To compare the performance of two charge models, RESP and QM-polarized ligand docking (QPLD), in MD simulations of the DNA/duocarmycin complex.
• To evaluate the accuracy of binding free energy (DeltaG(bind)) predictions using both models.
• To assess the stability of the simulated DNA/duocarmycin complex trajectories.

Main Methods:

• Performed 3-nanosecond unrestrained molecular dynamics (MD) simulations.
• Utilized two distinct charge models: the classic RESP charge model and the QM-polarized ligand docking (QPLD)-based charge model.
• Analyzed root-mean-square deviations (RMSDs) of trajectories and calculated binding free energies (DeltaG(bind)).

Main Results:

• The QPLD model demonstrated superior performance over the RESP model in MD simulations.
• Trajectories simulated using the QPLD model exhibited better stability.
• The QPLD model provided more accurate DeltaG(bind) estimations (-16.11 kcal/mol) compared to the RESP model (-10.05 kcal/mol).

Conclusions:

• The QM-polarized ligand docking (QPLD) charge model offers significant advantages for simulating DNA/duocarmycin interactions.
• QPLD-based MD simulations lead to more reliable predictions of binding affinity and complex stability.
• This study highlights the importance of advanced charge models in computational drug discovery.