Understanding noncovalent interactions: ligand binding energy and catalytic efficiency from ligand-induced reductions in motion within receptors and enzymes

Thermodynamic entropy of organic oxidation in the water environment: experimental evaluation compared to semi-empirical calculation

Residual organic matters in the secondary effluent are usually less biodegradable in terms of the total organic carbon content, and when discharged into a receiving water body, their further decomposition most likely mainly occurs due to chemical oxidation. Using this scenario, a semi-empirical method was previously developed to calculate the thermodynamic entropy of organic oxidation to quantitatively evaluate the impact of organic discharge on the water environment. In this study, the relationship between the entropy increase (ΔS_C) and excess organic mass (ΔTOC) was experimentally verified via combustion heat measurement using typical organic chemicals and mixtures. For individual organic chemicals, a linear relationship was detected between ΔS_C and ΔTOC with the same proportionality coefficient, 54.0 kJ/g, determined in the previous semi-empirical relationship. For the organic mixtures, a linear relationship was also identified; however, the proportionality coefficient was 69.2 kJ/g, indicating an approximately 28 % increase in the oxidation heat required to decompose the same organic mass. This increase in energy can likely be attributed to the synergistic effects of hydrogen bonding, hydrophobic interactions, π-π interactions, and van der Waals interactions between functional groups of different organic compounds. Intermolecular interactions may result in 17-32 % more dissociation energy for organic mixtures compared to the organic components' chemical structures. Because organics discharged into a water body are always a mixture of organic compounds, the proportionality coefficient obtained using organic mixtures should be adopted to modify the previously proposed semi-empirical equation.

Investigation of a Catenane with a Responsive Noncovalent Network: Mimicking Long-Range Responses in Proteins

We report a functional synthetic model for studying the noncovalent networks (NCNs) required for complex protein functions. The model [2]-catenane is self-assembled from dipeptide building blocks and contains an extensive network of hydrogen bonds and aromatic interactions. Perturbations to the catenane cause compensating changes in the NCNs structure and dynamics, resulting in long-distance changes reminiscent of a protein. Key findings include the notion that NCNs require regions of negative cooperativity, or "frustrated" noncovalent interactions, as a source of potential energy for driving the response. We refer to this potential energy as latent free energy and describe a mechanistic and energetic model for responsive systems.

Direct inference of site strength in basic solids upon CO2 adsorption: enthalpy-entropy compensation effects

The adsorption of CO2 coupled to calorimetry is a state-of-the-art technique for characterizing the basic properties of solids. In this paper, we show that the differential heat and entropy curves measured upon CO2 adsorption on a basic solid can be reasonably estimated from a single CO2 isotherm with no need for any independent heat (calorimetric) measurement. Our method relies on two important observations: (1) formulation of generalized F-H-TS thermodynamic isotherms, the former (F) being directly generated from the raw CO2 isotherms, and (2) the presence of unexpected enthalpy-entropy compensation effects upon CO2 adsorption linking the integral enthalpy and entropy of adsorption until saturation for different solids. Our thermodynamic method has been validated using a broad library of basic solids with variable site strength and heterogeneity. Finally, a new scale of basicity is proposed using the parameters fitted from the thermodynamic isotherm (free energy basis) as descriptors of basic strength. This method opens an avenue to the inference of site strength of basic solids without the need for expensive calorimeters.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

An Evolving Insight into Chiral H-Bond Catalyzed Aza-Henry Reactions: A Cooperative Role for Noncovalent Attractive Interactions Unveiled by Density Functional Theory

The role of cooperative effects arising from noncovalent attractive interactions as a vital factor governing stereoinduction in chiral H-bond catalyzed aza-Henry reactions is reported. Supporting this finding were density functional theory (DFT) calculations which revealed a shape and size dependency existed between the catalyst and substrates that when matched lead to high enantioselectivity, as reflected by favorable activation parameters. Associated with optimal catalyst and substrate pairing were a closed catalytic binding pocket and a synclinal orientation of the substrates that reinforced favorable stereoelectronic effects and dispersive type forces. Meanwhile, unfavorable steric interactions were found to be a dominant effect controlling diastereoselection.

Hydrocarbon constrained peptides - understanding preorganisation and binding affinity

The development of constrained peptides represents an emerging strategy to generate peptide based probes and hits for drug-discovery that address challenging protein-protein interactions (PPIs). In this manuscript we report on the use of a novel α-alkenylglycine derived amino acid to synthesise hydrocarbon constrained BH3-family sequences (BIM and BID). Our biophysical and structural analyses illustrate that whilst the introduction of the constraint increases the population of the bioactive α-helical conformation of the peptide in solution, it does not enhance the inhibitory potency against pro-apoptotic Bcl-x_L and Mcl-1 PPIs. SPR analyses indicate binding occurs via an induced fit mechanism whilst X-ray analyses illustrate none of the key interactions between the helix and protein are disturbed. The behaviour derives from enthalpy-entropy compensation which may be considered in terms of the ground state energies of the unbound constrained and unconstrained peptides; this has implications for the design of preorganised peptides to target protein-protein interactions.

Challenges and Opportunities in the Development of Aptamers for TNFα

RNA aptamers for tumor necrosis factor-alpha (TNFα), for which functionality was demonstrated in L929 cells, show only little affinity for the protein in vitro. Detailed investigation of the aptamer-protein interaction by surface plasmon resonance and quartz crystal microbalance analysis revealed that affinity is not the only crucial parameter for efficacy and functionality of those aptamers. Instead, the sensitive equilibrium of the monomeric and homotrimeric form of soluble TNFα decides on aptamer binding. Our results show that the field of application and the source of TNFα have to be carefully defined before selection of aptamer sequences.

Triangulating Nucleic Acid Conformations Using Multicolor Surface Energy Transfer

Optical ruler methods employing multiple fluorescent labels offer great potential for correlating distances among several sites, but are generally limited to interlabel distances under 10 nm and suffer from complications due to spectral overlap. Here we demonstrate a multicolor surface energy transfer (McSET) technique able to triangulate multiple points on a biopolymer, allowing for analysis of global structure in complex biomolecules. McSET couples the competitive energy transfer pathways of Förster Resonance Energy Transfer (FRET) with gold-nanoparticle mediated Surface Energy Transfer (SET) in order to correlate systematically labeled points on the structure at distances greater than 10 nm and with reduced spectral overlap. To demonstrate the McSET method, the structures of a linear B-DNA and a more complex folded RNA ribozyme were analyzed within the McSET mathematical framework. The improved multicolor optical ruler method takes advantage of the broad spectral range and distances achievable when using a gold nanoparticle as the lowest energy acceptor. The ability to report distance information simultaneously across multiple length scales, short-range (10-50 Å), mid-range (50-150 Å), and long-range (150-350 Å), distinguishes this approach from other multicolor energy transfer methods.

Enhancing binding affinity and selectivity through preorganization and cooperative enhancement of the receptor

When direct host–guest binding interactions are weakened by unfavorable solvent competition, guest-triggered intrareceptor interactions could be used to augment the binding.

An unusual CC⋯CO interaction in (Z)-3-[(4-halogenphenyl)amino]-2-cyanoprop-2-enoates

An unusual CC⋯CO interaction has been discovered in (Z)-3-[(4-halogenphenyl)amino]-2-cyanoprop-2-enoates and rationalized by the density functional theory calculations.

Efficiency parameters in artificial allosteric systems

It is shown that the conformational energy ΔG_C required for the formation of a suitable cavity or cleft for ligand binding in absence of an effector E is the essential factor for the efficiency of synthetic allosteric systems with positive cooperativity.

Assembly of Lock-and-Key Colloids Mediated by Polymeric Depletant

Polymer-mediated lock-and-key assembly via depletion attraction is purely a shape-recognition process without any molecular bonding. Since the depletion attraction relates to osmotic pressure and excluded volume, the binding tendency in a dispersion of lock-and-key colloids can be controlled by adjusting the characteristics of polymeric depletants. In this work, dissipative particle dynamics accounting for explicit solvents, polymers, and multiple lock-key pairs are performed to investigate the influences of the polymer concentration, chain length, solvent quality, and chain stiffness. As the polymer concentration associated with osmotic pressure is increased, the binding free energy (E(b)) between a lock-key pair rises linearly and the binding fraction (θ(LK)) in the dispersion grows sigmoidally. Moreover, the increases in the chain length, solvent quality, and chain stiffness lead to the expansion of the polymer size associated with excluded volume and thus both E(b) and θ(LK) rise accordingly. However, E(b) and θ(LK) grow to be insensitive to the chain length for long enough polymer coils but still can be enhanced if the polymer becomes rod-like. As the solvent quality is varied, θ(LK) can be dramatically altered, although the radius of gyration of polymers is slightly changed.

Probing zeolites by vibrational spectroscopies

This review addresses the most relevant aspects of vibrational spectroscopies (IR, Raman and INS) applied to zeolites and zeotype materials. Surface Brønsted and Lewis acidity and surface basicity are treated in detail. The role of probe molecules and the relevance of tuning both the proton affinity and the steric hindrance of the probe to fully understand and map the complex site population present inside microporous materials are critically discussed. A detailed description of the methods needed to precisely determine the IR absorption coefficients is given, making IR a quantitative technique. The thermodynamic parameters of the adsorption process that can be extracted from a variable-temperature IR study are described. Finally, cutting-edge space- and time-resolved experiments are reviewed. All aspects are discussed by reporting relevant examples. When available, the theoretical literature related to the reviewed experimental results is reported to support the interpretation of the vibrational spectra on an atomic level.

Globular and disordered-the non-identical twins in protein-protein interactions

In biology proteins from different structural classes interact across and within classes in ways that are optimized to achieve balanced functional outputs. The interactions between intrinsically disordered proteins (IDPs) and other proteins rely on changes in flexibility and this is seen as a strong determinant for their function. This has fostered the notion that IDP's bind with low affinity but high specificity. Here we have analyzed available detailed thermodynamic data for protein-protein interactions to put to the test if the thermodynamic profiles of IDP interactions differ from those of other protein-protein interactions. We find that ordered proteins and the disordered ones act as non-identical twins operating by similar principles but where the disordered proteins complexes are on average less stable by 2.5 kcal mol(-1).

Terpyridine-Cu(ii) targeting human telomeric DNA to produce highly stereospecific G-quadruplex DNA metalloenzyme

The cofactors commonly involved in natural enzymes have provided the inspiration for numerous advances in the creation of artificial metalloenzymes. Nevertheless, to design an appropriate cofactor for a given biomolecular scaffold or vice versa remains a challenge in developing efficient catalysts in biochemistry. Herein, we extend the idea of G-quadruplex-targeting anticancer drug design to construct a G-quadruplex DNA metalloenzyme. We found that a series of terpyridine-Cu(ii) complexes (CuLn) can serve as excellent cofactors to dock with human telemetric G-quadruplex DNA. The resulting G-quadruplex DNA metalloenzyme utilising CuL1 catalyzes an enantioselective Diels-Alder reaction with enantioselectivity of >99% enantiomeric excess and about 73-fold rate acceleration compared to CuL1 alone. The terpyridine-Cu(ii) complex cofactors demonstrate dual functions, both as an active site to perform catalysis and as a structural regulator to promote the folding of human telemetric G-quadruplex DNA towards excellent catalysts.

Enthalpy-driven nuclease-like activity and mechanism of peptide-chlorambucil conjugates

We report the results of attaching the anticancer drug chlorambucil (CLB) to two high-affinity DNA binding peptides: Met-Hyp-Arg-Lys-(Py)4-Lys-Arg-NH2 (HyM-10) and Gln-Hyp-Arg-Lys-(Py)4-Lys-Arg-NH2 (HyQ-10). These CLB-peptide conjugates cleave DNA very effectively and sequence-selectively without the use of chemicals, heat, or UV irradiation. Polyacrylamide gel electrophoresis identifies the sites where CLB-HyM-10 and CLB-HyQ-10 attack a complementary pair of 5'-(32)P-labeled duplexes derived from pBR322 in the absence of piperidine or other chemical additives. Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS) has confirmed the preferential cleavage sites as well as a novel stepwise cleavage mechanism of sequence-selective DNA cleavage. Resembling restriction endonucleases, the CLB-peptide conjugates appear to be capable of producing double strand DNA breaks. Circular dichroism studies show that CLB-HyM-10 and CLB-HyQ-10 induce significant local conformational changes in DNA via the minor groove, possibly with dimeric binding stoichiometry. The energetic basis of DNA binding by these conjugates has been investigated by isothermal titration calorimetry, revealing that the binding of both the peptides and their CLB conjugates is overwhelmingly enthalpy-driven. The maintenance of a conserved negative binding free energy in DNA-conjugate interactions is a crucial feature of the universal enthalpy-entropy compensation phenomenon. The strongly enthalpy-driven binding of CLB-peptide conjugates to preferred loci in DNA furnishes the required proximity effect to generate the observed nuclease-like sequence-selective cleavage.

Collisional unfolding of multiprotein complexes reveals cooperative stabilization upon ligand binding

Cooperative binding mechanisms are a common feature in biology, enabling a diverse range of protein-based molecular machines to regulate activities ranging from oxygen uptake to cellular membrane transport. Much, however, is not known about such cooperative binding mechanisms, including how such events typically add to the overall stability of such protein systems. Measurements of such cooperative stabilization events are challenging, as they require the separation and resolution of individual protein complex bound states within a mixture of potential stoichiometries to individually assess protein stabilities. Here, we report ion mobility-mass spectrometry results for the concanavalin A tetramer bound to a range of polysaccharide ligands. We use collision induced unfolding, a relatively new methodology that functions as a gas-phase analog of calorimetry experiments in solution, to individually assess the stabilities of concanavalin A bound states. By comparing the differences in activation voltage required to unfold different concanavalin A-ligand stoichiometries, we find evidence suggesting a cooperative stabilization of concanavalin A occurs upon binding most carbohydrate ligands. We critically evaluate this observation by assessing a broad range of ligands, evaluating the unfolding properties of multiple protein charge states, and by comparing our gas-phase results with those obtained from calorimetry experiments carried out in solution.

Ionic interactions. Subnanoscale hydrophobic modulation of salt bridges in aqueous media

Polar interactions such as electrostatic forces and hydrogen bonds play an essential role in biological molecular recognition. On a protein surface, polar interactions occur mostly in a hydrophobic environment because nonpolar amino acid residues cover ~75% of the protein surface. We report that ionic interactions on a hydrophobic surface are modulated by their subnanoscale distance to the surface. We developed a series of ionic head groups-appended self-assembled monolayers with C2, C6, C8, and C12 space-filling alkyl chains, which capture a dendritic guest via the formation of multiple salt bridges. The guest release upon protonolysis is progressively suppressed when its distance from the background hydrophobe changes from 1.2 (C2) to 0.2 (C12) nanometers, with an increase in salt bridge strength of ~3.9 kilocalories per mole.

Binding energies of tyrosine kinase inhibitors: Error assessment of computational methods for imatinib and nilotinib binding

The binding energies of imatinib and nilotinib to tyrosine kinase have been determined by quantum mechanical (QM) computations, and compared with literature binding energy studies using molecular mechanics (MM). The potential errors in the computational methods include these critical factors. Errors in X-ray structures such as structural distortions and steric clashes give unrealistically high van der Waals energies, and erroneous binding energies.MM optimization gives a very different configuration to the QM optimization for nilotinib, whereas the imatinib ion gives similar configurations. Solvation energies are a major component of the overall binding energy. The QM based solvent model (PCM/SMD) gives different values from those used in the implicit PBSA solvent MM models. A major error in inhibitor—kinase binding lies in the non-polar solvation terms. Solvent transfer free energies and the required empirical solvent accessible surface area factors for nilotinib and imatinib ion to give the transfer free energies have been reverse calculated. These values differ from those used in the MM PBSA studies.An intertwined desolvation—conformational binding selectivity process is a balance of thermodynamic desolvation and intramolecular conformational kinetic control.The configurational entropies (TΔS) are minor error sources.

Dominant Alcohol-Protein Interaction via Hydration-Enabled Enthalpy-Driven Binding Mechanism

Water plays an important role in weak associations of small drug molecules with proteins. Intense focus has been on binding-induced structural changes in the water network surrounding protein binding sites, especially their contributions to binding thermodynamics. However, water is also tightly coupled to protein conformations and dynamics, and so far little is known about the influence of water-protein interactions on ligand binding. Alcohols are a type of low-affinity drugs, and it remains unclear how water affects alcohol-protein interactions. Here, we present alcohol adsorption isotherms under controlled protein hydration using in situ NMR detection. As functions of hydration level, Gibbs free energy, enthalpy, and entropy of binding were determined from the temperature dependence of isotherms. Two types of alcohol binding were found. The dominant type is low-affinity nonspecific binding, which is strongly dependent on temperature and the level of hydration. At low hydration levels, this nonspecific binding only occurs above a threshold of alcohol vapor pressure. An increased hydration level reduces this threshold, with it finally disappearing at a hydration level of h ≈ 0.2 (g water/g protein), gradually shifting alcohol binding from an entropy-driven to an enthalpy-driven process. Water at charged and polar groups on the protein surface was found to be particularly important in enabling this binding. Although further increase in hydration has smaller effects on the changes of binding enthalpy and entropy, it results in a significant negative change in Gibbs free energy due to unmatched enthalpy-entropy compensation. These results show the crucial role of water-protein interplay in alcohol binding.

Novel α-L-arabinofuranosidase from Cellulomonas fimi ATCC 484 and its substrate-specificity analysis with the aid of computer

In the process of gene mining for novel α-L-arabinofuranosidases (AFs), the gene Celf_3321 from Cellulomonas fimi ATCC 484 encodes an AF, termed as AbfCelf, with potent activity, 19.4 U/mg under the optimum condition, pH 6.0 and 40 °C. AbfCelf can hydrolyze α-1,5-linked oligosaccharides, sugar beet arabinan, linear 1,5-α-arabinan, and wheat flour arabinoxylan, which is partly different from some previously well-characterized GH 51 AFs. The traditional substrate-specificity analysis for AFs is labor-consuming and money costing, because the substrates include over 30 kinds of various 4-nitrophenol (PNP)-glycosides, oligosaccharides, and polysaccharides. Hence, a preliminary structure and mechanism based method was applied for substrate-specificity analysis. The binding energy (ΔG, kcal/mol) obtained by docking suggested the reaction possibility and coincided with the experimental results. AbfA crystal 1QW9 was used to test the rationality of docking method in simulating the interaction between enzyme and substrate, as well the credibility of the substrate-specificity analysis method in silico.

Dynamics of nitric oxide synthase-calmodulin interactions at physiological calcium concentrations

The intracellular Ca²⁺ concentration is an important regulator of many cellular functions. The small acidic protein calmodulin (CaM) serves as a Ca²⁺ sensor and control element for many enzymes. Nitric oxide synthase (NOS) is one of the proteins that is activated by CaM and plays a major role in a number of key physiological and pathological processes. Previous studies have shown CaM to act like a switch that causes a conformational change in NOS to allow for the electron transfer between the reductase and oxygenase domains through a process that is thought to be highly dynamic. We have analyzed the structure and dynamics of complexes formed by peptides based on inducible NOS (iNOS) and endothelial NOS (eNOS) with CaM at Ca²⁺ concentrations that mimic the physiological basal (17 and 100 nM) and elevated levels (225 nM) found in mammalian cells using fluorescence techniques and nuclear magnetic resonance spectroscopy. The results show the CaM-NOS complexes have similar structures at physiological and fully saturated Ca²⁺ levels; however, their dynamics are remarkably different. At 225 nM Ca²⁺, the CaM-NOS complexes show overall an increase in backbone dynamics, when compared to the dynamics of the complexes at saturating Ca²⁺ concentrations. Specifically, the N-lobe of CaM in the CaM-iNOS complex displays a lower internal mobility (higher S²) and higher exchange protection compared to those of the CaM-eNOS complex. In contrast, the C-lobe of CaM in the CaM-eNOS complex is less dynamic. These results illustrate that structures of CaM-NOS complexes determined at saturated Ca²⁺ concentrations cannot provide a complete picture because the differences in intramolecular dynamics become visible only at physiological Ca²⁺ levels.

Insight into the structural stability of coumestrol with human estrogen receptor α and β subtypes: a combined approach involving docking and molecular dynamics simulation studies

Structural stability of coumestrol with human estrogen receptor α and β subtypes.

A QSAR study on the inhibition mechanism of matrix metalloproteinase-12 by arylsulfone analogs based on molecular orbital calculations

The inhibition mechanism of matrix metalloproteinase-12 by arylsulfone analogs is revealed using a comprehensive computational approach including docking simulations, molecular orbital calculations, and QSAR.

Insight into the effect of inhibitor resistant S130G mutant on physico-chemical properties of SHV type beta-lactamase: a molecular dynamics study

Bacterial resistance is a serious threat to human health. The production of β-lactamase, which inactivates β-lactams is most common cause of resistance to the β-lactam antibiotics. The Class A enzymes are most frequently encountered among the four β-lactamases in the clinic isolates. Mutations in class A β-lactamases play a crucial role in substrate and inhibitor specificity. SHV and TEM type are known to be most common class A β-lactamases. In the present study, we have analyzed the effect of inhibitor resistant S130G point mutation of SHV type Class-A β-lactamase using molecular dynamics and other in silico approaches. Our study involved the use of different in silico methods to investigate the affect of S130G point mutation on the major physico-chemical properties of SHV type class A β-lactamase. We have used molecular dynamics approach to compare the dynamic behaviour of native and S130G mutant form of SHV β-lactamase by analyzing different properties like root mean square deviation (RMSD), H-bond, Radius of gyration (Rg) and RMS fluctuation of mutation. The results clearly suggest notable loss in the stability of S130G mutant that may further lead to decrease in substrate specificity of SHV. Molecular docking further indicates that S130G mutation decreases the binding affinity of all the three inhibitors in clinical practice.

Dynamic changes in composition of extracts of natural products as monitored by in situ NMR

The direct in situ NMR observation and quantification, based on the aldehyde -CH chemical shift region, of the inter-conversion of secoiridoid derivatives due to temperature and solvent effects is demonstrated in complex extracts of natural products without prior isolation of the individual components. The equilibrium between the aldehyde hydrate form and the dialdehyde form of the oleuropein aglycon of an olive leaf aqueous extract in D(2)O was shown to be temperature dependent. The resulting thermodynamic values of the Van't Hoff plot with ΔH(o)  = -26.34 ± 1.00 kJ mol(-1) and TΔS° (298 K) = -24.70 ± 1.00 kJ mol(-1) demonstrate a significant entropy term which nearly compensates the effect of enthalpy at room temperature. The equilibrium between the two diastereomeric hemiacetal forms and the dialdehyde form of the oleuropein 6-O-β-d-glucopyranoside aglycon of an olive leaf aqueous extract in CD(3) OD was also shown to be strongly temperature dependent again because of the significant entropy term (TΔS° (298 K) = -26.50 ± 1.39 kJ mol(-1)) compared with that of the enthalpy term (ΔH(o)  = -36.64 ± 1.46 kJ mol(-1)). This is the first demonstration of the significant role of the entropy parameter in determining the equilibrium of chemical transformations in complex mixtures of natural products due to solvent and temperature effects.

Thermodynamics of binding of structurally similar ligands to histone deacetylase 8 sheds light on challenges in the rational design of potent and isozyme-selective inhibitors of the enzyme

Among the different histone deacetylase (HDAC) isozymes, HDAC8 is the most highly malleable enzyme, and it exhibits the potential to accommodate structurally diverse ligands (albeit with moderate binding affinities) in its active site pocket. To probe the molecular basis of this feature, we performed detailed thermodynamic studies of the binding of structurally similar ligands, which differed with respect to the “cap”, “linker”, and “metal-binding” regions of the suberoylanilide hydroxamic acid (SAHA) pharmacophore, to HDAC8. The experimental data revealed that although the enthalpic (ΔH°) and entropic (ΔS°) changes for the binding of individual SAHA analogues to HDAC8 were substantially different, their binding free energies (ΔG°) were markedly similar, conforming to a strong enthalpy–entropy compensation effect. This effect was further observed in the temperature-dependent thermodynamics of binding of all SAHA analogues to the enzyme. Notably, in contrast to other metalloenzymes, our isothermal titration calorimetry experiments (performed in different buffers of varying ionization enthalpies) suggest that depending on the ligand, its zinc-binding group may or may not be deprotonated upon the binding to HDAC8. Furthermore, the heat capacity changes (ΔC_p°) associated with the ligand binding to HDAC8 markedly differed from one SAHA analogue to the other, and such features could primarily be rationalized in light of the dynamic flexibility in the enzyme structure in conjunction with the reorganization of the active site resident water molecules. Arguments are presented that although the binding thermodynamic features described above would facilitate identification of weak to moderately tight-binding HDAC8 inhibitors (by a high-throughput and/or virtual screening of libraries of small molecules), they would pose major challenges for the structure-based rational design of highly potent and isozyme-selective inhibitors of human HDAC8.

The Molecular and Structural Basis of HBV-resistance to Nucleos(t)ide Analogs

Infection with hepatitis B virus (HBV) is a worldwide health problem. Chronic hepatitis B can lead to fibrosis, liver cirrhosis, and hepatocellular carcinoma (HCC). Management of the latter two conditions often requires liver transplantation. Treatment with conventional interferon or pegylated interferon alpha can clear the virus, but the rates are very low. The likelihood, however, of viral resistance to interferon is minimal. The main problems with this therapy are the frequency and severity of side effects. In contrast, nucleos(t)ide analogs (NAs) have significantly lower side effects, but require long term treatment as sustained virological response rates are extremely low. However, long term treatment with NAs increases the risk for the development of anti-viral drug resistance. Only by understanding the molecular basis of resistance and using agents with multiple sites of action can drugs be designed to optimally prevent the occurrence of HBV antiviral resistance.

A B3LYP and MP2(full) theoretical investigation into cooperativity effects, aromaticity and thermodynamic properties in the Na(+)⋯benzonitrile⋯H2O ternary complex

The cooperativity effects between H-bonding and Na(+)⋯π or Na(+)⋯σ interactions in Na(+)⋯benzonitrile⋯H2O complexes were investigated using the B3LYP and MP2(full) methods with 6-311++G(2d,p) and aug-cc-pVTZ basis sets. The thermodynamic cooperativity and the influence of this cooperativity on aromaticity was evaluated by nucleus-independent chemical shifts (NICS). The results showed that the influence of the Na(+)⋯σ or Na(+)⋯π interaction on the hydrogen bond is more pronounced than that of the latter on the former. The cooperativity effect appeared in the Na(+)⋯σ interaction complex while the anti-cooperativity effect tended to be in the Na(+)⋯π system. The change in enthalpy is the major factor driving cooperativity. Thermodynamic cooperativity is not in accordance with the cooperativity effect evaluated by the change of interaction energy. The ring aromaticity of is weakened while the bond dissociation energy (BDE) of the C-CN bond increases upon ternary complex formation. The cooperativity effect (E coop) correlates with R c (NICS(1)ternary/NICS(1)binary) and ΔΔδ (Δδ ternary - Δδ binary) involving the ring and C ≡ N bond, as well as R BDE(C-CN) [BDE(C-CN)ternary/BDE(C-CN)binary], respectively. AIM (atoms in molecules) analysis confirms the existence of cooperativity.