Biophysical optimization of a therapeutic protein by nonstandard mutagenesis: studies of an iodo-insulin derivative

Synthesis of a highly thermostable insulin by phenylalanine conjugation at B29 Lysine

Globally, millions of diabetic patients require daily life-saving insulin injections. Insulin heat-lability and fibrillation pose significant challenges, especially in parts of the world without ready access to uninterrupted refrigeration. Here, we have synthesized four human insulin analogs by conjugating ε-amine of B29 lysine of insulin with acetic acid, phenylacetic acid, alanine, and phenylalanine residues. Of these, phenylalanine-conjugated insulin, termed FHI, was the most stable under high temperature (65 °C), elevated salt stress (25 mM NaCl), and varying pH levels (ranging from highly acidic pH 1.6 to physiological pH 7.4). It resists fibrillation for a significantly longer duration with sustained biological activity in in vitro, ex vivo, and in vivo and displays prolonged stability over its native counterpart. We further unravel the critical interactions, such as additional aromatic π-π interactions and hydrogen bonding in FHI, that are notably absent in native insulin. These interactions confer enhanced structural stability of FHI and offer a promising solution to the challenges associated with insulin heat sensitivity.

4S-fluorination of ProB29 in insulin lispro slows fibril formation

Recombinant insulin is a life-saving therapeutic for millions of patients affected by diabetes mellitus. Standard mutagenesis has led to insulin variants with improved control of blood glucose; for instance, the fast-acting insulin lispro contains two point mutations that suppress dimer formation and expedite absorption. However, insulins undergo irreversible denaturation, a process accelerated for the insulin monomer. Here we replace ProB29 of insulin lispro with 4R-fluoroproline, 4S-fluoroproline, and 4,4-difluoroproline. All three fluorinated lispro variants reduce blood glucose in diabetic mice, exhibit similar secondary structure as measured by CD, and rapidly dissociate from the zinc- and resorcinol-bound hexamer upon dilution. Notably, however, we find that 4S-fluorination of ProB29 delays the formation of undesired insulin fibrils that can accumulate at the injection site in vivo and can complicate insulin production and storage. These results demonstrate how subtle molecular changes achieved through non-canonical amino acid mutagenesis can improve the stability of protein therapeutics.

Omniligase-1-Mediated Phage-Peptide Library Modification and Insulin Engineering

Chemical and enzymatic modifications of peptide-displayed libraries have been successfully employed to expand the phage display library. However, the requirement of specific epitopes and scaffolds has limited the scope of protein engineering using phage display. In this study, we present a novel approach utilizing omniligase-1-mediated selective and specific ligation on the phage pIII protein, offering a high conversion rate and compatibility with commercially available phage libraries. We applied this method to perform high-throughput engineering of insulin analogues with randomized B chain C-terminal regions. Insulin analogues with different B chain C-terminal segments were selected and exhibited biological activity equivalent to that of human insulin. Molecular dynamics studies of insulin analogues revealed a novel interaction between the insulin B27 residue and insulin receptor L1 domain. In summary, our findings highlight the potential of omniligase-1-mediated phage display in the development and screening of disulfide-rich peptides and proteins. This approach holds promise for the creation of novel insulin analogues with enhanced therapeutic properties and exhibits potential for the development of other therapeutic compounds.

Serine-mediated hydrazone ligation displaying insulin-like peptides on M13 phage pIII

Phage display has emerged as a tool for the discovery of therapeutic antibodies and proteins. However, the effective display and engineering of structurally complex proteins, such as insulin, pose significant challenges due to the sequence of insulin, which is composed of two peptide chains linked by three disulfide bonds. In this study, we developed a new approach for the display of insulin-like peptides on M13 phage pIII, employing N-terminal serine-mediated hydrazone ligation. The insulin-displaying phage retains the biological binding affinity of human insulin. To address the viability loss after ligation, we introduced a trypsin-cleavable spacer on pIII, enabling insulin-displayed phage library selection. This method offers a general pathway for the display of structurally complex proteins on pIII, enhancing the practicality of selecting chemically modified phage libraries and opening avenues for the engineering of new insulin analogs for the treatment of diabetes by using phage display.

Strategies for interference of insulin fibrillogenesis: challenges and advances

The discovery of insulin came up with very high hopes for diabetic patients. In the year 2021, the world celebrated the 100 th anniversary of the discovery of this vital hormone. However, external use of insulin is highly affected by its aggregating tendency that occurs during its manufacturing, transportation, and improper handling which ultimately leads its pharmaceutically and biologically ineffective form. In this review, we aim to discuss the various approaches used for decelerating insulin aggregation which results in the enhancement of its overall structural stability and usage. The approaches that are discussed are broadly classified as either a measure through excipient additions or by intrinsic modifications in the insulin native structure.

Altered islet prohormone processing: A cause or consequence of diabetes?

Peptide hormones are first produced as larger precursor prohormones that require endoproteolytic cleavage to liberate the mature hormones. A structurally conserved but functionally distinct family of nine prohormone convertase enzymes (PCs) are responsible for cleavage of protein precursors of which PC1/3 and PC2 are known to be exclusive to neuroendocrine cells and responsible for prohormone cleavage. Differential expression of PCs within tissues define prohormone processing; whereas glucagon is the major product liberated from proglucagon via PC2 in pancreatic α-cells, proglucagon is preferentially processed by PC1/3 in intestinal L cells to produce glucagon-like peptides 1 and 2 (GLP-1, GLP-2). Beyond our understanding of processing of islet prohormones in healthy islets, there is convincing evidence that proinsulin, proIAPP, and proglucagon processing is altered during prediabetes and diabetes. There is predictive value of elevated circulating proinsulin or proinsulin : C-peptide ratio for progression to type 2 diabetes and elevated proinsulin or proinsulin : C-peptide is predictive for development of type 1 diabetes in at risk groups. After onset of diabetes, patients have elevated circulating proinsulin and proIAPP and proinsulin may be an autoantigen in type 1 diabetes. Further, preclinical studies reveal that α-cells have altered proglucagon processing during diabetes leading to increased GLP-1 production. We conclude that despite strong associative data, current evidence is inconclusive on the potential causal role of impaired prohormone processing in diabetes, and suggest that future work should focus on resolving the question of whether altered prohormone processing is a causal driver or merely a consequence of diabetes pathology.

Evolution of insulin at the edge of foldability and its medical implications

Proteins have evolved to be foldable, and yet determinants of foldability may be inapparent once the native state is reached. Insight has emerged from studies of diseases of protein misfolding, exemplified by monogenic diabetes mellitus due to mutations in proinsulin leading to endoplasmic reticulum stress and β-cell death. Cellular foldability of human proinsulin requires an invariant Phe within a conserved crevice at the receptor-binding surface (position B24). Any substitution, even related aromatic residue TyrB24, impairs insulin biosynthesis and secretion. As a seeming paradox, a monomeric TyrB24 insulin analog exhibits a native-like structure in solution with only a modest decrement in stability. Packing of TyrB24 is similar to that of PheB24, adjoining core cystine B19-A20 to seal the core; the analog also exhibits native self-assembly. Although affinity for the insulin receptor is decreased ∼20-fold, biological activities in cells and rats were within the range of natural variation. Together, our findings suggest that the invariance of PheB24 among vertebrate insulins and insulin-like growth factors reflects an essential role in enabling efficient protein folding, trafficking, and secretion, a function that is inapparent in native structures. In particular, we envision that the para-hydroxyl group of TyrB24 hinders pairing of cystine B19-A20 in an obligatory on-pathway folding intermediate. The absence of genetic variation at B24 and other conserved sites near this disulfide bridge-excluded due to β-cell dysfunction-suggests that insulin has evolved to the edge of foldability. Nonrobustness of a protein's fitness landscape underlies both a rare monogenic syndrome and "diabesity" as a pandemic disease of civilization.

Bioenergetic functions in subpopulations of heart mitochondria are preserved in a non-obese type 2 diabetes rat model (Goto-Kakizaki)

A distinct bioenergetic impairment of heart mitochondrial subpopulations in diabetic cardiomyopathy is associated with obesity; however, many type 2 diabetic (T2DM) patients with high-risk for cardiovascular disease are not obese. In the absence of obesity, it is unclear whether bioenergetic function in the subpopulations of mitochondria is affected in heart with T2DM. To address this issue, a rat model of non-obese T2DM was used to study heart mitochondrial energy metabolism, measuring bioenergetics and enzyme activities of the electron transport chain (ETC). Oxidative phosphorylation in the presence of substrates for ETC and ETC activities in both populations of heart mitochondria in T2DM rats were unchanged. Despite the preservation of mitochondrial function, aconitase activity in T2DM heart was reduced, suggesting oxidative stress in mitochondria. Our study indicate that metabolic function of heart mitochondria is unchanged in the face of oxidative stress and point to a critical role of obesity in T2DM cardiomyopathy.

Reassessment of an Innovative Insulin Analogue Excludes Protracted Action yet Highlights the Distinction between External and Internal Diselenide Bridges

Long-acting insulin analogues represent the most prescribed class of therapeutic proteins. An innovative design strategy was recently proposed: diselenide substitution of an external disulfide bridge. This approach exploited the distinctive physicochemical properties of selenocysteine (U). Relative to wild type (WT), Se-insulin[C7U^A , C7U^B ] was reported to be protected from proteolysis by insulin-degrading enzyme (IDE), predicting prolonged activity. Because of this strategy's novelty and potential clinical importance, we sought to validate these findings and test their therapeutic utility in an animal model of diabetes mellitus. Surprisingly, the analogue did not exhibit enhanced stability, and its susceptibility to cleavage by either IDE or a canonical serine protease (glutamyl endopeptidase Glu-C) was similar to WT. Moreover, the analogue's pharmacodynamic profile in rats was not prolonged relative to a rapid-acting clinical analogue (insulin lispro). Although [C7U^A , C7U^B ] does not confer protracted action, nonetheless its comparison to internal diselenide bridges promises to provide broad biophysical insight.

"Register-shift" insulin analogs uncover constraints of proteotoxicity in protein evolution

Globular protein sequences encode not only functional structures (the native state) but also protein foldability, i.e. a conformational search that is both efficient and robustly minimizes misfolding. Studies of mutations associated with toxic misfolding have yielded insights into molecular determinants of protein foldability. Of particular interest are residues that are conserved yet dispensable in the native state. Here, we exploited the mutant proinsulin syndrome (a major cause of permanent neonatal-onset diabetes mellitus) to investigate whether toxic misfolding poses an evolutionary constraint. Our experiments focused on an invariant aromatic motif (Phe^B24-Phe^B25-Tyr^B26) with complementary roles in native self-assembly and receptor binding. A novel class of mutations provided evidence that insulin can bind to the insulin receptor (IR) in two different modes, distinguished by a "register shift" in this motif, as visualized by molecular dynamics (MD) simulations. Register-shift variants are active but defective in cellular foldability and exquisitely susceptible to fibrillation in vitro Indeed, expression of the corresponding proinsulin variant induced endoplasmic reticulum stress, a general feature of the mutant proinsulin syndrome. Although not present among vertebrate insulin and insulin-like sequences, a prototypical variant ([Gly^B24]insulin) was as potent as WT insulin in a rat model of diabetes. Although in MD simulations the shifted register of receptor engagement is compatible with the structure and allosteric reorganization of the IR-signaling complex, our results suggest that this binding mode is associated with toxic misfolding and so is disallowed in evolution. The implicit threat of proteotoxicity limits sequence variation among vertebrate insulins and insulin-like growth factors.

Biophysical Elucidation of Amyloid Fibrillation Inhibition and Prevention of Secondary Nucleation by Cholic Acid: An Unexplored Function of Cholic Acid

Protein misfolding and its deviant self-assembly to converge into amyloid fibrils is associated with the perturbation of cellular functions and thus with debilitating neurodegenerative diseases including Alzheimer's disease, Parkinson's disease, etc. A great deal of research has already been carried out to discover a potential amyloid inhibitor that can slow down, prevent, or remodel toxic amyloids. In the present study with the help of a combination of biophysical, imaging, and computational techniques, we investigated the mechanism of interaction of cholic acid (CA), a primary bile acid, with human insulin and Aβ-42 and found CA to be effective in inhibiting amyloid formation. From ThT data, we inferred that CA encumbers amyloid fibrillation up to 90% chiefly by targeting elongation of fibrils with an insignificant effect on lag time, while in the case of Aβ-42, CA stabilizes the peptide in its native state preventing its fibrillation. Strikingly upon adding initially at the secondary nucleation stage, CA also detained the progression/growth of insulin fibrils. CA is unable to prevent the conformational changes completely during fibrillation but tends to resist and maintain an α helical structure up to a significant extent at a primary nucleation stage while reducing the β sheet rich content at the secondary nucleation stage. Moreover, CA treated samples exhibited reduced cytotoxicity and different morphology. Furthermore, the results obtained after molecular docking indicated that CA is interacting with insulin via hydrogen bonds. For future research, this study can be considered as preliminary research for the development of CA, a metabolite of our body, as a potential therapeutic agent against Alzheimer's disease without even stimulating the immunological responses.

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.

Biosynthesis, structure, and folding of the insulin precursor protein

Insulin synthesis in pancreatic β-cells is initiated as preproinsulin. Prevailing glucose concentrations, which oscillate pre- and postprandially, exert major dynamic variation in preproinsulin biosynthesis. Accompanying upregulated translation of the insulin precursor includes elements of the endoplasmic reticulum (ER) translocation apparatus linked to successful orientation of the signal peptide, translocation and signal peptide cleavage of preproinsulin-all of which are necessary to initiate the pathway of proper proinsulin folding. Evolutionary pressures on the primary structure of proinsulin itself have preserved the efficiency of folding ("foldability"), and remarkably, these evolutionary pressures are distinct from those protecting the ultimate biological activity of insulin. Proinsulin foldability is manifest in the ER, in which the local environment is designed to assist in the overall load of proinsulin folding and to favour its disulphide bond formation (while limiting misfolding), all of which is closely tuned to ER stress response pathways that have complex (beneficial, as well as potentially damaging) effects on pancreatic β-cells. Proinsulin misfolding may occur as a consequence of exuberant proinsulin biosynthetic load in the ER, proinsulin coding sequence mutations, or genetic predispositions that lead to an altered ER folding environment. Proinsulin misfolding is a phenotype that is very much linked to deficient insulin production and diabetes, as is seen in a variety of contexts: rodent models bearing proinsulin-misfolding mutants, human patients with Mutant INS-gene-induced Diabetes of Youth (MIDY), animal models and human patients bearing mutations in critical ER resident proteins, and, quite possibly, in more common variety type 2 diabetes.

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.

Solution structure of an ultra-stable single-chain insulin analog connects protein dynamics to a novel mechanism of receptor binding

Domain-minimized insulin receptors (IRs) have enabled crystallographic analysis of insulin-bound "micro-receptors." In such structures, the C-terminal segment of the insulin B chain inserts between conserved IR domains, unmasking an invariant receptor-binding surface that spans both insulin A and B chains. This "open" conformation not only rationalizes the inactivity of single-chain insulin (SCI) analogs (in which the A and B chains are directly linked), but also suggests that connecting (C) domains of sufficient length will bind the IR. Here, we report the high-resolution solution structure and dynamics of such an active SCI. The hormone's closed-to-open transition is foreshadowed by segmental flexibility in the native state as probed by heteronuclear NMR spectroscopy and multiple conformer simulations of crystallographic protomers as described in the companion article. We propose a model of the SCI's IR-bound state based on molecular-dynamics simulations of a micro-receptor complex. In this model, a loop defined by the SCI's B and C domains encircles the C-terminal segment of the IR α-subunit. This binding mode predicts a conformational transition between an ultra-stable closed state (in the free hormone) and an active open state (on receptor binding). Optimization of this switch within an ultra-stable SCI promises to circumvent insulin's complex global cold chain. The analog's biphasic activity, which serendipitously resembles current premixed formulations of soluble insulin and microcrystalline suspension, may be of particular utility in the developing world.

An ultra-stable single-chain insulin analog resists thermal inactivation and exhibits biological signaling duration equivalent to the native protein

Thermal degradation of insulin complicates its delivery and use. Previous efforts to engineer ultra-stable analogs were confounded by prolonged cellular signaling in vivo, of unclear safety and complicating mealtime therapy. We therefore sought an ultra-stable analog whose potency and duration of action on intravenous bolus injection in diabetic rats are indistinguishable from wild-type (WT) insulin. Here, we describe the structure, function, and stability of such an analog, a 57-residue single-chain insulin (SCI) with multiple acidic substitutions. Cell-based studies revealed native-like signaling properties with negligible mitogenic activity. Its crystal structure, determined as a novel zinc-free hexamer at 2.8 Å, revealed a native insulin fold with incomplete or absent electron density in the C domain; complementary NMR studies are described in the accompanying article. The stability of the analog (ΔG_U 5.0(±0.1) kcal/mol at 25 °C) was greater than that of WT insulin (3.3(±0.1) kcal/mol). On gentle agitation, the SCI retained full activity for >140 days at 45 °C and >48 h at 75 °C. These findings indicate that marked resistance to thermal inactivation in vitro is compatible with native duration of activity in vivo Further, whereas WT insulin forms large and heterogeneous aggregates above the standard 0.6 mm pharmaceutical strength, perturbing the pharmacokinetic properties of concentrated formulations, dynamic light scattering, and size-exclusion chromatography revealed only limited SCI self-assembly and aggregation in the concentration range 1-7 mm Such a combination of favorable biophysical and biological properties suggests that SCIs could provide a global therapeutic platform without a cold chain.

4S-Hydroxylation of Insulin at ProB28 Accelerates Hexamer Dissociation and Delays Fibrillation

Daily injections of insulin provide lifesaving benefits to millions of diabetics. But currently available prandial insulins are suboptimal: The onset of action is delayed by slow dissociation of the insulin hexamer in the subcutaneous space, and insulin forms amyloid fibrils upon storage in solution. Here we show, through the use of noncanonical amino acid mutagenesis, that replacement of the proline residue at position 28 of the insulin B-chain (ProB28) by (4S)-hydroxyproline (Hzp) yields an active form of insulin that dissociates more rapidly, and fibrillates more slowly, than the wild-type protein. Crystal structures of dimeric and hexameric insulin preparations suggest that a hydrogen bond between the hydroxyl group of Hzp and a backbone amide carbonyl positioned across the dimer interface may be responsible for the altered behavior. The effects of hydroxylation are stereospecific; replacement of ProB28 by (4R)-hydroxyproline (Hyp) causes little change in the rates of fibrillation and hexamer disassociation. These results demonstrate a new approach that fuses the concepts of medicinal chemistry and protein design, and paves the way to further engineering of insulin and other therapeutic proteins.

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

Insulin, a protein critical for metabolic homeostasis, provides a classical model for protein design with application to human health. Recent efforts to improve its pharmaceutical formulation demonstrated that iodination of a conserved tyrosine (TyrB26) enhances key properties of a rapid-acting clinical analog. Moreover, the broad utility of halogens in medicinal chemistry has motivated the use of hybrid quantum- and molecular-mechanical methods to study proteins. Here, we (i) undertook quantitative atomistic simulations of 3-[iodo-TyrB26]insulin to predict its structural features, and (ii) tested these predictions by X-ray crystallography. Using an electrostatic model of the modified aromatic ring based on quantum chemistry, the calculations suggested that the analog, as a dimer and hexamer, exhibits subtle differences in aromatic-aromatic interactions at the dimer interface. Aromatic rings (TyrB16, PheB24, PheB25, 3-I-TyrB26, and their symmetry-related mates) at this interface adjust to enable packing of the hydrophobic iodine atoms within the core of each monomer. Strikingly, these features were observed in the crystal structure of a 3-[iodo-TyrB26]insulin analog (determined as an R6 zinc hexamer). Given that residues B24-B30 detach from the core on receptor binding, the environment of 3-I-TyrB26 in a receptor complex must differ from that in the free hormone. Based on the recent structure of a "micro-receptor" complex, we predict that 3-I-TyrB26 engages the receptor via directional halogen bonding and halogen-directed hydrogen bonding as follows: favorable electrostatic interactions exploiting, respectively, the halogen's electron-deficient σ-hole and electronegative equatorial band. Inspired by quantum chemistry and molecular dynamics, such "halogen engineering" promises to extend principles of medicinal chemistry to proteins.

A Toolkit to Fit Nonbonded Parameters from and for Condensed Phase Simulations

The quality of atomistic simulations depends decisively on the accuracy of the underlying energy function (force field). Of particular importance for condensed-phase properties are nonbonded interactions, including the electrostatic and Lennard-Jones terms. Permanent atomic multipoles (MTPs) are an extension to common point-charge (PC) representations in atomistic simulations. MTPs are commonly determined from and fitted to an ab initio Electrostatic Potential (ESP), and Lennard-Jones (LJ) parameters are obtained from comparison of experimental and computed observables using molecular dynamics (MD) simulations. For this a set of thermodynamic observables such as density, heat of vaporization, and hydration free energy is chosen, to which the parametrization is fitted. The current work introduces a comprehensive computing environment (Fitting Wizard (FW)) for optimizing nonbonded interactions for atomistic force fields of different qualities. The FW supports fitting of standard PC-based force fields and more physically motivated multipolar (MTP) force fields. A broader study including 20 molecules ranging from N-methyl-acetamide and benzene to halogenated benzenes, phenols, anilines, and pyridines yields a root mean squared deviation for hydration free energies of 0.36 kcal/mol over a range of 8 kcal/mol. It is furthermore shown that PC-based force fields are not necessarily inferior compared to MTP parametrizations depending on the molecule considered.

Impact of Quadrupolar Electrostatics on Atoms Adjacent to the Sigma-Hole in Condensed-Phase Simulations

Halogenation is one of the cases for which advanced molecular simulation methods are mandatory for quantitative and predictive studies. The present work provides a systematic investigation of the importance of higher-order multipoles on specific sites of halobenzenes, other than the halogen, for static and dynamic properties in condensed-phase simulations. For that purpose, solute-solvent interactions using point charge (PC), multipole (MTP), and hybrid point charge/multipole (HYB) electrostatic models are analyzed in regions of halogen bonding and extended to regions of π orbitals of phenyl carbons. Using molecular dynamics simulations and quantum chemical methods, it is found that the sigma-hole does not only affect the halogen and the carbon bound to it but its effect extends to the carbons adjacent to the CX bond. This effect increases with the magnitude of the positive potential of the sigma-hole. With the MTP and HYB3 models, all hydration free energies of the PhX compounds are reproduced within 0.1 kcal/mol. Analysis of pair distribution functions and hydration free energies of halogenated benzenes provides a microscopic explanation why "point charge"-based representations with off-site charges fail in reproducing thermodynamic properties of the sigma-hole. Application of the hybrid models to study protein-ligand binding demonstrates both their accuracy and computational efficiency.

Aromatic anchor at an invariant hormone-receptor interface: function of insulin residue B24 with application to protein design

Crystallographic studies of insulin bound to fragments of the insulin receptor have recently defined the topography of the primary hormone-receptor interface. Here, we have investigated the role of Phe(B24), an invariant aromatic anchor at this interface and site of a human mutation causing diabetes mellitus. An extensive set of B24 substitutions has been constructed and tested for effects on receptor binding. Although aromaticity has long been considered a key requirement at this position, Met(B24) was found to confer essentially native affinity and bioactivity. Molecular modeling suggests that this linear side chain can serve as an alternative hydrophobic anchor at the hormone-receptor interface. These findings motivated further substitution of Phe(B24) by cyclohexanylalanine (Cha), which contains a nonplanar aliphatic ring. Contrary to expectations, [Cha(B24)]insulin likewise exhibited high activity. Furthermore, its resistance to fibrillation and the rapid rate of hexamer disassembly, properties of potential therapeutic advantage, were enhanced. The crystal structure of the Cha(B24) analog, determined as an R6 zinc-stabilized hexamer at a resolution of 1.5 Å, closely resembles that of wild-type insulin. The nonplanar aliphatic ring exhibits two chair conformations with partial occupancies, each recapitulating the role of Phe(B24) at the dimer interface. Together, these studies have defined structural requirements of an anchor residue within the B24-binding pocket of the insulin receptor; similar molecular principles are likely to pertain to insulin-related growth factors. Our results highlight in particular the utility of nonaromatic side chains as probes of the B24 pocket and suggest that the nonstandard Cha side chain may have therapeutic utility.