Structure and Energetics of Ligand-Fluorine Interactions with Galectin-3 Backbone and Side-Chain Amides: Insight into Solvation Effects and Multipolar Interactions

Interplay of halogen bonding and solvation in protein-ligand binding

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

Nanopore surface engineering of molecular imprinted mesoporous organosilica for rapid and selective detection of L-thyroxine

To develop a biosensing platform for precise diagnosis and management of thyroid-related diseases, the sensitive and selective recognition and identification of L-thyroxine (T4), a thyroid hormone, remains challenging. We herein introduce T4-imprinted mesoporous organosilica (T4-IMO) for sensitive and specific detection of T4 via the sophisticated engineering of pore surfaces using additives with different polarities. The pore surface of T4-IMO emitting a stable fluorescence signal is simply modified by fixed additives. Additives embedded in the pore surface promote the rebinding response of T4 into the recognized cavities, subsequently sensitizing T4 detection. Notably, T4-IMO containing abundant fluorine elements on the pore surface shows a high affinity toward T4, remarkably boosting the rebinding capacity. In addition to good selectivity to T4, the "turn-off" fluorescent signal exhibits a linear relationship with the logarithm of T4 concentration in a range of 0-500 nM with a detection limit of 0.47 nM in synthetic urine samples. Our findings can establish an insightful strategy for the rational design of molecular-recognition-based sensor systems for the selective and sensitive detection of target analytes.

Ligand Sulfur Oxidation State Progressively Alters Galectin-3-Ligand Complex Conformations To Induce Affinity-Influencing Hydrogen Bonds

Galectins play biological roles in immune regulation and tumor progression. Ligands with high affinity for the shallow, hydrophilic galectin-3 ligand binding site rely primarily on a galactose core with appended aryltriazole moieties, making hydrophobic interactions and π-stacking. We designed and synthesized phenyl sulfone, sulfoxide, and sulfide-triazolyl thiogalactoside derivatives to create affinity-enhancing hydrogen bonds, hydrophobic and π-interactions. Crystal structures and thermodynamic analyses revealed that the sulfoxide and sulfone ligands form hydrogen bonds while retaining π-interactions, resulting in improved affinities and unique binding poses. The sulfoxide, bearing one hydrogen bond acceptor, leads to an affinity decrease compared to the sulfide, whereas the corresponding sulfone forms three hydrogen bonds, two directly with Asn and Arg side chains and one water-mediated to an Asp side chain, respectively, which alters the complex structure and increases affinity. These findings highlight that the sulfur oxidation state influences both the interaction thermodynamics and structure.

Intermolecular n→π* Interactions in Supramolecular Chemistry and Catalysis

The n→π* interactions describing attractive force between lone pairs (lps) of nucleophile and carbonyl or polarized unsaturated bonds have recently attracted growing attentions in various disciplines. So far, such non-covalent driving force are mainly concentrated to intramolecular systems. Intermolecular n→π* interactions in principle could produce fascinated supramolecular systems or facilitate organic reactions, however, they remain largely underexplored due to the very weak energy of individual interaction. This review attempts to give an overview of the challenging intermolecular n→π* interactions, much efforts emphasize the supramolecular systems, catalytic processes and spectroscopic measurements.

New therapeutic approaches against pulmonary fibrosis

Pulmonary fibrosis is the end-stage change of a large class of lung diseases characterized by the proliferation of fibroblasts and the accumulation of a large amount of extracellular matrix, accompanied by inflammatory damage and tissue structure destruction, which also shows the normal alveolar tissue is damaged and then abnormally repaired resulting in structural abnormalities (scarring). Pulmonary fibrosis has a serious impact on the respiratory function of the human body, and the clinical manifestation is progressive dyspnea. The incidence of pulmonary fibrosis-related diseases is increasing year by year, and no curative drugs have appeared so far. Nevertheless, research on pulmonary fibrosis have also increased in recent years, but there are no breakthrough results. Pathological changes of pulmonary fibrosis appear in the lungs of patients with coronavirus disease 2019 (COVID-19) that have not yet ended, and whether to improve the condition of patients with COVID-19 by means of the anti-fibrosis therapy, which are the questions we need to address now. This review systematically sheds light on the current state of research on fibrosis from multiple perspectives, hoping to provide some references for design and optimization of subsequent drugs and the selection of anti-fibrosis treatment plans and strategies.

An Automated Approach to Assess Relative Galectin-Glycan Affinity Following Glycan Microarray Analysis

Numerous studies have highlighted the utility of glycan microarray analysis for the elucidation of protein-glycan interactions. However, most current glycan microarray studies analyze glycan binding protein (GBP)-glycan interactions at a single protein concentration. While this approach provides useful information related to a GBP's overall binding capabilities, extrapolation of true glycan binding preferences using this method fails to account for printing variations or other factors that may confound relative binding. To overcome this limitation, we examined glycan array binding of three galectins over a range of concentrations to allow for a more complete assessment of binding preferences. This approach produced a richer data set than single concentration analysis and provided more accurate identification of true glycan binding preferences. However, while this approach can be highly informative, currently available data analysis approaches make it impractical to perform binding isotherms for each glycan present on currently available platforms following GBP evaluation. To overcome this limitation, we developed a method to directly optimize the efficiency of assessing association constants following multi-GBP concentration glycan array analysis. To this end, we developed programs that automatically analyze raw array data (kdMining) to generate output graphics (kaPlotting) following array analysis at multiple doses. These automatic programing methods reduced processing time from 32.8 h to 1.67 min. Taken together, these results demonstrate an effective approach to glycan array analysis that provides improved detail and efficiency when compared to previous methods.

Engineering the ligand specificity of the human galectin-1 by incorporation of tryptophan analogs

Galectin-1 is a β-galactoside-binding lectin with manifold biological functions. A single tryptophan residue (W68) in its carbohydrate binding site plays a major role in ligand binding and is highly conserved among galectins. To fine tune galectin-1 specificity, we introduced several non-canonical tryptophan analogs at this position of human galectin-1 and analyzed the resulting variants using glycan microarrays. Two variants containing 7-azatryptophan and 7-fluorotryptophan showed a reduced affinity for 3'-sulfated oligosaccharides. Their interaction with different ligands was further analyzed by fluorescence polarization competition assay. Using molecular modeling we provide structural clues that the change in affinities comes from modulated interactions and solvation patterns. Thus, we show that the introduction of subtle atomic mutations in the ligand binding site of galectin-1 is an attractive approach for fine-tuning its interactions with different ligands.

Toward Glycomaterials with Selectivity as Well as Affinity

Multivalent glycosylated materials (polymers, surfaces, and particles) often show high affinity toward carbohydrate binding proteins (e.g., lectins) due to the nonlinear enhancement from the cluster glycoside effect. This affinity gain has potential in applications from diagnostics, biosensors, and targeted delivery to anti-infectives and in an understanding of basic glycobiology. This perspective highlights the question of selectivity, which is less often addressed due to the reductionist nature of glycomaterials and the promiscuity of many lectins. The use of macromolecular features, including architecture, heterogeneous ligand display, and the installation of non-natural glycans, to address this challenge is discussed, and examples of selectivity gains are given.

Molecular mechanism of interspecies differences in the binding affinity of TD139 to Galectin-3

Galectin-3 (Gal-3), a β-galactoside-binding lectin, has been implicated in a plethora of pathological disorders including fibrosis, inflammation, cancer and metabolic diseases. TD139-a thio-digalactoside inhibitor developed by Galecto Biotech as a potential therapeutic for idiopathic pulmonary fibrosis-is the most advanced small-molecule Gal-3 inhibitor in clinical studies. It binds to human Gal-3 with high affinity but has lower affinity towards mouse and rat homologs, which is also manifested in the differential inhibition of Gal-3 function. Using biophysical methods and high-resolution X-ray co-crystal structures of TD139 and Gal-3 proteins, we demonstrate that a single amino acid change corresponding to A146 in human Gal-3 is sufficient for the observed reduction in the binding affinity of TD139 in rodents. Site-directed mutagenesis of A146V (in human Gal-3) and V160A (in mouse Gal-3) was sufficient to interchange the affinities, mainly by affecting the off rates of the inhibitor binding. In addition, molecular dynamics simulations of both wild-type and mutant structures revealed the sustained favorable noncovalent interactions between the fluorophenyl ring and the active site A146 (human Gal-3 and mouse V160A) that corroborate the finding from biophysical studies. Current findings have ramifications in the context of optimization of drug candidates against Gal-3.

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.

Molecular dynamics simulations elucidate oligosaccharide recognition pathways by galectin-3 at atomic resolution

The recognition of carbohydrates by lectins plays key roles in diverse cellular processes such as cellular adhesion, proliferation, and apoptosis, which makes it a therapeutic target of significance against cancers. One of the most functionally active lectins, galectin-3 is distinctively known for its specific binding affinity toward β-galactoside. However, despite the prevalence of high-resolution crystallographic structures, the mechanistic basis and more significantly, the dynamic process underlying carbohydrate recognition by galectin-3 are currently elusive. To this end, we employed extensive Molecular Dynamics simulations to unravel the complete binding event of human galectin-3 with its native natural ligand N-acetyllactosamine (LacNAc) at atomic precision. The simulation trajectory demonstrates that the oligosaccharide diffuses around the protein and eventually identifies and binds to the biologically designated binding site of galectin-3 in real time. The simulated bound pose correlates with the crystallographic pose with atomic-level accuracy and recapitulates the signature stabilizing galectin-3/oligosaccharide interactions. The recognition pathway also reveals a set of transient non-native ligand poses in its course to the receptor. Interestingly, kinetic analysis in combination with a residue-level picture revealed that the key to the efficacy of a more active structural variant of the LacNAc lay in the ligand's resilience against disassociation from galectin-3. By catching the ligand in the act of finding its target, our investigations elucidate the detailed recognition mechanism of the carbohydrate-binding domain of galectin-3 and underscore the importance of ligand-target binary complex residence time in understanding the structure-activity relationship of cognate ligands.

Computational Study of Potential Galectin-3 Inhibitors in the Treatment of COVID-19

Galectin-3 is a carbohydrate-binding protein and the most studied member of the galectin family. It regulates several functions throughout the body, among which are inflammation and post-injury remodelling. Recent studies have highlighted the similarity between Galectin-3's carbohydrate recognition domain and the so-called "galectin fold" present on the N-terminal domain of the S1 sub-unit of the SARS-CoV-2 spike protein. Sialic acids binding to the N-terminal domain of the Spike protein are known to be crucial for viral entry into humans, and the role of Galectin-3 as a mediator of lung fibrosis has long been the object of study since its levels have been found to be abnormally high in alveolar macrophages following lung injury. In this context, the discovery of a double inhibitor may both prevent viral entry and reduce post-infection pulmonary fibrosis. In this study, we use a database of 56 compounds, among which 37 have known experimental affinity with Galectin-3. We carry out virtual screening of this database with respect to Galectin-3 and Spike protein. Several ligands are found to exhibit promising binding affinity and interaction with the Spike protein's N-terminal domain as well as with Galectin-3. This finding strongly suggests that existing Galectin-3 inhibitors possess dual-binding capabilities to disrupt Spike-ACE2 interactions. Herein we identify the most promising inhibitors of Galectin-3 and Spike proteins, of which five emerge as potential dual effective inhibitors. Our preliminary results warrant further in vitro and in vivo testing of these putative inhibitors against SARS-CoV-2 with the hope of being able to halt the spread of the virus in the future.

Exploring ligand dynamics in protein crystal structures with ensemble refinement

Understanding the dynamics of ligands bound to proteins is an important task in medicinal chemistry and drug design. However, the dominant technique for determining protein-ligand structures, X-ray crystallography, does not fully account for dynamics and cannot accurately describe the movements of ligands in protein binding sites. In this article, an alternative method, ensemble refinement, is used on six protein-ligand complexes with the aim of understanding the conformational diversity of ligands in protein crystal structures. The results show that ensemble refinement sometimes indicates that the flexibility of parts of the ligand and some protein side chains is larger than that which can be described by a single conformation and atomic displacement parameters. However, since the electron-density maps are comparable and Rfree values are slightly increased, the original crystal structure is still a better model from a statistical point of view. On the other hand, it is shown that molecular-dynamics simulations and automatic generation of alternative conformations in crystallographic refinement confirm that the flexibility of these groups is larger than is observed in standard refinement. Moreover, the flexible groups in ensemble refinement coincide with groups that give high atomic displacement parameters or non-unity occupancy if optimized in standard refinement. Therefore, the conformational diversity indicated by ensemble refinement seems to be qualitatively correct, indicating that ensemble refinement can be an important complement to standard crystallographic refinement as a tool to discover which parts of crystal structures may show extensive flexibility and therefore are poorly described by a single conformation. However, the diversity of the ensembles is often exaggerated (probably partly owing to the rather poor force field employed) and the ensembles should not be trusted in detail.

Entropy-Entropy Compensation between the Protein, Ligand, and Solvent Degrees of Freedom Fine-Tunes Affinity in Ligand Binding to Galectin-3C

Molecular recognition is fundamental to biological signaling. A central question is how individual interactions between molecular moieties affect the thermodynamics of ligand binding to proteins and how these effects might propagate beyond the immediate neighborhood of the binding site. Here, we investigate this question by introducing minor changes in ligand structure and characterizing the effects of these on ligand affinity to the carbohydrate recognition domain of galectin-3, using a combination of isothermal titration calorimetry, X-ray crystallography, NMR relaxation, and computational approaches including molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). We studied a congeneric series of ligands with a fluorophenyl-triazole moiety, where the fluorine substituent varies between the ortho, meta, and para positions (denoted O, M, and P). The M and P ligands have similar affinities, whereas the O ligand has 3-fold lower affinity, reflecting differences in binding enthalpy and entropy. The results reveal surprising differences in conformational and solvation entropy among the three complexes. NMR backbone order parameters show that the O-bound protein has reduced conformational entropy compared to the M and P complexes. By contrast, the bound ligand is more flexible in the O complex, as determined by ¹⁹F NMR relaxation, ensemble-refined X-ray diffraction data, and MD simulations. Furthermore, GIST calculations indicate that the O-bound complex has less unfavorable solvation entropy compared to the other two complexes. Thus, the results indicate compensatory effects from ligand conformational entropy and water entropy, on the one hand, and protein conformational entropy, on the other hand. Taken together, these different contributions amount to entropy-entropy compensation among the system components involved in ligand binding to a target protein.

Artificial Intelligence-Based Application to Explore Inhibitors of Neurodegenerative Diseases

Neuroinflammation is a common factor in neurodegenerative diseases, and it has been demonstrated that galectin-3 activates microglia and astrocytes, leading to inflammation. This means that inhibition of galectin-3 may become a new strategy for the treatment of neurodegenerative diseases. Based on this motivation, the objective of this study is to explore an integrated new approach for finding lead compounds that inhibit galectin-3, by combining universal artificial intelligence algorithms with traditional drug screening methods. Based on molecular docking method, potential compounds with high binding affinity were screened out from Chinese medicine database. Manifold artificial intelligence algorithms were performed to validate the docking results and further screen compounds. Among all involved predictive methods, the deep learning-based algorithm made 500 modeling attempts, and the square correlation coefficient of the best trained model on the test sets was 0.9. The XGBoost model reached a square correlation coefficient of 0.97 and a mean square error of only 0.01. We switched to the ZINC database and performed the same experiment, the results showed that the compounds in the former database showed stronger affinity. Finally, we further verified through molecular dynamics simulation that the complex composed of the candidate ligand and the target protein showed stable binding within 100 ns of simulation time. In summary, combined with the application based on artificial intelligence algorithms, we unearthed the active ingredients 1,2-Dimethylbenzene and Typhic acid contained in Crataegus pinnatifida and Typha angustata might be the effective inhibitors of neurodegenerative diseases. The high prediction accuracy of the models shows that it has practical application value on small sample data sets such as drug screening.

Introducing affinity and selectivity into galectin-targeting nanoparticles with fluorinated glycan ligands

Galectins are potential biomarkers and therapeutic targets. However, galectins display broad affinity towards β-galactosides meaning glycan-based (nano)biosensors lack the required selectivity and affinity. Using a polymer-stabilized nanoparticle biosensing platform, we herein demonstrate that the specificity of immobilised lacto-N-biose towards galectins can be 'turned on/off' by using site-specific glycan fluorination and in some cases reversal of specificity can be achieved. The panel of fluoro-glycans were obtained by a chemoenzymatic approach, exploiting BiGalK and BiGalHexNAcP enzymes from Bifidobacterium infantis which are shown to tolerate fluorinated glycans, introducing structural diversity which would be very laborious by chemical methods alone. These results demonstrate that integrating non-natural, fluorinated glycans into nanomaterials can encode unprecedented selectivity with potential applications in biosensing.

Addressing the Structural Complexity of Fluorinated Glucose Analogues: Insight into Lipophilicities and Solvation Effects

In this work, we synthesized all mono-, di-, and trifluorinated glucopyranose analogues at positions C-2, C-3, C-4, and C-6. This systematic investigation allowed us to perform direct comparison of ¹⁹ F resonances of fluorinated glucose analogues and also to determine their lipophilicities. Compounds with a fluorine atom at C-6 are usually the most hydrophilic, whereas those with vicinal polyfluorinated motifs are the most lipophilic. Finally, the solvation energies of fluorinated glucose analogues were assessed for the first time by using density functional theory. This method allowed the log P prediction of fluoroglucose analogues, which was comparable to the C log P values obtained from various web-based programs.

Understanding the role of galectin inhibitors as potential candidates for SARS-CoV-2 spike protein: in silico studies

Galectin 3 have the potential to inhibit the SARS-CoV-2 spike protein. We validated the studies by docking, MD and MM/GBSA calculations.

3-Substituted 1-Naphthamidomethyl-C-galactosyls Interact with Two Unique Sub-sites for High-Affinity and High-Selectivity Inhibition of Galectin-3

The galectins are a family of galactose-binding proteins playing key roles in inflammatory processes and cancer. However, they are structurally very closely related, and discovery of highly selective inhibitors is challenging. In this work, we report the design of novel inhibitors binding to a subsite unique to galectin-3, which confers both high selectivity and affinity towards galectin-3. Olefin cross metathesis between allyl β-C-galactopyranosyl and 1-vinylnaphthalenes or acylation of aminomethyl β-C-galactopyranosyl with 1-naphthoic acid derivatives gave C-galactopyranosyls carrying 1-naphthamide structural elements that interacted favorably with a galectin-3 unique subsite according to molecular modeling and X-ray structural analysis of two inhibitor-galectin-3 complexes. Affinities were down to sub-µM and selectivities over galectin-1, 2, 4 N-terminal domain, 4 C-terminal domain, 7, 8 N-terminal domain, 9 N-terminal domain, and 9 C-terminal domain were high. These results show that high affinity and selectivity for a single galectin can be achieved by targeting unique subsites, which holds promise for further development of small and selective galectin inhibitors.

The Sweet-Side of Leukocytes: Galectins as Master Regulators of Neutrophil Function

Among responders to microbial invasion, neutrophils represent one of the earliest and perhaps most important factors that contribute to initial host defense. Effective neutrophil immunity requires their rapid mobilization to the site of infection, which requires efficient extravasation, activation, chemotaxis, phagocytosis, and eventual killing of potential microbial pathogens. Following pathogen elimination, neutrophils must be eliminated to prevent additional host injury and subsequent exacerbation of the inflammatory response. Galectins, expressed in nearly every tissue and regulated by unique sensitivity to oxidative and proteolytic inactivation, appear to influence nearly every aspect of neutrophil function. In this review, we will examine the impact of galectins on neutrophils, with a particular focus on the unique biochemical traits that allow galectin family members to spatially and temporally regulate neutrophil function.

Structure and Energetics of Ligand-Fluorine Interactions with Galectin-3 Backbone and Side-Chain Amides: Insight into Solvation Effects and Multipolar Interactions

Multipolar fluorine–amide interactions with backbone and side‐chain amides have been described as important for protein–ligand interactions and have been used to improve the potency of synthetic inhibitors. In this study, fluorine interactions within a well‐defined binding pocket on galectin‐3 were investigated systematically using phenyltriazolyl‐thiogalactosides fluorinated singly or multiply at various positions on the phenyl ring. X‐ray structures of the C‐terminal domain of galectin‐3 in complex with eight of these ligands revealed potential orthogonal fluorine–amide interactions with backbone amides and one with a side‐chain amide. The two interactions involving main‐chain amides seem to have a strong influence on affinity as determined by fluorescence anisotropy. In contrast, the interaction with the side‐chain amide did not influence affinity. Quantum mechanics calculations were used to analyze the relative contributions of these interactions to the binding energies. No clear correlation could be found between the relative energies of the fluorine–main‐chain amide interactions and the overall binding energy. Instead, dispersion and desolvation effects play a larger role. The results confirm that the contribution of fluorine–amide interactions to protein–ligand interactions cannot simply be predicted, on geometrical considerations alone, but require careful consideration of the energetic components.