Extending Halogen-based Medicinal Chemistry to Proteins: IODO-INSULIN AS A CASE STUDY

Hydrogen- and halogen bond synergy in the self-assembly of 3,5-dihalo-tyrosines: structural and theoretical insights

Halogenation, generally introduced on aromatic aminoacids, is becoming a key supramolecular tool in peptides. Herein, we report the crystal structures and DFT study of two bis-halogenated tyrosines showing the subtle...

Hydration dynamics and IR spectroscopy of 4-fluorophenol

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy. In this work, the structural dynamics and infrared spectroscopy of para-fluorophenol (F-PhOH) and phenol (PhOH) is investigated in the gas phase and in water using a combination of experiment and molecular dynamics (MD) simulations. The gas phase and solvent dynamics around F-PhOH and PhOH is characterized from atomistic simulations using empirical energy functions with point charges or multipoles for the electrostatics, Machine Learning (ML) based parametrizations and with full ab initio (QM) and mixed Quantum Mechanical/Molecular Mechanics (QM/MM) simulations with a particular focus on the CF- and OH-stretch region. The CF-stretch band is heavily mixed with other modes whereas the OH-stretch in solution displays a characteristic high-frequency peak around 3600 cm⁻¹ most likely associated with the –OH group of PhOH and F-PhOH together with a characteristic progression below 3000 cm⁻¹ due to coupling with water modes which is also reproduced by several of the simulations. Solvent and radial distribution functions indicate that the CF-site is largely hydrophobic except for simulations using point charges which renders them unsuited for correctly describing hydration and dynamics around fluorinated sites. The hydrophobic character of the CF-group is particularly relevant for applications in pharmaceutical chemistry with a focus on local hydration and interaction with the surrounding protein.

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy.

Structural Stability of Insulin Oligomers and Protein Association-Dissociation Processes: Free Energy Landscape and Universal Role of Water

Association and dissociation of proteins are important biochemical events. In this Feature Article, we analyze the available studies of these processes for insulin oligomers in aqueous solution. We focus on the solvation of the insulin monomer in water, stability and dissociation of its dimer, and structural integrity of the hexamer. The intricate role of water in solvation of the dimer- and hexamer-forming surfaces, in long-range interactions between the monomers and the stability of the oligomers, is discussed. Ten water molecules inside the central cavity stabilize the structure of the insulin hexamer. We discuss how different order parameters can be used to understand the dissociation of the insulin dimer. The calculation of the rate using a recently computed multidimensional free energy provides considerable insight into the interplay between protein and water dynamics.

A Reduced Generalized Force Field for Biological Halogen Bonds

The halogen bond (or X-bond) is a noncovalent interaction that is increasingly recognized as an important design tool for engineering protein-ligand interactions and controlling the structures of proteins and nucleic acids. In the past decade, there have been significant efforts to characterize the structure-energy relationships of this interaction in macromolecules. Progress in the computational modeling of X-bonds in biological molecules, however, has lagged behind these experimental studies, with most molecular mechanics/dynamics-based simulation methods not properly treating the properties of the X-bond. We had previously derived a force field for biological X-bonds (ffBXB) based on a set of potential energy functions that describe the anisotropic electrostatic and shape properties of halogens participating in X-bonds. Although fairly accurate for reproducing the energies within biomolecular systems, including X-bonds engineered into a DNA junction, the ffBXB with its seven variable parameters was considered to be too unwieldy for general applications. In the current study, we have generalized the ffBXB by reducing the number of variables to just one for each halogen type and show that this remaining electrostatic variable can be estimated for any new halogenated molecule through a standard restricted electrostatic potential calculation of atomic charges. In addition, we have generalized the ffBXB for both nucleic acids and proteins. As a proof of principle, we have parameterized this reduced and more general ffBXB against the AMBER force field. The resulting parameter set was shown to accurately recapitulate the quantum mechanical landscape and experimental interaction energies of X-bonds incorporated into DNA junction and T4 lysozyme model systems. Thus, this reduced and generalized ffBXB is more readily adaptable for incorporation into classical molecular mechanics/dynamics algorithms, including those commonly used to design inhibitors against therapeutic targets in medicinal chemistry and materials in biomolecular engineering.

Insights into the Sensitivity of Arginine Concentration to Preserve the Folded Form of Insulin Monomer under Thermal Stress

Arginine, although popularly known as aggregation suppressor additive, has been found to quench proteins' structure and function by destabilizing their conformations. Driven by such controversial evidence, in this work we performed a series of atomistic molecular dynamics simulations of insulin monomer, a biologically active hormone protein, in arginine solution of varying concentrations (0.5, 1, and 2 M) at ambient and elevated temperature (400 K) to explore the arginine concentration driven structure-based stability of the protein. Our study reveals that the flexibility of the protein's structure is dependent on the arginine concentration, and among all the used solutions, 2 M arginine, a "neutral crowder" that mimics the cellular environment, can preserve the native folded form of the protein at ambient temperature in an excellent manner. Further, while the protein unfolds at 400 K in pure water, this solution worked satisfactorily to preserve the protein's folded conformation more firmly than the other solutions. The replica-exchange MD of insulin in 2 M arginine solution further supports the fact. In this aspect an important issue in molecular pharmacology is to identify and recognize the physical origin of the stability of a protein, i.e, in this case, how arginine directs the conformational flexibility of the protein and preserves its native folded form. We identified that the exclusion of arginine from the protein surface increases the local structuration of water around the protein, thereby preserving its "biological water" layer, and makes the protein more hydrated at 2 M concentration as compared to the other arginine solutions. Additionally, our microscopic investigation on the interactions of the protein-solvation layer revealed that the structural heterogeneity of the protein surface, arising from the differential physicochemical nature of the amino acid residues, controls the favorable formation of sluggish water-arginine mixed solvation layer at higher arginine concentration that helps the protein to maintain its structural rigidity. Importantly, apart from the protein-solvent hydrogen-bonding interactions, the anion-pi interactions, established between the carboxyl group of arginine and the aromatic amino acid residues of insulin, were recognized to facilitate the protein to maintain its native folded form at the experimental temperatures.

Substitution of an Internal Disulfide Bridge with a Diselenide Enhances both Foldability and Stability of Human Insulin

Insulin analogues, mainstays in the modern treatment of diabetes mellitus, exemplify the utility of protein engineering in molecular pharmacology. Whereas chemical syntheses of the individual A and B chains were accomplished in the early 1960s, their combination to form native insulin remains inefficient because of competing disulfide pairing and aggregation. To overcome these limitations, we envisioned an alternative approach: pairwise substitution of cysteine residues with selenocysteine (Sec, U). To this end, CysA6 and CysA11 (which form the internal intrachain A6-A11 disulfide bridge) were each replaced with Sec. The A chain[C6U, C11U] variant was prepared by solid-phase peptide synthesis; while sulfitolysis of biosynthetic human insulin provided wild-type B chain-di-S-sulfonate. The presence of selenium atoms at these sites markedly enhanced the rate and fidelity of chain combination, thus solving a long-standing challenge in chemical insulin synthesis. The affinity of the Se-insulin analogue for the lectin-purified insulin receptor was indistinguishable from that of WT-insulin. Remarkably, the thermodynamic stability of the analogue at 25 °C, as inferred from guanidine denaturation studies, was augmented (ΔΔGu ≈0.8 kcal mol-1 ). In accordance with such enhanced stability, reductive unfolding of the Se-insulin analogue and resistance to enzymatic cleavage by Glu-C protease occurred four times more slowly than that of WT-insulin. 2D-NMR and X-ray crystallographic studies demonstrated a native-like three-dimensional structure in which the diselenide bridge was accommodated in the hydrophobic core without steric clash.

Halogen bonding in halocarbon-protein complexes and computational tools for rational drug design

Introduction: Halogens have a prominent role in drug design. Often used as a mean to improve ADME properties, they are also becoming a tool in protein-ligand recognition given their ability to form a non-covalent interaction, termed halogen bond, where halogens act as electrophilic species interacting with electron-rich partners. Rational drug design of halogen-bonding lead molecules requires an accurate description of halocarbon-protein complexes by computational tools though not all methods are able to tackle this non-covalent interaction. Areas covered: The authors present a review of computational methodologies that can be used to properly describe halogen bonds in the context of protein-ligand complexes, providing also insights on how these methods can be used in the context of computer-aided drug design. Expert opinion: Although in the last few years many computational tools, ranging from fast screening methods to the more expensive QM calculations, have been developed to tackle the halogen bonding phenomenon, they are not yet standard in the literature. This will eventually change as official software distributions are including support for halogen bonding in their methods. Tackling desolvation of halogenated species seems to be a good strategy to improve the accuracy of computational methods, that will be more commonly used prior to laboratory work in the future.

Novel protein science enabled by total chemical synthesis

Chemical synthesis is a well-established method for the preparation in the research laboratory of multiple-tens-of-milligram amounts of correctly folded, high purity protein molecules. Chemically synthesized proteins enable a broad spectrum of novel protein science. Racemic mixtures consisting of d-protein and l-protein enantiomers facilitate crystallization and determination of protein structures by X-ray diffraction. d-Proteins enable the systematic development of unnatural mirror image protein molecules that bind with high affinity to natural protein targets. The d-protein form of a therapeutic target can also be used to screen natural product libraries to identify novel small molecule leads for drug development. Proteins with novel polypeptide chain topologies including branched, circular, linear-loop, and interpenetrating polypeptide chains can be constructed by chemical synthesis. Medicinal chemistry can be applied to optimize the properties of therapeutic protein molecules. Chemical synthesis has been used to redesign glycoproteins and for the a priori design and construction of covalently constrained novel protein scaffolds not found in nature. Versatile and precise labeling of protein molecules by chemical synthesis facilitates effective application of advanced physical methods including multidimensional nuclear magnetic resonance and time-resolved FTIR for the elucidation of protein structure-activity relationships. The chemistries used for total synthesis of proteins have been adapted to making artificial molecular devices and protein-inspired nanomolecular constructs. Research to develop mirror image life in the laboratory is in its very earliest stages, based on the total chemical synthesis of d-protein forms of polymerase enzymes.

Enzyme kinetics from circular dichroism of insulin reveals mechanistic insights into the regulation of insulin-degrading enzyme

Insulin-degrading enzyme (IDE) is a zinc metalloprotease that selectively degrades biologically important substrates associated with type 2 diabetes and Alzheimer’s disease (AD). As such, IDE is an attractive target for therapeutic innovations. A major requirement is an understanding of how other molecules present in cells regulate the activity of the enzyme toward insulin, IDE’s most important physiologically relevant substrate. Previous kinetic studies of the IDE-dependent degradation of insulin in the presence of potential regulators have used iodinated insulin, a chemical modification that has been shown to alter the biological and biochemical properties of insulin. Here, we present a novel kinetic assay that takes advantage of the loss of helical circular dichroic signals of insulin with IDE-dependent degradation. As proof of concept, the resulting Michaelis–Menten kinetic constants accurately predict the known regulation of IDE by adenosine triphosphate (ATP). Intriguingly, we found that when Mg²⁺ is present with ATP, the regulation is abolished. The implication of this result for the development of preventative and therapeutic strategies for AD is discussed. We anticipate that the new assay presented here will lead to the identification of other small molecules that regulate the activity of IDE toward insulin.

Biomolecular Simulations of Halogen Bonds with a GROMOS Force Field

Halogen bonds (XBs) are non-covalent interactions in which halogens (X), acting as electrophiles, interact with Lewis bases. XBs are able to mediate protein-ligand recognition and therefore play an important role in rational drug design. In this context, the development of molecular modeling tools that can tackle XBs is paramount. XBs are predominantly explained by the existence of a positive region on the electrostatic potential of X named the σ-hole. Typically, with molecular mechanics force fields, this region is modeled using a charged extra point (EP) linked to X along the R-X covalent bond axis. In this work, we developed the first EP-based strategy for GROMOS force fields (specifically GROMOS 54A7) using bacteriophage T4 lysozyme in complex with both iodobenzene and iodopentafluorobenzene as a prototype system. Several EP parametrization schemes were tested by adding a virtual interaction site to ligand topologies retrieved from the Automated Topology Builder (ATB) and Repository. Contrary to previous approaches using other force fields, our analysis is based on the capability of each parametrization scheme to sample XBs during MD simulations. Our results indicate that the implementation of an EP at a distance from iodine corresponding to R_min provides a good qualitative description of XBs in MD simulations, supporting the compatibility of our approach with the GROMOS 54A7 force field.

Structure-based stabilization of insulin as a therapeutic protein assembly via enhanced aromatic-aromatic interactions

Key contributions to protein structure and stability are provided by weakly polar interactions, which arise from asymmetric electronic distributions within amino acids and peptide bonds. Of particular interest are aromatic side chains whose directional π-systems commonly stabilize protein interiors and interfaces. Here, we consider aromatic-aromatic interactions within a model protein assembly: the dimer interface of insulin. Semi-classical simulations of aromatic-aromatic interactions at this interface suggested that substitution of residue TyrB26 by Trp would preserve native structure while enhancing dimerization (and hence hexamer stability). The crystal structure of a [TrpB26]insulin analog (determined as a T3Rf3 zinc hexamer at a resolution of 2.25 Å) was observed to be essentially identical to that of WT insulin. Remarkably and yet in general accordance with theoretical expectations, spectroscopic studies demonstrated a 150-fold increase in the in vitro lifetime of the variant hexamer, a critical pharmacokinetic parameter influencing design of long-acting formulations. Functional studies in diabetic rats indeed revealed prolonged action following subcutaneous injection. The potency of the TrpB26-modified analog was equal to or greater than an unmodified control. Thus, exploiting a general quantum-chemical feature of protein structure and stability, our results exemplify a mechanism-based approach to the optimization of a therapeutic protein assembly.

Increasing Enzyme Stability and Activity through Hydrogen Bond-Enhanced Halogen Bonds

The construction of more stable proteins is important in biomolecular engineering, particularly in the design of biologics-based therapeutics. We show here that replacing the tyrosine at position 18 (Y18) of T4 lysozyme with the unnatural amino acid m-chlorotyrosine (^mClY) increases both the thermal stability (increasing the melting temperature by ∼1 °C and the melting enthalpy by 3 kcal/mol) and the enzymatic activity at elevated temperatures (15% higher than that of the parent enzyme at 40 °C) of this classic enzyme. The chlorine of ^mClY forms a halogen bond (XB) to the carbonyl oxygen of the peptide bond at glycine 28 (G28) in a tight loop near the active site. In this case, the XB potential of the typically weak XB donor Cl is shown from quantum chemical calculations to be significantly enhanced by polarization via an intramolecular hydrogen bond (HB) from the adjacent hydroxyl substituent of the tyrosyl side chain, resulting in a distinctive synergistic HB-enhanced XB (or HeX-B for short) interaction. The larger halogens (bromine and iodine) are not well accommodated within this same loop and, consequently, do not exhibit the effects on protein stability or function associated with the HeX-B interaction. Thus, we have for the first time demonstrated that an XB can be engineered to stabilize and increase the activity of an enzyme, with the increased stabilizing potential of the HeX-B further extending the application of halogenated amino acids in the design of more stable protein therapeutics.

From in silica to in silico: retention thermodynamics at solid–liquid interfaces

The dynamics of solvated molecules at the solid/liquid interface is essential for a molecular-level understanding for the solution thermodynamics in reversed phase liquid chromatography (RPLC).

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Molecular Mechanisms Underlying Solute Retention at Heterogeneous Interfaces

Despite considerable effort, a molecular-level understanding of the mechanisms governing adsorption/desorption in reversed-phase liquid chromatography is still lacking. This impedes rational design of columns and the development of reliable, computationally more efficient approaches to predict the selectivity of a particular column design. Using state-of-the art, validated force fields and free-energy simulations, the adsorption thermodynamics of benzene derivatives is investigated in atomistic detail and provides a quantitative microscopic understanding of retention when compared with experimental data. It is found that pure partitioning or pure adsorption is rather the exception than the rule. Typically, a pronounced ∼1 kcal/mol stabilization on the surface is accompanied by a broad trough indicative of partitioning before the probe molecule incorporates into the mobile phase. The present findings provide a quantitative and rational basis to develop improved effective, coarse-grained computational models and to design columns for specific applications.