Structure of human insulin monomer in water/acetonitrile solution

Here we present evidence that in water/acetonitrile solvent detailed structural and dynamic information can be obtained for important proteins that are naturally present as oligomers under native conditions. An NMR-derived human insulin monomer structure in H₂O/CD₃CN, 65/35 vol%, pH 3.6 is presented and compared with the available X-ray structure of a monomer that forms part of a hexamer (Acta Crystallogr. 2003 Sec. D59, 474) and with NMR structures in water and organic cosolvent. Detailed analysis using PFGSE NMR, temperature-dependent NMR, dilution experiments and CSI proves that the structure is monomeric in the concentration and temperature ranges 0.1–3 mM and 10–30°C, respectively. The presence of long-range interstrand NOEs, as found in the crystal structure of the monomer, provides the evidence for conservation of the tertiary structure. Starting from structures calculated by the program CYANA, two different molecular dynamics simulated annealing refinement protocols were applied, either using the program AMBER in vacuum (AMBER_VC), or including a generalized Born solvent model (AMBER_GB). Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Wojciech Bocian contributed equally to this work.     

NMR structure of a parallel-stranded DNA duplex at atomic resolution

DNA dodecamers have been designed with two cytosines on each end and intervening A and T stretches, such that the oligomers have fully complementary A:T base pairs when aligned in the parallel orientation. Spectroscopic (UV, CD and IR), NMR and molecular dynamics studies have shown that oligomers having the sequences d(CCATAATTTACC) and d(CCTATTAAATCC) form a parallel-stranded duplex when dissolved at 1:1 stoichiometry in aqueous solution. This is due to the C:C⁺ clamps on either end and extensive mismatches in the antiparallel orientation. The structure is stable at neutral and acidic pH. At higher temperatures, the duplex melts into single strands in a highly cooperative fashion. All adenine, cytosine and thymine nucleotides adopt the anti conformation with respect to the glycosidic bond. The A:T base pairs form reverse Watson–Crick base pairs. The duplex shows base stacking and NOEs between the base protons T(H6)/A(H8) and the sugar protons (H1′/H2′/H2″) of the preceding nucleotide, as has been observed in antiparallel duplexes. However, no NOEs are observed between base protons H2/H6/H8 of sequential nucleotides, though such NOEs are observed between T(CH₃) and A(H8). A three-dimensional structure of the parallel-stranded duplex at atomic resolution has been obtained using molecular dynamics simulations under NMR constraints. The simulated structures have torsional angles very similar to those found in B-DNA duplexes, but the base stacking and helicoid parameters are significantly different.     

Solution structure of the coiled-coil trimerization domain from lung surfactant protein D

Surfactant protein D (SP-D) is one of four known protein components of the pulmonary surfactant lining the lung alveoli. It is involved in immune and allergic responses. SP-D occurs as a tetramer of trimers. Trimerization is thought to be initiated by a coiled coil domain. We have determined the solution structure of a 64-residue peptide encompassing the coiled coil domain of human SP-D. As predicted, the domain forms a triple-helical parallel coiled coil. As with all symmetric oligomers, the structure calculation was complicated by the symmetry degeneracy in the NMR spectra. We used the symmetry-ADR (ambiguous distance restraint) structure calculation method to solve the structure. The results demonstrate that the leucine zipper region of SP-D is an autonomously folded domain. The structure is very similar to the independently determined X-ray crystal structure, differing mainly at a single residue, Tyr248. This residue is completely symmetric in the solution structure, and markedly asymmetric in the crystalline phase. This difference may be functionally important, as it affects the orientation of the antigenic surface presented by SP-D.     

Zinc–Ligand Interactions Modulate Assembly and Stability of the Insulin Hexamer – A Review

Zinc and calcium ions play important roles in the biosynthesis and storage of insulin. Insulin biosynthesis occurs within the β-cells of the pancreas via preproinsulin and proinsulin precursors. In the golgi apparatus, proinsulin is sequestered within Zn²⁺- and Ca²⁺-rich storage/secretory vesicles and assembled into a Zn²⁺ and Ca²⁺ containing hexameric species, (Zn²⁺)₂(Ca²⁺)(Proin)₆. In the vesicle, (Zn²⁺)₂(Ca²⁺)(Proin)₆ is converted to the insulin hexamer, (Zn²⁺)₂(Ca²⁺)(In)₆, by excision of the C-peptide through the action of proteolytic enzymes. The conversion of (Zn²⁺)₂(Ca²⁺)(Proin)₆to (Zn²⁺)₂(Ca²⁺)(In)₆ significantly lowers the solubility of the hexamer, causing crystallization within the vesicle. The (Zn²⁺)₂(Ca²⁺)(In)₆ hexamer is an allosteric protein that undergoes ligand-mediated interconversion among three global conformation states designated T₆, T₃R₃ and R₆. Two classes of allosteric sites have been identified; hydrophobic pockets (3 in T₃R₃ and 6 in R₆) that bind phenolic ligands, and anion sites (1 in T₃R₃ and 2 in R₆) that bind monovalent anions. The allosteric states differ widely with respect to the physical and chemical stability of the insulin subunits. Fusion of the vesicle with the plasma membrane results in the expulsion of the insulin crystals into the intercellular fluid. Dissolution of the crystals, dissociation of the hexamers to monomer and transport of monomers to the liver and other tissues then occurs via the blood stream. Insulin action then follows binding to the insulin receptors. The role of Zn²⁺ in the assembly, structure, allosteric properties, and dynamic behavior of the insulin hexamer will be discussed in relation to biological function.     

The solution structure of cytochrome c 6 from the green alga Monoraphidium braunii

 The three-dimensional structure in solution of the reduced form of cytochrome c ₆ from the green alga Monoraphidium braunii has been solved through NMR data. Cytochrome c ₆ acts as a small mono-heme electron carrier protein between the two membrane-embedded complexes cytochrome f and photosystem I. The structure was determined using 1278 relevant interproton NOEs out of 1776 assigned NOEs with distance geometry (DG) calculations which included 36 stereospecific assignments and 20 experimentally found angle constraints. The family of structures obtained from the DG calculations was subjected to energy minimization and molecular dynamics simulation using previously defined force field parameters for the heme and its ligands. In all stages of the calculations, the solution structure is well defined and similar to the now available X-ray structure. Received: 18 January 1996 / Accepted: 19 April 1996     

NMR structure of biosynthetic engineered human insulin monomer B31(Lys)-B32(Arg) in water/acetonitrile solution. Comparison with the solution structure of native human insulin monomer

A solution NMR-derived structure of a new long -acting, B31(Lys)-B32(Arg) (LysArg), engineered human insulin monomer, in H(2)O/CD(3)CN, 65/35 vol %, pH 3.6, is presented and compared with the available X-ray structure of a monomer that forms part of a hexamer (Smith, et al., Acta Crystallogr D 2003, 59, 474) and with NMR structure of human insulin in the same solvent (Bocian, et al., J Biomol NMR 2008, 40, 55-64). Detailed analysis using PFGSE NMR (Pulsed Field Gradient Spin Echo NMR) in dilution experiments and CSI analysis prove that the structure is monomeric in the concentration range 0.1-3 mM. The presence of long-range interstrand NOEs in a studied structure, relevant to the distances found in the crystal structure of the monomer, provides the evidence for conservation of the tertiary structure. Therefore the results suggest that this solvent system is a suitable medium for studying the native conformation of the protein, especially in situations (as found for insulins) in which extensive aggregation renders structure elucidations in water difficult or impossible. Starting from the structures calculated by the program CYANA, two different molecular dynamics (MD) simulated annealing refinement protocols were applied, either using the program AMBER in vacuum (AMBER_VC), or including a generalized Born solvent model (AMBER_GB). Here we present another independent evidence to the one presented recently by us (Bocian et al., J Biomol NMR 2008, 40, 55-64), that in water/acetonitrile solvent detailed structural and dynamic information can be obtained for important proteins that are naturally present as oligomers under native conditions. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 820-830, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.     

NMR triple-quantum filtered relaxation analysis of O-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 ¹⁷O-water was used to study the properties of water in insulin solutions at different Zn²⁺ 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, τ_sb, is 10 ns, the apparent number of strongly bound water molecules was 3 to 4 per insulin monomer. Addition of Zn²⁺ equivalent to 2 g-atoms per hexamer did not produce substantial increases in the overall ¹⁷O-water TTQF signal intensity and apparent fraction of bound water. The dramatic enhancement of the TTQF signals observed for samples with a Zn²⁺/hexamer ratio greater than 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.     

Three-dimensional structure of the weakly associated protein homodimer SeR13 using RDCs and paramagnetic surface mapping

The traditional NMR‐based method for determining oligomeric protein structure usually involves distinguishing and assigning intra‐ and intersubunit NOEs. This task becomes challenging when determining symmetric homo‐dimer structures because NOE cross‐peaks from a given pair of protons occur at the same position whether intra‐ or intersubunit in origin. While there are isotope‐filtering strategies for distinguishing intra from intermolecular NOE interactions in these cases, they are laborious and often prove ineffectual in cases of weak dimers, where observation of intermolecular NOEs is rare. Here, we present an efficient procedure for weak dimer structure determination based on residual dipolar couplings (RDCs), chemical shift changes upon dilution, and paramagnetic surface perturbations. This procedure is applied to the Northeast Structural Genomics Consortium protein target, SeR13, a negatively charged Staphylococcus epidermidis dimeric protein (K_d 3.4 ± 1.4 mM) composed of 86 amino acids. A structure determination for the monomeric form using traditional NMR methods is presented, followed by a dimer structure determination using docking under orientation constraints from RDCs data, and scoring under residue pair potentials and shape‐based predictions of RDCs. Validation using paramagnetic surface perturbation and chemical shift perturbation data acquired on sample dilution is also presented. The general utility of the dimer structure determination procedure and the possible relevance of SeR13 dimer formation are discussed.     

Observation of intersubunit NOEs in a dimeric P22 Mnt repressor mutant by a time-shared [15N, 13C] double half-filter technique

Summary A time-shared [¹⁵N, ¹³C] half-filter technique is presented, which can be used to study proton-proton NOEs between biomolecules. The filter is demonstrated in a 2D [¹⁵N, ¹³C] double filtered NOESY experiment of a dimeric Mnt repressor mutant consisting of completely [¹⁵N, ¹³C] labeled monomer and unlabeled monomer. The benefit of this combined [¹⁵N, ¹³C] half-filter is that a single NMR experiment can be designed that yields all NOE interactions between labeled and unlabeled protons ((¹³C, ¹⁴N/¹²C), (¹⁵N, ¹⁴N/¹²C), (¹²C, ¹⁵N/¹³C) and (¹⁴N, ¹⁵N/¹³C)) in the protein, where conventional half-filters would require at least three separate NMR experiments to obtain these NOEs. The intermonomer NOEs of the Mnt mutant confirmed the secondary structure of the DNA-binding domain as an antiparallel ribbon, formed from an N-terminal segment contributed by each monomer. Moreover, several intersubunit NOEs were characterized in the C-terminal part of the Mnt mutant for which no structural data is available yet.     

Backbone dynamics of the human CC-chemokine eotaxin

Eotaxin is a CC chemokine with potent chemoattractant activity towards eosinophils. ¹⁵N NMR relaxation data have been used to characterize the backbone dynamics of recombinant human eotaxin. ¹⁵N longitudinal (R₁) and transverse (R₂) auto relaxation rates, heteronuclear ¹H-¹⁵N steady-state NOEs, and transverse cross-relaxation rates (_xy) were obtained at 30 °C for all resolved backbone secondary amide groups using ¹ H-detected two-dimensional NMR experiments. Ratios of transverse auto and cross relaxation rates were used to identify NH groups influenced by slow conformational rearrangement. Relaxation data were fit to the extended model free dynamics formalism, yielding parameters describing axially symmetric molecular rotational diffusion and the internal dynamics of each NH group. The molecular rotational correlation time (_m) is 5.09±0.02 ns, indicating that eotaxin exists predominantly as a monomer under the conditions of the NMR study. The ratio of diffusion rates about unique and perpendicular axes (D/D) is 0.81±0.02. Residues with large amplitudes of subnanosecond motion are clustered in the N-terminal region (residues 1–19), the C-terminus (residues 68–73) and the loop connecting the first two -strands (residues 30–37). N-terminal flexibility appears to be conserved throughout the chemokine family and may have implications for the mechanism of chemokine receptor activation. Residues exhibiting significant dynamics on the microsecond–millisecond time scale are located close to the two conserved disulfide bonds, suggesting that these motions may be coupled to disulfide bond isomerization.     

Structural signatures of the complex formed between 3-nitro-4-hydroxybenzoate and the Zn(II)-substituted R(6) insulin hexamer

3-Nitro-4-hydroxybenzoate (3N4H) is a probe of the structure and dynamics of the metal-centered His B10 assembly sites of the insulin hexamer. Each His B10 site consists of a approximately 12 A-long cavity situated on the threefold symmetry axis. These sites play an important role in the storage and release of insulin in vivo. The allosteric behavior of the insulin hexamer is modulated by ligand binding to the His B10 zinc sites and to the phenolic pockets. Binding to these sites drives transitions among three allosteric states, designated T(6), T(3)R(3), and R(6). Although a wide variety of mono anions bind to the His B10 zinc sites of R(3), X-ray structures of ligands complexed to this site exist only for H(2)O, Cl(-), and SCN(-). This work combines one- and two-dimensional (1)H NMR and UV-Vis absorbance studies of the structure and dynamics of the 3N4H complex, which establish the following: (1). relative to the NMR time scale, 3N4H exchange between free and bound states is slow, while flipping among three equivalent orientations about the site threefold axis is fast; (2). binding of 3N4H perturbs resonances within the His B10 zinc site and generates NOEs between ligand resonances and the insulin C-alpha and side chain resonances of ValB2, AsnB3, LeuB6, and CysB7; and (3).3N4H exchange for other ligands is limited by a protein conformational transition. These results are consistent with coordination of the 3N4H carboxylate to the His B10 zinc ion and van der Waals interactions with Val B2, Asn B3, Leu B6, and Cys A7.     

Metabolomics for the Quality Assessment of Lycium chinense Fruits

Aibellin is a 20-residue peptide antibiotic that has been isolated from the fungus Verticimonosporium ellipticum. Sequence-specific assignment of the ¹H- and ¹³C-NMR signals of aibellin in a methanol solution was achieved by using the two-dimensional NMR technique. Furthermore, its secondary structure was characterized by circular dichroism (CD) and NOESY spectra. The observed NOEs, ³J_NHCαH coupling constants and amide hydrogen–deuterium (H–D) exchange rates show that the peptide consisted of two α-helices and a bent structure around a Pro-14 residue.     

Conformational sampling of bioactive conformers: a low-temperature NMR study of 15N-Leu–enkephalin

Conformational studies of enkephalins are hampered by their high flexibility which leads to mixtures of quasi-isoenergetic conformers in solution and makes NOEs very difficult to detect in NMR spectra. In order to improve the quality of the NMR data, Leu–enkephalin was synthesized with ¹⁵N-labelled uniformly on all amide nitrogens and examined in a viscous solvent medium at low temperature. HMQC NOESY spectra of the labelled Leu–enkephalin in a DMSO_d6/H₂O mixture at 275 K do show numerous NOEs, but these are not consistent with a single conformer and are only sufficient to describe the conformational state as a mixture of several conformers. Here a different approach to the structure–activity relationships of enkephalins is presented: it is possible to analyse the NMR data in terms of limiting canonical structures (i.e. β- and γ-turns) and finally to select only those consistent with the requirements of δ selective agonists and antagonists. This strategy results in the prediction of a family of conformers that may be useful in the design of new δ selective opioid peptides. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Rotational diffusion tensor of nucleic acids from 13C NMR relaxation

Rotational diffusion properties have been derived for the DNA dodecamer d(CGCGAATTCGCG)₂ from ¹³C R₁ and R₁ measurements on the C₁, C₃, and C₄ carbons in samples uniformly enriched in ¹³C. The narrow range of C-H bond vector orientations relative to the DNA axis make the analysis particularly sensitive to small structural deviations. As a result, the R₁/R₁ ratios are found to fit poorly to the crystal structures of this dodecamer, but well to a recent solution NMR structure, determined in liquid crystalline media, even though globally the structures are quite similar. A fit of the R₁/R₁ ratios to the solution structure is optimal for an axially symmetric rotational diffusion model, with a diffusion anisotropy, D_||/D, of 2.1±0.4, and an overall rotational correlation time, (2D_||+4D)^–1, of 3.35 ns at 35 °C in D₂O, in excellent agreement with values obtained from hydrodynamic modeling.     

Biophysical Optimization of a Therapeutic Protein by Nonstandard Mutagenesis: STUDIES OF AN IODO-INSULIN DERIVATIVE*

Insulin provides a model for the therapeutic application of protein engineering. A paradigm in molecular pharmacology was defined by design of rapid-acting insulin analogs for the prandial control of glycemia. Such analogs, a cornerstone of current diabetes regimens, exhibit accelerated subcutaneous absorption due to more rapid disassembly of oligomeric species relative to wild-type insulin. This strategy is limited by a molecular trade-off between accelerated disassembly and enhanced susceptibility to degradation. Here, we demonstrate that this trade-off may be circumvented by nonstandard mutagenesis. Our studies employed Lys^B28, Pro^B29-insulin (“lispro”) as a model prandial analog that is less thermodynamically stable and more susceptible to fibrillation than is wild-type insulin. We have discovered that substitution of an invariant tyrosine adjoining the engineered sites in lispro (Tyr^B26) by 3-iodo-Tyr (i) augments its thermodynamic stability (ΔΔG_u 0.5 ±0.2 kcal/mol), (ii) delays onset of fibrillation (lag time on gentle agitation at 37 °C was prolonged by 4-fold), (iii) enhances affinity for the insulin receptor (1.5 ± 0.1-fold), and (iv) preserves biological activity in a rat model of diabetes mellitus. ¹H NMR studies suggest that the bulky iodo-substituent packs within a nonpolar interchain crevice. Remarkably, the 3-iodo-Tyr^B26 modification stabilizes an oligomeric form of insulin pertinent to pharmaceutical formulation (the R₆ zinc hexamer) but preserves rapid disassembly of the oligomeric form pertinent to subcutaneous absorption (T₆ hexamer). By exploiting this allosteric switch, 3-iodo-Tyr^B26-lispro thus illustrates how a nonstandard amino acid substitution can mitigate the unfavorable biophysical properties of an engineered protein while retaining its advantages.     

RDC-assisted modeling of symmetric protein homo-oligomers

Protein oligomerization serves an important function in biological processes, yet solving structures of protein oligomers has always been a challenge. For solution NMR, the challenge arises both from the increased size of these systems and, in the case of homo-oligomers, from ambiguities in assignment of intra- as opposed to intersubunit NOEs. In this study, we present a residual dipolar coupling (RDC)-assisted method for constructing models of homo-oligomers with purely rotational symmetry. Utilizing the fact that one of the principal axes of the tensor describing the alignment needed for RDC measurement is always parallel to the oligomer symmetry axis, it is possible to greatly restrict possible models for the oligomer. Here, it is shown that, if the monomer structure is known, all allowed dimer models can be constructed using a grid search algorithm and evaluated based on RDC simulations and the quality of the interface between the subunits. Using the Bacillus subtilis protein YkuJ as an example, it is shown that the evaluation criteria based on just two sets of NH RDCs are very selective and can unambiguously produce a model in good agreement with an existing X-ray structure of YkuJ.     

Sequence-specific 1H NMR assignment and secondary structure of the Arc repressor of bacteriophage P22, as determined by two-dimensional 1H NMR spectroscopy

The Arc repressor of bacteriophage P22 is a DNA binding protein that does not belong to any of the known classes of such proteins. We have undertaken a 1H NMR study of the protein with the aim of elucidating its three-dimensional structure in solution and its mode of binding of operator DNA. Here we present the 1H nuclear magnetic resonance (NMR) assignments of all backbone protons and most of the side-chain protons of Arc repressor. Elements of secondary structure have been identified on the basis of networks of characteristic sequential and medium-range nuclear Overhauser enhancements (NOEs). Two alpha-helical regions have been found in the peptide regions 16-29 and 35-45. The ends of the helices could not yet be firmly established and could extend to residue 31 for the first helix and to residue 49 for the second. Immediately before the first helix, between residues 8 and 14, a region is present with beta-sheet characteristics dominated by a close proximity of the alpha-protons of residues 9 and 13. Because of the dimeric nature of the protein there are still two possible ways in which the NOEs in the beta-sheet region can be interpreted. If the NOEs are intramonomer, this requires a tight turn involving residues 10-12. Alternatively, if the NOEs are intermonomer, then and antiparallel beta-sheet would be implicated comprising two strands of different Arc monomers. While the data presently do not allow an unambiguous choice between these two possibilities, some evidence is discussed that favors the latter (beta-sheet between monomers).(ABSTRACT TRUNCATED AT 250 WORDS)     

Crystal structure of the N-terminal domain of linker L(R) and the assembly of cyanobacterial phycobilisome rods

Phycobilisomes are light‐harvesting supramolecular complexes in cyanobacteria and red algae. Linkers play a pivotal role in the assembly and energy transfer modulation of phycobilisomes. However, how linkers function remains unclear due to the lack of structural and biochemical studies of linkers, especially the N‐terminal domain of L_R (pfam00427). Here, we report the crystal structure of the pfam00427 domain of the linker L_R³⁰ from Synechocystis sp. PCC 6803 at 1.9 Å. The pfam00427 presents as a previously uncharacterized point symmetric six α‐helix bundle. To elucidate the binding style of pfam00427 in the C‐phycocyanin (C‐PC) (αβ)₆ hexamer, we fixed pfam00427 computationally into the C‐PC (αβ)₆ inner cavity using the program AutoDock. Combined with a conserved ‘C‐PC binding patch’ on pfam00427 identified, we arrived at a model for the pfam00427–C‐PC (αβ)₆ complex. This model was further optimized and evaluated as a reasonable result by a molecular dynamics simulation. In the resulting model, the pfam00427 domain is stably positioned in the central hole of the C‐PC trimer. Moreover, the L_RT (pfam01383) was docked into our pfam00427–C‐PC model to generate a complete phycobilisome rod in which the linkers join individual biliprotein hexamers.     

Solution structure of T4moC,the Rieske ferredoxin component of the toluene 4-monooxygenase complex

Toluene 4-monooxygenase, a four-protein complex from Pseudomonas mendocina KR1, catalyzes the NADH- and O₂-dependent hydroxylation of toluene to form p-cresol. The solution structure of the 112-amino-acid Rieske ferredoxin component, T4moC, was determined from 2D and 3D ¹H, ¹³C, and ¹⁵N NMR data. The structural model was refined through simulated annealing by molecular dynamics in torsion angle space with input from 1650 experimental restraints, including 1264 inter-proton distance restraints obtained from NOEs, 247 non-redundant intra-residue NOEs, 26 hydrogen bond restraints, and 113 dihedral angle (, ) restraints. The 20 calculated conformers that best satisfied the input restraints were submitted to refinement in explicit solvent to improve the stereochemical quality. With exclusion of ill-defined N- and C-terminal segments (Ser2; His111–Ser112) and residues near to the [2Fe-2S] cluster, the atomic root mean square deviation for the 20 conformers with respect to the mean coordinates was 1.09 Å for the backbone and 1.60 Å for all non-hydrogen atoms. The T4moC structure consists of 10 -strands arranged in the three anti-parallel -sheet topology observed in all Rieske [2Fe-2S] domain proteins. The S of Cys45 and Cys64 and the N¹ of His47 and His67 provide the ligands to the [2Fe-2S] cluster of T4moC. ¹H–¹⁵N HSQC measurements show that both His47-N² and His67-N² are protonated at the pH of the NMR experiments. Comparisons are made between the present NMR structure, previous paramagnetic NMR studies of T4moC, and the X-ray structures of other members of the Rieske protein family.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00775-004-0594-4     

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 ΔF 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.