Computational and structural evidence for neurotransmitter-mediated modulation of the oligomeric states of human insulin in storage granules.

Human insulin is a pivotal protein hormone controlling metabolism, growth, and aging and whose malfunctioning underlies diabetes, some cancers, and neurodegeneration. Despite its central position in human physiology, the in vivo oligomeric state and conformation of insulin in its storage granules in the pancreas are not known. In contrast, many in vitro structures of hexamers of this hormone are available and fall into three conformational states: T6, T3Rf3, and R6 As there is strong evidence for accumulation of neurotransmitters, such as serotonin and dopamine, in insulin storage granules in pancreatic β-cells, we probed by molecular dynamics (MD) and protein crystallography (PC) if these endogenous ligands affect and stabilize insulin oligomers. Parallel studies independently converged on the observation that serotonin binds well within the insulin hexamer (site I), stabilizing it in the T3R3 conformation. Both methods indicated serotonin binding on the hexamer surface (site III) as well. MD, but not PC, indicated that dopamine was also a good site III ligand. Some of the PC studies also included arginine, which may be abundant in insulin granules upon processing of pro-insulin, and stable T3R3 hexamers loaded with both serotonin and arginine were obtained. The MD and PC results were supported further by in solution spectroscopic studies with R-state-specific chromophore. Our results indicate that the T3R3 oligomer is a plausible insulin pancreatic storage form, resulting from its complex interplay with neurotransmitters, and pro-insulin processing products. These findings may have implications for clinical insulin formulations.

Introduction

Insulin is one of the key n class="Species">human protein hormones that is responsible for the maintenance of metabolic homeostasis, with an influence on cell proliferation and regulation of aging (1, 2). Defects in insulin bioavailability or impaired insulin receptor signaling lead to different pathological conditions such as diabetes (3–5), cancers (6–8), and Alzheimer's disease (9).

Insulin is a 51-amino acid protein consisting of two n class="Chemical">disulfide-linked chains (A1-A21, B1-B30) and exerts its functions through binding as a monomer to the (αβ)2 heterodimer tyrosine-kinase-type insulin receptor (IR) (10, 11). Insulin is produced from a single chain pro-insulin and stored in pancreatic β-cells in storage granules (termed large dense core vesicles (LDCVs)) from which it is released into the bloodstream in response to elevated blood glucose levels. The first 3D crystal structure of insulin was described by Hodgkin and coworkers in 1969 (12) in the form of its hexameric assembly (with three and two-fold rotational symmetry (i.e. trimer of dimers), obtained in the presence of Zn2+ ions. Two Zn2+ ions were identified in the hexamer on its 3-fold axis, being coordinated by 3 imidazole side chains of HisB10. This and subsequent similar findings resulted in a generally accepted paradigm that the hexamer is the storage form of insulin in LDCVs in pancreatic β-cells, whereas upon its release into the bloodstream it dissociates to monomers, which represent the active form of the enzyme (13). The pioneering work of D. C. Hodgkin was followed by 3D descriptions of hundreds of in vitro studied insulin hexamers, dimers, and monomers (for reviews, see for example Refs. 14–16). However, the actual in vivo storage form of insulin in pancreatic β-cells LDCVs is still not known, and it is extrapolated from its in vitro structures. Surprisingly, there are more advances into the very elusive nature of insulin/IR interactions (17, 18) than into the in vivo form of this hormone, investigation of which presents formidable experimental challenges.

The in vitro data showed that insulin hexamers can be grouped into three structurally distinct states/families: T6, T3Rf3, and R6 (15, 16), which differ by the conformation of the B1-B6 N termini of the hormone B-chain (see Fig. 1 and supplemental Fig. S1). In the T6 hexamers, the initial B1-B6 segments of the B chains are fully extended followed by B7-n class="Gene">B10 type II′ β-turns that go into invariant B9-B19 α-helices. In contrast, the B1-B6 segments acquire an α-helical conformation in the R-state, extending the B9-B19 α-helix (19). The insulin Rf3 state is similar to the α-helical R6 state; however, the long R-state B-chain helix is shortened to the B3-B19 segment, whereas the B1-B3 residues depart (“fray”) from the α-helical fold (20, 21). The TRf transition can be induced by an increase of, for example, the concentration of Zn2+, SCN−, and Cl− ions (22) or by the presence of phenolic derivatives (e.g. phenol or similar small cyclic alcohols) at low concentrations (23–25). The RfR (or full TR) transition can be accomplished by a further increase of the concentration of phenolic ligands (19, 24). The octahedral Zn2+ coordination is also a signature of the T-state (assured by three HisB10 imidazoles and three water molecules), whereas the Zn2+ tetrahedral coordination (three HisB10 and, for example, a Cl− ion) is typical for the Rf/R-states (e.g. Ref. 26). This coordination switch results from a smaller space above the Zn2+ ion in the R-state that is obstructed there by a longer B1-B8 helical part of the B9-B19-helix and a newly formed binding pocket for cyclic alcohols on the dimer interface, formed by the side chains of A6–7, A9–11, A16, B7, B11, and B5–7 (the so-called site I, or phenol main-binding site) However, variations of the Zn2+ coordination spheres are also observed frequently (27).

Figure 1.

General structural organization of the three main forms of insulin hexamer: T6 (A), T3Rf3 (B), and R6 (C) in the hexamer top view along its 3-fold symmetry axis. Insulin B chains are in white, A chains are in pink; the chains of one representative insulin dimer are indicated in A in light blue (B-chains) and yellow (A-chains). Some N and C termini of one dimer are also indicated in A and only B-chain N termini in B and C, with * corresponding to the symmetry-related monomer within a dimer. The B1-B8 segments of the B chains that contribute to the largest structural changes in TR transitions are in red (top of the hexamer) and in green (bottom of the hexamer). The Zn2+ ion is in blue. D–F, representative insulin monomers in T, Rf, and R-state, respectively, with a coloring code as in A–C. Typical, main ligand binding sites I and III in the R insulin form are indicated in C as S1 and S3 respectively.

General structural organization of the three main forms of insulin hexamer: T6 (A), T3Rf3 (B), and R6 (C) in the hexamer top view along its 3-fold symmetry axis. n class="Gene">Insulin B chains are in white, A chains are in pink; the chains of one representative insulin dimer are indicated in A in light blue (B-chains) and yellow (A-chains). Some N and C termini of one dimer are also indicated in A and only B-chain N termini in B and C, with * corresponding to the symmetry-related monomer within a dimer. The B1-B8 segments of the B chains that contribute to the largest structural changes in TR transitions are in red (top of the hexamer) and in green (bottom of the hexamer). The Zn2+ ion is in blue. D–F, representative insulin monomers in T, Rf, and R-state, respectively, with a coloring code as in A–C. Typical, main ligand binding sites I and III in the R insulin form are indicated in C as S1 and S3 respectively.

Although site I is the main binding cavity for all phenol-like ligands, they have also been identified in other regions of the hexamers; e.g. site II, formed by side chains B9, n class="Gene">B12, B16, B17, B9, B10, B13; site I/II formed by B10, B14, B16, B17, and B9 (28). Hexamer surface-exposed (between insulin dimers) phenolic site III (A14, A17) has been observed as well (29). The serendipitous character of the discovery of the phenol-stabilized R-state resulted from the bactericide-like applications of phenol in clinical formulations of insulin. Nevertheless, the R-like states of the hormone have clinical importance, as the Rf/R-trimers/hexamers are more stable (30). This is mostly due to the slower Zn2+-solvent exchange, which results in a higher hexamer-monomer dissociation constant. However, the physiological relevance of T/Rf/R states for an effective insulin/IR binding and for the storage form of the hormone (e.g. protection against proteolysis) is still unknown.

Here, we attempt to shed light on a possible conformation of the in vivo storage form of insulin in n class="Gene">insulin-containing granules in pancreatic β-cells. Although the exact chemical composition of LDCVs still awaits full characterization, they contain Zn2+ ions and phenolic neurotransmitters such as dopamine and serotonin, which are involved in regulation of exocytosis of the granules and insulin release (31–33). Therefore, we hypothesize that both neurotransmitters can bind to, and affect the structure of insulin hexamers, assembled in the presence of zinc cations. We probe these effects in parallel via molecular dynamics simulations (MD) and protein X-ray crystallography. These were further supplemented by spectroscopic studies with the use of a sensitive chromophore, 4-hydroxy-3-nitrobenzoic acid (4H3N), which specifically binds to R-state insulin as indicated by a red shift of its absorption spectrum. Moreover, as arginine can accumulate in LDCVs upon processing of pro-insulin at its two Arg-rich sites (34, 35) and as arginine-rich polypeptides (spermine/spermidine) are used in some insulin injections formulations (36, 37), we also probe the structural effect of this amino acid on insulin oligomeric assembly.

Both the MD and X-ray serotonin-related studies point to a stable n class="Gene">insulin hexamer-serotonin-site I complex and its dynamic association with site III. Both methodologies excluded dopamine binding in site I. MD simulations also indicate preferential dopamine over serotonin binding to site III, whereas solid-state X-ray analyses yields only serotonin-site III complexes. In addition, crystal structures of serotonin-insulin hexamers loaded with arginine have been obtained that underline a possible more complex, physiological role of insulin-arginine interactions.

Results concerning site I complexes are corroborated by n class="Chemical">4H3N spectroscopic experiments, which provide apparent binding constants for interaction of insulin hexamers with phenolic ligands and arginine. Our findings thus provide a novel insight into the storage mechanism of insulin in pancreatic β-cells, also relevant for novel insulin formulations for clinical applications.

Results

MD simulations of serotonin, dopamine, and phenol binding in site I

Here, we address the question of whether two neurotransmitters, serotonin and n class="Chemical">dopamine, can substitute phenol in its insulin R6 hexamer binding site I. MD simulations were performed for the R6 insulin hexamer, where site I-bound phenol was systematically replaced with neurotransmitters (each ligand in two steric-clash-free orientations; see Fig. 2). As expected, the MD simulations yielded the benchmark phenol-bound insulin R6-state hexamer (InsPheR6) with a low root mean square deviation (r.m.s.d.) of ∼1.6–1.8 Å from the reference InsPheR6 structure, indicating that the system does not deviate significantly from the NMR structure during the simulation (Fig. 3).

Figure 2.

Starting conformations of phenol (

Figure 3.

r.m.s.d. of the protein backbone (heavy atoms) from the insulin R Labels 1 and 2 correspond to the two starting orientations of the neurotransmitter molecules.

Both serotonin (InsSerR6) and n class="Chemical">dopamine (InsDopR6) complexes with ligands orientation (1) behaved similarly to InsPheR6, with the backbone r.m.s.d. just above 2 Å. However, the starting orientation (2) of the ligands led to higher r.m.s.d. of ∼2.8 Å (Fig. 2) indicating somewhat larger structural changes. A more direct insight into the behavior of the ligand in site I can be obtained from a direct comparison of their r.m.s.d. (Fig. 4).

Figure 4.

r.m.s.d. from simulations of all phenolic ligands (serotonin or dopamine in two orientations) from the phenol-binding pockets in the insulin R Each line corresponds to the r.m.s.d. of one phenolic ligand from its starting position with respect to the O atom of CysA6 and the N atom of CysA11 involved in a hydrogen bond. A zero or small value of r.m.s.d. means a strong ligand, R6 hexamer hydrogen bond, whereas a significant increase of r.m.s.d. indicates breaking of this bond.

r.m.s.d. from simulations of all phenolic ligands (n class="Chemical">serotonin or dopamine in two orientations) from the phenol-binding pockets in the insulin R Each line corresponds to the r.m.s.d. of one phenolic ligand from its starting position with respect to the O atom of CysA6 and the N atom of CysA11 involved in a hydrogen bond. A zero or small value of r.m.s.d. means a strong ligand, R6 hexamer hydrogen bond, whereas a significant increase of r.m.s.d. indicates breaking of this bond.

The MD simulations confirm that phenol molecules are essentially fixed in site I with a minimal dynamic behavior. Two of the six n class="Chemical">phenols transiently break HBs (both with CysA6 and CysA11); however, they remain in site I, eventually reestablishing the original geometry of binding. Similarly, all but one of the serotonin molecules in the starting conformation (1) stay in site I (only one serotonin ligand temporarily leaves this cavity). However, serotonin in starting geometry (2) as well as dopamine in both starting orientations, break most of the crucial hydrogen bonds. Serotonin, but not dopamine, thus essentially mimics phenol interactions observed in InsPheR6 with its charged aminoethyl-group-forming HB with backbone carbonyl CO of CysA11, fitting well into the binding site I.

The strength of binding (i.e. the binding free energy and the corresponding dissociation constant) of phenol and the two neurotransmitters was calculated using the thermodynamic integn class="Species">ration method (see Table 1). The free energy calculations confirmed and quantified the pattern observed in direct MD simulations, indicating a similarly strong binding of phenol and serotonin to site I (K of 5.4 × 10−4 m and 8.1 × 10−4 m, respectively) as well as demonstrating a lack of any stable conformation of dopamine in this cavity (reflected in a positive binding free energy).

Table 1

Standard free energies of binding of phenol, dopamine, and serotonin, molecules to the binding site I (phenolic pocket) together with the corresponding dissociation constants

Ligand	ΔGb°	Kd
	kcal/mol	[M]
Phenol	−4.49 ± 1.55	5.4 × 10−4
Dopamine	1.10 ± 1.73	6.3
Serotonin	−4.24 ± 1.87	8.1 × 10−4

Standard free energies of binding of phenol, n class="Chemical">dopamine, and serotonin, molecules to the binding site I (phenolic pocket) together with the corresponding dissociation constants

MD simulations of serotonin and dopamine probing the hexamer surface sites

Initial simulations of the R6 insulin hexamer with six n class="Chemical">dopamine molecules located in the phenolic pockets indicated that dopamine does not bind strongly to the phenolic pocket. At the same time, several dopamine molecules left the phenolic pocket to bind to the hitherto unexplored site III at the surface of the insulin hexamer. Therefore, an additional set of calculations was performed to characterize the binding of neurotransmitters to the surface of the insulin R6 hexamer (for simulation details, see the supplemental information). These simulations showed that dopamine binds at three equivalent binding sites in between the symmetry-related dimers (Fig. 5, top), which correspond to the previously sporadically observed site III. Analogous MD simulations for serotonin showed a similar mode of binding in site III but with a spatially looser neurotransmitter distribution, which reflects its weaker binding to this site compared with dopamine (Fig. 5, bottom). Finally, the binding of neurotransmitters at the surface sites has only a minor effect on the structure of the hexamer.

Figure 5.

Dopamine ( The A/B-chains are in blue and yellow, respectively; the same isodensity value (∼50× the bulk concentration) was used for both neurotransmitters. Black circles depict the site III binding pockets.

Dopamine ( The A/B-chains are in blue and yellow, respectively; the same isodensity value (∼50× the bulk concentn class="Species">ration) was used for both neurotransmitters. Black circles depict the site III binding pockets.

Both dopamine and n class="Chemical">serotonin bind in site III in a pocket between GluA17 and Tyr14 from neighboring dimers (Fig. 6). The stronger binding of dopamine over serotonin shown by free energy calculations (see below) results from its interactions with TyrA14, a hydrogen bond of the phenolic OH to GluA17, and a salt bridge between the aminoethyl group and the carboxylic group of GluA17. Similar interactions are present for serotonin; however, the steric fit is not as good as for dopamine. Dopamine, serotonin, and phenol binding in site III were quantified by evaluation of the binding free energies and the corresponding dissociation constants (Fig. 7 and Table 2) using the umbrella sampling method, with symmetry and volume entropy corrections. In contrast to site I, dopamine is the strongest binder to site III (K 4.38 × 10−3 m), with serotonin and phenol showing a much weaker affinity for this surface site (K of 2.13 × 10−1 m and 2.59 × 10−1 m, respectively).

Figure 6.

Serotonin ( Detailed structures of amino acids involved in binding (GluA17 and TyrA14) are also shown.

Figure 7.

Free energy profiles of phenolic ligands entering the surface binding site III.

Table 2

Standard free energies of binding of phenol, dopamine, and serotonin molecules to the binding site III and the corresponding dissociation constants

W0 is free energy difference directly from umbrella sampling calculations, 〈‹ΔGsymm〉 is hexamer symmetry contribution, ΔGvol is volume entropy contribution (with respect to the standard state at 1 m), ΔG° is overall standard free energy of binding, and K is the corresponding dissociation constant.

Ligand in the binding site III	W0	〈‹ΔGsymm〉	ΔGvol	ΔGb°	Kd
	kcal/mol	kcal/mol	kcal/mol	kcal/mol	m
Phenol	−1.83 ± 0.28	−0.36	1.38 ± 0.22	−0.81 ± 0.36	2.59 × 10−1
Dopamine	−6.23 ± 0.67	−0.36	3.2 ± 0.2	−3.24 ± 0.70	4.38 × 10−3
Serotonin	−2.65 ± 0.55	−0.36	2.0 ± 0.1	−0.92 ± 0.56	2.13 × 10−1

Standard free energies of binding of phenol, n class="Chemical">dopamine, and serotonin molecules to the binding site III and the corresponding dissociation constants

Insulin-serotonin crystal complex

The crystal structure of insulin grown in the presence of n class="Chemical">serotonin and Zn2+ revealed a T3R3 hexamer (referred to here as InsSerT3R3), with six neurotransmitters and two Zn2+ ions per hexamer (Fig. 8). The minimum serotonin concentration that still yielded complex crystals was within the 35–40 mm range. The asymmetric unit of this crystal contains 14 molecules of insulin, hence providing reliable, independent multicopy structural evidence. The hexamer quaternary structure is indeed T3R3 as the B1-B3 N termini acquire here a more R-like (fully α-helical) fold rather than an Rf, “frayed,” conformation observed in T3Rf3 oligomers.

Figure 8.

Serotonin binding sites in insulin InsSerT Sites I and III are indicated; some N and C termini of A/B-chains (T-state: T, R-state: R) are shown. Serotonin in atom-colored coded, with C-atoms in yellow.

Serotonin binding sites in n class="Gene">insulin InsSerT Sites I and III are indicated; some N and C termini of A/B-chains (T-state: T, R-state: R) are shown. Serotonin in atom-colored coded, with C-atoms in yellow.

Coordination of both Zn2+ ions is tetrahedral regardless of the R or T protein environment with axial Cl− ligands. The overall fold of InsSerT3R3 is very similar to other T3Rf3 hexamers. For example, the r.m.s.d. of InsSerT3R3 from complexes with n class="Chemical">phenol (PDB ID 1MPJ Ref. 29), Tylenol (PDB ID 1TYL; Ref. 23), paraben (3 MTH; Ref. 29), and 4-hydroxy-benzamide (PDB ID 1BEN; Ref. 25) are 0.8911 Å, 0.9886 Å, 0.9105 Å, and 1.1570 Å, respectively. Also, this hexamer is not significantly different from phenol-free T3Rf3 hexamers induced by ions such as SCN− (PDB ID 2TCI, r.m.s.d. 0.9919 Å; Ref. 29) and Cl− (PDB ID 1G7A, r.m.s.d. 0.7765 Å (Ref. 38).

Serotonin occupies six sites in InsSerT3R3: three n class="Chemical">phenol “main” sites I and three hexamer surface sites III. It is tethered (Fig. 9) into site I via (typical for phenolic ligands) HBs of its OH group to -CO of CysA6 (range of 2.41–2.43 Å) and to -NH of CysA11 (2.93–3.17 Å); a weaker HB -OH contact to -CO of SerA9 (3.38–3.41 Å) further stabilizes this phenolic anchor. A significant HisB5-Nϵ2·π pyrrole HB (3.46–3.56 Å) contributes to the immobilization of the serotonin indole ring, which is assisted further by more subtle van der Waals (hydrophobic) interactions, especially with LeuB17 (3.45–3.59 Å). The aminoethyl side chain of serotonin points toward the surface of the hexamer, i.e. into the space/gap between two 3-fold related dimers. The -NζH amino group of this side chain is hydrogen-bonded to Oϵ1 of GluB21 (2.55–2.75 Å) and -CO of CysA11 (2.72–2.79 Å).

Figure 9.

Serotonin binding sites: Dashed lines, HB; *R and *T, correspond to symmetry-related dimer and T/R state of the monomer, respectively. Strong HisB5-Nϵ2·π-pyrrole-center contact is also shown.

Serotonin binding sites: Dashed lines, HB; *R and *T, correspond to symmetry-related dimer and T/R state of the monomer, respectively. Strong n class="Chemical">HisB5-Nϵ2·π-pyrrole-center contact is also shown.

In contrast to site I, the serotonin-binding mode to site III is n class="Species">rather dynamic, indicating its mobility on the InsSerT3R3 surface. Although serotonin is pincered there in a hydrophobic clamp of TyrA14 (π·π interactions ∼3.5 Å) and the side chain of LeuA13 (from a symmetry related dimer), the ligands' electron density is generally less well defined there. The serotonin phenolic OH group makes hydrogen bond to the Oϵ atoms of GluA17 (2.96–3.44 Å), and it can be stabilized by water-mediated HB to the hydroxyl of TyrA14. The HBs of the GluA17 carboxylic side chain to the guanidinium group of ArgB22 (2.75–2.82 Å) contributes further to the structural stability of the serotonin-OH environment. However, it is evident that the aromatic moieties of Tyr14 and serotonin can slide parallel to each other, indicating further adaptability of this mode of binding. Relative flexibility of site III is also seen in different modes of the HB network to the -Nζ group of the serotonin aminoethyl side chain. It is engaged primarily by HBs to the carboxyl side chain of GluA17 (2.65–2.76 Å) and via a water molecule to ArgB22, but this network of HBs is prone to disruption due to the flexibility of the serotonin surface-exposed side chain. In summary, the GluA14/ArgB22-pair are the 2-fold symmetry-related providers of HBs to serotonin (one pair to the ligand's OH group, the other one to the end of its side chain) in site III, but these HB networks are not fully symmetrical in overall geometry and strength due to a non-symmetrical ligand in this site. The swing of the other (not involved in serotonin binding) TyrA14 side chain, away from site III, breaks down the 2-fold symmetry of this interface even further. Finally, we note that all our dopamine/insulin co-crystallizations were unsuccessful, likely due to the oxidation of the ligand, despite the inclusion of some anti-oxidation agents in the media.

Insulin-serotonin-arginine crystal complex

A similar set of crystallization conditions yielded crystals of insulin in a ternary complex with n class="Chemical">serotonin and arginine (referred to here as InsSerArgT3R3). Arginine- and Zn2+-containing crystallizations of insulin were carried out in the presence and absence of the previously established optimum serotonin concentration (40 mm). Serotonin-arginine-containing solutions yielded several morphologically different, but crystallographically very isomorphous, crystal forms. Despite their similarity, two of these forms (f1, f2) (referred to as InsSerArgf1-T3R3 and InsSerArgf2-T3R3) are reported here due to a dynamic (i.e. with high level of disorder) nature of Arg/insulin binding, hence the need for more independent structural evidence about the nature of these interactions.

Both InsSerArg complexes appear in the T3R3 hexamer state, with all n class="Chemical">Zn2+ ions in a tetrahedral coordination with Cl− as an axial ion. They are very similar to the InsSerT3R3 complex, with overall r.m.s.d. values of 0.2193/0.2349 Å between these structures.

Despite the presence of arginine, the n class="Chemical">serotonin modes of binding in sites I are practically the same as in the Ins-Ser complex (Fig. 10). However, the site III shows more significant structural variety in binding of this neurotransmitter. Some serotonins in site III are bound in a fashion observed in InsSerT3R3; however, some of these ligands are flipped in both InsSerArgT3R3 crystal forms, i.e. the -OH and aminoethyl groups switch their positions. As both of them benefit from a similar HB network of symmetrical GluA14/ArgB22 tandems, the overall HB-connectivity is conserved here as well, with the exception of the Nζ atom of the aminoethyl side chain. It forms in some sites III a new HB with the hydroxyl group of TyrA14 (2.67 Å) besides maintaining a HB to the side chain of GluA14 as well.

Figure 10.

Serotonin and arginine binding sites in insulin InsSerArgT

Serotonin and n class="Chemical">arginine binding sites in insulin InsSerArgT

It has to be stressed that several serotonins in site III and n class="Chemical">arginine molecules (see below) in the InsSerArgT3R3 complex have been refined with half-occupancies and, in some cases, with more arbitrary modeling due to their partial definition in the electron density maps. This reflects the dynamic and mobile character of some ligand/hormone interactions in InsSerArgT3R3 hexamers.

Arginine molecules occupy the “T-state-half” of the trimer (Fig. 11), with their Cα moieties being part of the T-state trimer surface (on the level with the Cl− ion), whereas their side chains point toward the hexamer core, relatively parallel to its central axis. The aminocarboxy-Cα groups of n class="Chemical">Arg have been modeled in both InsSerArgT3R3 forms with some flipped/alternative conformations around the Cα atom.

Figure 11.

Arginine binding sites in the insulin InsSerArgT Labeling and color-coding is as in Fig. 8, with Cl− ion as yellow.

Arginine binding sites in the n class="Gene">insulin InsSerArgT Labeling and color-coding is as in Fig. 8, with Cl− ion as yellow.

The dominant stabilizations of the Arg Cα-end result from HBs of the carboxyl group to the NH of n class="Chemical">CysB7 (2.76–3.31 Å) and amino acid amine group to CO of HisB5 (2.86–3.18 Å). However, it is possible to model this part of arginine with alternative conformations, a likely effect of a dynamic character of these ligands at the wider, and more solvent-accessible, T3 side of the hexamer.

In contrast, the HBs of the Arg n class="Chemical">guanidinium groups are much better defined, especially in the InsSerArgf2T3R3 form. One of their main signatures is the extensive HB network with both GluB13 and HisB10 residues. Here, the GluB13 Oϵ1/Oϵ2 atom forms HBs to Arg Nη1 (2.11–2.33 Å) and Nδ1 of HisB10 (2.82–3.14 Å), whereas the second Nη2 atom of the guanidinium group is hydrogen-bonded to CO of HisB10 as well (2.71–2.82 Å). The last guanidinium Nϵ1 atom locks this group by HB to the hydroxyl Oγ of SerB9. Some of these HBs are broken in the InsSerArgf1T3R3 form, also indicating the increased mobility of this environment, especially some shifting of the arginine parallel to the hexamer 3-fold axis.

Determination of ligand K values by solution 4H3N assay

The interactions of insulin hexamers with n class="Chemical">phenol, serotonin, dopamine, and arginine in solution were investigated and quantified by spectroscopic studies, with the 4H3N chromophore, which binds exclusively to HisB10 sites in the R-state only, which is associated with the red shift of its absorption spectrum. This allows for the determination of the apparent binding constants (K) of selected ligands to insulin hexamers and estimation of the values of their maximum specific binding (Bmax), which can be considered here as a measure of the amount of the R state induced by the ligand. In addition, the Hill coefficients (h), which indicate the scale of cooperativity (likely very complex) between ligand and/or insulin its hexameric state upon ligand binding, were evaluated as well.

First, we determined the change in 4H3N absorption upon titn class="Species">ration of selected ligand to T6 insulin hexamer preincubated with this chromophore. Fitting of all curves (Fig. 12) gave K, Bmax, and h values (Table 3) for all types of experiments. The data from Table 3 are also presented in a bar plot in supplemental Fig. S2. The 4H3N spectroscopic data indicate that phenol is the strongest binder to insulin hexamers, about 4× stronger than serotonin (K values of 0.85 mm and 3.34 mm, respectively; see Table 3). In contrast, titration of the same T6 insulin hexamers by arginine and dopamine did not induced any changes in the 4H3N absorption spectra (data not shown), indicating a lack of binding to site I and, subsequently, inability to induce T6→T3R3/R6 transitions.

Figure 12.

Titration of insulin hexamers and 4H3N with phenol (A) or serotonin (B). Titration of insulin hexamers and 4H3N preincubated with 5 mm arginine with phenol (C) or serotonin (D). Titration of insulin hexamers and 4H3N preincubated with 5 mm dopamine with phenol (E) or serotonin (F). All measured binding curves are shown.

Table 3

Values of

Ligand	Kd ± S.D.	Bmax± S.D.	h ± S.D.
	mm
Phenol (n = 4)	0.86 ± 0.13	0.17 ± 0.01	2.22 ± 0.30
Phenol and 5 mm Arg (n = 3)	0.75 ± 0.01	0.10 ± 0.00	2.07 ± 0.19
Phenol and 5 mm dopamine (n = 3)	0.64 ± 0.06	0.13 ± 0.01	1.68 ± 0.18
Serotonin (n = 3)	3.34 ± 0.52	0.11 ± 0.08	1.63 ± 0.10
Serotonin and 5 mm Arg (n = 3)	3.12 ± 0.31	0.06 ± 0.00	2.64 ± 0.27
Serotonin and 5 mm dopamine (n = 3)	3.93 ± 1.28	0.05 ± 0.04	7.12 ± 7.50

Titration of n class="Gene">insulin hexamers and 4H3N with phenol (A) or serotonin (B). Titration of insulin hexamers and 4H3N preincubated with 5 mm arginine with phenol (C) or serotonin (D). Titration of insulin hexamers and 4H3N preincubated with 5 mm dopamine with phenol (E) or serotonin (F). All measured binding curves are shown.

The addition of phenol or n class="Chemical">serotonin to insulin/4H3N that was also preincubated with 5 mm Arg resulted, however, in different patterns of binding curves (Fig. 12, C and D), compared with sole ligands-insulin/4H3N titrations (Fig. 12, A and B). Similar trends were observed for phenol and serotonin titrations of insulin/4H3N preincubated with 5 mm dopamine (Fig. 12 E and 11F). It seems that in both cases the preincubation of insulin/4H3N with arginine or dopamine (which are not hormone Site I binders on their own) lowers the phenol K to 0.75 and 0.64 mm, respectively. In contrast, arginine and dopamine do not change the binding affinity of serotonin to insulin hexamers (Table 3, Fig. 12A).

The higher Bmax value for phenol than for n class="Chemical">serotonin (0.17 versus 0.11) may also indicate that serotonin is not able to fully induce the R6 state of the hexamers. Note that preincubation of insulin/4H3N with arginine or dopamine lowers the Bmax values of phenol or serotonin binding to insulin (Fig. 12B). The Hill coefficient of phenol binding to insulin hexamers (Fig. 12C) is also higher (h = 2.2) than that for serotonin (1.6), indicating a lower degree of cooperativity for serotonin-insulin interactions. Interestingly, the presence of arginine or dopamine lowers cooperativity (i.e. the h value) of phenol-insulin binding but increases the h value for serotonin interactions with the hormone. This effect seems to be even more noticeable in the case of dopamine (h = 7.1), but these data are accompanied by a high experimental error (see Table 3 and supplemental Fig. S2).

Discussion

Our parallel and independent MD simulations, X-ray crystallography, and ligand binding experiments provide a concise picture showing that the insulin hexamer can accommodate and be affected by ligands that are present in the hormone storage n class="Chemical">LDCV granules in pancreatic β-cells. These ligands can be of endogenous origin, such as neurotransmitters (serotonin or dopamine) involved in regulation of insulin secretion, or they can be byproducts of pro-hormone/pro-insulin processing (e.g. arginine) that accumulate in the granules.

All the above methods indicated that serotonin is a good n class="Chemical">phenol-compatible ligand for insulin hexamer site I, whereas dopamine is not a site I binder. Serotonin-site I binding engages its phenolic OH in “phenol-canonic” HBs to CO of CysA6 and NH of CysA11. This is despite an ∼29° tilt of its indole plane with respect to the phenol benzene ring, a movement of the aromatic moiety that is also observed (to a different degree) in Tylenol (23), paraben (29), and benzamide (25) complexes. However, the side chains of the latter phenolic ligands are oriented toward the central channel of the hexamer, whereas the aminoethyl group of serotonin points into the opposite direction, i.e. toward the surface of the hexamer. This results in further enhancement of serotonin stability in site I by fixing its side chain via HBs with the CO of CysA11 and carboxyl of GluB21. The fitting of serotonin in site I is completed with a strong HisB5-Nϵ2·π-pyrrole contact. The firmness of serotonin-site I interactions, which is propagated likely toward the N termini of the B-helix, combined also with a very good overall fit of serotonin into site I may be behind the induction of the R (i.e. not Rf) conformation of the T3R3 hexamer. This is rather unique, as this type of insulin quaternary arrangement occurs mainly in the T3Rf3 state. T3R3 hexamers were observed only within a R6→T3R3 retrograde conformational change (39–41) upon evaporation of phenol from R6 insulin crystals, which is a process that cannot occur for the non-volatile serotonin ligand. Therefore, the stability of the SerInsT3R3 complex may indicate that this hormone's conformer is indeed close to its putative native form in pancreatic LDCV (as postulated, for example, by Dunn (26)). Moreover, larger Bmax and Hill coefficient h values (obtained from 4H3N spectroscopy studies) for phenol as compared with serotonin suggest that the latter, despite its “perfect” fit into site I, cannot fully induce the insulin R6 state (as also confirmed by the lack of R6 serotonin hexamer crystals even at a very high ligand concentration). These findings further corroborate the potential physiological aspect of the T3R3 insulin conformer as the only possible oligomeric species in the presence of serotonin. Although MD simulations did not exclude the presence of a stable SerInsR6 conformer, the R6 state was postulated in the MD work a priori (due to its thermodynamic superiority in the presence of phenol ligands) to identify the potential ligand-binding sites rather than to investigate the impact of a particular ligand on the state of the insulin hexamer.

The lack of stable site I-dopamine complexes was demonstn class="Species">rated by all methods employed in this study. It seems that this is not due to a steric ligand-cavity incompatibility, as the dopamine dispersion contribution to free energy of binding in the MD simulations is more negative compared with phenol. Rather, as the superposition of dopamine with phenol or serotonin ligands in site I suggests, the lack of dopamine binding is due to an electrostatic mismatch, namely its inability to accommodate the charged aminoethyl side chain in the site I cavity. The smaller aromatic scaffold of dopamine thus does not allow its side chain to reach the hexamer surface and form HBs observed in serotonin complexes.

As the 4H3N spectroscopic applicability was limited here to the environment of n class="Chemical">HisB10 and, hence, the site I, the evidence for other insulin site-ligand interactions came solely from MD simulations and X-ray crystallography. Both showed serotonin binding to surface site III, whereas only MD calculations predicted the possibility of a stable site III-dopamine complex (note, however, that dopamine-insulin co-crystallization may be hampered by a rapid chemical degradation/oxidation of dopamine). Interestingly, MD calculations indicated superior dopamine binding to site III, with its K ∼50 times better than for phenol or serotonin. MD-predicted dopamine dominance at site III is interesting, as these simulations show its binding in a similar mode to serotonin-site III complex, which is observed both in MD simulations and crystal structures. Both neurotransmitters engage TyrA14 and GluA17/ArgB22 side chains from 2-fold symmetry-related dimers in the network of HB and π·π-stacking interactions. Moreover, crystal structures of InsSerT3R3 also underline the dynamical character of site III mode of binding, which involves flipping of the serotonin indole ring (i.e. swapping positions of -OH and the aminoethyl side chain) and a certain variability of the TyrA14-OH/GluA17/ArgB22 HB patterns.

The impact of these ligands, which are natural components of pancreatic β-cells insulin-storing n class="Chemical">LDCVs, was also extended from neurotransmitters to probing of arginine-insulin hexamer interactions. It was postulated earlier that arginine can accumulate in these granules upon processing of pro-insulin at its two sites (Arg-31–Arg-32 and Lys-64–Arg-65) (34, 35), which are cleaved by specific convertases during maturation of this hormone (42, 43). Here, the X-ray crystallography revealed that arginine can penetrate serotonin-containing hexamers (InsSerArgT3R3) without significant changes to their InsSerT3R3 state. Arginine is accommodated there in channels within the T3 trimer, parallel to the hexamer 3-fold axis and opposite to R3 trimer sites I filled by serotonin, which leads to filling the voids around Zn2+ and His-B10. Binding of arginine in InsSerArgT3R3 is dynamic with high mobility and alternative, slightly different ligand conformations. The flexibility of arginine is concentrated mostly around its -CαNH2COO− moieties, which form weak alternative HB to main chain groups of CysB7 and HisB5, also exhibiting some flipping around the Cα atom. In contrast, the guanidinium groups of the arginines are firmly engaged in extensive networks of HBs, which link them with side chains of GluB13, the side- and main-chain of HisB10, and the side chain of SerB9. This simultaneous binding of the two very different ligands to the T3R3 insulin hexamer polarizes this oligomer even further into two structurally distinct trimers, Arg3T3 and Ser6R3. Although Arg binding in the ArgT3 trimer has a notable impact on the stability of the usually mobile GluB13 side chains, the solvent-exposed flexible -CαNH2COO− moieties may contribute to the disorder of the N-terminal parts of B-helices, which is, however, untraceable up to the HisB5 site in both arginine-containing and arginine-free serotonin complexes. In contrast, the GluB13 side chains are fully disordered in the Ser6R3 trimer, whereas the serotonins in sites I facilitate fixing of the B1-B19 helices in the R-state.

The structural partition of Arg3T3 and n class="Chemical">Ser6R3 is more complex. First, although the site I pocket is mainly formed by two R-monomers, the T-state monomer also contributes here by LeuB17 and GluB21 toward serotonin aromatic rings and side-chain binding, respectively. Secondly, serotonin surface site III also mediates the structural interface between T3 and R3, as this binding site (formed by TyrA14/GluA17) is provided by the T- and R-monomers from the 2-fold symmetry related dimers.

Interestingly, 4H3N spectroscopic data indicate that the presence of n class="Chemical">arginine increases (∼13%) phenol (but not serotonin) affinity for site I and also increases the cooperative character of serotonin binding (as demonstrated by a higher Hill h factor). These structural and spectroscopic data provide evidence for the T3-R3 trimers cross-talk and, as postulated previously, for a heterotropic allostery within the insulin hexamer (22, 30, 44–47). However, the exact molecular detail behind this phenomenon has been elusive, likely involving a network of propagating longer range interactions, which fit the Seydoux, Malhotra, and Bernhard (SMB) cooperativity model (48). The stabilizing effect of guanidinium groups on GluB13-HisB10 and Zn2+ may be one of the initial steps in this process, priming the other half of the hexamer for more effective serotonin binding (i.e. an element of the so-called half-site reactivity). The stability of InsSerArgT3R3 and InsSerT3R3 hexamers in the presence of physiological ligands, i.e. trapping of these oligomers in a stable T3R3 state, agrees also with the finding that conformational fluctuations of the T3 trimer needed for the appearance of site I become constrained upon formation of T3R3 (49). This can be considered as one of the features of the negative cooperativity effect in the insulin hexamer, predicted by the SMB model as well.

The present 4H3N spectroscopic data also suggest a n class="Chemical">dopamine-lowering effect of phenol's K. As there is no evidence for dopamine site I binding, this allosteric effect may result from longer-range cross-talk of the dopamine-site III complex with the environment of site I. On the other hand, dopamine does not seem to significantly change the serotonin K, and therefore, some “competitive” crowding effect (e.g. dynamic occlusion of site I) between these two neurotransmitters cannot be excluded. It may also mean that dopamine-site III binding facilitates the diffusion of a smaller ligand (i.e. phenol) into site I, whereas the increased stability of the insulin R3-like trimer (via new HBs and interactions in site III) may obtrude its effect on binding of a larger ligand (such as serotonin).

It is worth mentioning that several tubular electron densities are also observed in InsSerArgT3R3 on the top of the T3 trimer. These are linked with the modeled n class="Chemical">arginine ligands, suggesting the presence of more of these molecules on the hexamer surface, in the form of a hydrogen-binding network of ligands. Some of this electron density is in contact with the T3 Zn2+ ion. Therefore, a contribution of the arginine -COO− group to the tetrahedral/monodentate coordination of Zn2+ (observed also in the insulin R6 structure; Ref. 41) cannot be excluded either. This direct arginine-Zn2+ interaction would hamper the T3→R3 transition even further, thus reinforcing the stability of the T3R3 hexamer.

Arginine is widely used in vitro as a nonspecific protein-folding stabilizing agent (49, 50), and so its presence in InsSern class="Chemical">ArgT3R3 structures may reflect some of its physiologically relevant roles as a ubiquitous insulin ligand. First, it may serve as an unspecific, granules-abundant, insulin-folding stabilizer and hexamer half-life enhancer. Second, it may act as a more specific endogenous ligand that modulates the hexamer cooperativity, also enhancing the stability of its particular T3R3 storage form. This more specific role of arginine may result, among others, from its guanidinium counterion-like neutralization of the repulsive system of the GluB13 side chains in the hexamer core. Moreover, a simultaneous arginine-mediated cross-linking of GluB13 side chains with Zn2+-binding HisB10, leading to stabilization of this structurally key region, further exemplifies the ligand-like role of this amino acid.

It has to be stressed that insulin crystallization in the granules can be also affected by the intragranular n class="Chemical">Zn2+ concentrations, which are modulated by some Zn2+ transporters (ZnT), especially by the disease-associated ZnT8 variants (51). However, the focus here was on insulin/neurotransmitter interactions, hence the experimental Zn2+ concentrations (0.4–8 mm) were maintained within typical insulin in vitro crystallization ranges. Interestingly, insulin crystallized here only in the form of 2Zn2+ hexamers, regardless the very high (25:1) Zn2+-hormone molar excess in the crystal-yielding conditions.

There are ongoing efforts for a better insight into the size and content of insulin storage granules (e.g. Ref. 52). Their exocytosis can also be morphologically and mechanistically quite heterogeneous (so-called “kiss-and-run” phenomenon), releasing varying amounts and forms of granule content and on different time scales (53). Therefore, the final amount and form of n class="Gene">insulin discharged to the circulation can be a very complex and multifactorial process; whether insulin-neurotransmitter-arginine hexamers are one of its variables remains to be seen.

Conclusions

We have demonstrated here that n class="Gene">insulin oligomeric forms in storage granules in β-pancreatic cells could be regulated by certain endogenous components of these vesicles, such as serotonin, dopamine, and arginine. They are able to (i) shift the insulin oligomeric equilibrium toward the T3R3 state, which therefore may be considered as the insulin storage form in pancreatic β-cells, (ii) affect and modulate the allostery of the insulin hexamer, (iii) may provide folding stability and protection for the insulin hexamer, and consequently (iv), may have a direct role in the modulation of insulin release from β-cells in vivo. Moreover, these ligands can act in a synergistic, heterotropic fashion, pointing to a co-operative, complex, and dynamic nature of the interactions within the insulin hexamer.

Therefore, the importance of the neurotransmitters for insulin beta-cell biology can be either cumulative or even independent. n class="Chemical">Serotonin/dopamine are involved in the regulation of exocytosis of the granules, whereas their impact on insulin storage can be seen as a serendipitous side effect of the evolution of the hormone within the β-cell environment and its optimization for a particular, species-specific, physiological profile.

It is also tempting to consider serotonin/n class="Chemical">dopamine-human insulin interactions as molecular fingerprints of insulin evolution in the animal kingdom, as insulin-like hormones are expressed in neurons in many invertebrates (54). The newly discovered (and still puzzling) role of insulin in the central nervous system and the emergence of neurodegeneration-linked so-called type III diabetes thus places our findings in a much wider physiological context (55, 56).

In summary, we show here that the insulin hexamer can act as a macromolecular sponge, keen to bind and assimilate a variety of physiological ligands and ions. Tn class="Chemical">his suggests that clinical exploration and use of insulin-hexamer-stabilizing ligands, driven by steering of the pharmacokinetics of this hormone, may mirror its physiological canonical properties, rooted also in insulin evolution. Although the present work is aimed at shedding new light on the physiological storage state of this hormone, it may also open ways to novel approaches in medical formulations of insulin, toward its more desired closer to in vivo profile during subcutaneous injections administered in diabetes.

Experimental procedures

Molecular dynamics simulations

All MD simulations were performed using the AMBER 14 program (57, 58). In all simulations the AMBER ff03 protein force field (59) and SPC/E n class="Chemical">water model were used (60). Production simulations were performed in the isothermal-isobaric ensemble at ambient conditions of T = 300 K and p = 1 atm using the Berendsen barostat and thermostat (61). The only exception was thermodynamic integration calculations, where the temperature was controlled by a Langevin thermostat with a reference temperature of T = 300 K and a collision frequency of 5 ps−1 in order to avoid problems of non-ergodicity when the ligand is fully decoupled from its environment. 3D periodic boundary conditions were applied with a non-bonded interaction cutoff of 9 Å. The long-range electrostatic interactions were accounted for using the particle mesh Ewald method (62) using a cubic spline interpolation. The density of the charge grid was 64 × 64 × 64, and the direct sum threshold was 10−5. Van der Waals interactions beyond the cutoff were treated using the continuum model correction for energy and pressure. All bonds containing hydrogen atoms were constrained using the SHAKE algorithm (63). A time step of 2 fs was employed.

Mildly acidic conditions (pH ∼ 5.5) of insulin storage granules would generally tend to favor a fully protonated side chain of n class="Chemical">HisB5 (64). However, due to zinc coordination the Nϵ2 atom of HisB10 is actually deprotonated (65, 66). In addition, the six GluB13 side chains in the middle region of the insulin hexamer were considered as deprotonated. The overall charge of the protein is, therefore, −6. Initial simulations indicated diffusion of Na+ cation/cations into the middle region. At least one Na+ cation was always present in the middle region of the hexamer being placed to the middle region from the start of the simulation.

Substitution of the site I-bound phenol by the neurotransmitters n class="Chemical">serotonin or dopamine was also investigated. The R6-state insulin complex with six phenols in sites I (referred to here as InsPheR6) was used as a reference structure (PDB ID 1AIY) due to a higher stability of this conformer and homogenous saturation of insulin hexamer with the same ligand (67). The starting structures for simulations were obtained by exchanging phenol molecules in site I by either dopamine or serotonin molecules with two initial conformations for each ligand. The ligands were built from the phenol core (1 oxygen atom and 6 carbon atoms), the structure of which was available from X-ray structure). Subsequent restrained minimization eliminated any potential steric clashes. Starting geometries of phenolic ligands inside the phenolic pockets are depicted in Fig. 2.

Phenol, n class="Chemical">serotonin, and dopamine were assigned parameters from ff03 AMBER force field using the ANTECHAMBER package (68). Both serotonin and dopamine have a charge of +1 e at pH 5.5. Partial atomic charges were obtained by the RESP (Restrained Electrostatic Potential Fit) (83) method, calculated at the HF/6–31G+ level, and they are listed in supplemental Table S2, with the explanation in the supplemental Figure S5. These calculations were performed using the Gaussian 09 package (69).

The B10-Zn interaction potential had to be re-parametrized in order to account (at least partially) for the electronic polarization and charge transfer effects and thus to reproduce experimental data. The charges on zinc ions and on HisB10 were modified according to the results from ab initio calculations on small model systems to +1.5 e and +0.1677 e, respectively. These were obtained using the RESP/NPA (Natural Population Analysis) (84) analysis employing the B3LYP/aug-cc-pvtz level of theory (85, 86). In a similar spirit, to account for electronic polarization effects in a mean-field way (70–73), we rescaled bulk Na+ and Cl− (including the Na+ cation located in the middle of the hexamer) ionic charges by a factor of 0.75.

Each of the liganded insulin R6 hexamers was immersed into a unit cell containing 9000 n class="Gene">SPC/E water molecules, with Na+/Cl− ions added to acquire overall electroneutrality with no excess of salt present. After preparation, the energy of each of the systems was minimized using 5000 steps of the steepest descent method, where the protein and phenolic ligands were restrained with a harmonic potential. The systems were then subjected to 200 ps of isothermal-isochoric molecular dynamics, where the temperature was slowly raised from 10 to 300 K. This was followed by 1.2 ns of an isothermal-isobaric equilibration, which led to an equilibrated cell size of ∼69 × 68 × 66 Å3. Systems were assessed as equilibrated by monitoring temperature, cell size, density, and root mean square displacement of the protein. After equilibration, the production runs were propagated for 600 ns.

In addition to direct MD simulations free energy calculations employing the thermodynamic integration method were performed to determine apparent binding constants K of n class="Chemical">phenol and the two neurotransmitters to the R6 hexamer. The free energy cycles and further details of these simulations are presented in the supplemental Figs. S6 and S7, and Table S3. Similar free energy calculations for evaluation of the K values to ligand binding sites III were performed using the umbrella sampling methods (for more details see the supplemental information, and Table S4). Here, to stabilize the hexamer R6 state, site I was filled with phenol before the search for dopamine site III binding mode, as this neurotransmitter did not bind to site I in MD simulations.

X-ray crystallography

Crystallizations of all insulin complexes reported here were performed with the in-house n class="Gene">insulin crystallization screens that cover most of the previously reported crystal growth parameters. Crystallization conditions, data collection, refinement, and models statistics as well as PDB codes are provided in supplemental Table S1. All crystals were directly flash-cooled in liquid N2. X-ray data were collected at 100 K and processed by xia2 (74), and model building and refinement (F>0σF) were performed by COOT (75) and the CCP4 suite of programs (76). Crystal structures were solved by Molrep (77) with the B1-B6-truncated insulin hexamer, hexamer-derived dimer, and insulin monomer as a model (based on PDB ID 1MSO; Ref. 78) and refined by Refmac 5.8 (79). Examples of electron density maps are shown in supplemental Figs. S3 and S4. The figures were made using CCP4mg (80). For structural comparisons, all insulin hexamer structures were superimposed in COOT by the SSM fit option. For comparison of the site I ligand-binding modes, the relevant dimers were superimposed by Cα matching of the B9–19 helix for one of the monomers by the LSQ option in COOT.

Determination of ligand K values by a solution 4H3N assay

Solution studies of the interactions of phenol, n class="Chemical">serotonin, dopamine, and arginine with porcine insulin were performed following the protocols by Huang et al. (81), Bloom et al. (82), and Huus et al. (46). An anionic ligand 4H3N was used as a sensitive chromophoric probe to determine ligand binding curves by monitoring UV-visible absorption. 4H3N binds only in the insulin R-state (in the vicinity of HisB10), which red-shifts its absorption spectrum. The samples were dissolved in 10 mm potassium phosphate buffer at pH 7.4. In all measurements porcine hexameric Zn2+-insulin (two atoms of zinc per one hexamer) and 4H3N were used at 0.6 mm and 0.225 mm concentrations, respectively. Insulin was preincubated with 4H3N in a total volume 0.6 ml for 20 min at room temperature, after which a particular ligand (phenol, serotonin, dopamine, or arginine) was added in a minimum volume (0.6–3 μl) of 10 mm potassium phosphate buffer. After 15 min of reverse transcription equilibration, the UV-visible absorbance spectra were collected on a Perkin, Lambda 25 UV-visible Spectrometer. The spectra were taken in the range of 300–550 nm, and the difference spectra at λmax of 444 nm were obtained by subtracting the absorbance of ligand-free 4H3N from the absorbance of ligand-bound 4H3N. The binding curves were analyzed using a method of non-linear regression and a fitting program considering one-site specific binding with a Hill slope using GraphPad Prism 5.0. The final K, Bmax, and h values with standard errors were calculated from at least three independently determined binding curves for each system.

Author contributions

V. P. carried out the MD simulations. C. M. V., M. K., and T. R. G, carried out the crystallographic experiments. K. K. and P. H. carried out the solution 4H3N assays. J. P. T. collected the crystallographic data. A. M. B., P. J., J. J., and L. Z. conceived the study, designed the experiments, and analyzed the data. A. M. B., P. J., and J. J. wrote the paper. All authors discussed the results and commented on the manuscript.