Search

The infrared spectroscopy and dynamics of -CO labels in wild type and mutant insulin monomer and dimer are characterized from molecular dynamics simulations using validated force fields. It is found that the spectroscopy of monomeric and dimeric forms in the region of the amide-I vibration differs for residues B24-B26 and D24-D26, which are involved in dimerization of the hormone. Also, the spectroscopic signatures change for mutations at position B24 from phenylalanine, which is conserved in many organisms and is known to play a central role in insulin aggregation, to alanine or glycine. Using three different methods to determine the frequency trajectories (solving the nuclear Schrödinger equation on an effective 1-dimensional potential energy curve, using instantaneous normal modes, and using parametrized frequency maps) leads to the same overall conclusions. The spectroscopic response of monomeric WT and mutant insulin differs from that of their respective dimers, and the spectroscopy of the two monomers in the dimer is also not identical. For the WT and F24A and F24G monomers, spectroscopic shifts are found to be ∼20 cm -1 for residues (B24-B26) located at the dimerization interface. Although the crystal structure of the dimer is that of a symmetric homodimer, dynamically the two monomers are not equivalent on the nanosecond time scale. Together with earlier work on the thermodynamic stability of the WT and the same mutants, it is concluded that combining computational and experimental infrared spectroscopy provides a potentially powerful way to characterize the aggregation state and dimerization energy of modified insulins.

The monomer-dimer equilibrium for insulin is one of the essential steps in forming the receptor-binding competent monomeric form of the hormone. Despite this importance, the thermodynamic stability, in particular for modified insulins, is quite poorly understood, in part, due to experimental difficulties. This work explores one- and two-dimensional infrared spectroscopy in the range of the amide-I band for the hydrated monomeric and dimeric wild-type hormone. It is found that for the monomer the frequency fluctuation correlation function (FFCF) and the one-dimensional infrared spectra are position sensitive. The spectra of the -CO probes at the dimerization interface (residues Phe24, Phe25, and Tyr26) split and indicate an asymmetry despite the overall (formal) point symmetry of the dimer structure. Also, the decay times of the FFCF for the same -CO probe in the monomer and the dimer can differ by up to 1 order of magnitude, for example, for residue PheB24, which is solvent exposed for the monomer but at the interface for the dimer. The spectroscopic shifts correlate approximately with the average number of hydration waters and the magnitude of the FFCF at time zero. However, this correlation is only qualitative due to the heterogeneous and highly dynamical environment. Based on density functional theory calculations, the dominant contribution for solvent-exposed -CO is found to arise from the surrounding water (∼75%), whereas the protein environment contributes considerably less. The results suggest that infrared spectroscopy is a positionally sensitive probe of insulin dimerization, in particular in conjunction with isotopic labeling of the probe.

Insulin dimerization and aggregation play important roles in the endogenous delivery of the hormone. One of the important residues at the insulin dimer interface is Phe B24 , which is an invariant aromatic anchor that packs toward its own monomer inside a hydrophobic cavity formed by Val B12 , Leu B15 , Tyr B16 , Cys B19 , and Tyr B26 . Using molecular dynamics and free-energy simulations within explicit solvent, the structural and dynamical consequences of mutations of Phe at position B24 to glycine (Gly), alanine (Ala), and d-Ala and the des-PheB25 variant are quantified. Consistent with experiments, it is found that the Gly and Ala modifications lead to insulin dimers with reduced stability by 4 and 5 kcal/mol from thermodynamic integration and 4 and 8 kcal/mol from results using molecular mechanics-generalized Born surface area, respectively. Given the experimental difficulties to quantify the thermodynamic stability of modified insulin dimers, such computations provide a valuable complement. Interestingly, the Gly mutant exists as a strongly and a weakly interacting dimer. Analysis of the molecular dynamics simulations shows that this can be explained by water molecules that replace direct monomer-monomer H-bonding contacts at the dimerization interface involving residues B24 to B26. It is concluded that such solvent molecules play an essential role and must be included in future insulin dimerization studies.

Insulin, a protein critical for metabolic homeostasis, provides a classical model for protein design with application to human health. Recent efforts to improve its pharmaceutical formulation demonstrated that iodination of a conserved tyrosine (Tyr B26 ) enhances key properties of a rapid-acting clinical analog. Moreover, the broad utility of halogens in medicinal chemistry has motivated the use of hybrid quantum- and molecular-mechanical methods to study proteins. Here, we (i) undertook quantitative atomistic simulations of 3-[iodo-Tyr B26 ]insulin to predict its structural features, and (ii) tested these predictions by X-ray crystallography. Using an electrostatic model of the modified aromatic ring based on quantum chemistry, the calculations suggested that the analog, as a dimer and hexamer, exhibits subtle differences in aromatic-aromatic interactions at the dimer interface. Aromatic rings (Tyr B16 , Phe B24 , Phe B25 , 3-I-Tyr B26 , and their symmetry-related mates) at this interface adjust to enable packing of the hydrophobic iodine atoms within the core of each monomer. Strikingly, these features were observed in the crystal structure of a 3-[iodo-Tyr B26 ]insulin analog (determined as an R 6 zinc hexamer). Given that residues B24-B30 detach from the core on receptor binding, the environment of 3-I-Tyr B26 in a receptor complex must differ from that in the free hormone. Based on the recent structure of a "micro-receptor" complex, we predict that 3-I-Tyr B26 engages the receptor via directional halogen bonding and halogen-directed hydrogen bonding as follows: favorable electrostatic interactions exploiting, respectively, the halogen's electron-deficient σ-hole and electronegative equatorial band. Inspired by quantum chemistry and molecular dynamics, such "halogen engineering" promises to extend principles of medicinal chemistry to proteins., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

A new computational approach is proposed to probe the importance of residue side chains for the stability of a protein fold. Computational mutations to estimate protein stability (CMEPS) is based on the notion that the binding free energy corresponding to the complexation of a given side chain, considered as a "pseudo-ligand" of the wild type protein, reflects the importance of this side chain to the thermodynamic stability of the protein. The contribution of a particular side chain to the folding energy is estimated according to the molecular mechanics-generalized born surface area MM-GBSA approach, using a single molecular dynamics simulation trajectory of the wild type protein. CMEPS is a first principles method which does not contain any adjustable parameter that could be fitted to experimental data. The approach is first validated for Barnase and the B1 domain of protein L, for which a correlation coefficient R = 0.73, between experimental and CMEPS calculated DeltaDeltaG values, is found and then applied to the insulin monomer. In the present application, CMEPS replaces each amino acid by an alanine residue. Therefore, most mutations lead to cavities in the protein. From this the change in stability can be correlated with increased cavity volume. For insulin, this correlation is very similar compared with data previously analyzed for T4 lysozyme from an experiment for buried apolar side chains. There, the increased cavity volume has been related to the hydrophobic effect. However, since CMEPS uses the energetics in terms of electrostatic and van der Waals interactions (and not the hydrophobic effect which is difficult to relate to physical interactions), it is possible to study the effect of mutations of polar and solvent accessible side chains. According to CMEPS, residues Leu A16, Tyr A19, Leu B11, Leu B15, and Arg B22 are most important for the stability of the monomeric insulin fold. This is in agreement with experimental data. As a consequence, mutation of these residues may lead to misfolded and inactive insulin analogues., (Copyright 2006 Wiley Periodicals, Inc.)

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.

A calculation of the binding free energy for the dimerization of insulin has been performed using the molecular mechanics-generalized Born surface area approach. The calculated absolute binding free energy is -11.9 kcal/mol, in approximate agreement with the experimental value of -7.2 kcal/mol. The results show that the dimerization is mainly due to nonpolar interactions. The role of the hydrogen bonds between the 2 monomers appears to give the direction of the interactions. A per-atom decomposition of the binding free energy has been performed to identify the residues contributing most to the self association free energy. Residues B24-B26 are found to make the largest favorable contributions to the dimerization. Other residues situated at the interface between the 2 monomers were found to make favorable but smaller contributions to the dimerization: Tyr B16, Val B12, and Pro B28, and to an even lesser extent, Gly B23. The energy decomposition on a per-residue basis is in agreement with experimental alanine scanning data. The results obtained from a single trajectory (i.e., the dimer trajectory is also used for the monomer analysis) and 2 trajectories (i.e., separate trajectories are used for the monomer and dimer) are similar., ((c) 2005 Wiley-Liss, Inc.)

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.

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., (Copyright 2004 Wiley-Liss, Inc.)