Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins

Although molecular docking has evolved dramatically over the years, its application to glycosaminoglycans (GAGs) has remained challenging because of their intrinsic flexibility, highly anionic character and rather ill-defined site of binding on proteins. GAGs have been treated as either fully “rigid...

Full description

Bibliographic Details
Main Authors: Samuel G. Holmes, Umesh R. Desai
Format: Article
Language:English
Published: MDPI AG 2023-11-01
Series:Biomolecules
Subjects:
Online Access:https://www.mdpi.com/2218-273X/13/11/1633
_version_ 1797460015304409088
author Samuel G. Holmes
Umesh R. Desai
author_facet Samuel G. Holmes
Umesh R. Desai
author_sort Samuel G. Holmes
collection DOAJ
description Although molecular docking has evolved dramatically over the years, its application to glycosaminoglycans (GAGs) has remained challenging because of their intrinsic flexibility, highly anionic character and rather ill-defined site of binding on proteins. GAGs have been treated as either fully “rigid” or fully “flexible” in molecular docking. We reasoned that an intermediate semi-rigid docking (SRD) protocol may be better for the recapitulation of native heparin/heparan sulfate (Hp/HS) topologies. Herein, we study 18 Hp/HS–protein co-complexes containing chains from disaccharide to decasaccharide using genetic algorithm-based docking with rigid, semi-rigid, and flexible docking protocols. Our work reveals that rigid and semi-rigid protocols recapitulate native poses for longer chains (5→10 mers) significantly better than the flexible protocol, while 2→4-mer poses are better predicted using the semi-rigid approach. More importantly, the semi-rigid docking protocol is likely to perform better when no crystal structure information is available. We also present a new parameter for parsing selective versus non-selective GAG–protein systems, which relies on two computational parameters including consistency of binding (i.e., RMSD) and docking score (i.e., GOLD Score). The new semi-rigid protocol in combination with the new computational parameter is expected to be particularly useful in high-throughput screening of GAG sequences for identifying promising druggable targets as well as drug-like Hp/HS sequences.
first_indexed 2024-03-09T16:59:12Z
format Article
id doaj.art-fef52af94f8c48c1ab9668d1431ddb1c
institution Directory Open Access Journal
issn 2218-273X
language English
last_indexed 2024-03-09T16:59:12Z
publishDate 2023-11-01
publisher MDPI AG
record_format Article
series Biomolecules
spelling doaj.art-fef52af94f8c48c1ab9668d1431ddb1c2023-11-24T14:32:02ZengMDPI AGBiomolecules2218-273X2023-11-011311163310.3390/biom13111633Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto ProteinsSamuel G. Holmes0Umesh R. Desai1Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USADepartment of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USAAlthough molecular docking has evolved dramatically over the years, its application to glycosaminoglycans (GAGs) has remained challenging because of their intrinsic flexibility, highly anionic character and rather ill-defined site of binding on proteins. GAGs have been treated as either fully “rigid” or fully “flexible” in molecular docking. We reasoned that an intermediate semi-rigid docking (SRD) protocol may be better for the recapitulation of native heparin/heparan sulfate (Hp/HS) topologies. Herein, we study 18 Hp/HS–protein co-complexes containing chains from disaccharide to decasaccharide using genetic algorithm-based docking with rigid, semi-rigid, and flexible docking protocols. Our work reveals that rigid and semi-rigid protocols recapitulate native poses for longer chains (5→10 mers) significantly better than the flexible protocol, while 2→4-mer poses are better predicted using the semi-rigid approach. More importantly, the semi-rigid docking protocol is likely to perform better when no crystal structure information is available. We also present a new parameter for parsing selective versus non-selective GAG–protein systems, which relies on two computational parameters including consistency of binding (i.e., RMSD) and docking score (i.e., GOLD Score). The new semi-rigid protocol in combination with the new computational parameter is expected to be particularly useful in high-throughput screening of GAG sequences for identifying promising druggable targets as well as drug-like Hp/HS sequences.https://www.mdpi.com/2218-273X/13/11/1633heparin/heparan sulfatemolecular dockingglycosaminoglycansknowledge-based docking
spellingShingle Samuel G. Holmes
Umesh R. Desai
Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins
Biomolecules
heparin/heparan sulfate
molecular docking
glycosaminoglycans
knowledge-based docking
title Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins
title_full Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins
title_fullStr Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins
title_full_unstemmed Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins
title_short Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins
title_sort assessing genetic algorithm based docking protocols for prediction of heparin oligosaccharide binding geometries onto proteins
topic heparin/heparan sulfate
molecular docking
glycosaminoglycans
knowledge-based docking
url https://www.mdpi.com/2218-273X/13/11/1633
work_keys_str_mv AT samuelgholmes assessinggeneticalgorithmbaseddockingprotocolsforpredictionofheparinoligosaccharidebindinggeometriesontoproteins
AT umeshrdesai assessinggeneticalgorithmbaseddockingprotocolsforpredictionofheparinoligosaccharidebindinggeometriesontoproteins