Tofacitinib and peficitinib inhibitors of Janus kinase for autoimmune disease treatment: a quantum biochemistry approach.

Autoimmune inflammatory diseases, such as rheumatoid arthritis (RA) and ulcerative colitis, are associated with an uncontrolled production of cytokines leading to the pronounced inflammatory response of these disorders. Their therapy is currently focused on the inhibition of cytokine receptors, such as the Janus kinase (JAK) protein family. Tofacitinib and peficitinib are JAK inhibitors that have been recently approved to treat rheumatoid arthritis. In this study, an in-depth analysis was carried out through quantum biochemistry to understand the interactions involved in the complexes formed by JAK1 and tofacitinib or peficitinib. Computational analyses provided new insights into the binding mechanisms between tofacitinib or peficitinib and JAK1. The essential amino acid residues that support the complex are also identified and reported. Additionally, we report new interactions, such as van der Waals; hydrogen bonds; and alkyl, pi-alkyl, and pi-sulfur forces, that stabilize the complexes. The computational results revealed that peficitinib presents a similar affinity to JAK1 compared to tofacitinib based on their interaction energies.

Unveiling fructose and glucose binding to human serum albumin: fluorescence measurements and docking, molecular dynamics and quantum biochemistry computations.

This research examines the interaction between human serum albumin (HSA) and various sugar forms (ß-D-fructofuranose (FRC), α-D-glucopyranose (GLC), Keto-D-fructose (FRO), Aldehydo-D-glucose (GLO), and modified Aldehydo-D-glucose (GLOm)) using fluorescent spectroscopy, molecular docking simulations, molecular dynamics, protein conformational clusters (EnGens), molecular fractionation with conjugate caps (MFCC) and quantum biochemistry analysis. We analyze molecular and quantum aspects, uncovering interaction energies between sugar atoms and amino acids. Total interaction energy considers protein fragmentation, energetic decomposition, and interaction energy from a bottom-up perspective. Molecular dynamics reveal that unmodified Aldehydo-D-glucose (GLO) escapes HSA binding sites, explaining gradual glycation. We pioneer studying HSA's binding mechanism with glucose and fructose in a 1:1 ratio using long molecular dynamics simulations. Results suggest the transitional GLOm form has a higher Sudlow I site propensity than unmodified glucose, crucial for K195 glycation. FRO and GLOm interaction tendencies move toward a deeper FA7 cavity, near its center. This approach effectively elucidates small molecule binding mechanisms, consistent with previous experimental results.Communicated by Ramaswamy H. Sarma.

Noncovalent binding of carbofuran to acetylcholinesterase from Homo sapiens, Danio rerio, Apis mellifera and Caenorhabditis elegans: Homology modelling, molecular docking and dynamics, and quantum biochemistry description.

Although various regulatory agencies have banned or severely restricted the use of carbofuran (CAR), recent reports indicate the presence of CAR residues in both cultivated and wild areas. This pesticide is a potent inhibitor of acetylcholinesterase (AChE), which acts by preventing the hydrolysis of acetylcholine (ACh). Given the critical role of AChE::ACh in the proper functioning of the nervous system, we thought it appropriate to investigate the binding of CAR to AChEs from Homo sapiens, Danio rerio, Apis mellifera, and Caenorhabditis elegans using homology modelling, molecular docking, molecular dynamics, and quantum biochemistry. Molecular docking and dynamics results indicated peculiar structural behavior in each AChE::CAR system. Quantum biochemistry results showed similar affinities for all complexes, confirming the description of carbofuran as a broad-spectrum pesticide, and have a limited correlation with IC50 values. We found the following decreasing affinity order of AChE species: H. sapiens > A. mellifera > C. elegans > D. rerio. The computational results suggest that CAR occupies different pockets in the AChEs studied. In addition, our results showed that CAR binds to hsAChE and ceAChE in a very similar manner: it has high affinities for the same subsites in both species and forms hydrogen bonds with residues (hsTYR124 and ceTRP107) occupying homologous positions in the peripheral site. This suggests that this nematode is a potential model to evaluate the toxicity of carbamates, even though the sequence identity between them is only 41 %. Interestingly, we also observed that the catalytic histidines of drAChE and amAChE exhibited favorable contacts with carbofuran, suggesting that the non-covalent binding of carbofuran to these proteins may promote faster carbamylation rates than the binding modes to human and worm acetylcholinesterases. Our computational results provide a better understanding of the binding mechanisms in these complexes, as well as new insights into the mechanism of carbamylation.

In silico description of the adsorption of cell signaling pathway proteins ovalbumin, glutathione, LC3, TLR4, ASC PYCARD, PI3K and NF-Kß on 7.0 nm gold nanoparticles: obtaining their Lennard-Jones-like potentials through docking and molecular mechanics.

The impact of vaccination on the world's population is difficult to calculate. For developing different types of vaccines, adjuvants are substances added to vaccines to increase the magnitude and durability of the immune response and the effectiveness of the vaccine. This work explores the potential use of spherical gold nanoparticles (AuNPs) as adjuvants. Thus, we employed docking techniques and molecular mechanics to describe how a AuNP 7.0 nm in diameter interacts with cell signaling pathway proteins. Initially, we used X-ray crystallization data of the proteins ovalbumin, glutathione, LC3, TLR4, ASC PYCARD, PI3K, and NF-Kß to study the adsorption with an AuNP through molecular docking. Therefore, interaction energies were obtained for the AuNP complexes and individual proteins, as well as the AuNP and OVA complex (AuNP@OVA) with each cellular protein, respectively. Results showed that AuNPs had the highest affinity for OVA individually, followed by glutathione, ASC PYCARD domain, LC3, PI3K, NF-Kß, and TLR4. Furthermore, when evaluating the AuNP@OVA complex, glutathione showed a greater affinity with more potent interaction energy when compared to the other studied systems.

Pharmacokinetic Profile Evaluation of Novel Combretastatin Derivative, LASSBio-1920, as a Promising Colorectal Anticancer Agent.

LASSBio-1920 was synthesized due to the poor solubility of its natural precursor, combretastatin A4 (CA4). The cytotoxic potential of the compound against human colorectal cancer cells (HCT-116) and non-small cell lung cancer cells (PC-9) was evaluated, yielding IC50 values of 0.06 and 0.07 µM, respectively. Its mechanism of action was analyzed by microscopy and flow cytometry, where LASSBio-1920 was found to induce apoptosis. Molecular docking simulations and the enzymatic inhibition study with wild-type (wt) EGFR indicated enzyme-substrate interactions similar to other tyrosine kinase inhibitors. We suggest that LASSBio-1920 is metabolized by O-demethylation and NADPH generation. LASSBio-1920 demonstrated excellent absorption in the gastrointestinal tract and high central nervous system (CNS) permeability. The pharmacokinetic parameters obtained by predictions indicated that the compound presents zero-order kinetics and, in a human module simulation, accumulates in the liver, heart, gut, and spleen. The pharmacokinetic parameters obtained will serve as the basis to initiate in vivo studies regarding LASSBio-1920's antitumor potential.

Losartan as an ACE inhibitor: a description of the mechanism of action through quantum biochemistry.

Losartan (LST) is a potent and selective angiotensin II (Ang II) type 1 (AT1) receptor antagonist widely used in the treatment of hypertension. The formation of Ang II is catalyzed by the angiotensin I-converting enzyme (ACE) through proteolytic cleavage of angiotensin I (Ang I), which is involved in the control of blood pressure. Despite the vast literature on the relationship of losartan with the renin-angiotensin system (RAS), the actions of losartan on the sACE enzyme are so far poorly understood. In view of this, we investigated how losartan can interact with the sACE enzyme to block its activity and intracellular signaling. After performing docking assays following quantum biochemistry calculations using losartan and sACE crystallographic data, we report that their interaction results reveal a new mechanism of action with important implications for understanding its effects on hypertension.

Characterization of the binding interaction between atrazine and human serum albumin: Fluorescence spectroscopy, molecular dynamics and quantum biochemistry.

Atrazine (ATR), one of the most used herbicides worldwide, causes persistent contamination of water and soil due to its high resistance to degradation. ATR is associated with low fertility and increased risk of prostate cancer in humans, as well as birth defects, low birth weight and premature delivery. Describing ATR binding to human serum albumin (HSA) is clinically relevant to future studies about pharmacokinetics, pharmacodynamics and toxicity of ATR, as albumin is the most abundant carrier protein in plasma and binds important small biological molecules. In this work we characterize, for the first time, the binding of ATR to HSA by using fluorescence spectroscopy and performing simulations using molecular docking, classical molecular dynamics and quantum biochemistry based on density functional theory (DFT). We determine the most likely binding sites of ATR to HSA, highlighting the fatty acid binding site FA8 (located between subdomains IA-IB-IIA and IIB-IIIA-IIIB) as the most important one, and evaluate each nearby amino acid residue contribution to the binding interactions explaining the fluorescence quenching due to ATR complexation with HSA. The stabilization of the ATR/FA8 complex was also aided by the interaction between the atrazine ring and SER454 (hydrogen bond) and LEU481(alkyl interaction).

Gallic acid leads to cell death of Candida albicans by the apoptosis mechanism.

Aim: To evaluate the antifungal activity of gallic acid (GA) against the strains of Candida spp. resistant to fluconazole and to determine its mechanism of action. Materials & methods: Antifungal activity was evaluated using the broth microdilution and flow cytometry techniques. Results: GA presented minimum inhibitory concentrations ranging from 16 to 72 µg/ml, causing alterations of the membrane integrity and mitochondrial transmembrane potential, production of reactive oxygen species and externalization of phosphatidylserine. Conclusion: GA has potential antifungal activity against Candida spp.

ACE2-derived peptides interact with the RBD domain of SARS-CoV-2 spike glycoprotein, disrupting the interaction with the human ACE2 receptor.

Vaccines could be the solution to the current SARS-CoV-2 outbreak. However, some studies have shown that the immunological memory only lasts three months. Thus, it is imperative to develop pharmacological treatments to cope with COVID-19. Here, the in silico approach by molecular docking, dynamic simulations and quantum biochemistry revealed that ACE2-derived peptides strongly interact with the SARS-CoV-2 RBD domain of spike glycoprotein (S-RBD). ACE2-Dev-PepI, ACE2-Dev-PepII, ACE2-Dev-PepIII and ACE2-Dev-PepIV complexed with S-RBD provoked alterations in the 3D structure of S-RBD, leading to disruption of the correct interaction with the ACE2 receptor, a pivotal step for SARS-CoV-2 infection. This wrong interaction between S-RBD and ACE2 could inhibit the entry of SARS-CoV-2 in cells, and thus virus replication and the establishment of COVID-19 disease. Therefore, we suggest that ACE2-derived peptides can interfere with recognition of ACE2 in human cells by SARS-CoV-2 in vivo. Bioinformatic prediction showed that these peptides have no toxicity or allergenic potential. By using ACE2-derived peptides against SARS-CoV-2, this study points to opportunities for further in vivo research on these peptides, seeking to discover new drugs and entirely new perspectives to treat COVID-19.Communicated by Ramaswamy H. Sarma.

Quantum biochemistry, molecular docking, and dynamics simulation revealed synthetic peptides induced conformational changes affecting the topology of the catalytic site of SARS-CoV-2 main protease.

The recent outbreak caused by SARS-CoV-2 continues to threat and take many lives all over the world. The lack of an efficient pharmacological treatments are serious problems to be faced by scientists and medical staffs worldwide. In this work, an in silico approach based on the combination of molecular docking, dynamics simulations, and quantum biochemistry revealed that the synthetic peptides RcAlb-PepI, PepGAT, and PepKAA, strongly interact with the main protease (Mpro) a pivotal protein for SARS-CoV-2 replication. Although not binding to the proteolytic site of SARS-CoV-2 Mpro, RcAlb-PepI, PepGAT, and PepKAA interact with other protein domain and allosterically altered the protease topology. Indeed, such peptide-SARS-CoV-2 Mpro complexes provoked dramatic alterations in the three-dimensional structure of Mpro leading to area and volume shrinkage of the proteolytic site, which could affect the protease activity and thus the virus replication. Based on these findings, it is suggested that RcAlb-PepI, PepGAT, and PepKAA could interfere with SARS-CoV-2 Mpro role in vivo. Also, unlike other antiviral drugs, these peptides have no toxicity to human cells. This pioneering in silico investigation opens up opportunity for further in vivo research on these peptides, towards discovering new drugs and entirely new perspectives to treat COVID-19.Communicated by Ramaswamy H. Sarma.

Computational approach, scanning electron and fluorescence microscopies revealed insights into the action mechanisms of anticandidal peptide Mo-CBP3-PepIII.

AIMS: The Candida genus is composed of opportunistic pathogens that threaten public health. Given the increase in resistance to current drugs, it is necessary to develop new drugs to treat infections by these pathogens. Antimicrobial peptides are promising alternative molecules with low cost, broad action spectrum and low resistance induction. This study aimed to clarify the action mechanisms of synthetic peptides against Candida albicans. MAIN METHODS: The mode of action of the anticandidal peptides Mo-CBP3-PepIII were analyzed through molecular dynamics and quantum biochemistry methods against Exo-ß-1,3-glucanase (EXG), vital to cell wall metabolism. Furthermore, scanning electron (SEM) and fluorescence (FM) microscopies were employed to corroborate the in silico data. KEY FINDINGS: Mo-CBP3-PepIII strongly interacted with EXG (-122.2 kcal mol-1) at the active site, higher than the commercial inhibitor pepstatin. Also, molecular dynamics revealed the insertion of Mo-CBP3-PepIII into the yeast membrane. SEM analyses revealed that Mo-CBP3-PepIII induced cracks and scars of the cell wall and FM analyses confirmed the pore formation on the Candida membrane. SIGNIFICANCE: Mo-CBP3-PepIII has strong potential as a new drug with a broad spectrum of action, given its different mode of action compared to conventional drugs.

Molecular insight on the binding of stevia glycosides to bovine serum albumin.

The interaction of the steviol and its glycosides (SG), steviolbioside, and rebaudioside A, with bovine serum albumin (BSA) was studied by absorption and fluorescence spectroscopy techniques alongside molecular docking. The stevia derivatives quenched the fluorescence of BSA by a dynamic quenching mechanism, indicating the interaction between the stevia derivatives and BSA. The binding constant (Kb) of steviol was 100-1000-fold higher than those of SG. The stevia derivative/BSA binding reaction was spontaneous and involved the formation of hydrogen bonds and van der Waals interactions between steviol and steviolbioside with BSA, and water reorganization around the rebaudioside A/BSA complex. Molecular docking pointed out the FA1 and FA9 binding sites of BSA as the probable binding sites of steviol and SG, respectively. In conclusion, steviol enhanced hydrophobicity and small size compared to SG may favor its binding to BSA. As steviol and its glycosides share binding sites on BSA with free fatty acids and drugs, they may be competitively displaced from plasma albumin under various physiological states or disease conditions. These findings are clinically relevant and provide an insight into the pharmacokinetics and pharmacodynamics of the stevia glycosides.

Quantum biochemistry in cancer immunotherapy: New insights about CTLA-4/ipilimumab and design of ipilimumab-derived peptides with high potential in cancer treatment.

Cancer is a group of diseases involving disordered growth of abnormal cells with the potential to invade and spread to other parts of the body. Today, immunotherapy is the most efficient treatment, with fewer side effects. Notably, the employment of monoclonal antibodies to inhibit checkpoint proteins, such as CTLA-4, has caused much excitement among cancer immunotherapy researchers. Thus, in-depth analysis through quantum biochemistry and molecular dynamics simulations was performed to understand the complex formed by ipilimumab and its target CTLA-4. Our computational results provide a better understanding of the binding mechanisms and new insights about the CTLA-4: ipilimumab interaction, identifying essential amino acid residues to support the complex. Additionally, we report new interactions such as aromatic-aromatic, aromatic-sulfur, and cation-pi interactions to stabilize the CTLA-4:ipilimumab complex. Finally, quantum biochemistry analyses reveal the most important amino acid residues involved in the CTLA-4:ipilimumab interface, which were used to design synthetic peptides to inhibit CTLA-4. The computational results presented here provide a better understanding of the CTLA-4:ipilimumab binding mechanisms, and can support the development of alternative antibody-based drugs with high relevance in cancer immunotherapy.

Crystal structure and specific location of a germin-like protein with proteolytic activity from Thevetia peruviana.

Peruvianin-I is a cysteine peptidase (EC 3.4.22) purified from Thevetia peruviana. Previous studies have shown that it is the only germin-like protein (GLP) with proteolytic activity described so far. In this work, the X-ray crystal structure of peruvianin-I was determined to a resolution of 2.15 Å (PDB accession number: 6ORM) and its specific location was evaluated by different assays. Its overall structure shows an arrangement composed of a homohexamer (a trimer of dimers) where each monomer exhibits a typical ß-barrel fold and two glycosylation sites (Asn55 and Asn144). Analysis of its active site confirmed the absence of essential amino acids for typical oxalate oxidase activity of GLPs. Details of the active site and molecular docking results, using a specific cysteine peptidase inhibitor (iodoacetamide), were used to discuss a plausible mechanism for proteolytic activity of peruvianin-I. Histological analyses showed that T. peruviana has articulated anastomosing laticifers, i.e., rows of cells which merge to form continuous tubes throughout its green organs. Moreover, peruvianin-I was detected exclusively in the latex. Because latex peptidases have been described as defensive molecules against insects, we hypothesize that peruvianin-I contributes to protect T. peruviana plants against herbivory.

Antitumor Potential of the Isoflavonoids (+)- and (-)-2,3,9-Trimethoxypterocarpan: Mechanism-of-Action Studies.

Synthetically derived samples of (+)-(6aS,11aS)-2,3,9-trimethoxypterocarpan [(+)-1] and its enantiomer [(-)-1], both of which are examples of naturally occurring isoflavonoids, were evaluated, together with the corresponding racemate, as cytotoxic agents against the HL-60, HCT-116, OVCAR-8, and SF-295 tumor cell lines. As a result it was established that compound (+)-1 was particularly active with OVCAR-8 cells being the most sensitive and responding in a dose-dependent manner. A study of cell viability and drug-induced morphological changes revealed the compound causes cell death through a mechanism characteristic of apoptosis. Finally, a computational study of the interactions of compound (+)-1 and (S)-monastrol, an established, synthetically derived, potent, and cell-permeant inhibitor of mitosis, with the kinesin-type protein Eg5 revealed that both bind to this receptor in a similar manner. Significantly, compound (+)-1 binds with greater affinity, an effect attributed to the presence of the associated methoxy groups.

Ribosomal RNA-Aminoglycoside Hygromycin B Interaction Energy Calculation within a Density Functional Theory Framework.

We intend to investigate the drug-binding energy of each nucleotide inside the aminoglycoside hygromycin B (hygB) binding site of 30S ribosomal RNA (rRNA) subunit by using the molecular fractionation with conjugate caps (MFCC) strategy based on the density functional theory (DFT), considering the functional LDA/PWC, OBS, and the dielectric constant parametrization. Aminoglycosides are bactericidal antibiotics that have high affinity to the prokaryotic rRNA, inhibiting the synthesis of proteins by acting on the main stages of the translation mechanism, whereas binding to rRNA 16S, a component of the 30S ribosomal subunit in prokaryotes. The identification of the nucleotides presenting the most negative binding energies allows us to stabilize hygB in a suitable binding pocket of the 30S ribosomal subunit. In addition, it should be highlighted that mutations in these residues may probably lead to resistance to ribosome-targeting antibiotics. Quantum calculations of aminoglycoside hygromycin B-ribosome complex might contribute to further quantum studies with antibiotics like macrolides and other aminoglycosides.

Explaining urokinase type plasminogen activator inhibition by amino-5-hydroxybenzimidazole and two naphthamidine-based compounds through quantum biochemistry.

Urokinase plasminogen activator (uPA) is a biomarker and therapeutic target for several cancer types whose inhibition has been shown to slow tumor growth and metastasis. In this work, crystallographic data of uPA complexed with distinct ligands (PDB id: 1SQA, 1SQO, and 1FV9) were used to perform quantum biochemistry calculations based on the framework of density functional theory (DFT) and within the molecular fractionation with conjugated caps (MFCC) scheme. Our calculations revealed a total energy interaction of -107.30, -99.5, and -35.30 kcal mol-1 for two naphthamidine-based compounds (Ul1 and UI2) and 2-amino-5-hydroxybenzimidazole (172), respectively, which are in good agreement with known inhibitory experiments. Residues Asp189, Ser190, Cys191-Cys220, Gln192, Trp 215, Gly216, and Gly219 were identified as the main interacting amino acid residues with interaction energy contributions lower than -4.0 kcal mol-1 for uPA/UI1 and UPA/UI2 complexes. In the case of compound 172, our calculations have shown that the most important interactions occur with residues Asp189, Cys191-Cys220, and Ser190. Our results highlight the relevance of the protonation state of ligands and residues and that the naphthamidine scaffold of UI1 and UI2 is the main determinant of their potency, followed by their aminopyrimidine substitution. Altogether, the results of this work contribute to the understanding of the uPA binding mechanisms of the inhibitory compounds Ul1 and 172, stimulating the use of quantum biochemistry theoretical approaches for the development of new uPA inhibitors as new medicines for cancer treatment.

Cloning of cDNA sequences encoding cowpea (Vigna unguiculata) vicilins: Computational simulations suggest a binding mode of cowpea vicilins to chitin oligomers.

Vicilins are 7S globulins which constitute the major seed storage proteins in leguminous species. Variant vicilins showing differential binding affinities for chitin have been implicated in the resistance and susceptibility of cowpea to the bruchid Callosobruchus maculatus. These proteins are members of the cupin superfamily, which includes a wide variety of enzymes and non-catalytic seed storage proteins. The cupin fold does not share similarity with any known chitin-biding domain. Therefore, it is poorly understood how these storage proteins bind to chitin. In this work, partial cDNA sequences encoding ß-vignin, the major component of cowpea vicilins, were obtained from developing seeds. Three-dimensional molecular models of ß-vignin showed the characteristic cupin fold and computational simulations revealed that each vicilin trimer contained 3 chitin-binding sites. Interaction models showed that chito-oligosaccharides bound to ß-vignin were stabilized mainly by hydrogen bonds, a common structural feature of typical carbohydrate-binding proteins. Furthermore, many of the residues involved in the chitin-binding sites of ß-vignin are conserved in other 7S globulins. These results support previous experimental evidences on the ability of vicilin-like proteins from cowpea and other leguminous species to bind in vitro to chitin as well as in vivo to chitinous structures of larval C. maculatus midgut.

The vibrational properties of the bee-killer imidacloprid insecticide: A molecular description.

The chemical imidacloprid belongs to the neonicotinoids insecticide class, widely used for insect pest control mainly for crop protection. However, imidacloprid is a non-selective agrochemical to the insects and it is able to kill the most important pollinators, the bees. The high toxicity of imidacloprid requires controlled release and continuous monitoring. For this purpose, high performance liquid chromatography (HPLC) is usually employed; infrared and Raman spectroscopy, however, are simple and viable techniques that can be adapted to portable devices for field application. In this communication, state-of-the-art quantum level simulations were used to predict the infrared and Raman spectra of the most stable conformer of imidacloprid. Four molecular geometries were investigated in vacuum and solvated within the Density Functional Theory (DFT) approach employing the hybrid meta functional M06-2X and the hybrid functional B3LYP. The M062X/PCM model proved to be the best to predict structural features, while the values of harmonic vibrational frequencies were predicted more accurately using the B3LYP functional.

cDNA cloning, molecular modeling and docking calculations of L-type lectins from Swartzia simplex var. grandiflora (Leguminosae, Papilionoideae), a member of the tribe Swartzieae.

The genus Swartzia is a member of the tribe Swartzieae, whose genera constitute the living descendants of one of the early branches of the papilionoid legumes. Legume lectins comprise one of the main families of structurally and evolutionarily related carbohydrate-binding proteins of plant origin. However, these proteins have been poorly investigated in Swartzia and to date, only the lectin from S. laevicarpa seeds (SLL) has been purified. Moreover, no sequence information is known from lectins of any member of the tribe Swartzieae. In the present study, partial cDNA sequences encoding L-type lectins were obtained from developing seeds of S. simplex var. grandiflora. The amino acid sequences of the S. simplex grandiflora lectins (SSGLs) were only averagely related to the known primary structures of legume lectins, with sequence identities not greater than 50-52%. The SSGL sequences were more related to amino acid sequences of papilionoid lectins from members of the tribes Sophoreae and Dalbergieae and from the Cladratis and Vataireoid clades, which constitute with other taxa, the first branching lineages of the subfamily Papilionoideae. The three-dimensional structures of 2 representative SSGLs (SSGL-A and SSGL-E) were predicted by homology modeling using templates that exhibit the characteristic ß-sandwich fold of the L-type lectins. Molecular docking calculations predicted that SSGL-A is able to interact with D-galactose, N-acetyl-D-galactosamine and α-lactose, whereas SSGL-E is probably a non-functional lectin due to 2 mutations in the carbohydrate-binding site. Using molecular dynamics simulations followed by density functional theory calculations, the binding free energies of the interaction of SSGL-A with GalNAc and α-lactose were estimated as -31.7 and -47.5 kcal/mol, respectively. These findings gave insights about the carbohydrate-binding specificity of SLL, which binds to immobilized lactose but is not retained in a matrix containing D-GalNAc as ligand.