Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)2 has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by 13C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Plant growth in elevated CO2 alters mitochondrial number and chloroplast fine structure

With increasing interest in the effects of elevated atmospheric CO2 on plant growth and the global carbon balance, there is a need for greater understanding of how plants respond to variations in atmospheric partial pressure of CO2. Our research shows that elevated CO2 produces significant fine structural changes in major cellular organelles that appear to be an important component of the metabolic responses of plants to this global change. Nine species (representing seven plant families) in several experimental facilities with different CO2-dosing technologies were examined. Growth in elevated CO2 increased numbers of mitochondria per unit cell area by 1.3–2.4 times the number in control plants grown in lower CO2 and produced a statistically significant increase in the amount of chloroplast stroma (nonappressed) thylakoid membranes compared with those in lower CO2 treatments. There was no observable change in size of the mitochondria. However, in contrast to the CO2 effect on mitochondrial number, elevated CO2 promoted a decrease in the rate of mass-based dark respiration. These changes may reflect a major shift in plant metabolism and energy balance that may help to explain enhanced plant productivity in response to elevated atmospheric CO2 concentrations.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., blast and fasta validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.

Solution structure of the 2-amino-1- methyl-6-phenylimidazo[4,5-b]pyridine C8-deoxyguanosine adduct in duplex DNA

The carcinogenic heterocyclic amine (HA) 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) is formed during the cooking of various meats. To enable structure/activity studies aimed at understanding how DNA damaged by a member of the HA class of compounds can ultimately lead to cancer, we have determined the first solution structure of an 11-mer duplex containing the C8-dG adduct formed by reaction with N-acetoxy-PhIP. A slow conformational exchange is observed in which the PhIP ligand either intercalates into the DNA helix by denaturing and displacing the modified base pair (main form) or is located outside the helix in a minimally perturbed B-DNA duplex (minor form). In the main base-displaced intercalation structure, the minor groove is widened, and the major groove is compressed at the lesion site because of the location of the bulky PhIP-N-methyl and phenyl ring in the minor groove; this distortion causes significant bending of the helix. The PhIP phenyl ring interacts with the phosphodiester-sugar ring backbone of the complementary strand and its fast rotation with respect to the intercalated imidazopyridine ring causes substantial distortions at this site, such as unwinding and bulging-out of the strand. The glycosidic torsion angle of the [PhIP]dG residue is syn, and the displaced guanine base is directed toward the 3′ end of the modified strand. This study contributes, to our knowledge, the first structural information on the biologically relevant HA class to a growing body of knowledge about how conformational similarities and differences for a variety of types of lesions can influence protein interactions and ultimately biological outcome.

Interplay of structure and disorder in cochaperonin mobile loops.

Protein-protein interactions typically are characterized by highly specific interfaces that mediate binding with precisely tuned affinities. Binding of the Escherichia coli cochaperonin GroES to chaperonin GroEL is mediated, at least in part, by a mobile polypeptide loop in GroES that becomes immobilized in the GroEL/GroES/nucleotide complex. The bacteriophage T4 cochaperonin Gp31 possesses a similar highly flexible polypeptide loop in a region of the protein that shows low, but significant, amino acid similarity with GroES and other cochaperonins. When bound to GroEL, a synthetic peptide representing the mobile loop of either GroES or Gp31 adopts a characteristic bulged hairpin conformation as determined by transferred nuclear Overhauser effects in NMR spectra. Thermodynamic considerations suggest that flexible disorder in the cochaperonin mobile loops moderates their affinity for GroEL to facilitate cycles of chaperonin-mediated protein folding.

Structure, tissue distribution, and chromosomal localization of the prepronociceptin gene.

Nociceptin (orphanin FQ), the newly discovered natural agonist of opioid receptor-like (ORL1) receptor, is a neuropeptide that is endowed with pronociceptive activity in vivo. Nociceptin is derived from a larger precursor, prepronociceptin (PPNOC), whose human, mouse, and rat genes we have now isolated. The PPNOC gene is highly conserved in the three species and displays organizational features that are strikingly similar to those of the genes of preproenkephalin, preprodynorphin, and preproopiomelanocortin, the precursors to endogenous opioid peptides, suggesting the four genes belong to the same family-i.e., have a common evolutionary origin. The PPNOC gene encodes a single copy of nociceptin as well as of other peptides whose sequence is strictly conserved across murine and human species; hence it is likely to be neurophysiologically significant. Northern blot analysis shows that the PPNOC gene is predominantly transcribed in the central nervous system (brain and spinal cord) and, albeit weakly, in the ovary, the sole peripheral organ expressing the gene. By using a radiation hybrid cell line panel, the PPNOC gene was mapped to the short arm of human chromosome 8 (8p21), between sequence-tagged site markers WI-5833 and WI-1172, in close proximity of the locus encoding the neurofilament light chain NEFL. Analysis of yeast artificial chromosome clones belonging to the WC8.4 contig covering the 8p21 region did not allow to detect the presence of the gene on these yeast artificial chromosomes, suggesting a gap in the coverage within this contig.

Identification of pentosidine as a native structure for advanced glycation end products in beta-2-microglobulin-containing amyloid fibrils in patients with dialysis-related amyloidosis.

beta-2-Microglobulin (beta-2m) is a major constituent of amyloid fibrils in patients with dialysis-related amyloidosis (DRA). Recently, we found that the pigmented and fluorescent adducts formed nonenzymatically between sugar and protein, known as advanced glycation end products (AGEs), were present in beta-2m-containing amyloid fibrils, suggesting the possible involvement of AGE-modified beta-2m in bone and joint destruction in DRA. As an extension of our search for the native structure of AGEs in beta-2m of patients with DRA, the present study focused on pentosidine, a fluorescent cross-linked glycoxidation product. Determination by both HPLC assay and competitive ELISA demonstrated a significant amount of pentosidine in amyloid-fibril beta-2m from long-term hemodialysis patients with DRA, and the acidic isoform of beta-2m in the serum and urine of hemodialysis patients. A further immunohistochemical study revealed the positive immunostaining for pentosidine and immunoreactive AGEs and beta-2m in macrophage-infiltrated amyloid deposits of long-term hemodialysis patients with DRA. These findings implicate a potential link of glycoxidation products in long-lived beta-2m-containing amyloid fibrils to the pathogenesis of DRA.

Structure of the N intermediate of bacteriorhodopsin revealed by x-ray diffraction.

X-ray diffraction experiments revealed the structure of the N photointermediate of bacteriorhodopsin. Since the retinal Schiff base is reprotonated from Asp-96 during the M to N transition in the photocycle, and Asp-96 is reprotonated during the lifetime of the N intermediate, or immediately after, N is a key intermediate for understanding the light-driven proton pump. The N intermediate accumulates in large amounts during continuous illumination of the F171C mutant at pH 7 and 5 degrees Celsius. Small but significant changes of the structure were detected in the x-ray diffraction profile under these conditions. The changes were reversible and reproducible. The difference Fourier map indicates that the major change occurs near helix F. The observed diffraction changes between N and the original state were essentially identical to the diffraction changes reported for the M intermediate of the D96N mutant of bacteriorhodopsin. Thus, we find that the protein conformations of the M and N intermediates of the photocycle are essentially the same, in spite of the fact that in M the Schiff base is unprotonated and in N it is protonated. The observed structural change near helix F will increase access of the Schiff base and Asp-96 to the cytoplasmic surface and facilitate the proton transfer events that begin with the decay of the M state.

Bathorhodopsin structure in the room-temperature rhodopsin photosequence: picosecond time-resolved coherent anti-Stokes Raman scattering.

Structural changes in the retinal chromophore during the formation of the bathorhodopsin intermediate (bathoRT) in the room-temperature rhodopsin (RhRT) photosequence (i.e., vision) are examined using picosecond time-resolved coherent anti-Stokes Raman scattering. Specifically, the retinal structure assignable to bathoRT following 8-ps excitation of RhRT is measured via vibrational Raman spectroscopy at a 200-ps time delay where the only intermediate present is bathoRT. Significant differences are observed between the C=C stretching frequencies of the retinal chromophore at low temperature where bathorhodopsin is stabilized and at room temperature where bathorhodopsin is a transient species in the RhRT photosequence. These vibrational data are discussed in terms of the formation of bathoRT, an important step in the energy storage/transduction mechanism of RhRT.

Crystal structure of the complex of a catalytic antibody Fab fragment with a transition state analog: structural similarities in esterase-like catalytic antibodies.

The x-ray structure of the complex of a catalytic antibody Fab fragment with a phosphonate transition-state analog has been determined. The antibody (CNJ206) catalyzes the hydrolysis of p-nitrophenyl esters with significant rate enhancement and substrate specificity. Comparison of this structure with that of the uncomplexed Fab fragment suggests hapten-induced conformational changes: the shape of the combining site changes from a shallow groove in the uncomplexed Fab to a deep pocket where the hapten is buried. Three hydrogen-bond donors appear to stabilize the charged phosphonate group of the hapten: two NH groups of the heavy (H) chain complementarity-determining region 3 (H3 CDR) polypeptide chain and the side-chain of histidine-H35 in the H chain (His-H35) in the H1 CDR. The combining site shows striking structural similarities to that of antibody 17E8, which also has esterase activity. Both catalytic antibody ("abzyme") structures suggest that oxyanion stabilization plays a significant role in their rate acceleration. Additional catalytic groups that improve efficiency are not necessarily induced by the eliciting hapten; these groups may occur because of the variability in the combining sites of different monoclonal antibodies that bind to the same hapten.

Structure determination of murine mitochondrial carbonic anhydrase V at 2.45-A resolution: implications for catalytic proton transfer and inhibitor design.

The three-dimensional structure of murine mitochondrial carbonic anhydrase V has been determined and refined at 2.45-A resolution (crystallographic R factor = 0.187). Significant structural differences unique to the active site of carbonic anhydrase V are responsible for differences in the mechanism of catalytic proton transfer as compared with other carbonic anhydrase isozymes. In the prototypical isozyme, carbonic anhydrase II, catalytic proton transfer occurs via the shuttle group His-64; carbonic anhydrase V has Tyr-64, which is not an efficient proton shuttle due in part to the bulky adjacent side chain of Phe-65. Based on analysis of the structure of carbonic anhydrase V, we speculate that Tyr-131 may participate in proton transfer due to its proximity to zinc-bound solvent, its solvent accessibility, and its electrostatic environment in the protein structure. Finally, the design of isozyme-specific inhibitors is discussed in view of the complex between carbonic anhydrase V and acetazolamide, a transition-state analogue. Such inhibitors may be physiologically important in the regulation of blood glucose levels.

Genetic and comparative analyses reveal an alternative secondary structure in the region of nt 912 of Escherichia coli 16S rRNA.

Mutations at position 912 of Escherichia coli 16S rRNA result in two notable phenotypes. The C-->U transition confers resistance to streptomycin, a translational-error-inducing antibiotic, while a C-->G transversion causes marked retardation of cell growth rate. Starting with the slow-growing G912 mutant, random mutagenesis was used to isolate a second site mutation that restored growth nearly to the wild-type rate. The second site mutation was identified as a G-->C transversion at position 885 in 16S rRNA. Cells containing the G912 mutation had an increased doubling time, abnormal sucrose gradient ribosome/subunit profile, increased sensitivity to spectinomycin, dependence upon streptomycin for growth in the presence of spectinomycin, and slower translation rate, whereas cells with the G912/C885 double mutation were similar to wild type in these assays. Comparative analysis showed there was significant covariation between positions 912 and 885. Thus the second-site suppressor analysis, the functional assays, and the comparative data suggest that the interaction between nt 912 and nt 885 is conserved and necessary for normal ribosome function. Furthermore, the comparative data suggest that the interaction extends to include G885-G886-G887 pairing with C912-U911-C910. An alternative secondary structure element for the central domain of 16S rRNA is proposed.

Tertiary structure of an amyloid immunoglobulin light chain protein: a proposed model for amyloid fibril formation.

An immunoglobulin light chain protein was isolated from the urine of an individual (BRE) with systemic amyloidosis. Complete amino acid sequence of the variable region of the light chain (VL) protein established it as a kappa I, which when compared with other kappa I amyloid associated proteins had unique residues, including Ile-34, Leu-40, and Tyr-71. To study the tertiary structure, BRE VL was expressed in Escherichia coli by using a PCR product amplified from the patient BRE's bone marrow DNA. The PCR product was ligated into pCZ11, a thermal-inducible replication vector. Recombinant BRE VL was isolated, purified to homogeneity, and crystallized by using ammonium sulfate as the precipitant. Two crystal forms were obtained. In crystal form I the BRE VL kappa domain crystallizes as a dimer with unit cell constants isomorphous to previously published kappa protein structures. Comparison with a nonamyloid VL kappa domain from patient REI, identified significant differences in position of residues in the hypervariable segments plus variations in framework region (FR) segments 40-46 (FR2) and 66-67 (FR3). In addition, positional differences can be seen along the two types of local diads, corresponding to the monomer-monomer and dimer-dimer interfaces. From the packing diagram, a model for the amyloid light chain (AL) fibril is proposed based on a pseudohexagonal spiral structure with a rise of approximately the width of two dimers per 360 degree turn. This spiral structure could be consistent with the dimensions of amyloid fibrils as determined by electron microscopy.