Transcriptional activation by protein-induced DNA bending: evidence for a DNA structural transmission model.

Integration host factor (IHF) is a DNA-bending protein that binds to an upstream activating sequence (UAS1) and, on a negatively supercoiled DNA template, activates transcription from the ilvPG promoter of the ilvG-MEDA operon of Escherichia coli. The transcriptional initiation site of the ilvGMEDA operon is located 92 bp downstream of UAS1. Activation is still observed when the orientation of the upstream IHF binding site is reversed. This manipulation places the IHF binding site on the opposite face of the DNA helix, directs the IHF-induced DNA bend in the opposite direction, and presents the opposite face of the nonsymmetrical, heterodimeric, IHF molecule to the downstream RNA polymerase. Lymphoid enhancer-binding factor, LEF-1, is a DNA-bending, lymphoid-specific, mammalian transcription factor that shares no amino acid sequence similarity with IHF. When the IHF site in UAS1 is replaced with a LEF-1 site, LEF-1 activates transcription from the downstream ilvPG promoter in E. coli as well as it is activated by its natural activator, IHF. These results suggest that specific interactions between IHF and RNA polymerase are not required for activation. The results of DNA structural studies show that IHF forms a protein-DNA complex in the UAS1 region that, in the absence of RNA polymerase, alters the structure of the DNA helix in the -10 hexanucleotide region of the downstream ilvPG promoter. The results of in vitro abortive transcription assays show that IIIF also increases the apparent rate of RNA polymerase isomerization from a closed to an open complex. We suggest, therefore, that IHF activates transcription by forming a higher-order protein-DNA complex in the UAS1 region that structurally alters the DNA helix in a way that facilitates open complex formation at the downstream ilvPG promoter site.

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 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.

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.

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.

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.

A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

Quantitative fine-structural analysis of olfactory cortical synapses

To determine the extent to which hippocampal synapses are typical of those found in other cortical regions, we have carried out a quantitative analysis of olfactory cortical excitatory synapses, reconstructed from serial electron micrograph sections of mouse brain, and have compared these new observations with previously obtained data from hippocampus. Both superficial and deep layer I olfactory cortical synapses were studied. Although individual synapses in each of the areas—CA1 hippocampus, olfactory cortical layer Ia, olfactory cortical area Ib—might plausibly have been found in any of the other areas, the average characteristics of the three synapse populations are distinct. Olfactory cortical synapses in both layers are, on average, about 2.5 times larger than their hippocampal counterparts. The layer Ia olfactory cortical synapses have fewer synaptic vesicles than do the layer Ib synapses, but the absolute number of vesicles docked to the active zone in the layer Ia olfactory cortical synapses is about equal to the docked vesicle number in the smaller hippocampal synapses. As would be predicted from studies on hippocampus that relate paired-pulse facilitation to the number of docked vesicles, the synapses in layer 1a exhibit facilitation, whereas the ones in layer 1b do not. Although hippocampal synapses provide as a good model system for central synapses in general, we conclude that significant differences in the average structure of synapses from one cortical region to another exist, and this means that generalizations based on a single synapse type must be made with caution.

A model for the structure of the P domains in the subtilisin-like prohormone convertases

The proprotein convertases are a family of at least seven calcium-dependent endoproteases that process a wide variety of precursor proteins in the secretory pathway. All members of this family possess an N-terminal proregion, a subtilisin-like catalytic module, and an additional downstream well-conserved region of ≈150 amino acid residues, the P domain, which is not found in any other subtilase. The pro and catalytic domains cannot be expressed in the absence of the P domains; their thermodynamic instability may be attributable to the presence of large numbers of negatively charged Glu and Asp side chains in the substrate binding region for recognition of multibasic residue cleavage sites. Based on secondary structure predictions, we here propose that the P domains consist of 8-stranded β-barrels with well-organized inner hydrophobic cores, and therefore are independently folded components of the proprotein convertases. We hypothesize further that the P domains are integrated through strong hydrophobic interactions with the catalytic domains, conferring structural stability and regulating the properties and activity of the convertases. A molecular model of these interdomain interactions is proposed in this report.

Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Chrysanthemum chlorotic mottle viroid: Unusual structural properties of a subgroup of self-cleaving viroids with hammerhead ribozymes

The causal agent of chrysanthemum chlorotic mottle (CChM) disease has been identified, cloned, and sequenced. It is a viroid RNA (CChMVd) of 398–399 nucleotides. In vitro transcripts with the complete CChMVd sequence were infectious and induced the typical symptoms of the CChM disease. CChMVd can form hammerhead structures in both polarity strands. Plus and minus monomeric CChMVd RNAs self-cleaved during in vitro transcription and after purification as predicted by these structures, which are stable and most probably act as single hammerhead structures as in peach latent mosaic viroid (PLMVd), but not in avocado sunblotch viroid (ASBVd). Moreover, the plus CChMVd hammerhead structure also appears to be active in vivo, because the 5′ terminus of the linear plus CChMVd RNA isolated from infected tissue is that predicted by the corresponding hammerhead ribozyme. Both hammerhead structures of CChMVd display some peculiarities: the plus self-cleaving domain has an unpaired A after the conserved A9 residue, and the minus one has an unusually long helix II. The most stable secondary structure predicted for CChMVd is a branched conformation that does not fulfill the rod-like or quasi-rod-like model proposed for the in vitro structure of most viroids with the exception of PLMVd, whose proposed secondary structure of lowest free energy also is branched. The unusual conformation of CChMVd and PLMVd is supported by their insolubility in 2 M LiCl, in contrast to ASBVd and a series of representative non-self-cleaving viroids that are soluble under the same high salt conditions. These results support the classification of self-cleaving viroids into two subgroups, one formed by ASBVd and the other one by PLMVd and CChMVd.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

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.

Tissue phenotype depends on reciprocal interactions between the extracellular matrix and the structural organization of the nucleus

What determines the nuclear organization within a cell and whether this organization itself can impose cellular function within a tissue remains unknown. To explore the relationship between nuclear organization and tissue architecture and function, we used a model of human mammary epithelial cell acinar morphogenesis. When cultured within a reconstituted basement membrane (rBM), HMT-3522 cells form polarized and growth-arrested tissue-like acini with a central lumen and deposit an endogenous BM. We show that rBM-induced morphogenesis is accompanied by relocalization of the nuclear matrix proteins NuMA, splicing factor SRm160, and cell cycle regulator Rb. These proteins had distinct distribution patterns specific for proliferation, growth arrest, and acini formation, whereas the distribution of the nuclear lamina protein, lamin B, remained unchanged. NuMA relocalized to foci, which coalesced into larger assemblies as morphogenesis progressed. Perturbation of histone acetylation in the acini by trichostatin A treatment altered chromatin structure, disrupted NuMA foci, and induced cell proliferation. Moreover, treatment of transiently permeabilized acini with a NuMA antibody led to the disruption of NuMA foci, alteration of histone acetylation, activation of metalloproteases, and breakdown of the endogenous BM. These results experimentally demonstrate a dynamic interaction between the extracellular matrix, nuclear organization, and tissue phenotype. They further show that rather than passively reflecting changes in gene expression, nuclear organization itself can modulate the cellular and tissue phenotype.

Estimating the probability for a protein to have a new fold: A statistical computational model

Structural genomics aims to solve a large number of protein structures that represent the protein space. Currently an exhaustive solution for all structures seems prohibitively expensive, so the challenge is to define a relatively small set of proteins with new, currently unknown folds. This paper presents a method that assigns each protein with a probability of having an unsolved fold. The method makes extensive use of protomap, a sequence-based classification, and scop, a structure-based classification. According to protomap, the protein space encodes the relationship among proteins as a graph whose vertices correspond to 13,354 clusters of proteins. A representative fold for a cluster with at least one solved protein is determined after superposition of all scop (release 1.37) folds onto protomap clusters. Distances within the protomap graph are computed from each representative fold to the neighboring folds. The distribution of these distances is used to create a statistical model for distances among those folds that are already known and those that have yet to be discovered. The distribution of distances for solved/unsolved proteins is significantly different. This difference makes it possible to use Bayes' rule to derive a statistical estimate that any protein has a yet undetermined fold. Proteins that score the highest probability to represent a new fold constitute the target list for structural determination. Our predicted probabilities for unsolved proteins correlate very well with the proportion of new folds among recently solved structures (new scop 1.39 records) that are disjoint from our original training set.