Ulrika Brath

Symmetric Halogen Bonding Is Preferred in Solution

Halogen bonding is a recently rediscovered secondary interaction that shows potential to become a complementary molecular tool to hydrogen bonding in rational drug design and in material sciences. Whereas hydrogen bond symmetry has been the subject of systematic studies for decades, the understanding of the analogous three-center halogen bonds is yet in its infancy. The isotopic perturbation of equilibrium (IPE) technique with (13)C NMR detection was applied to regioselectively deuterated pyridine complexes to investigate the symmetry of [N-I-N](+) and [N-Br-N](+) halogen bonding in solution. Preference for a symmetric arrangement was observed for both a freely adjustable and for a conformationally restricted [N-X-N](+) model system, as also confirmed by computation on the DFT level. A closely attached counterion is shown to be compatible with the preferred symmetric arrangement. The experimental observations and computational predictions reveal a high energetic gain upon formation of symmetric, three-center four-electron halogen bonding. Whereas hydrogen bonds are generally asymmetric in solution and symmetric in the crystalline state, the analogous bromine and iodine centered halogen bonds prefer symmetric arrangement in solution.

Differential Responses of the Backbone and Side-Chain Conformational Dynamics in FKBP12 upon Binding the Transition-State Analog FK506: Implications for Transition-State Stabilization and Target Protein Recognition

FKBP12 serves a dual role as a peptidyl-prolyl cis-trans isomerase and as a modulator of several cell signaling pathways. The macrolide FK506 is a transition-state analog of the catalyzed reaction and displaces FKBP12 from its natural target proteins. We compared the conformational exchange dynamics of the backbone and methyl-bearing side chains of FKBP12 in the free and FK506-bound states using NMR relaxation-dispersion experiments. Our results show that the free enzyme exchanges between the ground state and an excited state that resembles the ligand-bound state or Michaelis complex. In FK506-bound FKBP12, the backbone is confined to a single conformation, while conformational exchange prevails for many methyl groups. The residual side-chain dynamics in the transition-state analog-bound state suggests that the transition-state ensemble involves multiple conformations, a finding that challenges the long-standing concept of conformational restriction in the transition-state complex. Furthermore, exchange between alternative conformations is observed in the bound state for an extended network of methyl groups that includes locations remote from the active site. Several of these locations are known to be important for interactions with cellular target proteins, including calcineurin and the ryanodine receptor, suggesting that the conformational heterogeneity might play a role in the promiscuous binding of FKBP12 to different targets.

Solvent effects on halogen bond symmetry

The symmetric arrangement of the iodine and bromine centred 3-center–4-electron halogen bond is revealed to remain preferred in a polar, aprotic solvent environment. Acetonitrile is unable to compete with pyridine for halogen bonding; however, its polarity weakly modulates the energy of the interaction and influences IPE-NMR experiments.

Paramagnetic Ligand Tagging To Identify Protein Binding Sites

Transient biomolecular interactions are the cornerstones of the cellular machinery. The identification of the binding sites for low affinity molecular encounters is essential for the development of high affinity pharmaceuticals from weakly binding leads but is hindered by the lack of robust methodologies for characterization of weakly binding complexes. We introduce a paramagnetic ligand tagging approach that enables localization of low affinity protein–ligand binding clefts by detection and analysis of intermolecular protein NMR pseudocontact shifts, which are invoked by the covalent attachment of a paramagnetic lanthanoid chelating tag to the ligand of interest. The methodology is corroborated by identification of the low millimolar volatile anesthetic interaction site of the calcium sensor protein calmodulin. It presents an efficient route to binding site localization for low affinity complexes and is applicable to rapid screening of protein–ligand systems with varying binding affinity.

Mapping the sevoflurane-binding sites of calmodulin.

General anesthetics, with sevoflurane (SF) being the first choice inhalational anesthetic agent, provide reversible, broad depressor effects on the nervous system yet have a narrow margin of safety. As characterization of low‐affinity binding interactions of volatile substances is exceptionally challenging with the existing methods, none of the numerous cellular targets proposed as chief protagonists in anesthesia could yet be confirmed. The recognition that most critical functions modulated by volatile anesthetics are under the control of intracellular Ca2+ concentration, which in turn is primarily regulated by calmodulin (CaM), motivated us for characterization of the SF–CaM interaction. Solution NMR (Nuclear Magnetic Resonance) spectroscopy was used to identify SF‐binding sites using chemical shift displacement, NOESY and heteronuclear Overhauser enhancement spectroscopy (HOESY) experiments. Binding affinities were measured using ITC (isothermal titration calorimetry). SF binds to both lobes of (Ca2+)4‐CaM with low mmol/L affinity whereas no interaction was observed in the absence of Ca2+. SF does not affect the calcium binding of CaM. The structurally closely related SF and isoflurane are shown to bind to the same clefts. The SF‐binding clefts overlap with the binding sites of physiologically relevant ion channels and bioactive small molecules, but the binding affinity suggests it could only interfere with very weak CaM targets.

The impact of interchain hydrogen bonding on β‐hairpin stability is readily predicted by molecular dynamics simulation

Peptides are frequently used model systems for protein folding. They are also gaining increased importance as therapeutics. Here, the ability of molecular dynamics (MD) simulation for describing the structure and dynamics of β‐hairpin peptides was investigated, with special attention given to the impact of a single interstrand sidechain to sidechain interaction. The MD trajectories were compared to structural information gained from solution NMR. By assigning frames from restraint‐free MD simulations to an intuitive hydrogen bond on/off pattern, folding ratios and folding pathways were predicted. The computed molecular model successfully reproduces the folding ratios determined by NMR, indicating that MD simulation may be straightforwardly used as a screening tool in β‐hairpin design.

Assessing the Ability of Spectroscopic Methods to Determine the Difference in the Folding Propensities of Highly Similar β-Hairpins

We have evaluated the ability of nuclear magnetic resonance (NMR) and circular dichroism (CD) spectroscopies to describe the difference in the folding propensities of two structurally highly similar cyclic β-hairpins, comparing the outcome to that of molecular dynamics simulations. NAMFIS-type NMR ensemble analysis and CD spectroscopy were observed to accurately describe the consequence of altering a single interaction site, whereas a single-site 13C NMR chemical shift melting curve-based technique was not.