Improving protein-ligand docking results using the Semiempirical quantum mechanics: testing on the PDBbind 2016 core
Zainab Mohebbinia1, Rohoullah Firouzi1, Mohammad Hossein Karimi-Jafari2
1Department of Physical Chemistry, Chemistry and Chemical Engineering Research Center of Iran, Tehran, Iran.
Journal of Biomolecular Structure & Dynamics
|January 2, 2024
Summary
This study optimizes molecular docking poses using semi-empirical quantum mechanics (SQM) for drug design. Combining AutoDock Vina and AutoDock4 offers improved accuracy in predicting ligand binding and aids virtual screening.
Area of Science:
- Computational chemistry
- Drug discovery
- Molecular modeling
Background:
- Molecular docking is crucial for predicting ligand binding in drug design.
- Accurate pose prediction and affinity ranking are essential for virtual screening.
- Refinement of docked poses is necessary to remove steric clashes and improve accuracy.
Purpose of the Study:
- To evaluate the effectiveness of semi-empirical quantum mechanics (SQM) method PM7 for optimizing molecular docking poses.
- To compare the sampling and ranking capabilities of AutoDock4 and AutoDock Vina.
- To propose an improved protocol for structure-based virtual screening by combining both docking programs.
Main Methods:
- Geometry optimization of top-ranked docking poses using the PM7 SQM method.
- Utilizing the PDBbind core set (2016) for high-quality protein-ligand complexes.
- Comparative analysis of AutoDock4 and AutoDock Vina performance in pose generation and ranking.
Main Results:
- Docking pose optimization using PM7 enhances prediction accuracy and refines docked complexes by eliminating clashes.
- AutoDock Vina demonstrates superior sampling power for accurate ligand pose generation (RMSD ≤ 2.0 Å).
- AutoDock4 shows better ranking power compared to AutoDock Vina.
Conclusions:
- Geometry optimization with SQM methods like PM7 is beneficial for improving docking pose accuracy.
- A hybrid approach combining AutoDock Vina's sampling and AutoDock4's ranking provides a robust protocol for virtual screening.
- The proposed protocol enhances the reliability of structure-based virtual screening in drug design.
Related Concept Videos
Ligand Binding Sites
12.9K
Proteins are dynamic macromolecules that carry out a wide variety of essential processes; however, the activities of most proteins depend on their interactions with other molecules or ions, known as ligands.
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
12.9K
The Equilibrium Binding Constant and Binding Strength
12.9K
The equilibrium binding constant (Kb) quantifies the strength of a protein-ligand interaction. Kb can be calculated as follows when the reaction is at equilibrium:
12.9K
Conserved Binding Sites
4.2K
Many proteins’ biological role depends on their interactions with their ligands, small molecules that bind to specific locations on the protein known as ligand-binding sites. Ligand-binding sites are often conserved among homologous proteins as these sites are critical for protein function.
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
4.2K
Protein-protein Interfaces
12.5K
Many proteins form complexes to carry out their functions, making protein-protein interactions (PPIs) essential for an organism's survival. Most PPIs are stabilized by numerous weak noncovalent chemical forces. The physical shape of the interfaces determines the way two proteins interact. Many globular proteins have closely-matching shapes on their surfaces, which form a large number of weak bonds. Additionally, many PPIs occur between two helices or between a surface cleft and a...
12.5K
Ligand Binding and Linkage
4.8K
Allosteric proteins have more than one ligand binding site; the binding of a ligand to any of these sites influences the binding of ligands to the other sites. When a protein is allosteric, its binding sites are called coupled or linked. In the case of enzymes, the site that binds to the substrate is known as the active site and the other site is known as the regulatory site. When a ligand binds to the regulatory site, this leads to conformational changes in the protein that can influence...
4.8K
Protein Organization
6.5K
Proteins are polymers of amino acid residues. They are versatile and responsible for different cellular functions, including DNA replication, molecular transport, catalysis, and structural support. Proteins have a hierarchical structure comprising at least three levels of organization: primary, secondary, and tertiary structure. Some large proteins have a quaternary structure where individual protein subunits are linked together.
The primary structure of a protein is its amino acid sequence....
The primary structure of a protein is its amino acid sequence....
6.5K


