Alanine scanning mutagenesis of insulin

Molecular engineering of insulin for recombinant expression in yeast

Since the first administration of insulin to a person with diabetes in 1922, scientific contributions from academia and industry have improved insulin therapy and access. The pharmaceutical need for insulin is now more than 40 tons annually, half of which is produced by recombinant secretory expression in Saccharomyces cerevisiae. We discuss how, in this yeast species, adaptation of insulin precursors by removable structural elements is pivotal for efficient secretory expression. The technologies reviewed have been implemented at industrial scale and are seminal for the supply of human insulin and insulin analogues to people with diabetes now and in the future. Engineering of a target protein with removable structural elements may provide a general approach to yield optimisation.

Clues to the Design of Aggregation-Resistant Insulin from Proline Scanning of Highly Amyloidogenic Peptides Derived from the N-Terminal Segment of the A-Chain

Insulin aggregation poses a significant problem in pharmacology and medicine as it occurs during prolonged storage of the hormone and in vivo at insulin injection sites. We have recently shown that dominant forces driving the self-assembly of insulin fibrils are likely to arise from intermolecular interactions involving the N-terminal segment of the A-chain (ACC1-13). Here, we study how proline substitutions within the pilot GIVEQ sequence of this fragment affect its propensity to aggregate in both neutral and acidic environments. In a reasonable agreement with in silico prediction based on the Cordax algorithm, proline substitutions at positions 3, 4, and 5 turn out to be very effective in preventing aggregation according to thioflavin T-fluorescence-based kinetic assay, infrared spectroscopy, and atomic force microscopy (AFM). Since the valine and glutamate side chains within this segment are strongly involved in the interactions with the insulin receptor, we have focused on the possible implications of the Q → P substitution for insulin's stability and interactions with the receptor. To this end, comparative molecular dynamics (MD) simulations of the Q5P mutant and wild-type insulin were carried out for both free and receptor-bound (site 1) monomers. The results point to a mild destabilization of the mutant vis à vis the wild-type monomer, as well as partial preservation of key contacts in the complex between Q5P insulin and the receptor. We discuss the implications of these findings in the context of the design of aggregation-resistant insulin analogues retaining hormonal activity.

Engineering of human tryptophan hydroxylase 2 for efficient synthesis of 5-hydroxytryptophan

L-Tryptophan hydroxylation catalyzed by tryptophan hydroxylase (TPH) presents a promising method for synthesizing 5-hydroxytryptophan (5-HTP), yet the limited activity of wild-type human TPH2 restricts its application. A high-activity mutant, MT10 (H318E/H323E), was developed through semi-rational active site saturation testing (CAST) of wild-type TPH2, exhibiting a 2.85-fold increase in kcat/Km over the wild type, thus enhancing catalytic efficiency. Two biotransformation systems were developed, including an in vitro one-pot system and a Whole-Cell Catalysis System (WCCS). In the WCCS, MT10 achieved a conversion rate of only 31.5 % within 32 h. In the one-pot reaction, MT10 converted 50 mM L-tryptophan to 44.5 mM 5-HTP within 8 h, achieving an 89 % conversion rate, outperforming the M1 (NΔ143/CΔ26) variant. Molecular dynamics simulations indicated enhanced interactions of MT10 with the substrate, suggesting improved binding affinity and system stability. This study offers an effective approach for the efficient production of 5-HTP.

pan-MHC and cross-Species Prediction of T Cell Receptor-Antigen Binding

Profiling the binding of T cell receptors (TCRs) of T cells to antigenic peptides presented by MHC proteins is one of the most important unsolved problems in modern immunology. Experimental methods to probe TCR-antigen interactions are slow, labor-intensive, costly, and yield moderate throughput. To address this problem, we developed pMTnet-omni, an Artificial Intelligence (AI) system based on hybrid protein sequence and structure information, to predict the pairing of TCRs of αβ T cells with peptide-MHC complexes (pMHCs). pMTnet-omni is capable of handling peptides presented by both class I and II pMHCs, and capable of handling both human and mouse TCR-pMHC pairs, through information sharing enabled this hybrid design. pMTnet-omni achieves a high overall Area Under the Curve of Receiver Operator Characteristics (AUROC) of 0.888, which surpasses competing tools by a large margin. We showed that pMTnet-omni can distinguish binding affinity of TCRs with similar sequences. Across a range of datasets from various biological contexts, pMTnet-omni characterized the longitudinal evolution and spatial heterogeneity of TCR-pMHC interactions and their functional impact. We successfully developed a biomarker based on pMTnet-omni for predicting immune-related adverse events of immune checkpoint inhibitor (ICI) treatment in a cohort of 57 ICI-treated patients. pMTnet-omni represents a major advance towards developing a clinically usable AI system for TCR-pMHC pairing prediction that can aid the design and implementation of TCR-based immunotherapeutics.

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.

Insulin receptor Arg717 and IGF-1 receptor Arg704 play a key role in ligand binding and in receptor activation

The insulin receptor (IR, with its isoforms IR-A and IR-B) and the insulin-like growth factor 1 receptor (IGF-1R) are related tyrosine kinase receptors. Recently, the portfolio of solved hormone-receptor structures has grown extensively thanks to advancements in cryo-electron microscopy. However, the dynamics of how these receptors transition between their inactive and active state are yet to be fully understood. The C-terminal part of the alpha subunit (αCT) of the receptors is indispensable for the formation of the hormone-binding site. We mutated the αCT residues Arg717 and His710 of IR-A and Arg704 and His697 of IGF-1R. We then measured the saturation binding curves of ligands on the mutated receptors and their ability to become activated. Mutations of Arg704 and His697 to Ala in IGF-1R decreased the binding of IGF-1. Moreover, the number of binding sites for IGF-1 on the His697 IGF-1R mutant was reduced to one-half, demonstrating the presence of two binding sites. Both mutations of Arg717 and His710 to Ala in IR-A inactivated the receptor. We have proved that Arg717 is important for the binding of insulin to its receptor, which suggests that Arg717 is a key residue for the transition to the active conformation.

A stepwise activation model for the insulin receptor

The binding of insulin to the insulin receptor (IR) triggers a cascade of receptor conformational changes and autophosphorylation, leading to the activation of metabolic and mitogenic pathways. Recent advances in the structural and functional analyses of IR have revealed the conformations of the extracellular domains of the IR in inactive and fully activated states. However, the early activation mechanisms of this receptor remain poorly understood. The structures of partially activated IR in complex with aptamers provide clues for understanding the initial activation mechanism. In this review, we discuss the structural and functional features of IR complexed with various ligands and propose a model to explain the sequential activation mechanism. Moreover, we discuss the structures of IR complexed with biased agonists that selectively activate metabolic pathways and provide insights into the design of selective agonists and their clinical implications.

The Activation Mechanism of the Insulin Receptor: A Structural Perspective

The insulin receptor (IR) is a type II receptor tyrosine kinase that plays essential roles in metabolism, growth, and proliferation. Dysregulation of IR signaling is linked to many human diseases, such as diabetes and cancers. The resolution revolution in cryo-electron microscopy has led to the determination of several structures of IR with different numbers of bound insulin molecules in recent years, which have tremendously improved our understanding of how IR is activated by insulin. Here, we review the insulin-induced activation mechanism of IR, including (a) the detailed binding modes and functions of insulin at site 1 and site 2 and (b) the insulin-induced structural transitions that are required for IR activation. We highlight several other key aspects of the activation and regulation of IR signaling and discuss the remaining gaps in our understanding of the IR activation mechanism and potential avenues of future research.

Targeting the insulin receptor with hormone and peptide dimers

Insulin is a key hormone involved in the regulation of overall energetic homeostasis of the organism. The dimeric character of the receptor for insulin evokes ideas about its activation or inhibition with peptide dimers that could either trigger or block the structural transition of the insulin receptor, leading to its activation. Herewith, we present the chemical engineering and biological characterization of several series of insulin dimers or dimers of specific peptides that should be able to bind receptors for insulin or insulin growth factor 1. The hormones or peptides in the dimers were interconnected with different linkers, consisting of triazole moieties and 3, 6, 8, 11, or 23 polyethylene glycol units. The prepared dimers were weaker in binding to insulin receptors than human insulin. However, some of the insulin dimers showed preferential binding specificity toward the isoform A of the insulin receptor, and the insulin dimers also stimulated the insulin receptor more strongly than would be consistent with their binding affinities. Our results suggest that designing insulin dimers may be a promising strategy for modulating the ability of the hormone to activate the receptor or to alter its specificity toward insulin receptor isoforms.

Mutations at hypothetical binding site 2 in insulin and insulin-like growth factors 1 and 2

Elucidating how insulin and the related insulin-like growth factors 1 and 2 (IGF-1 and IGF-2) bind to their cellular receptors (IR and IGF-1R) and how the receptors are activated has been the holy grail for generations of scientists. However, deciphering the 3D structure of tyrosine kinase receptors and their hormone-bound complexes has been complicated by the flexible and dimeric nature of the receptors and the dynamic nature of their interaction with hormones. Therefore, mutagenesis of hormones and kinetic studies first became an important tool for studying receptor interactions. It was suggested that hormones could bind to receptors through two binding sites on the hormone surface called site 1 and site 2. A breakthrough in knowledge came with the solution of cryoelectron microscopy (cryoEM) structures of hormone-receptor complexes. In this chapter, we document in detail the mutagenesis of insulin, IGF-1, and IGF-2 with emphasis on modifications of the hypothetical binding site 2 in the hormones, and we discuss the results of structure-activity studies in light of recent cryoEM structures of hormone complexes with IR and IGF-1R.

Cryo-electron Microscopic Analysis of Single-Pass Transmembrane Receptors

Single-pass transmembrane receptors (SPTMRs) represent a diverse group of integral membrane proteins that are involved in many essential cellular processes, including signal transduction, cell adhesion, and transmembrane transport of materials. Dysregulation of the SPTMRs is linked with many human diseases. Despite extensive efforts in past decades, the mechanisms of action of the SPTMRs remain incompletely understood. One major hurdle is the lack of structures of the full-length SPTMRs in different functional states. Such structural information is difficult to obtain by traditional structural biology methods such as X-ray crystallography and nuclear magnetic resonance (NMR). The recent rapid development of single-particle cryo-electron microscopy (cryo-EM) has led to an exponential surge in the number of high-resolution structures of integral membrane proteins, including SPTMRs. Cryo-EM structures of SPTMRs solved in the past few years have tremendously improved our understanding of how SPTMRs function. In this review, we will highlight these progresses in the structural studies of SPTMRs by single-particle cryo-EM, analyze important structural details of each protein involved, and discuss their implications on the underlying mechanisms. Finally, we also briefly discuss remaining challenges and exciting opportunities in the field.

[The insulin receptor discovery is 50 years old - A review of achieved progress]

The isolation of insulin from the pancreas and its purification to a degree permitting its safe administration to type 1 diabetic patients were accomplished 100 years ago at the University of Toronto by Banting, Best, Collip and McLeod and constitute undeniably one of the major medical therapeutic revolutions, recognized by the attribution of the 1923 Nobel Prize in Physiology or Medicine to Banting and McLeod. The clinical spin off was immediate as well as the internationalization of insulin's commercial production. The outcomes regarding basic research were much slower, in particular regarding the molecular mechanisms of insulin action on its target cells. It took almost a half-century before the determination of the tri-dimensional structure of insulin in 1969 and the characterization of its cell receptor in 1970-1971. The demonstration that the insulin receptor is in fact an enzyme named tyrosine kinase came in the years 1982-1985, and the crystal structure of the intracellular kinase domain 10 years later. The crystal structure of the first intracellular kinase substrate (IRS-1) in 1991 paved the way for the elucidation of the intracellular signalling pathways but it took 15 more years to obtain the complete crystal structure of the extracellular receptor domain (without insulin) in 2006. Since then, the determination of the structure of the whole insulin-receptor complex in both the inactive and activated states has made considerable progress, not least due to recent improvement in the resolution power of cryo-electron microscopy. I will here review the steps in the development of the concept of hormone receptor, and of our knowledge of the structure and molecular mechanism of activation of the insulin receptor.

Molecular Aspects of Insulin Aggregation and Various Therapeutic Interventions

Protein aggregation leading to the formation of amyloid fibrils has various adverse effects on human health ranging from fatigue and numbness to organ failure and death in extreme cases. Insulin, a peptide hormone commonly used to treat diabetes, undergoes aggregation at the site of repeated injections in diabetic patients as well as during its industrial production and transport. The reduced bioavailability of insulin due to aggregation hinders the proper control of glucose levels in diabetic patients. Thus, it is necessary to develop rational approaches for inhibiting insulin aggregation, which in turn requires a detailed understanding of the mechanism of fibrillation. Given the relative simplicity of insulin and ease of access, insulin has also served as a model system for studying amyloids. Approaches to inhibit insulin aggregation have included the use of natural molecules, synthetic peptides or small molecules, and bacterial chaperone machinery. This review focuses on insulin aggregation with an emphasis on its mechanism, the structural features of insulin fibrils, and the reported inhibitors that act at different stages in the aggregation pathway. We discuss molecules that can serve as leads for improved inhibitors for use in commercial insulin formulations. We also discuss the aggregation propensity of fast- and slow-acting insulin biosimilars, commonly administered to diabetic patients. The development of better insulin aggregation inhibitors and insights into their mechanism of action will not only aid diabetic therapies, but also enhance our knowledge of protein amyloidosis.

Synergistic activation of the insulin receptor via two distinct sites

Insulin receptor (IR) signaling controls multiple facets of animal physiology. Maximally four insulins bind to IR at two distinct sites, termed site-1 and site-2. However, the precise functional roles of each binding event during IR activation remain unresolved. Here, we showed that IR incompletely saturated with insulin predominantly forms an asymmetric conformation and exhibits partial activation. IR with one insulin bound adopts a Γ-shaped conformation. IR with two insulins bound assumes a Ƭ-shaped conformation. One insulin binds at site-1 and another simultaneously contacts both site-1 and site-2 in the Ƭ-shaped IR dimer. We further show that concurrent binding of four insulins to sites-1 and -2 prevents the formation of asymmetric IR and promotes the T-shaped symmetric, fully active state. Collectively, our results demonstrate how the synergistic binding of multiple insulins promotes optimal IR activation.

Synthesis and Characterization of Phenylboronic Acid-Modified Insulin With Glucose-Dependent Solubility

Glucose-responsive insulin represents a promising approach to regulate blood glucose levels. We previously showed that attaching two fluorophenylboronic acid (FPBA) residues to the C-terminal B chain of insulin glargine led to glucose-dependent solubility. Herein, we demonstrated that relocating FPBA from B chain to A chain increased the baseline solubility without affecting its potency. Furthermore, increasing the number of FPBA groups led to increased glucose-dependent solubility.

Identification and Affinity Determination of Protein-Antibody and Protein-Aptamer Epitopes by Biosensor-Mass Spectrometry Combination

Analytical methods for molecular characterization of diagnostic or therapeutic targets have recently gained high interest. This review summarizes the combination of mass spectrometry and surface plasmon resonance (SPR) biosensor analysis for identification and affinity determination of protein interactions with antibodies and DNA-aptamers. The binding constant (KD) of a protein-antibody complex is first determined by immobilizing an antibody or DNA-aptamer on an SPR chip. A proteolytic peptide mixture is then applied to the chip, and following removal of unbound material by washing, the epitope(s) peptide(s) are eluted and identified by MALDI-MS. The SPR-MS combination was applied to a wide range of affinity pairs. Distinct epitope peptides were identified for the cardiac biomarker myoglobin (MG) both from monoclonal and polyclonal antibodies, and binding constants determined for equine and human MG provided molecular assessment of cross immunoreactivities. Mass spectrometric epitope identifications were obtained for linear, as well as for assembled ("conformational") antibody epitopes, e.g., for the polypeptide chemokine Interleukin-8. Immobilization using protein G substantially improved surface fixation and antibody stabilities for epitope identification and affinity determination. Moreover, epitopes were successfully determined for polyclonal antibodies from biological material, such as from patient antisera upon enzyme replacement therapy of lysosomal diseases. The SPR-MS combination was also successfully applied to identify linear and assembled epitopes for DNA-aptamer interaction complexes of the tumor diagnostic protein C-Met. In summary, the SPR-MS combination has been established as a powerful molecular tool for identification of protein interaction epitopes.

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.

Understanding insulin and its receptor from their three-dimensional structures

BACKGROUND

Insulin's discovery 100 years ago and its ongoing use since that time to treat diabetes belies the molecular complexity of its structure and that of its receptor. Advances in single-particle cryo-electron microscopy have over the past three years revolutionized our understanding of the atomic detail of insulin-receptor interactions.

SCOPE OF REVIEW

This review describes the three-dimensional structure of insulin and its receptor and details on how they interact. This review also highlights the current gaps in our structural understanding of the system.

MAJOR CONCLUSIONS

A near-complete picture has been obtained of the hormone receptor interactions, providing new insights into the kinetics of the interactions and necessitating a revision of the extant two-site cross-linking model of hormone receptor engagement. How insulin initially engages the receptor and the receptor's traversed trajectory as it undergoes conformational changes associated with activation remain areas for future investigation.

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.

Commemorating insulin's centennial: engineering insulin pharmacology towards physiology

The life-saving discovery of insulin in Toronto in 1921 is one of the most impactful achievements in medical history, at the time being hailed as a miracle treatment for diabetes. The insulin molecule itself, however, is poorly amenable as a pharmacological intervention, and the formidable challenge of optimizing insulin therapy has been ongoing for a century. We review early academic insights into insulin structure and its relation to self-association and receptor binding, as well as recombinant biotechnology, which have all been seminal for drug design. Recent developments have focused on combining genetic and chemical engineering with pharmaceutical optimization to generate ultra-rapid and ultra-long-acting, tissue-selective, or orally delivered insulin analogs. We further discuss these developments and propose that future scientific efforts in molecular engineering include realizing the dream of glucose-responsive insulin delivery.

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.

The Chemical Synthesis of Insulin: An Enduring Challenge

The pancreatic peptide hormone insulin, first discovered exactly 100 years ago, is essential for glycemic control and is used as a therapeutic for the treatment of type 1 and, increasingly, type 2 diabetes. With a worsening global diabetes epidemic and its significant health budget imposition, there is a great demand for new analogues possessing improved physical and functional properties. However, the chemical synthesis of insulin's intricate 51-amino acid, two-chain, three-disulfide bond structure, together with the poor physicochemical properties of both the individual chains and the hormone itself, has long represented a major challenge to organic chemists. This review provides a timely overview of the past efforts to chemically assemble this fascinating hormone using an array of strategies to enable both correct folding of the two chains and selective formation of disulfide bonds. These methods not only have contributed to general peptide synthesis chemistry and enabled access to the greatly growing numbers of insulin-like and cystine-rich peptides but also, today, enable the production of insulin at the synthetic efficiency levels of recombinant DNA expression methods. They have led to the production of a myriad of novel analogues with optimized structural and functional features and of the feasibility for their industrial manufacture.

Modulation of co-stimulatory signal from CD2-CD58 proteins by a grafted peptide

Peptides were designed to inhibit the protein-protein interaction of CD2 and CD58 to modulate the immune response. This work involved the design and synthesis of eight different peptides by replacing each amino acid residue in peptide 6 with alanine as well as grafting the peptide to the sunflower trypsin-inhibitor framework. From the alanine scanning studies, mutation at position 2 of the peptide was shown to result in increased potency to inhibit cell adhesion interactions. The most potent peptide from the alanine scanning was further studied for its detailed three-dimensional structure and binding to CD58 protein using surface plasmon resonance and flow cytometry. This peptide was used to graft to the sunflower trypsin inhibitor to improve the stability of the peptide. The grafted peptide, SFTI-a1, was further studied for its potency as well as its thermal, chemical, and enzymatic stability. The grafted peptide exhibited improved activity compared to our previously grafted peptide and was stable against thermal and enzymatic degradation.

Modulation of Insulin Sensitivity by Insulin-Degrading Enzyme

Insulin-degrading enzyme (IDE) is a highly conserved and ubiquitously expressed metalloprotease that degrades insulin and several other intermediate-size peptides. For many decades, IDE had been assumed to be involved primarily in hepatic insulin clearance, a key process that regulates availability of circulating insulin levels for peripheral tissues. Emerging evidence, however, suggests that IDE has several other important physiological functions relevant to glucose and insulin homeostasis, including the regulation of insulin secretion from pancreatic β-cells. Investigation of mice with tissue-specific genetic deletion of Ide in the liver and pancreatic β-cells (L-IDE-KO and B-IDE-KO mice, respectively) has revealed additional roles for IDE in the regulation of hepatic insulin action and sensitivity. In this review, we discuss current knowledge about IDE’s function as a regulator of insulin secretion and hepatic insulin sensitivity, both evaluating the classical view of IDE as an insulin protease and also exploring evidence for several non-proteolytic functions. Insulin proteostasis and insulin sensitivity have both been highlighted as targets controlling blood sugar levels in type 2 diabetes, so a clearer understanding the physiological functions of IDE in pancreas and liver could led to the development of novel therapeutics for the treatment of this disease.

Biological behaviors of mutant proinsulin contribute to the phenotypic spectrum of diabetes associated with insulin gene mutations

Insulin gene mutation is the second most common cause of neonatal diabetes (NDM). It is also one of the genes involved in maturity-onset diabetes of the young (MODY). We aim to investigate molecular behaviors of different INS gene variants that may correlate with the clinical spectrum of diabetes phenotypes. In this study, we concentrated on two previously uncharacterized MODY-causing mutants, proinsulin-p.Gly44Arg [G(B20)R] and p.Pro52Leu [P(B28)L] (a novel mutant identified in one French family), and an NDM causing proinsulin-p.(Cys96Tyr) [C(A7)Y]. We find that these proinsulin mutants exhibit impaired oxidative folding in the endoplasmic reticulum (ER) with blocked ER export, ER stress, and apoptosis. Importantly, the proinsulin mutants formed abnormal intermolecular disulfide bonds that not only involved the mutant proinsulin, but also the co-expressed WT-proinsulin, forming misfolded disulfide-linked proinsulin complexes. This impaired the intracellular trafficking of WT-proinsulin and limited the production of bioactive mature insulin. Notably, although all three mutants presented with similar defects in folding, trafficking, and dominant negative behavior, the degrees of these defects appeared to be different. Specifically, compared to MODY mutants G(B20)R and P(B28)L that partially affected folding and trafficking of co-expressed WT-proinsulin, the NDM mutant C(A7)Y resulted in an almost complete blockade of the ER export of WT-proinsulin, decreasing insulin production, inducing more severe ER stress and apoptosis. We thus demonstrate that differences in cell biological behaviors among different proinsulin mutants correlate with the spectrum of diabetes phenotypes caused by the different INS gene mutations.

Uncovering the basis of protein-protein interaction specificity with a combinatorially complete library

Protein-protein interaction specificity is often encoded at the primary sequence level. However, the contributions of individual residues to specificity are usually poorly understood and often obscured by mutational robustness, sequence degeneracy, and epistasis. Using bacterial toxin-antitoxin systems as a model, we screened a combinatorially complete library of antitoxin variants at three key positions against two toxins. This library enabled us to measure the effect of individual substitutions on specificity in hundreds of genetic backgrounds. These distributions allow inferences about the general nature of interface residues in promoting specificity. We find that positive and negative contributions to specificity are neither inherently coupled nor mutually exclusive. Further, a wild-type antitoxin appears optimized for specificity as no substitutions improve discrimination between cognate and non-cognate partners. By comparing crystal structures of paralogous complexes, we provide a rationale for our observations. Collectively, this work provides a generalizable approach to understanding the logic of molecular recognition.

Understanding IGF-II Action through Insights into Receptor Binding and Activation

The insulin-like growth factor (IGF) system regulates metabolic and mitogenic signaling through an intricate network of related receptors and hormones. IGF-II is one of several hormones within this system that primarily regulates mitogenic functions and is especially important during fetal growth and development. IGF-II is also found to be overexpressed in several cancer types, promoting growth and survival. It is also unique in the IGF system as it acts through both IGF-1R and insulin receptor isoform A (IR-A). Despite this, IGF-II is the least investigated ligand of the IGF system. This review will explore recent developments in IGF-II research including a structure of IGF-II bound to IGF-1R determined using cryo-electron microscopy (cryoEM). Comparisons are made with the structures of insulin and IGF-I bound to their cognate receptors. Finally discussed are outstanding questions in the mechanism of action of IGF-II with the goal of developing antagonists of IGF action in cancer.

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.

Deep2Full: Evaluating strategies for selecting the minimal mutational experiments for optimal computational predictions of deep mutational scan outcomes

Performing a complete deep mutational scan with all single point mutations may not be practical, and may not even be required, especially if predictive computational models can be developed. Computational models are however naive to cellular response in the myriads of assay-conditions. In a realistic paradigm of assay context-aware predictive hybrid models that combine minimal experimental data from deep mutational scans with structure, sequence information and computational models, we define and evaluate different strategies for choosing this minimal set. We evaluated the trivial strategy of a systematic reduction in the number of mutational studies from 85% to 15%, along with several others about the choice of the types of mutations such as random versus site-directed with the same 15% data completeness. Interestingly, the predictive capabilities by training on a random set of mutations and using a systematic substitution of all amino acids to alanine, asparagine and histidine (ANH) were comparable. Another strategy we explored, augmenting the training data with measurements of the same mutants at multiple assay conditions, did not improve the prediction quality. For the six proteins we analyzed, the bin-wise error in prediction is optimal when 50-100 mutations per bin are used in training the computational model, suggesting that good prediction quality may be achieved with a library of 500-1000 mutations.

Calculating the absolute binding free energy of the insulin dimer in an explicit solvent

Insulin is a significant hormone in the regulation of glucose level in the blood. Its monomers bind to each other to form dimers or hexamers through a complex process. To study the binding of the insulin dimer, we first calculate its absolute binding free energy by the steered molecular dynamics method and the confinement method based on a fictitious thermodynamic cycle. After considering some special correction terms, the final calculated binding free energy at 298 K is −8.97 ± 1.41 kcal mol⁻¹, which is close to the experimental value of −7.2 ± 0.8 kcal mol⁻¹. Furthermore, we discuss the important residue–residue interactions between the insulin monomers, including hydrophobic interactions, π–π interactions and hydrogen bond interactions. The analysis reveals five key residues, Vla^B12, Tyr^B16, Phe^B24, Phe^B25, and Tyr^B26, for the dimerization of the insulin. We also perform MM-PBSA calculations for the wild-type dimer and some mutants and study the roles of the key residues by the change of the binding energy of the insulin dimer.

In this paper, we calculate the absolute binding free energy of an insulin dimer by steered MD method. The result of −8.97 kcal mol⁻¹ is close to the experimental value −7.2 kcal mol⁻¹. We also analyze the residue–residue interactions.

Mutations at hypothetical binding site 2 in insulin and insulin-like growth factors 1 and 2 result in receptor- and hormone-specific responses

Information on how insulin and insulin-like growth factors 1 and 2 (IGF-1 and -2) activate insulin receptors (IR-A and -B) and the IGF-1 receptor (IGF-1R) is crucial for understanding the difference in the biological activities of these peptide hormones. Cryo-EM studies have revealed that insulin uses its binding sites 1 and 2 to interact with IR-A and have identified several critical residues in binding site 2. However, mutagenesis studies suggest that Ile-A10, Ser-A12, Leu-A13, and Glu-A17 also belong to insulin's site 2. Here, to resolve this discrepancy, we mutated these insulin residues and the equivalent residues in IGFs. Our findings revealed that equivalent mutations in the hormones can result in differential biological effects and that these effects can be receptor-specific. We noted that the insulin positions A10 and A17 are important for its binding to IR-A and IR-B and IGF-1R and that A13 is important only for IR-A and IR-B binding. The IGF-1/IGF-2 positions 51/50 and 54/53 did not appear to play critical roles in receptor binding, but mutations at IGF-1 position 58 and IGF-2 position 57 affected the binding. We propose that IGF-1 Glu-58 interacts with IGF-1R Arg-704 and belongs to IGF-1 site 1, a finding supported by the NMR structure of the less active Asp-58–IGF-1 variant. Computational analyses indicated that the aforementioned mutations can affect internal insulin dynamics and inhibit adoption of a receptor-bound conformation, important for binding to receptor site 1. We provide a molecular model and alternative hypotheses for how the mutated insulin residues affect activity.

Activation mechanism of the insulin receptor revealed by cryo-EM structure of the fully liganded receptor-ligand complex

Insulin signaling controls metabolic homeostasis. Here, we report the cryo-EM structure of full-length insulin receptor (IR) and insulin complex in the active state. This structure unexpectedly reveals that maximally four insulins can bind the 'T'-shaped IR dimer at four distinct sites related by 2-fold symmetry. Insulins 1 and 1' bind to sites 1 and 1', formed by L1 of one IR protomer and α-CT and FnIII-1 of the other. Insulins 2 and 2' bind to sites 2 and 2' on FnIII-1 of each protomer. Mutagenesis and cellular assays show that both sites 1 and 2 are required for optimal insulin binding and IR activation. We further identify a homotypic FnIII-2-FnIII-2 interaction in mediating the dimerization of membrane proximal domains in the active IR dimer. Our results indicate that binding of multiple insulins at two distinct types of sites disrupts the autoinhibited apo-IR dimer and stabilizes the active dimer.

A Disulfide Scan of Insulin by [3 + 1] Methodology Exhibits Site-Specific Influence on Bioactivity

Insulin is the principal hormone involved in the regulation of metabolism and has served a seminal role in the treatment of diabetes. Building upon advances in insulin synthetic methodology, we have developed a straightforward route to novel insulins containing a fourth disulfide bond in a [3 + 1] fashion establishing the first disulfide scan of the hormone. All the targeted analogs accommodated the constraint to demonstrate an unexpected conformational flexibility of native insulin. The bioactivity was established for the constrained (4-DS) and unconstrained (3-DS) analogs by in vitro methods, and extended to in vivo study for select peptides. We also identified residue B10 as a preferred anchor to introduce a tether that would regulate insulin bioactivity. We believe that the described [3 + 1] methodology might constitute the preferred approach for performing similar disulfide scanning in peptides that contain multiple disulfides.

Probing the Differential Dynamics of the Monomeric and Dimeric Insulin from Amide-I IR Spectroscopy

The monomer-dimer equilibrium for insulin is one of the essential steps in forming the receptor-binding competent monomeric form of the hormone. Despite this importance, the thermodynamic stability, in particular for modified insulins, is quite poorly understood, in part, due to experimental difficulties. This work explores one- and two-dimensional infrared spectroscopy in the range of the amide-I band for the hydrated monomeric and dimeric wild-type hormone. It is found that for the monomer the frequency fluctuation correlation function (FFCF) and the one-dimensional infrared spectra are position sensitive. The spectra of the -CO probes at the dimerization interface (residues Phe24, Phe25, and Tyr26) split and indicate an asymmetry despite the overall (formal) point symmetry of the dimer structure. Also, the decay times of the FFCF for the same -CO probe in the monomer and the dimer can differ by up to 1 order of magnitude, for example, for residue PheB24, which is solvent exposed for the monomer but at the interface for the dimer. The spectroscopic shifts correlate approximately with the average number of hydration waters and the magnitude of the FFCF at time zero. However, this correlation is only qualitative due to the heterogeneous and highly dynamical environment. Based on density functional theory calculations, the dominant contribution for solvent-exposed -CO is found to arise from the surrounding water (∼75%), whereas the protein environment contributes considerably less. The results suggest that infrared spectroscopy is a positionally sensitive probe of insulin dimerization, in particular in conjunction with isotopic labeling of the probe.

Mapping protein selectivity landscapes using multi-target selective screening and next-generation sequencing of combinatorial libraries

Characterizing the binding selectivity landscape of interacting proteins is crucial both for elucidating the underlying mechanisms of their interaction and for developing selective inhibitors. However, current mapping methods are laborious and cannot provide a sufficiently comprehensive description of the landscape. Here, we introduce a novel and efficient strategy for comprehensively mapping the binding landscape of proteins using a combination of experimental multi-target selective library screening and in silico next-generation sequencing analysis. We map the binding landscape of a non-selective trypsin inhibitor, the amyloid protein precursor inhibitor (APPI), to each of the four human serine proteases (kallikrein-6, mesotrypsin, and anionic and cationic trypsins). We then use this map to dissect and improve the affinity and selectivity of APPI variants toward each of the four proteases. Our strategy can be used as a platform for the development of a new generation of target-selective probes and therapeutic agents based on selective protein-protein interactions.

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.

Equilibrium Ensembles for Insulin Folding from Bias-Exchange Metadynamics

Earliest events in the aggregation process, such as single molecule reconfiguration, are extremely important and the most difficult to characterize in experiments. To this end, we have used well-tempered bias exchange metadynamics simulations to determine the equilibrium ensembles of an insulin molecule under amyloidogenic conditions of low pH and high temperature. A bin-based clustering method that uses statistics accumulated in bias exchange metadynamics trajectories was employed to construct a detailed thermodynamic and kinetic model of insulin folding. The highest lifetime, lowest free-energy ensemble identified consisted of native conformations adopted by a folded insulin monomer in solution, namely, the R-, the R^f-, and the T-states of insulin. The lowest free-energy structure had a root mean square deviation of only 0.15 nm from native x-ray structure. The second longest-lived metastable state was an unfolded, compact monomer with little similarity to the native structure. We have identified three additional long-lived, metastable states from the bin-based model. We then carried out an exhaustive structural characterization of metastable states on the basis of tertiary contact maps and per-residue accessible surface areas. We have also determined the lowest free-energy path between two longest-lived metastable states and confirm earlier findings of non-two-state folding for insulin through a folding intermediate. The ensemble containing the monomeric intermediate retained 58% of native hydrophobic contacts, however, accompanied by a complete loss of native secondary structure. We have discussed the relative importance of nativelike versus nonnative tertiary contacts for the folding transition. We also provide a simple measure to determine the importance of an individual residue for folding transition. Finally, we have compared and contrasted this intermediate with experimental data obtained in spectroscopic, crystallographic, and calorimetric measurements during early stages of insulin aggregation. We have also determined stability of monomeric insulin by incubation at a very low concentration to isolate protein-protein interaction effects.

Molecular Physiology of Freeze Tolerance in Vertebrates

Freeze tolerance is an amazing winter survival strategy used by various amphibians and reptiles living in seasonally cold environments. These animals may spend weeks or months with up to ∼65% of their total body water frozen as extracellular ice and no physiological vital signs, and yet after thawing they return to normal life within a few hours. Two main principles of animal freeze tolerance have received much attention: the production of high concentrations of organic osmolytes (glucose, glycerol, urea among amphibians) that protect the intracellular environment, and the control of ice within the body (the first putative ice-binding protein in a frog was recently identified), but many other strategies of biochemical adaptation also contribute to freezing survival. Discussed herein are recent advances in our understanding of amphibian and reptile freeze tolerance with a focus on cell preservation strategies (chaperones, antioxidants, damage defense mechanisms), membrane transporters for water and cryoprotectants, energy metabolism, gene/protein adaptations, and the regulatory control of freeze-responsive hypometabolism at multiple levels (epigenetic regulation of DNA, microRNA action, cell signaling and transcription factor regulation, cell cycle control, and anti-apoptosis). All are providing a much more complete picture of life in the frozen state.

Effects of localized interactions and surface properties on stability of protein-based therapeutics

OBJECTIVES

Protein-based therapeutics garner significant attention because of exquisite specificity and limited side effects and are now being used to accomplish targeted delivery of small-molecule drugs. This review identifies and highlights individual chemical attributes and categorizes how site-specific changes affect protein stability based on published high-resolution molecular analyses.

KEY FINDINGS

Because it is challenging to determine the mechanisms by which the stability of large, complex molecules is altered and data are sparse, smaller, therapeutic proteins (insulin, erythropoietin, interferons) are examined alongside antibody data. Integrating this large pool of information with the limited available studies on antibodies reveals common mechanisms by which specific alterations affect protein structure and stability.

SUMMARY

Physical and chemical stability of therapeutic proteins and antibody drug conjugates (ADCs) is of critical importance because insufficient stability prevents molecules from making it to market. Individual moieties on/near the surface of proteins have substantial influence on structure and stability. Seemingly small, superficial modification may have far-reaching consequences on structure, conformational dynamics, and solubility of the protein, and hence physical stability of the molecule. Chemical modifications, whether spontaneous (e.g. oxidation, deamidation) or intentional, as with ADCs, may adversely impact stability by disrupting local surface properties or higher order protein structure.

Probing Receptor Specificity by Sampling the Conformational Space of the Insulin-like Growth Factor II C-domain

Insulin and insulin-like growth factors I and II are closely related protein hormones. Their distinct evolution has resulted in different yet overlapping biological functions with insulin becoming a key regulator of metabolism, whereas insulin-like growth factors (IGF)-I/II are major growth factors. Insulin and IGFs cross-bind with different affinities to closely related insulin receptor isoforms A and B (IR-A and IR-B) and insulin-like growth factor type I receptor (IGF-1R). Identification of structural determinants in IGFs and insulin that trigger their specific signaling pathways is of increasing importance in designing receptor-specific analogs with potential therapeutic applications. Here, we developed a straightforward protocol for production of recombinant IGF-II and prepared six IGF-II analogs with IGF-I-like mutations. All modified molecules exhibit significantly reduced affinity toward IR-A, particularly the analogs with a Pro-Gln insertion in the C-domain. Moreover, one of the analogs has enhanced binding affinity for IGF-1R due to a synergistic effect of the Pro-Gln insertion and S29N point mutation. Consequently, this analog has almost a 10-fold higher IGF-1R/IR-A binding specificity in comparison with native IGF-II. The established IGF-II purification protocol allowed for cost-effective isotope labeling required for a detailed NMR structural characterization of IGF-II analogs that revealed a link between the altered binding behavior of selected analogs and conformational rearrangement of their C-domains.

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 Methods and Protein Synthesis for Definition of Autoantibody Epitopes

Epitope mapping is the process of experimentally identifying the binding sites, or "epitopes," of antibodies on their target antigens. Understanding the antibody-epitope interaction provides a basis for the rational design of potential preventative vaccines. Islet autoantibodies are currently the best available biomarkers for predicting future type 1 diabetes. These include autoantibodies to the islet beta cell proteins, insulin and the tyrosine phosphatase islet antigen-2 (IA-2) which selectively bind to a small number of dominant epitopes associated with increased risk of disease progression. The major epitope regions of insulin and IA-2 autoantibodies have been identified, but need to be mapped more precisely. In order to characterize these epitopes more accurately, this article describes the methods of cloning and mutagenesis of insulin and IA-2 and subsequent purification of the proteins that can be tested in displacement analysis and used to monitor immune responses, in vivo, to native and mutated proteins in a humanized mouse model carrying the high-risk HLA class II susceptibility haplotype DRB1*04-DQ8.

Design, synthesis and biological evaluation of novel peptide MC2 analogues from Momordica charantia as potential anti-diabetic agents

Three series of Momordica charantia (MC)2 analogues were designed, synthesized and evaluated for their anti-hyperglycaemic effects. Alanine scanning focusing on the peptide MC2 indicated the importance of the residues proline (Pro)(3), serine (Ser)(6), isoleucine (Ile)(7) and Ser(10) for anti-hyperglycaemic effects. Among the first series of MC2 analogues, peptide I-4 exhibited a better anti-hyperglycaemic effect and was chosen for further modification. A further two series of conformationally constrained analogues were designed by scanning the residues Pro(3), Ser(6), Ile(7), and Ser(10) with an i - (i + 2) lactam bridge consisting of a glutamic acid-Xaa-lysine (Glu-Xaa-Lys) scaffold and a diproline fragment. By screening in normal mice and mice with diabetes mellitus, peptides II-1, II-2 and III-3 showed a significant improvement in anti-hyperglycaemic and anti-oxidative activities compared with I-4. These data suggest that II-1, II-2 and III-3 could be candidates for future treatment of diabetes mellitus.

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.