Exploring the Origins of Enthalpy–Entropy Compensation by Calorimetric Studies of Cyclodextrin Complexes

Interactions of a Novel Anthracycline with Oligonucleotide DNA and Cyclodextrins in an Aqueous Environment

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

ITC analysis of polydisperse systems: Unravelling the impact of sample heterogeneity

Binding interactions often involve heterogeneous samples displaying a distribution of binding sites that vary in affinity and binding enthalpy. Examples include biological samples like proteins and chemically produced samples like modified cyclodextrins. Experimental studies often ignore sample heterogeneity and treat the system as an interaction of two homogeneous species, i.e. a chemically well-defined ligand binding to one type of site. The present study explores, by simulations and experiments, the impact of heterogeneity in isothermal titration calorimetry (ITC) setups where one of the binding components is heterogeneous. It is found that the standard single-site model, based on the assumption of two homogeneous binding components, provides excellent fits to simulated ITC data when the binding free energy is normally distributed and all sites have similar binding enthalpies. In such cases, heterogeneity can easily go undetected but leads to underestimated binding constants. Heterogeneity in the binding enthalpy is a bigger problem and may result in enthalpograms of increased complexity that are likely to be misinterpreted as two-site binding or other complex binding models. Finally, it is shown that heterogeneity can account for previously observed experimental anomalies. All simulations are accessible in Google Colab for readers to experiment with the simulation parameters.

The temperature-dependence of host-guest binding thermodynamics: experimental and simulation studies

The thermodynamic parameters of host-guest binding can be used to describe, understand, and predict molecular recognition events in aqueous systems. However, interpreting binding thermodynamics remains challenging, even for these relatively simple molecules, as they are determined by both direct and solvent-mediated host-guest interactions. In this contribution, we focus on the contributions of water to binding by studying binding thermodynamics, both experimentally and computationally, for a series of nearly rigid, electrically neutral host-guest systems and report the temperature-dependent thermodynamic binding contributions ΔGb(T), ΔHb(T), ΔSb(T), and ΔCp,b. Combining isothermal titration calorimetry (ITC) measurements with molecular dynamics (MD) simulations, we provide insight into the binding forces at play for the macrocyclic hosts cucurbit[n]uril (CBn, n = 7-8) and β-cyclodextrin (β-CD) with a range of guest molecules. We find consistently negative changes in heat capacity on binding (ΔCp,b) for all systems studied herein - as well as for literature host-guest systems - indicating increased enthalpic driving forces for binding at higher temperatures. We ascribe these trends to solvation effects, as the solvent properties of water deteriorate as temperature rises. Unlike the entropic and enthalpic contributions to binding, with their differing signs and magnitudes for the classical and non-classical hydrophobic effect, heat capacity changes appear to be a unifying and more general feature of host-guest complex formation in water. This work has implications for understanding protein-ligand interactions and other complex systems in aqueous environments.

The role of loop dynamics in the prediction of ligand-protein binding enthalpy

The enthalpic and entropic components of ligand-protein binding free energy reflect the interactions and dynamics between ligand and protein. Despite decades of study, our understanding and hence our ability to predict these individual components remains poor. In recent years, there has been substantial effort and success in the prediction of relative and absolute binding free energies, but the prediction of the enthalpic (and entropic) contributions in biomolecular systems remains challenging. Indeed, it is not even clear what kind of performance in terms of accuracy could currently be obtained for such systems. It is, however, relatively straight-forward to compute the enthalpy of binding. We thus evaluated the performance of absolute enthalpy of binding calculations using molecular dynamics simulation for ten inhibitors against a member of the bromodomain family, BRD4-1, against isothermal titration calorimetry data. Initial calculations, with the AMBER force-field showed good agreement with experiment (R² = 0.60) and surprisingly good accuracy with an average of root-mean-square error (RMSE) = 2.49 kcal mol^-1. Of the ten predictions, three were obvious outliers that were all over-predicted compared to experiment. Analysis of various simulation factors, including parameterization, buffer concentration and conformational dynamics, revealed that the behaviour of a loop (the ZA loop on the periphery of the binding site) strongly dictates the enthalpic prediction. Consistent with previous observations, the loop exists in two distinct conformational states and by considering one or the other or both states, the prediction for the three outliers can be improved dramatically to the point where the R² = 0.95 and the accuracy in terms of RMSE improves to 0.90 kcal mol^-1. However, performance across force-fields is not consistent: if OPLS and CHARMM are used, different outliers are observed and the correlation with the ZA loop behaviour is not recapitulated, likely reflecting parameterization as a confounding problem. The results provide a benchmark standard for future study and comparison.

Molecular Recognition in an Aqueous Medium Using Water-Soluble Prismarene Hosts

The synthesis of water-soluble prism[n]arenes (n = 5 and 6) bearing anionic carboxylato groups on the rims is described. The binding properties of this novel class of water-soluble hosts are studied by nuclear magnetic resonance and calorimetry. The complexation of singly and doubly charged ammonium guests with the more rigid pentamer is enthalpically driven by secondary interactions, while the thermodynamic fingerprint for the larger hexamer reveals driving forces that strongly depend on the guest charge and/or size.

Tuning DNA Supramolecular Polymers by the Addition of Small, Functionalized Nucleobase Mimics

Nucleobase mimicking small molecules able to reconfigure DNA are a recently discovered strategy that promises to extend the structural and functional diversity of nucleic acids. However, only simple, unfunctionalized molecules such as cyanuric acid and melamine have so far been used in this approach. In this work, we show that the addition of substituted cyanuric acid molecules can successfully program polyadenine strands to assemble into supramolecular fibers. Unlike conventional DNA nanostructure functionalization, which typically end-labels DNA strands, our approach incorporates functional groups into DNA with high density using small molecules and results in new DNA triple helices coated with alkylamine or alcohol units that grow into micrometer-long fibers. We find that small changes in the small molecule functional group can result in large structural and energetic variation in the overall assembly. A combination of circular dichroism, atomic force microscopy, molecular dynamics simulations, and a new thermodynamic method, transient equilibrium mapping, elucidated the molecular factors behind these large changes. In particular, we identify substantial DNA sugar and phosphate group deformations to accommodate a hydrogen bond between the phosphate and the small-molecule functional groups, as well as a critical chain length of the functional group which switches this interaction from intra- to interfiber. These parameters allow the controlled formation of hierarchical, hybrid DNA assemblies simply through the addition and variation of small, functionalized molecules.

Thermodynamic Studies of Interactions between Sertraline Hydrochloride and Randomly Methylated β-Cyclodextrin Molecules Supported by Circular Dichroism Spectroscopy and Molecular Docking Results

The interaction between sertraline hydrochloride (SRT) and randomly methylated β-cyclodextrin (RMβCD) molecules have been investigated at 298.15 K under atmospheric pressure. The method used-Isothermal Titration Calorimetry (ITC) enabled to determine values of the thermodynamic functions like the enthalpy (ΔH), the entropy (ΔS) and the Gibbs free energy (ΔG) of binding for the examined system. Moreover, the stoichiometry coefficient of binding (n) and binding/association constant (K) value have been calculated from the experimental results. The obtained outcome was compared with the data from the literature for other non-ionic βCD derivatives interacting with SRT and the enthalpy-entropy compensation were observed and interpreted. Furthermore, the connection of RMβCD with SRT was characterized by circular dichroism spectroscopy (CD) and complexes of βCD derivatives with SRT were characterized through the computational studies with the use of molecular docking (MD).

Cyclodextrin Complexation as a Way of Increasing the Aqueous Solubility and Stability of Carvedilol

We studied the effect of several CDs on carvedilol's solubility and chemical stability in various aqueous media. Our present results show that it is possible to achieve a carvedilol concentration of 5 mg/mL (12.3 mM) in the presence of 5 eq of γCD or RAMEB in an aqueous medium with an acceptable acid pH (between 3.5 and 4.7). Carvedilol formed 1:1 inclusion complexes but those with RAMEB appear to be stronger (K = 317 M-1 at 298 K) than that with γCD (K = 225 M-1 at 298 K). The complexation of carvedilol by RAMEB significantly increased the drug's photochemical stability in aqueous solution. These results might constitute a first step towards the development of a novel oral formulation of carvedilol.

Interaction of Polyelectrolytes with Proteins: Quantifying the Role of Water

A theoretical model is presented for the free energy ΔGb of complex formation between a highly charged polyelectrolyte and a protein. The model introduced here comprises both the effect of released counterions and the uptake or release of water molecules during complex formation. The resulting expression for ΔGb is hence capable of describing the dependence of ΔGb on temperature as well as on the concentration of salt in the system: An increase of the salt concentration in the solution increases the activity of the ions and counterion release becomes less effective for binding. On the other hand, an increased salt concentration leads to the decrease of the activity of water in bulk. Hence, release of water molecules during complex formation will be more advantageous and lead to an increase of the magnitude of ΔGb and the binding constant. It is furthermore demonstrated that the release or uptake of water molecules is the origin of the marked enthalpy-entropy cancellation observed during complex formation of polyelectrolytes with proteins. The comparison with experimental data on complex formation between a synthetic (sulfated dendritic polyglycerol) and natural polyelectrolytes (DNA; heparin) with proteins shows full agreement with theory.

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.

A physicochemical, thermodynamical, structural and computational evaluation of kynurenic acid/cyclodextrin complexes

In this work, the interaction between Kynurenic acid (KYNA) and several natural and modified cyclodextrins (CDs) is carried out. Among all the CD tested, HPβ-CD showed the strongest complexation constant (K_F), with a value of 270.94 ± 29.80 M^-1. Between natural (α- and β-) CDs, the complex of KYNA with β-CD was the most efficient. The inclusion complex of KYNA with CDs showed a strong influence of pH and temperature. The K_F value decreased at high pH values, when the pK_a was passed. Moreover, an increase of the temperature caused a decrease in the K_F values. The thermodynamic parameters of the complexation (ΔH°, ΔS° and ΔG°) were studied with negative entropy, enthalpy and spontaneity of the process at 25 °C. Moreover, the inclusion complex was also characterized using FTIR and TGA. Finally, molecular docking calculations provided different interactions and their influence in the complexation constant.

Understanding the Interaction of Polyelectrolyte Architectures with Proteins and Biosystems

The counterions neutralizing the charges on polyelectrolytes such as DNA or heparin may dissociate in water and greatly influence the interaction of such polyelectrolytes with biomolecules, particularly proteins. In this Review we give an overview of studies on the interaction of proteins with polyelectrolytes and how this knowledge can be used for medical applications. Counterion release was identified as the main driving force for the binding of proteins to polyelectrolytes: Patches of positive charge become multivalent counterions of the polyelectrolyte and lead to the release of counterions from the polyelectrolyte and a concomitant increase in entropy. This is shown from investigations on the interaction of proteins with natural and synthetic polyelectrolytes. Special emphasis is paid to sulfated dendritic polyglycerols (dPGS). The Review demonstrates that we are moving to a better understanding of charge-charge interactions in systems of biological relevance. Research along these lines will aid and promote the design of synthetic polyelectrolytes for medical applications.

Evaluating the Performance of Water Models with Host-Guest Force Fields in Binding Enthalpy Calculations for Cucurbit[7]uril-Guest Systems

Computational prediction of thermodynamic components with computational methods has become increasingly routine in computer-aided drug design. Although there has been significant recent effort and improvements in the calculation of free energy, the prediction of enthalpy (and entropy) remains underexplored. Furthermore, there has been relatively little work reported so far that attempts to comparatively assess how well different force fields and water models perform in conjunction with each other. Here, we report a comprehensive assessment of force fields and water models using host-guest systems that mimic many features of protein-ligand systems. These systems are computationally inexpensive, possibly because of their small size compared to protein-ligand systems. We present absolute enthalpy calculations using the multibox approach on a set of 25 cucurbit[7]uril-guest pairs. Eight water models were considered (TIP3P, TIP4P, TIP4P-Ew, SPC, SPC/E, OPC, TIP5P, Bind3P), along with five force fields commonly used in the literature (GAFFv1, GAFFv2, CGenFF, Parsley, and SwissParam). We observe that host-guest binding enthalpies are strongly sensitive to the selection of force field and water model. In terms of water models, we find that TIP3P and its derivative Bind3P are the best performing models for this particular host-guest system. The performance is generally better for aliphatic compounds than for aromatic ones, suggesting that aromaticity remains a difficult property to include accurately in these simple force fields.

Adamantane Functionalized Poly(2-oxazoline)s with Broadly Tunable LCST-Behavior by Molecular Recognition

Smart or adaptive materials often utilize stimuli-responsive polymers, which undergo a phase transition in response to a given stimulus. So far, various stimuli have been used to enable the modulation of drug release profiles, cell-interactive behavior, and optical and mechanical properties. In this respect, molecular recognition is a powerful tool to fine-tune the stimuli-responsive behavior due to its high specificity. Within this contribution, a poly(2-oxazoline) copolymer bearing adamantane side chains was synthesized via triazabicyclodecene-catalyzed amidation of the ester side chains of a poly(2-ethyl-2-oxazoline-stat-2-methoxycarbonylpropyl-2-oxazoline) statistical copolymer. Subsequent complexation of the pendant adamantane groups with sub-stoichiometric amounts (0-1 equivalents) of hydroxypropyl β-cyclodextrin or β-cyclodextrin enabled accurate tuning of its lower critical solution temperature (LCST) over an exceptionally wide temperature range, spanning from 30 °C to 56 °C. Furthermore, the sharp thermal transitions display minimal hysteresis, suggesting a reversible phase transition of the complexed polymer chains (i.e., the β-cyclodextrin host collapses together with the polymers) and a minimal influence by the temperature on the supramolecular association. Analysis of the association constant of the polymer with hydroxypropyl β-cyclodextrin via 1H NMR spectroscopy suggests that the selection of the macrocyclic host and rational polymer design can have a profound influence on the observed thermal transitions.

Comment on "Physicochemical stimuli as tuning parameters to modulate the structure and stability of nanostructured lipid carriers and release kinetics of encapsulated antileprosy drugs" by R. Kanwar, M. Gradzielski, S. Prevost, G. Kaur, M. S. Appavou and S. K. Mehta, J. Mater. Chem. B, 2019, , 6539

In a recent article [R. Kanwar et al., J. Mater. Chem. B, 2019, 7(42), 6539-6555], the authors characterized the interactions between drug-loaded nanostructured lipid carriers and bovine serum albumin using thermodynamics. They found that the interactions are spontaneous and driven by entropy. In this present paper, we report our analysis of these results in terms of equilibrium thermodynamics to show that the binding reactions exhibit enthalpy-entropy compensation. Our findings may prove useful for designing nanostructured lipid carriers.