An MD View of Ligand Binding.

Area of Science:

• Biochemistry
• Structural Biology
• Computational Chemistry

Background:

• Protein-ligand complexes are typically characterized by specific bonding interactions in crystal structures.
• Understanding the dynamic behavior of these complexes is crucial for drug discovery and protein function studies.

Purpose of the Study:

• To investigate the dynamic behavior of a representative protein-ligand complex using molecular dynamics (MD) simulations.
• To compare the MD behavior and system energies of complexes starting from crystal structures versus docked structures, including native and analog ligands.

Main Methods:

• One-microsecond molecular dynamics simulations were performed on a protein-ligand complex.
• Simulations were initiated from crystal structures and structures with docked native ligand and ligand analogs.
• Molecular Mechanics with the Generalized Born Surface Area (MM/GBSA) method was used to calculate system energies.

Main Results:

• MD behaviors and calculated system energies were similar when simulations started from crystal or docked native ligand structures, though replicate simulations showed variability.
• Interatomic contact analysis revealed that some crystal contacts were rarely sampled, others intermittently, and new persistent contacts were formed.
• Docking of non-native ligand analogs resulted in similar dynamic behaviors and, in some cases, comparable calculated energies to the native ligand.

Conclusions:

• Ligands and proteins exhibit dynamic cooperation in forming binding interactions.
• The binding site can accommodate novel ligand analogs, maintaining similar interaction dynamics and energies.
• MD simulations provide insights into the plasticity and adaptability of protein-ligand interactions beyond static crystal structures.