Gene Duplication and the Evolution of Hemoglobin Isoform Differentiation in Birds

The majority of bird species co-express two functionally distinct hemoglobin (Hb) isoforms in definitive erythrocytes as follows: HbA (the major adult Hb isoform, with α-chain subunits encoded by the α^A-globin gene) and HbD (the minor adult Hb isoform, with α-chain subunits encoded by the α^D-globin gene). The α^D-globin gene originated via tandem duplication of an embryonic α-like globin gene in the stem lineage of tetrapod vertebrates, which suggests the possibility that functional differentiation between the HbA and HbD isoforms may be attributable to a retained ancestral character state in HbD that harkens back to a primordial, embryonic function. To investigate this possibility, we conducted a combined analysis of protein biochemistry and sequence evolution to characterize the structural and functional basis of Hb isoform differentiation in birds. Functional experiments involving purified HbA and HbD isoforms from 11 different bird species revealed that HbD is characterized by a consistently higher O₂ affinity in the presence of allosteric effectors such as organic phosphates and Cl⁻ ions. In the case of both HbA and HbD, analyses of oxygenation properties under the two-state Monod-Wyman-Changeux allosteric model revealed that the pH dependence of Hb-O₂ affinity stems primarily from changes in the O₂ association constant of deoxy (T-state)-Hb. Ancestral sequence reconstructions revealed that the amino acid substitutions that distinguish the adult-expressed Hb isoforms are not attributable to the retention of an ancestral (pre-duplication) character state in the α^D-globin gene that is shared with the embryonic α-like globin gene.     

Post-translational Transformation of Methionine to Aspartate Is Catalyzed by Heme Iron and Driven by Peroxide: A NOVEL SUBUNIT-SPECIFIC MECHANISM IN HEMOGLOBIN*

A pathogenic V67M mutation occurs at the E11 helical position within the heme pockets of variant human fetal and adult hemoglobins (Hb). Subsequent post-translational modification of Met to Asp was reported in γ subunits of human fetal Hb Toms River (γ67(E11)Val → Met) and β subunits of adult Hb (HbA) Bristol-Alesha (β67(E11)Val → Met) that were associated with hemolytic anemia. Using kinetic, proteomic, and crystal structural analysis, we were able to show that the Met → Asp transformation involves heme cycling through its oxoferryl state in the recombinant versions of both proteins. The conversion to Met and Asp enhanced the spontaneous autoxidation of the mutants relative to wild-type HbA and human fetal Hb, and the levels of Asp were elevated with increasing levels of hydrogen peroxide (H₂O₂). Using H₂¹⁸O₂, we verified incorporation of ¹⁸O into the Asp carboxyl side chain confirming the role of H₂O₂ in the oxidation of the Met side chain. Under similar experimental conditions, there was no conversion to Asp at the αMet(E11) position in the corresponding HbA Evans (α62(E11)Val → Met). The crystal structures of the three recombinant Met(E11) mutants revealed similar thioether side chain orientations. However, as in the solution experiments, autoxidation of the Hb mutant crystals leads to electron density maps indicative of Asp(E11) formation in β subunits but not in α subunits. This novel post-translational modification highlights the nonequivalence of human Hb α, β, and γ subunits with respect to redox reactivity and may have direct implications to α/β hemoglobinopathies and design of oxidatively stable Hb-based oxygen therapeutics.     

Autoxidation and Oxygen Binding Properties of Recombinant Hemoglobins with Substitutions at the αVal-62 or βVal-67 Position of the Distal Heme Pocket

The E11 valine in the distal heme pocket of either the α- or β-subunit of human adult hemoglobin (Hb A) was replaced by leucine, isoleucine, or phenylalanine. Recombinant proteins were expressed in Escherichia coli and purified for structural and functional studies. ¹H NMR spectra were obtained for the CO and deoxy forms of Hb A and the mutants. The mutations did not disturb the α₁β₂ interface in either form, whereas the H-bond between αHis-103 and βGln-131 in the α₁β₁ interfaces of the deoxy α-subunit mutants was weakened. Localized structural changes in the mutated heme pocket were detected for the CO form of recombinant Hb (rHb) (αV62F), rHb (βV67I), and rHb (βV67F) compared with Hb A. In the deoxy form the proximal histidyl residue in the β-subunit of rHb (βV67F) has been altered. Furthermore, the interactions between the porphyrin ring and heme pocket residues have been perturbed in rHb (αV62I), rHb (αV62F), and rHb (βV67F). Functionally, the oxygen binding affinity (P₅₀), cooperativity (n₅₀), and the alkaline Bohr Effect of the three α-subunit mutants and rHb (βV67L) are similar to those of Hb A. rHb (βV67I) and rHb (βV67F) exhibit low and high oxygen affinity, respectively. rHb (βV67F) has P₅₀ values lower that those reported for rHb (αL29F), a B10 mutant studied previously in our laboratory (Wiltrout, M. E., Giovannelli, J. L., Simplaceanu, V., Lukin, J. A., Ho, N. T., and Ho, C. (2005) Biochemistry 44, 7207–7217). These E11 mutations do not slow down the autoxidation and azide-induced oxidation rates of the recombinant proteins. Results from this study provide new insights into the roles of E11 mutants in the structure-function relationship in hemoglobin.     

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73–96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27–30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.     

α-Hemoglobin-stabilizing Protein (AHSP) Perturbs the Proximal Heme Pocket of Oxy-α-hemoglobin and Weakens the Iron-Oxygen Bond

α-Hemoglobin (αHb)-stabilizing protein (AHSP) is a molecular chaperone that assists hemoglobin assembly. AHSP induces changes in αHb heme coordination, but how these changes are facilitated by interactions at the αHb·AHSP interface is not well understood. To address this question we have used NMR, x-ray absorption spectroscopy, and ligand binding measurements to probe αHb conformational changes induced by AHSP binding. NMR chemical shift analyses of free CO-αHb and CO-αHb·AHSP indicated that the seven helical elements of the native αHb structure are retained and that the heme Fe(II) remains coordinated to the proximal His-87 side chain. However, chemical shift differences revealed alterations of the F, G, and H helices and the heme pocket of CO-αHb bound to AHSP. Comparisons of iron-ligand geometry using extended x-ray absorption fine structure spectroscopy showed that AHSP binding induces a small 0.03 Å lengthening of the Fe-O₂ bond, explaining previous reports that AHSP decreases αHb O₂ affinity roughly 4-fold and promotes autooxidation due primarily to a 3–4-fold increase in the rate of O₂ dissociation. Pro-30 mutations diminished NMR chemical shift changes in the proximal heme pocket, restored normal O₂ dissociation rate and equilibrium constants, and reduced O₂-αHb autooxidation rates. Thus, the contacts mediated by Pro-30 in wild-type AHSP promote αHb autooxidation by introducing strain into the proximal heme pocket. As a chaperone, AHSP facilitates rapid assembly of αHb into Hb when βHb is abundant but diverts αHb to a redox resistant holding state when βHb is limiting.     

Sickle Cell Hemoglobin with Mutation at αHis-50 Has Improved Solubility

The unliganded tetrameric Hb S has axial and lateral contacts with neighbors and can polymerize in solution. Novel recombinants of Hb S with single amino acid substitutions at the putative axial (recombinant Hb (rHb) (βE6V/αH20R) and rHb (βE6V/αH20Q)) or lateral (rHb (βE6V/αH50Q)) or double amino acid substitutions at both the putative axial and lateral (rHb (βE6V/αH20R/αH50Q) and rHb (βE6V/αH20Q/αH50Q)) contact sites were expressed in Escherichia coli and purified for structural and functional studies. The ¹H NMR spectra of the CO and deoxy forms of these mutants indicate that substitutions at either αHis-20 or αHis-50 do not change the subunit interfaces or the heme pockets of the proteins. The double mutants show only slight structural alteration in the β-heme pockets. All mutants have similar cooperativity (n₅₀), alkaline Bohr effect, and autoxidation rate as Hb S. The oxygen binding affinity (P₅₀) of the single mutants is comparable with that of Hb S. The double mutants bind oxygen with slightly higher affinity than Hb S under the acidic conditions. In high salt, rHb (βE6V/αH20R) is the only mutant that has a shorter delay time of polymerization and forms polymers more readily than Hb S with a dextran-Csat value of 1.86 ± 0.20 g/dl. Hb S, rHb (βE6V/αH20Q), rHb (βE6V/αH50Q), rHb (βE6V/αH20R/αH50Q), and rHb (βE6V/αH20Q/αH50Q) have dextran-Csat values of 2.95 ± 0.10, 3.04 ± 0.17, 11.78 ± 0.59, 7.11 ± 0.66, and 10.89 ± 0.83 g/dl, respectively. rHb (βE6V/αH20Q/αH50Q) is even more stable than Hb S under elevated temperature (60 °C).     

Subunit δ Is the Key Player for Assembly of the H+-translocating Unit of Escherichia coli FOF1 ATP Synthase

The ATP synthase (F_OF₁) of Escherichia coli couples the translocation of protons across the cytoplasmic membrane to the synthesis or hydrolysis of ATP. This nanomotor is composed of the rotor c₁₀γϵ and the stator ab₂α₃β₃δ. To study the assembly of this multimeric enzyme complex consisting of membrane-integral as well as peripheral hydrophilic subunits, we combined nearest neighbor analyses by intermolecular disulfide bond formation or purification of partially assembled F_OF₁ complexes by affinity chromatography with the use of mutants synthesizing different sets of F_OF₁ subunits. Together with a time-delayed in vivo assembly system, the results demonstrate that F_OF₁ is assembled in a modular way via subcomplexes, thereby preventing the formation of a functional H⁺-translocating unit as intermediate product. Surprisingly, during the biogenesis of F_OF₁, F₁ subunit δ is the key player in generating stable F_O. Subunit δ serves as clamp between ab₂ and c₁₀α₃β₃γϵ and guarantees that the open H⁺ channel is concomitantly assembled within coupled F_OF₁ to maintain the low membrane proton permeability essential for viability, a general prerequisite for the assembly of multimeric H⁺-translocating enzymes.     

An Origin of Cooperative Oxygen Binding of Human Adult Hemoglobin: Different Roles of the α and β Subunits in the α2β2 Tetramer

Human hemoglobin (Hb), which is an α₂β₂ tetramer and binds four O₂ molecules, changes its O₂-affinity from low to high as an increase of bound O₂, that is characterized by ‘cooperativity’. This property is indispensable for its function of O₂ transfer from a lung to tissues and is accounted for in terms of T/R quaternary structure change, assuming the presence of a strain on the Fe-histidine (His) bond in the T state caused by the formation of hydrogen bonds at the subunit interfaces. However, the difference between the α and β subunits has been neglected. To investigate the different roles of the Fe-His(F8) bonds in the α and β subunits, we investigated cavity mutant Hbs in which the Fe-His(F8) in either α or β subunits was replaced by Fe-imidazole and F8-glycine. Thus, in cavity mutant Hbs, the movement of Fe upon O₂-binding is detached from the movement of the F-helix, which is supposed to play a role of communication. Recombinant Hb (rHb)(αH87G), in which only the Fe-His in the α subunits is replaced by Fe-imidazole, showed a biphasic O₂-binding with no cooperativity, indicating the coexistence of two independent hemes with different O₂-affinities. In contrast, rHb(βH92G), in which only the Fe-His in the β subunits is replaced by Fe-imidazole, gave a simple high-affinity O₂-binding curve with no cooperativity. Resonance Raman, ¹H NMR, and near-UV circular dichroism measurements revealed that the quaternary structure change did not occur upon O₂-binding to rHb(αH87G), but it did partially occur with O₂-binding to rHb(βH92G). The quaternary structure of rHb(αH87G) appears to be frozen in T while its tertiary structure is changeable. Thus, the absence of the Fe-His bond in the α subunit inhibits the T to R quaternary structure change upon O₂-binding, but its absence in the β subunit simply enhances the O₂-affinity of α subunit.     

Subunit Composition Determines the Single Channel Kinetics of the Epithelial Sodium Channel

We have further characterized at the single channel level the properties of epithelial sodium channels formed by coexpression of α with either wild-type β or γ subunits and α with carboxy-terminal truncated β (β_T) or γ (γ_T) subunits in Xenopus laevis oocytes. αβ and αβ_T channels (9.6 and 8.7 pS, respectively, with 150 mM Li⁺) were found to be constitutively open. Only upon inclusion of 1 μM amiloride in the pipette solution could channel activity be resolved; both channel types had short open and closed times. Mean channel open probability (P _o) for αβ was 0.54 and for αβ_T was 0.50. In comparison, αγ and αγ_T channels exhibited different kinetics: αγ channels (6.7 pS in Li⁺) had either long open times with short closings, resulting in a high P _o (0.78), or short openings with long closed times, resulting in a low P _o (0.16). The mean P _o for all αγ channels was 0.48. αγ_T (6.6 pS in Li⁺) behaved as a single population of channels with distinct kinetics: mean open time of 1.2 s and closed time of 0.4 s, with a mean P _o of 0.6, similar to that of αγ. Inclusion of 0.1 μM amiloride in the pipette solution reduced the mean open time of αγ_T to 151 ms without significantly altering the closed time. We also examined the kinetics of amiloride block of αβ, αβ_T (1 μM amiloride), and αγ_T (0.1 μM amiloride) channels. αβ and αβ_T had similar blocking and unblocking rate constants, whereas the unblocking rate constant for αγ_T was 10-fold slower than αβ_T. Our results indicate that subunit composition of ENaC is a main determinant of P _o. In addition, channel kinetics and P _o are not altered by carboxy-terminal deletion in the β subunit, whereas a similar deletion in the γ subunit affects channel kinetics but not P _o.     

Differences in Hematological Traits between High- and Low-Altitude Lizards (Genus Phrynocephalus)

Phrynocephalus erythrurus (Lacertilia: Agamidae) is considered to be the highest living reptile in the world (about 4500-5000 m above sea level), whereas Phrynocephalus przewalskii inhabits low altitudes (about 1000-1500 m above sea level). Here, we report the differences in hematological traits between these two different Phrynocephalus species. Compared with P. przewalskii, the results indicated that P. erythrurus own higher oxygen carrying capacity by increasing red blood cell count (RBC), hemoglobin concentration ([Hb]) and hematocrit (Hct) and these elevations could promote oxygen carrying capacity without disadvantage of high viscosity. The lower partial pressure of oxygen in arterial blood (PaO₂) of P. erythrurus did not cause the secondary alkalosis, which may be attributed to an efficient pulmonary system for oxygen (O₂) loading. The elevated blood-O₂ affinity in P. erythrurus may be achieved by increasing intrinsic O₂ affinity of isoHbs and balancing the independent effects of potential heterotropic ligands. We detected one α-globin gene and three β-globin genes with 1 and 33 amino acid substitutions between these two species, respectively. Molecular dynamics simulation results showed that amino acids substitutions in β-globin chains could lead to the elimination of hydrogen bonds in T-state Hb models of P. erythrurus. Based on the present data, we suggest that P. erythrurus have evolved an efficient oxygen transport system under the unremitting hypobaric hypoxia.     

Regulation of Endothelial Nitric-oxide Synthase (NOS) S-Glutathionylation by Neuronal NOS: EVIDENCE OF A FUNCTIONAL INTERACTION BETWEEN MYOCARDIAL CONSTITUTIVE NOS ISOFORMS*

Myocardial constitutive No production depends on the activity of both endothelial and neuronal NOS (eNOS and nNOS, respectively). Stimulation of myocardial β₃-adrenergic receptor (β₃-AR) produces a negative inotropic effect that is dependent on eNOS. We evaluated whether nNOS also plays a role in β₃-AR signaling and found that the β₃-AR-mediated reduction in cell shortening and [Ca²⁺]_i transient amplitude was abolished both in eNOS^−/− and nNOS^−/− left ventricular (LV) myocytes and in wild type LV myocytes after nNOS inhibition with S-methyl-l-thiocitrulline. LV superoxide (O₂^˙̄) production was increased in nNOS^−/− mice and reduced by l-N^ω-nitroarginine methyl ester (l-NAME), indicating uncoupling of eNOS activity. eNOS S-glutathionylation and Ser-1177 phosphorylation were significantly increased in nNOS^−/− myocytes, whereas myocardial tetrahydrobiopterin, eNOS Thr-495 phosphorylation, and arginase activity did not differ between genotypes. Although inhibitors of xanthine oxidoreductase (XOR) or NOX2 NADPH oxidase caused a similar reduction in myocardial O₂^˙̄, only XOR inhibition reduced eNOS S-glutathionylation and Ser-1177 phosphorylation and restored both eNOS coupled activity and the negative inotropic and [Ca²⁺]_i transient response to β₃-AR stimulation in nNOS^−/− mice. In summary, our data show that increased O₂^˙̄ production by XOR selectively uncouples eNOS activity and abolishes the negative inotropic effect of β₃-AR stimulation in nNOS^−/− myocytes. These findings provide unequivocal evidence of a functional interaction between the myocardial constitutive NOS isoforms and indicate that aspects of the myocardial phenotype of nNOS^−/− mice result from disruption of eNOS signaling.     

Selectivity of Digitalis Glycosides for Isoforms of Human Na,K-ATPase

There are four isoforms of the α subunit (α1–4) and three isoforms of the β subunit (β1–3) of Na,K-ATPase, with distinct tissue-specific distribution and physiological functions. α2 is thought to play a key role in cardiac and smooth muscle contraction and be an important target of cardiac glycosides. An α2-selective cardiac glycoside could provide important insights into physiological and pharmacological properties of α2. The isoform selectivity of a large number of cardiac glycosides has been assessed utilizing α1β1, α2β1, and α3β1 isoforms of human Na,K-ATPase expressed in Pichia pastoris and the purified detergent-soluble isoform proteins. Binding affinities of the digitalis glycosides, digoxin, β-methyl digoxin, and digitoxin show moderate but highly significant selectivity (up to 4-fold) for α2/α3 over α1 (K_D α1 > α2 = α3). By contrast, ouabain shows moderate selectivity (≈2.5-fold) for α1 over α2 (K_D α1 ≤ α3 < α2). Binding affinities for the three isoforms of digoxigenin, digitoxigenin, and all other aglycones tested are indistinguishable (K_D α1 = α3 = α2), showing that the sugar determines isoform selectivity. Selectivity patterns for inhibition of Na,K-ATPase activity of the purified isoform proteins are consistent with binding selectivities, modified somewhat by different affinities of K⁺ ions for antagonizing cardiac glycoside binding on the three isoforms. The mechanistic insight on the role of the sugars is strongly supported by a recent structure of Na,K-ATPase with bound ouabain, which implies that aglycones of cardiac glycosides cannot discriminate between isoforms. In conclusion, several digitalis glycosides, but not ouabain, are moderately α2-selective. This supports a major role of α2 in cardiac contraction and cardiotonic effects of digitalis glycosides.     

Importance of the Voltage Dependence of Cardiac Na/K ATPase Isozymes

Cardiac cells express more than one isoform of the Na, K-ATPase (NKA), the heteromeric enzyme that creates the Na⁺ and K⁺ gradients across the plasmalemma. Cardiac isozymes contain one catalytic α-subunit isoform (α1, α2, or α3) associated with an auxiliary β-subunit isoform (β1 or β2). Past studies using biochemical approaches have revealed minor kinetic differences between isozymes formed by different α-β isoform combinations; these results make it difficult to understand the physiological requirement for multiple isoforms. In intact cells, however, NKA enzymes operate in a more complex environment, which includes a substantial transmembrane potential. We evaluated the voltage dependence of human cardiac NKA isozymes expressed in Xenopus oocytes, and of native NKA isozymes in rat ventricular myocytes, using normal mammalian physiological concentrations of Na⁺_o and K⁺_o. We demonstrate that although α1 and α3 pumps are functional at all physiologically relevant voltages, α2β1 pumps and α2β2 pumps are inhibited by ∼75% and ∼95%, respectively, at resting membrane potentials, and only activate appreciably upon depolarization. Furthermore, phospholemman (FXYD1) inhibits pump function without significantly altering the pump’s voltage dependence. Our observations provide a simple explanation for the physiological relevance of the α2 subunit (∼20% of total α subunits in rat ventricle): they act as a reserve and are recruited into action for extra pumping during the long-lasting cardiac action potential, where most of the Na⁺ entry occurs. This strong voltage dependence of α2 pumps also helps explain how cardiotonic steroids, which block NKA pumps, can be a beneficial treatment for heart failure: by only inhibiting the α2 pumps, they selectively reduce NKA activity during the cardiac action potential, leading to an increase in systolic Ca²⁺, due to reduced extrusion through the Na/Ca exchanger, without affecting resting Na⁺ and Ca²⁺ concentrations.     

Differential α4(+)/(?)β2 Agonist-binding Site Contributions to α4β2 Nicotinic Acetylcholine Receptor Function within and between Isoforms

Two α4β2 nicotinic acetylcholine receptor (α4β2-nAChR) isoforms exist with (α4)₂(β2)₃ and (α4)₃(β2)₂ subunit stoichiometries and high versus low agonist sensitivities (HS and LS), respectively. Both isoforms contain a pair of α4(+)/(−)β2 agonist-binding sites. The LS isoform also contains a unique α4(+)/(−)α4 site with lower agonist affinity than the α4(+)/(−)β2 sites. However, the relative roles of the conserved α4(+)/(−)β2 agonist-binding sites in and between the isoforms have not been studied. We used a fully linked subunit concatemeric nAChR approach to express pure populations of HS or LS isoform α4β2*-nAChR. This approach also allowed us to mutate individual subunit interfaces, or combinations thereof, on each isoform background. We used this approach to systematically mutate a triplet of β2 subunit (−)-face E-loop residues to their non-conserved α4 subunit counterparts or vice versa (β2HQT and α4VFL, respectively). Mutant-nAChR constructs (and unmodified controls) were expressed in Xenopus oocytes. Acetylcholine concentration-response curves and maximum function were measured using two-electrode voltage clamp electrophysiology. Surface expression was measured with ¹²⁵I-mAb 295 binding and was used to define function/nAChR. If the α4(+)/(−)β2 sites contribute equally to function, making identical β2HQT substitutions at either site should produce similar functional outcomes. Instead, highly differential outcomes within the HS isoform, and between the two isoforms, were observed. In contrast, α4VFL mutation effects were very similar in all positions of both isoforms. Our results indicate that the identity of subunits neighboring the otherwise equivalent α4(+)/(−)β2 agonist sites modifies their contributions to nAChR activation and that E-loop residues are an important contributor to this neighbor effect.     

Stabilization of the ADP/Metaphosphate Intermediate during ATP Hydrolysis in Pre-power Stroke Myosin: QUANTITATIVE ANATOMY OF AN ENZYME*

It has been proposed recently that ATP hydrolysis in ATPase enzymes proceeds via an initial intermediate in which the dissociated γ-phosphate of ATP is bound in the protein as a metaphosphate (P_γO₃⁻). A combined quantum/classical analysis of this dissociated nucleotide state inside myosin provides a quantitative understanding of how the enzyme stabilizes this unusual metaphosphate. Indeed, in vacuum, the energy of the ADP³⁻·P_γO₃⁻·Mg²⁺ complex is much higher than that of the undissociated ATP⁴⁻. The protein brings it to a surprisingly low value. Energy decomposition reveals how much each interaction in the protein stabilizes the metaphosphate state; backbone peptides of the P-loop contribute 50% of the stabilization energy, and the side chain of Lys-185⁺ contributes 25%. This can be explained by the fact that these groups make strong favorable interactions with the α- and β-phosphates, thus favoring the charge distribution of the metaphosphate state over that of the ATP state. Further stabilization (16%) is achieved by a hydrogen bond between the backbone C=O of Ser-237 (on loop Switch-1) and a water molecule perfectly positioned to attack the P_γO₃⁻ in the subsequent hydrolysis step. The planar and singly negative P_γO₃⁻ is a much better target for the subsequent nucleophilic attack by a negatively charged OH⁻ than the tetrahedral and doubly negative P_γO₄²⁻ group of ATP. Therefore, we argue that the present mechanism of metaphosphate stabilization is common to the large family of nucleotide-hydrolyzing enzymes. Methodologically, this work presents a computational approach that allows us to obtain a truly quantitative conception of enzymatic strategy.     

N-terminal Isoforms of the Large-conductance Ca2+-activated K+ Channel Are Differentially Modulated by the Auxiliary β1-Subunit

The large-conductance Ca²⁺-activated K⁺ (BK_Ca) channel is essential for maintaining the membrane in a hyperpolarized state, thereby regulating neuronal excitability, smooth muscle contraction, and secretion. The BK_Ca α-subunit has three predicted initiation codons that generate proteins with N-terminal ends starting with the amino acid sequences MANG, MSSN, or MDAL. Because the N-terminal region and first transmembrane domain of the α-subunit are required for modulation by auxiliary β1-subunits, we examined whether β1 differentially modulates the N-terminal BK_Ca α-subunit isoforms. In the absence of β1, all isoforms had similar single-channel conductances and voltage-dependent activation. However, whereas β1 did not modulate the voltage-activation curve of MSSN, β1 induced a significant leftward shift of the voltage activation curves of both the MDAL and MANG isoforms. These shifts, of which the MDAL was larger, occurred at both 10 μm and 100 μm Ca²⁺. The β1-subunit increased the open dwell times of all three isoforms and decreased the closed dwell times of MANG and MDAL but increased the closed dwell times of MSSN. The distinct modulation of voltage activation by the β1-subunit may be due to the differential effect of β1 on burst duration and interburst intervals observed among these isoforms. Additionally, we observed that the related β2-subunit induced comparable leftward shifts in the voltage-activation curves of all three isoforms, indicating that the differential modulation of these isoforms was specific to β1. These findings suggest that the relative expression of the N-terminal isoforms can fine-tune BK_Ca channel activity in cells, highlighting a novel mechanism of BK_Ca channel regulation.     

Characteristics of an R-Phycoerythrin with Two γ Subunits Prepared from Red Macroalga Polysiphonia urceolata

An R-phycoerythrin (R-PE) was isolated by gel filtrations on Sepharose CL-4B and Sephadex G-150 from the phycobiliprotein extract of the marine red macroalga Polysiphonia urceolata Grev and further purified by ion exchange chromatography on DEAE-Sepharose Fast Flow. The purified R-PE showed three absorption peaks at 498 nm, 538 nm, 566 nm and one fluorescent emission maximum at 577 nm. Although the R-PE showed a single band on the examination by native PAGE, it exhibited two very close bands at pH about 4.7 in native isoelectric focusing (IEF). Polypeptide analysis of the R-PE demonstrated that it contained four chromophore-carrying subunits, α^18.2, β^20.6, γ^31.6 (γ''), γ^34.6 (γ), and no colorless polypeptide; its subunit composition was 6α^18.2:6β^20.6:1 γ^31.6:2γ^34.6. The α and β subunits were distributed within a acidic pH range from 5.0 to 6.0 in denaturing IEF and the γ subunits were in a basic pH range from 7.6 to 8.1. These results reveal that the prepared R-PE may exist in two hexamers of γ (αβ)₃ γ (αβ)₃γ'' and γ (αβ)₃ γ''(αβ)₃ γ and that the R-PE participate in the rod domain assembly of P. urceolata phycobilisomes by stacking each of its trimer (αβ)₃ face-to-face with the aid of one γ subunit (γ or γ'').     

Mechanisms of Attenuation of Pulmonary V’O2 Slow Component in Humans after Prolonged Endurance Training

In this study we have examined the effect of prolonged endurance training program on the pulmonary oxygen uptake (V’O₂) kinetics during heavy-intensity cycling-exercise and its impact on maximal cycling and running performance. Twelve healthy, physically active men (mean±SD: age 22.33±1.44 years, V’O_2peak 3198±458 mL ∙ min^-1) performed an endurance training composed mainly of moderate-intensity cycling, lasting 20 weeks. Training resulted in a decrease (by ~5%, P = 0.027) in V’O₂ during prior low-intensity exercise (20 W) and in shortening of τ_p of the V’O₂ on-kinetics (30.1±5.9 s vs. 25.4±1.5 s, P = 0.007) during subsequent heavy-intensity cycling. This was accompanied by a decrease of the slow component of V’O₂ on-kinetics by 49% (P = 0.001) and a decrease in the end-exercise V’O₂ by ~5% (P = 0.005). An increase (P = 0.02) in the vascular endothelial growth factor receptor 2 mRNA level and a tendency (P = 0.06) to higher capillary-to-fiber ratio in the vastus lateralis muscle were found after training (n = 11). No significant effect of training on the V’O_2peak was found (P = 0.12). However, the power output reached at the lactate threshold increased by 19% (P = 0.01). The power output obtained at the V’O_2peak increased by 14% (P = 0.003) and the time of 1,500-m performance decreased by 5% (P = 0.001). Computer modeling of the skeletal muscle bioenergetic system suggests that the training-induced decrease in the slow component of V’O₂ on-kinetics found in the present study is mainly caused by two factors: an intensification of the each-step activation (ESA) of oxidative phosphorylation (OXPHOS) complexes after training and decrease in the ‘‘additional” ATP usage rising gradually during heavy-intensity exercise.     

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θ_sil=∼0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima''s D=0.316 vs D=−0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (π_tot=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4N_e generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role.