Structural basis for ligand recognition in a mushroom lectin: solvent structure as specificity predictor

New biochemical insights into the dynamics of water molecules at the GMP or IMP binding site of human GMPR enzyme: A molecular dynamics study

Human GMP reductase (hGMPR) enzyme is involved in a cellular metabolic pathway, converting GMP into IMP, and also it is an important target for anti-leukemic agents. Present computational investigations explain dynamical behavior of water molecules during the conformational transition process from GMP to IMP using molecular dynamics simulations. Residues at substrate-binding site of cancerous protein (PDB Id. 2C6Q) are mostly more dynamic in nature than the normal protein (PDB Id. 2BLE). Nineteen conserved water molecules are identified at the GMP/IMP binding site and are classified as (i) conserved stable dynamic and (ii) infrequent dynamic. Water molecules W11, W14, and W16 are classified as conserved stable dynamic due to their immobile character, whereas remaining water molecules (W1, W2, W3, W4, W5, W7, W8, W9, W10, W12, W13, W15, W17, W18, and W19) are infrequent with dynamic nature. Entrance or displacement of these infrequent water molecules at GMP/IMP sites may occur due to forward and backward movement of reference residues involving ligands. Four water molecules of hGMPR-I and nine water molecules of hGMPR-II are observed in repetitive transitions from GMP to IMP pathway, which indicates discrimination between two isoforms of hGMPRs. Water molecules in cancerous protein are more dynamic and unstable compared to normal protein. These water molecules execute rare dynamical events at GMP binding site and could assist in detailed understanding of conformational transitions that influence the hGMPR's biological functionality. The present study should be of interest to the experimental community engaged in leukemia research and drug discovery for CML cancer.

Investigation of structural analogs of hydroxychloroquine for SARS-CoV-2 main protease (Mpro): A computational drug discovery study

The main protease (Mpro) is the key enzyme of nCOVID-19 and plays a decisive role that makes it an attractive drug target. Multiple analysis of crystal structures reveals the presence of W1, W2, and W3 water locations in the active site pocket of Mpro; W1 and W2 are unstable and are weakly bonded with protein in comparison to W3 of Mpro-native. So, we adopt the water displacement method to occupy W1 or W2 sites by triggering HCQ or its analogs to inactivate the enzyme. Virtual screening is employed to find out best analogs of HCQ, molecular docking is used for water displacement from catalytic region of Mpro, and finally, MD simulations are conducted for validation of these findings. The docking study reveals that W1 and W2 are occupied by respective atoms of ZINC28706440 whereas W2 by HCQ and indacaterol. Finally, MD results demonstrate (i) HCQ occupies W1 and W2 positions, but its analogs (indacaterol and ZINC28706440) are inadequate to retain either W1 or W2 (ii) His41 and Asp187 are stabilized by W3 in Mpro-native and His41, Cys145 and HCQ by W7 in ZINC28706440, and W4, W5, and W6 make water mediated bridge between indacaterol with His41. The structural, dynamical, and thermodynamic (WFP and J value) profiling parameters suggest that W3, W4, and W7 are prominent in their corresponding positions in comparison with W5 and W6. The final results conclude that ZINC28706440 may act as a best analog of HCQ with acceptable physico-chemical and toxicological scores and may further be synthesized for experimental validation.

Structural and Dynamical Impact of Water Molecules at Substrate- or Product-Binding Sites in Human GMPR Enzyme: A Study by Molecular Dynamics Simulations

Human guanosine monophosphate reductase (hGMPR) enzyme maintains the intracellular balance between adenine and guanine nucleotide pools, and it is an excellent target for the design of isoform-specific antileukemic agents. In the present study, we have investigated solvation properties of substrate GMP or product inosine-5'-monophosphate (IMP)-binding pocket of hGMPR by employing molecular dynamics simulations on conformations A (substrate GMP), B [substrate GMP with cofactor nicotinamide adenine dinucleotide phosphate (NDP)], C (product IMP with cofactor NDP), and D (product IMP). Nineteen water sites are identified precisely; they are responsible for the catalytic activity of this site, control structural and dynamical integrity, and electronic consequences of GMP or IMP in the binding site of hGMPR. The water sites of category-1 (W1, W4, W5, W6, W13, and W15) in normal protein and category-2 (W2, W3, W7, W8, W10, W17, and W18) in cancerous protein are unique and stabilize the guanosine or inosine group of GMP or IMP for participation in the enzymatic reaction, whereas the remaining water centers either stabilize pentose sugar ribose or the phosphate group of GMP or IMP. Furthermore, water sites of category-4 (W11, W14, and W16) appear to be conserved in all conformations during the entire simulation. The GMP-binding site in cancerous protein 2C6Q is significantly expanded, and its dynamics are very different from normal protein 2BLE. Furthermore, unique interactions of GMP(N1)···W2···Asp129/Asn158, IMP(N1)···W3···Glu289, and IMP(O6)···W10···Ser270 might be used in a water mimic drug design for hGMPR-II. In this context, water finding probability, relative interaction energy (J) associated with water site W, entropy, and topologies of these three water sites are thermodynamically acceptable for the water displacement method by the modified ligand. Hence, their positions in the catalytic pocket may also facilitate future drug discovery for chronic myelogenous leukemia by the design of appropriately oriented chemical groups that may displace these water molecules to mimic their structural, electronic, and thermodynamic properties.

A C-type lectin (OmCTL) in Onychostoma macrolepis: Binding ability to LPS, PGN and agglutinating activity against bacteria

C-type lectins (CTLs) are calcium-dependent carbohydrate-binding proteins that mainly bind to carbohydrate-based or other ligands to mediate cell adhesion, recognize pathogens, and play important roles in the immune system. In the present study, a novel C-type lectin (OmCTL) isolated from Onychostoma macrolepis was investigated. The open reading frame of OmCTL comprises 468 bp, encoding a 155 amino acid polypeptide with an 18 amino acid putative signaling peptide. The predicted primary OmCTL structure contains a signal peptide, a single carbohydrate recognition domain (CRD) and an EPN/WND motif required for carbohydrate-binding specificity. Using tissue expression pattern analysis, OmCTL has been shownto be highly expressed in the liver, and is also detected in other tissues. OmCTL was significantly upregulated in the liver and spleen following infection with Aeromonas hydrophila, suggesting its involvement in immune response. The recombinant OmCTL protein (rOmCTL) agglutinated two gram-negative bacteria, Escherichia coli and A. hydrophila, in vitro in the presence of Ca²⁺, showing that it is a typical Ca²⁺-dependent carbohydrate-binding protein.Furthermore, rOmCTL purified from E. coli BL21 (DE3) strongly bound to LPS and PGN, as well as all tested bacteria in a Ca²⁺-independent manner. These results indicate that OmCTL plays a central role in the innate immune response and as a pattern recognition receptor that recognizes diverse pathogens among O. macrolepis.

Characterization and expression of galectin-3 in grass carp (Ctenopharyngodon idella)

Galectins are members of evolutionary conserved lectin family and play important roles in the innate and adaptive immunity of both vertebrates and invertebrates. Galectin-3 is the only chimera galectin with one C-terminal carbohydrate recognition domain (CRD) connected to the N-terminal end. Here, a galectin-3 (named CiGal3) from grass carp was identified and characterized, which encodes polypeptides 362 amino acids with a predicted molecular mass of 36.45 kDa and theoretical isoelectric point of 4.91. The sugar binding motifs involved in carbohydrate binding activity (H-N-R, V-N and W--E-R) were detected in CRD. In comparison to other species, CiGal3 showed the highest similarity and identity to Cyprinus carpio (95.3% sequence similarity and 92.5% sequence identity). The subcellular localization of CiGal3 was distributed in the cytoplasm and nucleus of transfected cells. The CiGal3 transcripts were ubiquitously expressed in all checked tissues and highly expressed in immune tissues. In addition, the expression of CiGal3 in liver and spleen was induced post grass carp reovirus (GCRV), lipopolysaccharide (LPS), and polyinosinic:polycytidylic acid (poly I:C) challenge. These results suggest that CiGal3 plays a vital role in the immune system.

Seaweed Secondary Metabolites with Beneficial Health Effects: An Overview of Successes in In Vivo Studies and Clinical Trials

Macroalgae are increasingly viewed as a source of secondary metabolites with great potential for the development of new drugs. In this development, in vitro studies are only the first step in a long process, while in vivo studies and clinical trials are the most revealing stages of the true potential and limitations that a given metabolite may have as a new drug. This literature review aims to give a critical overview of the secondary metabolites that reveal the most interesting results in these two steps. Phlorotannins show great pharmaceutical potential in in vivo models and, among the several examples, the anti-dyslipidemia activity of dieckol must be highlighted because it was more effective than lovastatin in an in vivo model. The IRLIIVLMPILMA tridecapeptide that exhibits an in vivo level of activity similar to the hypotensive clinical drug captopril should still be stressed, as well as griffithsin which showed such stunning results over a variety of animal models and which will probably move onto clinical trials soon. Regarding clinical trials, studies with pure algal metabolites are scarce, limited to those carried out with kahalalide F and fucoxanthin. The majority of clinical trials currently aim to ascertain the effect of algae consumption, as extracts or fractions, on obesity and diabetes.

An efficient use of X-ray information, homology modeling, molecular dynamics and knowledge-based docking techniques to predict protein-monosaccharide complexes

Unraveling the structure of lectin-carbohydrate complexes is vital for understanding key biological recognition processes and development of glycomimetic drugs. Molecular Docking application to predict them is challenging due to their low affinity, hydrophilic nature and ligand conformational diversity. In the last decade several strategies, such as the inclusion of glycan conformation specific scoring functions or our developed solvent-site biased method, have improved carbohydrate docking performance but significant challenges remain, in particular, those related to receptor conformational diversity. In the present work we have analyzed conventional and solvent-site biased autodock4 performance concerning receptor conformational diversity as derived from different crystal structures (apo and holo), Molecular Dynamics snapshots and Homology-based models, for 14 different lectin-monosaccharide complexes. Our results show that both conventional and biased docking yield accurate lectin-monosaccharide complexes, starting from either apo or homology-based structures, even when only moderate (45%) sequence identity templates are available. An essential element for success is a proper combination of a middle-sized (10-100 structures) conformational ensemble, derived either from Molecular dynamics or multiple homology model building. Consistent with our previous works, results show that solvent-site biased methods improve overall performance, but that results are still highly system dependent. Finally, our results also show that docking can select the correct receptor structure within the ensemble, underscoring the relevance of joint evaluation of both ligand pose and receptor conformation.

Computational screening of known broad-spectrum antiviral small organic molecules for potential influenza HA stem inhibitors

Background

With the emergence of new influenza virus strains that are resistant to current inhibitors such as oseltamivir (anti-neuraminidase (NA)) and amantadine (anti-M2 proton channel), influenza A viruses continue to be a serious threat to the public health worldwide. With this in view, there is a persistent need for the development of broader and more effective vaccines and therapeutics. Identification of broadly neutralizing antibodies (bNAbs) that recognize relatively invariant structures ‎on influenza haemagglutinin (HA) stem has invigorated efforts to develop universal influenza vaccines.

Aim

The current computational study is designed to identify potential flavonoid inhibitors that bind to the contact epitopes of HA stem that are targeted by broadly neutralizing antibodies (bNAb).

Method

In this study, we utilized the three-dimensional crystallographic structure of different HA subtypes (H1, H2, H5, H3, and H7) in complex with bNAb to screen for potential broadly reactive influenza inhibitors. We performed Quantitative Structure-Activity and Relationship (QSAR) for 100 natural compounds known for their antiviral activity and performed molecular docking using AutoDock 4.2 suite. Furthermore, we conducted virtual screening of 1413 bioassay hit compounds by using virtual lab bench CLC Drug Discovery.

Results

The results showed 18 lead flavonoids with strong binding abilities to bNAb epitopes of various HA subtypes. These 18 broadly reactive compounds exhibited significant interactions with an average of seven Hbonds, docking energy of -22.43 kcal·mol⁻¹, and minimum interaction ‎ energy of -4.65 kcal·mol⁻¹, with functional contact residues. Procyanidin depicted strong interactions with group 1 HAs, whereas both sorbitol and procyanidin exhibited significant interactions with group 2 HAs.

Conclusion

Using in silico docking analysis, we identified 18 bioactive flavonoids with potential strong binding cababilities to influenza HA-stems of various subtypes, which are the target for bNAb. The virtual screened bioassay hit compounds depicted a high number of Hbonds but low interaction and docking values compared to antiviral flavonoids. Using structure-based design and nanotechnology-based approaches, identified molecules could be modified to generate next generation anti-influenza drugs.

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

Synthesis of β-galactosylamides as ligands of the peanut lectin. Insights into the recognition process

The synthesis of mono and divalent β-galactosylamides linked to a hydroxylated chain having a C2 symmetry axis derived from l-tartaric anhydride is reported. Reference compounds devoid of hydroxyl groups in the linker were also prepared from β-galactosylamine and succinic anhydride. After functionalization with an alkynyl residue, the resulting building blocks were grafted onto different azide-equipped scaffolds through the copper catalyzed azide-alkyne cycloaddition. Thus, a family of structurally related mono and divalent β-N-galactopyranosylamides was obtained and fully characterized. The binding affinities of the ligands towards the model lectin PNA were measured by the enzyme-linked lectin assay (ELLA). The IC₅₀ values were significantly higher than that of galactose but the presence of hydroxyl groups in the aglycone chain improved lectin recognition. Docking and molecular dynamics experiments were in accordance with the hypothesis that a hydroxyl group properly disposed in the linker could mimic the Glc O3 in the recognition process. On the other hand, divalent presentation of the ligands led to lectin affinity enhancements.

A galectin from Hyriopsis cumingii involved in the innate immune response against to pathogenic microorganism and its expression profiling during pearl sac formation

Hyriopsis cumingii is the most important freshwater pearl mussel cultured in China. The operation for implantation is one necessary technical step for pearl culture. However, implantation-induced trauma results in a series of immune responses and can enable the invasion of pathogenic microbes. Lectin proteins are found widely in nature and play important roles in innate immunity. Galectins are members of the lectin superfamily and are characterized by one or several carbohydrate recognition domains (CRDs) that produce multiple sugar binding sites on the protein. Here we cloned and characterized the H. cumingii galectin gene HcGal1, which encodes a 312 amino acid galectin protein. The HcGal1 transcript was detected in all tested H. cumingii tissues and showed higher expression specifically in immune tissues. The significant upregulation of HcGal1 expression was observed after challenging the mussel with lipopolysaccharide or Gram-negative and Gram-positive bacteria. After implantation, significant downregulation of the HcGal1 transcript was noted in the mantle, hemocytes, and pearl sac in the acute-stress stage (0-24 h) and the stage of wound healing and pearl-sac formation (24 h-7 d). In addition, significant upregulation of HcGal1 expression was observed in the liver in the stage of wound healing and pearl-sac formation. In the pearl-secretion stage (7-35 d), the HcGal1 transcript levels returned to normal in all tested tissues. We also show that recombinantly expressed and purified HcGal1 can agglutinate some Gram-negative and Gram-positive bacteria. In addition, in vivo experiments showed that the recombinant protein HcGal1 could promote phagocytosis by hemocytes. Our data suggest that HcGal1 plays a role in innate immune responses involved in pathogen recognition and wound healing.

Purification and characterization of a mucin specific mycelial lectin from Aspergillus gorakhpurensis: application for mitogenic and antimicrobial activity

BACKGROUND

Lectins are carbohydrate binding proteins or glycoproteins that bind reversibly to specific carbohydrates present on the apposing cells, which are responsible for their ability to agglutinate red blood cells, lymphocytes, fibroblasts, etc. Interest in lectins has been intensified due to their carbohydrate specificity as they can be valuable reagents for the investigation of cell surface sugars, purification and characterization of glycoproteins. The present study reports the purification, characterization and evaluation of mitogenic and antimicrobial potential of a mycelial lectin from Aspergillus gorakhpurensis.

METHODS

Affinity chromatography on mucin-sepharose column was carried out for purification of Aspergillus gorakhpurensis lectin. The lectin was characterized for physico-chemical parameters. Mitogenic potential of the lectin was evaluated against splenocytes of Swiss albino mice by MTT assay. Antimicrobial activity of the purified lectin has also been evaluated by disc diffusion assay.

RESULTS

Single-step affinity purification resulted in 18.6-fold purification of the mycelial lectin. The molecular mass of the lectin was found to be 70 kDa and it was composed of two subunits of 34.8 kDa as determined by gel filtration chromatography, SDS-PAGE and MALDI-TOF analysis. pH optima of the lectin was found to be 6.5-9.5, while optimum temperature for lectin activity was 20-30 °C. Lectin was stable within a pH range of 7.0-10.5 and showed fair thermostability. EDTA did not affect lectin activity whereas it was found susceptible to the denaturants tested. MTT assay revealed strong mitogenic potential of A. gorakhpurensis lectin at a concentration upto 150 µg/mL. Antimicrobial activity assay showed its potent antibacterial activity against Bacillus cereus, Staphylococcous aureus and Escherichia coli and marginal antifungal activity against Saccharomyces cerevisiae.

CONCLUSION

This is the first report on the mitogenic and antimicrobial potential of Aspergillus gorakhpurensis lectin. The results will provide useful guidelines for further research in clinical applications of this lectin.

Recent advances in employing molecular modelling to determine the specificity of glycan-binding proteins

Impressive improvements in docking performance can be achieved by applying energy bonuses to poses in which glycan hydroxyl groups occupy positions otherwise preferred by bound waters. In addition, inclusion of glycosidic conformational energies allows unlikely glycan conformations to be appropriately penalized. A method for predicting the binding specificity of glycan-binding proteins has been developed, which is based on grafting glycan branches onto a minimal binding determinant in the binding site. Grafting can be used either to screen virtual libraries of glycans, such as the known glycome, or to identify docked poses of minimal binding determinants that are consistent with specificity data. The reviewed advances allow accurate modelling of carbohydrate-protein 3D co-complexes, but challenges remain in ranking the affinity of congeners.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

The allosteric modulation of thyroxine-binding globulin affinity is entropy driven

BACKGROUND

Thyroxine-binding globulin (TBG) is a non-inhibitory member of the serpin family of proteins whose main structural element is the reactive center loop (RCL), that, upon cleavage by proteases, is inserted into the protein core adopting a β-strand conformation (stressed to relaxed transition, S-to-R). After S-to-R transition thyroxine (T4) affinity decreases. However, crystallographic studies in the presence or absence of the hormone in different states are unable to show significant differences in the structure and interactions of the binding site. Experimental results also suggest the existence of several S states (differing in the number of inserted RCL residues), associated with a differential affinity.

METHODS

To shed light into the molecular basis that regulates T4 affinity according to the degree of RCL insertion in TBG, we performed extended molecular dynamics simulations combined with several thermodynamic analysis of the T4 binding to TBG in three different S states, and in the R state.

RESULTS

Our results show that, despite T4 binding in the protein by similar interactions in all states, a good correlation between the degree of RCL insertion and the binding affinity, driven by a change in TBG conformational entropy, was observed.

CONCLUSION

TBG allosteric regulation is entropy driven. The presence of multiple S states may allow more efficient T4 release due to protease activity.

GENERAL SIGNIFICANCE

The presented results are clear examples of how computer simulation methods can reveal the thermodynamic basis of allosteric effects, and provide a general framework for understanding serpin allosteric affinity regulation.

Solvent structure improves docking prediction in lectin-carbohydrate complexes

Recognition and complex formation between proteins and carbohydrates is a key issue in many important biological processes. Determination of the three-dimensional structure of such complexes is thus most relevant, but particularly challenging because of their usually low binding affinity. In silico docking methods have a long-standing tradition in predicting protein-ligand complexes, and allow a potentially fast exploration of a number of possible protein-carbohydrate complex structures. However, determining which of these predicted complexes represents the correct structure is not always straightforward. In this work, we present a modification of the scoring function provided by AutoDock4, a widely used docking software, on the basis of analysis of the solvent structure adjacent to the protein surface, as derived from molecular dynamics simulations, that allows the definition and characterization of regions with higher water occupancy than the bulk solvent, called water sites. They mimic the interaction held between the carbohydrate -OH groups and the protein. We used this information for an improved docking method in relation to its capacity to correctly predict the protein-carbohydrate complexes for a number of tested proteins, whose ligands range in size from mono- to tetrasaccharide. Our results show that the presented method significantly improves the docking predictions. The resulting solvent-structure-biased docking protocol, therefore, appears as a powerful tool for the design and optimization of development of glycomimetic drugs, while providing new insights into protein-carbohydrate interactions. Moreover, the achieved improvement also underscores the relevance of the solvent structure to the protein carbohydrate recognition process.

Combining molecular docking and molecular dynamics to predict the binding modes of flavonoid derivatives with the neuraminidase of the 2009 H1N1 influenza A virus

Control of flavonoid derivatives inhibitors release through the inhibition of neuraminidase has been identified as a potential target for the treatment of H1N1 influenza disease. We have employed molecular dynamics simulation techniques to optimize the 2009 H1N1 influenza neuraminidase X-ray crystal structure. Molecular docking of the compounds revealed the possible binding mode. Our molecular dynamics simulations combined with the solvated interaction energies technique was applied to predict the docking models of the inhibitors in the binding pocket of the H1N1 influenza neuraminidase. In the simulations, the correlation of the predicted and experimental binding free energies of all 20 flavonoid derivatives inhibitors is satisfactory, as indicated by R² = 0.75.