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

Se-Glargine: Chemical Synthesis of a Basal Insulin Analogue Stabilized by an Internal Diselenide Bridge

Insulin has long provided a model for studies of protein folding and stability, enabling enhanced treatment of diabetes mellitus via analogue design. We describe the chemical synthesis of a basal insulin analogue stabilized by substitution of an internal cystine (A6-A11) by a diselenide bridge. The studies focused on insulin glargine (formulated as Lantus® and Toujeo®; Sanofi). Prepared at pH 4 in the presence of zinc ions, glargine exhibits a shifted isoelectric point due to a basic B chain extension (ArgB31 -ArgB32 ). Subcutaneous injection leads to pH-dependent precipitation of a long-lived depot. Pairwise substitution of CysA6 and CysA11 by selenocysteine was effected by solid-phase peptide synthesis; the modified A chain also contained substitution of AsnA21 by Gly, circumventing acid-catalyzed deamidation. Although chain combination of native glargine yielded negligible product, in accordance with previous synthetic studies, the pairwise selenocysteine substitution partially rescued this reaction: substantial product was obtained through repeated combination, yielding a stabilized insulin analogue. This strategy thus exploited both (a) the unique redox properties of selenocysteine in protein folding and (b) favorable packing of an internal diselenide bridge in the native state, once achieved. Such rational optimization of protein folding and stability may be generalizable to diverse disulfide-stabilized proteins of therapeutic interest.

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.

The three-dimensional structure of insulin and its receptor

Insulin is a peptide hormone essential for maintaining normal blood glucose levels. Individuals unable to secrete sufficient insulin or not able to respond properly to insulin develop diabetes. Since the discovery of insulin its structure and function has been intensively studied with the aim to develop effective diabetes treatments. The three-dimensional crystal structure of this 51 amino acid peptide paved the way for discoveries, outlined in this review, of determinants important for receptor binding and hormone stability that have been instrumental in development of insulin analogs used in the clinic today. Important for the future development of effective diabetes treatments will be a detailed understanding of the insulin receptor structure and function. Determination of the three-dimensional structure of the insulin receptor, a receptor tyrosine kinase, proved challenging but with the recent advent of high-resolution cryo-electron microscopy significant progress has been made. There are now >40 structures of the insulin:insulin receptor complex deposited in the Protein Data Bank. From these structures we have a detailed picture of how insulin binds and activates the receptor. Still lacking are details of the initial binding events and the exact sequence of structural changes within the receptor and insulin. In this review, the focus will be on the most recent structural studies of insulin:insulin receptor complexes and how they have contributed to the current understanding of insulin receptor activation and signaling outcome. Molecular mechanisms underlying insulin receptor signaling bias emerging from the latest structures are described.

Determinants of IGF-II influencing stability, receptor binding and activation

Insulin like growth factor II (IGF-II) is involved in metabolic and mitogenic signalling in mammalian cells and plays important roles in normal fetal development and postnatal growth. It is structurally similar to insulin and binds not only with high affinity to the type 1 insulin-like growth factor receptor (IGF-1R) but also to the insulin receptor isoform A (IR-A). As IGF-II expression is commonly upregulated in cancer and its signalling promotes cancer cell survival, an antagonist that blocks IGF-II action without perturbing insulin signalling would be invaluable. The high degree of structural homology between the IR and IGF-1R makes selectively targeting either receptor in the treatment of IGF-II-dependent cancers very challenging. However, there are sequence differences between insulin and IGF-II that convey receptor selectivity and influence binding affinity and signalling outcome. Insulin residue YB16 is a key residue involved in maintaining insulin stability, dimer formation and IR binding. Mutation of this residue to glutamine (as found in IGF-II) results in reduced binding affinity. In this study we sought to determine if the equivalent residue Q18 in IGF-II plays a similar role. We show through site-directed mutagenesis of Q18 that this residue contributes to IGF-II structural integrity, selectivity of IGF-1R/IR binding, but surprisingly does not influence IR-A signalling activation. These findings provide insights into a unique IGF-II residue that can influence receptor binding specificity whilst having little influence on signalling outcome.

Symmetric and asymmetric receptor conformation continuum induced by a new insulin

Cone snail venoms contain a wide variety of bioactive peptides, including insulin-like molecules with distinct structural features, binding modes and biochemical properties. Here, we report an active humanized cone snail venom insulin with an elongated A chain and a truncated B chain, and use cryo-electron microscopy (cryo-EM) and protein engineering to elucidate its interactions with the human insulin receptor (IR) ectodomain. We reveal how an extended A chain can compensate for deletion of B-chain residues, which are essential for activity of human insulin but also compromise therapeutic utility by delaying dissolution from the site of subcutaneous injection. This finding suggests approaches to developing improved therapeutic insulins. Curiously, the receptor displays a continuum of conformations from the symmetric state to a highly asymmetric low-abundance structure that displays coordination of a single humanized venom insulin using elements from both of the previously characterized site 1 and site 2 interactions.

Peptide Model of the Mutant Proinsulin Syndrome. I. Design and Clinical Correlation

The mutant proinsulin syndrome is a monogenic cause of diabetes mellitus due to toxic misfolding of insulin's biosynthetic precursor. Also designated mutant INS-gene induced diabetes of the young (MIDY), this syndrome defines molecular determinants of foldability in the endoplasmic reticulum (ER) of β-cells. Here, we describe a peptide model of a key proinsulin folding intermediate and variants containing representative clinical mutations; the latter perturb invariant core sites in native proinsulin (LeuB15→Pro, LeuA16→Pro, and PheB24→Ser). The studies exploited a 49-residue single-chain synthetic precursor (designated DesDi), previously shown to optimize in vitro efficiency of disulfide pairing. Parent and variant peptides contain a single disulfide bridge (cystine B19-A20) to provide a model of proinsulin's first oxidative folding intermediate. The peptides were characterized by circular dichroism and redox stability in relation to effects of the mutations on (a) in vitro foldability of the corresponding insulin analogs and (b) ER stress induced in cell culture on expression of the corresponding variant proinsulins. Striking correlations were observed between peptide biophysical properties, degree of ER stress and age of diabetes onset (neonatal or adolescent). Our findings suggest that age of onset reflects the extent to which nascent structure is destabilized in proinsulin's putative folding nucleus. We envisage that such peptide models will enable high-resolution structural studies of key folding determinants and in turn permit molecular dissection of phenotype-genotype relationships in this monogenic diabetes syndrome. Our companion study (next article in this issue) employs two-dimensional heteronuclear NMR spectroscopy to define site-specific perturbations in the variant peptides.

Structural Lessons From the Mutant Proinsulin Syndrome

Insight into folding mechanisms of proinsulin has been provided by analysis of dominant diabetes-associated mutations in the human insulin gene (INS). Such mutations cause pancreatic β-cell dysfunction due to toxic misfolding of a mutant proinsulin and impairment in trans of wild-type insulin secretion. Anticipated by the "Akita" mouse (a classical model of monogenic diabetes mellitus; DM), this syndrome illustrates the paradigm endoreticulum (ER) stress leading to intracellular proteotoxicity. Diverse clinical mutations directly or indirectly perturb native disulfide pairing leading to protein misfolding and aberrant aggregation. Although most introduce or remove a cysteine (Cys; leading in either case to an unpaired thiol group), non-Cys-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the hormone's evolution has been constrained not only by structure-function relationships, but also by the susceptibility of its single-chain precursor to impaired foldability. An intriguing hypothesis posits that INS overexpression in response to peripheral insulin resistance likewise leads to chronic ER stress and β-cell dysfunction in the natural history of non-syndromic Type 2 DM. Cryptic contributions of conserved residues to folding efficiency, as uncovered by rare genetic variants, define molecular links between biophysical principles and the emerging paradigm of Darwinian medicine: Biosynthesis of proinsulin at the edge of non-foldability provides a key determinant of "diabesity" as a pandemic disease of civilization.

Insertion of a synthetic switch into insulin provides metabolite-dependent regulation of hormone-receptor activation

Insulin-signaling requires conformational change: whereas the free hormone and its receptor each adopt autoinhibited conformations, their binding leads to structural reorganization. To test the functional coupling between insulin's "hinge opening" and receptor activation, we inserted an artificial ligand-dependent switch into the insulin molecule. Ligand-binding disrupts an internal tether designed to stabilize the hormone's native closed and inactive conformation, thereby enabling productive receptor engagement. This scheme exploited a diol sensor (meta-fluoro-phenylboronic acid at GlyA1) and internal diol (3,4-dihydroxybenzoate at LysB28). The sensor recognizes monosaccharides (fructose > glucose). Studies of insulin-signaling in human hepatoma-derived cells (HepG2) demonstrated fructose-dependent receptor autophosphorylation leading to appropriate downstream signaling events, including a specific kinase cascade and metabolic gene regulation (gluconeogenesis and lipogenesis). Addition of glucose (an isomeric ligand with negligible sensor affinity) did not activate the hormone. Similarly, metabolite-regulated signaling was not observed in control studies of 1) an unmodified insulin analog or 2) an analog containing a diol sensor without internal tethering. Although secondary structure (as probed by circular dichroism) was unaffected by ligand-binding, heteronuclear NMR studies revealed subtle local and nonlocal monosaccharide-dependent changes in structure. Insertion of a synthetic switch into insulin has thus demonstrated coupling between hinge-opening and allosteric holoreceptor signaling. In addition to this foundational finding, our results provide proof of principle for design of a mechanism-based metabolite-responsive insulin. In particular, replacement of the present fructose sensor by an analogous glucose sensor may enable translational development of a "smart" insulin analog to mitigate hypoglycemic risk in diabetes therapy.

Diabetes mellitus due to toxic misfolding of proinsulin variants

BACKGROUND

Dominant mutations in the human insulin gene (INS) lead to pancreatic β-cell dysfunction and diabetes mellitus (DM) due to toxic misfolding of a mutant proinsulin. Analogous to a classical mouse model of monogenic DM ("Akita"), this syndrome highlights the susceptibility of β-cells to endoreticulum (ER) stress due to protein misfolding and aberrant aggregation.

SCOPE OF REVIEW

Diverse clinical mutations directly or indirectly perturb native disulfide pairing. Whereas most introduce or remove a cysteine (Cys; leading in either case to an unpaired thiol group), non-Cys-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the hormone's evolution has been constrained not only by structure-function relationships but also by the susceptibility of its single-chain precursor to impaired foldability. An intriguing hypothesis posits that INS overexpression in response to peripheral insulin resistance likewise leads to chronic ER stress and β-cell dysfunction in the natural history of nonsyndromic Type 2 DM.

MAJOR CONCLUSIONS

Cryptic contributions of conserved residues to folding efficiency, as uncovered by rare genetic variants, define molecular links between biophysical principles and the emerging paradigm of Darwinian medicine: Biosynthesis of proinsulin at the edge of nonfoldability provides a key determinant of "diabesity" as a pandemic disease of civilization.

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.

A structurally minimized yet fully active insulin based on cone-snail venom insulin principles

Human insulin and its current therapeutic analogs all show propensity, albeit varyingly, to self-associate into dimers and hexamers, which delays their onset of action and makes blood glucose management difficult for people with diabetes. Recently, we described a monomeric, insulin-like peptide in cone-snail venom with moderate human insulin-like bioactivity. Here, with insights from structural biology studies, we report the development of mini-Ins-a human des-octapeptide insulin analog-as a structurally minimal, full-potency insulin. Mini-Ins is monomeric and, despite the lack of the canonical B-chain C-terminal octapeptide, has similar receptor binding affinity to human insulin. Four mutations compensate for the lack of contacts normally made by the octapeptide. Mini-Ins also has similar in vitro insulin signaling and in vivo bioactivities to human insulin. The full bioactivity of mini-Ins demonstrates the dispensability of the PheB24-PheB25-TyrB26 aromatic triplet and opens a new direction for therapeutic insulin development.

"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.

A thing of beauty: Structure and function of insulin's "aromatic triplet"

The classical crystal structure of insulin was determined in 1969 by D.C. Hodgkin et al. following a 35-year program of research. This structure depicted a hexamer remarkable for its self-assembly as a zinc-coordinated trimer of dimer. Prominent at the dimer interface was an "aromatic triplet" of conserved residues at consecutive positions in the B chain: PheB24 , PheB25 and TyrB26 . The elegance of this interface inspired the Oxford team to poetry: "A thing of beauty is a joy forever" (John Keats as quoted by Blundell, T.L., et al. Advances in Protein Chemistry 26:279-286 [1972]). Here, we revisit this aromatic triplet in light of recent advances in the structural biology of insulin bound as a monomer to fragments of the insulin receptor. Such co-crystal structures have defined how these side chains pack at the primary hormone-binding surface of the receptor ectodomain. On receptor binding, the B-chain β-strand (residues B24-B28) containing the aromatic triplet detaches from the α-helical core of the hormone. Whereas TyrB26 lies at the periphery of the receptor interface and may functionally be replaced by a diverse set of substitutions, PheB24 and PheB25 engage invariant elements of receptor domains L1 and αCT. These critical contacts were anticipated by the discovery of diabetes-associated mutations at these positions by Donald Steiner et al. at the University of Chicago. Conservation of PheB24 , PheB25 and TyrB26 among vertebrate insulins reflects the striking confluence of structure-based evolutionary constraints: foldability, protective self-assembly and hormonal activity.

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.

Computational study of the activity, dynamics, energetics and conformations of insulin analogues using molecular dynamics simulations: Application to hyperinsulinemia and the critical residue B26

Due to the increasing prevalence of diabetes, finding therapeutic analogues for insulin has become an urgent issue. While many experimental studies have been performed towards this end, they have limited scope to examine all aspects of the effect of a mutation. Computational studies can help to overcome these limitations, however, relatively few studies that focus on insulin analogues have been performed to date. Here, we present a comprehensive computational study of insulin analogues-three mutant insulins that have been identified with hyperinsulinemia and three mutations on the critical B26 residue that exhibit similar binding affinity to the insulin receptor-using molecular dynamics simulations with the aim of predicting how mutations of insulin affect its activity, dynamics, energetics and conformations. The time evolution of the conformers is studied in long simulations. The probability density function and potential of mean force calculations are performed on each insulin analogue to unravel the effect of mutations on the dynamics and energetics of insulin activation. Our conformational study can decrypt the key features and molecular mechanisms that are responsible for an enhanced or reduced activity of an insulin analogue. We find two key results: 1) hyperinsulinemia may be due to the drastically reduced activity (and binding affinity) of the mutant insulins. 2) Y26_BS and Y26_BE are promising therapeutic candidates for insulin as they are more active than WT-insulin. The analysis in this work can be readily applied to any set of mutations on insulin to guide development of more effective therapeutic analogues.

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.

Elucidating the Activation Mechanism of the Insulin-Family Proteins with Molecular Dynamics Simulations

The insulin-family proteins bind to their own receptors, but insulin-like growth factor II (IGF-II) can also bind to the A isoform of the insulin receptor (IR-A), activating unique and alternative signaling pathways from those of insulin. Although extensive studies of insulin have revealed that its activation is associated with the opening of the B chain-C terminal (BC-CT), the activation mechanism of the insulin-like growth factors (IGFs) still remains unknown. Here, we present the first comprehensive study of the insulin-family proteins comparing their activation process and mechanism using molecular dynamics simulations to reveal new insights into their specificity to the insulin receptor. We have found that all the proteins appear to exhibit similar stochastic dynamics in their conformational change to an active state. For the IGFs, our simulations show that activation involves two opening locations: the opening of the BC-CT section away from the core, similar to insulin; and the additional opening of the BC-CT section away from the C domain. Furthermore, we have found that these two openings occur simultaneously in IGF-I, but not in IGF-II, where they can occur independently. This suggests that the BC-CT section and the C domain behave as a unified domain in IGF-I, but as two independent domains in IGF-II during the activation process, implying that the IGFs undergo different activation mechanisms for receptor binding. The probabilities of the active and inactive states of the proteins suggest that IGF-II is hyperactive compared to IGF-I. The hinge residue and the hydrophobic interactions in the core are found to play a critical role in the stability and activity of IGFs. Overall, our simulations have elucidated the crucial differences and similarities in the activation mechanisms of the insulin-family proteins, providing new insights into the molecular mechanisms responsible for the observed differences between IGF-I and IGF-II in receptor binding.

Contribution of TyrB26 to the Function and Stability of Insulin: STRUCTURE-ACTIVITY RELATIONSHIPS AT A CONSERVED HORMONE-RECEPTOR INTERFACE

Crystallographic studies of insulin bound to receptor domains have defined the primary hormone-receptor interface. We investigated the role of Tyr(B26), a conserved aromatic residue at this interface. To probe the evolutionary basis for such conservation, we constructed 18 variants at B26. Surprisingly, non-aromatic polar or charged side chains (such as Glu, Ser, or ornithine (Orn)) conferred high activity, whereas the weakest-binding analogs contained Val, Ile, and Leu substitutions. Modeling of variant complexes suggested that the B26 side chains pack within a shallow depression at the solvent-exposed periphery of the interface. This interface would disfavor large aliphatic side chains. The analogs with highest activity exhibited reduced thermodynamic stability and heightened susceptibility to fibrillation. Perturbed self-assembly was also demonstrated in studies of the charged variants (Orn and Glu); indeed, the Glu(B26) analog exhibited aberrant aggregation in either the presence or absence of zinc ions. Thus, although Tyr(B26) is part of insulin's receptor-binding surface, our results suggest that its conservation has been enjoined by the aromatic ring's contributions to native stability and self-assembly. We envisage that such classical structural relationships reflect the implicit threat of toxic misfolding (rather than hormonal function at the receptor level) as a general evolutionary determinant of extant protein sequences.

Molecular characterization of insulin from squamate reptiles reveals sequence diversity and possible adaptive evolution

The Squamata are the most adaptive and prosperous group among ectothermic amniotes, reptiles, due to their species-richness and geographically wide habitat. Although the molecular mechanisms underlying their prosperity remain largely unknown, unique features have been reported from hormones that regulate energy metabolism. Insulin, a central anabolic hormone, is one such hormone, as its roles and effectiveness in regulation of blood glucose levels remain to be examined in squamates. In the present study, cDNAs coding for insulin were isolated from multiple species that represent various groups of squamates. The deduced amino acid sequences showed a high degree of divergence, with four lineages showing obviously higher number of amino acid substitutions than most of vertebrates, from teleosts to mammals. Among 18 sites presented to comprise the two receptor binding surfaces (one with 12 sites and the other with 6 sites), substitutions were observed in 13 sites. Among them was the substitution of HisB10, which results in the loss of the ability to hexamerize. Furthermore, three of these substitutions were reported to increase mitogenicity in human analogues. These substitutions were also reported from insulin of hystricomorph rodents and agnathan fishes, whose mitogenic potency have been shown to be increased. The estimated value of the non-synonymous-to-synonymous substitution ratio (ω) for the Squamata clade was larger than those of the other reptiles and aves. Even higher values were estimated for several lineages among squamates. These results, together with the regulatory mechanisms of digestion and nutrient assimilation in squamates, suggested a possible adaptive process through the molecular evolution of squamate INS. Further studies on the roles of insulin, in relation to the physiological and ecological traits of squamate species, will provide an insight into the molecular mechanisms that have led to the adaptivity and prosperity of squamates.

Molecular Dynamics Simulations of Insulin: Elucidating the Conformational Changes that Enable Its Binding

A sequence of complex conformational changes is required for insulin to bind to the insulin receptor. Recent experimental evidence points to the B chain C-terminal (BC-CT) as the location of these changes in insulin. Here, we present molecular dynamics simulations of insulin that reveal new insights into the structural changes occurring in the BC-CT. We find three key results: 1) The opening of the BC-CT is inherently stochastic and progresses through an open and then a “wide-open” conformation—the wide-open conformation is essential for receptor binding, but occurs only rarely. 2) The BC-CT opens with a zipper-like mechanism, with a hinge at the Phe24 residue, and is maintained in the dominant closed/inactive state by hydrophobic interactions of the neighboring Tyr26, the critical residue where opening of the BC-CT (activation of insulin) is initiated. 3) The mutation Y26N is a potential candidate as a therapeutic insulin analogue. Overall, our results suggest that the binding of insulin to its receptor is a highly dynamic and stochastic process, where initial docking occurs in an open conformation and full binding is facilitated through interactions of insulin receptor residues with insulin in its wide-open conformation.