Assembly and dissociation of human insulin and LysB28ProB29-insulin hexamers: a comparison study

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

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

Kinetics of Phenol Escape from the Insulin R6 Hexamer

Therapeutic preparations of insulin often contain phenolic molecules, which can impact both pharmacokinetics and shelf life. Thus, understanding the interactions of insulin and phenolic molecules can aid in designing improved therapeutics. In this study, we use molecular dynamics to investigate phenol release from the insulin hexamer. Leveraging recent advances in methods for analyzing molecular dynamics data, we expand on existing simulation studies to identify and quantitatively characterize six phenol binding/unbinding pathways for wild-type and A10 Ile → Val and B13 Glu → Gln mutant insulins. A number of these pathways involve large-scale opening of the primary escape channel, suggesting that the hexamer is much more dynamic than previously appreciated. We show that phenol unbinding is a multipathway process, with no single pathway representing more than 50% of the reactive current and all pathways representing at least 10%. We use the mutant simulations to show how the contributions of specific pathways can be rationally manipulated. Predicting the net effects of mutations is more challenging because the kinetics depend on all of the pathways, demanding quantitatively accurate simulations and experiments.

Structural principles of insulin formulation and analog design: A century of innovation

BACKGROUND

The discovery of insulin in 1921 and its near-immediate clinical use initiated a century of innovation. Advances extended across a broad front, from the stabilization of animal insulin formulations to the frontiers of synthetic peptide chemistry, and in turn, from the advent of recombinant DNA manufacturing to structure-based protein analog design. In each case, a creative interplay was observed between pharmaceutical applications and then-emerging principles of protein science; indeed, translational objectives contributed to a growing molecular understanding of protein structure, aggregation and misfolding.

SCOPE OF REVIEW

Pioneering crystallographic analyses-beginning with Hodgkin's solving of the 2-Zn insulin hexamer-elucidated general features of protein self-assembly, including zinc coordination and the allosteric transmission of conformational change. Crystallization of insulin was exploited both as a step in manufacturing and as a means of obtaining protracted action. Forty years ago, the confluence of recombinant human insulin with techniques for site-directed mutagenesis initiated the present era of insulin analogs. Variant or modified insulins were developed that exhibit improved prandial or basal pharmacokinetic (PK) properties. Encouraged by clinical trials demonstrating the long-term importance of glycemic control, regimens based on such analogs sought to resemble daily patterns of endogenous β-cell secretion more closely, ideally with reduced risk of hypoglycemia.

MAJOR CONCLUSIONS

Next-generation insulin analog design seeks to explore new frontiers, including glucose-responsive insulins, organ-selective analogs and biased agonists tailored to address yet-unmet clinical needs. In the coming decade, we envision ever more powerful scientific synergies at the interface of structural biology, molecular physiology and therapeutics.

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.

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

Investigation of factors that cause insulin precipitation and/or amyloid formation in insulin formulations

Background

Multiple daily subcutaneous injections (MDSIs) are mainly used for formulating an insulin therapy for diabetic patients; however, they also cause insulin-derived amyloidosis (IDA) and lead to poor glycemic control. In addition, for the continuous subcutaneous insulin infusion system (CSII), precipitation frequently causes catheter occlusion and, if the precipitate in the formulations is amyloid, the injection of the insoluble amyloid into the subcutaneous tissue leads to IDA. The aim of this study was to conduct in vitro experiments and present a situation where insulin formulations cause precipitation and amyloid formation.

Methods

Humulin®R and NovoRapid® were used as model formulations for MDSIs and CSII, respectively. The generation of the precipitation was evaluated by measuring turbidity, and amyloid formation was evaluated by using Thioflavin T. Humulin®R was mixed with saline buffer solutions and glucose solutions to evaluate the effect of dilution. In addition, we created an experimental system to consider the effect of the time course of condition changes, and investigated the effects of insulin concentration, m-cresol existence, and pH change on the generation of the precipitate and amyloid in the formulation.

Results

In both the original and diluted formulations, physical stimulation resulted in the formation of a precipitate, which in most cases was an amyloid. The amyloid was likely to be formed at a near neutral pH. On the contrary, although a precipitate formed when the pH was decreased to near the isoelectric point, this precipitate was not an amyloid. Further decreases in pH resulted in the formation of amyloids, suggesting that both the positive and negative charged states of insulin tended to form amyloids. The formulation additive m-cresol suppressed amyloid formation. When additives were removed from the formulation, the amyloid-containing gel was formed in the field of substance exchange.

Conclusions

To consider changes in conditions that may occur for insulin formulations, the relationship between the formation of precipitates and amyloids was demonstrated in vitro by using insulin formulations. From the in vitro study, m-cresol was shown to have an inhibitory effect on amyloid formation.

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.

Insight into human insulin aggregation revisited using NMR derived translational diffusion parameters

The NMR derived translational diffusion coefficients were performed on unlabeled and uniformly labeled ¹³C,¹⁵N human insulin in water, both in neat, with zinc ions only, and in pharmaceutical formulation, containing only m-cresol as phenolic ligand, glycerol and zinc ions. The results show the dominant role of the pH parameter and the concentration on aggregation. The diffusion coefficient Dav was used for monitoring the overall average state of oligomeric ensemble in solution. The analysis of the experimental data of diffusion measurements, using the direct exponential curve resolution algorithm (DECRA) allows suggesting the two main components of the oligomeric ensemble. The 3D HSQC-iDOSY, (diffusion ordered HSQC) experiments performed on ¹³C, ¹⁵N-fully labeled insulin at the two pH values, 4 and 7.5, allow for the first time a more detailed experimental observation of individual components in the ensemble. The discussion involves earlier static and dynamic laser light scattering experiments and recent NMR derived translational diffusion results. The results bring new informations concerning the preparation of pharmaceutical formulation and in particular a role of Zn²⁺ ions. They also will enable better understanding and unifying the results of studies on insulin misfolding effects performed in solution by diverse physicochemical methods at different pH and concentration.

Rapid-Acting and Human Insulins: Hexamer Dissociation Kinetics upon Dilution of the Pharmaceutical Formulation

Purpose

Comparison of the dissociation kinetics of rapid-acting insulins lispro, aspart, glulisine and human insulin under physiologically relevant conditions.

Methods

Dissociation kinetics after dilution were monitored directly in terms of the average molecular mass using combined static and dynamic light scattering. Changes in tertiary structure were detected by near-UV circular dichroism.

Results

Glulisine forms compact hexamers in formulation even in the absence of Zn²⁺. Upon severe dilution, these rapidly dissociate into monomers in less than 10 s. In contrast, in formulations of lispro and aspart, the presence of Zn²⁺ and phenolic compounds is essential for formation of compact R6 hexamers. These slowly dissociate in times ranging from seconds to one hour depending on the concentration of phenolic additives. The disadvantage of the long dissociation times of lispro and aspart can be diminished by a rapid depletion of the concentration of phenolic additives independent of the insulin dilution. This is especially important in conditions similar to those after subcutaneous injection, where only minor dilution of the insulins occurs.

Conclusion

Knowledge of the diverging dissociation mechanisms of lispro and aspart compared to glulisine will be helpful for optimizing formulation conditions of rapid-acting insulins.

Electronic supplementary material

The online version of this article (doi:10.1007/s11095-017-2233-0) contains supplementary material, which is available to authorized users.

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

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

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.

Theoretical and computational studies of peptides and receptors of the insulin family

Synergistic interactions among peptides and receptors of the insulin family are required for glucose homeostasis, normal cellular growth and development, proliferation, differentiation and other metabolic processes. The peptides of the insulin family are disulfide-linked single or dual-chain proteins, while receptors are ligand-activated transmembrane glycoproteins of the receptor tyrosine kinase (RTK) superfamily. Binding of ligands to the extracellular domains of receptors is known to initiate signaling via activation of intracellular kinase domains. While the structure of insulin has been known since 1969, recent decades have seen remarkable progress on the structural biology of apo and liganded receptor fragments. Here, we review how this useful structural information (on ligands and receptors) has enabled large-scale atomically-resolved simulations to elucidate the conformational dynamics of these biomolecules. Particularly, applications of molecular dynamics (MD) and Monte Carlo (MC) simulation methods are discussed in various contexts, including studies of isolated ligands, apo-receptors, ligand/receptor complexes and intracellular kinase domains. The review concludes with a brief overview and future outlook for modeling and computational studies in this family of proteins.

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

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

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

Insulin provides a model for the therapeutic application of protein engineering. A paradigm in molecular pharmacology was defined by design of rapid-acting insulin analogs for the prandial control of glycemia. Such analogs, a cornerstone of current diabetes regimens, exhibit accelerated subcutaneous absorption due to more rapid disassembly of oligomeric species relative to wild-type insulin. This strategy is limited by a molecular trade-off between accelerated disassembly and enhanced susceptibility to degradation. Here, we demonstrate that this trade-off may be circumvented by nonstandard mutagenesis. Our studies employed Lys(B28), Pro(B29)-insulin ("lispro") as a model prandial analog that is less thermodynamically stable and more susceptible to fibrillation than is wild-type insulin. We have discovered that substitution of an invariant tyrosine adjoining the engineered sites in lispro (Tyr(B26)) by 3-iodo-Tyr (i) augments its thermodynamic stability (ΔΔGu 0.5 ± 0.2 kcal/mol), (ii) delays onset of fibrillation (lag time on gentle agitation at 37 °C was prolonged by 4-fold), (iii) enhances affinity for the insulin receptor (1.5 ± 0.1-fold), and (iv) preserves biological activity in a rat model of diabetes mellitus. (1)H NMR studies suggest that the bulky iodo-substituent packs within a nonpolar interchain crevice. Remarkably, the 3-iodo-Tyr(B26) modification stabilizes an oligomeric form of insulin pertinent to pharmaceutical formulation (the R6 zinc hexamer) but preserves rapid disassembly of the oligomeric form pertinent to subcutaneous absorption (T6 hexamer). By exploiting this allosteric switch, 3-iodo-Tyr(B26)-lispro thus illustrates how a nonstandard amino acid substitution can mitigate the unfavorable biophysical properties of an engineered protein while retaining its advantages.

Effects of phenol and meta-cresol depletion on insulin analog stability at physiological temperature

The stability of three commercial "fast-acting" insulin analogs, insulin lispro, insulin aspart, and insulin glulisine, was studied at various concentrations of phenolic preservatives (phenol and/or meta-cresol) during 9 days of incubation at 37 °C. The analysis by both size-exclusion and reversed-phase chromatography showed degradation of lispro and aspart that was inversely dependent on the concentration of phenolic preservatives. Insulin glulisine was much more stable than the other analogs and showed minimal degradation even in the absence of phenolic preservatives. With sedimentation velocity ultracentrifugation, we determined the preservatives' effect on the insulins' self-assembly. When depleted of preservatives, insulin glulisine dissociates from higher molecular weight species into a number of intermediate molecular weight species, in between monomer and hexamer, whereas insulin aspart and insulin lispro dissociate into monomers and dimers. Decreased stability of insulin lispro and insulin aspart seems to be because of the extent of dissociation when depleted of preservative. Insulin glulisine's dissociation to intermediate molecular weight species appears to help minimize its degradation during incubation at 37 °C.

Thermodynamic contributions to the stability of the insulin hexamer

The insulin hexamer is resistant to degradation and fibrillation, which makes it an important quaternary structure for its in vivo storage in Zn(2+)- and Ca(2+)-rich vesicles in the pancreas and for pharmaceutical formulations. In addition to the two Zn(2+) ions that are required for its formation, three other species, Zn-coordinating anions (e.g., Cl(-)), Ca(2+), and phenols (e.g., resorcinol), bind to the hexamer and affect the subunit conformation and stability. The contributions of these four species to the thermodynamics of insulin unfolding have been quantified by differential scanning calorimetry and thermal unfolding measurements to determine the extent and nature of their stabilization of the insulin hexamer. Both Zn(2+) and resorcinol make a significant enthalpic contribution, while Ca(2+) primarily affects the protein heat capacity (solvation) by its interactions in the central cation-binding cavity, which is modulated by the surrounding subunit conformations. Coordinating anions have a negligible effect on the stability of the hexamer, even though subunits shift to an alternate conformation when these anions bind to the Zn(2+) ions. Finally, Zn(2+) in excess of the two that are required to form the hexamer further stabilizes the protein by additional enthalpic contributions.

Insulin aspart pharmacokinetics: an assessment of its variability and underlying mechanisms

BACKGROUND

Insulin aspart (IAsp) is used by many diabetics as a meal-time insulin to control post-prandial glucose levels. As is the case with many other insulin types, the pharmacokinetics (PK), and consequently the pharmacodynamics (PD), is associated with clinical variability, both between and within individuals. The present article identifies the main physiological mechanisms that govern the PK of IAsp following subcutaneous administration and quantifies them in terms of their contribution to the overall variability.

MATERIAL AND METHODS

CT scanning data from Thomsen et al. (2012) are used to investigate and quantify the properties of the subcutaneous depot. Data from Brange et al. (1990) are used to determine the effects of insulin chemistry in subcutis on the absorption rate. Intravenous (i.v.) bolus and infusion PK data for human insulin are used to understand and quantify the systemic distribution and elimination (Pørksen et al., 1997; Sjöstrand et al., 2002). PK and PD profiles for type 1 diabetics from Chen et al. (2005) are analyzed to demonstrate the effects of IAsp antibodies in terms of bound and unbound insulin. PK profiles from Thorisdottir et al. (2009) and Ma et al. (2012b) are analyzed in the nonlinear mixed effects software Monolix® to determine the presence and effects of the mechanisms described in this article.

RESULTS

The distribution of IAsp in the subcutaneous depot show an initial dilution of approximately a factor of two in a single experiment. Injected insulin hexamers exist in a chemical equilibrium with monomers and dimers, which depends strongly on the degree of dilution in subcutis, the presence of auxiliary substances, and a variety of other factors. Sensitivity to the initial dilution in subcutis can thus be a cause of some of the variability. Temporal variations in the PK are explained by variations in the subcutaneous blood flow. IAsp antibodies are found to be a large contributor to the variability of total insulin PK in a study by Chen et al. (2005), since only the free fraction is eliminated via the receptors. The contribution of these and other sources of variability to the total variability is quantified via a population PK analysis and two recent clinical studies (Thorisdottir et al., 2009; Ma et al., 2012b), which support the presence and significance of the identified mechanisms.

CONCLUSIONS

IAsp antibody binding, oligomeric transitions in subcutis, and blood flow dependent variations in absorption rate seem to dominate the PK variability of IAsp. It may be possible via e.g. formulation design to reduce some of these variability factors.

Phenomena of insulin peak fronting in size exclusion chromatography and strategies to reduce fronting

Insulin peak fronting in size exclusion chromatography (SEC) results in more than 10% yield loss in the production of insulin. The goal of this study is to understand the mechanisms of peak fronting and to develop strategies to reduce fronting and increase insulin yield. Chromatography experiments ruled out pressure surge, viscous fingering, and adsorption as the cause for peak fronting. Theoretical analysis based on a general rate model indicated that reversible dimerization is the major cause for peak fronting and reducing the dimerization equilibrium constant is the most effective method for reducing fronting. Two strategies were developed and tested to reduce the degree of dimer formation. The first strategy was to use 0.1N acetic acid as the presaturant and eluent. The second strategy was to use 0.8 or 2.8N acetic acid in 20vol.% denatured ethanol as the mobile phase. The experimental results showed that both strategies can reduce insulin peak fronting in SEC, maintain desired product purity, and increase insulin yield to higher than 98%.

Insulin Self-association: Effects on Lung Disposition Kinetics in the Airways of the Isolated Perfused Rat Lung (IPRL)

PURPOSE

To characterize the kinetic dependence of pulmonary absorption and metabolism of insulin and lispro on the magnitude of their hexameric association.

METHODS

Hexamer content by weight percent (%Hex) in various insulin-zinc and lispro-zinc solutions were determined by quantitative centrifugal ultrafiltration and zinc titration with terpyridine (QCUF-ZTT). Each of the solutions (0.1 ml) was then administered into the airways of the IPRL of normal and experimental diabetic animals. Rate constants were determined for lung absorption (k (a)) and non-absorptive loss (k (nal); comprising mucociliary clearance and metabolism).

RESULTS

%Hex in administered solutions ranged from 3.3 to 94.4%. Data analysis showed excellent correlations between the values for k (a) or k (nal) and %Hex, irrespective of insulin type, concentration, solution pH or ionic strength. The values for k (a) decreased (0.22 --> 0.05 h(-1)) with increasing %Hex, as did values for k (nal). At %Hex in administered solutions >/=50%, values for k (nal) approached estimates for the rate constant for mucociliary clearance, implying that lung metabolism occurred primarily with monomeric insulin. There were no differences in insulin disposition kinetics between lungs taken from experimental diabetic and sham-control animals.

CONCLUSIONS

The kinetics of pulmonary insulin disposition depended on the magnitude of molecular self-association. Dissociated forms of insulin (dimers or monomers) in the dosing solution showed higher rates than hexamers for both lung absorption and metabolism.

Insulin glulisine complementing basal insulins: a review of structure and activity

The advancement in protein engineering offers targeted development of insulin analogs that display either faster absorption kinetics or longer time-action profiles compared with human insulin and, therefore, more closely mimic endogenous insulin secretion. Insulin glulisine (3(B)Lys29(B) Glu-human insulin) is a new fast-acting analog that provides absorption and onset of action more rapidly with a shorter duration of action compared with regular human insulin, and thus better resembles physiologic mealtime insulin requirements. Insulin glulisine has been designed to exhibit intrinsic stability while maintaining rapid deployment of insulin monomers. Pharmacokinetic and pharmacodynamic profiling of insulin glulisine in healthy subjects and patients with type 1 and type 2 diabetes not only confirms the rapid absorption and fast action of insulin glulisine compared with human insulin, but also provides evidence that the unique drug formulation may offer additional benefits. Insulin glulisine complements insulin glargine (21(A)-Gly30(Ba)-L-Arg-30(Bb)-L-Arg-human insulin), the first long-acting basal insulin analog that displays a smoothed time-action profile with a 24-h duration of action. Together these analogs offer patients a more physiologic approach to insulin replacement.

Theoretical analysis of the effects of reversible dimerization in size exclusion chromatography

Reversible dimer formation in size exclusion chromatography (SEC) can cause peak splitting, merging, tailing, and fronting. Such behavior can be predicted by the association rate and the dissociation rate relative to the convection rate. Slow association and dissociation result in separated monomer and dimer peaks. Fast association and slow dissociation result in one single dimer peak. Slow association and fast dissociation result in one single monomer peak. Intermediate association and dissociation result in a merged, broad peak with either fronting when monomers dominate or tailing when dimers dominate. A diagram based on the two relative rates is generated to predict general peak shape and retention behavior in SEC.

Molecular variants and derivatives of insulin for improved glycemic control in diabetes

Insulin is a historic molecule. It presents many first instances, such as the first protein to be fully sequenced, one of the first proteins to be crystallized in pure form, one among the early proteins whose structure was investigated using X-ray crystallography, the first protein to be chemically synthesized and the first Biotech drug. Therefore, the development of insulin in the early years is intricately intertwined with the progress in molecular and structural biology. In recent years, development of a range of insulin analogs has led to better control of glucose levels, thus preventing secondary complications and improving the quality of life in diabetic patients. Such analogs were obtained by modification of the native insulin sequence. They vary with regard to their pharmacokinetic profile, stability, tissue specificity and mode of administration. In addition, alterations involving incorporation of various chemical moieties in insulin and its co-crystallization with insoluble derivatives are used to modulate the time-action profile of the drug. This article traces the development of molecular variants and derivatives of insulin. It discusses future directions for further improvement in their properties to produce still better insulin therapeutics for tight glycemic control.

In situ probing of insulin aggregation in chromatography effluents with spectroturbidimetry

Probing protein aggregation in situ is quite important for analyzing and developing chromatographic protein purification processes. A spectroturbidimetry method with a photodiode array detector is developed and tested for probing insulin aggregation in solution and determining the aggregation number, n(m). All aggregates examined are in the Rayleigh light scattering regime, where the turbidity between 400 and 350 nm is proportional to lambda(-4). Insulin at 25 degrees C in 3.5 N acetic acid is mainly monomeric (non-aggregated). At 25 degrees C and lower acetic acid concentrations, from 0.1 to 1 N, the average insulin aggregation number n(m) ranges from 2.9 to 1.6. Aggregates, with n(m) = 2-3, are found in 2.6 N acetic acid with 20 vol% acetonitrile. In 0.8 N acetic acid with 20 vol% denatured ethanol, n(m) = 1.2. At 4 degrees C, as acetic acid concentration decreases from 3.5 to 0.1 N, n(m) decreases from 2.4 to 1.8. In 2.8 N acetic acid with 20 vol% denatured ethanol at 4 degrees C, insulin exists mainly in monomer form. In situ probing of size exclusion chromatography, SEC, effluents in 3.5 N acetic acid at 4 degrees C shows n(m) = 1.6 at the fronting portion (a mixture of monomers and dimers or other oligomers) and n(m) = 1.1 (mostly monomers) at the tailing portion of the main peak. In another example, for LysPro-insulin in reversed phase chromatography at 4 degrees C, complex elution patterns and broad peaks are due to substantial aggregation. For a linear gradient of acetonitrile from 10 to 60 vol% at 4 degrees C, n(m) ranges from 2.2 to 12, in order of elution. For a linear gradient of ethanol from 30 to 50 vol% at 4 degrees C, n(m) ranges from 14 to 27, in order of elution. Analytical HPLC results at 25 degrees C imply that the aggregates are reversible.

Ligand Binding and Thermostability of Different Allosteric States of the Insulin Zinc−Hexamer

The influence of ligand binding and conformation state on the thermostability of hexameric zinc-insulin was studied by differential scanning calorimetry (DSC). The insulin hexamer exists in equilibrium between the forms T6, T3R3, and R6. Phenolic ligands induce and stabilize the T3R3- and R6-states which are further stabilized by binding of certain anions that do not stabilize the T6-state. It was shown that the thermostability of the resorcinol-stabilized R6-state was significantly higher than that of the T6-state. Further analysis showed that phenol- and m-cresol-stabilized R6-hexamer loses three ligands before reaching the unfolding temperature and hence unfolds from the T3R3-state. The relative affinity of the four tested anionic ligands was found, by DSC, to be thiocyanate > or = 4-hydroxy-3-nitrobenzoate >> p-aminobenzoate >> chloride. The results correlate with other methods and demonstrate that DSC provides a general and useful method of evaluation of both phenolic and anionic ligand binding to insulin without the use of probes or other alterations of the system of interest. However, it is a prerequisite that the binding is strong enough to saturate the binding sites at temperatures around the unfolding transition.

MD simulation of protein-ligand interaction: formation and dissociation of an insulin-phenol complex

Complexes of proteins with small ligands are of utmost importance in biochemistry, and therefore equilibria, formation, and decay have been investigated extensively by means of biochemical and biophysical methods. Theoretical studies of the molecular dynamics of such systems in solution are restricted to 10 ns, i.e., to fast processes. Only recently new theoretical methods have been developed not to observe the process in real time, but to explore its pathway(s) through the energy landscape. From the profiles of free energy, equilibrium and kinetic quantities can be determined using transition-state theory. This study is dedicated to the pharmacologically relevant insulin-phenol complex. The distance of the center of mass chosen as a reaction coordinate allows a reasonable description over most of the pathway. The analysis is facilitated by analytical expressions we recently derived for distance-type reaction coordinates. Only the sudden onset of rotations at the very release of the ligand cannot be parameterized by a distance. They obviously require a particular treatment. Like a preliminary study on a peptide, the present case emphasizes the contribution of internal friction inside a protein, which can be computed from simulation data. The calculated equilibrium constant and the friction-corrected rates agree well with experimental data.

Structural studies of a crystalline insulin analog complex with protamine by atomic force microscopy

Crystallographic studies of insulin-protamine complexes, such as neutral protamine Hagedorn (NPH) insulin, have been hampered by high crystal solvent content, small crystal dimensions, and extensive disorder in the protamine molecules. We report herein in situ tapping mode atomic force microscopy (TMAFM) studies of crystalline neutral protamine Lys(B28)Pro(B29) (NPL), a complex of Lys(B28)Pro(B29) insulin, in which the C-terminal prolyl and lysyl residues of human insulin are inverted, and protamine that is used as an intermediate time-action therapy for treating insulin-dependent diabetes. Tapping mode AFM performed at 6 degrees C on bipyramidally tipped tetragonal rod-shaped NPL crystals revealed large micron-sized islands separated by 44-A tall steps. Lattice images obtained by in situ TMAFM phase and height imaging on these islands were consistent with the arrangement of individual insulin-protamine complexes on the P4(1)2(1)2 (110) crystal plane of NPH, based on a low-resolution x-ray diffraction structure of NPH, arguing that the NPH and NPL insulins are isostructural. Superposition of the height and phase images indicated that tip-sample adhesion was larger in the interstices between NPL complexes in the (110) crystal plane than over the individual complexes. These results demonstrate the utility of low-temperature TMAFM height and phase imaging for the structural characterization of biomolecular complexes.

Pulmonary insulin administration using the AERx system: physiological and physicochemical factors influencing insulin effectiveness in healthy fasting subjects

BACKGROUND

Orally inhaled insulin may provide a convenient and effective therapy for prandial glucose control in patients with diabetes. This study evaluated the influence of formulation pH and concentration and different respiratory maneuvers on pharmacokinetic and pharmacodynamic properties of inhaled insulin.

METHODS

Three, open-label crossover studies in a total of 23 healthy subjects were conducted in which the safety, pharmacokinetics, and pharmacodynamics of insulin inhalation were compared to subcutaneous (SC) injection into the abdomen of commercially available regular insulin. A novel, aerosol generating system (AERx Diabetes Management System, Aradigm Corporation, Hayward, CA) was used to deliver aqueous insulin bolus aerosols to the lower respiratory tract from formulations at pH 3.5 or 7.4 and concentrations of U250 (250 U/mL) or U500 (500 U/mL).

RESULTS

Time to maximum insulin concentration in serum (Tmax) after SC dosing occurred approximately 50-60 minutes with the time to minimum plasma glucose concentration (i.e., maximum hypoglycemic effect), (TGmin), occurring later, at around 100-120 minutes. In contrast, pulmonary delivery led to a significantly earlier Tmax (7-20 minutes) and TGmin (60-70 minutes), parameters that were shown to be largely unaffected by changing the pH or concentration of the insulin. However, investigation of changes in inhaled volume (achieved by different programming of the AERx system) for administration of the same sized aerosol bolus revealed significant effects. Significantly slower absorption and time to peak hypoglycemic activity occurred when aerosol delivery of insulin occurred during a shallow (approximately 40% vital capacity) as opposed to a deep (approximately 80% vital capacity) inspiration. In addition, it was shown that serum concentration of insulin increased immediately after a series of forced expiraratory maneuvers 30 minutes after inhaled delivery.

CONCLUSIONS

Pulmonary delivery of aqueous bolus aerosols of insulin in healthy subjects resulted in rapid absorption with an associated hypoglycemic effect quicker than is achieved after subcutaneous dosing of regular insulin. Inhaled insulin pharmacokinetics and pharmacodynamics were independent of formulation variables (pH, concentration) but affected by certain respiratory maneuvers.

Preparation of a microcrystalline suspension formulation of Lys(B28)Pro(B29)-human insulin with ultralente properties

The monomeric analogue, Lys(B28)Pro(B29)-human insulin (LysPro), has been crystallized using similar conditions employed to prepare extended-acting insulin ultralente formulations. In the presence of zinc ions, sodium acetate and sodium chloride, but without phenolic preservative, LysPro surprisingly forms small rhombohedral crystals with similar morphology to human insulin ultralente crystals with a mean particle size of 20 +/- 1 microm. X-ray powder diffraction studies on the LysPro crystals prior to dilution in ultralente vehicle ([NaCl] = 1.2 M) revealed the presence of T(3)R(3)(f) hexamers. Consistent with human insulin ultralente preparations, LysPro crystals formulated as an ultralente suspension ([NaCl] = 0. 12 M) contain T(6) hexamers indicating that a conformational change occurs in the hexamer units of the crystals upon dilution of the salt concentration. The pharmacological properties of subcutaneously administered ultralente LysPro (ULP) were compared to ultralente human insulin (UHI) using a conscious dog model (n = 5) with glucose levels clamped at basal. There were no statistically significant differences between the kinetic and dynamic responses of ULP compared to UHI [C(max) (ng/mL): 3.58 +/- 0.76, ULP and 3.61 +/- 0. 66, UHI; T(max) (min): 226 +/- 30, ULP and 185 +/- 42, UHI; R(max) (mg/kg min): 11.2 +/- 1.9, ULP and 13.3 +/- 2.0, UHI; and T(Rmax) (min): 336 +/- 11, ULP and 285 +/- 57, UHI]. Although the Pro to Lys sequence inversion destabilizes insulin self-assembly and greatly alters the time action of soluble LysPro preparations, this modification has now been found neither to prevent the formation of ultralente crystals in the absence of phenolics nor to compromise the protracted activity of the insulin analogue suspension.

Insulin analogs with improved pharmacokinetic profiles

The aim of insulin replacement therapy is to normalize blood glucose in order to reduce the complications of diabetes. The pharmacokinetics of the traditional insulin preparations, however, do not match the profiles of physiological insulin secretion. The introduction of the rDNA technology 20 years ago opened new ways to create insulin analogs with altered properties. Fast-acting analogs are based on the idea that an insulin with less tendency to self-association than human insulin would be more readily absorbed into the systemic circulation. Protracted-acting analogs have been created to mimic the slow, steady rate of insulin secretion in the fasting state. The present paper provides a historical review of the efforts to change the physicochemical and pharmacological properties of insulin in order to improve insulin therapy. The available clinical studies of the new insulins are surveyed and show, together with modeling results, that new strategies for optimal basal-bolus treatment are required for utilization of the new fast-acting analogs.

Site-specific insulin conjugates with enhanced stability and extended action profile

Two different hydrophilic moieties, carboxyl derivatives of monosaccharidic (Glc, Gal, Man, Fuc) glycosides and methoxypolyethylene glycols of varying MW, were covalently attached to the insulin GlyA1, PheB1 and/or LysB29 amino groups (seven possible derivatives), and resulting insulin conjugates purified to homogeneity. In vivo bioactivity in rats of most derivatives was preserved while disubstituted PEG-insulins showed decreased potency. Only site-specific modification of PheB1 amino group with either moiety resulted in pronouncedly increased resistance of insulin to fibrillation, indicating that the B-chain N-terminus of the insulin molecule is mechanistically involved in the fibrillation process. Immunogenicity in vivo and in vitro of monoglycosylated insulins was comparable to that of insulin, diglycosylated insulins showed immunogenicity enhancement. Immunogenicity of PEG-insulins was significantly suppressed. PheB1-glycosylated insulins administered subcutaneously in dogs displayed extended action profiles, the most effective being PheB1-galactosylated insulin, resembling the pharmacodynamic response of intermediate-acting insulin preparations. The pharmacokinetic parameters of these insulin derivatives were not significantly different from that of insulin even though absorption and residence time and clearance were increased, providing some explanation for prolonged action profile. Lectin-specific binding as a retardation basis is not likely involved. In support of this, subcutaneously administered PheB1-PEG(600)-insulin showed an even more protracted action profile, suggesting that the basis of retardation is physical and nonspecific. This implies that by increasing PEG chain MW, further delay/prolongation of action can be achieved to yield new soluble basal insulin substitutes with potential clinical applications.

Self-association properties of monomeric insulin analogs under formulation conditions

PURPOSE

The purpose of the current study was to investigate the effects of two important excipients, zinc and m-cresol, on the self-association properties of a series of monomeric insulin analogs. In this way, the effects on formulation behavior of individual amino acid substitutions in the C-terminal region of the insulin B-chain could be compared.

METHODS

The self-association of ten insulin analogs was monitored by equilibrium and velocity analytical ultracentrifugation under three different conditions: (i) in neutral buffer alone; (ii) in neutral buffer containing zinc ion; and (iii) in neutral buffer containing both zinc ion and phenolic preservative (a typical condition for insulin formulations). The self-association properties of these analogs were compared to those of human insulin and the rapid-acting insulin analog Lys(B28)Pro(B29)-human insulin.

RESULTS

The analogs in the current study exhibited a wide range of association properties when examined in neutral buffer alone or in neutral buffer containing zinc ion. However, all of these analogs had association properties similar to human insulin in the presence of both zinc and m-cresol. Under these formulation conditions each analog had an apparent sedimentation coefficient of s* = 2.9-3.1 S, which corresponds to the insulin hexamer.

CONCLUSIONS

Analogs with changes in the B27-B29 region of human insulin form soluble hexamers in the presence of both zinc and m-cresol, and m-cresol binding overrides the otherwise destabilizing effects of these mutations on self assembly.

Atomic force microscopy of crystalline insulins: the influence of sequence variation on crystallization and interfacial structure

The self-association of proteins is influenced by amino acid sequence, molecular conformation, and the presence of molecular additives. In the presence of phenolic additives, LysB28ProB29 insulin, in which the C-terminal prolyl and lysyl residues of wild-type human insulin have been inverted, can be crystallized into forms resembling those of wild-type insulins in which the protein exists as zinc-complexed hexamers organized into well-defined layers. We describe herein tapping-mode atomic force microscopy (TMAFM) studies of single crystals of rhombohedral (R3) LysB28ProB29 that reveal the influence of sequence variation on hexamer-hexamer association at the surface of actively growing crystals. Molecular scale lattice images of these crystals were acquired in situ under growth conditions, enabling simultaneous identification of the rhombohedral LysB28ProB29 crystal form, its orientation, and its dynamic growth characteristics. The ability to obtain crystallographic parameters on multiple crystal faces with TMAFM confirmed that bovine and porcine insulins grown under these conditions crystallized into the same space group as LysB28ProB29 (R3), enabling direct comparison of crystal growth behavior and the influence of sequence variation. Real-time TMAFM revealed hexamer vacancies on the (001) terraces of LysB28ProB29, and more rounded dislocation noses and larger terrace widths for actively growing screw dislocations compared to wild-type bovine and porcine insulin crystals under identical conditions. This behavior is consistent with weaker interhexamer attachment energies for LysB28ProB29 at active growth sites. Comparison of the single crystal x-ray structures of wild-type insulins and LysB28ProB29 suggests that differences in protein conformation at the hexamer-hexamer interface and accompanying changes in interhexamer bonding are responsible for this behavior. These studies demonstrate that subtle changes in molecular conformation due to a single sequence inversion in a region critical for insulin self-association can have a significant effect on the crystallization of proteins.

Preparation and characterization of a cocrystalline suspension of [LysB28,ProB29]-human insulin analogue

Soluble preparations of [LysB28,ProB29]-human insulin analogue (LysPro) exhibit more rapid absorption than human insulin upon subcutaneous injection. Biphasic mixtures of LysPro and intermediate-acting insulin suspensions could provide advantages over current preparations for the treatment of diabetes. To prepare biphasic mixtures of LysPro, a suspension formulation of the analogue is required. We have devised a method for crystallizing LysPro with the basic peptide protamine yielding neutral protamine LysPro (NPL) suspension. The crystallization conditions are strongly dependent on the precipitation procedure and temperature. Using various techniques, the crystalline and suspension characteristics of NPL are found to be similar to human insulin (neutral protamine Hagedorn, NPH) (8:1 molar ratio insulin:protamine, rod-shaped crystals, particle size of 4.0-6.0 microns, and Point of Zero Charge at 6.0-7.0). Using a dog model with NPL or NPH injected subcutaneously and glucose levels clamped at basal, NPL was found to have kinetic and dynamic responses analogous to human insulin NPH [Cmax (maximal insulin or LysPro concentration, ng/mL) of 2.61 +/- 0.22, NPL; 2.58 +/- 0.36, NPH, attained at Tmax (min) of 93 +/- 22, NPL; 145 +/- 33 NPH, and Rmax (maximal rate of glucose infusion, mg/kg min) of 10.8 +/- 1.2, NPL; 13.2 +/- 1.9, NPH, attained at TRmax (min) of 277 +/- 58, NPL; 265 +/- 38, NPH]. There are no statistically significant differences between the insulin curves or the glucose responses. These results provide insight into the mechanism of action of NPH suspensions and the relationship to duration of action. Furthermore, the formulation of a suspension of LysPro having an intermediate time-action makes possible the preparation of stable biphasic mixtures containing LysPro and NPL.