The platelet protein kinase C substrate pleckstrin binds directly to SDPR protein

Pleckstrin is a modular platelet protein consisting of N- and C-terminal pleckstrin homology (PH) domains, a central dishevelled egl10 and pleckstrin (DEP) domain and a phosphorylation region. Following agonist-induced platelet stimulation, dimeric pleckstrin translocates to the plasma membrane, is phosphorylated and then monomerizes. A recent study found that pleckstrin null platelets from a knockout mouse have a defect in granule secretion, actin polymerization and aggregation. However, the mechanism of pleckstrin signaling for this function is unknown. Our recent studies have led to the identification of a novel pleckstrin-binding protein, serum deprivation response protein (SDPR), by co-immunoprecipitation, GST-pulldowns and nanospray quadruple time of flight mass spectrometry. We show that this interaction occurs directly through N-terminal sequences of pleckstrin. Both pleckstrin and SDPR are phosphorylated by protein kinase C (PKC), but the interaction between pleckstrin and SDPR was shown to be independent of PKC inhibition or activation. These results suggest that SDPR may facilitate the translocation of nonphosphorylated pleckstrin to the plasma membrane in conjunction with phosphoinositides that bind to the C-terminal PH domain. After binding of pleckstrin to the plasma membrane, its phosphorylation by PKC exerts downstream effects on platelet aggregation/secretion.

Lobe specific Ca2+-calmodulin nano-domain in neuronal spines: a single molecule level analysis.

Calmodulin (CaM) is a ubiquitous Ca(2+) buffer and second messenger that affects cellular function as diverse as cardiac excitability, synaptic plasticity, and gene transcription. In CA1 pyramidal neurons, CaM regulates two opposing Ca(2+)-dependent processes that underlie memory formation: long-term potentiation (LTP) and long-term depression (LTD). Induction of LTP and LTD require activation of Ca(2+)-CaM-dependent enzymes: Ca(2+)/CaM-dependent kinase II (CaMKII) and calcineurin, respectively. Yet, it remains unclear as to how Ca(2+) and CaM produce these two opposing effects, LTP and LTD. CaM binds 4 Ca(2+) ions: two in its N-terminal lobe and two in its C-terminal lobe. Experimental studies have shown that the N- and C-terminal lobes of CaM have different binding kinetics toward Ca(2+) and its downstream targets. This may suggest that each lobe of CaM differentially responds to Ca(2+) signal patterns. Here, we use a novel event-driven particle-based Monte Carlo simulation and statistical point pattern analysis to explore the spatial and temporal dynamics of lobe-specific Ca(2+)-CaM interaction at the single molecule level. We show that the N-lobe of CaM, but not the C-lobe, exhibits a nano-scale domain of activation that is highly sensitive to the location of Ca(2+) channels, and to the microscopic injection rate of Ca(2+) ions. We also demonstrate that Ca(2+) saturation takes place via two different pathways depending on the Ca(2+) injection rate, one dominated by the N-terminal lobe, and the other one by the C-terminal lobe. Taken together, these results suggest that the two lobes of CaM function as distinct Ca(2+) sensors that can differentially transduce Ca(2+) influx to downstream targets. We discuss a possible role of the N-terminal lobe-specific Ca(2+)-CaM nano-domain in CaMKII activation required for the induction of synaptic plasticity.

Nonclassical mechanisms of progesterone action in the brain: I. Protein kinase C activation in the hypothalamus of female rats.

The modulation of gene regulation by progesterone (P) and its classical intracellular regulation by progestin receptors in the brain, resulting in alterations in physiology and behavior has been well studied. The mechanisms mediating the short latency effects of P are less well understood. Recent studies have revealed rapid nonclassical signaling action of P involving the activation of intracellular signaling pathways. We explored the involvement of protein kinase C (PKC) in P-induced rapid signaling in the ventromedial nucleus of the hypothalamus (VMN) and preoptic area (POA) of the rat brain. Both the Ca2+-independent (basal) PKC activity representing the activation of PKC by the in vivo treatments and the Ca+2-dependent (total) PKC activity assayed in the presence of exogenous cofactors in vitro were determined. A comparison of the two activities demonstrated the strength and temporal status of PKC regulation by steroid hormones in vivo. P treatment resulted in a rapid increase in basal PKC activity in the VMN but not the POA. Estradiol benzoate priming augmented P-initiated increase in PKC basal activity in both the VMN and POA. These increases were inhibited by intracerebroventricular administration of a PKC inhibitor administered 30 min prior to P. The total PKC activity remained unchanged demonstrating maximal PKC activation within 30 min in the VMN. In contrast, P regulation in the POA significantly attenuated total PKC activity +/- estradiol benzoate priming. These rapid changes in P-initiated PKC activity were not due to changes in PKC protein levels or phosphorylation status.

Multiple site phosphorylation of tyrosine hydroxylase: A potential role for protein kinase C in the regulation of catecholamine synthesis

Tyrosine hydroxylase (TH), the initial and rate limiting enzyme in the catecholaminergic biosynthetic pathway, is phosphorylated on multiple serine residues by multiple protein kinases. Although it has been demonstrated that many protein kinases are capable of phosphorylating and activating TH in vitro, it is less clear which protein kinases participate in the physiological regulation of catecholamine synthesis in situ. These studies were designed to determine if protein kinase C (PK-C) plays such a regulatory role.^ Stimulation of intact bovine adrenal chromaffin cells with phorbol esters results in stimulation of catecholamine synthesis, tyrosine hydroxylase phosphorylation and activation. These responses are both time and concentration dependent, and are specific for those phorbol ester analogues which activate PK-C. RP-HPLC analysis of TH tryptic phosphopeptides indicate that PK-C phosphorylates TH on three putative sites. One of these (pepetide 6) is the same as that phosphorylated by both cAMP-dependent protein kinase (PK-A) and calcium/calmodulin-dependent protein kinase (CaM-K). However, two of these sites (peptides 4 and 7) are unique, and, to date, have not been shown to be phosphorylated by any other protein kinase. These peptides correspond to those which are phosphorylated with a slow time course in response to stimulation of chromaffin cells with the natural agonist acetylcholine. The activation of TH produced by PK-C is most closely correlated with the phosphorylation of peptide 6. But, as evident from pH profiles of tyrosine hydroxylase activity, phosphorylation of peptides 4 and 7 affect the expression of the activation produced by phosphorylation of peptide 6.^ These data support a role for PK-C in the control of TH activity, and suggest a two stage model for the physiological regulation of catecholamine synthesis by phosphorylation in response to cholinergic stimulation. An initial fast response, which appears to be mediated by CaM-K, and a slower, sustained response which appears to be mediated by PK-C. In addition, the multiple site phosphorylation of TH provides a mechanism whereby the regulation of catecholamine synthesis appears to be under the control of multiple protein kinases, and allows for the convergence of multiple, diverse physiological and biochemical signals. ^

Fatty acylated proteins and the protein kinase C signaling pathway

Numerous proteins in intracellular signaling pathways are known to be covalently modified by long chain fatty acids. The objective of this project was to identify potentially novel components of the protein kinase C signaling pathway by virtue of their fatty acylation. A 64 kDa palmitoylated protein (p64) was identified that became deacylated following stimulation of quiescent cells with serum, FGF, or PDBu, suggesting that stimulus-dependent deacylation might alter interactions between p64 and other membrane/cytoskeletal components. A myristoylated protein of 68 kDa observed during these studies was identified as the "80K" PKC substrate. This protein was acylated cotranslationally with myristate through an amide linkage. The majority of the 80K protein was tightly associated with the plasma membrane, with approximately 20% in the cytosol. Although phosphorylation of the membrane-bound and soluble forms of the protein was increased 6-fold in response to PDBu, no changes in the subcellular distribution or myristoylation of the protein were observed. A cDNA encoding the murine form of this protein was cloned, and its deduced amino acid sequence revealed the presence of an N-terminal myristoylation consensus and five potential sites for phosphorylation by PKC. A mutant in which the N-terminal glycine residue was changed to alanine was no longer a substrate for NMT and consequently lost its membrane-binding potential. However, its ability to be phosphorylated in response to purified growth factors and phorbol esters was unimpaired. These results indicate that the myristoylated N-terminus of the 80K protein is required for its association with the plasma membrane, and that the cytoplasmic form of the protein can be phosphorylated independently of the membrane-bound form. Mutants of PKC were constructed in which the regulatory domain was removed and replaced by the N-terminus of the 80K or Al proteins. Unexpectedly, both the myristoylated and nonmyristoylated fusion proteins were tightly associated with the nuclear envelope. Further deletion analyses mapped nuclear targeting signals to the hinge region and a portion of the catalytic domain of PKC, explaining the ability of PKC to be translocated to the nucleus in response to certain stimuli. ^

Functional and Evolutionary Analysis of the CASPARIAN STRIP MEMBRANE DOMAIN PROTEIN Family

CASPARIAN STRIP MEMBRANE DOMAIN PROTEINS (CASPs) are four-membrane-span proteins that mediate the deposition of Casparian strips in the endodermis by recruiting the lignin polymerization machinery. CASPs show high stability in their membrane domain, which presents all the hallmarks of a membrane scaffold. Here, we characterized the large family of CASP-like (CASPL) proteins. CASPLs were found in all major divisions of land plants as well as in green algae; homologs outside of the plant kingdom were identified as members of the MARVEL protein family. When ectopically expressed in the endodermis, most CASPLs were able to integrate the CASP membrane domain, which suggests that CASPLs share with CASPs the propensity to form transmembrane scaffolds. Extracellular loops are not necessary for generating the scaffold, since CASP1 was still able to localize correctly when either one of the extracellular loops was deleted. The CASP first extracellular loop was found conserved in euphyllophytes but absent in plants lacking Casparian strips, an observation that may contribute to the study of Casparian strip and root evolution. In Arabidopsis (Arabidopsis thaliana), CASPL showed specific expression in a variety of cell types, such as trichomes, abscission zone cells, peripheral root cap cells, and xylem pole pericycle cells.

Maximal acceleration of calcium release refractoriness by β-adrenergic stimulation requires dual activation of protein kinase A and CaMKII in mouse ventricular myocytes

Time-dependent refractoriness of calcium (Ca2+) release in cardiac myocytes is an important factor in determining whether pro-arrhythmic release patterns develop. At the subcellular level of the Ca2+ spark, recent studies have suggested that recovery of spark amplitude is controlled by local sarcoplasmic reticulum (SR) refilling whereas refractoriness of spark triggering depends on both refilling and the sensitivity of the ryanodine receptor (RyR) release channels that produce sparks. Here we studied regulation of Ca2+ spark refractoriness in mouse ventricular myocytes by examining how β-adrenergic stimulation influenced sequences of Ca2+ sparks originating from individual RyR clusters. Our protocol allowed us to separately measure recovery of spark amplitude and delays between successive sparks, and data were interpreted quantitatively through simulations with a stochastic mathematical model. We found that, compared with spark sequences measured under control conditions: (1) β-adrenergic stimulation with isoproterenol accelerated spark amplitude recovery and decreased spark-to-spark delays; (2) activating protein kinase A (PKA) with forskolin accelerated amplitude recovery but did not affect spark-to-spark delays; (3) inhibiting PKA with H89 retarded amplitude recovery and increased spark- to-spark delays; (4) preventing phosphorylation of the RyR at serine 2808 with a knock-in mouse prevented the decrease in spark-to-spark delays seen with β-adrenergic stimulation; (5) inhibiting either PKA or Ca2+/calmodulin-dependent protein kinase II (CaMKII) during β-adrenergic stimulation prevented the decrease in spark-to-spark delays seen) without inhibition. The results suggest that activation of either PKA or CaMKII is sufficient to speed SR refilling, but activation of both kinases appears necessary to observe increased RyR sensitivity. The data provide novel insight into β-adrenergic regulation of Ca2+ release refractoriness in mouse myocytes.

Identification of Novel Death-Associated Protein Kinase 2 Interaction Partners by Proteomic Screening Coupled with Bimolecular Fluorescence Complementation.

Death-associated protein kinase 2 (DAPK2) is a Ca(2+)/calmodulin-dependent Ser/Thr kinase that possesses tumor-suppressive functions and regulates programmed cell death, autophagy, oxidative stress, hematopoiesis, and motility. As only few binding partners of DAPK2 have been determined, the molecular mechanisms governing these biological functions are largely unknown. We report the identification of 180 potential DAPK2 interaction partners by affinity purification-coupled mass spectrometry, 12 of which are known DAPK binding proteins. A small subset of established and potential binding proteins detected in this screen was further investigated by bimolecular fluorescence complementation (BiFC) assays, a method to visualize protein interactions in living cells. These experiments revealed that α-actinin-1 and 14-3-3-β are novel DAPK2 binding partners. The interaction of DAPK2 with α-actinin-1 was localized at the plasma membrane, resulting in massive membrane blebbing and reduced cellular motility, whereas the interaction of DAPK2 with 14-3-3-β was localized to the cytoplasm, with no impact on blebbing, motility, or viability. Our results therefore suggest that DAPK2 effector functions are influenced by the protein's subcellular localization and highlight the utility of combining mass spectrometry screening with bimolecular fluorescence complementation to identify and characterize novel protein-protein interactions.

Murine coronavirus ubiquitin-like domain is important for papain-like protease stability and viral pathogenesis

Ubiquitin-like domains (Ubls) now are recognized as common elements adjacent to viral and cellular proteases; however, their function is unclear. Structural studies of the papain-like protease (PLP) domains of coronaviruses (CoVs) revealed an adjacent Ubl domain in severe acute respiratory syndrome CoV, Middle East respiratory syndrome CoV, and the murine CoV, mouse hepatitis virus (MHV). Here, we tested the effect of altering the Ubl adjacent to PLP2 of MHV on enzyme activity, viral replication, and pathogenesis. Using deletion and substitution approaches, we identified sites within the Ubl domain, residues 785 to 787 of nonstructural protein 3, which negatively affect protease activity, and valine residues 785 and 787, which negatively affect deubiquitinating activity. Using reverse genetics, we engineered Ubl mutant viruses and found that AM2 (V787S) and AM3 (V785S) viruses replicate efficiently at 37°C but generate smaller plaques than wild-type (WT) virus, and AM2 is defective for replication at higher temperatures. To evaluate the effect of the mutation on protease activity, we purified WT and Ubl mutant PLP2 and found that the proteases exhibit similar specific activities at 25°C. However, the thermal stability of the Ubl mutant PLP2 was significantly reduced at 30°C, thereby reducing the total enzymatic activity. To determine if the destabilizing mutation affects viral pathogenesis, we infected C57BL/6 mice with WT or AM2 virus and found that the mutant virus is highly attenuated, yet it replicates sufficiently to elicit protective immunity. These studies revealed that modulating the Ubl domain adjacent to the PLP reduces protease stability and viral pathogenesis, revealing a novel approach to coronavirus attenuation. IMPORTANCE Introducing mutations into a protein or virus can have either direct or indirect effects on function. We asked if changes in the Ubl domain, a conserved domain adjacent to the coronavirus papain-like protease, altered the viral protease activity or affected viral replication or pathogenesis. Our studies using purified wild-type and Ubl mutant proteases revealed that mutations in the viral Ubl domain destabilize and inactivate the adjacent viral protease. Furthermore, we show that a CoV encoding the mutant Ubl domain is unable to replicate at high temperature or cause lethal disease in mice. Our results identify the coronavirus Ubl domain as a novel modulator of viral protease stability and reveal manipulating the Ubl domain as a new approach for attenuating coronavirus replication and pathogenesis.

RNA-Activated Protein Kinase, PKR

The double-stranded RNA (dsRNA) activated protein kinase, PKR, is one of the several enzymes induced by interferons and a key molecule mediating the antiviral effects of interferons. PKR contain an N-terminal, double-stranded RNA binding domain (dsRBD), which has two tandem copies of the motifs (dsRBM I and dsRBM II). Upon binding to viral dsRNA, PKR is activated via autophosphorylation. Activated PKR has several substrates; one of the examples is eukaryotic translation initiation factor 2 (eIF2a). The phosphorylation of eIF2a leads to the termination of cell growth by inhibiting protein synthesis in response to viral infection. The objective of this project was to characterize the dsRBM I and define the dsRNA binding using biophysical methods. First, the dsRBM I gene was cloned from a pET-28b to a pET-11a expression plasmid. N-terminal poly-histidine tags on pET-28b are for affinity purification; however, these tags can alter the structure and function of proteins, thus the gene of dsRBM I was transferred into the plasmid without tags (pET-11a) and expressed as a native protein. The dsRBM I was transformed into and expressed by Rosetta DE3plyS expression cells. Purification was done by FPLC using a Sepharose IEX ion exchange followed by Heparin affinity column; yielding pure protein was assayed by PAGE. Analytical Ultracentrifugation, Sedimentation Velocity, was used to characterize free solution association state and hydrodynamic properties of the protein. The slight decrease in S-value with concentration is due to the hydrodynamic non-ideality. No self association was observed. The obtained molecule weight was 10,079 Da. The calculated sedimentation constant at zero concentration at 20°C in water was 1.23 and its friction coefficient was 3.575 ´ 10-8. The frictional ratio of sphere and dsRBM I became 1.30. Therefore, dsRBM I must be non-globular and more asymmetric shape. Isolated dsRBM I exhibits the same tertiary fold as compared to context in the full domain but it exhibited weaker binding affinity than full domain to a 20 bp dsRNA. However, when the conditions allowed for its saturation, dsRBM I to 20 bp dsRNA has similar stoichiometry as full dsRBD.

Heregulin -beta 1 directly associates with RhoA via the immunoglobulin -like domain thereby modulating breast cancer metastasis-related properties

Heregulins constitute a family of growth factors belonging to the epidermal growth factor (EGF) family. Breast cancers that overexpress specific members of the EGF receptor family (EGFR, ErbB2, ErbB3, ErbB4) have increased metastatic potential, and Heregulin-β1 (HRGβ1), a ligand for ErbB3 and ErbB4, has also been shown to induce metastasis-related properties in breast cancer cells in vitro. The secreted form of the HRGβ1 is composed of five distinct structural domains, including the N-terminal domain, an immunoglobulin-like domain (IgG-like), a glycosylation domain, an EGF-like domain, and a β1-specific domain. Of these, the EGF-like domain is well characterized for its function in metastasis-related properties as well as its structure. However, the contributions of the other HRGβ1 domains in breast cancer metastasis remains unclear. ^ To investigate this, HRGβ1 proteins with targeted domain deletions were purified and subjected to assays for metastasis-related properties, including aggregation, invasion, activation of EGFR family members, and motility of breast cancer cells. These assays showed that retaining the EGF-like domain of HRGβ1 is important for activation of EGFRs. Interestingly, the HRGβ1 protein lacking the IgG-like domain (NGEB) led to a decrease in breast cancer cell motility, indicating the IgG-like domain modulates cell motility, an important step in cancer metastasis. ^ To understand the underlying mechanisms, I performed protein sequence and structural analysis of HRGβ1 and identified that the IgG-like domain of HRGβ1 shares sequence homology and three-dimensional structural similarity with the IgG-like domain of TRIO. TRIO is a cytoplasmic protein that directly associates with RhoA, a GTPase involved in cell reorganization and cell motility. Therefore, I hypothesized that HRGβ1 may translocate inside the breast cancer cells through receptor mediated endocytosis and bind to RhoA via its IgG-like domain. I show wild type HRGβ1 but not NGEB binds RhoA in vitro and in vivo, leading to RhoA activation. Inhibition of HRG-β1 internalization via endocytosis disrupted HRGβ1 binding to RhoA. Additionally, breast cancer cell motility induced by HRG-β1 is reduced after treatment with inhibitors to both endocytosis and RhoA function, similar to levels seen with NGEB treatment. ^ Thus, in addition to the well-known role of HRGβ1 as an extracellular stimulator of the EGFR family members, HRGβ1 also functions within the cell as a binding partner and activator of RhoA to modulate cancer cell motility. ^

Activating mutations in the extracellular domain of the fibroblast growth factor receptor 2 function by disruption of the disulfide bond in the third immunoglobulin-like domain

Multiple human skeletal and craniosynostosis disorders, including Crouzon, Pfeiffer, Jackson–Weiss, and Apert syndromes, result from numerous point mutations in the extracellular region of fibroblast growth factor receptor 2 (FGFR2). Many of these mutations create a free cysteine residue that potentially leads to abnormal disulfide bond formation and receptor activation; however, for noncysteine mutations, the mechanism of receptor activation remains unclear. We examined the effect of two of these mutations, W290G and T341P, on receptor dimerization and activation. These mutations resulted in cellular transformation when expressed as FGFR2/Neu chimeric receptors. Additionally, in full-length FGFR2, the mutations induced receptor dimerization and elevated levels of tyrosine kinase activity. Interestingly, transformation by the chimeric receptors, dimerization, and enhanced kinase activity were all abolished if either the W290G or the T341P mutation was expressed in conjunction with mutations that eliminate the disulfide bond in the third immunoglobulin-like domain (Ig-3). These results demonstrate a requirement for the Ig-3 cysteine residues in the activation of FGFR2 by noncysteine mutations. Molecular modeling also reveals that noncysteine mutations may activate FGFR2 by altering the conformation of the Ig-3 domain near the disulfide bond, preventing the formation of an intramolecular bond. This allows the unbonded cysteine residues to participate in intermolecular disulfide bonding, resulting in constitutive activation of the receptor.

Protein kinase mutants of human ATR increase sensitivity to UV and ionizing radiation and abrogate cell cycle checkpoint control

In fission yeast, the rad3 gene product plays a critical role in sensing DNA structure defects and activating damage response pathways. A structural homologue of rad3 in humans (ATR) has been identified based on sequence similarity in the protein kinase domain. General information regarding ATR expression, protein kinase activity, and cellular localization is known, but its function in human cells remains undetermined. In the current study, the ATR protein was examined by gel filtration of protein extracts and was found to exist predominantly as part of a large protein complex. A kinase-inactivated form of the ATR gene was prepared by site-directed mutagenesis and was used in transfection experiments to probe the function of this complex. Introduction of this kinase-dead ATR into a normal fibroblast cell line, an ATM-deficient fibroblast line derived from a patient with ataxia–telangiectasia, or a p53 mutant cell line all resulted in significant losses in cell viability. Clones expressing the kinase-dead ATR displayed increased sensitivity to x-rays and UV and a loss of checkpoint control. We conclude that ATR functions as a critical part of a protein complex that mediates responses to ionizing and UV radiation in human cells. These responses include effects on cell viability and cell cycle checkpoint control.

Activation of protein kinase C by tyrosine phosphorylation in response to H2O2

Protein kinase C (PKC) isoforms, α, βI, and γ of cPKC subgroup, δ and ɛ of nPKC subgroup, and ζ of aPKC subgroup, were tyrosine phosphorylated in COS-7 cells in response to H2O2. These isoforms isolated from the H2O2-treated cells showed enhanced enzyme activity to various extents. The enzymes, PKC α and δ, recovered from the cells were independent of lipid cofactors for their catalytic activity. Analysis of mutated molecules of PKC δ showed that tyrosine residues, which are conserved in the catalytic domain of the PKC family, are critical for PKC activation induced by H2O2. These results suggest that PKC isoforms can be activated through tyrosine phosphorylation in a manner unrelated to receptor-coupled hydrolysis of inositol phospholipids.

Three members of a novel small gene-family from Arabidopsis thaliana able to complement functionally an Escherichia coli mutant defective in PAPS reductase activity encode proteins with a thioredoxin-like domain and “APS reductase” activity

Three different cDNAs, Prh-19, Prh-26, and Prh-43 [3′-phosphoadenosine-5′-phosphosulfate (PAPS) reductase homolog], have been isolated by complementation of an Escherichia coli cysH mutant, defective in PAPS reductase activity, to prototrophy with an Arabidopsis thaliana cDNA library in the expression vector λYES. Sequence analysis of the cDNAs revealed continuous open reading frames encoding polypeptides of 465, 458, and 453 amino acids, with calculated molecular masses of 51.3, 50.5, and 50.4 kDa, respectively, that have strong homology with fungal, yeast, and bacterial PAPS reductases. However, unlike microbial PAPS reductases, each PRH protein has an N-terminal extension, characteristic of a plastid transit peptide, and a C-terminal extension that has amino acid and deduced three-dimensional homology to thioredoxin proteins. Adenosine 5′-phosphosulfate (APS) was shown to be a much more efficient substrate than PAPS when the activity of the PRH proteins was tested by their ability to convert 35S-labeled substrate to acid-volatile 35S-sulfite. We speculate that the thioredoxin-like domain is involved in catalytic function, and that the PRH proteins may function as novel “APS reductase” enzymes. Southern hybridization analysis showed the presence of a small multigene family in the Arabidopsis genome. RNA blot hybridization with gene-specific probes revealed for each gene the presence of a transcript of ≈1.85 kb in leaves, stems, and roots that increased on sulfate starvation. To our knowledge, this is the first report of the cloning and characterization of plant genes that encode proteins with APS reductase activity and supports the suggestion that APS can be utilized directly, without activation to PAPS, as an intermediary substrate in reductive sulfate assimilation.