Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

Branching in the sequential folding pathway of cytochrome c

Previous results indicate that the folding pathways of cytochrome c and other proteins progressively build the target native protein in a predetermined stepwise manner by the sequential formation and association of native-like foldon units. The present work used native state hydrogen exchange methods to investigate a structural anomaly in cytochrome c results that suggested the concerted folding of two segments that have little structural relationship in the native protein. The results show that the two segments, an 18-residue omega loop and a 10-residue helix, are able to unfold and refold independently, which allows a branch point in the folding pathway. The pathway that emerges assembles native-like foldon units in a linear sequential manner when prior native-like structure can template a single subsequent foldon, and optional pathway branching is seen when prior structure is able to support the folding of two different foldons.     

Downhill versus Barrier-Limited Folding of BBL 1: Energetic and Structural Perturbation Effects upon Protonation of a Histidine of Unusually Low pKa

A dispersion of melting temperatures at pH 5.3 for individual residues of the BBL protein domain has been adduced as evidence for barrier-free downhill folding. Other members of the peripheral subunit domain family fold cooperatively at pH 7. To search for possible causes of anomalies in BBL's denaturation behavior, we measured the pH titration of individual residues by heteronuclear NMR. At 298 K, the pK_a of His142 was close to that of free histidine at 6.47 ± 0.04, while that of the more buried His166 was highly perturbed at 5.39 ± 0.02. Protonation of His166 is thus energetically unfavorable and destabilizes the protein by ∼ 1.5 kcal/mol. Changes in C^α secondary shifts at pH 5.3 showed a decrease in helicity of the C-terminus of helix 2, where His166 is located, which was accompanied by a measured decrease of 1.1 ± 0.2 kcal/mol in stability from pH 7 to 5.3. Protonation of His166 perturbs, therefore, the structure of BBL. Only ∼ 1% of the structurally perturbed state will be present at the biologically relevant pH 7.6. Experiments at pH 5.3 report on a near-equal mixture of the two different native states. Further, at this pH, small changes of pH and pK_a induced by changes in temperature will have near-maximal effects on pH-dependent conformational equilibria and on propagation of experimental error. Accordingly, conventional barrier-limited folding predicts some dispersion of measured thermal unfolding curves of individual residues at pH 5.3.     

How cytochrome c folds, and why: submolecular foldon units and their stepwise sequential stabilization

Native state hydrogen exchange experiments have shown that the cytochrome c (Cyt c) protein consists of five cooperative folding-unfolding units, called foldons. These are named, in the order of increasing unfolding free energy, the nested-Yellow, Red, Yellow, Green, and Blue foldons. Previous results suggest that these units unfold in a stepwise sequential way so that each higher energy partially unfolded form includes all of the previously unfolded lower free energy units. If this is so, then selectively destabilizing any given foldon should equally destabilize each subsequent unfolding step above it in the unfolding ladder but leave the lower ones before it unaffected. To perform this test, we introduced the mutation Glu62Gly, which deletes a salt link in the Yellow unit and destabilizes the protein by 0.8 kcal/mol. Native state hydrogen exchange and other experiments show that the stability of the Yellow unit and the states above it in the free energy ladder are destabilized by about the same amount while the lower lying states are unaffected. These results help to confirm the sequential stepwise nature of the Cyt c unfolding pathway and therefore a similar refolding pathway. The steps in the pathway are dictated by the concerted folding-unfolding property of the individual unit foldons; the order of steps is determined by the sequential stabilization of progressively added foldons in the native context. Much related information for Cyt c strongly conforms with this mechanism. Its generality is supported by available information for other proteins.     

Increased folding stability of TEM-1 beta-lactamase by in vitro selection

In vitro selections of stabilized proteins lead to more robust enzymes and, at the same time, yield novel insights into the principles of protein stability. We employed Proside, a method of in vitro selection, to find stabilized variants of TEM-1 β-lactamase from Escherichia coli. Proside links the increased protease resistance of stabilized proteins to the infectivity of a filamentous phage. Several libraries of TEM-1 β-lactamase variants were generated by error-prone PCR, and variants with increased protease resistance were obtained by raising temperature or guanidinium chloride concentration during proteolytic selections. Despite the small size of phage libraries, several strongly stabilizing mutations could be obtained, and a manual combination of the best shifted the profiles for thermal unfolding and temperature-dependent inactivation of β-lactamase by almost 20 °C to a higher temperature. The wild-type protein unfolds in two stages: from the native state via an intermediate of the molten-globule type to the unfolded form. In the course of the selections, the native protein was stabilized by 27 kJ mol^− 1 relative to the intermediate and the cooperativity of unfolding was strongly increased. Three of our stabilizing replacements (M182T, A224V, and R275L) had been identified independently in naturally occurring β-lactamase variants with extended substrate spectrum. In these variants, they acted as global suppressors of destabilizations caused by the mutations in the active site. The comparison between the crystal structure of our best variant and the crystal structure of the wild-type protein indicates that most of the selected mutations optimize helices and their packing. The stabilization by the E147G substitution is remarkable. It removes steric strain that originates from an overly tight packing of two helices in the wild-type protein. Such unfavorable van der Waals repulsions are not easily identified in crystal structures or by computational approaches, but they strongly reduce the conformational stability of a protein.     

Thermodynamic and structural analysis of homodimeric proteins: Model of β-lactoglobulin

The energetics of protein homo-oligomerization was analyzed in detail with the application of a general thermodynamic model. We have studied the thermodynamic aspects of protein-protein interaction employing β-lactoglobulin A from bovine milk at pH = 6.7 where the protein is mainly in its dimeric form. We performed differential calorimetric scans at different total protein concentration and the resulting thermograms were analyzed with the thermodynamic model for oligomeric proteins previously developed. The thermodynamic model employed, allowed the prediction of the sign of the enthalpy of dimerization, the analysis of complex calorimetric profiles without transitions baselines subtraction and the obtainment of the thermodynamic parameters from the unfolding and the association processes and the compared with association parameters obtained with Isothermal Titration Calorimetry performed at different temperatures. The dissociation and unfolding reactions were also monitored by Fourier-transform infrared spectroscopy and the results indicated that the dimer of β-lactoglobulin (N₂) reversibly dissociates into monomeric units (N) which are structurally distinguishable by changes in their infrared absorbance spectra upon heating. Hence, it is proposed that β-lactoglobulin follows the conformational path induced by temperature:N₂ ? 2N ? 2D. The general model was validated with these results indicating that it can be employed in the study of the thermodynamics of other homo-oligomeric protein systems.     

NMR analysis of partially folded states and persistent structure in the alpha subunit of tryptophan synthase: implications for the equilibrium folding mechanism of a 29-kDa TIM barrel protein

Structural insights into the equilibrium folding mechanism of the alpha subunit of tryptophan synthase (αTS) from Escherichia coli, a (βα)₈ TIM barrel protein, were obtained with a pair of complementary nuclear magnetic resonance (NMR) spectroscopic techniques. The secondary structures of rare high-energy partially folded states were probed by native-state hydrogen-exchange NMR analysis of main-chain amide hydrogens. 2D heteronuclear single quantum coherence NMR analysis of several ¹⁵N-labeled nonpolar amino acids was used to probe the side chains involved in stabilizing a highly denatured intermediate that is devoid of secondary structure. The dynamic broadening of a subset of isoleucine and leucine side chains and the absence of protection against exchange showed that the highest energy folded state on the free-energy landscape is stabilized by a hydrophobic cluster lacking stable secondary structure. The core of this cluster, centered near the N-terminus of αTS, serves as a nucleus for the stabilization of what appears to be nonnative secondary structure in a marginally stable intermediate. The progressive decrease in protection against exchange from this nucleus toward both termini and from the N-termini to the C-termini of several β-strands is best described by an ensemble of weakly coupled conformers. Comparison with previous data strongly suggests that this ensemble corresponds to a marginally stable off-pathway intermediate that arises in the first few milliseconds of folding and persists under equilibrium conditions. A second, more stable intermediate, which has an intact β-barrel and a frayed α-helical shell, coexists with this marginally stable species. The conversion of the more stable intermediate to the native state of αTS entails the formation of a stable helical shell and completes the acquisition of the tertiary structure.     

Protein hydrogen exchange mechanism: local fluctuations

Experiments were done to study the dynamic structural motions that determine protein hydrogen exchange (HX) behavior. The replacement of a solvent-exposed lysine residue with glycine (Lys8Gly) in a helix of recombinant cytochrome c does not perturb the native structure, but it entropically potentiates main-chain flexibility and thus can promote local distortional motions and large-scale unfolding. The mutation accelerates amide hydrogen exchange of the mutated residue by about 50-fold, neighboring residues in the same helix by less, and residues elsewhere in the protein not at all, except for Leu98, which registers the change in global stability. The pattern of HX changes shows that the coupled structural distortions that dominate exchange can be several residues in extent, but they expose to exchange only one amide NH at a time. This "local fluctuation" mode of hydrogen exchange may be generally recognized by disparate near-neighbor rates and a low dependence on destabilants (denaturant, temperature, pressure). In contrast, concerted unfolding reactions expose multiple neighboring amide NHs with very similar computed protection factors, and they show marked destabilant sensitivity. In both modes, ionic hydrogen exchange catalysts attack from the bulk solvent without diffusing through the protein matrix.     

The cDNA sequence of three hemocyanin subunits from the garden snail Helix lucorum

Hemocyanins are blue copper containing respiratory proteins residing in the hemolymph of many molluscs and arthropods. They can have different molecular masses and quaternary structures. Moreover, several molluscan hemocyanins are isolated with one, two or three isoforms occurring as decameric, didecameric, multidecameric or tubule aggregates. We could recently isolate three different hemocyanin isopolypeptides from the hemolymph of the garden snail Helix lucorum (HlH). These three structural subunits were named α_D-HlH, α_N-HlH and β-HlH. We have cloned and sequenced their cDNA which is the first result ever reported for three isoforms of a molluscan hemocyanin. Whereas the complete gene sequence of α_D-HlH and β-HlH was obtained, including the 5′ and 3′ UTR, 180 bp of the 5′ end and around 900 bp at the 3′ end are missing for the third subunit. The subunits α_D-HlH and β-HlH comprise a signal sequence of 19 amino acids plus a polypeptide of 3409 and 3414 amino acids, respectively. We could determine 3031 residues of the α_N-HLH subunit. Sequence comparison with other molluscan hemocyanins shows that α_D-HlH is more related to Aplysia californicum hemocyanin than to each of its own isopolypeptides. The structural subunits comprise 8 different functional units (FUs: a, b, c, d, e, f, g, h) and each functional unit possesses a highly conserved copper-A and copper-B site for reversible oxygen binding. Potential N-glycosylation sites are present in all three structural subunits. We confirmed that all three different isoforms are effectively produced and secreted in the hemolymph of H. lucorum by analyzing a tryptic digest of the purified native hemocyanin by MALDI-TOF and LC-FTICR mass spectrometry.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Protein hydrogen exchange: testing current models

To investigate the determinants of protein hydrogen exchange (HX), HX rates of most of the backbone amide hydrogens of Staphylococcal nuclease were measured by NMR methods. A modified analysis was used to improve accuracy for the faster hydrogens. HX rates of both near surface and well buried hydrogens are spread over more than 7 orders of magnitude. These results were compared with previous hypotheses for HX rate determination. Contrary to a common assumption, proximity to the surface of the native protein does not usually produce fast exchange. The slow HX rates for unprotected surface hydrogens are not well explained by local electrostatic field. The ability of buried hydrogens to exchange is not explained by a solvent penetration mechanism. The exchange rates of structurally protected hydrogens are not well predicted by algorithms that depend only on local interactions or only on transient unfolding reactions. These observations identify some of the present difficulties of HX rate prediction and suggest the need for returning to a detailed hydrogen by hydrogen analysis to examine the bases of structure-rate relationships, as described in the companion paper (Skinner et al., Protein Sci 2012;21:996-1005).     

Crystal structure of glycoside hydrolase family 78 alpha-L-Rhamnosidase from Bacillus sp. GL1

α-L-Rhamnosidase (EC 3.2.1.40) catalyzes the hydrolytic release of rhamnose from polysaccharides and glycosides. Bacillus sp. GL1 α-L-rhamnosidase (RhaB), a member of glycoside hydrolase (GH) family 78, is responsible for degrading the bacterial biofilm gellan, and also functions as a debittering agent for citrus fruit in the food and beverage industries through the release of rhamnose from plant glycoside, naringin. The X-ray crystal structure of RhaB was determined by single-wavelength anomalous diffraction using a selenomethionine derivative and refined at 1.9 Å resolution with a final R-factor of 18.2%. As is seen in the homodimeric form of the active enzyme, the structure of RhaB in crystal packing is a homodimer containing 1908 amino acids (residues 3-956), 43 glycerol molecules, four calcium ions, and 1755 water molecules. The overall structure consists of five domains, four of which are β-sandwich structures designated as domains N, D1, D2, and C, and an (α/α)₆-barrel structure designated as domain A. Structural comparison by DALI showed that RhaB shares its highest level of structural similarity with chitobiose phosphorylase (Z score of 25.3). The structure of RhaB in complex with the reaction product rhamnose (inhibitor constant, K_i = 1.8 mM) was also determined and refined at 2.1 Å with a final R-factor of 19.5%. Rhamnose is bound to the deep cleft of the (α/α)₆-barrel domain, as is seen in the clan-L GHs. Several negatively charged residues, such as Asp567, Glu572, Asp579, and Glu841, conserved in GH family 78 enzymes, interact with rhamnose, and RhaB mutants of these residues have drastically reduced enzyme activity, indicating that the residues are crucial for enzyme catalysis and/or substrate binding. To our knowledge, this is the first report on the determination of the crystal structure of α-L-rhamnosidase and identification of its clan-L (α/α)₆-barrel as a catalytic domain.     

Crystal Structure of an Essential Enzyme in Seed Starch Degradation: Barley Limit Dextrinase in Complex with Cyclodextrins

Barley limit dextrinase [Hordeum vulgare limit dextrinase (HvLD)] catalyzes the hydrolysis of α-1,6 glucosidic linkages in limit dextrins. This activity plays a role in starch degradation during germination and presumably in starch biosynthesis during grain filling. The crystal structures of HvLD in complex with the competitive inhibitors α-cyclodextrin (CD) and β-CD are solved and refined to 2.5 Å and 2.1 Å, respectively, and are the first structures of a limit dextrinase. HvLD belongs to glycoside hydrolase 13 family and is composed of four domains: an immunoglobulin-like N-terminal eight-stranded β-sandwich domain, a six-stranded β-sandwich domain belonging to the carbohydrate binding module 48 family, a catalytic (β/α)₈-like barrel domain that lacks α-helix 5, and a C-terminal eight-stranded β-sandwich domain of unknown function. The CDs are bound at the active site occupying carbohydrate binding subsites + 1 and + 2. A glycerol and three water molecules mimic a glucose residue at subsite − 1, thereby identifying residues involved in catalysis. The bulky Met440, a unique residue at its position among α-1,6 acting enzymes, obstructs subsite − 4. The steric hindrance observed is proposed to affect substrate specificity and to cause a low activity of HvLD towards amylopectin. An extended loop (Asp513-Asn520) between β5 and β6 of the catalytic domain also seems to influence substrate specificity and to give HvLD a higher affinity for α-CD than pullulanases. The crystal structures additionally provide new insight into cation sites and the concerted action of the battery of hydrolytic enzymes in starch degradation.     

The role of a disulfide bridge in the stability and folding kinetics of Arabidopsis thaliana cytochrome c6A

Cytochrome c_6A is a eukaryotic member of the Class I cytochrome c family possessing a high structural homology with photosynthetic cytochrome c₆ from cyanobacteria, but structurally and functionally distinct through the presence of a disulfide bond and a heme mid-point redox potential of + 71 mV (vs normal hydrogen electrode). The disulfide bond is part of a loop insertion peptide that forms a cap-like structure on top of the core α-helical fold. We have investigated the contribution of the disulfide bond to thermodynamic stability and (un)folding kinetics in cytochrome c_6A from Arabidopsis thaliana by making comparison with a photosynthetic cytochrome c₆ from Phormidium laminosum and through a mutant in which the Cys residues have been replaced with Ser residues (C67/73S). We find that the disulfide bond makes a significant contribution to overall stability in both the ferric and ferrous heme states. Both cytochromes c_6A and c₆ fold rapidly at neutral pH through an on-pathway intermediate. The unfolding rate for the C67/73S variant is significantly increased indicating that the formation of this region occurs late in the folding pathway. We conclude that the disulfide bridge in cytochrome c_6A acts as a conformational restraint in both the folding intermediate and native state of the protein and that it likely serves a structural rather than a previously proposed catalytic role.     

Mutational analysis of the stability of the H2A and H2B histone monomers

The eukaryotic histone heterodimer H2A-H2B folds through an obligatory dimeric intermediate that forms in a nearly diffusion-limited association reaction in the stopped-flow dead time. It is unclear whether there is partial folding of the isolated monomers before association. To address the possible contributions of structure in the monomers to the rapid association, we characterized H2A and H2B monomers in the absence of their heterodimeric partner. By far-UV circular dichroism, the H2A and H2B monomers are 15% and 31% helical, respectively—significantly less than observed in X-ray crystal structures. Acrylamide quenching of the intrinsic Tyr fluorescence was indicative of tertiary structure. The H2A and H2B monomers exhibit free energies of unfolding of 2.5 and 2.9 kcal mol^− 1, respectively; at 10 μM, the sum of the stability of the monomers is ∼ 60% of the stability of the native dimer. The helical content, stability, and m values indicate that H2B has a more stable, compact structure than H2A. The monomer m values are larger than expected for the extended histone fold motif, suggesting that the monomers adopt an overly collapsed structure. Stopped-flow refolding—initiated from urea-denatured monomers or the partially folded monomers populated at low denaturant concentrations—yielded essentially identical rates, indicating that monomer folding is productive in the rapid association and folding of the heterodimer. A series of Ala and Gly mutations were introduced into H2A and H2B to probe the importance of helix propensity on the structure and stability of the monomers. The mutational studies show that the central α-helix of the histone fold, which makes extensive intermonomer contacts, is structured in H2B but only partially folded in H2A.     

Crystal structure of α-galactosidase from Lactobacillus acidophilus NCFM: insight into tetramer formation and substrate binding

Lactobacillus acidophilus NCFM is a probiotic bacterium known for its beneficial effects on human health. The importance of α-galactosidases (α-Gals) for growth of probiotic organisms on oligosaccharides of the raffinose family present in many foods is increasingly recognized. Here, the crystal structure of α-Gal from L. acidophilus NCFM (LaMel36A) of glycoside hydrolase (GH) family 36 (GH36) is determined by single-wavelength anomalous dispersion. In addition, a 1.58-Å-resolution crystallographic complex with α-d-galactose at substrate binding subsite − 1 was determined. LaMel36A has a large N-terminal twisted β-sandwich domain, connected by a long α-helix to the catalytic (β/α)₈-barrel domain, and a C-terminal β-sheet domain. Four identical monomers form a tightly packed tetramer where three monomers contribute to the structural integrity of the active site in each monomer. Structural comparison of LaMel36A with the monomeric Thermotoga maritima α-Gal (TmGal36A) reveals that O2 of α-d-galactose in LaMel36A interacts with a backbone nitrogen in a glycine-rich loop of the catalytic domain, whereas the corresponding atom in TmGal36A is from a tryptophan side chain belonging to the N-terminal domain. Thus, two distinctly different structural motifs participate in substrate recognition. The tetrameric LaMel36A furthermore has a much deeper active site than the monomeric TmGal36A, which possibly modulates substrate specificity. Sequence analysis of GH36, inspired by the observed structural differences, results in four distinct subgroups having clearly different active-site sequence motifs. This novel subdivision incorporates functional and architectural features and may aid further biochemical and structural analyses within GH36.     

Intimate view of a kinetic protein folding intermediate: residue-resolved structure, interactions, stability, folding and unfolding rates, homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.     

The regulatory factor SipA is a highly stable β-II class protein with a SH3 fold

The small regulator SipA, interacts with the ATP-binding domain of non-bleaching sensor histidine kinase (NblS), the most conserved histidine kinase in cyanobacteria. NblS regulates photosynthesis and acclimation to a variety of environmental conditions. We show here that SipA is a highly stable protein in a wide pH range, with a thermal denaturation midpoint of 345 K. Circular dichroism and 1D ¹H NMR spectroscopies, as well as modelling, suggest that SipA is a β-II class protein, with short strands followed by turns and long random-coil polypeptide patches, matching the SH3 fold. The experimentally determined m-value and the heat capacity change upon thermal unfolding (ΔC_p) closely agreed with the corresponding theoretical values predicted from the structural model, further supporting its accuracy.     

A new β-chain haemoglobin variant with increased oxygen affinity: Hb Roma [β115(g17)Ala → Val]

Background

Haemoglobin Roma [β115(G17)Ala → Val] is a new adult haemoglobin variant found in a patient presenting a mild hypochromia and microcytosis. We studied this previously uncharacterised variant in order to evaluate the effect on the structural and funcional properties of the Ala → Val substitution at the α1β1 interface.

Methods and results

The variant chain was identified by direct DNA sequencing of the β-globin gene, which revealed a GCC → GTC mutation in codon 115. This mutation was confirmed by mass spectrometric analysis of the tetramers and peptides. The oxygen-binding properties of the haemoglobin A/haemoglobin Roma mixture, in which the variant makes up 25% of the haemoglobins, showed a significant increase in oxygen affinity with respect to normal haemoglobin A, both in the absence and presence of 2,3-bisphosphoglycerate. The role of the βG17 position, situated at the α₁β₁ interface, has been examined using computational models of haemoglobin Roma and other known βG17 variants, in comparison with normal haemoglobin A.

Conclusions

This study suggests that the β115(G17)Ala → Val substitution at the α1β1 interface is responsible for increased oxygen affinity and mild destabilisation of the haemoglobin Roma.

General significance

An amino acid substitution at the G17 position of the α1β1 interface may result in stabilisation of the high affinity R-state of the haemoglobin molecule.