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

A λ-Dynamics Investigation of Insulin Wakayama and Other A3 Variant Binding Affinities to the Insulin Receptor

Insulin Wakayama is a clinical insulin variant where a conserved valine at the third residue on insulin's A chain (ValA3) is replaced with a leucine (LeuA3), weakening insulin receptor (IR) binding by 140-500-fold. This severe impact on binding from a subtle modification has posed an intriguing problem for decades. Although experimental investigations of natural and unnatural A3 mutations have highlighted the sensitivity of insulin-IR binding at this site, atomistic explanations of these binding trends have remained elusive. We investigate this problem computationally using λ-dynamics free energy calculations to model structural changes in response to perturbations of the ValA3 side chain and to calculate associated relative changes in binding free energy (ΔΔGbind). The Wakayama LeuA3 mutation and seven other A3 substitutions were studied in this work. The calculated ΔΔGbind results showed high agreement compared to experimental binding potencies with a Pearson correlation of 0.88 and a mean unsigned error of 0.68 kcal/mol. Extensive structural analyses of λ-dynamics trajectories revealed that critical interactions were disrupted between insulin and the insulin receptor as a result of the A3 mutations. This investigation also quantifies the effect that adding an A3 Cδ atom or losing an A3 Cγ atom has on insulin's binding affinity to the IR. Thus, λ-dynamics was able to successfully model the effects of mutations to insulin's A3 side chain on its protein-protein interactions with the IR and shed new light on a decades-old mystery: the exquisite sensitivity of hormone-receptor binding to a subtle modification of an invariant insulin residue.

Single-chain insulin analogs threaded by the insulin receptor αCT domain

Insulin is a mainstay of therapy for diabetes mellitus, yet its thermal stability complicates global transportation and storage. Cold-chain transport, coupled with optimized formulation and materials, prevents to some degree nucleation of amyloid and hence inactivation of hormonal activity. These issues hence motivate the design of analogs with increased stability, with a promising approach being single-chain insulins (SCIs), whose C domains (foreshortened relative to proinsulin) resemble those of the single-chain growth factors (IGFs). We have previously demonstrated that optimized SCIs can exhibit native-like hormonal activity with enhanced thermal stability and marked resistance to fibrillation. Here, we describe the crystal structure of an ultrastable SCI (C-domain length 6; sequence EEGPRR) bound to modules of the insulin receptor (IR) ectodomain (N-terminal α-subunit domains L1-CR and C-terminal αCT peptide; "microreceptor" [μIR]). The structure of the SCI-μIR complex, stabilized by an Fv module, was determined using diffraction data to a resolution of 2.6 Å. Remarkably, the αCT peptide (IR-A isoform) "threads" through a gap between the flexible C domain and the insulin core. To explore such threading, we undertook molecular dynamics simulations to 1) compare threaded with unthreaded binding modes and 2) evaluate effects of C-domain length on these alternate modes. The simulations (employing both conventional and enhanced sampling simulations) provide evidence that very short linkers (C-domain length of -1) would limit gap opening in the SCI and so impair threading. We envisage that analogous threading occurs in the intact SCI-IR complex-rationalizing why minimal C-domain lengths block complete activity-and might be exploited to design novel receptor-isoform-specific analogs.

Progress in Simulation Studies of Insulin Structure and Function

Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.