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

Structural models of viral insulin-like peptides and their analogs

The insulin receptor (IR), the insulin-like growth factor-1 receptor (IGF1R), and the insulin/IGF1 hybrid receptors (hybR) are homologous transmembrane receptors. The peptide ligands, insulin and IGF1, exhibit significant structural homology and can bind to each receptor via site-1 and site-2 residues with distinct affinities. The variants of the Iridoviridae virus family show capability in expressing single-chain insulin/IGF1 like proteins, termed viral insulin-like peptides (VILPs), which can stimulate receptors from the insulin family. The sequences of VILPs lacking the central C-domain (dcVILPs) are known, but their structures in unbound and receptor-bound states have not been resolved to date. We report all-atom structural models of three dcVILPs (dcGIV, dcSGIV, and dcLCDV1) and their complexes with the receptors (μIR, μIGF1R, and μhybR), and probed the peptide/receptor interactions in each system using all-atom molecular dynamics (MD) simulations. Based on the nonbonded interaction energies computed between each residue of peptides (insulin and dcVILPs) and the receptors, we provide details on residues establishing significant interactions. The observed site-1 insulin/μIR interactions are consistent with previous experimental studies, and a residue-level comparison of interactions of peptides (insulin and dcVILPs) with the receptors revealed that, due to sequence differences, dcVILPs also establish some interactions distinct from those between insulin and IR. We also designed insulin analogs and report enhanced interactions between some analogs and the receptors.

Structures and interactions of insulin-like peptides from cone snail venom

The venomous insulin-like peptides released by certain cone snails stimulate hypoglycemic shock to immobilize fish and catch the prey. Compared to human insulin (hIns), the cone snail insulins (Con-Ins) are typically monomeric and shorter in sequence, yet they exhibit moderate hIns-like biological activity. We have modeled six variants of Con-Ins (G3, K1, K2, T1A, T1B, and T2) and carried out explicit-solvent molecular dynamics (MD) simulations of eight types of insulins, two with known structures (hIns and Con-Ins-G1) and six Con-Ins with modeled structures, to characterize key residues of each insulin that interact with the truncated human insulin receptor (μIR). We show that each insulin/μIR complex is stable during explicit-solvent MD simulations and hIns interactions indicate the highest affinity for the "site 1" of IR. The residue contact maps reveal that each insulin preferably interacts with the αCT peptide than the L1 domain of IR. Through analysis of the average nonbonded interaction energy contribution of every residue of each insulin for the μIR, we probe the residues establishing favorable interactions with the receptor. We compared the interaction energy of each residue of every Con-Ins to the μIR and observed that γ-carboxylated glutamate (Gla), His, Thr, Tyr, Tyr/His, and Asn in Con-Ins are favorable substitutions for GluA4, AsnA21, ValB12, LeuB15, GlyB20, and ArgB22 in hIns, respectively. The identified insulin analogs, although lacking the last eight residues of the B-chain of hIns, bind strongly to μIR. Our findings are potentially useful in designing potent fast-acting therapeutic insulin.

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.

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.

All-Atom Structural Models of the Transmembrane Domains of Insulin and Type 1 Insulin-Like Growth Factor Receptors

The receptor tyrosine kinase superfamily comprises many cell-surface receptors including the insulin receptor (IR) and type 1 insulin-like growth factor receptor (IGF1R) that are constitutively homodimeric transmembrane glycoproteins. Therefore, these receptors require ligand-triggered domain rearrangements rather than receptor dimerization for activation. Specifically, binding of peptide ligands to receptor ectodomains transduces signals across the transmembrane domains for trans-autophosphorylation in cytoplasmic kinase domains. The molecular details of these processes are poorly understood in part due to the absence of structures of full-length receptors. Using MD simulations and enhanced conformational sampling algorithms, we present all-atom structural models of peptides containing 51 residues from the transmembrane and juxtamembrane regions of IR and IGF1R. In our models, the transmembrane regions of both receptors adopt helical conformations with kinks at Pro961 (IR) and Pro941 (IGF1R), but the C-terminal residues corresponding to the juxtamembrane region of each receptor adopt unfolded and flexible conformations in IR as opposed to a helix in IGF1R. We also observe that the N-terminal residues in IR form a kinked-helix sitting at the membrane-solvent interface, while homologous residues in IGF1R are unfolded and flexible. These conformational differences result in a larger tilt-angle of the membrane-embedded helix in IGF1R in comparison to IR to compensate for interactions with water molecules at the membrane-solvent interfaces. Our metastable/stable states for the transmembrane domain of IR, observed in a lipid bilayer, are consistent with a known NMR structure of this domain determined in detergent micelles, and similar states in IGF1R are consistent with a previously reported model of the dimerized transmembrane domains of IGF1R. Our all-atom structural models suggest potentially unique structural organization of kinase domains in each receptor.