Chromatin structure mapping in Saccharomyces cerevisiae in vivo with DNase I

Most methods for assessment of chromatin structure involve chemical or nuclease damage to DNA followed by analysis of distribution and susceptibility of cutting sites. The agents used generally do not permeate cells, making nuclear isolation mandatory. In vivo mapping strategies might allow detection of labile constituents and/or structures that are lost when chromatin is swollen in isolated nuclei at low ionic strengths. DNase I has been the most widely used enzyme to detect chromatin sites where DNA is active in transcription, replication or recombination. We have introduced the bovine DNase I gene into yeast under control of a galactose-responsive promoter. Expression of the nuclease leads to DNA degradation and cell death. Shorter exposure to the active enzyme allows mapping of chromatin structure in whole cells without isolation of nuclei. The validity and efficacy of the strategy are demonstrated by footprinting a labile repressor bound to its operator. Investigation of the inter-nucleosome linker regions in several types of repressed domains has revealed different degrees of protection in cells, relative to isolated nuclei.

Global properties of the mapping between local amino acid sequence and local structure in proteins.

Local protein structure prediction efforts have consistently failed to exceed approximately 70% accuracy. We characterize the degeneracy of the mapping from local sequence to local structure responsible for this failure by investigating the extent to which similar sequence segments found in different proteins adopt similar three-dimensional structures. Sequence segments 3-15 residues in length from 154 different protein families are partitioned into neighborhoods containing segments with similar sequences using cluster analysis. The consistency of the sequence-to-structure mapping is assessed by comparing the local structures adopted by sequence segments in the same neighborhood in proteins of known structure. In the 154 families, 45% and 28% of the positions occur in neighborhoods in which one and two local structures predominate, respectively. The sequence patterns that characterize the neighborhoods in the first class probably include virtually all of the short sequence motifs in proteins that consistently occur in a particular local structure. These patterns, many of which occur in transitions between secondary structural elements, are an interesting combination of previously studied and novel motifs. The identification of sequence patterns that consistently occur in one or a small number of local structures in proteins should contribute to the prediction of protein structure from sequence.

Cosmic microwave background theory

A long-standing goal of theorists has been to constrain cosmological parameters that define the structure formation theory from cosmic microwave background (CMB) anisotropy experiments and large-scale structure (LSS) observations. The status and future promise of this enterprise is described. Current band-powers in ℓ-space are consistent with a ΔT flat in frequency and broadly follow inflation-based expectations. That the levels are ∼(10−5)2 provides strong support for the gravitational instability theory, while the Far Infrared Absolute Spectrophotometer (FIRAS) constraints on energy injection rule out cosmic explosions as a dominant source of LSS. Band-powers at ℓ ≳ 100 suggest that the universe could not have re-ionized too early. To get the LSS of Cosmic Background Explorer (COBE)-normalized fluctuations right provides encouraging support that the initial fluctuation spectrum was not far off the scale invariant form that inflation models prefer: e.g., for tilted Λ cold dark matter sequences of fixed 13-Gyr age (with the Hubble constant H0 marginalized), ns = 1.17 ± 0.3 for Differential Microwave Radiometer (DMR) only; 1.15 ± 0.08 for DMR plus the SK95 experiment; 1.00 ± 0.04 for DMR plus all smaller angle experiments; 1.00 ± 0.05 when LSS constraints are included as well. The CMB alone currently gives weak constraints on Λ and moderate constraints on Ωtot, but theoretical forecasts of future long duration balloon and satellite experiments are shown which predict percent-level accuracy among a large fraction of the 10+ parameters characterizing the cosmic structure formation theory, at least if it is an inflation variant.

Three-dimensional modeling of and ligand docking to vitamin D receptor ligand binding domain

The ligand binding domain of the human vitamin D receptor (VDR) was modeled based on the crystal structure of the retinoic acid receptor. The ligand binding pocket of our VDR model is spacious at the helix 11 site and confined at the β-turn site. The ligand 1α,25-dihydroxyvitamin D3 was assumed to be anchored in the ligand binding pocket with its side chain heading to helix 11 (site 2) and the A-ring toward the β-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1α- and 25-hydroxyl groups of 1α,25-dihydroxyvitamin D3 were identified and confirmed by mutational analysis: the 1α-hydroxyl group is forming pincer-type hydrogen bonds with S237 and R274 and the 25-hydroxyl group is interacting with H397. Docking potential for various ligands to the VDR model was examined, and the results are in good agreement with our previous three-dimensional structure-function theory.

A unified theory of carcinogenesis based on order–disorder transitions in DNA structure as studied in the human ovary and breast

Fourier transform-infrared/statistics models demonstrate that the malignant transformation of morphologically normal human ovarian and breast tissues involves the creation of a high degree of structural modification (disorder) in DNA, before restoration of order in distant metastases. Order–disorder transitions were revealed by methods including principal components analysis of infrared spectra in which DNA samples were represented by points in two-dimensional space. Differences between the geometric sizes of clusters of points and between their locations revealed the magnitude of the order–disorder transitions. Infrared spectra provided evidence for the types of structural changes involved. Normal ovarian DNAs formed a tight cluster comparable to that of normal human blood leukocytes. The DNAs of ovarian primary carcinomas, including those that had given rise to metastases, had a high degree of disorder, whereas the DNAs of distant metastases from ovarian carcinomas were relatively ordered. However, the spectra of the metastases were more diverse than those of normal ovarian DNAs in regions assigned to base vibrations, implying increased genetic changes. DNAs of normal female breasts were substantially disordered (e.g., compared with the human blood leukocytes) as were those of the primary carcinomas, whether or not they had metastasized. The DNAs of distant breast cancer metastases were relatively ordered. These findings evoke a unified theory of carcinogenesis in which the creation of disorder in the DNA structure is an obligatory process followed by the selection of ordered, mutated DNA forms that ultimately give rise to metastases.

Vascular endothelial growth factor: Crystal structure and functional mapping of the kinase domain receptor binding site

Vascular endothelial growth factor (VEGF) is a homodimeric member of the cystine knot family of growth factors, with limited sequence homology to platelet-derived growth factor (PDGF) and transforming growth factor β2 (TGF-β). We have determined its crystal structure at a resolution of 2.5 Å, and identified its kinase domain receptor (KDR) binding site using mutational analysis. Overall, the VEGF monomer resembles that of PDGF, but its N-terminal segment is helical rather than extended. The dimerization mode of VEGF is similar to that of PDGF and very different from that of TGF-β. Mutational analysis of VEGF reveals that symmetrical binding sites for KDR are located at each pole of the VEGF homodimer. Each site contains two functional “hot spots” composed of binding determinants presented across the subunit interface. The two most important determinants are located within the largest hot spot on a short, three-stranded sheet that is conserved in PDGF and TGF-β. Functional analysis of the binding epitopes for two receptor-blocking antibodies reveal different binding determinants near each of the KDR binding hot spots.

Structure of the gene for porcine peptide antibiotic PR-39, a cathelin gene family member: comparative mapping of the locus for the human peptide antibiotic FALL-39.

PR-39 is a porcine 39-aa peptide antibiotic composed of 49% proline and 24% arginine, with an activity against Gram-negative bacteria comparable to that of tetracycline. In Escherichia coli, it inhibits DNA and protein synthesis. PR-39 was originally isolated from pig small intestine, but subsequent cDNA cloning showed that the gene is expressed in the bone marrow. The open reading frame of the clone showed that PR-39 is made as 173-aa precursor whose proregion belongs to the cathelin family. The PR39 gene, which is rather compact and spans only 1784 bp has now been sequenced. The coding information is split into four exons. The first exon contains the signal sequence of 29 residues and the first 37 residues of the cathelin propart. Exons 2 and 3 contain only cathelin information, while exon 4 codes for the four C-terminal cathelin residues and the mature PR-39 peptide extended by three residues. The sequenced upstream region (1183 bp) contains four potential recognition sites for NF-IL6 and three for APRF, transcription factors known to regulate genes for both cytokines and acute phase response factors. Genomic hybridizations revealed a fairly high level of restriction fragment length polymorphism and indicated that there are at least two copies of the PR39 gene in the pig genome. PR39 was mapped to pig chromosome 13 by linkage and in situ hybridization mapping. The gene for the human peptide antibiotic FALL-39 (also a member of the cathelin family) was mapped to human chromosome 3, which is homologous to pig chromosome 13.

Dynamic mapping at the laminar level of odor-elicited responses in rat olfactory bulb by functional MRI

We have applied functional MRI (fMRI) based on blood oxygenation level-dependent (BOLD) image-contrast to map odor-elicited olfactory responses at the laminar level in the rat olfactory bulb (OB) elicited by iso-amyl acetate (10−2 dilution of saturated vapor) with spatial and temporal resolutions of 220×220×1,000 μm and 36 s. The laminar structure of the OB was clearly depicted by high-resolution in vivo anatomical MRI with spatial resolution of 110×110×1,000 μm. In repeated BOLD fMRI measurements, highly significant (P < 0.001) foci were located in the outer layers of both OBs. The occurrence of focal OB activity within a domain at the level of individual glomeruli or groups of glomeruli was corroborated on an intra- and inter-animal basis under anesthetized conditions with this noninvasive method. The dynamic studies demonstrated that the odor-elicited BOLD activations were highly reproducible on a time scale of minutes, whereas over tens of minutes the activations sometimes varied slowly. We found large BOLD signal (ΔS/S = 10–30%) arising from the olfactory nerve layer, which is devoid of synapses and composed of unmyelinated fibers and glial cells. Our results support previous studies with other methods showing that odors elicit activity within glomerular layer domains in the mammalian OB, and extend the analysis to shorter time periods at the level of individual glomeruli or groups of glomeruli. With further improvement, BOLD fMRI should be ideal for systematic analysis of the functional significance of individual glomeruli in olfactory information encoding and of spatiotemporal processing within the olfactory system.

Linkage disequilibrium mapping in isolated populations: The example of Finland revisited

Linkage disequilibrium analysis can provide high resolution in the mapping of disease genes because it incorporates information on recombinations that have occurred during the entire period from the mutational event to the present. A circumstance particularly favorable for high-resolution mapping is when a single founding mutation segregates in an isolated population. We review here the population structure of Finland in which a small founder population some 100 generations ago has expanded into 5.1 million people today. Among the 30-odd autosomal recessive disorders that are more prevalent in Finland than elsewhere, several appear to have segregated for this entire period in the “panmictic” southern Finnish population. Linkage disequilibrium analysis has allowed precise mapping and determination of genetic distances at the 0.1-cM level in several of these disorders. Estimates of genetic distance have proven accurate, but previous calculations of the confidence intervals were too small because sampling variation was ignored. In the north and east of Finland the population can be viewed as having been “founded” only after 1500. Disease mutations that have undergone such a founding bottleneck only 20 or so generations ago exhibit linkage disequilibrium and haplotype sharing over long genetic distances (5–15 cM). These features have been successfully exploited in the mapping and cloning of many genes. We review the statistical issues of fine mapping by linkage disequilibrium and suggest that improved methodologies may be necessary to map diseases of complex etiology that may have arisen from multiple founding mutations.

A physical basis for protein secondary structure

A physical theory of protein secondary structure is proposed and tested by performing exceedingly simple Monte Carlo simulations. In essence, secondary structure propensities are predominantly a consequence of two competing local effects, one favoring hydrogen bond formation in helices and turns, the other opposing the attendant reduction in sidechain conformational entropy on helix and turn formation. These sequence specific biases are densely dispersed throughout the unfolded polypeptide chain, where they serve to preorganize the folding process and largely, but imperfectly, anticipate the native secondary structure.

Localization of the N terminus of hepatitis B virus capsid protein by peptide-based difference mapping from cryoelectron microscopy

Recently, cryoelectron microscopy of isolated macromolecular complexes has advanced to resolutions below 10 Å, enabling direct visualization of α-helical secondary structure. To help correlate such density maps with the amino acid sequences of the component proteins, we advocate peptide-based difference mapping, i.e., insertion of peptides, ≈10 residues long, at targeted points in the sequence and visualization of these peptides as bulk labels in cryoelectron microscopy-derived difference maps. As proof of principle, we have appended an extraneous octapeptide at the N terminus of hepatitis B virus capsid protein and determined its location on the capsid surface by difference imaging at 11 Å resolution. Hepatitis B virus capsids are icosahedral particles, ≈300 Å in diameter, made up of T-shaped dimers (subunit Mr, 16–21 kDa, depending on construct). The stems of the Ts protrude outward as spikes, whereas the crosspieces pack to form the contiguous shell. The two N termini per dimer reside on either side of the spike-stem, at the level at which it enters the shell. This location is consistent with formation of the known intramolecular disulfide bond between the cysteines at positions 61 and −7 (in the residual propeptide) in the “e-antigen” form of the capsid protein and has implications for why this clinically important antigen remains unassembled in vivo.

Propagating structure of Alzheimer’s β-amyloid(10–35) is parallel β-sheet with residues in exact register

The pathognomonic plaques of Alzheimer’s disease are composed primarily of the 39- to 43-aa β-amyloid (Aβ) peptide. Crosslinking of Aβ peptides by tissue transglutaminase (tTg) indicates that Gln15 of one peptide is proximate to Lys16 of another in aggregated Aβ. Here we report how the fibril structure is resolved by mapping interstrand distances in this core region of the Aβ peptide chain with solid-state NMR. Isotopic substitution provides the source points for measuring distances in aggregated Aβ. Peptides containing a single carbonyl 13C label at Gln15, Lys16, Leu17, or Val18 were synthesized and evaluated by NMR dipolar recoupling methods for the measurement of interpeptide distances to a resolution of 0.2 Å. Analysis of these data establish that this central core of Aβ consists of a parallel β-sheet structure in which identical residues on adjacent chains are aligned directly, i.e., in register. Our data, in conjunction with existing structural data, establish that the Aβ fibril is a hydrogen-bonded, parallel β-sheet defining the long axis of the Aβ fibril propagation.

Mapping the interface between calmodulin and MARCKS-related protein by fluorescence spectroscopy

MARCKS-related protein (MRP) is a myristoylated protein kinase C substrate that binds calmodulin (CaM) with nanomolar affinity. To obtain structural information on this protein, we have engineered 10 tryptophan residues between positions 89 and 104 in the effector domain, a 24-residue-long amphipathic segment that mediates binding of MRP to CaM. We show that the effector domain is in a polar environment in free MRP, suggesting exposure to water, in agreement with a rod-shaped structure of the protein. The effector domain participates in the binding of MRP to CaM, as judged by the dramatic changes observed in the fluorescent properties of the mutants on complex formation. Intermolecular quenching of the fluorescence emission of the tryptophan residues in MRP by selenomethionine residues engineered in CaM reveals that the N-terminal side of the effector domain contacts the C-terminal domain of CaM, whereas the C-terminal side of the effector domain contacts the N-terminal domain of CaM. Finally, a comparison of the fluorescent properties of the myristoylated and unmyristoylated forms of a construct in which a tryptophan residue was introduced at position 4 close to the myristoylated N terminus of MRP suggests that the lipid moiety is also involved in the interaction of MRP with CaM.

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.