Use of intrinsic modes in biology: Examples of indicial response of pulmonary blood pressure to ± step hypoxia

Recently, a new method to analyze biological nonstationary stochastic variables has been presented. The method is especially suitable to analyze the variation of one biological variable with respect to changes of another variable. Here, it is illustrated by the change of the pulmonary blood pressure in response to a step change of oxygen concentration in the gas that an animal breathes. The pressure signal is resolved into the sum of a set of oscillatory intrinsic mode functions, which have zero “local mean,” and a final nonoscillatory mode. With this device, we obtain a set of “mean trends,” each of which represents a “mean” in a definitive sense, and together they represent the mean trend systematically with different degrees of oscillatory content. Correspondingly, the oscillatory content of the signal about any mean trend can be represented by a set of partial sums of intrinsic mode functions. When the concept of “indicial response function” is used to describe the change of one variable in response to a step change of another variable, we now have a set of indicial response functions of the mean trends and another set of indicial response functions to describe the energy or intensity of oscillations about each mean trend. Each of these can be represented by an analytic function whose coefficients can be determined by a least-squares curve-fitting procedure. In this way, experimental results are stated sharply by analytic functions.

Light-induced exposure of the cytoplasmic end of transmembrane helix seven in rhodopsin

A key step in signal transduction in the visual cell is the light-induced conformational change of rhodopsin that triggers the binding and activation of the guanine nucleotide-binding protein. Site-directed mAbs against bovine rhodopsin were produced and used to detect and characterize these conformational changes upon light activation. Among several antibodies that bound exclusively to the light-activated state, an antibody (IgG subclass) with the highest affinity (Ka ≈ 6 × 10−9 M) was further purified and characterized. The epitope of this antibody was mapped to the amino acid sequence 304–311. This epitope extends from the central region to the cytoplasmic end of the seventh transmembrane helix and incorporates a part of a highly conserved NPXXY motif, a critical region for signaling and agonist-induced internalization of several biogenic amine and peptide receptors. In the dark state, no binding of the antibody to rhodopsin was detected. Accessibility of the epitope to the antibody correlated with formation of the metarhodopsin II photointermediate and was reduced significantly at the metarhodopsin III intermediate. Further, incubation of the antigen–antibody complex with 11-cis-retinal failed to regenerate the native rhodopsin chromophore. These results suggest significant and reversible conformational changes in close proximity to the cytoplasmic end of the seventh transmembrane helix of rhodopsin that might be important for folding and signaling.

Yeast chorismate mutase in the R state: Simulations of the active site

The isomerization of chorismate to prephenate by chorismate mutase in the biosynthetic pathway that forms Tyr and Phe involves C5—O (ether) bond cleavage and C1—C9 bond formation in a Claisen rearrangement. Development of negative charge on the ether oxygen, stabilized by Lys-168 and Glu-246, is inferred from the structure of a complex with a transition state analogue (TSA) and from the pH-rate profile of the enzyme and the E246Q mutant. These studies imply a protonated Glu-246 well above pH 7. Here, several 500-ps molecular dynamics simulations test the stability of enzyme–TSA complexes by using a solvated system with stochastic boundary conditions. The simulated systems are (i) protonated Glu-246 (stable), (ii) deprotonated Glu-246 (unstable), (iii) deprotonated Glu-246 plus one H2O between Glu-246 and the ether oxygen (unstable), (iv) the E246Q mutant (stable), and (v) addition of OH− between protonated Glu-246 and the ether oxygen. In (v), a local conformational change of Lys-168 displaced the OH− into the solvent region, suggesting a possible rate-determining step that precedes the catalytic step. In a 500-ps simulation of the enzyme complexed with the reactant chorismate or the product prephenate, no water molecule remained near the oxygen of the ligand. Calculations using the linearized Poisson–Boltzmann equation show that the effective pKa of Glu-246 is shifted from 5.8 to 8.1 as the negative charge on the ether oxygen of the TSA is changed from −0.56 electron to −0.9 electron. Altogether, these results support retention of a proton on Glu-246 to high pH and the absence of a water molecule in the catalytic steps.

The anticoagulant activation of antithrombin by heparin

Antithrombin, a plasma serpin, is relatively inactive as an inhibitor of the coagulation proteases until it binds to the heparan side chains that line the microvasculature. The binding specifically occurs to a core pentasaccharide present both in the heparans and in their therapeutic derivative heparin. The accompanying conformational change of antithrombin is revealed in a 2.9-Å structure of a dimer of latent and active antithrombins, each in complex with the high-affinity pentasaccharide. Inhibitory activation results from a shift in the main sheet of the molecule from a partially six-stranded to a five-stranded form, with extrusion of the reactive center loop to give a more exposed orientation. There is a tilting and elongation of helix D with the formation of a 2-turn helix P between the C and D helices. Concomitant conformational changes at the heparin binding site explain both the initial tight binding of antithrombin to the heparans and the subsequent release of the antithrombin–protease complex into the circulation. The pentasaccharide binds by hydrogen bonding of its sulfates and carboxylates to Arg-129 and Lys-125 in the D-helix, to Arg-46 and Arg-47 in the A-helix, to Lys-114 and Glu-113 in the P-helix, and to Lys-11 and Arg-13 in a cleft formed by the amino terminus. This clear definition of the binding site will provide a structural basis for developing heparin analogues that are more specific toward their intended target antithrombin and therefore less likely to exhibit side effects.

Tyrosine Phosphorylation Regulates Cell Cycle-dependent Nuclear Localization of Cdc48p

Cdc48p from Saccharomyces cerevisiae and its highly conserved mammalian homologue VCP (valosin-containing protein) are ATPases with essential functions in cell division and homotypic fusion of endoplasmic reticulum vesicles. Both are mainly attached to the endoplasmic reticulum, but relocalize in a cell cycle-dependent manner: Cdc48p enters the nucleus during late G1; VCP aggregates at the centrosome during mitosis. The nuclear import signal sequence of Cdc48p was localized near the amino terminus and its function demonstrated by mutagenesis. The nuclear import is regulated by a cell cycle-dependent phosphorylation of a tyrosine residue near the carboxy terminus. Two-hybrid studies indicate that the phosphorylation results in a conformational change of the protein, exposing the nuclear import signal sequence previously masked by a stretch of acidic residues.

TGF-β–induced Phosphorylation of Smad3 Regulates Its Interaction with Coactivator p300/CREB-binding Protein

Smads are intermediate effector proteins that transduce the TGF-β signal from the plasma membrane to the nucleus, where they participate in transactivation of downstream target genes. We have shown previously that coactivators p300/CREB-binding protein are involved in TGF-β–mediated transactivation of two Cdk inhibitor genes, p21 and p15. Here we examined the possibility that Smads function to regulate transcription by directly interacting with p300/CREB-binding protein. We show that Smad3 can interact with a C-terminal fragment of p300 in a temporal and phosphorylation-dependent manner. TGF-β–mediated phosphorylation of Smad3 potentiates the association between Smad3 and p300, likely because of an induced conformational change that removes the autoinhibitory interaction between the N- and C-terminal domains of Smad3. Consistent with a role for p300 in the transcription regulation of multiple genes, overexpression of a Smad3 C-terminal fragment causes a general squelching effect on multiple TGF-β–responsive reporter constructs. The adenoviral oncoprotein E1A can partially block Smad-dependent transcriptional activation by directly competing for binding to p300. Taken together, these findings define a new role for phosphorylation of Smad3: in addition to facilitating complex formation with Smad4 and promoting nuclear translocation, the phosphorylation-induced conformational change of Smad3 modulates its interaction with coactivators, leading to transcriptional regulation.

Molecular heterochrony in the early development of Drosophila

Heterochrony, the relative change of developmental timing, is one of the major modes of macroevolutionary change; it identifies temporally disassociated units of developmental evolution. Here, we report the results of a fine-scale temporal study for the expression of the developmental gene hairy and morphological development in three species of Drosophila, D. melanogaster, D. simulans, and D. pseudoobscura. The results suggest that between and among closely related species, temporal displacement of ontogenetic trajectory is detected even at the earliest stage of development. Overall, D. simulans shows the earliest expression, followed by D. melanogaster, and then by D. pseudoobscura. Setting D. melanogaster as the standard, we find the approximate time to full expression is accelerated by 13 min, 48 s in D. simulans and retarded by 24 min in D. pseudoobscura. Morphologically, again with D. melanogaster setting the standard, initiation of cellularization is faster in D. simulans by 15 min, 42 s; and initiation of morphogenesis is faster in D. simulans by 18 min, 7 s. These results seem to be consistent with the finding that the approximate time to full expression of hairy is accelerated by 13 min, 48 s in D. simulans. On the other hand, the same morphological events are delayed by 5 min, 32 s, and by 11 min, 32 s, respectively, in D. pseudoobscura. These delays are small, compared with the 24-min delay in full expression. The timing changes, in total, seem consistent with continuous phyletic evolution of temporal trajectories. Finally, we speculate that epigenetic interactions of hairy expression timing and cell-cycle timing may have led to morphological differences in the terminal system of the larvae.

Nonlinear indicial response of complex nonstationary oscillations as pulmonary hypertension responding to step hypoxia

This paper is devoted to the quantization of the degree of nonlinearity of the relationship between two biological variables when one of the variables is a complex nonstationary oscillatory signal. An example of the situation is the indicial responses of pulmonary blood pressure (P) to step changes of oxygen tension (ΔpO2) in the breathing gas. For a step change of ΔpO2 beginning at time t1, the pulmonary blood pressure is a nonlinear function of time and ΔpO2, which can be written as P(t-t1 | ΔpO2). An effective method does not exist to examine the nonlinear function P(t-t1 | ΔpO2). A systematic approach is proposed here. The definitions of mean trends and oscillations about the means are the keys. With these keys a practical method of calculation is devised. We fit the mean trends of blood pressure with analytic functions of time, whose nonlinearity with respect to the oxygen level is clarified here. The associated oscillations about the mean can be transformed into Hilbert spectrum. An integration of the square of the Hilbert spectrum over frequency yields a measure of oscillatory energy, which is also a function of time, whose mean trends can be expressed by analytic functions. The degree of nonlinearity of the oscillatory energy with respect to the oxygen level also is clarified here. Theoretical extension of the experimental nonlinear indicial functions to arbitrary history of hypoxia is proposed. Application of the results to tissue remodeling and tissue engineering of blood vessels is discussed.

A cyclin-dependent kinase-activating kinase regulates differentiation of root initial cells in Arabidopsis

Cell division and differentiation continue throughout the plant life cycle without significant loss of control. However, little is known about the mechanisms that allow the continuous development of meristems. Cell division is controlled by a family of cyclin-dependent kinases (CDKs). CDK-activating kinases (CAKs) are known to phosphorylate and activate almost all CDKs and thus may have a crucial role in controlling CDK activities in each cell of the meristems. Here, we show that overexpression of sense or antisense gene for Cak1At in Arabidopsis by using the glucocorticoid-mediated transcriptional induction system resulted in a reduction of CDK activities. After 14–24 h of glucocorticoid treatment, starch granules appeared in columellar initials in the root meristem, and cortical initials were periclinally divided into cortical and endodermal cells. Accumulation of the cyclin∷β-glucuronidase fusion protein ceased after 72 h of glucocorticoid treatment. Our results indicate that a change of Cak1At activity leads to differentiation of initial cells, followed by cessation of cell division. Therefore, we propose that differentiation of initial cells is controlled by Cak1At but is maintained independent of cell division.

A scrapie-like unfolding intermediate of the prion protein domain PrP(121–231) induced by acidic pH

The infectious agent of transmissible spongiform encephalopathies is believed to consist of an oligomeric isoform, PrPSc, of the monomeric cellular prion protein, PrPC. The conversion of PrPC to PrPSc is characterized by a decrease in α-helical structure, an increase in β-sheet content, and the formation of PrPSc amyloid. Whereas the N-terminal part of PrPC comprising residues 23–120 is flexibly disordered, its C-terminal part, PrP(121–231), forms a globular domain with three α-helices and a small β-sheet. Because the segment of residues 90–231 is protease-resistant in PrPSc, it is most likely structured in the PrPSc form. The conformational change of the segment containing residues 90–120 thus constitutes the minimal structural difference between PrPC and a PrPSc monomer. To test whether PrP(121–231) is also capable to undergo conformational transitions, we analyzed its urea-dependent unfolding transitions at neutral and acidic pH. We identified an equilibrium unfolding intermediate of PrP(121–231) that is exclusively populated at acidic pH and shows spectral characteristics of a β-sheet protein. The intermediate is in rapid equilibrium with native PrP(121–231), significantly populated in the absence of urea at pH 4.0, and may have important implications for the presumed formation of PrPSc during endocytosis.

Photochemically induced nuclear spin polarization in reaction centers of photosystem II observed by 13C-solid-state NMR reveals a strongly asymmetric electronic structure of the P680.+ primary donor chlorophyll

We report 13C magic angle spinning NMR observation of photochemically induced dynamic nuclear spin polarization (photo- CIDNP) in the reaction center (RC) of photosystem II (PS2). The light-enhanced NMR signals of the natural abundance 13C provide information on the electronic structure of the primary electron donor P680 (chlorophyll a molecules absorbing around 680 nm) and on the pz spin density pattern in its oxidized form, P680⨥. Most centerband signals can be attributed to a single chlorophyll a (Chl a) cofactor that has little interaction with other pigments. The chemical shift anisotropy of the most intense signals is characteristic for aromatic carbon atoms. The data reveal a pronounced asymmetry of the electronic spin density distribution within the P680⨥. PS2 shows only a single broad and intense emissive signal, which is assigned to both the C-10 and C-15 methine carbon atoms. The spin density appears shifted toward ring III. This shift is remarkable, because, for monomeric Chl a radical cations in solution, the region of highest spin density is around ring II. It leads to a first hypothesis as to how the planet can provide itself with the chemical potential to split water and generate an oxygen atmosphere using the Chl a macroaromatic cycle. A local electrostatic field close to ring III can polarize the electronic charge and associated spin density and increase the redox potential of P680 by stabilizing the highest occupied molecular orbital, without a major change of color. This field could be produced, e.g., by protonation of the keto group of ring V. Finally, the radical cation electronic structure in PS2 is different from that in the bacterial RC, which shows at least four emissive centerbands, indicating a symmetric spin density distribution over the entire bacteriochlorophyll macrocycle.

Glucose-regulated interaction of a regulatory subunit of protein phosphatase 1 with the Snf1 protein kinase in Saccharomyces cerevisiae

The Snf1 protein kinase family has been conserved in eukaryotes. In the yeast Saccharomyces cerevisiae, Snf1 is essential for transcription of glucose-repressed genes in response to glucose starvation. The direct interaction between Snf1 and its activating subunit, Snf4, within the kinase complex is regulated by the glucose signal. Glucose inhibition of the Snf1-Snf4 interaction depends on protein phosphatase 1 and its targeting subunit, Reg1. Here we show that Reg1 interacts with the Snf1 catalytic domain in the two-hybrid system. This interaction increases in response to glucose limitation and requires the conserved threonine in the activation loop of the kinase, a putative phosphorylation site. The inhibitory effect of Reg1 appears to require the Snf1 regulatory domain because a reg1Δ mutation no longer relieves glucose repression of transcription when Snf1 function is provided by the isolated catalytic domain. Finally, we show that abolishing the Snf1 catalytic activity by mutation of the ATP-binding site causes elevated, constitutive interaction with Reg1, indicating that Snf1 negatively regulates its own interaction with Reg1. We propose a model in which protein phosphatase 1, targeted by Reg1, facilitates the conformational change of the kinase complex from its active state to the autoinhibited state.

Structure of the recombinant full-length hamster prion protein PrP(29–231): The N terminus is highly flexible

The prion diseases seem to be caused by a conformational change of the prion protein (PrP) from the benign cellular form PrPC to the infectious scrapie form PrPSc; thus, detailed information about PrP structure may provide essential insights into the mechanism by which these diseases develop. In this study, the secondary structure of the recombinant Syrian hamster PrP of residues 29–231 [PrP(29–231)] is investigated by multidimensional heteronuclear NMR. Chemical shift index analysis and nuclear Overhauser effect data show that PrP(29–231) contains three helices and possibly one short β-strand. Most striking is the random-coil nature of chemical shifts for residues 30–124 in the full-length PrP. Although the secondary structure elements are similar to those found in mouse PrP fragment PrP(121–231), the secondary structure boundaries of PrP(29–231) are different from those in mouse PrP(121–231) but similar to those found in the structure of Syrian hamster PrP(90–231). Comparison of resonance assignments of PrP(29–231) and PrP(90–231) indicates that there may be transient interactions between the additional residues and the structured core. Backbone dynamics studies done by using the heteronuclear [1H]-15N nuclear Overhauser effect indicate that almost half of PrP(29–231), residues 29–124, is highly flexible. This plastic region could feature in the conversion of PrPC to PrPSc by template-assisted formation of β-structure.

Locking in alternate conformations of the integrin αLβ2 I domain with disulfide bonds reveals functional relationships among integrin domains

We used integrin αLβ2 heterodimers containing I domains locked open (active) or closed (inactive) with disulfide bonds to investigate regulatory interactions among domains in integrins. mAbs to the αL I domain and β2 I-like domain inhibit adhesion of wild-type αLβ2 to intercellular adhesion molecule-1. However, with αLβ2 containing a locked open I domain, mAbs to the I domain were subdivided into subsets (i) that did not inhibit, and thus appear to inhibit by favoring the closed conformation, and (ii) that did inhibit, and thus appear to bind to the ligand binding site. Furthermore, αLβ2 containing a locked open I domain was completely resistant to inhibition by mAbs to the β2 I-like domain, but became fully susceptible to inhibition after disulfide reduction with DTT. This finding suggests that the I-like domain indirectly contributes to ligand binding by regulating opening of the I domain in wild-type αLβ2. Conversely, locking the I domain closed partially restrained conformational change of the I-like domain by Mn2+, as measured with mAb m24, which we map here to the β2 I-like domain. By contrast, locking the I domain closed or open did not affect constitutive or Mn2+-induced exposure of the KIM127 epitope in the β2 stalk region. Furthermore, locked open I domains, in αLβ2 complexes or expressed in isolation on the cell surface, bound to intercellular adhesion molecule-1 equivalently in Mg2+ and Mn2+. These results suggest that Mn2+ activates αLβ2 by binding to a site other than the I domain, most likely the I-like domain of β2.

In vitro assembly of phytochrome B apoprotein with synthetic analogs of the phytochrome chromophore

Phytochrome B (PhyB), one of the major photosensory chromoproteins in plants, mediates a variety of light-responsive developmental processes in a photoreversible manner. To analyze the structural requirements of the chromophore for the spectral properties of PhyB, we have designed and chemically synthesized 20 analogs of the linear tetrapyrrole (bilin) chromophore and reconstituted them with PhyB apoprotein (PHYB). The A-ring acts mainly as the anchor for ligation to PHYB, because the modification of the side chains at the C2 and C3 positions did not significantly influence the formation or difference spectra of adducts. In contrast, the side chains of the B- and C-rings are crucial to position the chromophore properly in the chromophore pocket of PHYB and for photoreversible spectral changes. The side-chain structure of the D-ring is required for the photoreversible spectral change of the adducts. When methyl and ethyl groups at the C17 and C18 positions are replaced with an n-propyl, n-pentyl, or n-octyl group, respectively, the photoreversible spectral change of the adducts depends on the length of the side chains. From these studies, we conclude that each pyrrole ring of the linear tetrapyrrole chromophore plays a different role in chromophore assembly and the photochromic properties of PhyB.