Molecular Determinants of the Kinetic Binding Properties of Antihistamines at the Histamine H1 Receptors

The binding affinity of ligands for their receptors is determined by their kinetic and thermodynamic binding properties. Kinetic analyses of the rate constants of association and dissociation (kon and koff, respectively) of antihistamines have suggested that second-generation antihistamines have a long duration of action owing to the long residence time (1/koff) at the H1 receptors. In this study, we examined the relationship between the kinetic and thermodynamic binding properties of antihistamines, followed by an evaluation of the structural determinants responsible for their kinetic binding properties using quantitative structure–activity relationship (QSAR) analyses. We found that whereas the binding enthalpy and entropy might contribute to the increase and decrease, respectively, in the koff values, there was no significant relationship with the kon values. QSAR analyses indicated that kon and koff values could be determined by the descriptors FASA_H (water-accessible surface area of all hydrophobic atoms divided by total water-accessible surface area) and vsurf_CW2 (a 3D molecular field descriptor weighted by capacity factor 2, the ratio of the hydrophilic surface to the total molecular surface), respectively. These findings provide further insight into the mechanisms by which the kinetic binding properties of antihistamines are regulated by their thermodynamic binding forces and physicochemical properties.

Differential Regulation of Bilastine Affinity for Human Histamine H1 Receptors by Lys 179 and Lys 191 via Its Binding Enthalpy and Entropy

Bilastine, a zwitterionic second-generation antihistamine containing a carboxyl group, has higher selectivity for H1 receptors than first-generation antihistamines. Ligand-receptor docking simulations have suggested that the electrostatic interaction between the carboxyl group of second-generation antihistamines and the amino group of Lys179ECL2 and Lys1915.39 of human H1 receptors might contribute to increased affinity of these antihistamines to H1 receptors. In this study, we evaluated the roles of Lys179ECL2 and Lys1915.39 in regulating the electrostatic and hydrophobic binding of bilastine to H1 receptors by thermodynamic analyses. The binding enthalpy and entropy of bilastine were estimated from the van ’t Hoff equation using the dissociation constants. These constants were obtained from the displacement curves against the binding of [3H] mepyramine to membrane preparations of Chinese hamster ovary cells expressing wild-type human H1 receptors and their Lys179ECL2 or Lys1915.39 mutants to alanine at various temperatures. We found that the binding of bilastine to wild-type H1 receptors occurred by enthalpy-dependent binding forces and, more dominantly, entropy-dependent binding forces. The mutation of Lys179ECL2 and Lys1915.39 to alanine reduced the affinity of bilastine to H1 receptors by reducing enthalpy- and entropy-dependent binding forces, respectively. These results suggest that Lys179ECL2 and Lys1915.39 differentially contribute to the increased binding affinity to bilastine via electrostatic and hydrophobic binding forces.

Does Antibody Stabilize the Ligand Binding in GP120 of HIV-1 Envelope Protein? Evidence from MD Simulation

CD4-mimetic HIV-1 entry inhibitors are small sized molecules which imitate similar conformational flexibility, in gp120, to the CD4 receptor. However, the mechanism of the conformational flexibility instigated by these small sized inhibitors is little known. Likewise, the effect of the antibody on the function of these inhibitors is also less studied. In this study, we present a thorough inspection of the mechanism of the conformational flexibility induced by a CD4-mimetic inhibitor, NBD-557, using Molecular Dynamics Simulations and free energy calculations. Our result shows the functional importance of Asn425 in substrate induced conformational dynamics in gp120. The MD simulations of Asn425Gly mutant provide a less dynamic gp120 in the presence of NBD-557 without incapacitating the binding enthalpy of NBD-557. The MD simulations of complexes with the antibody clearly show the enhanced affinity of NBD-557 due to the presence of the antibody, which is in good agreement with experimental Isothermal Titration Calorimetry results (Biochemistry2006, 45, 10973–10980).

SARS-CoV-2 Fusion Peptide has a Greater Membrane Perturbating Effect than SARS-CoV with Highly Specific Dependence on Ca2+

Coronaviruses are a major infectious disease threat, and include the zoonotic-origin human pathogens SARS-CoV-2, SARS-CoV, and MERS-CoV (SARS-2, SARS-1, and MERS). Entry of coronaviruses into host cells is mediated by the spike (S) protein. In our previous ESR studies, the local membrane ordering effect of the fusion peptide (FP) of various viral glycoproteins including the S of SARS-1 and MERS has been consistently observed. We previously determined that the sequence immediately downstream from the S2’ cleavage site is the bona fide SARS-1 FP. In this study, we used sequence alignment to identify the SARS-2 FP, and studied its membrane ordering effect. Although there are only three residue difference, SARS-2 FP induces even greater membrane ordering than SARS-1 FP, possibly due to its greater hydrophobicity. This may be a reason that SARS-2 is better able to infect host cells. In addition, the membrane binding enthalpy for SARS-2 is greater. Both the membrane ordering of SARS-2 and SARS-1 FPs are dependent on Ca2+, but that of SARS-2 shows a greater response to the presence of Ca2+. Both FPs bind two Ca2+ ions as does SARS-1 FP, but the two Ca2+ binding sites of SARS-2 exhibit greater cooperativity. This Ca2+ dependence by the SARS-2 FP is very ion-specific. These results show that Ca2+ is an important regulator that interacts with the SARS-2 FP and thus plays a significant role in SARS-2 viral entry. This could lead to therapeutic solutions that either target the FP-calcium interaction or block the Ca2+ channel.

Experimental Characterization of the Association of Nine Novel Cyclodextrin Derivatives with Two Guest Compounds

We investigate the binding of native β-cyclodextrin (β-CD) and eight novel β-CD derivatives with two different guest compounds, using isothermal calorimetry (ITC) and 2D NOESY NMR. In all cases, the stoichiometry is 1:1 and binding is exothermic. Overall, modifications at the 3’ position of β-CD, which is at the secondary face, weaken binding by several kJ/mol relative to native β-CD, while modifications at the 6’ position (primary face) maintain or somewhat reduce the binding affinity. The variations in binding enthalpy are larger than the variations in binding free energy, so entropy-enthalpy compensation is observed. Characterization of the bound conformations with NOESY NMR shows that the polar groups of the guests may be situated at either face, depending on the host molecule, and, in some cases, both orientations are populated. The present results were used in the SAMPL7 blinded prediction challenge whose results are detailed in the same special issue of JCAMD.

Experimental Characterization of the Association of Nine Novel Cyclodextrin Derivatives with Two Guest Compounds

We investigate the binding of native β-cyclodextrin (β-CD) and eight novel β-CD derivatives with two different guest compounds, using isothermal calorimetry (ITC) and 2D NOESY NMR. In all cases, the stoichiometry is 1:1 and binding is exothermic. Overall, modifications at the 3’ position of β-CD, which is at the secondary face, weaken binding by several kJ/mol relative to native β-CD, while modifications at the 6’ position (primary face) maintain or somewhat reduce the binding affinity. The variations in binding enthalpy are larger than the variations in binding free energy, so entropy-enthalpy compensation is observed. Characterization of the bound conformations with NOESY NMR shows that the polar groups of the guests may be situated at either face, depending on the host molecule, and, in some cases, both orientations are populated. The present results were used in the SAMPL7 blinded prediction challenge whose results are detailed in the same special issue of JCAMD.

Dapsone is not a Pharmacodynamic Lead Compound for its Aryl Derivatives

Background: The relatedness between the linear equations of thermodynamics and QSAR was studied thanks to the recently elucidated crystal structure complexes between sulfonamide pterin conjugates and dihydropteroate synthase (DHPS) together with a published set of thirty- six synthetic dapsone derivatives with their reported entropy-driven activity data. Only a few congeners were slightly better than dapsone. Objective : Our study aimed at demonstrating the applicability of thermodynamic QSAR and to shed light on the mechanistic aspects of sulfone binding to DHPS. Methods: To this end ligand docking to DHPS, quantum mechanical properties, 2D- and 3D-QSAR as well as Principle Component Analysis (PCA) were carried out. Results: The short aryl substituents of the docked pterin-sulfa conjugates were outward oriented into the solvent space without interacting with target residues which explains why binding enthalpy (ΔH) did not correlate with potency. PCA revealed how chemically informative descriptors are evenly loaded on the first three PCs (interpreted as ΔG, ΔH and ΔS), while chemically cryptic ones reflected higher dimensional (complex) loadings. Conclusions: It is safe to utter that synthesis efforts to introduce short side chains for aryl derivatization of the dapsone scaffold have failed in the past. On theoretical grounds we provide computed evidence why dapsone is not a pharmacodynamic lead for drug profiling because enthalpic terms do not change significantly at the moment of ligand binding to target.

Simulation of Intrinsic Defects and Cd Site Occupation in LaIn3 and LuIn3

In previous work, perturbed angular correlation spectroscopy (PAC) was used to determine jump rates of 111Cd, the daughter of the 111In radiotracer, in the series of phases RIn3 (R = rare-earth element) through nuclear quadrupole relaxation. Greater relaxation, indicating faster Cd jump rates, was observed in heavy rare-earths for compositions more deficient in indium, as would be expected for diffusion mediated by vacancies on the In sublattice. On the other hand, greater relaxation was observed for light rare-earths (R = La, Ce, and Pr) for compositions with excess indium, suggesting Cd diffusion is mediated there by a different mechanism. In this work, computer simulations were carried out to better understand the nature of the relaxation observed for the light rare-earths and the origin of the change in behavior across the rare-earth series. As a first step, formation enthalpies of intrinsic defects were calculated using density functional theory (DFT) for series end-members LaIn3 and LuIn3. Both compounds were found to exhibit Schottky thermal disorder. Additional DFT simulations show that the binding enthalpy between In-and R-vacancies is larger in LaIn3 than in LuIn3, suggesting that diffusion in LaIn3 might be mediated by divacancies. Site enthalpies of Cd also were calculated, and it was found more favorable energetically for Cd to occupy the In sublattice than the R sublattice in both end-member phases.