Aleksandra Maršavelski

The Quantum Nature of Drug-Receptor Interactions: Deuteration Changes Binding Affinities for Histamine Receptor Ligands

In this article we report a combined experimental and computational study concerning the effects of deuteration on the binding of histamine and two other histaminergic agonists to 3H-tiotidine-labeled histamine H2 receptor in neonatal rat astrocytes. Binding affinities were measured by displacing radiolabeled tiotidine from H2 receptor binding sites present on cultured neonatal rat astrocytes. Quantum-chemical calculations were performed by employing the empirical quantization of nuclear motion within a cluster model of the receptor binding site extracted from the homology model of the entire H2 receptor. Structure of H2 receptor built by homology modelling is attached in the supporting information (S1 Table) Experiments clearly demonstrate that deuteration affects the binding by increasing the affinity for histamine and reducing it for 2-methylhistamine, while basically leaving it unchanged for 4-methylhistamine. Ab initio quantum-chemical calculations on the cluster system extracted from the homology H2 model along with the implicit quantization of the acidic N–H and O–H bonds demonstrate that these changes in the binding can be rationalized by the altered strength of the hydrogen bonding upon deuteration known as the Ubbelohde effect. Our computational analysis also reveals a new mechanism of histamine binding, which underlines an important role of Tyr250 residue. The present work is, to our best knowledge, the first study of nuclear quantum effects on ligand receptor binding. The ligand H/D substitution is relevant for therapy in the context of perdeuterated and thus more stable drugs that are expected to enter therapeutic practice in the near future. Moreover, presented approach may contribute towards understanding receptor activation, while a distant goal remains in silico discrimination between agonists and antagonists based on the receptor structure.

Supramolecular Ionic-Liquid Gels with High Ionic Conductivity.

Supramolecular ionogels were prepared by the gelation of room-temperature ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIm][BF4 ]) with (S,S)-bis(leucinol)oxalamide. Remarkably, the ionic conductivity of solutions and ionogels with low gelator concentrations is higher than that of neat [BMIm][BF4 ]. On the basis of molecular dynamics simulations and quantum mechanical calculations, the origin of this phenomenon is attributed to the higher affinity of gelator molecules towards [BF4 ](-) ions, which reduces the electrostatic attraction between [BMIm](+) and [BF4 ](-) and thus increases their mobility. With increasing gelator concentration, the ionic conductivity decreases due to the formation of a denser gelator matrix, which hinders the pathways for ionic transport. However, even for very dense ionogels, this decrease is less than one order of magnitude relative to neat [BMIm][BF4 ], and thus they can be classified as highly conductive materials with strong potential for application as functional electrolytes.

15N NMR Spectroscopy, X-ray and Neutron Diffraction, Quantum-Chemical Calculations, and UV/vis-Spectrophotometric Titrations as Complementary Techniques for the Analysis of Pyridine-Supported Bicyclic Guanidine Superbases

Pyridine substituted with one and two bicyclic guanidine groups has been studied as a potential source of superbases. 2-{hpp}C5H4N (I) and 2,6-{hpp}2C5H3N (II) (hppH = 1,3,4,6,7,8-hexahydro-2H-pyrimido[1,2-a]pyrimidine) were protonated using [HNEt3][BPh4] to afford [I-H][BPh4] (1a), [II-H][BPh4] (2), and [II-H2][BPh4]2 (3). Solution-state (1)H and (15)N NMR spectroscopy shows a symmetrical cation in 2, indicating a facile proton-exchange process in solution. Solid-state (15)N NMR data differentiates between the two groups, indicating a mixed guanidine/guanidinium. X-ray diffraction data are consistent with protonation at the imine nitrogen, confirmed for 1a by single-crystal neutron diffraction. The crystal structure of 1a shows association of two [I-H](+) cations within a cage of [BPh4](-) anions. Computational analysis performed in the gas phase and in MeCN solution shows that the free energy barrier to transfer a proton between imino centers in [II-H](+) is 1 order of magnitude lower in MeCN than in the gas phase. The results provide evidence that linking hpp groups with the pyridyl group stabilizes the protonation center, thereby increasing the intrinsic basicity in the gas phase, while the bulk prevents efficient cation solvation, resulting in diminished pKa(MeCN) values. Spectrophotometrically measured pKa values are in excellent agreement with calculated values and confirm that I and II are superbases in solution.

What a Difference a Methyl Group Makes: The Selectivity of Monoamine Oxidase B Towards Histamine and N-Methylhistamine

Monoamine oxidase (MAO) enzymes catalyze the degradation of a very broad range of biogenic and dietary amines including many neurotransmitters in the brain, whose imbalance is extensively linked with the biochemical pathology of various neurological disorders. Although sharing around 70 % sequence identity, both MAO A and B isoforms differ in substrate affinities and inhibitor sensitivities. Inhibitors that act on MAO A are used to treat depression, due to their ability to raise serotonin concentrations, whereas MAO B inhibitors decrease dopamine degradation and improve motor control in patients with Parkinson disease. Despite this functional importance, the factors affecting MAO selectivity are poorly understood. Here, we used a combination of molecular dynamics (MD) simulations, molecular mechanics with Poisson-Boltzmann and surface area solvation (MM-PBSA) binding free energy evaluations, and quantum mechanical (QM) cluster calculations to address the unexpected, yet challenging MAO B selectivity for N-methylhistamine (NMH) over histamine (HIS), differing only in a single methyl group distant from the reactive ethylamino center. This study shows that a dominant selectivity contribution is offered by a lower activation free energy for NMH by 2.6 kcal mol-1 , in excellent agreement with the experimental ΔΔG≠EXP =1.4 kcal mol-1 , together with a more favorable reaction exergonicity and active-site binding. This study also confirms the hydrophobic nature of the MAO B active site and underlines the important role of Ile199, Leu171, and Leu328 in properly orienting substrates for the reaction.

Isolation, biochemical characterization and anti-bacterial activity of BPIFA2 protein

OBJECTIVE Human BPIFA2 (parotid secretory protein) is a ubiquitous soluble salivary protein, which belongs to the PLUNC family of proteins. Having sequence similarity to bactericidal/permeability-increasing protein and lipopolysaccharide-binding protein, PLUNC proteins are probably involved in local antibacterial response at mucosal sites, such as oral cavity. The aim of the study was to isolate and characterize human BPIFA2. DESIGN In this paper, we report one-step affinity chromatography method for BPIFA2 purification from whole human saliva. The isolated BPIFA2 was identified by trypsin mass fingerprinting and characterized by electrophoretic methods. Antibacterial activity of BPIFA2 against model microorganism Pseudomonas aeruginosa was shown in minimum inhibitory concentration and time kill study assays. RESULTS The protein showed microheterogeneity, both in molecular weight and pI value. BPIFA2 inhibited the growth of P. aeruginosa in microgram concentration range determined by minimum inhibitory concentration assay. In the time kill study, 32μg/mL BPIFA2 showed clear bactericidal activity and did not cause any aggregation of bacteria. CONCLUSION Affinity chromatography is well suited for isolation of functional BPIFA2 with a potent bactericidal activity against P. aeruginosa.

Catalytic Dyad in the SGNH Hydrolase Superfamily: In-depth Insight into Structural Parameters Tuning the Catalytic Process of Extracellular Lipase from Streptomyces rimosus

SrLip is an extracellular enzyme from Streptomyces rimosus (Q93MW7) exhibiting lipase, phospholipase, esterase, thioesterase, and tweenase activities. The structure of SrLip is one of a very few lipases, among the 3D-structures of the SGNH superfamily of hydrolases, structurally characterized by synchrotron diffraction data at 1.75 Å resolution (PDB: 5MAL ). Its crystal structure was determined by molecular replacement using a homology model based on the crystal structure of phospholipase A1 from Streptomyces albidoflavus (PDB: 4HYQ ). The structure reveals the Rossmann-like 3-layer αβα sandwich fold typical of the SGNH superfamily stabilized by three disulfide bonds. The active site shows a catalytic dyad involving Ser10 and His216 with Ser10-OγH···NεHis216, His216-NδH···O═C-Ser214, and Gly54-NH···Oγ-Ser10 hydrogen bonds essential for the catalysis; the carbonyl oxygen of the Ser214 main chain acts as a hydrogen bond acceptor ensuring the orientation of the His216 imidazole ring suitable for a proton transfer. Molecular dynamics simulations of the apoenzyme and its complex with p-nitrophenyl caprylate were used to probe the positioning of the substrate ester group within the active site and its aliphatic chain within the binding site. Quantum-mechanical calculations at the DFT level revealed the precise molecular mechanism of the SrLip catalytic activity, demonstrating that the overall hydrolysis is a two-step process with acylation as the rate-limiting step associated with the activation free energy of ΔG⧧ENZ = 17.9 kcal mol-1, being in reasonable agreement with the experimental value of 14.5 kcal mol-1, thus providing strong support in favor of the proposed catalytic mechanism based on a dyad.

Empirical Valence Bond Simulations Suggest a Direct Hydride Transfer Mechanism for Human Diamine Oxidase

Diamine oxidase (DAO) is an enzyme involved in the regulation of cell proliferation and the immune response. This enzyme performs oxidative deamination in the catabolism of biogenic amines, including, among others, histamine, putrescine, spermidine, and spermine. The mechanistic details underlying the reductive half-reaction of the DAO-catalyzed oxidative deamination which leads to the reduced enzyme cofactor and the aldehyde product are, however, still under debate. The catalytic mechanism was proposed to involve a prototropic shift from the substrate–Schiff base to the product–Schiff base, which includes the rate-limiting cleavage of the Cα–H bond by the conserved catalytic aspartate. Our detailed mechanistic study, performed using a combined quantum chemical cluster approach with empirical valence bond simulations, suggests that the rate-limiting cleavage of the Cα–H bond involves direct hydride transfer to the topaquinone cofactor—a mechanism that does not involve the formation of a Schiff base. Additional investigation of the D373E and D373N variants supported the hypothesis that the conserved catalytic aspartate is indeed essential for the reaction; however, it does not appear to serve as the catalytic base, as previously suggested. Rather, the electrostatic contributions of the most significant residues (including D373), together with the proximity of the Cu2+ cation to the reaction site, lower the activation barrier to drive the chemical reaction.

A novel antimicrobial target—expanded and revisited mode of action of pantothenamides

Pantothenamides are analogs of pantothenic acid (vitamin B5), which is a natural precursor of coenzyme A (CoA). It has been shown that these compounds, predominantly N-substituted pantothenamides, possess antimicrobial activity against various pathogenic bacteria such as E. coli and S. aureus. It is widely accepted that these compounds act through combined inhibition of coenzyme A and fatty acid synthesis. However, the precise mechanism of action remains unrevealed. Here is reported the identification of a novel target of pantothenamides, never considered before. Molecular dynamics simulations together with free energy calculations reveal that the hydrophobic pocket of the acyl carrier protein (ACP) binds N-pentylpantothenamide. Consequently, the sequestration of the acyl chain attached to the Ppant prosthetic arm is defunct since the inhibitor occupies the hydrophobic core of the ACP. Thus, the acyl chain remains solvent-exposed and susceptible to hydrolysis. Moreover, the ACP with N-pentylpantothenamide bound could change its chain-flipping ability as well as its interaction propensity towards downstream enzyme partners of the fatty acid synthesis pathway, which could result in the suppression of the fatty acid synthesis rate.

Insight into recognition between aminoacyl carrier protein and its binding partner

Acyl carrier-protein (ACP) is one of the most promiscuous protein in terms of protein-protein interactions within the cell. Common to both, primary and secondary cell metabolic pathways, ACP acts as a central acceptor of various intermediate products of fatty acid, polyketide and nonribosomal peptide synthesis pathways. Thus, its quite puzzling how ACP selects correct protein partner among many possible upstream and downstream binding partners. To address this question we choose recently described protein-protein complex formed between aminoacyl carrier protein from Bradyrhizobium japonicum (Bj CP) and Bj Gly:CP ligase 1. We performed molecular dynamics simulations, MM-GBSA binding free energy calculations, alanine scanning mutagenesis, and multiple sequence alignment to get deeper insight into relevant protein-protein interactions between three forms of the aminoacyl carrier protein and their partner protein Bj Gly:CP ligase 1. Our result are in line with experimental findings and indicate that correct protein-protein communication is regulated through specific non-polar contribution of the interacting amino acid residues placed at the interface between aminoacyl carrier protein and partner Bj Gly:CP ligase 1. In silico alanine scanning mutagenesis along with multiple sequence alignment aided in classifying interacting residues as universally required for aminoacyl carrier protein binding and those that determine the selectivity.