The agonistic binding site at the histamine H2 receptor. I. Theoretical investigations of histamine binding to an oligopeptide mimicking a part of the fifth transmembrane α-helix

Summary Mutation studies on the histamine H₂ receptor were reported by Gantz et al. [J. Biol. Chem., 267 (1992) 20840], which indicate that both the mutation of the fifth transmembrane Asp¹⁸⁶ (to Ala¹⁸⁶) alone or in combination with Thr¹⁹⁰ (to Ala¹⁹⁰) maintained, albeit partially, the cAMP response to histamine. Recently, we have shown that histamine binds to the histamine H₂ receptor as a monocation in its proximal tautomeric form, and, moreover, we suggested that a proton is donated from the receptor towards the tele-position of the agonist, thereby triggering the biological effect [Nederkoorn et al., J. Mol. Graph., 12 (1994) 242; Eriks et al., Mol. Pharmacol., 44 (1993) 886]. These findings result in a close resemblance with the catalytic triad (consisting of Ser, His and Asp) found in serine proteases. Thr¹⁹⁰ resembles a triad's serine residue closely, and could also act as a proton donor. However, the mutation of Thr¹⁹⁰ to Ala¹⁹⁰ — the latter is unable to function as a proton donor — does not completely abolish the agonistic cAMP response. At the fifth transmembrane -helix of the histamine H₂ receptor near the extracellular surface, another amino acid is present, i.e. Tyr¹⁸², so an alternative couple of amino acids, Tyr¹⁸² and Asp¹⁸⁶, could constitute the histamine binding site at the fifth -helix instead of the (mutated) couple Asp¹⁸⁶ and Thr¹⁹⁰. In the first part of our present study, this hypothesis is investigated with the aid of an oligopeptide with an -helical backbone, which represents a part of the fifth transmembrane helix. Both molecular mechanics and ab initio data lead to the conclusion that the Tyr¹⁸²/Asp¹⁸⁶ couple is most likely to act as the binding site for the imidazole ring present in histamine.     

“In silico” study of the binding of two novel antagonists to the nociceptin receptor

Antagonists of the nociceptin receptor (NOP) are raising interest for their possible clinical use as antidepressant drugs. Recently, the structure of NOP in complex with some piperidine-based antagonists has been revealed by X-ray crystallography. In this study, a multi-flexible docking (MF-docking) procedure, i.e. docking to multiple receptor conformations extracted by preliminary molecular dynamics trajectories, together with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations have been carried out to provide the binding mode of two novel NOP antagonists, one of them selective (BTRX-246040, formerly named LY-2940094) and one non selective (AT-076), i.e. able to inactivate NOP as well as the classical µ- k- and δ-opioid receptors (MOP KOP and DOP). According to our results, the pivotal role of residue D130^3,32 (upper indexes are Ballesteros–Weinstein notations) is analogous to that enlighten by the already known X-ray structures of opioid receptors: binding of the molecules are predicted to require a slight readjustment of the hydrophobic pocket (residues Y131^3,33, M134^3,36, I219^5,43, Q280^6,52 and V283^6,55) in the orthosteric site of NOP, accommodating either the pyridine–pyrazole (BTRX-246040) or the isoquinoline (AT-076) moiety of the ligand, in turn allowing the protonated piperidine nitrogen to maximize interaction (salt-bridge) with residue D130^3,32 of the NOP, and the aromatic head to be sandwiched in optimal π-stacking between Y131^3,33 and M134^3,36. The QM/MM optimization after the MF-docking procedure has provided the more likely conformations for the binding to the NOP receptor of BTRX-246040 and AT-076, based on different pharmacophores and exhibiting different selectivity profiles. While the high selectivity for NOP of BTRX-246040 can be explained by interactions with NOP specific residues, the lack of selectivity of AT-076 could be associated to its ability to penetrate into the deep hydrophobic pocket of NOP, while retaining a conformation very similar to the one assumed by the antagonist JDTic into the K-opioid receptor. The proposed binding geometries fit better the binding pocket environment providing clues for experimental studies aimed to design selective or multifunctional opioid drugs.     

Free energy calculations offer insights into the influence of receptor flexibility on ligand–receptor binding affinities

Docking algorithms for computer-aided drug discovery and design often ignore or restrain the flexibility of the receptor, which may lead to a loss of accuracy of the relative free enthalpies of binding. In order to evaluate the contribution of receptor flexibility to relative binding free enthalpies, two host–guest systems have been examined: inclusion complexes of α-cyclodextrin (αCD) with 1-chlorobenzene (ClBn), 1-bromobenzene (BrBn) and toluene (MeBn), and complexes of DNA with the minor-groove binding ligands netropsin (Net) and distamycin (Dist). Molecular dynamics simulations and free energy calculations reveal that restraining of the flexibility of the receptor can have a significant influence on the estimated relative ligand–receptor binding affinities as well as on the predicted structures of the biomolecular complexes. The influence is particularly pronounced in the case of flexible receptors such as DNA, where a 50% contribution of DNA flexibility towards the relative ligand–DNA binding affinities is observed. The differences in the free enthalpy of binding do not arise only from the changes in ligand–DNA interactions but also from changes in ligand–solvent interactions as well as from the loss of DNA configurational entropy upon restraining.     

Importance of inter-residue interactions in ligand—receptor binding

In the present study, the role of inter-residue interactions in ligand binding and the ligand—receptor interactions were examined. Computational chemistry methods of ligand docking and molecular dynamics simulations were used to study the binding of β-funaltrexamine (β-FNA) and N-methyl-β-funaltrexamine (N-methyl-β-FNA) to μ- and κ-opioid receptors and to the μ-receptor with Lys303^6.58Glu mutation. It was found that inter-residue interactions Lys233^5.39—Glu303^6.58 in the mutant receptor and Lys227^5.39—Asp223^5.35 in the κ-receptor are more likely to prevent covalent bond formation between β-FNA and the receptor than the ligand-receptor interactions. This emphasizes the importance of inter-residue interactions in ligand binding as well as the effects of point-mutations.     

Chemoenzymatic synthesis and Fcγ receptor binding of homogeneous glycoforms of antibody Fc domain. Presence of a bisecting sugar moiety enhances the affinity of Fc to FcγIIIa receptor

Structurally well-defined IgG-Fc glycoforms are highly demanded for understanding the effects of glycosylation on an antibody's effector functions. We report in this paper chemoenzymatic synthesis and Fcγ receptor binding of an array of homogeneous IgG-Fc glycoforms. The chemoenzymatic approach consists of the chemical synthesis of defined N-glycan oxazolines as donor substrates, the expression of the Fc domain in a CHO cell line in the presence of an α-mannosidase inhibitor kifunensine, and an endoglycosidase-catalyzed glycosylation of the deglycosylated Fc domain (GlcNAc-Fc homodimer) with the synthetic glycan oxazolines. The enzyme from Arthrobacter protophormiae (Endo-A) was found to be remarkably efficient to take various modified N-glycan core oxazolines, including the bisecting sugar-containing derivatives, for Fc glycosylation remodeling, resulting in the formation of the corresponding homogeneous Fc glycoforms. Nevertheless, neither Endo-A nor the Mucor hiemalis endoglycosidase mutants (EndoM-N175A and EndoM-N175Q) were able to transfer full-length complex-type N-glycan to the Fc domain, implicating the limitations of these two enzymes in Fc glycosylation remodeling. Surface plasmon resonance (SPR) binding studies with the synthetic IgG-Fc glycoforms unambiguously proved that the presence of a bisecting GlcNAc moiety could significantly enhance the binding of Fc to FcγRIIIa, the activating Fcγ receptor, independent of Fc core-fucosylation. Interestingly, the Fc glycoforms carrying an unusual bisecting sugar moiety such as a mannose or a LacNAc moiety also demonstrated enhanced affinity to FcγRIIIa. On the orther hand, the presence of a bisecting GlcNAc or core-fucosylation had little effect on the affinity of Fc to the inhibitory Fcγ receptor, FcγRIIb. Our experimental data also showed that the α-linked mannose residues in the pentasaccharide Man3GlcNAc2 core was essential to maintain a high affinity of Fc to both FcγRIIIa and FcγRIIb. The synthetic homogeneous Fc glycoforms thus provide a useful tool for elucidating how a fine Fc N-glycan structure precisely affects the function of the Fc domain.     

Electrospray mass spectrometry for the direct accurate mass measurement of ligands in complex with the retinoid X receptor α ligand binding domain

Accurate mass measurements are often used in the structural determination of unknown compounds of low molecular mass (i.e., below approximately 500 Da). Recently, it has been shown that accurate mass measurements also can be made on small denatured proteins (i.e., M(r), approximately 17,000) to confirm their amino acid composition and identify the presence of isoforms. In the current report, we present nondenaturing electrospray (ES) mass spectrometry data on the direct accurate mass measurement of ligands in complex with the retinoid X receptor ligand binding domain (RXR LBD; M(r) 31,370.92). Average mass errors were below 0.198 Da, 6.3 ppm (standard deviation [SD], 0.146; n = 10) for low-affinity fatty acid agonists analyzed in complex with the RXR LBD. Protein consumption was less than 15 pmol, with fatty acid ligands present at concentrations corresponding to their median effective concentration value (low micromolar, determined in transfection assays). Although determination of fatty acid mass was only sufficiently accurate to give nominal mass values, measurements were of sufficient accuracy to assign fatty acid chain length, degree of unsaturation, or cyclization. Using 17beta-estradiol as a control, the ability to observe specific ligand binding is shown for both high- and low-affinity RXRalpha agonists. In addition, binding of a novel synthetic receptor agonist XCT0315908 to the RXRalpha LBD is reported. This compound showed a high degree of complex formation, and the receptor-ligand complex could be mass measured with an average mass error of -0.024 Da, 0.8 ppm (SD, 0.092; n = 9). Thus, specific binding of both nanomolar and micromolar affinity ligands to a nuclear receptor LBD can be directly observed using nondenaturing ES mass spectrometry and accurate mass measurements additionally can be made on intact complexes in the same experiment. This methodology also is applicable when ligands are present as components of mixtures.     

Effects of additives on the aggregation of surfactants: binding domination

The effects of electrolytes, alcohols, and urea on the aggregation of SDS, CTAB, and TritonX-100 at 25℃ have been investigated by fluorescence probing of pyrene. Both electrolytes and alcoholsreduce the critical micelle concentration (CMC) of the ionic surfactants, while the effect of the former ismore pronounced. It is shown that the effects of electrolytes mainly depend on the concentrationsand especially the valence of the opposite charge ions, and only slightly depend on the same charge ionsin respect of ion aggregate of micelle. The logarithm of CMC is not linearly correlated with theconcentrations of the counter ion or the electrolytes. The results are rationalized in terms of Hartley'smodel. Propanol increases the CMC of TX-100, while electrolytes and urea do not. In all the threekinds of surfactant micelles the excitation spectrum of pyrene slightly red-shifts (ca. 4 nm) from thatin water, but is not affected by the additives. The micropolarity of the environment in which pyrenemolecule resides in SDS micelle decreases with the increase of the concentrations of electrolytes. Thisis not the case when alcohols and urea were added to SDS or to TX-100. It is suggested that theaddition of electrolytes would result in more orderly orientation of SDS molecules. It is the bindingstrength of the counter ions that dominates the effects of additives on the aggregation of surfactants.     

Quantitative analysis of cation binding to the adenosine nucleotides using the variable ionic strength method: validation of the Debye-Hückel-Onsager theory of electrophoresis in the absence of counterion binding

The free solution mobilities of the adenosine nucleotides 5'-adenosine triphosphate (ATP), 5'-adenosine diphosphate (ADP), 5'-adenosine monophosphate (AMP), and 3'-5'-cyclic AMP (cAMP) have been measured in diethylmalonate buffers containing a wide variety of monovalent cations. The mobilities of all nucleotides increase gradually with the increase in intrinsic conductivity of the cation in the BGE. However, at a given conductivity, the mobilities observed for ATP, ADP, and AMP in BGEs containing alkali metal ions and other cations are lower than these observed in BGEs containing tetraalkylammonium ions. Since the mobility of cAMP is independent of the cation in the BGE, the results suggest that the relatively low mobilities observed for ATP, ADP, and AMP in BGEs containing cations other than a tetraalkylammonium ion are due to cation binding, reducing the effective net charge of the nucleotide and thereby reducing the observed mobility. To measure the binding quantitatively, the mobilities of the nucleotides were measured as a function of ionic strength. The mobilities of ATP, ADP, and AMP decrease nonlinearly with the square root of ionic strength (I(1/2)) in BGEs containing an alkali metal ion or Tris(+). By contrast, the mobilities decrease linearly with I(1/2) in BGEs containing a nonbinding quaternary ammonium ion, as expected from Debye-Hückel-Onsager (DHO) theory. The mobility of cAMP, a nonbinding analyte, decreases linearly with I(1/2), regardless of the cation in the BGE. Hence, a nonlinear decrease of the mobility of an analyte with I(1/2) appears to be a hallmark of counterion binding. The curved mobility profiles observed for ATP, ADP, and AMP in BGEs containing an alkali metal ion or Tris(+) were analyzed by nonlinear curve fitting, using difference mobility profiles to correct for the effect of the physical properties of BGE on the observed mobilities. The calculated apparent dissociation constants range from 22 to 344 mM, depending on the particular cation-nucleotide pair. Similar values have been obtained by other investigators, using different methods. Interestingly, Tris(+) and Li(+) bind to the adenosine nucleotides with approximately equal affinities, suggesting that positively charged Tris(+) buffer ions can compete with alkali metal ions in Tris-buffered solutions.     

Variations in binding among several agonists at two stoichiometries of the neuronal, α4β2 nicotinic receptor

Drug-receptor binding interactions of four agonists, ACh, nicotine, and the smoking cessation compounds varenicline (Chantix) and cytisine (Tabex), have been evaluated at both the 2:3 and 3:2 stoichiometries of the α4β2 nicotinic acetylcholine receptor (nAChR). Previous studies have established that unnatural amino acid mutagenesis can probe three key binding interactions at the nAChR: a cation-π interaction, and two hydrogen-bonding interactions to the protein backbone of the receptor. We find that all drugs make a cation-π interaction to TrpB of the receptor. All drugs except ACh, which lacks an N(+)H group, make a hydrogen bond to a backbone carbonyl, and ACh and nicotine behave similarly in acting as a hydrogen-bond acceptor. However, varenicline is not a hydrogen-bond acceptor to the backbone NH that interacts strongly with the other three compounds considered. In addition, we see interesting variations in hydrogen bonding interactions with cytisine that provide a rationalization for the stoichiometry selectivity seen with this compound.     

Conformational enantiomerization and estrogen receptor α binding of anti-cancer drug tamoxifen and its derivatives

The anticancer drug tamoxifen (TAM) displays two chiral vinyl propeller structures, which interconvert so rapidly that the process is undetectable on the NMR time scale. In the present work, the enantiomerization processes were investigated with molecular modeling techniques. The threshold mechanisms probed at the different rings were shown to be identical, i.e., involving a synchronous three-ring flip, with a correlated rotation of the rings. In order to reveal the pharmacological profiles of the two chiral forms, we performed structural studies on the ligand binding domain of estrogen receptor α (ERα LBD) and associated ligands. The enantiomers, with opposite torsional twist, were found to be discriminated by ERα. For TAM and its main metabolites, the effects of the stereoselectivity of ERα are overcome by the low energy cost for helical inversion between the two torsional enantiomers, estimated to be ～3 kcal/mol.     

Quantitative structure-activity relationship of compounds binding to estrogen receptor β based on heuristic method

Estrogen compounds may pose a serious threat to the health of humans and wildlife.The estrogen receptor (ER) exists as two subtypes,ERβ and ERβ.Compounds might have different relative affinities and binding modes for ERβ and ERβ.In this study,the heuristic method was performed on 31 compounds binding to ERβ to select 5 variances most related to the activity (LogRBA) from 1524 variances,which were then employed to develop the best model with the significant correlation and the best predictive power (r2 = 0.8...     

Influences of the Hydrophobicity of the Heme－binding Pocket on the Properties and Functions of Cytochrome b5 Mutants

The mutation sites of the four mutants F35Y, P40V, V45E and V45Y of cytochrome b5 are located at the edge of the heme-binding pocket. The solvent accessible areas of the “pocket inte-rior“ of the four mutants and the wild-type cytochrome b5 have been calculated based on their crystal structures at high resolu-tion. The change in the hydrophobicity of the heme-binding pocket resulting from the mutation can be quantitatively de-scribed using the difference of the solvent accessible area of the “pocket interior“ of each mutant from that of the wild-type cy-tochrome b5. The influences of the hydrophobicity of the heme-binding pocket on the protein stability and redox potential are discussed.     

Theoretical studies on the binding energy of β-sheet models

In this paper, B3LYP and MP2 methods are used to investigate the binding energy of seventeen antiparallel and parallel β-sheet models. The results indicate that the binding energy obtained from B3LYP calculations is weaker than that obtained from MP2 calculations but the relative binding energy yielded by B3LYP is almost the same as that by MP2. For the antiparallel β-sheets in which two N-H⋯O=C hydrogen bonds can form either a large hydrogen-bonded ring or a small hydrogen-bonded ring, the binding energy increases obviously when one large ring unit is added, whereas it only changes slightly when one small ring unit is added because of the secondary electrostatic repulsive interaction existing in the small ring unit which is estimated to be about 20 kJ/mol. For the parallel β-sheet models, the binding energy increases almost exactly linearly with the increase of the chain length.     

Neutral C–H bond vs. electron pair of N(sp~2): A binding site effect study of macrocycle anion receptor

To evaluate the effect of neutral C–H bond or electron pair of nitrogen atom with sp2hybridization(N(sp2)) involving into the same chemical environment for anion binding, two analogous tetracationic imidazolium macrocycles, namely cyclo[2](2,6-bis-(1H-imidazol-1-yl)pyridine) [2](1,3-dimethylenebenzene)(14+), and cyclo[2](2,6-bis-(1H-imidazol-1-yl)pyridine)[2](2,6-di methylenepyridine)(24+)were studied in detail as small inorganic anion receptors. The guest anions with different shapes are Cl,N3, NO3, and H2PO4. The host–guest interactions were characterized via1 H NMR spectroscopy,electrospray ionization mass spectrometry(ESI-MS) and single crystal X-ray crystallography. The results implied that macrocyclic hosts with similar backbone but two distinct binding sites(14+with neutral C–H vs. 24+with N(sp2)) vary markedly in their response to anions, including the binding modes and association constants. The finding will serve to the construction of new anion receptors, even improve insights into the anion binding process in biology.     

Computational prediction of structure-activity relationships for the binding of aminocyclitols to β-glucocerebrosidase

Glucocerebrosidase (GCase, acid β-Glucosidase) hydrolyzes the sphingolipid glucosylceramide into glucose and ceramide. Mutations in this enzyme lead to a lipid metabolism disorder known as Gaucher disease. The design of competitive inhibitors of GCase is a promising field of research for the design of pharmacological chaperones as new therapeutic agents. Using a series of recently reported molecules with experimental binding affinities for GCase in the nanomolar to micromolar range, we here report an extensive theoretical analysis of their binding mode. On the basis of molecular docking, molecular dynamics, and binding free energy calculations using the linear interaction energy method (LIE), we provide details on the molecular interactions supporting ligand binding in the different families of compounds. The applicability of other computational approaches, such as the COMBINE methodology, is also investigated. The results show the robustness of the standard parametrization of the LIE method, which reproduces the experimental affinities with a mean unsigned error of 0.7 kcal/mol. Several structure-activity relationships are established using the computational models here provided, including the identification of hot spot residues in the binding site. The models derived are envisaged as important tools in ligand-design programs for GCase inhibitors.     

A chemoenzymatic approach to the synthesis of the stereoisomers of a β-adrenergic receptor antagonist

The four stereoisomers of Δ²-isoxazoline 2, a β-adrenergic receptor antagonist structurally related to Falintolol 1, were prepared by an enzyme-catalyzed kinetic resolution of the unsaturated secondary alcohol (±)-7 followed by its cycloaddition to pyruvonitrile oxide. Through this strategy, diastereomeric aminoalcohols (+)-2a/(−)-2b and (−)-2a/(+)-2b were obtained in 99 and 92% enantiomeric excess, respectively. The absolute configuration to the target compounds was assigned via chemical correlation to the enantiomers of epoxides 4a and 4b, whose stereochemistry had been previously established.     

Prediction of ligand binding affinity using a multiple-conformations-multiple-protonation scheme: application to estrogen receptor α

A fast method that can predict the binding affinities of chemicals to a target protein with a high degree of accuracy will be very useful in drug design and regulatory science. We have been developing a scoring function for affinity prediction, which can be applied to extensive protein systems, and also trying to generate a prediction scheme that specializes in each target protein, with as high a predictive power as possible. In this study, we have constructed a prediction scheme with target-specific scores for estimating ligand-binding affinities to human estrogen receptor α (ERα), considering the major conformational change between agonist- and antagonist-bound forms and the change in protonation states of histidine at the ligand-binding site. The generated scheme calibrated with fewer training compounds (23 for the agonist-bound form, 17 for the antagonist-bound form) demonstrated good predictive power (a predictive r(2) of 0.83 for 154 validation compounds); this was also true for compounds with frameworks that were quite different from those of the training compounds. Our prediction scheme will be useful in drug development targeting ERα and in primary screening of endocrine disruptors, and provides a successful method of affinity prediction considering the major conformational changes in a protein.     

Simulation analysis of the cellulase Cel7A carbohydrate binding module on the surface of the cellulose Iβ

The Family 7 cellobiohydrolase (Cel7A) from Trichoderma reesei consists of a carbohydrate-binding module (CBM) joined by a linker to a catalytic domain. Cellulose hydrolysis is limited by the accessibility of Cel7A to crystalline substrates, which is perceived to be primarily mediated by the CBM. Here, the binding of CBM to the cellulose Iβ fiber is characterized by combined Brownian dynamics (BD) and molecular dynamics (MD) simulations. The results confirm that CBM prefers to dock to the hydrophobic than to the hydrophilic fiber faces. Both electrostatic (ES) and van der Waals (VDW) interactions are required for achieving the observed binding preference. The VDW interactions play a more important role in stabilizing the CBM-fiber binding, whereas the ES interactions contribute through the formation of a number of hydrogen bonds between the CBM and the fiber. At long distances, an ES steering effect is also observed that tends to align the CBM in an antiparallel manner relative to the fiber axis. Furthermore, the MD results reveal hindered diffusion of the CBM on all fiber surfaces. The binding of the CBM to the hydrophobic surfaces is found to involve partial dewetting at the CBM-fiber interface coupled with local structural arrangements of the protein. The present simulation results complement and rationalize a large body of previous work and provide detailed insights into the mechanism of the CBM-cellulose fiber interactions.