The impact of crystallization conditions on structure-based drug design: A case study on the methylene blue/acetylcholinesterase complex

The impact of molecular variants, crystallization conditions and the space group on ligand-protein complexes: a case study on bacterial phosphotriesterase

A bacterial phosphotriesterase was employed as an experimental paradigm to examine the effects of multiple factors, such as the molecular constructs, the ligands used during protein expression and purification, the crystallization conditions and the space group, on the visualization of molecular complexes of ligands with a target enzyme. In this case, the ligands used were organophosphates that are fragments of the nerve agents and insecticides on which the enzyme acts as a bioscavenger. 12 crystal structures of various phosphotriesterase constructs obtained by directed evolution were analyzed, with resolutions of up to 1.38 Å. Both apo forms and holo forms, complexed with the organophosphate ligands, were studied. Crystals obtained from three different crystallization conditions, crystallized in four space groups, with and without N-terminal tags, were utilized to investigate the impact of these factors on visualizing the organophosphate complexes of the enzyme. The study revealed that the tags used for protein expression can lodge in the active site and hinder ligand binding. Furthermore, the space group in which the protein crystallizes can significantly impact the visualization of bound ligands. It was also observed that the crystallization precipitants can compete with, and even preclude, ligand binding, leading to false positives or to the incorrect identification of lead drug candidates. One of the co-crystallization conditions enabled the definition of the spaces that accommodate the substituents attached to the P atom of several products of organophosphate substrates after detachment of the leaving group. The crystal structures of the complexes of phosphotriesterase with the organophosphate products reveal similar short interaction distances of the two partially charged O atoms of the P-O bonds with the exposed β-Zn2+ ion and the buried α-Zn2+ ion. This suggests that both Zn2+ ions have a role in stabilizing the transition state for substrate hydrolysis. Overall, this study provides valuable insights into the challenges and considerations involved in studying the crystal structures of ligand-protein complexes, highlighting the importance of careful experimental design and rigorous data analysis in ensuring the accuracy and reliability of the resulting phosphotriesterase-organophosphate structures.

Study on the active ingredients of Shenghui decoction inhibiting acetylcholinesterase based on molecular docking and molecular dynamics simulation

We aim to investigate the mechanism and effective components of Shenghui decoction (SHD), which has been shown to inhibit acetylcholinesterase (AChE) through molecular docking (MD) and molecular dynamics simulation (MDS). The effective ingredients in SHD were collected through the TCMSP database and literature review. All components were docked with AChE using CDOCKER. Receptor ligand interaction analysis was performed for the optimal ligand. Two simulation models (model I and II) containing AChE and acetylcholine (ACh) were constructed, in which model II contained the best-docked ligand. Perform 90ns MDS on 2 models. After the simulation, the distance between ACh and AChE peripheral active sites were calculated in both models. The root mean square deviation (RMSD) curve of ligand and receptor, the radius of gyration (Rog) of the receptor, the distance between ligand center and binding site center, and the binding energy of ligand and receptor were calculated in model II. 98 ingredients of SHD were collected, and the best ligand was Tumulosic acid. The residues that form conventional hydrogen bonds between AChE and Tumulosic acid include Tyr132 and Glu201. MDS revealed that ACh could bind to AChE active site in model I. In model II, ACh cannot bind to the binding cavity because the ligand occupies the active site. The RMSD of AChE and Tumulosic acid tends to be stable, the Rog curve of AChE is relatively stable, and the distance between ligand and binding cavity does not fluctuate greatly, indicating that the structure of receptor and ligand is relatively stable. The binding energy of AChE and Tumulosic acid was -24.14 ± 2.46 kcal/mol. SHD contains many effective ingredients that may inhibit AChE activity. Tumulosic acid can occupy the binding site to prevent ACh from entering the chemical domain, thus exerting AChE inhibitory effect.

Conjugates of Methylene Blue with Cycloalkaneindoles as New Multifunctional Agents for Potential Treatment of Neurodegenerative Disease

The development of multi-target-directed ligands (MTDLs) would provide effective therapy of neurodegenerative diseases (ND) with complex and nonclear pathogenesis. A promising method to create such potential drugs is combining neuroactive pharmacophoric groups acting on different biotargets involved in the pathogenesis of ND. We developed a synthetic algorithm for the conjugation of indole derivatives and methylene blue (MB), which are pharmacophoric ligands that act on the key stages of pathogenesis. We synthesized hybrid structures and performed a comprehensive screening for a specific set of biotargets participating in the pathogenesis of ND (i.e., cholinesterases, NMDA receptor, mitochondria, and microtubules assembly). The results of the screening study enabled us to find two lead compounds (4h and 4i) which effectively inhibited cholinesterases and bound to the AChE PAS, possessed antioxidant activity, and stimulated the assembly of microtubules. One of them (4i) exhibited activity as a ligand for the ifenprodil-specific site of the NMDA receptor. In addition, this lead compound was able to bypass the inhibition of complex I and prevent calcium-induced mitochondrial depolarization, suggesting a neuroprotective property that was confirmed using a cellular calcium overload model of neurodegeneration. Thus, these new MB-cycloalkaneindole conjugates constitute a promising class of compounds for the development of multitarget neuroprotective drugs which simultaneously act on several targets, thereby providing cognitive stimulating, neuroprotective, and disease-modifying effects.

Structure of Klebsiella pneumoniae adenosine monophosphate nucleosidase

Klebsiella pneumoniae is a bacterial pathogen that is increasingly responsible for hospital-acquired pneumonia and sepsis. Progressive development of antibiotic resistance has led to higher mortality rates and creates a need for novel treatments. Because of the essential role that nucleotides play in many bacterial processes, enzymes involved in purine and pyrimidine metabolism and transport are ideal targets for the development of novel antibiotics. Herein we describe the structure of K. pneumoniae adenosine monophosphate nucleosidase (KpAmn), a purine salvage enzyme unique to bacteria, as determined by cryoelectron microscopy. The data detail a well conserved fold with a hexameric overall structure and clear density for the putative active site residues. Comparison to the crystal structures of homologous prokaryotic proteins confirms the presence of many of the conserved structural features of this protein yet reveals differences in distal loops in the absence of crystal contacts. This first cryo-EM structure of an Amn enzyme provides a basis for future structure-guided drug development and extends the accuracy of structural characterization of this family of proteins beyond this clinically relevant organism.

Characterizing the potentially neuronal acetylcholinesterase reactivity toward chiral pyraclofos: Enantioselective insights from spectroscopy, in silico docking, molecular dynamics simulation and per-residue energy decomposition studies

Chiral organophosphorus agents are distributed ubiquitously in the environment, but the neuroactivity of these asymmetric chemicals to humans remains uncertain. This scenario was to explore the stereoselective neurobiological response of human acetylcholinesterase (AChE) to chiral pyraclofos at the enantiomeric scale, and then decipher the microscopic basis of enantioselective neurotoxicity of pyraclofos enantiomers. The results indicated that (R)-/(S)-pyraclofos can form the bioconjugates with AChE with a stoichiometric ratio of 1:1, but the neuronal affinity of (R)-pyraclofos (K = 6.31 × 10⁴ M^-1) with AChE was larger than that of (S)-pyraclofos (K = 1.86 × 10⁴ M^-1), and significant enantioselectivity was existed in the biochemical reaction. The modes of neurobiological action revealed that pyraclofos enantiomers were situated at the substrate binding domain, and the strength of the overall noncovalent bonds between (S)-pyraclofos and the residues was weaker than that of (R)-pyraclofos, resulting in the high inhibitory effect of (R)-pyraclofos toward the activity of AChE. Dynamic enantioselective biointeractions illustrated that the intervention of inherent conformational flexibility in the AChE-(R)-pyraclofos was greater than that of the AChE-(S)-pyraclofos, which arises from the big spatial displacement and the conformational flip of the binding domain composed of the residues Thr-64～Asn-89, Gly-122～Asp-134, and Thr-436～Tyr-449. Energy decomposition exhibited that the Gibbs free energies of the AChE-(R)-/(S)-pyraclofos were ΔG° = －37.4/－30.2 kJ mol^-1, respectively, and the disparity comes from the electrostatic energy during the stereoselective neurochemical reactions. Quantitative conformational analysis further confirmed the atomic-scale computational chemistry conclusions, and the perturbation of (S)-pyraclofos on the AChE's ordered conformation was lower than that of (R)-pyraclofos, which is germane to the interaction energies of the crucial residues, e.g. Tyr-124, Tyr-337, Asp-74, Trp-86, and Tyr-119. Evidently, this attempt will contribute mechanistic information to uncovering the neurobiological effects of chiral organophosphates on the body.

α-tocopherol, a slow-binding inhibitor of acetylcholinesterase

Acetylcholinesterase (AChE) is reversibly inhibited by α-tocopherol (α-T). Steady state kinetic analysis shows that α-T is a mixed slow-binding inhibitor of type A of human enzyme (Kci = 0.49 μM; Kui = 1.6 μM) with a residence time of 2 min on target. Molecular dynamics (MD) simulations support this mechanism, and indicate that α-T first forms multiple non-specific interactions with AChE surface near the gorge entrance, then binds to the peripheral side with alkylene chain slowly sliding down the gorge, inducing no significant conformational change. α-T slightly modulates the progressive inhibition of AChE by the cyclic organophosphorus, cresyl saligenylphosphate, accelerating the fast pseudo-first order process of phosphorylation. A moderate accelerating effect of α-T on phosphorylation by paraoxon was also observed after pre-incubation of AChE in the presence of α-T. This accelerating effect of α-T on ex vivo paraoxon-induced diaphragm muscle weakness was also observed. The effect of α-T on AChE phosphylation was interpreted in light of molecular modeling results. From all results it is clear that α-T does not protect AChE against phosphylation by organophosphorus.

Correcting cis-trans-transgressions in macromolecular structure models

Many macromolecular X-ray and cryo-EM structure models deposited in the PDB contain biologically relevant small molecule ligands with unsaturated fatty acid acyl chains, whose cis-trans stereochemistry is incorrect. The molecules are either not properly defined in their stereochemical restraint files, or the proper stereochemistry is neglected during model building. Often, the same molecules appear in deposited models in both isomeric configurations, one of which is almost always incorrect, and the use of the same moiety (HET) identifier and restraint files in model refinement is wrong. We present case studies of frequently occurring molecules and a compilation of identified cases of C-C=C-C cis-trans geometry in the deposited structure models. Full listings of cis/trans torsion angles are provided for models with commonly occurring molecules to assist identification and correction of cis-trans errors and prevent inadvertent use of incorrect models. Caveats for users, advice for modellers and suggestions for remediation efforts with a simple but effective restraint file modification are provided.

Small Glycols Discover Cryptic Pockets on Proteins for Fragment-Based Approaches

Cryptic pockets are visible in ligand-bound protein structures but are occluded in unbound structures. Utilizing these pockets in fragment-based drug-design provides an attractive option for proteins not tractable by classical binding sites. However, owing to their hidden nature, they are difficult to identify. Here, we show that small glycols find cryptic pockets on a diverse set of proteins. Initial crystallography experiments serendipitously revealed the ability of ethylene glycol, a small glycol, to identify a cryptic pocket on the W6A mutant of the RBSX protein (RBSX-W6A). Explicit-solvent molecular dynamics (MD) simulations of RBSX-W6A with the exposed state of the cryptic pocket (ethylene glycol removed) revealed closure of the pocket reiterating that the exposed state of cryptic pockets in general are unstable in the absence of ligands. Also, no change in the pocket was observed for simulations of RBSX-W6A with the occluded state of the cryptic pocket, suggesting that water molecules are not able to open the cryptic pocket. "Cryptic-pocket finding" potential of small glycols was then supported and generalized through additional crystallography experiments, explicit-cosolvent MD simulations, and protein data set construction and analysis. The cryptic pocket on RBSX-W6A was found again upon repeating the crystallography experiments with another small glycol, propylene glycol. Use of ethylene glycol as a probe molecule in cosolvent MD simulations led to the enhanced sampling of the exposed state of experimentally observed cryptic sites on a test set of two proteins (Niemann-Pick C2, Interleukin-2). Further, analyses of protein structures with validated cryptic sites showed that ethylene glycol molecules bind to sites on proteins (Bcl-xL, G-actin, myosin II, and glutamate receptor 2), which become apparent upon binding of biologically relevant ligands. Our study thus suggests potential application of the small glycols in experimental and computational fragment-based approaches to identify cryptic pockets in apparently undruggable and/or difficult targets, making these proteins amenable to drug-design strategies.

Estimating the ligand-binding affinity via λ-dependent umbrella sampling simulations

The umbrella sampling (US) approach has been demonstrated to be a very efficient method for estimating the ligand-binding affinity. However, most of the calculated values overestimate experimental ones that are probably caused by the inaccurate representation of the interaction between the ligand and the surrounding molecules. The issue can be resolved via the implementation aspects of λ-alteration simulation into the US approach, which we call the λ-dependent umbrella sampling (λUS) scheme. In particular, the electrostatic and van der Waals interactions were simultaneously changed by using the coupling parameter λ during λUS simulations. The mean value of obtained results, ∆ G US λ = 0.20 = - 11.59 ± 1.51 kcal mol^-1 , is in good fitting to the mean value of respective experiments, ∆G_EXP =  - 11.26 ± 0.89 kcal mol^-1 . Moreover, the correlation between the proposed approach and experiment is quite good with a value of R US λ = 0.20 = 0.82 ± 0.10 . The λUS scheme significantly enhances the calculated accuracy since the RMSE of the proposed scheme is smaller than traditional US simulations, RMSE US λ = 0.20 = 2.99 ± 0.82 kcal mol^-1 versus RMSE US λ = 0.00 = 5.48 ± 0.81 kcal mol^-1 . Furthermore, the precision is increased since the computed error via λUS approach, δ US λ = 0.20 = 1.51 kcal mol^-1 , was smaller than those of the US simulation, δ US λ = 0.00 = 1.78 kcal mol^-1 . Overall, the proposed approach perhaps provides an efficient way to accurately and precisely estimate the ligand-binding free energy.

Controlled dehydration, structural flexibility and gadolinium MRI contrast compound binding in the human plasma glycoprotein afamin

Afamin, which is a human blood plasma glycoprotein, a putative multifunctional transporter of hydrophobic molecules and a marker for metabolic syndrome, poses multiple challenges for crystallographic structure determination, both practically and in analysis of the models. Several hundred crystals were analysed, and an unusual variability in cell volume and difficulty in solving the structure despite an ∼34% sequence identity with nonglycosylated human serum albumin indicated that the molecule exhibits variable and context-sensitive packing, despite the simplified glycosylation in insect cell-expressed recombinant afamin. Controlled dehydration of the crystals was able to stabilize the orthorhombic crystal form, reducing the number of molecules in the asymmetric unit from the monoclinic form and changing the conformational state of the protein. An iterative strategy using fully automatic experiments available on MASSIF-1 was used to quickly determine the optimal protocol to achieve the phase transition, which should be readily applicable to many types of sample. The study also highlights the drawback of using a single crystallographic structure model for computational modelling purposes given that the conformational state of the binding sites and the electron density in the binding site, which is likely to result from PEGs, greatly varies between models. This also holds for the analysis of nonspecific low-affinity ligands, where often a variety of fragments with similar uncertainty can be modelled, inviting interpretative bias. As a promiscuous transporter, afamin also seems to bind gadoteridol, a magnetic resonance imaging contrast compound, in at least two sites. One pair of gadoteridol molecules is located near the human albumin Sudlow site, and a second gadoteridol molecule is located at an intermolecular site in proximity to domain IA. The data from the co-crystals support modern metrics of data quality in the context of the information that can be gleaned from data sets that would be abandoned on classical measures.

A new crystal form of human acetylcholinesterase for exploratory room-temperature crystallography studies

Structure-guided design of novel pharmacologically active molecules relies at least in part on functionally relevant accuracy of macromolecular structures for template based drug design. Currently, about 95% of all macromolecular X-ray structures available in the PDB (Protein Data Bank) were obtained from diffraction experiments at low, cryogenic temperatures. However, it is known that functionally relevant conformations of both macromolecules and pharmacological ligands can differ at higher, physiological temperatures. We describe in this article development and properties of new human acetylcholinesterase (AChE) crystals of space group P3₁ and a new unit cell, amenable for room-temperature X-ray diffraction studies. We co-crystallized hAChE in P3₁ unit cell with the reversible inhibitor 9-aminoacridine that binds at the base of the active center gorge in addition to inhibitors that span the full length of the gorge, donepezil (Aricept, E2020) and AChE specific inhibitor BW284c51. Their new low temperature P3₁ space group structures appear similar to those previously obtained in the different P3₁21 unit cell. Successful solution of the new room temperature 3.2 Å resolution structure of BW284c51*hAChE complex from large P3₁ crystals enables us to proceed with studying room temperature structures of lower affinity complexes, such as oxime reactivators bound to hAChE, where temperature-related conformational diversity could be expected in both oxime and hAChE, which could lead to better informed structure-based design under conditions approaching physiological temperature.

Productive reorientation of a bound oxime reactivator revealed in room temperature X-ray structures of native and VX-inhibited human acetylcholinesterase

Exposure to organophosphorus compounds (OPs) may be fatal if untreated, and a clear and present danger posed by nerve agent OPs has become palpable in recent years. OPs inactivate acetylcholinesterase (AChE) by covalently modifying its catalytic serine. Inhibited AChE cannot hydrolyze the neurotransmitter acetylcholine leading to its build-up at the cholinergic synapses and creating an acute cholinergic crisis. Current antidotes, including oxime reactivators that attack the OP-AChE conjugate to free the active enzyme, are inefficient. Better reactivators are sought, but their design is hampered by a conformationally rigid portrait of AChE extracted exclusively from 100K X-ray crystallography and scarcity of structural knowledge on human AChE (hAChE). Here, we present room temperature X-ray structures of native and VX-phosphonylated hAChE with an imidazole-based oxime reactivator, RS-170B. We discovered that inhibition with VX triggers substantial conformational changes in bound RS-170B from a "nonproductive" pose (the reactive aldoxime group points away from the VX-bound serine) in the reactivator-only complex to a "semi-productive" orientation in the VX-modified complex. This observation, supported by concurrent molecular simulations, suggested that the narrow active-site gorge of hAChE may be significantly more dynamic than previously thought, allowing RS-170B to reorient inside the gorge. Furthermore, we found that small molecules can bind in the choline-binding site hindering approach to the phosphorous of VX-bound serine. Our results provide structural and mechanistic perspectives on the reactivation of OP-inhibited hAChE and demonstrate that structural studies at physiologically relevant temperatures can deliver previously overlooked insights applicable for designing next-generation antidotes.

Generation of Reactive Oxygen Species by Photosensitizers and their Modes of Action on Proteins

In this review, we first survey the mechanisms underlying the chemical modification of amino acid residues in proteins by singlet oxygen elicited by photosensitizers. Singlet oxygen has the capacity to cause widespread chemical damage to cellular proteins. Its use in photodynamic therapy of tumors thus requires the development of methodologies for specific addressing of the photosensitizer to malignant cells while sparing normal tissue. We describe three targeting paradigms for achieving this objective. The first involves the use of a photosensitizer with a high affinity for its target protein; in this case, the photosensitizer is methylene blue for acetylcholinesterase. The second paradigm involves the use of the hydrophobic photosensitizer hypericin, which has the capacity to interact selectively with partially unfolded forms of proteins, including nascent species in rapidly dividing or virus-infected and cancer cells, acting preferentially at membrane interfaces. In this case, partially unfolded molten globule species of acetylcholinesterase serve as the model system. In the third paradigm, the photodynamic approach takes advantage of a general approach in 'state-of-the-art' chemotherapy, by coupling the photosensitizer emodin to a specific peptide hormone, GnRH, which recognizes malignant cells via specific GnRH receptors on their surface.

Kinetic and structural studies on the interactions of Torpedo californica acetylcholinesterase with two donepezil-like rigid analogues

Acetylcholinesterase inhibitors were introduced for the symptomatic treatment of Alzheimer’s disease (AD). Among the currently approved inhibitors, donepezil (DNP) is one of the most preferred choices in AD therapy. The X-ray crystal structures of Torpedo californica AChE in complex with two novel rigid DNP-like analogs, compounds 1 and 2, have been determined. Kinetic studies indicated that compounds 1 and 2 show a mixed-type inhibition against TcAChE, with K_i values of 11.12 ± 2.88 and 29.86 ± 1.12 nM, respectively. The DNP rigidification results in a likely entropy-enthalpy compensation with solvation effects contributing primarily to AChE binding affinity. Molecular docking evidenced the molecular basis for the binding of compounds 1 and 2 to the active site of β-secretase-1. Overall, these simplified DNP derivatives may represent new structural templates for the design of lead compounds for a more effective therapeutic strategy against AD by foreseeing a dual AChE and BACE-1 inhibitory activity.

Detect, correct, retract: How to manage incorrect structural models

The massive technical and computational progress of biomolecular crystallography has generated some adverse side effects. Most crystal structure models, produced by crystallographers or well-trained structural biologists, constitute useful sources of information, but occasional extreme outliers remind us that the process of structure determination is not fail-safe. The occurrence of severe errors or gross misinterpretations raises fundamental questions: Why do such aberrations emerge in the first place? How did they evade the sophisticated validation procedures which often produce clear and dire warnings, and why were severe errors not noticed by the depositors themselves, their supervisors, referees and editors? Once detected, what can be done to either correct, improve or eliminate such models? How do incorrect models affect the underlying claims or biomedical hypotheses they were intended, but failed, to support? What is the long-range effect of the propagation of such errors? And finally, what mechanisms can be envisioned to restore the validity of the scientific record and, if necessary, retract publications that are clearly invalidated by the lack of experimental evidence? We suggest that cognitive bias and flawed epistemology are likely at the root of the problem. By using examples from the published literature and from public repositories such as the Protein Data Bank, we provide case summaries to guide correction or improvement of structural models. When strong claims are unsustainable because of a deficient crystallographic model, removal of such a model and even retraction of the affected publication are necessary to restore the integrity of the scientific record.

Understanding inhibitor resistance in Mps1 kinase through novel biophysical assays and structures

Monopolar spindle 1 (Mps1/TTK) is a protein kinase essential in mitotic checkpoint signaling, preventing anaphase until all chromosomes are properly attached to spindle microtubules. Mps1 has emerged as a potential target for cancer therapy, and a variety of compounds have been developed to inhibit its kinase activity. Mutations in the catalytic domain of Mps1 that give rise to inhibitor resistance, but retain catalytic activity and do not display cross-resistance to other Mps1 inhibitors, have been described. Here we characterize the interactions of two such mutants, Mps1 C604Y and C604W, which raise resistance to two closely related compounds, NMS-P715 and its derivative Cpd-5, but not to the well characterized Mps1 inhibitor, reversine. We show that estimates of the IC₅₀ (employing a novel specific and efficient assay that utilizes a fluorescently labeled substrate) and the binding affinity (K_D ) indicate that, in both mutants, Cpd-5 should be better tolerated than the closely related NMS-P715. To gain further insight, we determined the crystal structure of the Mps1 kinase mutants bound to Cpd-5 and NMS-P715 and compared the binding modes of Cpd-5, NMS-P715, and reversine. The difference in steric hindrance between Tyr/Trp⁶⁰⁴ and the trifluoromethoxy moiety of NMS-P715, the methoxy moiety of Cpd-5, and complete absence of such a group in reversine, account for differences we observe in vitro Our analysis enforces the notion that inhibitors targeting Mps1 drug-resistant mutations can emerge as a feasible intervention strategy based on existing scaffolds, if the clinical need arises.

Recent developments in structural studies on acetylcholinesterase

This review focuses on several recent developments concerning structure-function relationships in vertebrate acetylcholinesterase. These include studies on high-resolution structures of human acetylcholinesterase and its complexes; the first crystal structure of a snake venom acetylcholinesterase, in which open and closed states of the 'back door' are visualized; a powerful algorithm for redesigning proteins for enhanced expression in prokaryotic systems, as applied to human acetylcholinesterase, which has hitherto been an intractable target; in situ implementation of 'click chemistry' in crystalline acetylcholinesterase, which yields novel insights into the steric and dynamic changes involved in the reaction within the active-site gorge; and a study that demonstrates the effect of crystallization conditions on ligand alignment within a protein complex, in this case the methylene blue-Torpedo californica acetylcholinesterase complex, which highlights the relevance of the precipitant employed to structure-based drug design. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

Molecular dynamics simulations reveal structural insights into inhibitor binding modes and functionality in human Group IIA phospholipase A2

Human Group IIA phospholipase A₂ (hGIIA) promotes inflammation in immune-mediated pathologies by regulating the arachidonic acid pathway through both catalysis-dependent and -independent mechanisms. The hGIIA crystal structure, both alone and inhibitor-bound, together with structures of closely related snake-venom-derived secreted phospholipase enzymes has been well described. However, differentiation of biological and nonbiological contacts and the relevance of structures determined from snake venom enzymes to human enzymes are not clear. We employed molecular dynamics (MD) and docking approaches to understand the binding of inhibitors that selectively or nonselectively block the catalysis-independent mechanism of hGIIA. Our results indicate that hGIIA behaves as a monomer in the solution environment rather than a dimer arrangement that is in the asymmetric unit of some crystal structures. The binding mode of a nonselective inhibitor, KH064, was validated by a combination of the experimental electron density and MD simulations. The binding mode of the selective pentapeptide inhibitor FLSYK to hGIIA was stipulated to be different to that of the snake venom phospholipases A₂ of Daboia russelli pulchella (svPLA₂ ). Our data suggest that the application of MD approaches to crystal structure data is beneficial in evaluating the robustness of conclusions drawn based on crystal structure data alone. Proteins 2017; 85:827-842. © 2016 Wiley Periodicals, Inc.