Dual role of Lys206-Lys296 interaction in human transferrin N-lobe: iron-release trigger and anion-binding site.

The unique structural feature of the dilysine (Lys206-Lys296) pair in the transferrin N-lobe (hTF/2N) has been postulated to serve a special function in the release of iron from the protein. These two lysines, which are located in opposite domains, hydrogen bond to each other in the iron-containing hTF/2N at neutral pH but are far apart in the apo-form of the protein. It has been proposed that charge repulsion resulting from the protonation of the dilysines at lower pH may be the trigger to open the cleft and facilitate iron release. The fact that the dilysine pair is positively charged and resides in a location close to the metal-binding center has also led to the suggestion that the dilysine pair is an anion-binding site for chelators. The present report provides comprehensive evidence to confirm that the dilysine pair plays this dual role in modulating release of iron. When either of the lysines is mutated to glutamate or glutamine or when both are mutated to glutamate, release of iron is much slower compared to the wild-type protein. This is due to the fact that the driving force for cleft opening is absent in the mutants or is converted to a lock-like interaction (in the case of the K206E and K296E mutants). Direct titration of the apo-proteins with anions as well as anion-dependent iron release studies show that the dilysine pair is part of an active anion-binding site which exists with the Lys296-Tyr188 interaction as a core. At this site, Lys296 serves as the primary anion-binding residue and Tyr188 is the main reporter for electronic spectral change, with smaller contributions from Lys206, Tyr85, and Tyr95. In iron-loaded hTF/2N, anion binding becomes invisible as monitored by UV-vis difference spectra since the spectral reporters Tyr188 and Tyr95 are bound to iron. Our data strongly support the hypothesis that the apo-hTF/2N exists in equilibrium between the open and closed conformations, because only in the closed form is Lys296 in direct contact with Tyr188. The current findings bring together observations, ideas, and experimental data from a large number of previous studies and shed further light on the detailed mechanism of iron release from the transferrin N-lobe. In iron-containing hTF/2N, Lys296 may still function as a target to introduce an anion (or a chelator) near to the iron-binding center. When the pH is lowered, the protonation of carbonate (synergistic anion for metal binding) and then the dilysine pair form the driving force to loosen the cleft, exposing iron; the nearby anion (or chelator) then binds to the iron and releases it from the protein.     

Toxic and Physiological Metal Uptake and Release by Human Serum Transferrin

An atomistic understanding of metal transport in the human body is critical to anticipate the side effects of metal-based therapeutics and holds promise for new drugs and drug delivery designs. Human serum transferrin (hTF) is a central part of the transport processes because of its ubiquitous ferrying of physiological Fe(III) and other transition metals to tightly controlled parts of the body. There is an atomistic mechanism for the uptake process with Fe(III), but not for the release process, or for other metals. This study provides initial insight into these processes for a range of transition metals—Ti(IV), Co(III), Fe(III), Ga(III), Cr(III), Fe(II), Zn(II)—through fully atomistic, extensive quantum mechanical/discrete molecular dynamics sampling and provides, to our knowledge, a new technique we developed to calculate relative binding affinities between metal cations and the protein. It identifies protonation of Tyr188 as a trigger for metal release rather than protonation of Lys206 or Lys296. The study identifies the difficulty of metal release from hTF as potentially related to cytotoxicity. Simulations identify a few critical interactions that stabilize the metal binding site in a flexible, nuanced manner.     

A QM/MM study of the complexes formed by aluminum and iron with serum transferrin at neutral and acidic pH

Serum transferrin (sTf) transports iron in serum and internalizes in cells via receptor mediated endocytosis. Additionally, sTf has been identified as the predominant aluminum carrier in serum. Some questions remain unclear about the exact mechanism for the metal release or whether the aluminum and iron show the same binding mode during the entire process. In the present work, simulation techniques at quantum and atomic levels have been employed in order to gain access into a molecular level understanding of the metal-bound sTf complex, and to describe the binding of Al(III) and Fe(III) ions to sTf. First, hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations were carried out in order to analyze the dynamics of the aluminum-loaded complex, taking into account the different pH conditions in blood and into the cell. Moreover, the complexes formed by transferrin with Al(III) and Fe(III) were optimized with high level density functional theory (DFT)/MM methods. All these results indicate that the interaction mode of Al(III) and Fe(III) with sTf change upon different pH conditions, and that the coordination of Al(III) and Fe(III) is not equivalent during the metal intake, transport and release processes. Our results emphasize the importance of the pH on the metal binding and release mechanism and suggest that Al(III) can follow the iron pathway to get access into cells, although once there, it may show a different binding mode, leading to a different mechanism for its release.     

Mutation of the iron ligand His 249 to Glu in the N-lobe of human transferrin abolishes the dilysine "trigger" but does not significantly affect iron release

Serum transferrin is the major iron transport protein in humans. Its function depends on its ability to bind iron with very high affinity, yet to release this bound iron at the lower intracellular pH. Possible explanations for the release of iron from transferrin at low pH include protonation of a histidine ligand and the existence of a pH-sensitive "trigger" involving a hydrogen-bonded pair of lysines in the N-lobe of transferrin. We have determined the crystal structure of the His249Glu mutant of the N-lobe half-molecule of human transferrin and compared its iron-binding properties with those of the wild-type protein and other mutants. The crystal structure, determined at 2.4 A resolution (R-factor 19.8%, R(free) 29.4%), shows that Glu 249 is directly bound to iron, in place of the His ligand, and that a local movement of Lys 296 has broken the dilysine interaction. Despite the loss of this potentially pH-sensitive interaction, the H249E mutant is only slightly more acid-stable than wild-type and releases iron slightly faster. We conclude that the loss of the dilysine interaction does make the protein more acid stable but that this is counterbalanced by the replacement of a neutral ligand (His) by a negatively charged one (Glu), thus disrupting the electroneutrality of the binding site.     

Crystal structures and iron release properties of mutants (K206A and K296A) that abolish the dilysine interaction in the N-lobe of human transferrin

Human transferrin (Tf) is responsible for the binding and transport of iron in the bloodstream of vertebrates. Delivery of this bound iron to cells occurs by a process of receptor-mediated endocytosis during which Tf releases its iron at the reduced endosomal pH of approximately 5.6. Iron release from Tf involves a large conformational change in which the two domains that enclose the binding site in each lobe move apart. We have examined the role of two lysines, Lys206 and Lys296, that form a hydrogen-bonded pair close to the N-lobe binding site of human Tf and have been proposed to form a pH-sensitive trigger for iron release. We report high-resolution crystal structures for the K206A and K296A mutants of the N-lobe half-molecule of Tf, hTf/2N, and quantitative iron release data on these mutants and the double mutant K206A/K296A. The refined crystal structures (for K206A, R = 19.6% and R(free) = 23.7%; for K296A, R= 21.2% and R(free) = 29.5%) reveal a highly conserved hydrogen bonding network in the dilysine pair region that appears to be maintained even when individual hydrogen bonding groups change. The iron release data show that the mutants retain iron to a pH 1 unit lower than the pH limit of wild type hTf/2N, and release iron much more slowly as a result of the loss of the dilysine interaction. Added chloride ions are shown to accelerate iron release close to the pH at which iron is naturally lost and the closed structure becomes destabilized, and to retard it at higher pH.     

Crystal structure and iron-binding properties of the R210K mutant of the N-lobe of human lactoferrin: implications for iron release from transferrins

Lactoferrin (Lf) and serum transferrin (Tf) combine high-affinity iron binding with an ability to release this iron at reduced pH. Lf, however, retains iron to significantly lower pH than Tf, giving the two proteins distinct functional roles. In this paper, we compared the iron-release profiles for human Lf, Tf, and their N-lobe half-molecules Lf(N) and Tf(N) and showed that half of the difference in iron retention at low pH ( approximately 1.3 pH units) results from interlobe interactions in Lf. To probe factors intrinsic to the N-lobes, we further examined the specific role of two basic residues that are proposed to form a pH-sensitive dilysine trigger for iron release in the N-lobe of Tf [Dewan, J. C., Mikami, B., Hirose, M., and Sacchettini, J. C. (1993) Biochemistry 32, 11963-11968] by mutating Arg 210 to Lys in the N-lobe half-molecule Lf(N). The R210K mutant was expressed, purified, and crystallized, and its crystal structure was determined and refined at 2.0-A resolution to a final R factor (R(free)) of 19.8% (25.0%). The structure showed that Lys 210 and Lys 301 in R210K do not form a dilysine interaction like that between Lys 206 and Lys 296 in human Tf. The R210K mutant retained iron to lower pH than Tf(N), consistent with the absence of the dilysine interaction but released iron at approximately 0.7 pH units higher than Lf(N). We conclude that (i) the ability of Lf to retain iron to significantly lower pH than Tf is due equally to interlobe interactions and to the absence in Lfs of an interaction analogous to the dilysine pair in Tfs, even when two lysines are present at the corresponding sequence positions, and (ii) an appropriately positioned basic residue (Arg 210 in human Lf) modulates iron release by inhibiting protonation of the N-lobe iron ligands, specifically His 253.     

Protein-Protein Docking with Dynamic Residue Protonation States

Protein-protein interactions depend on a host of environmental factors. Local pH conditions influence the interactions through the protonation states of the ionizable residues that can change upon binding. In this work, we present a pH-sensitive docking approach, pHDock, that can sample side-chain protonation states of five ionizable residues (Asp, Glu, His, Tyr, Lys) on-the-fly during the docking simulation. pHDock produces successful local docking funnels in approximately half (79/161) the protein complexes, including 19 cases where standard RosettaDock fails. pHDock also performs better than the two control cases comprising docking at pH 7.0 or using fixed, predetermined protonation states. On average, the top-ranked pHDock structures have lower interface RMSDs and recover more native interface residue-residue contacts and hydrogen bonds compared to RosettaDock. Addition of backbone flexibility using a computationally-generated conformational ensemble further improves native contact and hydrogen bond recovery in the top-ranked structures. Although pHDock is designed to improve docking, it also successfully predicts a large pH-dependent binding affinity change in the Fc–FcRn complex, suggesting that it can be exploited to improve affinity predictions. The approaches in the study contribute to the goal of structural simulations of whole-cell protein-protein interactions including all the environmental factors, and they can be further expanded for pH-sensitive protein design.     

A computational study of the open and closed forms of the N-lobe human serum transferrin apoprotein

Human serum transferrin tightly binds ferric ions in the blood stream but is able to release them in cells by a process involving receptor-mediated endocytosis and decrease in pH. Iron binding and release are accompanied by a large conformation change. In this study, we investigate theoretically the open and closed forms of the N-lobe human serum transferrin apoprotein by performing pKa calculations and molecular dynamics and free-energy simulations. In agreement with the hypothesis based on the x-ray crystal structures, our calculations show that there is a shift in the pKa values of the lysines forming the dilysine trigger when the conformation changes. We argue, however, that simple electrostatic repulsion between the lysines is not sufficient to trigger domain opening and, instead, propose an alternative explanation for the dilysine-trigger effect. Analysis of the molecular dynamics and free-energy results indicate that the open form is more mobile than the closed form and is much more stable at pH 5.3, in large part due to entropic effects. Despite a lower free energy, the dynamics simulation of the open form shows that it is flexible enough to sample conformations that are consistent with iron binding.     

Mutations at the histidine 249 ligand profoundly alter the spectral and iron-binding properties of human serum transferrin N-lobe

Human serum transferrin is an iron-binding and -transport protein which carries iron from the blood stream into various cells. Iron is held in two deep clefts located in the N- and C-lobes by coordinating to four amino acid ligands, Asp 63, Tyr 95, Tyr 188, and His 249 (N-lobe numbering), and to two oxygens from carbonate. We have previously reported the effect on the iron-binding properties of the N-lobe following mutation of the ligands Asp 63, Tyr 95, and Tyr 188. Here we report the profound functional changes which result from mutating His 249 to Ala, Glu, or Gln. The results are consistent with studies done in lactoferrin which showed that the histidine ligand is critical for the stability of the iron-binding site [H. Nicholson, B. F. Anderson, T. Bland, S. C. Shewry, J. W. Tweedie, and E. N. Baker (1997) Biochemistry 36, 341-346]. In the mutant H249A, the histidine ligand is disabled, resulting in a dramatic reduction in the kinetic stability of the protein toward loss of iron. The H249E mutant releases iron three times faster than wild-type protein but shows significant changes in both EPR spectra and the binding of anion. This appears to be the net effect of the metal ligand substitution from a neutral histidine residue to a negative glutamate residue and the disruption of the "dilysine trigger" [MacGillivray, R. T. A., Bewley, M. C., Smith, C. A., He, Q.-Y., Mason, A. B., Woodworth, R. C., and Baker, E. N. (2000) Biochemistry 39, 1211-1216]. In the H249Q mutant, Gln 249 appears not to directly contact the iron, given the similarity in the spectroscopic properties and the lability of iron release of this mutant to the H249A mutant. Further evidence for this idea is provided by the preference of both the H249A and H249Q mutants for nitrilotriacetate rather than carbonate in binding iron, probably because NTA is able to provide a third ligation partner. An intermediate species has been identified during the kinetic interconversion between the NTA and carbonate complexes of the H249A mutant. Thus, mutation of the His 249 residue does not abolish iron binding to the transferrin N-lobe but leads to the appearance of novel iron-binding sites of varying structure and stability.     

Release of ADP from the catalytic subunit of protein kinase A: a molecular dynamics simulation study

Substrate phosphorylation by cAMP-dependent-protein kinase A (protein kinase A, PKA) has been studied extensively. Phosphoryl transfer was found to be fast, whereas ADP release was found to be the slow, rate-limiting step. There is also evidence that ADP release may be preceded by a partially rate-limiting conformational change. However, the atomic details of the conformational change and the mode of ADP release are difficult to obtain experimentally. In this work, we studied ADP release from PKA by carrying out molecular dynamics simulations with different pulling forces applied to the ligand. The detailed ADP release pathway and the associated conformational changes were analyzed. The ADP release process was found to involve a swinging motion with the phosphate of ADP anchored to the Gly-rich loop, so that the more buried adenine base and ribose ring came out before the phosphate. In contrast to the common belief that a hinge-bending motion was responsible for the opening of the ligand-binding cleft, our simulations showed that the small lobe exhibited a large amplitude "rocking" motion when the ligand came out. The largest conformational change of the protein was observed at about the first quarter time point along the release pathway. Two prominent intermediate states were observed in the release process.     

"Dilysine trigger" in transferrins probed by mutagenesis of lactoferrin: crystal structures of the R210G,R210E,and R210L mutants of human lactoferrin

The mammalian iron-binding proteins lactoferrin (Lf) and transferrin (Tf) bind iron very tightly, but reversibly. Despite homologous structures and essentially identical iron binding sites, Tf begins to release iron at pH 6.0, whereas Lf retains iron to pH approximately 3.5. This difference in iron retention gives the two proteins different biological roles. Two lysine residues, Lys 206 and Lys 296, which form a hydrogen-bonded dilysine pair in human Tf, have been shown to strongly influence iron release from the N-lobe. The equivalent residues in human Lf are Arg 210 and Lys 301, and we have here mutated Arg 210 in the N-lobe half-molecule of human lactoferrin, Lf(N), to probe its role in iron release. The Lf(N) mutants R210G, R210E, and R210L were expressed, purified, and crystallized, and their crystal structures were determined and refined at resolutions of 1.95 A (R210G), 2.2 A (R210E), and 2.0 A (R210L). The overall structures are very similar to that of wild-type Lf(N), but with small differences in domain orientations. In each of the mutants, however, Lys 301 (equivalent to Lys 296 in Tf) changes its conformation to fill the space occupied by Arg 210 Neta2 in wild-type Lf(N), interacting with the two tyrosine ligands Tyr 92 and Tyr 192. By comparison with other Lf and Tf structures, we conclude that Lys 301 (or Lys 296 in Tf) only occupies this site when residue 210 (206 in Tf) is nonpositive (neutral as in R210G and R210L or negative as in R210E). Thus, Lys 206 in the Tf dilysine pair is identified as having a depressed pK(a). Three specific sites are variably occupied by polar groups in the Lf mutants and other Lf and Tf proteins, and when coupled with iron-release data, these give new insights into the factors that most influence iron retention at low pH.     

Conformational transitions of uracil transporter UraA from Escherichia coli: a molecular simulation study

The Escherichia coli uracil/H + symporter UraA, known as the representative nucleobase/cation symporter 2(NCS2) protein, gets involved in several crucial physiological processes for most living organisms on Earth, such as the uptake of nucleobases and transport of vitamin C. Some experiments proposed a working model to explain proton-coupling and uracil transporting process of UraA on the basis of the crystal structure of NCS2 protein, but the details of conformational changes remained unknown. Thus, in order to make clear conformational changes caused by the protonation and deprotonation process of some conserved proton-coupled residues, the molecular dynamics simulation was used to study the conformation of UraA complexes in different protonation states. The results demonstrated that the protonation of residue Glu241 and Glu290 resulted in the whole conformational transition from the inward-open to the outward-open state. It can be concluded that Glu290 was crucial in a network of hydrogen-bonds in the middle of the core domain involving another essential residue, mainly including tyr288 in TM8, Tyr342, Ser338 in TM12, and the network of hydrogen-bonds was the key to maintain the stability of conformation. Protonation of Glu290 affects the stability of network of H-bond and changed the domains TM3 TM10 TM12. Thus, Glu290 may play a vital role as a ‘proton trigger’ that affects spatial structural of amino and residues near substrate binding side leading to an outward-open conformation transition.     

Rational design of a conformation-switchable Ca2+- and Tb3+-binding protein without the use of multiple coupled metal-binding sites

Ca2+, as a messenger of signal transduction, regulates numerous target molecules via Ca2+-induced conformational changes. Investigation into the determinants for Ca2+-induced conformational change is often impeded by cooperativity between multiple metal-binding sites or protein oligomerization in naturally occurring proteins. To dissect the relative contributions of key determinants for Ca2+-dependent conformational changes, we report the design of a single-site Ca2+-binding protein (CD2.trigger) created by altering charged residues at an electrostatically sensitive location on the surface of the host protein rat Cluster of Differentiation 2 (CD2).CD2.trigger binds to Tb3+ and Ca2+ with dissociation constants of 0.3 +/- 0.1 and 90 +/- 25 microM, respectively. This protein is largely unfolded in the absence of metal ions at physiological pH, but Tb3+ or Ca2+ binding results in folding of the native-like conformation. Neutralization of the charged coordination residues, either by mutation or protonation, similarly induces folding of the protein. The control of a major conformational change by a single Ca2+ ion, achieved on a protein designed without reliance on sequence similarity to known Ca2+-dependent proteins and coupled metal-binding sites, represents an important step in the design of trigger proteins.     

Insights into conformational regulation of PfMATE transporter from Pyrococcus furiosus induced by alternating protonation state of Asp41 residue: A molecular dynamics simulation study

BackgroundMultidrug and toxic compound extrusion (MATE) family transporters induce multiple-drug resistance (MDR) of bacterial pathogens and cancer cells, thus causing critical reductions in the therapeutic efficacies of antibiotics and anti-cancer drugs. Unfortunately, to date, the details and intrinsic reason about conformational regulation mechanism of MATE transporters remain elusive.MethodIn this work, molecular dynamics (MD) simulations were conducted to explore the conformational regulation mechanism of PfMATE transporter from Pyrococcus furiosus based on different protonation state of Asp41. Two (MD) simulation systems were investigated: a system with protonation of Asp41 and a system without protonation of Asp41, which were named by D184(H)D41(H) system and D184(H) system, respectively.Results and conclusionsFirstly, MD simulation results indicate that conformational changes mainly happen in extracellular regions of PfMATE protein. Further analysis reveals that PfMATE protein experiences different motion mode and forms different conformation based on different protonation state of Asp41. In the D184(H)D41(H) system, PfMATE experiences an opening motion and forms a more outward-open conformation. As for the D184(H) system, the protein has an anticlockwise rotational motion with the channel axis of protein and the more outward-open conformation does not appear. It can be inferred that protonation of Asp41 is essential for conformational regulation of PfMATE during transporting substrates.General significanceThese findings provide intrinsic information for understanding the conformational regulation mechanism of PfMATE and will be very meaningful to explore the MDR mechanism of PfMATE further.     

Probing the Periplasmic-Open State of Lactose Permease in Response to Sugar Binding and Proton Translocation

Based on the crystal structure of lactose permease (LacY) open to the cytoplasm, a hybrid molecular simulation approach with self-guided Langevin dynamics is used to describe conformational changes that lead to a periplasmic-open state. This hybrid approach consists of implicit (IM) and explicit (EX) membrane simulations and requires self-guided Langevin dynamics to enhance protein motions during the IM simulations. The pore radius of the lumen increases by 3.5 Å on the periplasmic side and decreases by 2.5 Å on the cytoplasmic side (relative to the crystal structure), suggesting a lumen that is fully open to the periplasm to allow for extracellular sugar transport and closed to the cytoplasm. Based on our simulations, the mechanism that triggers this conformational change to the periplasmic-open state is the protonation of Glu269 and binding of the disaccharide. Then, helix packing is destabilized by breaking of several side chains involved in hydrogen bonding (Asn245, Ser41, Glu374, Lys42, and Gln242). For the periplasmic-open conformations obtained from our simulations, helix-helix distances agree well with experimental measurements using double electron-electron resonance, fluorescence resonance energy transfer, and varying sized cross-linkers. The periplasmic-open conformations are also in compliance with various substrate accessibility/reactivity measurements that indicate an opening of the protein lumen on the periplasmic side on sugar binding. The comparison with these measurements suggests a possible incomplete closure of the cytoplasmic half in our simulations. However, the closure is sufficient to prevent the disaccharide from transporting to the cytoplasm, which is in accordance with the well-established alternating access model. Ser53, Gln60, and Phe354 are determined to be important in sugar transport during the periplasmic-open stage of the sugar transport cycle and the sugar is found to undergo an orientational change in order to escape the protein lumen.     

The pH-induced release of iron from transferrin investigated with a continuum electrostatic model.

A reduction in pH induces the release of iron from transferrin in a process that involves a conformational change in the protein from a closed to an open form. Experimental evidence suggests that there must be changes in the protonation states of certain, as yet not clearly identified, residues in the protein accompanying this conformational change. Such changes in protonation states of residues and the consequent changes in electrostatic interactions are assumed to play a large part in the mechanism of release of iron from transferrin. Using the x-ray crystal structures of human ferri- and apo-lactoferrin, we calculated the pKa values of the titratable residues in both the closed (iron-loaded) and open (iron-free) conformations with a continuum electrostatic model. With the knowledge of a residue's pKa value, its most probable protonation state at any specified pH may be determined. The preliminary results presented here are in good agreement with the experimental observation that the binding of ferric iron and the synergistic anion bicarbonate/carbonate results in the release of approximately three H+ ions. It is suggested that the release of these three H+ ions may be accounted for, in most part, by the deprotonation of the bicarbonate and residues Tyr-92, Lys-243, Lys-282, and Lys-285 together with the protonation of residues Asp-217 and Lys-277.     

Structure-function relationships of the human bitter taste receptor hTAS2R1: insights from molecular modeling studies

Bitter taste receptors (T2Rs) belong to G-protein-coupled receptors (GPCRs). Despite extensive studies, the precise mechanisms of GPCR activation are still poorly understood. In this study, the models of the human bitter taste receptor hTAS2R1 alone and in complex with various ligands were constructed on the basis of template-based modeling and molecular docking. Then these models were subjected to all-atom molecular dynamics (MD) simulations in explicit lipid bilayers. The binding pocket of hTAS2R1 is mainly formed by transmembrane helix (TM) III, TM V, TM VI, and TM VII. Most of the residues contributing to ligand binding are positionally conserved comparing with other hTAS2Rs. By comparing the final conformations obtained by extensive MD simulations, we identified the changes in the transmembrane helices and the intra- and extracellular loops, which were expected to initiate the activation of the receptor. The intracellular loop II (ICL2) and TM III were found to play prominent roles in the process of activation. We proposed that a set of interactions between the aromatic Phe115 in the middle of ICL2 and three residues (Tyr103, Lys106, and Val107) at the cytoplasmic end of TM III may serve as a conformational switch of hTAS2R1 activation. All of the residues involved in the switch are highly conserved among T2Rs. This indicates that the control switch we proposed may be universal in T2Rs. Besides, our results also suggest that the formation of a short helical segment in ICL2 may be necessary for the activation of hTAS2R1.     

Molecular mechanics simulations of a conformational rearrangement of D-xylose in the active site of D-xylose isomerase.

A proposed reaction mechanism for the enzyme D-xylose isomerase involves the ring opening of the cyclic substrate with a subsequent conformational rearrangement to an extended open-chain form. Restrained energy minimization was used to simulate the rearrangement. In the ring-opening step, the substrate energy function was gradually altered from a cyclic to an open-chain form, with energy minimization after each change. The protein/sugar contact energy did not increase significantly during the process, showing that there was no steric hindrance to ring opening. The conformational rearrangement involves an alteration in the coordination of the substrate to metal ion [1], which was induced by gradually changing restraints on metal/ligand distances. By allowing varying amounts of flexibility in the protein and examining a simplified model system, the interactions of the sugar with metal ion [1] and its immediate ligands were found to be the most important contributors to the energy barrier for the change. Only small changes in the positions of protein atoms were required. The energy barrier to the rearrangement was estimated to be less than the Arrhenius activation energy for the enzymatic reaction. This is in accordance with experimental indications that the isomerization step is rate determining.     

Ligand variation in the transferrin family: the crystal structure of the H249Q mutant of the human transferrin N-lobe as a model for iron binding in insect transferrins

Proteins of the transferrin (Tf) family play a central role in iron homeostasis in vertebrates. In vertebrate Tfs, the four iron-binding ligands, 1 Asp, 2 Tyr, and 1 His, are invariant in both lobes of these bilobal proteins. In contrast, there are striking variations in the Tfs that have been characterized from insect species; in three of them, sequence changes in the C-lobe binding site render it nonfunctional, and in all of them the His ligand in the N-lobe site is changed to Gln. Surprisingly, mutagenesis of the histidine ligand, His249, to glutamine in the N-lobe half-molecule of human Tf (hTf/2N) shows that iron binding is destabilized and suggests that Gln249 does not bind to iron. We have determined the crystal structure of the H249Q mutant of hTf/2N and refined it at 1.85 A resolution (R = 0.221, R(free) = 0.246). The structure reveals that Gln249 does coordinate to iron, albeit with a lengthened Fe-Oepsilon1 bond of 2.34 A. In every other respect, the protein structure is unchanged from wild-type. Examination of insect Tf sequences shows that the K206.K296 dilysine pair, which aids iron release from the N-lobes of vertebrate Tfs, is not present in the insect proteins. We conclude that substitution of Gln for His does destabilize iron binding, but in the insect Tfs this is compensated by the loss of the dilysine interaction. The combination of a His ligand with the dilysine pair in vertebrate Tfs may have been a later evolutionary development that gives more sophisticated pH-mediated control of iron release from the N-lobe of transferrins.