Spectroscopic and molecular dynamics simulation studies of the interaction of insulin with glucose

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.

Tracking Insulin Glycation in Real Time by Time-Resolved Emission Spectroscopy

The application of time-resolved fluorescence sensing to the study of heterogenic biomolecular systems remains challenging because of the complexity of the resulting photophysics. Measuring the time-resolved emission spectroscopy (TRES) spectra can provide a more informative alternative to the modeling of the fluorescence decay that is currently employed. Here, we demonstrate this approach by monitoring real-time changes in intrinsic insulin fluorescence by TRES as a straightforward probe to directly measure kinetics of insulin aggregation and glycation. Our findings hold promise for monitoring the storage of insulin and its application in the control of diabetes and may support the development of more effective therapeutics against amyloidosis.

Rapid Non-Enzymatic Glycation of the Insulin Receptor under Hyperglycemic Conditions Inhibits Insulin Binding In Vitro: Implications for Insulin Resistance

The causes of insulin resistance are not well-understood in either type 1 or type 2 diabetes. Insulin (INS) is known to undergo rapid non-enzymatic covalent conjugation to glucose or other sugars (glycation). Because the insulin receptor (IR) has INS-like regions associated with both glucose and INS binding, we hypothesize that hyperglycemic conditions may rapidly glycate the IR, chronically interfering with INS binding. IR peptides were synthesized spanning IR- associated INS-binding regions. Glycation rates of peptides under hyperglycemic conditions were followed over six days using matrix assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. INS conjugated to horse-radish peroxidase was used to determine INS binding to IR peptides in glycated and non-glycated forms. Several IR peptides were glycated up to 14% within days of exposure to 20-60 mM glucose. Rates of IR-peptide glycation were comparable to those of insulin. Glycation of four IR peptides significantly inhibits INS binding to them. Glycation of intact IR also decreases INS binding by about a third, although it was not possible to confirm the glycation sites on the intact IR. Glycation of the IR may therefore provide a mechanism by which INS resistance develops in diabetes. Demonstration of glycation of intact IR in vivo is needed.

Vitamin E induces regular structure and stability of human insulin, more intense than vitamin D3

Changes in human environment and lifestyle over the last century have caused a dramatic increase in the occurrence of diabetes. Research of past decades illustrated that vitamin D and E have a key role in the improvement of diabetes by reducing oxidative stress, protein glycosylation, insulin resistance and also improving beta cell function. Binding properties and conformational changes of human insulin upon interaction with vitamins D₃ and E (α-tocopherol) were investigated by spectroscopy, differential scanning calorimetry (DSC) and molecular dynamic simulation. Tyrosine fluorescence quenching studies indicates changes in the human insulin conformation in the presence of vitamins. Binding constants of vitamins D₃ and E for human insulin were determined to be 2.7 and 1.5 (×10^-5M^-1) and the corresponding average numbers of binding sites were determined to be 1.3 and 1.2, respectively. Far- and near-UV circular dichroism studies showed that vitamin E can significantly change the secondary and tertiary structures of human insulin via an increase in the content of α-helix structure. Results of DSC showed that both vitamins D₃ and E stabilize the structure of human insulin. Molecular dynamic simulation results indicated that vitamin D₃ decreases the helical and strand structural contents of human insulin, but vitamin E stabilizes more regular secondary structures such as helical and strand structural contents as shown by experimental results.

Effect of the Freezing Step in the Stability and Bioactivity of Protein-Loaded PLGA Nanoparticles Upon Lyophilization

PURPOSE

The freezing step in lyophilization is the most determinant for the quality of biopharmaceutics. Using insulin as model of therapeutic protein, our aim was to evaluate the freezing effect in the stability and bioactivity of insulin-loaded PLGA nanoparticles. The performance of trehalose, sucrose and sorbitol as cryoprotectants was evaluated.

METHODS

Cryoprotectants were co-encapsulated with insulin into PLGA nanoparticles and lyophilized using an optimized cycle with freezing at -80°C, in liquid nitrogen, or ramped cooling at -40°C. Upon lyophilization, the stability of protein structure and in vivo bioactivity were assessed.

RESULTS

Insulin was co-encapsulated with cryoprotectants resulting in particles of 243-394 nm, zeta potential of -32 to -35 mV, and an association efficiency above 90%. The cryoprotectants were crucial to mitigate the freezing stresses and better stabilize the protein. The insulin structure maintenance was evident and close to 90%. Trehalose co-encapsulated insulin-loaded PLGA nanoparticles demonstrated enhanced hypoglycemic effect, comparatively to nanoparticles without cryoprotectant and added with trehalose, due to a superior insulin stabilization and bioactivity.

CONCLUSIONS

The freezing process may be detrimental to the structure of protein loaded into nanoparticles, with negative consequences to bioactivity. The co-encapsulation of cryoprotectants mitigated the freezing stresses with benefits to protein bioactivity.

Effects of several quinones on insulin aggregation

Protein misfolding and aggregation are associated with more than twenty diseases, such as neurodegenerative diseases and metabolic diseases. The amyloid oligomers and fibrils may induce cell membrane disruption and lead to cell apoptosis. A great number of studies have focused on discovery of amyloid inhibitors which may prevent or treat amyloidosis diseases. Polyphenols have been extensively studied as a class of amyloid inhibitors, with several polyphenols under clinical trials as anti-neurodegenerative drugs. As oxidative intermediates of natural polyphenols, quinones widely exist in medicinal plants or food. In this study, we used insulin as an amyloid model to test the anti-amyloid effects of four simple quinones and four natural anthraquinone derivatives from rhubarb, a traditional herbal medicine used for treating Alzheimer's disease. Our results demonstrated that all eight quinones show inhibitory effects to different extent on insulin oligomerization, especially for 1,4-benzoquinone and 1,4-naphthoquinone. Significantly attenuated oligomerization, reduced amount of amyloid fibrils and reduced hemolysis levels were found after quinones treatments, indicating quinones may inhibit insulin from forming toxic oligomeric species. The results suggest a potential action of native anthraquinone derivatives in preventing protein misfolding diseases, the quinone skeleton may thus be further explored for designing effective anti-amyloidosis compounds.

Insulin: a model system for nanomedicine?

One of the predominant aims of insulin therapy for diabetes is to appropriately mimic physiological insulin secretion levels and their correlation with glucose concentration in healthy individuals. This report outlines current methods and their limitations in glycemic control and their possible relationship to insufficient knowledge about the structure and dynamics of the insulin hormone itself. Based on recent experimental and computational work, a possible approach to less-invasive insulin administration is sketched.

Chemical analysis of solid-state irradiated human insulin

PURPOSE

To study the chemical modifications induced upon irradiation of solid human insulin at radiosterilization doses and investigate the influence of the absorbed dose on radiolysis.

MATERIALS AND METHODS

Volatile radiolytic products were monitored by gas chromatography coupled with mass spectrometry (GC-MS) and non-volatile products by two different high performance liquid chromatography (HPLC) methods: the formation of higher molecular weight proteins was assessed by size exclusion liquid chromatography whereas assays for related compounds and chemical potency tests were carried out using reverse-phase HPLC-UV. Conformational changes were investigated by measurements of circular dichroism.

RESULTS

After gamma irradiation at 10 kGy, the recovery of insulin was 96.8%; higher molecular weight proteins accounted for 0.35% (relative peak area) and other related compounds (including A21 desamido insulin) represented 1.29%. No major structural changes and no volatile radiolytic compounds were detected.

CONCLUSION

Human insulin samples irradiated in the solid-state at 10 kGy (gamma rays) and 14 kGy (electron-beam) meet the European Pharmacopoeia requirements and can be considered as quite stable towards radiation from a chemical analysis viewpoint.

An insulin-sensing sugar-based fluorescent hydrogel

We have prepared a small library of amphiphiles, each comprising a polar carbohydrate head group attached through an N-terminal amino acid to a nonpolar pyrene tail group. One of these derivatives is sensitive to the presence of insulin in aqueous media.

A Comparison of the Dynamic Behavior of Monomeric and Dimeric Insulin Shows Structural Rearrangements in the Active Monomer

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.

Investigation of glucose binding sites on insulin

Possible insulin binding sites for D-glucose have been investigated theoretically by docking and molecular dynamics (MD) simulations. Two different docking programs for small molecules were used; Multiple Copy Simultaneous Search (MCSS) and Solvation Energy for Exhaustive Docking (SEED) programs. The configurations resulting from the MCSS search were evaluated with a scoring function developed to estimate the binding free energy. SEED calculations were performed using various values for the dielectric constant of the solute. It is found that scores emphasizing non-polar interactions gave a preferential binding site in agreement with that inferred from recent fluorescence and NMR NOESY experiments. The calculated binding affinity of -1.4 to -3.5 kcal/mol is within the measured range of -2.0 +/- 0.5 kcal/mol. The validity of the binding site is suggested by the dynamical stability of the bound glucose when examined with MD simulations with explicit solvent. Alternative binding sites were found in the simulations and their relative stabilities were estimated. The motions of the bound glucose during molecular dynamics simulations are correlated with the motions of the insulin side chains that are in contact with it and with larger scale insulin motions. These results raise the question of whether glucose binding to insulin could play a role in its activity. The results establish the complementarity of molecular dynamics simulations and normal mode analyses with the search for binding sites proposed with small molecule docking programs.