Conformational analysis by nuclear magnetic resonance: insulin

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.

Ligand escape pathways and (un)binding free energy calculations for the hexameric insulin-phenol complex

Cooperative binding of phenolic species to insulin hexamers is known to stabilize pharmaceutical preparations of the hormone. Phenol dissociation is rapid on hexamer dissolution timescales, and phenol unbinding upon dilution is likely the first step in the conversion of (pharmaceutical) hexameric insulin to the active monomeric form upon injection. However, a clear understanding of the determinants of the rates of phenol unbinding remains obscure, chiefly because residues implicated in phenol exchange as determined by NMR are not all associated with likely unbinding routes suggested by the best-resolved hexamer structures. We apply random acceleration molecular dynamics simulation to identify potential escape routes of phenol from hydrophobic cavities in the hexameric insulin-phenol complex. We find three major pathways, which provide new insights into (un)binding mechanisms for phenol. We identify several residues directly participating in escape events that serve to resolve ambiguities from recent NMR experiments. Reaction coordinates for dissociation of phenol are developed based on these exit pathways. Potentials of mean force along the reaction coordinate for each pathway are resolved using multiple independent steered molecular dynamics simulations with second-order cumulant expansion of Jarzynski's equality. Our results for DeltaF agree reasonably well within the range of known experimental and previous simulation magnitudes of this quantity. Based on structural analysis and energetic barriers for each pathway, we suggest a plausible preferred mechanism of phenolic exchange that differs from previous mechanisms. Several weakly-bound metastable states are also observed for the first time in the phenol dissociation reaction.

H NMR Investigation of Thermally Triggered Insulin Release from Poly(N-isopropylacrylamide) Microgels

We describe investigations of insulin release from thermoresponsive microgels using variable temperature (1)H NMR. Microgel particles composed of poly(N-isopropylacrylamide) were loaded with the peptide via a swelling technique, and this method was compared to simple equilibrium partitioning. Variable temperature (1)H NMR studies suggest that the swelling loading method results in enhanced entrapment of the peptide versus equilibrium partitioning. A centrifugation-loading assay supports this finding. Pseudo-temperature jump (1)H NMR measurements suggest that the insulin release rate is partially decoupled from microgel collapse. These types of direct release investigations could prove to be useful methods in the future design of controlled macromolecule drug delivery devices.

Binding of phenol to R6 insulin hexamers

Small amounts of phenolic compounds are being used as preservatives in pharmaceutical insulin preparations. It has been shown previously that these compounds bind to specific sites on the insulin hexamer and act as allosteric effectors, inducing a transformation of the T6 hexamer to the R6 hexamer, via a T3R3 intermediate. In this article, the crystal structures of eight different insulin derivatives, all in the phenol-containing R6 form, are analyzed with respect to their phenol-binding sites. While six phenol molecules are normally bound per insulin hexamer, one of the engineered insulins appears to contain only three phenols but yet exists in an R6 conformation. This observation provides additional evidence for an inherent nonequivalence of the two trimers in the insulin hexamer. The unusual observation of a seventh phenol molecule bound to the hexamer of crystalline A21Gly-B31,B32Arg2 insulin (HOE 901), a long-acting derivative currently undergoing phase III clinical trials, provides a partial explanation for its protracted activity.

NMR triple-quantum filtered relaxation analysis of 17O-water in insulin solutions: an insight into the aggregation of insulin and the properties of its bound water

Transverse triple-quantum filtered NMR spectroscopy (TTQF) of 17O-water was used to study the properties of water in insulin solutions at different Zn2+ concentrations and pH values. It was established that strongly bound water molecules are already present in Zn-free insulin. On the assumption that the effective correlation time of a strongly bound water molecule, tau sb, is 10 ns, the apparent number of strongly bound water molecules was approximately 3 to 4 per insulin monomer. Addition of Zn2+ equivalent to approximately 2 g-atoms per hexamer did not produce substantial increases in the overall 17O-water TTQF signal intensity and apparent fraction of bound water. The dramatic enhancement of the TTQF signals observed for samples with a Zn2+/hexamer ratio greater than approximately 2:1 could be attributed to the increase in correlation time of the strongly bound water, due to the formation of higher-order oligomers of the protein.

Distinction of structural reorganisation and ligand binding in the T<==>R transition of insulin on the basis of allosteric models

Two allosteric models are presented for the T<==>R transition of insulin hexamers in the presence of phenolic ligands which are based on existing experimental information. The transition mainly involves residues 1-8 of the B-chain, i.e. 15% of the molecule, which are extended in the T- and helical in the R-state. The main facts to be accounted for are: 1) the transition is undergone trimer-wise; 2) the transition of the second trimer is disadvantaged compared to the first one; 3) the subunits of a trimer undergo transition in a cooperative process; 4) binding sites for phenolic ligands only exist in R3 trimers; 5) ligands shift the equilibrium by arresting the R-state; 6) the ligand is accommodated in a pocket made up between two adjacent subunits; 7) binding one ligand molecule extends the lifetime of the two other binding sites of a trimer; 8) only ligand-free trimers can undergo transitions. The two models allowed for CD spectroscopic titrations of zinc and cobalt insulin with phenol and m-cresol to be assessed in terms of structural reorganisation and ligand binding, and for the respective standard free energy differences to be calculated. delta G degrees for the reorganisation of the first timer in zinc-insulin is about 8 kJ/mol, and for that of the second trimer, 21kJ/mol. The corresponding values for cobalt-insulin are 12 and 24 kJ/mol, respectively. For both zinc- and cobalt-insulin, the delta G degrees for phenol and m-cresol binding is about -18 kJ/mol. Both models are equally compatible with the titration data.

The solution structure of a monomeric insulin. A two-dimensional 1H-NMR study of des-(B26-B30)-insulin in combination with distance geometry and restrained molecular dynamics

The solution conformation of des-(B26-B30)-insulin (DPI) has been investigated by 1H-NMR spectroscopy. A set of 250 approximate interproton distance restraints, derived from two-dimensional nuclear Overhauser enhancement spectra, were used as the basis of a structure determination using distance geometry (DG) and distance-bound driven dynamics (DDD). Sixteen DG structures were optimized using energy minimization (EM) and submitted to short 5-ps restrained molecular dynamics (RMD) simulations. A further refinement of the DDD structure with the lowest distance errors was done by energy minimization, a prolonged RMD simulation in vacuo and a time-averaged RMD simulation. An average structure was obtained from a trajectory generated during 20-ps RMD. The final structure was compared with the des-(B26-B30)-insulin crystal structure refined by molecular dynamics and the 2-Zn crystal structure of porcine insulin. This comparison shows that the overall structure of des-(B26-B30)-insulin is retained in solution with respect to the crystal structures with a high flexibility at the N-terminal part of the A chain and at the N-terminal and C-terminal parts of the B chain. In the RMD run a high mobility of Gly A1, Asn A21 and of the side chain of Phe B25 is noticed. One of the conformations adopted by des-(B26-B30)-insulin in solution is similar to that of molecule 1 (Chinese nomenclature) in the crystal structure of porcine insulin.

Two-dimensional NMR studies on des-pentapeptide-insulin. Proton resonance assignments and secondary structure analysis

The shortened analogue of insulin, des-(B26-B30)-pentapeptide insulin, has been characterized by two-dimensional 1H NMR. The 1H resonance assignments and the secondary structure in water solution are discussed The results indicate that the secondary structure in solution is very similar to that reported for the crystalline state. A high flexibility of both A and B chains is observed. Of the two conformations seen in the 2-Zn insulin crystals and indicated as molecules 1 and 2 (Chinese nomenclature), the structure of the analogue is more similar to that of molecule 1.

Sequence-specific 1H-NMR assignments for the aromatic region of several biologically active, monomeric insulins including native human insulin

The aromatic region of the 1H-FT-NMR spectrum of the biologically fully-potent, monomeric human insulin mutant, B9 Ser----Asp, B27 Thr----Glu has been investigated in D2O. At 1 to 5 mM concentrations, this mutant insulin is monomeric above pH 7.5. Coupling and amino acid classification of all aromatic signals is established via a combination of homonuclear one- and two-dimensional methods, including COSY, multiple quantum filters, selective spin decoupling and pH titrations. By comparisons with other insulin mutants and with chemically modified native insulins, all resonances in the aromatic region are given sequence-specific assignments without any reliance on the various crystal structures reported for insulin. These comparisons also give the sequence-specific assignments of most of the aromatic resonances of the mutant insulins B16 Tyr----Glu, B27 Thr----Glu and B25 Phe----Asp and the chemically modified species des-(B23-B30) insulin and monoiodo-Tyr A14 insulin. Chemical dispersion of the assigned resonances, ring current perturbations and comparisons at high pH have made possible the assignment of the aromatic resonances of human insulin, and these studies indicate that the major structural features of the human insulin monomer (including those critical to biological function) are also present in the monomeric mutant.

Changes in secondary structure follow the dissociation of human insulin hexamers: a circular dichroism study

Vacuum UV circular dichroism spectra measured down to 178 nm for hexameric 2-zinc human insulin, zinc-free human insulin, and the two engineered and biologically active monomeric mutants, [B/S9D] and [B/S9D,T27E] human insulin, show significant differences. The secondary structure analysis of the 2-zinc human insulin (T6) in neutral solution was determined: 57% helix, 1% beta-strand, 18% turn, and 24% random coil. This is very close to the corresponding crystal structure showing that the solution and solid structures are similar. The secondary structure of the monomer shows a 10-15% increase in antiparallel beta-structure and a corresponding reduction in random coil structure. These structural changes are consistent with an independent analysis of the corresponding difference spectra. The advantage of secondary structure analyses of difference spectra is that the contribution of odd spectral features stemming mainly from side chain chromophores is minimized and the sensitivity of the analyses improved. Analysis of the CD spectra of T6 2-zinc, zinc-free human insulin and monomeric mutant insulin by singular value decomposition indicates that the secondary structure changes following the dissociation of hexamers into dimers and monomers are two-state processes.

Phenol stabilizes more helix in a new symmetrical zinc insulin hexamer

SINCE insulin was first shown by Scott to crystallize in the presence of zinc ions in 1934, a variety of Zn-containing insulin crystals have been grown. The structures of insulin in the related rhombohedral crystals of 2Zn-insulin and 4Zn-insulin have been solved and reveal that the molecule is a hexamer, organized as three dimers, each containing a 2-fold symmetry axis and held together by Zn ions. In 2Zn-insulin the hexamer is nearly symmetrical with the two axial Zn ions and the two molecules of the dimer related closely by a local 2-fold axis. But in 4Zn-insulin the two molecules in the dimer differ remarkably, creating an asymmetric 4Zn-hexamer in which one trimer is essentially equivalent to that in 2Zn-insulin and the other is different by virtue of an additional stretch of N-terminal helix between residues B1 and B8 (refs 6, 7). We report here the structure of a new symmetrical hexamer, in which all six molecules have the B1-B8 helix seen in 4Zn-insulin. Phenol molecules, found bonding specifically to each molecule, evidently stabilize this new helical conformation.

1H NMR studies of insulin: histidine residues, metal binding, and dissociation in alkaline solution

The shifts of the H2 histidine B5 and B10 resonances of 2-Zn insulin hexamer were followed in 2H2O by 1H NMR spectroscopy at 270 MHz from pH 9.85 to 7. The two resonances present at high pH, previously assigned to H2 histidine B5 and B10 residues, moved slightly downfield and split into four resonances at pH 8.95 and also at pH 7. By use of a paramagnetic broadening probe (Mn2+) and the addition of Zn2+ to metal-free insulin, it was deduced that the four resonances arose from histidines B10 and B5 in two different magnetic environments, probably either bound to Zn2+ or not bound to Zn2+. The pK' values of the B5 and B10 histidines were determined in 60% 2H2O-40% dioxan, in which insulin was soluble throughout the pH range, to be 7.1 and 6.8, respectively at 37 degrees C. Studies at higher pH indicated that at a concentration level suitable for 1H NMR (approximately 1 mM) at 37 degrees C in 2H2O the 2-Zn hexamer was largely dissociated to dimer at pH 10.3 and to monomer at pH 10.8. Addition of paramagnetic shift probe Ni2+ to metal-free insulin caused changes to the spectrum similar to those produced on addition of diamagnetic Zn2+. Addition of Co2+ gave a different result, but there was no paramagnetic shift of the H2 histidine B10 resonance, probably because of rapid exchange at the binding site. Addition of Cd2+ and of Cd2+ and Ca2+ produced changes that were similar to each other but were different from those observed on addition of Zn2+, probably due to the binding of Cd2+ and Ca2+ at glutamate B13.

1H n.m.r. studies of insulin. Assignment of resonances and properties of tyrosine residues

The assignment of the aromatic 1H n.m.r. resonances of the four tyrosine residues of bovine 2-zinc insulin is reported, based on double resonance techniques, use of Hahn spin echo pulse sequences and examination of specific derivatives nitrated at tyrosines A14 and A19 as well as des-(B26-B30)-insulin. Titration curves of the four tyrosine residues show that residues A14 and B16 have normal pK' values of 10.3-10.6 in solution, consistent with their accessibility to solvent in monomer and dimer in the crystal. Tyrosine residues A19 and B26 have pK' values of 11.4 and exhibit other features in their titration curves that are consistent with limited accessibility to solvent and a nonpolar environment. The meta protons of residues B16 and B26 both observe the titration of a nearby tyrosine residue, probably A19. Interpretation of the n.m.r. data obtained in solution is consistent with the crystallographic data for the monomer and dimer obtained on insulin crystals [Blundell, Dodson, Hodgkin & Mercola (1972) Adv. Protein Chem. 26, 279-402].

A solution equivalent of the 2Zn----4Zn transformation of insulin in the crystal

Circular dichroic spectroscopy clearly reveals a solvent-induced conformational change of insulin in the presence of zinc ions. The spectral change corresponds to an increase in helix content. The transition observed in solution is an equivalent of the 2Zn----4Zn insulin transformation in the crystal. This is inferred from a series of observations. (1) The spectral effects are compatible with the refolding of the B-chain N-terminus into a helix known from crystal studies. (2) The spectral effects are induced by the very same conditions which are known to induce the 2Zn----4Zn insulin transformation in the crystal (i.e. threshold concentrations of NaCl, KSCN, NaI, for example). (3) They fail to be induced by the same conditions that fail to induce the crystal transformation (e.g. Ni2+ instead of Zn2+). It is concluded that the potential to undergo the transition resides in the hexamer since neither insulin dimers nor monomeric des-pentapeptideB26-30-insulin respond detectably to high halide concentration. Secondly the ability of zinc ions to accommodate tetrahedral coordination allows the transition which is not permitted by other divalent metal ions. Thirdly the transition is independent of the off-axial tetrahedral zinc coordination sites since it occurs in [AlaB5]insulin which lacks the B5 histidine necessary for their formation. A symmetrically rearranged hexamer thus appears possible with two tetrahedrally coordinated zinc ions on the threefold axis; this is consistent with the observation that in native insulin two zinc ions per hexamer are sufficient to produce the full spectral effect. The amount of additional helix derived from the circular dichroic spectral change, however, cannot settle whether the transition comprises only three or all six of the subunits to yield a symmetrical hexamer. Finally the transformation in solution evidently still occurs in an intramolecularly A1-B29-cross-linked insulin in spite of the partially reduced flexibility.

Photochemically induced nuclear polarization study of the accessibility of tyrosines in insulin

The accessible tyrosines of bovine insulin were studied by the photochemically induced dynamic nuclear polarization (photo-CIDNP) method. Tyrosine 1H nuclear polarization is observed in acidic, neutral, and basic solutions at all concentrations studied, in the absence of added salts as well as in the presence of 0.05-0.1 M chloride or phosphate. At pH 2.1 in the presence of chloride, at concentrations of 640 microM and above, most of the nuclear polarization at delta 6.82 originates from one group of tyrosines. On the basis of the crystallographic model, these are assumed to be the A14 tyrosines. We explored the possibility of a genuine concentration dependence of the photo-CIDNP intensity of insulin due to aggregation. In order to discern between such effects and trivial kinetic effects traceable to the optical irradiation method, the effects of concentration changes on polarization were examined in three apparently nonassociating trypsin inhibitor proteins. In insulin, the intensity of Tyr-A 14 polarization changes slowly at concentrations above 1 mM, suggesting that these residues are similarly accessible in all association states. At insulin concentrations below 320 microM, additional tyrosine emission signals were observed. These signals are probably due to B16 and B26 tyrosines of monomers. Polarization transfer effects from Tyr-A14 are evident in the tetramer and hexamer. Enhanced absorption effects in the two histidines (B5 and B10) of the insulin monomer were observed at pH 10 in the presence of 0.1 M phosphate.