Interplay between Conformational Entropy and Solvation Entropy in Protein-Ligand Binding

Effect of phenolic-hydroxy-group incorporation on the biological activity of a simplified aplysiatoxin analog with an (R)-(-)-carvone-based core

We synthesized a phenolic hydroxy group-bearing version (1) of a simplified analog of aplysiatoxin comprising a carvone-based conformation-controlling unit. Thereafter, we evaluated its antiproliferative activity against human cancer cell lines and its binding affinity to protein kinase C (PKC) isozymes. The antiproliferative activity and PKC-binding ability increased with the introduction of the phenolic hydroxy group. The results of molecular dynamics simulations and subsequent relative binding free-energy calculations conducted using an alchemical transformation procedure showed that the phenolic hydroxy group in 1 could form a hydrogen bond with a phospholipid and the PKC. The former hydrogen bonding formation facilitated the partitioning of the compound from water to the phospholipid membrane and the latter compensated for the loss of hydrogen bond with the phospholipid upon binding to the PKC. This information may facilitate the development of rational design methods for PKC ligands with additional hydrogen bonding groups.

Influence of Selective Deoxyfluorination on the Molecular Structure of Type-2 N-Acetyllactosamine

N-Acetyllactosamine is a common saccharide motif found in various biologically active glycans. This motif usually works as a backbone for additional modifications and thus significantly influences glycan conformational behavior and biological activity. In this work, we have investigated the type-2 N-acetyllactosamine scaffold using the complete series of its monodeoxyfluorinated analogs. These glycomimetics have been studied by molecular mechanics, quantum mechanics, X-ray crystallography, and various NMR techniques, which have provided a comprehensive and complete insight into the role of individual hydroxyl groups in the conformational behavior and lipophilicity of N-acetyllactosamine.

Thermodynamic consequences of stapling side-chains on a peptide ligand using a lactam-bridge: A theoretical study on anti-angiogenic peptides targeting VEGF

Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide-based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side-chain lactam-bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side-chain lactam-bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all-atom molecular dynamics simulations. It is found that introducing a side-chain lactam-bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand-target binding upon imposing a conformational constraint into a flexible ligand.

19F Fast Magic-Angle Spinning NMR Spectroscopy on Microcrystalline Complexes of Fluorinated Ligands and the Carbohydrate Recognition Domain of Galectin-3

Structural characterization of protein-ligand binding interfaces at atomic resolution is essential for improving the design of specific and potent inhibitors. Herein, we explored fast 19F- and 1H-detected magic angle spinning NMR spectroscopy to investigate the interaction between two fluorinated ligand diastereomers with the microcrystalline galectin-3 carbohydrate recognition domain. The detailed environment around the fluorine atoms was mapped by 2D 13C-19F and 1H-19F dipolar correlation experiments and permitted characterization of the binding interface. Our results demonstrate that 19F MAS NMR is a powerful tool for detailed characterization of protein-ligand interfaces and protein interactions at the atomic level.

Selection of Galectin-Binding Ligands from Synthetic Glycopeptide Libraries

Galectins, a class of carbohydrate-binding proteins, play a crucial role in various physiological and disease processes. Therefore, the identification of ligands that efficiently bind these proteins could potentially lead to the development of new therapeutic compounds. In this study, we present a method that involves screening synthetic click glycopeptide libraries to identify lectin-binding ligands with low micromolar affinity. Our methodology, initially optimized using Concanavalin A, was subsequently applied to identify binders for the therapeutically relevant galectin 1. Binding affinities were assessed using various methods and showed that the selected glycopeptides exhibited enhanced binding potency to the target lectins compared to the starting sugar moieties. This approach offers an alternative means of discovering galectin-binding ligands as well as other carbohydrate-binding proteins, which are considered important therapeutic targets.

Emerging Conformational-Analysis Protocols from the RTCONF55-16K Reaction Thermochemistry Conformational Benchmark Set

RTCONF55-16K is a new, reactive conformational data set based on cost-efficient methods to assess different conformational analysis protocols. Our reference calculations underpinned the accuracy of the CENSO (Grimme et al. J. Phys. Chem. A, 2021, 125, 4039) procedure and resulted in alternative recipes with different cost-accuracy compromises. Our general-purpose and economical protocols (CENSO-light and zero, respectively) were found to be 10-30 times faster than the original algorithm, adding only 0.4-0.7 kcal/mol absolute error to the relative free energy estimates.

Insect cuticle-inspired design of sustainably sourced composite bioplastics with enhanced strength, toughness and stretch-strengthening behavior

Insect cuticles that are mainly made of chitin, chitosan and proteins provide insects with rigid, stretchable and robust skins to defend harsh external environment. The insect cuticle therefore provides inspiration for engineering biomaterials with outstanding mechanical properties but also sustainability and biocompatibility. We herein propose a design of high-performance and sustainable bioplastics via introducing CPAP3-A1, a major structural protein in insect cuticles, to specifically bind to chitosan. Simply mixing 10w/w% bioengineered CPAP3-A1 protein with chitosan enables the formation of plastics-like, sustainably sourced chitosan/CPAP3-A1 composites with significantly enhanced strength (∼90 MPa) and toughness (∼20 MJ m -3), outperforming previous chitosan-based composites and most synthetic petroleum-based plastics. Remarkably, these bioplastics exhibit a stretch-strengthening behavior similar to the training living muscles. Mechanistic investigation reveals that the introduction of CPAP3-A1 induce chitosan chains to assemble into a more coarsened fibrous network with increased crystallinity and reinforcement effect, but also enable energy dissipation via reversible chitosan-protein interactions. Further uniaxial stretch facilitates network re-orientation and increases chitosan crystallinity and mechanical anisotropy, thereby resulting in stretch-strengthening behavior. In general, this study provides an insect-cuticle inspired design of high-performance bioplastics that may serve as sustainable and bio-friendly materials for a wide range of engineering and biomedical application potentials.

A high-resolution structural characterization and physicochemical study of how a peptoid binds to an oncoprotein MDM2

Peptoids are a promising drug modality targeting disease-related proteins, but how a peptoid engages in protein binding is poorly understood. This is primarily due to a lack of high-resolution peptoid-protein complex structures and systematic physicochemical studies. Here, we present the first crystal structure of a peptoid bound to a protein, providing high-resolution structural information about how a peptoid binds to a protein. We previously reported a rigid peptoid, oligo(N-substituted alanine) (oligo-NSA), and developed an oligo-NSA-type peptoid that binds to MDM2. X-ray crystallographic analysis of the peptoid bound to MDM2 showed that the peptoid recognizes the MDM2 surface predominantly through the interaction of the N-substituents, while the main chain acts as a scaffold. Additionally, conformational, thermodynamic, and kinetic analysis of the peptoid and its derivatives with a less rigid main chain revealed that rigidification of the peptoid main chain contributes to improving the protein binding affinity. This improvement is thermodynamically attributed to an increased magnitude of the binding enthalpy change, and kinetically to an increased association rate and decreased dissociation rate. This study provides invaluable insights into the design of protein-targeting peptoids.

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.

Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

Discovery of Novel Inhibitors of BRD4 for Treating Prostate Cancer: A Comprehensive Case Study for Considering Water Networks in Virtual Screening and Drug Design

Androgen receptor (AR) is the primary target for treating prostate cancer (PCa), which inevitably progresses due to drug-resistant mutations. Bromodomain-containing protein 4 (BRD4) has been a new potential drug target for PCa treatment. Herein, we report the rational design and discovery of novel BRD4 inhibitors through computer-aided drug design (CADD), and a hit compound SQ-1 (IC50 = 676 nM) was identified by structure-based virtual screening (SBVS) with the conserved water network. To optimize the structure of SQ-1, the free energy landscape was constructed, and the binding mechanism was explored by characterizing the water profile and the dissociation mechanism. Finally, the compound SQ-17 with improved inhibitory activity (IC50 < 100 nM) was discovered, which showed potent antiproliferative activity against LNCaP. These data highlighted a successful attempt to identify and optimize a small molecule by comprehensive CADD application and provided essential clues for developing novel therapeutics for PCa treatment.

Small-Molecule Conformer Generators: Evaluation of Traditional Methods and AI Models on High-Quality Data Sets

Small-molecule conformer generation (SMCG) is an extremely important task in both ligand- and structure-based computer-aided drug design, especially during the hit discovery phase. Recently, a multitude of artificial intelligence (AI) models tailored for SMCG have emerged. Despite developers typically furnishing performance evaluation data upon releasing their AI models, a comprehensive and equitable performance comparison between AI models and conventional methods is still lacking. In this study, we curated a new benchmarking data set comprising 3354 high-quality ligand bioactive conformations. Subsequently, we conducted a systematic assessment of the performance of four widely adopted traditional methods (i.e., ConfGenX, Conformator, OMEGA, and RDKit ETKDG) and five AI models (i.e., ConfGF, DMCG, GeoDiff, GeoMol, and torsional diffusion) in the tasks of reproducing bioactive and low-energy conformations of small molecules. In the former task, the AI models have no advantage, particularly with a maximum ensemble size of 1. Even the best-performing AI model GeoMol is still worse than any of the tested traditional methods. Conversely, in the latter task, the torsional diffusion model shows obvious advantages, surpassing the best-performing traditional method ConfGenX by 26.09 and 12.97% on the COV-R and COV-P metrics, respectively. Furthermore, the influence of force field-based fine-tuning on the quality of the generated conformers was also discussed. Finally, a user-friendly Web server called fastSMCG was developed to enable researchers to rapidly and flexibly generate small-molecule conformers using both traditional and AI methods. We anticipate that our work will offer valuable practical assistance to the scientific community in this field.

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.

1,2,3-Triazoles in Biomolecular Crystallography: A Geometrical Data-Mining Approach

The 1,2,3-triazole scaffold has become very attractive to identify new chemical entities in drug discovery projects. Despite the widespread use of click chemistry to synthesize numerous 123Ts, there are few drugs on the market that incorporate this scaffold as a substructure. To investigate the true potential of 123Ts in protein-ligand interactions, we examined the noncovalent interactions between the 1,2,3-triazole ring and amino acids in protein-ligand cocrystals using a geometrical approach. For this purpose, we constructed a nonredundant database of 220 PDB IDs from available 123T-protein cocrystal structures. Subsequently, using the Protein Ligand Interaction Profiler web platform (PLIP), we determined whether 1,2,3-triazoles primarily act as linkers or if they can be considered interactive scaffolds. We then manually analyzed the geometrical descriptors from 333 interactions between 1,4-disubstituted 123T rings and amino acid residues in proteins. This study demonstrates that 1,2,3-triazoles exhibit diverse preferred interactions with amino acids, which contribute to protein-ligand binding.

Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation

The sulfation patterns and degree of sulfation of chondroitin sulfate (CS), an important class of glycosaminoglycans (GAG), and their interactions with chemokines are accountable for various diseases. To realize the underlying mechanism of such complex biological phenomena at a molecular level and their application in rational drug design, a study on conformations and dynamics of CSs is necessary. To explore this, in this study, we performed a series of atomistic molecular dynamics (MD) simulations with different sulfated variants of octadecasaccharide CS, like CS-C, CS-E, and CS-T, in their free forms and when bound to the protein chemokine CXCL8 dimer in an aqueous medium. The calculated binding free energy of CSs with the CXCL8 dimer is favorable, and the degree of sulfation favors the complexation process further with prominent hydrophobic and hydrogen-bonded interactions. We find that the recognition is associated with the configurational entropy loss of the CS molecules as calculated from the Gaussian mixture approach, which supports that the degree of sulfation regulates the process. Cluster analysis through the k-means algorithm and end-to-end distance measurement revealed that although the free CS molecules adopted linear conformations, the nonlinear conformations during binding with protein were noted. Adaptation of nonlinear forms in the bound forms is noteworthy for the less-sulfated CS-C and CS-E. Apart from favorable 4C1 conformations, the occasional appearance of skew-boat forms from the free-energy map of ring pucker for the GlcUA unit was observed, which remains unaffected by the sulfation. We find that during recognition, the average relaxation time of intra-CS and inter-CS-CXCL8 hydrogen bonds (HBs) is about a magnitude lesser than that of CS-water HBs, most prominent on the involvement of higher sulfated CS-T analogues. The translational motion of surrounded water molecules in CSs exhibited sublinear diffusion, and the degree of sublinearity increases around the heavily sulfated molecules due to the hindrance created by them as well as the presence of the chemokine and exhibited markedly slow heterogeneous diffusion.

Conformational study into N-alkyl-N'-aryl ureas to inform drug discovery

Ureas are an important functional group in small molecule drugs as well as having wider applications in organic chemistry. Understanding of their conformation is of critical importance for rational design of urea-containing bioactive compounds. Whilst the conformational preferences of biaryl ureas have been extensively studied, very little attention has been paid to alkylated analogues. We carried out a systematic study of N-aryl (phenyl and pyridyl)-N'-cyclopentyl ureas with differing N-methylation patterns using Well Tempered Metadynamics at a semi-empirical level in implicit water (GBSA) using Well-Tempered Metadynamics to generate their conformational free-energy landscapes. Geometries and energetics of the most relevant configurations are further refined using DFT level of theory. Validation for the computation was obtained by synthesis of all 8 analogues followed by conformational studies by X-ray crystallography and NMR. These findings reveal that the methylation pattern significantly affects the conformational preference of the system. Most notably, N-phenyl-N'-cyclopentyl urea is shown to adopt both the trans-trans, and cis-trans conformations with equal energy and that the cis-trans conformation can be significantly stabilised by the presence of an internal hydrogen bond to the N'-hydrogen. This study will be of utility for the design of N-alkyl-N'-aryl ureas as drug candidates.

Improving binding entropy by higher ligand symmetry? - A case study with human matriptase

Understanding different contributions to the binding entropy of ligands is of utmost interest to better predict affinity and the thermodynamic binding profiles of protein-ligand interactions and to develop new strategies for ligand optimization. To these means, the largely neglected effects of introducing higher ligand symmetry, thereby reducing the number of energetically distinguishable binding modes on binding entropy using the human matriptase as a model system, were investigated. A set of new trivalent phloroglucinol-based inhibitors that address the roughly symmetric binding site of the enzyme was designed, synthesized, and subjected to isothermal titration calorimetry. These highly symmetric ligands that can adopt multiple indistinguishable binding modes exhibited high entropy-driven affinity in line with affinity-change predictions.

Selectively Modified Lactose and N-Acetyllactosamine Analogs at Three Key Positions to Afford Effective Galectin-3 Ligands

Galectins constitute a family of galactose-binding lectins overly expressed in the tumor microenvironment as well as in innate and adaptive immune cells, in inflammatory diseases. Lactose ((β-D-galactopyranosyl)-(1→4)-β-D-glucopyranose, Lac) and N-Acetyllactosamine (2-acetamido-2-deoxy-4-O-β-D-galactopyranosyl-D-glucopyranose, LacNAc) have been widely exploited as ligands for a wide range of galectins, sometimes with modest selectivity. Even though several chemical modifications at single positions of the sugar rings have been applied to these ligands, very few examples combined the simultaneous modifications at key positions known to increase both affinity and selectivity. We report herein combined modifications at the anomeric position, C-2, and O-3' of each of the two sugars, resulting in a 3'-O-sulfated LacNAc analog having a Kd of 14.7 µM against human Gal-3 as measured by isothermal titration calorimetry (ITC). This represents a six-fold increase in affinity when compared to methyl β-D-lactoside having a Kd of 91 µM. The three best compounds contained sulfate groups at the O-3' position of the galactoside moieties, which were perfectly in line with the observed highly cationic character of the human Gal-3 binding site shown by the co-crystal of one of the best candidates of the LacNAc series.

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

Total free energy analysis of fully hydrated proteins

The total free energy of a hydrated biomolecule and its corresponding decomposition of energy and entropy provides detailed information about regions of thermodynamic stability or instability. The free energies of four hydrated globular proteins with different net charges are calculated from a molecular dynamics simulation, with the energy coming from the system Hamiltonian and entropy using multiscale cell correlation. Water is found to be most stable around anionic residues, intermediate around cationic and polar residues, and least stable near hydrophobic residues, especially when more buried, with stability displaying moderate entropy-enthalpy compensation. Conversely, anionic residues in the proteins are energetically destabilized relative to singly solvated amino acids, while trends for other residues are less clear-cut. Almost all residues lose intraresidue entropy when in the protein, enthalpy changes are negative on average but may be positive or negative, and the resulting overall stability is moderate for some proteins and negligible for others. The free energy of water around single amino acids is found to closely match existing hydrophobicity scales. Regarding the effect of secondary structure, water is slightly more stable around loops, of intermediate stability around β strands and turns, and least stable around helices. An interesting asymmetry observed is that cationic residues stabilize a residue when bonded to its N-terminal side but destabilize it when on the C-terminal side, with a weaker reversed trend for anionic residues.

Harnessing the cyclization strategy for new drug discovery

The design of new ligands with high affinity and specificity against the targets of interest has been a central focus in drug discovery. As one of the most commonly used methods in drug discovery, the cyclization represents a feasible strategy to identify new lead compounds by increasing structural novelty, scaffold diversity and complexity. Such strategy could also be potentially used for the follow-on drug discovery without patent infringement. In recent years, the cyclization strategy has witnessed great success in the discovery of new lead compounds against different targets for treating various diseases. Herein, we first briefly summarize the use of the cyclization strategy in the discovery of new small-molecule lead compounds, including the proteolysis targeting chimeras (PROTAC) molecules. Particularly, we focus on four main strategies including fused ring cyclization, chain cyclization, spirocyclization and macrocyclization and highlight the use of the cyclization strategy in lead generation. Finally, the challenges including the synthetic intractability, relatively poor pharmacokinetics (PK) profiles and the absence of the structural information for rational structure-based cyclization are also briefly discussed. We hope this review, not exhaustive, could provide a timely overview on the cyclization strategy for the discovery of new lead compounds.

Decoding molecular recognition of inhibitors targeting HDAC2 via molecular dynamics simulations and configurational entropy estimation

Molecular recognition by enzymes is a complicated process involving thermodynamic energies governing protein-ligand interactions. In order to aid the estimation of inhibitory activity of compounds targeting an enzyme, several computational methods can be employed to dissect this intermolecular contact. Herein, we report a structural dynamics investigation of an epigenetic enzyme HDAC2 in differentiating its binding to various inhibitors within the sub-sites of its active site. Molecular dynamics (MD) simulation was employed to elucidate the intermolecular interactions as well as the dynamics behavior of ligand binding. MD trajectories of five distinct HDAC2-inhibitor complexes reveal that compounds lacking adequate contacts with the opening rim of the active site possess high fluctuation along the cap portion, thus weakening the overall affinity. Key intermolecular interactions determining the effective binding of inhibitors include hydrogen bonds with Gly154, Asp181, and Tyr308; hydrophobic interactions between Phe155/Phe210 and the linker region; and a pi-stacking with Arg39 at the foot pocket. Decomposition of the binding free energy calculated per-residue by MM/PBSA also indicates that the interactions within the internal foot pocket, especially with residues Met35, Leu144, Gly305, and Gly306, can contribute significantly to the ligand binding. Additionally, configurational entropy of the binding was estimated and compared to the scale of the binding free energy in order to assess its contribution to the binding and to differentiate various ligand partners. It was found that the levels of entropic contribution are comparable among a set of structurally similar carbamide ligands, while it is greatly different for the set of unrelated ligands, ranging from 2.75 to 16.38 kcal/mol for the five inhibitors examined. These findings exemplify the importance of assessing molecular dynamics as well as estimating the entropic contribution in evaluating the ligand binding mechanism.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Determination of immunogenic epitopes in major house dust mite allergen, Der p 2, via nanoallergens

BACKGROUND

Despite the high prevalence of allergic asthma, currently, avoidance of the responsible allergens, which is nearly impossible for allergens such as house dust mite (HDM), remains among the most effective treatment. Consequently, determination of the immunogenic epitopes of allergens will aid in developing a better understanding of the condition for diagnostic and therapeutic purposes. Current methods of epitope identification, however, only evaluate immunoglobulin E-epitope binding interactions, which is not directly related to epitope immunogenicity.

OBJECTIVE

To determine and rank the immunogenicity of the epitopes of major HDM allergen, Der p 2.

METHODS

We performed degranulation assays with RBL-SX38 cells primed using patient plasma and challenged with nanoallergens which multivalently displayed epitopes to study the relative immunogenicity of various epitopes of Der p 2. Nanoallergens were used to evaluate epitopes individually or in combination.

RESULTS

When evaluated using 3 patient samples, 3 epitopes in 2 distal regions of Der p 2 were identified as highly immunogenic when presented in combination, whereas no individual epitope triggered relevant degranulation. One of the epitopes (69-DPNACHYMKCPLVKGQQY-86) was identified to be cooperatively immunogenic when combined with other epitopes.

CONCLUSION

Our study highlights the importance of conformational epitopes in HDM-related allergies. This study also provides further evidence of the versatility of nanoallergens and their value for functional characterization of allergy epitopes, by ranking the Der p 2 epitopes according to immunogenicity. We believe that nanoallergens, by aiding in identification and understanding of immunogenic epitopes, will provide a better understanding of the manifestation of the allergic condition and potentially aid in developing new treatments.

Comparison of Grand Canonical and Conventional Molecular Dynamics Simulation Methods for Protein-Bound Water Networks

Water molecules play important roles in all biochemical processes. Therefore, it is of key importance to obtain information of the structure, dynamics, and thermodynamics of water molecules around proteins. Numerous computational methods have been suggested with this aim. In this study, we compare the performance of conventional and grand-canonical Monte Carlo (GCMC) molecular dynamics (MD) simulations to sample the water structure, as well GCMC and grid-based inhomogeneous solvation theory (GIST) to describe the energetics of the water network. They are evaluated on two proteins: the buried ligand-binding site of a ferritin dimer and the solvent-exposed binding site of galectin-3. We show that GCMC/MD simulations significantly speed up the sampling and equilibration of water molecules in the buried binding site, thereby making the results more similar for simulations started from different states. Both GCMC/MD and conventional MD reproduce crystal-water molecules reasonably for the buried binding site. GIST analyses are normally based on restrained MD simulations. This improves the precision of the calculated energies, but the restraints also significantly affect both absolute and relative energies. Solvation free energies for individual water molecules calculated with and without restraints show a good correlation, but with large quantitative differences. Finally, we note that the solvation free energies calculated with GIST are ∼5 times larger than those estimated by GCMC owing to differences in the reference state.

Exploring Glycan Binding Specificity of Odorranalectin by Alanine Scanning Library

Fluorescently labelled alanine scan analogues of odorranalectin (OL), a cyclic peptide that exhibits lectin like properties, were screened for binding BSA-conjugated monosaccharides using an enzyme-linked lectin assay (ELLA). Results revealed that Lys⁵, Phe⁷, Tyr⁹, Gly¹², Leu¹⁴, and Thr¹⁷ were crucial for binding BSA-L-fucose, BSA-D-galactose and BSA-N-acetyl-D-galactosamine. Notably, Ala substitution of Ser³, Pro⁴, and Val¹³ resulted in higher binding affinities compared to the native OL. The obtained data also indicated that Arg⁸ plays an important role in differentiation of binding for BSA-L-fucose/D-galactose from BSA-N-acetyl-D-galactosamine. The thermodynamics of binding of the selected alanine analogues was evaluated by isothermal titration calorimetry. Low to moderate binding affinities were determined for the tetravalent MUC1 glycopeptide and asialofetuin, respectively, and high for the fucose rich polysaccharide, fucoidan. The thermodynamic profile of interactions with asialofetuin exhibits shift to an entropy-driven mechanism compared to the fucoidan, which displayed an enthalpyentropy compensation, typically associated with the carbohydratelectin recognition process.

Spatially Resolved Hydration Thermodynamics in Biomolecular Systems

Water is essential for the structure, dynamics, energetics, and thus the function of biomolecules. It is a formidable challenge to elicit, in microscopic detail, the role of the solvation-related driving forces of biomolecular processes, such as the enthalpy and entropy contributions to the underlying free-energy landscape. In this Perspective, we discuss recent developments and applications of computational methods that provide a spatially resolved map of hydration thermodynamics in biomolecular systems and thus yield atomic-level insights to guide the interpretation of experimental observations. An emphasis is on the challenge of quantifying the hydration entropy, which requires characterization of both the motions of the biomolecules and of the water molecules in their surrounding.

Ligand binding remodels protein side-chain conformational heterogeneity

While protein conformational heterogeneity plays an important role in many aspects of biological function, including ligand binding, its impact has been difficult to quantify. Macromolecular X-ray diffraction is commonly interpreted with a static structure, but it can provide information on both the anharmonic and harmonic contributions to conformational heterogeneity. Here, through multiconformer modeling of time- and space-averaged electron density, we measure conformational heterogeneity of 743 stringently matched pairs of crystallographic datasets that reflect unbound/apo and ligand-bound/holo states. When comparing the conformational heterogeneity of side chains, we observe that when binding site residues become more rigid upon ligand binding, distant residues tend to become more flexible, especially in non-solvent-exposed regions. Among ligand properties, we observe increased protein flexibility as the number of hydrogen bonds decreases and relative hydrophobicity increases. Across a series of 13 inhibitor-bound structures of CDK2, we find that conformational heterogeneity is correlated with inhibitor features and identify how conformational changes propagate differences in conformational heterogeneity away from the binding site. Collectively, our findings agree with models emerging from nuclear magnetic resonance studies suggesting that residual side-chain entropy can modulate affinity and point to the need to integrate both static conformational changes and conformational heterogeneity in models of ligand binding.

P2X3-selective mechanism of Gefapixant, a drug candidate for the treatment of refractory chronic cough

•

The mechanism by which Gefapixant/AF-219 selectively acts on the P2X3 receptor is unclear.

•

The negative allosteric site of AF-219 at P2X3 is also a potent allosteric site for other P2X subtypes.

•

The selectivity of AF-219 for P2X3 is determined by the accessibility of binding site and the internal shape of this pocket.

•

The finding will provide new perspectives for drug design against P2X3-mediated diseases such as RCC.

Gefapixant/AF-219, a selective inhibitor of the P2X3 receptor, is the first new drug other than dextromethorphan to be approved for the treatment of refractory chronic cough (RCC) in nearly 60 years. To date, seven P2X subtypes (P2X1-7) activated by extracellular ATP have been cloned, and subtype selectivity of P2X inhibitors is a prerequisite for reducing side effects. We previously identified the site and mechanism of action of Gefapixant/AF-219 on the P2X3 receptor, which occupies a pocket consisting of the left flipper (LF) and lower body (LB) domains. However, the mechanism by which AF-219 selectively acts on the P2X3 receptor is unknown. Here, we combined mutagenesis, chimera construction, molecular simulations, covalent occupation and chemical synthesis, and find that the negative allosteric site of AF-219 at P2X3 is also present in other P2X subtypes, at least for P2X1, P2X2 and P2X4. By constructing each chimera of AF-219 sensitive P2X3 and insensitive P2X2 subtypes, the insensitive P2X2 subtype was made to acquire the inhibitory properties of AF-219 and AF-353, an analog of AF-219 with higher affinity. Our results suggest that the selectivity of AF-219/AF-353 for P2X3 over the other P2X subtypes is determined by a combination of the accessibility of P2X3 binding site and the internal shape of this pocket, a finding that could provide new perspectives for drug design against P2X3-mediated diseases such as RCC, idiopathic pulmonary fibrosis, hypertension and overactive bladder disorder.

Increased Molecular Flexibility Widens the Gap between K i and K d values in Screening for Retinoid X Receptor Modulators

Screening for small-molecule modulators targeting a particular receptor is frequently based on measurement of K _d, i.e., the binding constant between the receptor and the compound of interest. However, K _d values also reflect binding at receptor protein sites other than the modulatory site. We designed derivatives of retinoid X receptor (RXR) antagonist CBTF-EE (1) with modifications that altered their conformational flexibility. Compounds 6a,b and 7a,b showed quite similar K _d values, but 7a,b exhibited 10-fold higher K _i values than those of 6a,b. Further, 6a,b showed potent RXR-antagonistic activity, while 7a,b were inactive. These results suggest that increased conformational flexibility promotes binding at nontarget receptor sites. In this situation, conventional determination of K _d is less effective for screening purposes than the determination of K _i using a ligand that binds specifically to the site regulating transcriptional activity. Thus, the use of K _i values for orthosteric ligands may increase the hit rate in screening active regulatory molecules.

Double Mutant of Chymotrypsin Inhibitor 2 Stabilized through Increased Conformational Entropy

The conformational heterogeneity of a folded protein can affect not only its function but also stability and folding. We recently discovered and characterized a stabilized double mutant (L49I/I57V) of the protein CI2 and showed that state-of-the-art prediction methods could not predict the increased stability relative to the wild-type protein. Here, we have examined whether changed native-state dynamics, and resulting entropy changes, can explain the stability changes in the double mutant protein, as well as the two single mutant forms. We have combined NMR relaxation measurements of the ps-ns dynamics of amide groups in the backbone and the methyl groups in the side chains with molecular dynamics simulations to quantify the native-state dynamics. The NMR experiments reveal that the mutations have different effects on the conformational flexibility of CI2: a reduction in conformational dynamics (and entropy estimated from this) of the native state of the L49I variant correlates with its decreased stability, while increased dynamics of the I57V and L49I/I57V variants correlates with their increased stability. These findings suggest that explicitly accounting for changes in native-state entropy might be needed to improve the predictions of the effect of mutations on protein stability.

Toward a Computational NMR Procedure for Modeling Dipeptide Side-Chain Conformation

Theoretical relationships between the vicinal spin–spin coupling constants (SSCCs) and the χ₁ torsion angles have been studied to predict the conformations of protein side chains. An efficient computational procedure is developed to obtain the conformation of dipeptides through theoretical and experimental SSCCs, Karplus equations, and quantum chemistry methods, and it is applied to three aliphatic hydrophobic residues (Val, Leu, and Ile). Three models are proposed: unimodal-static, trimodal-static-stepped, and trimodal-static-trigonal, where the most important factors are incorporated (coupled nuclei, nature and orientation of the substituents, and local geometric properties). Our results are validated by comparison with NMR and X-ray empirical data described in the literature, obtaining successful results on the 29 residues considered. Using out trimodal residue treatment, it is possible to detect and resolve residues with a simple conformation and those with two or three staggered conformers. In four residues, a deeper analysis explains that they do not have a unique conformation and that the population of each conformation plays an important role.

Flexible CDOCKER: Hybrid Searching Algorithm and Scoring Function with Side Chain Conformational Entropy

The binding of small-molecule ligands to protein or nucleic acid targets is important to numerous biological processes. Accurate prediction of the binding modes between a ligand and a macromolecule is of fundamental importance in structure-based structure-function exploration. When multiple ligands with different sizes are docked to a target receptor, it is reasonable to assume that the residues in the binding pocket may adopt alternative conformations upon interacting with the different ligands. In addition, it has been suggested that the entropic contribution to binding can be important. However, only a few attempts to include the side chain conformational entropy upon binding within the application of flexible receptor docking methodology exist. Here, we propose a new physics-based scoring function that includes both enthalpic and entropic contributions upon binding by considering the conformational variability of the flexible side chains within the ensemble of docked poses. We also describe a novel hybrid searching algorithm that combines both molecular dynamics (MD)-based simulated annealing and genetic algorithm crossovers to address the enhanced sampling of the increased search space. We demonstrate improved accuracy in flexible cross-docking experiments compared with rigid cross-docking. We test our developments by considering five protein targets, thrombin, dihydrofolate reductase(DHFR), T4 L99A, T4 L99A/M102Q, and PDE10A, which belong to different enzyme classes with different binding pocket environments, as a representative set of diverse ligands and receptors. Each target contains dozens of different ligands bound to the same binding pocket. We also demonstrate that this flexible docking algorithm may be applicable to RNA docking with a representative riboswitch example. Our findings show significant improvements in top ranking accuracy across this set, with the largest improvement relative to rigid, 23.64%, occurring for ligands binding to DHFR. We then evaluate the ability to identify lead compounds among a large chemical space for the proposed flexible receptor docking algorithm using a subset of the DUD-E containing receptor targets MCR, GCR, and ANDR. We demonstrate that our new algorithms show improved performance in modeling flexible binding site residues compared to DOCK. Finally, we select the T4 L99A and T4 L99A/M102Q decoy sets, containing dozens of binders and experimentally validated nonbinders, to test our approach in distinguishing binders from nonbinders. We illustrate that our new algorithms for searching and scoring have superior performance to rigid receptor CDOCKER as well as AutoDock Vina. Finally, we suggest that flexible CDOCKER is sufficiently fast to be utilized in high-throughput docking screens in the context of hierarchical approaches.

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.

A Step toward NRF2-DNA Interaction Inhibitors by Fragment-Based NMR Methods

The NRF2 transcription factor is a key regulator in cellular oxidative stress response, and acts as a tumor suppressor. Aberrant activation of NRF2 has been implicated in promoting chemo-resistance, tumor growth, and metastasis by activating its downstream target genes. Hence, inhibition of NRF2 promises to be an attractive therapeutic strategy to suppress cell proliferation and enhance cell apoptosis in cancer. Direct targeting of NRF2 with small-molecules to discover protein-DNA interaction inhibitors is challenging as it is a largely intrinsically disordered protein. To discover molecules that bind to NRF2 at the DNA binding interface, we performed an NMR-based fragment screen against its DNA-binding domain. We discovered several weakly binding fragment hits that bind to a region overlapping with the DNA binding site. Using SAR by catalogue we developed an initial structure-activity relationship for the most interesting initial hit series. By combining NMR chemical shift perturbations and data-driven docking, binding poses which agreed with NMR information and the observed SAR were elucidated. The herein discovered NRF2 hits and proposed binding modes form the basis for future structure-based optimization campaigns on this important but to date 'undrugged' cancer driver.

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.

Observation of Sub-Microsecond Protein Methyl-Side Chain Dynamics by Nanoparticle-Assisted NMR Spin Relaxation

Amino-acid side-chain properties in proteins are key determinants of protein function. NMR spin relaxation of side chains is an important source of information about local protein dynamics and flexibility. However, traditional solution NMR relaxation methods are most sensitive to sub-nanosecond dynamics lacking information on slower ns-μs time-scale motions. Nanoparticle-assisted NMR spin relaxation (NASR) of methyl-side chains is introduced here as a window into these ns-μs dynamics. NASR utilizes the transient and nonspecific interactions between folded proteins and slowly tumbling spherical nanoparticles (NPs), whereby the increase of the relaxation rates reflects motions on time scales from ps all the way to the overall tumbling correlation time of the NPs ranging from hundreds of ns to μs. The observed motional amplitude of each methyl group can then be expressed by a model-free NASR S² order parameter. The method is demonstrated for ²H-relaxation of CH₂D methyl moieties and cross-correlated relaxation of CH₃ groups for proteins Im7 and ubiquitin in the presence of anionic silica-nanoparticles. Both types of relaxation experiments, dominated by either quadrupolar or dipolar interactions, yield highly consistent results. Im7 shows additional dynamics on the intermediate time scales taking place in a functionally important loop, whereas ubiquitin visits the majority of its conformational substates on the sub-ns time scale. These experimental observations are in good agreement with 4-10 μs all-atom molecular dynamics trajectories. NASR probes side-chain dynamics on a much wider range of motional time scales than previously possible, thereby providing new insights into the interplay between protein structure, dynamics, and molecular interactions that govern protein function.

A Computational Perspective on Molecular Recognition by Galectins

This article presents an overview of recent computational studies dedicated to the analysis of binding between galectins and small-molecule ligands. We first present a summary of the most popular simulation techniques adopted for calculating binding poses and binding energies, and then discuss relevant examples reported in the literature for the three main classes of galectins (dimeric, tandem and chimera). We show that simulation of galectin-ligand interactions is a mature field which has proven invaluable for completing and unraveling experimental observations. Future perspectives to further improve the accuracy and cost-effectiveness of existing computational approaches will involve the development of new schemes to account for solvation and entropy effects, which represent the main current limitations to the accuracy of computational results.

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.

Energy-entropy method using multiscale cell correlation to calculate binding free energies in the SAMPL8 host-guest challenge

Free energy drives a wide range of molecular processes such as solvation, binding, chemical reactions and conformational change. Given the central importance of binding, a wide range of methods exist to calculate it, whether based on scoring functions, machine-learning, classical or electronic structure methods, alchemy, or explicit evaluation of energy and entropy. Here we present a new energy-entropy (EE) method to calculate the host-guest binding free energy directly from molecular dynamics (MD) simulation. Entropy is evaluated using Multiscale Cell Correlation (MCC) which uses force and torque covariance and contacts at two different length scales. The method is tested on a series of seven host-guest complexes in the SAMPL8 (Statistical Assessment of the Modeling of Proteins and Ligands) "Drugs of Abuse" Blind Challenge. The EE-MCC binding free energies are found to agree with experiment with an average error of 0.9 kcal mol-1. MCC makes clear the origin of the entropy changes, showing that the large loss of positional, orientational, and to a lesser extent conformational entropy of each binding guest is compensated for by a gain in orientational entropy of water released to bulk, combined with smaller decreases in vibrational entropy of the host, guest and contacting water.

Structural Dynamics from NMR Relaxation by SRLS Analysis: Local Geometry, Potential Energy Landscapes, and Spectral Densities

We have developed the two-body coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics by nuclear magnetic resonance (NMR) relaxation. The rotators are represented by diffusion tensors D₁ for overall protein tumbling and D₂ for locally ordered probe motion. D₁ and D₂ are coupled dynamically by a potential, u, typically given by linear combinations of the Wigner functions D₀₀² and (D₀₂² + D_0-2²). Until now, our SRLS analyses provided the tensors, D₁ and D₂, the potential, u, and the geometric link between SRLS and NMR. Here we enhance this description by also examining the SRLS spectral densities obtained by solving the SRLS Smoluchowski equation. In addition, we show that the form of u specified above complies with two NMR-detected potential energy landscapes representing preferential ordering along N-H or C^α-C^α. Pictorial illustrations thereof are provided. The extended SRLS analysis is applied to ¹⁵N-H relaxation from the carbohydrate recognition domain of galectin-3 (Gal3C) in complex with two diastereomeric ligands, S and R. We find that D₂ is isotropic with a principal value, D₂, of 10¹⁰ s^-1 on average, and it is faster in the strands β₃, β₅, and β₈. The potential, u, is strong (∼20 kT); it is slightly rhombic when N-H is the main ordering axis and highly rhombic when C^α-C^α is the main ordering axis. Gal3C-S exhibits primarily preferential ordering along C^α-C^α; Gal3C-R exhibits both types of ordering. The binding-associated polypeptide chain segment of Gal3C-S is homogeneous, whereas that of Gal3C-R is diversified, with regard to D₂ and ordering preference. We associate these features with the previously determined diminished binding constant of Gal3C-R in comparison with Gal3C-S. Thus, the present study enhances the SRLS analysis, in general, and provides new insights into the dynamic structure and binding properties of Gal3C-S and Gal3C-R, in particular.

Structural Insights into the Molecular Recognition Mechanism of the Cancer and Pathogenic Epitope, LacdiNAc by Immune-Related Lectins

Interactions of glycan-specific epitopes to human lectin receptors represent novel immune checkpoints for investigating cancer and infection diseases. By employing a multidisciplinary approach that combines isothermal titration calorimetry, NMR spectroscopy, molecular dynamics simulations, and X-ray crystallography, we investigated the molecular determinants that govern the recognition of the tumour and pathogenic glycobiomarker LacdiNAc (GalNAcβ1-4GlcNAc, LDN), including their comparison with the ubiquitous LacNAc epitope (Galβ1-4GlcNAc, LN), by two human immune-related lectins, galectin-3 (hGal-3) and the macrophage galactose C-type lectin (hMGL). A different mechanism of binding and interactions was observed for the hGal-3/LDN and hMGL/LDN complexes, which explains the remarkable difference in the binding specificity of LDN and LN by these two lectins. The new structural clues reported herein are fundamental for the chemical design of mimetics targeting hGal-3/hMGL recognition process.

Dissecting the Conformational Free Energy of a Small Peptide in Solution

The free energy surface of a small peptide was analyzed based on an unbiased microsecond molecular dynamics simulation. The peptide sampled disordered conformational ensembles of distinct compactness, and its free energy was decomposed into separate contributions from the intramolecular potential energy, conformational entropy, and solvation free energy. The latter was further broken down into enthalpic and entropic contributions due to peptide-water and water-water interactions. This decomposition was enabled by a generalized linear response relation between the peptide-water interaction energy and the solvation free energy, which was empirically parametrized by explicit solvation free energy calculations for representative peptide conformations. This full dissection of the peptide free energy identifies individual contributions that stabilize and destabilize compact and extended peptide conformational ensembles and reveals the origin of a free energy barrier associated with transitions between them.

Applications of isothermal titration calorimetry in pure and applied research from 2016 to 2020

The last 5 years have seen a series of advances in the application of isothermal titration microcalorimetry (ITC) and interpretation of ITC data. ITC has played an invaluable role in understanding multiprotein complex formation including proteolysis-targeting chimeras (PROTACS), and mitochondrial autophagy receptor Nix interaction with LC3 and GABARAP. It has also helped elucidate complex allosteric communication in protein complexes like trp RNA-binding attenuation protein (TRAP) complex. Advances in kinetics analysis have enabled the calculation of kinetic rate constants from pre-existing ITC data sets. Diverse strategies have also been developed to study enzyme kinetics and enzyme-inhibitor interactions. ITC has also been applied to study small molecule solvent and solute interactions involved in extraction, separation, and purification applications including liquid-liquid separation and extractive distillation. Diverse applications of ITC have been developed from the analysis of protein instability at different temperatures, determination of enzyme kinetics in suspensions of living cells to the adsorption of uremic toxins from aqueous streams.

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.

Ligand Strain and Its Conformational Complexity Is a Major Factor in the Binding of Cyclic Dinucleotides to STING Protein

STING (stimulator of interferon genes) is a key regulator of innate immunity that has recently been recognized as a promising drug target. STING is activated by cyclic dinucleotides (CDNs) which eventually leads to expression of type I interferons and other cytokines. Factors underlying the affinity of various CDN analogues are poorly understood. Herein, we correlate structural biology, isothermal calorimetry (ITC) and computational modeling to elucidate factors contributing to binding of six CDNs-three pairs of natural (ribo) and fluorinated (2'-fluororibo) 3',3'-CDNs. X-ray structural analyses of six {STING:CDN} complexes did not offer any explanation for the different affinities of the studied ligands. ITC showed entropy/enthalpy compensation up to 25 kcal mol-1 for this set of similar ligands. The higher affinities of fluorinated analogues are explained with help of computational methods by smaller loss of entropy upon binding and by smaller strain (free) energy.