Comparative studies on the dynamics of crosslinked 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.

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

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

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.

Growth factor receptors: structure, mechanism, and drug discovery

The focus of this review is the relationship between the three-dimensional structure of ligands of the various members of the growth factor receptor tyrosine kinase superfamily and their interaction with the cognate receptor. Particular attention is given to the transforming growth factor-alpha, epidermal growth factor (EGF); nerve growth factor, neurotrophin; and insulin-like growth factor-1 (IGF-1), insulin systems since these have been extensively studied in recent years. The three receptor types, which bind these ligands, are the epidermal growth factor receptor family (erb B receptors), the neurotrophin or Trk receptor family, and IGF-1/insulin receptors, respectively, and represent three distinct members of the tyrosine kinase superfamily. For each of these, formation of the ligand-receptor complex initiates the signal transduction cascade through autophosphorylation by the intracellular tyrosine kinase domain. The extracellular portion of the receptor that contains the ligand binding domain in these systems varies significantly in organization in each case. For the EGF receptor system, ligand binding induces homo- and heterodimerization of the receptor leading to activation of the intracellular kinase. For the Trk receptor system, homodimerization of receptors has been shown to occur, although a second receptor, p75, is also required for high affinity binding of neurotrophins and for enhanced sensitivity of tyrosine kinase activation at low ligand concentrations. The IGF-1 and insulin receptors exist as covalent cross-linked dimers where each monomer is composed of two subunits. The aim of this review is also to discuss the mechanism of ligand-receptor interaction for each of these cases; however, since no structural information is yet available for the ligand-receptor complex, the discussion will largely be centered on the molecular requirements of ligand binding. As these receptors are activated through the ligand binding site on the extracellular domain, this represents a possible target for pharmacological intervention by inhibition or stimulation of this portion of the receptor. Thus from a drug design perspective, the focus of this review is to discuss progress in the development of agonists or antagonists of the ligand for these receptors.

Mini-proinsulin and mini-IGF-I: homologous protein sequences encoding non-homologous structures

Protein minimization highlights essential determinants of structure and function. Minimal models of proinsulin and insulin-like growth factor I contain homologous A and B domains as single-chain analogues. Such models (designated mini-proinsulin and mini-IGF-I) have attracted wide interest due to their native foldability but complete absence of biological activity. The crystal structure of mini-proinsulin, determined as a T3R3 hexamer, is similar to that of the native insulin hexamer. Here, we describe the solution structure of a monomeric mini-proinsulin under physiologic conditions and compare this structure to that of the corresponding two-chain analogue. The two proteins each contain substitutions in the B-chain (HisB10-->Asp and ProB28-->Asp) designed to destabilize self-association by electrostatic repulsion; the proteins differ by the presence or absence of a peptide bond between LysB29 and GlyA1. The structures are essentially identical, resembling in each case the T-state crystallographic protomer. Differences are observed near the site of cross-linking: the adjoining A1-A8 alpha-helix (variable among crystal structures) is less well-ordered in mini-proinsulin than in the two-chain variant. The single-chain analogue is not completely inactive: its affinity for the insulin receptor is 1500-fold lower than that of the two-chain analogue. Moreover, at saturating concentrations mini-proinsulin retains the ability to stimulate lipogenesis in adipocytes (native biological potency). These results suggest that a change in the conformation of insulin, as tethered by the B29-A1 peptide bond, optimizes affinity but is not integral to the mechanism of transmembrane signaling. Surprisingly, the tertiary structure of mini-proinsulin differs from that of mini-IGF-I (main-chain rms deviation 4.5 A) despite strict conservation of non-polar residues in their respective hydrophobic cores (side-chain rms deviation 4.9 A). Three-dimensional profile scores suggest that the two structures each provide acceptable templates for threading of insulin-like sequences. Mini-proinsulin and mini-IGF-I thus provide examples of homologous protein sequences encoding non-homologous structures.

Inactive conformation of an insulin despite its wild-type sequence

The peptide group between residues B24 and B25 of insulin was replaced by an ester bond. This modification only in the backbone was meant to eliminate a structurally important H-bond between the amide proton of B25 and the carbonyl oxygen of A19, and consequently to enhance detachment of the C-terminal B-chain from the body of the molecule, exposing the underlying A-chain. According to a model derived from the effects of side-chain substitutions, main-chain shortening, and crosslinking, this conformational change is prerequisite for receptor binding. Contrary to the expectation that increased flexibility would increase receptor binding and activity, depsi-insulin ([B24-B25 CO-O]insulin) has turned out be only 3-4% potent. In search of an explanation for this observation, the solution structure of depsi-insulin was determined by two-dimensional 1H-NMR spectroscopy. It was found that the loss of the B25-A19 H-bond does not entail detachment of the C-terminal B-chain. On the contrary, it is overcompensated by a gain in hydrophobic interaction achieved by insertion of the Phe B25 side chain into the molecule's core. This is possible because of increased rotational freedom in the backbone owing to the ester bond. Distortion of the B20-B23 turn and an altered direction of the distal B-chain are consequences that also affect self-association. The exceptional position of the B25 side chain is thus the key feature of the depsi-insulin structure. Being buried in the interior, it is not available for guiding the interaction with the receptor, a crucial role attributed to it by the model. This seems to be the main reason why the structure of depsi-insulin represents an inactive conformation.