Structural model of cytochrome b559 in photosystem II based on a mutant with genetically fused subunits

Photosystem II is a reaction center protein complex located in photosynthetic membranes of plants, algae, and cyanobacteria. Using light energy, photosystem II catalyzes the oxidation of water and the reduction of plastoquinone, resulting in the release of molecular oxygen. A key component of photosystem II is cytochrome b559, a membrane-embedded heme protein with an unknown function. The cytochrome is unusual in that a heme links two separate polypeptide subunits, α and β, either as a heterodimer (αβ) or as two homodimers (α2 and β2). To determine the structural organization of cytochrome b559 in the membrane, we used site-directed mutagenesis to fuse the coding regions of the two respective genes in the cyanobacterium Synechocystis sp. PCC 6803. In this construction, the C terminus of the α subunit (9 kDa) is attached to the N terminus of the β subunit (5 kDa) to form a 14-kDa αβ fusion protein that is predicted to have two membrane-spanning α-helices with antiparallel orientations. Cells containing the αβ fusion protein grow photoautotrophically and assemble functional photosystem II complexes. Optical spectroscopy shows that the αβ fusion protein binds heme and is incorporated into photosystem II. These data support a structural model of cytochrome b559 in which one heme is coordinated to an α2 homodimer and a second heme is coordinated to a β2 homodimer. In this model, each photosystem II complex contains two cytochrome b559 hemes, with the α2 heme located near the stromal side of the membrane and the β2 heme located near the lumenal side.

A structural model for a homeotic protein-extradenticle-DNA complex accounts for the choice of HOX protein in the heterodimer.

The genes of the homeotic complex (HOX) encode DNA binding homeodomain proteins that control developmental fates by differentially regulating the transcription of downstream target genes. Despite their unique in vivo functions, disparate HOX proteins often bind to very similar DNA sequences in vitro. Thus, a critical question is how HOX proteins select the correct sets of target genes in vivo. The homeodomain proteins encoded by the Drosophila extradenticle gene and its mammalian homologues, the pbx genes, contribute to HOX specificity by cooperatively binding to DNA with HOX proteins. For example, the HOX protein labial cooperatively binds with extradenticle protein to a 20-bp oligonucleotide that is sufficient to direct a labial-like expression pattern in Drosophila embryos. Here we have analyzed the protein-DNA interactions that are important for forming the labial-extradenticle-DNA complex. The data suggest a model in which labial and extradenticle, separated by only 4 bp, bind this DNA as a heterodimer in a head-to-tail orientation. We have confirmed several aspects of this model by characterizing extradenticle-HOX binding to mutant oligonucleotides. Most importantly, mutations in base pairs predicted to contact the HOX N-terminal arm resulted in a change in HOX preference in the heterodimer, from labial to Ultrabithorax. These results demonstrate that extradenticle prefers to bind cooperatively with different HOX proteins depending on subtle differences in the heterodimer binding site.

Correlation of the structural and functional domains in the membrane protein Vpu from HIV-1

Vpu is an 81-residue membrane protein encoded by the HIV-1 genome. NMR experiments show that the protein folds into two distinct domains, a transmembrane hydrophobic helix and a cytoplasmic domain with two in-plane amphipathic α-helices separated by a linker region. Resonances in one-dimensional solid-state NMR spectra of uniformly 15N labeled Vpu are clearly segregated into two bands at chemical shift frequencies associated with NH bonds in a transmembrane α-helix, perpendicular to the membrane surface, and with NH bonds in the cytoplasmic helices parallel to the membrane surface. Solid-state NMR spectra of truncated Vpu2–51 (residues 2–51), which contains the transmembrane α-helix and the first amphipathic helix of the cytoplasmic domain, and of a construct Vpu28–81 (residues 28–81), which contains only the cytoplasmic domain, support this structural model of Vpu in the membrane. Full-length Vpu (residues 2–81) forms discrete ion-conducting channels of heterogeneous conductance in lipid bilayers. The most frequent conductances were 22 ± 3 pS and 12 ± 3 pS in 0.5 M KCl and 29 ± 3 pS and 12 ± 3 pS in 0.5 M NaCl. In agreement with the structural model, truncated Vpu2–51, which has the transmembrane helix, forms discrete channels in lipid bilayers, whereas the cytoplasmic domain Vpu28–81, which lacks the transmembrane helix, does not. This finding shows that the channel activity is associated with the transmembrane helical domain. The pattern of channel activity is characteristic of the self-assembly of conductive oligomers in the membrane and is compatible with the structural and functional findings.

Photoactive yellow protein: A structural prototype for the three-dimensional fold of the PAS domain superfamily

PAS domains are found in diverse proteins throughout all three kingdoms of life, where they apparently function in sensing and signal transduction. Although a wealth of useful sequence and functional information has become recently available, these data have not been integrated into a three-dimensional (3D) framework. The very early evolutionary development and diverse functions of PAS domains have made sequence analysis and modeling of this protein superfamily challenging. Limited sequence similarities between the ∼50-residue PAS repeats and one region of the bacterial blue-light photosensor photoactive yellow protein (PYP), for which ground-state and light-activated crystallographic structures have been determined to high resolution, originally were identified in sequence searches using consensus sequence probes from PAS-containing proteins. Here, we found that by changing a few residues particular to PYP function, the modified PYP sequence probe also could select PAS protein sequences. By mapping a typical ∼150-residue PAS domain sequence onto the entire crystallographic structure of PYP, we show that the PAS sequence similarities and differences are consistent with a shared 3D fold (the PAS/PYP module) with obvious potential for a ligand-binding cavity. Thus, PYP appears to prototypically exhibit all the major structural and functional features characteristic of the PAS domain superfamily: the shared PAS/PYP modular domain fold of ∼125–150 residues, a sensor function often linked to ligand or cofactor (chromophore) binding, and signal transduction capability governed by heterodimeric assembly (to the downstream partner of PYP). This 3D PAS/PYP module provides a structural model to guide experimental testing of hypotheses regarding ligand-binding, dimerization, and signal transduction.

α1-Antitrypsin Portland, a bioengineered serpin highly selective for furin: Application as an antipathogenic agent

The important role of furin in the proteolytic activation of many pathogenic molecules has made this endoprotease a target for the development of potent and selective antiproteolytic agents. Here, we demonstrate the utility of the protein-based inhibitor α1-antitrypsin Portland (α1-PDX) as an antipathogenic agent that can be used prophylactically to block furin-dependent cell killing by Pseudomonas exotoxin A. Biochemical analysis of the specificity of a bacterially expressed His- and FLAG-tagged α1-PDX (α1-PDX/hf) revealed the selectivity of the α1-PDX/hf reactive site loop for furin (Ki, 600 pM) but not for other proprotein convertase family members or other unrelated endoproteases. Kinetic studies show that α1-PDX/hf inhibits furin by a slow tight-binding mechanism characteristic of serpin molecules and functions as a suicide substrate inhibitor. Once bound to furin’s active site, α1-PDX/hf partitions with equal probability to undergo proteolysis by furin at the C-terminal side of the reactive center -Arg355-Ile-Pro-Arg358-↓ or to form a kinetically trapped SDS-stable complex with the enzyme. This partitioning between the complex-forming and proteolytic pathways contributes to the ability of α1-PDX/hf to differentially inhibit members of the proprotein convertase family. Finally, we propose a structural model of the α1-PDX-reactive site loop that explains the high degree of enzyme selectivity of this serpin and which can be used to generate small molecule furin inhibitors.

The solution structure of the N-terminal domain of α2-macroglobulin receptor-associated protein

The three-dimensional structure of the N-terminal domain (residues 18–112) of α2-macroglobulin receptor-associated protein (RAP) has been determined by NMR spectroscopy. The structure consists of three helices composed of residues 23–34, 39–65, and 73–88. The three helices are arranged in an up-down-up antiparallel topology. The C-terminal 20 residues were shown not to be in a well defined conformation. A structural model for the binding of RAP to the family of low-density lipoprotein receptors is proposed. It defines a role in binding for both the unordered C terminus and the structural scaffold of the core structure. Pathogenic epitopes for the rat disease Heymann nephritis, an experimental model of human membranous glomerulonephritis, have been identified in RAP and in the large endocytic receptor gp330/megalin. Here we provide the three-dimensional structure of the pathogenic epitope in RAP. The amino acid residues known to form the epitope are in a helix–loop–helix conformation, and from the structure it is possible to rationalize the published results obtained from studies of fragments of the N-terminal domain.

Comparative analyses of the secondary structures of synthetic and intracellular yeast MFA2 mRNAs

The overall folded (global) structure of mRNA may be critical to translation and turnover control mechanisms, but it has received little experimental attention. Presented here is a comparative analysis of the basic features of the global secondary structure of a synthetic mRNA and the same intracellular eukaryotic mRNA by dimethyl sulfate (DMS) structure probing. Synthetic MFA2 mRNA of Saccharomyces cerevisiae first was examined by using both enzymes and chemical reagents to determine single-stranded and hybridized regions; RNAs with and without a poly(A) tail were compared. A folding pattern was obtained with the aid of the mfold program package that identified the model that best satisfied the probing data. A long-range structural interaction involving the 5′ and 3′ untranslated regions and causing a juxtaposition of the ends of the RNA, was examined further by a useful technique involving oligo(dT)-cellulose chromatography and antisense oligonucleotides. DMS chemical probing of A and C nucleotides of intracellular MFA2 mRNA was then done. The modification data support a very similar intracellular structure. When low reactivity of A and C residues is found in the synthetic RNA, ≈70% of the same sites are relatively more resistant to DMS modification in vivo. A slightly higher sensitivity to DMS is found in vivo for some of the A and C nucleotides predicted to be hybridized from the synthetic structural model. With this small mRNA, the translation process and mRNA-binding proteins do not block DMS modifications, and all A and C nucleotides are modified the same or more strongly than with the synthetic RNA.

Two amino acid substitutions convert a guanylyl cyclase, RetGC-1, into an adenylyl cyclase

Guanylyl cyclases (GCs) and adenylyl cyclases (ACs) have fundamental roles in a wide range of cellular processes. Whereas GCs use GTP as a substrate to form cGMP, ACs catalyze the analogous conversion of ATP to cAMP. Previously, a model based on the structure of adenylate cyclase was used to predict the structure of the nucleotide-binding pocket of a membrane guanylyl cyclase, RetGC-1. Based on this model, we replaced specific amino acids in the guanine-binding pocket of GC with their counterparts from AC. A change of two amino acids, E925K together with C995D, is sufficient to completely alter the nucleotide specificity from GTP to ATP. These experiments strongly validate the AC-derived RetGC-1 structural model and functionally confirm the role of these residues in nucleotide discrimination.

An important 2′-OH group for an RNA–protein interaction

We have investigated the role of 2′-OH groups in the specific interaction between the acceptor stem of Escherichia coli tRNACys and cysteine-tRNA synthetase. This interaction provides for the high aminoacylation specificity observed for cysteine-tRNA synthetase. A synthetic RNA microhelix that recapitulates the sequence of the acceptor stem was used as a substrate and variants containing systematic replacement of the 2′-OH by 2′-deoxy or 2′-O-methyl groups were tested. Except for position U73, all substitutions had little effect on aminoacylation. Interestingly, the deoxy substitution at position U73 had no effect on aminoacylation, but the 2′-O-methyl substitution decreased aminoacylation by 10-fold and addition of the even bulkier 2′-O-propyl group decreased aminoacylation by another 2-fold. The lack of an effect by the deoxy substitution suggests that the hydrogen bonding potential of the 2′-OH at position U73 is unimportant for aminoacylation. The decrease in activity upon alkyl substitution suggests that the 2′-OH group instead provides a monitor of the steric environment during the RNA–synthetase interaction. The steric role was confirmed in the context of a reconstituted tRNA and is consistent with the observation that the U73 base is the single most important determinant for aminoacylation and therefore is a site that is likely to be in close contact with cysteine-tRNA synthetase. A steric role is supported by an NMR-based structural model of the acceptor stem, together with biochemical studies of a closely related microhelix. This role suggests that the U73 binding site for cysteine-tRNA synthetase is sterically optimized to accommodate a 2′-OH group in the backbone, but that the hydroxyl group itself is not involved in specific hydrogen bonding interactions.

The Sm domain is an ancient RNA-binding motif with oligo(U) specificity

Sm and Sm-like proteins are members of a family of small proteins that is widespread throughout eukaryotic kingdoms. These proteins form heteromers with one another and bind, as heteromeric complexes, to various RNAs, recognizing primarily short U-rich stretches. Interestingly, completion of several genome projects revealed that archaea also contain genes that may encode Sm-like proteins. Herein, we studied the properties of one Sm-like protein derived from the archaebacterium Archaeoglobus fulgidus and overexpressed in Escherichia coli. This single small protein closely reflects the properties of an Sm or Sm-like protein heteromer. It binds to RNA with a high specificity for oligo(U), and assembles onto the RNA to form a complex that exhibits, as judged by electron microscopy, a ring-like structure similar to the ones observed with the Sm core ribonucleoprotein and the like Sm (LSm) protein heteromer. Importantly, multivariate statistical analysis of negative-stain electron-microscopic images revealed a sevenfold symmetry for the observed ring structure, indicating that the proteins form a homoheptamer. These results support the structural model of the Sm proteins derived from crystallographic studies on Sm heterodimers and demonstrate that the Sm protein family evolved from a single ancestor that was present before the eukaryotic and archaeal kingdoms separated.

Monitoring the GAP catalyzed H-Ras GTPase reaction at atomic resolution in real time

The molecular reaction mechanism of the GTPase-activating protein (GAP)-catalyzed GTP hydrolysis by Ras was investigated by time resolved Fourier transform infrared (FTIR) difference spectroscopy using caged GTP (P3-1-(2-nitro)phenylethyl guanosine 5′-O-triphosphate) as photolabile trigger. This approach provides the complete GTPase reaction pathway with time resolution of milliseconds at the atomic level. Up to now, one structural model of the GAP⋅Ras⋅GDP⋅AlFx transition state analog is known, which represents a “snap shot” along the reaction-pathway. As now revealed, binding of GAP to Ras⋅GTP shifts negative charge from the γ to β phosphate. Such a shift was already identified by FTIR in GTP because of Ras binding and is now shown to be enhanced by GAP binding. Because the charge distribution of the GAP⋅Ras⋅GTP complex thus resembles a more dissociative-like transition state and is more like that in GDP, the activation free energy is reduced. An intermediate is observed on the reaction pathway that appears when the bond between β and γ phosphate is cleaved. In the intermediate, the released Pi is strongly bound to the protein and surprisingly shows bands typical of those seen for phosphorylated enzyme intermediates. All these results provide a mechanistic picture that is different from the intrinsic GTPase reaction of Ras. FTIR analysis reveals the release of Pi from the protein complex as the rate-limiting step for the GAP-catalyzed reaction. The approach presented allows the study not only of single proteins but of protein–protein interactions without intrinsic chromophores, in the non-crystalline state, in real time at the atomic level.

X-ray structure of clotting factor IXa: active site and module structure related to Xase activity and hemophilia B.

Hereditary deficiency of factor IXa (fIXa), a key enzyme in blood coagulation, causes hemophilia B, a severe X chromosome-linked bleeding disorder afflicting 1 in 30,000 males; clinical studies have identified nearly 500 deleterious variants. The x-ray structure of porcine fIXa described here shows the atomic origins of the disease, while the spatial distribution of mutation sites suggests a structural model for factor X activation by phospholipid-bound fIXa and cofactor VIIIa. The 3.0-A-resolution diffraction data clearly show the structures of the serine proteinase module and the two preceding epidermal growth factor (EGF)-like modules; the N-terminal Gla module is partially disordered. The catalytic module, with covalent inhibitor D-Phe-1I-Pro-2I-Arg-3I chloromethyl ketone, most closely resembles fXa but differs significantly at several positions. Particularly noteworthy is the strained conformation of Glu-388, a residue strictly conserved in known fIXa sequences but conserved as Gly among other trypsin-like serine proteinases. Flexibility apparent in electron density together with modeling studies suggests that this may cause incomplete active site formation, even after zymogen, and hence the low catalytic activity of fIXa. The principal axes of the oblong EGF-like domains define an angle of 110 degrees, stabilized by a strictly conserved and fIX-specific interdomain salt bridge. The disorder of the Gla module, whose hydrophobic helix is apparent in electron density, can be attributed to the absence of calcium in the crystals; we have modeled the Gla module in its calcium form by using prothrombin fragment 1. The arched module arrangement agrees with fluorescence energy transfer experiments. Most hemophilic mutation sites of surface fIX residues occur on the concave surface of the bent molecule and suggest a plausible model for the membrane-bound ternary fIXa-FVIIIa-fX complex structure: fIXa and an equivalently arranged fX arch across an underlying fVIIIa subdomain from opposite sides; the stabilizing fVIIIa interactions force the catalytic modules together, completing fIXa active site formation and catalytic enhancement.

A Polymer Model for the Structural Organization of Chromatin Loops and Minibands in Interphase Chromosomes

A quantitative model of interphase chromosome higher-order structure is presented based on the isochore model of the genome and results obtained in the field of copolymer research. G1 chromosomes are approximated in the model as multiblock copolymers of the 30-nm chromatin fiber, which alternately contain two types of 0.5- to 1-Mbp blocks (R and G minibands) differing in GC content and DNA-bound proteins. A G1 chromosome forms a single-chain string of loop clusters (micelles), with each loop ∼1–2 Mbp in size. The number of ∼20 loops per micelle was estimated from the dependence of geometrical versus genomic distances between two points on a G1 chromosome. The greater degree of chromatin extension in R versus G minibands and a difference in the replication time for these minibands (early S phase for R versus late S phase for G) are explained in this model as a result of the location of R minibands at micelle cores and G minibands at loop apices. The estimated number of micelles per nucleus is close to the observed number of replication clusters at the onset of S phase. A relationship between chromosomal and nuclear sizes for several types of higher eukaryotic cells (insects, plants, and mammals) is well described through the micelle structure of interphase chromosomes. For yeast cells, this relationship is described by a linear coil configuration of chromosomes.

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Role of hydrophobic interactions and desolvation in determining the structural properties of a model alpha beta peptide.

Model AB, a 20-amino acid peptide that was designed to adopt an alpha beta tertiary structure stabilized by hydrophobic interactions between residues in adjacent helical and extended segments, exhibited large pKa shifts of several ionizable groups and slow hydrogen/deuterium exchange rates of nearly all the peptide amide groups [Butcher, D. J., Bruch, M. D. & Moe, G. T. (1995) Biopolymers 36, 109-120]. These properties, which depend on structure and hydration, are commonly observed in larger proteins but are quite unusual for small peptides. To identify which of several possible features of the peptide design are most important in determining these properties, several closely related analogs of Model AB were characterized by CD and NMR spectroscopy. The results show that hydrophobic interactions between adjacent helical and extended segments are structure-determining and have the additional effect of altering water-peptide interactions over much of the peptide surface. These results may have important implications for understanding mechanisms of protein folding and for the design of independently folding peptides.