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.

A local electrostatic change is the cause of the large-scale protein conformation shift in bacteriorhodopsin

During light-driven proton transport bacteriorhodopsin shuttles between two protein conformations. A large-scale structural change similar to that in the photochemical cycle is produced in the D85N mutant upon raising the pH, even without illumination. We report here that (i) the pKa values for the change in crystallographic parameters and for deprotonation of the retinal Schiff base are the same, (ii) the retinal isomeric configuration is nearly unaffected by the protein conformation, and (iii) preventing rotation of the C13—C14 double bond by replacing the retinal with an all-trans locked analogue makes little difference to the Schiff base pKa. We conclude that the direct cause of the conformational shift is destabilization of the structure upon loss of interaction of the positively charged Schiff base with anionic residues that form its counter-ion.

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.

Global air-sea flux of CO2: An estimate based on measurements of sea–air pCO2 difference

Approximately 250,000 measurements made for the pCO2 difference between surface water and the marine atmosphere, ΔpCO2, have been assembled for the global oceans. Observations made in the equatorial Pacific during El Nino events have been excluded from the data set. These observations are mapped on the global 4° × 5° grid for a single virtual calendar year (chosen arbitrarily to be 1990) representing a non-El Nino year. Monthly global distributions of ΔpCO2 have been constructed using an interpolation method based on a lateral advection–diffusion transport equation. The net flux of CO2 across the sea surface has been computed using ΔpCO2 distributions and CO2 gas transfer coefficients across sea surface. The annual net uptake flux of CO2 by the global oceans thus estimated ranges from 0.60 to 1.34 Gt-C⋅yr−1 depending on different formulations used for wind speed dependence on the gas transfer coefficient. These estimates are subject to an error of up to 75% resulting from the numerical interpolation method used to estimate the distribution of ΔpCO2 over the global oceans. Temperate and polar oceans of the both hemispheres are the major sinks for atmospheric CO2, whereas the equatorial oceans are the major sources for CO2. The Atlantic Ocean is the most important CO2 sink, providing about 60% of the global ocean uptake, while the Pacific Ocean is neutral because of its equatorial source flux being balanced by the sink flux of the temperate oceans. The Indian and Southern Oceans take up about 20% each.

RIGS (repeat-induced gene silencing) in Arabidopsis is transcriptional and alters chromatin configuration.

We have previously reported repeat-induced gene silencing (RIGS) in Arabidopsis, in which transgene expression may be silenced epigenetically when repeated sequences are present. Among an allelic series of lines comprising a primary transformant and various recombinant progeny carrying different numbers of drug resistance gene copies at the same locus, silencing was found to depend strictly on repeated sequences and to correlate with an absence of steady-state mRNA. We now report characterization, in nuclei isolated from the same transgenic lines, of gene expression by nuclear run-on assay and of chromatin structure by nuclease protection assay. We find that silencing is correlated with absence of run-on transcripts, indicating that expression is silenced at the level of transcription. We find further that silencing is also correlated with increased resistance to both DNase I and micrococcal nuclease, indicating that the silenced state reflects a change in chromatin configuration. We propose that silencing results when a locally paired region of homologous repeated nucleotide sequences is flanked by unpaired heterologous DNA, which leads chromatin to adopt a local configuration that is difficult to transcribe, and possibly akin to heterochromatin.

A rat genetic map constructed by representational difference analysis markers with suitability for large-scale typing.

Representational difference analysis (RDA) was applied to isolate chromosomal markers in the rat. Four series of RDA [restriction enzymes, BamHI and HindIII; subtraction of ACI/N (ACI) amplicon from BUF/Nac (BUF) amplicon and vice versa] yielded 131 polymorphic markers; 125 of these markers were mapped to all chromosomes except for chromosome X. This was done by using a mapping panel of 105 ACI x BUF F2 rats. To complement the relative paucity of chromosomal markers in the rat, genetically directed RDA, which allows isolation of polymorphic markers in the specific chromosomal region, was performed. By changing the F2 driver-DNA allele frequency around the region, four markers were isolated from the D1Ncc1 locus. Twenty-five of 27 RDA markers were informative regarding the dot blot analysis of amplicons, hybridizing only with tester amplicons. Dot blot analysis at a high density per unit of area made it possible to process a large number of samples. Quantitative trait loci can now be mapped in the rat genome by processing a large number of samples with RDA markers and then by isolating markers close to the loci of interest by genetically directed RDA.

Isolation and characterization of a pseudoautosomal region-specific genetic marker in C57BL/6 mice using genomic representational difference analysis.

Representational difference analysis was used to identify strain-specific differences in the pseudoautosomal region (PAR) of mouse X and Y chromosomes. One second generation (C57BL/6 x Mus spretus) x Mus spretus interspecific backcross male carrying the C57BL/6 (B6) PAR was used for tester DNA. DNA from five backcross males from the same generation that were M. spretus-type for the PAR was pooled for the driver. A cloned probe designated B6-38 was recovered that is B6-specific in Southern analysis. Analysis of genomic DNA from several inbred strains of laboratory mice and diverse Mus species and subspecies identified a characteristic Pst I pattern of fragment sizes that is present only in the C57BL family of strains. Hybridization was observed with sequences in DBA/2J and to a limited extent with Mus musculus (PWK strain) and Mus castaneus DNA. No hybridization was observed in DNA of different Mus species, M. spretus, M. hortulanus, and M. caroli. Genetic analyses of B6-38 was conducted using C57BL congenic males that carry M. spretus alleles for distal X chromosome loci and the PAR and outcrosses of heterozygous congenic females with M. spretus. These analyses demonstrated that the B6-38 sequences were inherited with both the X and Y chromosome. B6-38 sequences were genetically mapped as a locus within the PAR using two interspecific backcrosses. The locus defined by B6-38 is designated DXYRp1. Preliminary analyses of recombination between the distal X chromosome gene amelogenin (Amg) and the PAR loci for either TelXY or sex chromosome association (Sxa) suggest that the locus DXYRp1 maps to the distal portion of the PAR.

A maximum likelihood algorithm for genome mapping of cytogenetic loci from meiotic configuration data.

Frequencies of meiotic configurations in cytogenetic stocks are dependent on chiasma frequencies in segments defined by centromeres, breakpoints, and telomeres. The expectation maximization algorithm is proposed as a general method to perform maximum likelihood estimations of the chiasma frequencies in the intervals between such locations. The estimates can be translated via mapping functions into genetic maps of cytogenetic landmarks. One set of observational data was analyzed to exemplify application of these methods, results of which were largely concordant with other comparable data. The method was also tested by Monte Carlo simulation of frequencies of meiotic configurations from a monotelodisomic translocation heterozygote, assuming six different sample sizes. The estimate averages were always close to the values given initially to the parameters. The maximum likelihood estimation procedures can be extended readily to other kinds of cytogenetic stocks and allow the pooling of diverse cytogenetic data to collectively estimate lengths of segments, arms, and chromosomes.

Identification by representational difference analysis of a homozygous deletion in pancreatic carcinoma that lies within the BRCA2 region.

Homozygous deletions have been central to the discovery of several tumor-suppressor genes, but their finding has often been either serendipitous or the result of a directed search. A recently described technique [Lisitsyn, N., Lisitsyn, N. & Wigler, M. (1993) Science 259, 946-951] held out the potential to efficiently discover such events in an unbiased manner. Here we present the application of the representational difference analysis (RDA) to the study of cancer. We cloned two DNA fragments that identified a homozygous deletion in a human pancreatic adenocarcinoma, mapping to a 1-centimorgan region at chromosome 13q12.3 flanked by the markers D13S171 and D13S260. Interestingly, this lies within the 6-centimorgan region recently identified as the BRCA2 locus of heritable breast cancer susceptibility. This suggests that the same gene may be involved in multiple tumor types and that its function is that of a tumor suppressor rather than that of a dominant oncogene.