Survey of human mitochondrial diseases using new genomic/proteomic tools.

BACKGROUND: We have constructed Bayesian prior-based, amino-acid sequence profiles for the complete yeast mitochondrial proteome and used them to develop methods for identifying and characterizing the context of protein mutations that give rise to human mitochondrial diseases. (Bayesian priors are conditional probabilities that allow the estimation of the likelihood of an event - such as an amino-acid substitution - on the basis of prior occurrences of similar events.) Because these profiles can assemble sets of taxonomically very diverse homologs, they enable identification of the structurally and/or functionally most critical sites in the proteins on the basis of the degree of sequence conservation. These profiles can also find distant homologs with determined three-dimensional structures that aid in the interpretation of effects of missense mutations. RESULTS: This survey reports such an analysis for 15 missense mutations, one insertion and three deletions involved in Leber's hereditary optic neuropathy, Leigh syndrome, mitochondrial neurogastrointestinal encephalomyopathy, Mohr-Tranebjaerg syndrome, iron-storage disorders related to Friedreich's ataxia, and hereditary spastic paraplegia. We present structural correlations for seven of the mutations. CONCLUSIONS: Of the 19 mutations analyzed, 14 involved changes in very highly conserved parts of the affected proteins. Five out of seven structural correlations provided reasonable explanations for the malfunctions. As additional genetic and structural data become available, this methodology can be extended. It has the potential for assisting in identifying new disease-related genes. Furthermore, profiles with structural homologs can generate mechanistic hypotheses concerning the underlying biochemical processes - and why they break down as a result of the mutations.

Why do we need a three-dimensional architecture of the ligand-binding domain of the nuclear 1alpha,25-dihydroxyvitamin D(3) receptor?

Highly specific binding of 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2)D(3)) by vitamin D receptor (VDR), a nuclear transcriptional factor, activates a genomic mechanism that is manifested in the multiple biologic properties of 1alpha,25(OH)(2)D(3). Numerous synthetic analogs of 1alpha,25(OH)(2)D(3) have been employed to study the interaction between 1alpha,25(OH)(2)D(3) and VDR, and to identify structural markers in 1alpha,25(OH)(2)D(3) that are important for VDR-binding. On the other hand the three-dimensional structure of VDR remained elusive till very recently. In the present study we employed affinity labeling (by 1alpha,25-dihydroxyvitamin D(3)-3-bromoacetate, 1alpha,25(OH)(2)D(3)-3-BE) of VDR to identify C(288) as the anchoring residue for the 3-hydroxyl group of 1alpha,25(OH)(2)D(3) inside the ligand-binding domain of VDR (VDR-LBD). In addition we carried out mutation/hormone-binding analyses to determine the importance of M(284) and W(286) toward hormone binding. We incorporated this information with the three-dimensional structure of the LBD of progesterone receptor to develop a homology-extension model of VDR-LBD. This model identified several amino acid residues as ligand-contact points inside the LBD. Mutational and hormone-binding analyses of these residues verified the structure-functional authenticity of this model, in comparison with the crystal structure of VDR, bound to 1alpha,25(OH)(2)D(3).

Molecular modeling, affinity labeling, and site-directed mutagenesis define the key points of interaction between the ligand-binding domain of the vitamin D nuclear receptor and 1 alpha,25-dihydroxyvitamin D3.

We have combined molecular modeling and classical structure-function techniques to define the interactions between the ligand-binding domain (LBD) of the vitamin D nuclear receptor (VDR) and its natural ligand, 1alpha,25-dihydroxyvitamin D(3) [1alpha,25-(OH)(2)D(3)]. The affinity analogue 1alpha,25-(OH)(2)D(3)-3-bromoacetate exclusively labeled Cys-288 in the VDR-LBD. Mutation of C288 to glycine abolished this affinity labeling, whereas the VDR-LBD mutants C337G and C369G (other conserved cysteines in the VDR-LBD) were labeled similarly to the wild-type protein. These results revealed that the A-ring 3-OH group docks next to C288 in the binding pocket. We further mutated M284 and W286 (separately creating M284A, M284S, W286A, and W286F) and caused severe loss of ligand binding, indicating the crucial role played by the contiguous segment between M284 and C288. Alignment of the VDR-LBD sequence with the sequences of nuclear receptor LBDs of known 3-D structure positioned M284 and W286 in the presumed beta-hairpin of the molecule, thereby identifying it as the region contacting the A-ring of 1alpha, 25-(OH)(2)D(3). From the multiple sequence alignment, we developed a homologous extension model of the VDR-LBD. The model has a canonical nuclear receptor fold with helices H1-H12 and a single beta hairpin but lacks the long insert (residues 161-221) between H2 and H3. We docked the alpha-conformation of the A-ring into the binding pocket first so as to incorporate the above-noted interacting residues. The model predicts hydrogen bonding contacts between ligand and protein at S237 and D299 as well as at the site of the natural mutation R274L. Mutation of S237 or D299 to alanine largely abolished ligand binding, whereas changing K302, a nonligand-contacting residue, to alanine left binding unaffected. In the "activation" helix 12, the model places V418 closest to the ligand, and, consistent with this prediction, the mutation V418S abolished ligand binding. The studies together have enabled us to identify 1alpha,25-(OH)(2)D(3)-binding motifs in the ligand-binding pocket of VDR.

In silico detection of control signals: mRNA 3'-end-processing sequences in diverse species.

We have investigated mRNA 3'-end-processing signals in each of six eukaryotic species (yeast, rice, arabidopsis, fruitfly, mouse, and human) through the analysis of more than 20,000 3'-expressed sequence tags. The use and conservation of the canonical AAUAAA element vary widely among the six species and are especially weak in plants and yeast. Even in the animal species, the AAUAAA signal does not appear to be as universal as indicated by previous studies. The abundance of single-base variants of AAUAAA correlates with their measured processing efficiencies. As found previously, the plant polyadenylation signals are more similar to those of yeast than to those of animals, with both common content and arrangement of the signal elements. In all species examined, the complete polyadenylation signal appears to consist of an aggregate of multiple elements. In light of these and previous results, we present a broadened concept of 3'-end-processing signals in which no single exact sequence element is universally required for processing. Rather, the total efficiency is a function of all elements and, importantly, an inefficient word in one element can be compensated for by strong words in other elements. These complex patterns indicate that effective tools to identify 3'-end-processing signals will require more than consensus sequence identification.

Vitamin D receptor interacts with DnaK/heat shock protein 70: identification of DnaK interaction site on vitamin D receptor.

Vitamin D receptor (VDR) regulates the expression of vitamin D-dependent genes upon binding to its cognate ligand, 1alpha, 25-dihydroxyvitamin D3 (1,25(OH)2D3). This process represents a complex interaction of ligand-bound VDR with nuclear proteins like retinoid X receptor, nuclear accessory factors, and regulatory elements of the target gene. Expression of full-length VDR in Escherichia coli revealed that VDR binds DnaK, a member of heat-shock protein (Hsp) family, with high affinity. By systematic N-terminal truncation of VDR, the interaction site of DnaK on VDR was localized within a 17-amino-acid segment (105-122) representing the "hinge region" between the DNA-binding and hormone-binding domains of VDR. The putative DnaK-binding site was further localized between residues 105 to 109 of VDR by using binding-energy-minimization studies. The interaction of DnaK with VDR did not influence the binding of 1,25(OH)2D3 or nuclear accessory factor(s) to VDR. Furthermore, bovine brain Hsp 70, similar to DnaK, interacted with VDR-ligand-binding domain (105-427). These results suggest that DnaK/Hsp 70 may interact with VDR prior to the activation of the latter by 1,25(OH)2D3-binding.

Genomic detection of new yeast pre-mRNA 3'-end-processing signals.

To investigate Saccharomyces cerevisiae 3'-end-processing signals, a set of 1352 unique pre-mRNA 3'-end-processing sites, corresponding to 861 different genes, was identified by alignment of expressed sequence tag sequences with the complete yeast genome. Nucleotide word frequencies in the vicinity of the cleavage sites were analyzed to reveal the signal element features. In addition to previously recognized processing signals, two previously uncharacterized components of the 3'-end-processing signal sequence were discovered, specifically a predominance of U-rich sequences located on either side of the cleavage site. One of these, the downstream U-rich signal, provides a further link between the 3'-end-processing mechanisms of yeast and higher eukaryotes. Analysis of the complete set of 3'-end-processing sites by means of a discrimination function supports a 'contextual' model in which the sum total effectiveness of the signals in all four elements determines whether or not processing occurs.

Homology model for the ligand-binding domain of the human estrogen receptor.

We have modeled the ligand-binding domain (LBD) of the human estrogen receptor protein (hER) by homology to the known crystal structure of the LBD of the alpha isoform of human retinoate-X receptor (hRX). Alignment of hER with members of the nuclear receptor superfamily defined probable secondary structures which we used to constrain backbone torsion angles and hydrogen bonds. From published studies we identified key interactions between hER and estradiol to use to dock the hormone in its ligand-binding pocket. Since the hRX crystal structure corresponds to the unliganded form of the LBD, we adopted the "mousetrap" mechanism proposed by Renaud et al to predict the structure of the E2-bound hER. Refinement by molecular dynamics and energy minimization gave a model which matches well the known facts about the estradiol phamacophore. It also provides a possible explanation for how hER discriminates between estradiol and testosterone.

Synthesis and characterization of a 29-amino acid residue DNA-binding peptide derived from alpha/beta-type small, acid-soluble spore proteins (SASP) of bacteria.

A 29-amino acid residue peptide (SASP-peptide) derived from the sequence of the putative DNA-contacting portion at the carboxyl terminus of an alpha/beta-type small, acid-soluble spore protein (SASP) of Bacillus subtilis has been synthesized by automated solid-phase methods and tested for its ability to interact with DNA. Circular dichroism (CD) spectroscopy reveals an interaction between this SASP-peptide and DNA, both by an increase in alpha-helix content of the peptide (which alone has a mostly random conformation) and by enhancement of the 275-nm CD band of the DNA. In contrast to results with intact alpha/beta-type SASP, however, the peptide does not induce a B----A conformational transition in the DNA. The SASP-peptide also binds to poly(dG).poly(dC) and protects this polynucleotide against DNase I digestion and UV light-induced cytosine dimer formation, parallel to findings made previously with native alpha/beta-type SASP. The results confirm the hypothesis that the carboxyl-terminal region of the alpha/beta-type SASP directly contacts DNA and possesses some, but not all, of the functional characteristics of the intact molecule.

Binding of small acid-soluble spore proteins from Bacillus subtilis changes the conformation of DNA from B to A.

Small acid-soluble spore proteins (SASPs) appear 3-4 hr after the onset of sporulation in Gram-positive bacteria and constitute up to 20% of the protein of mature spores. Previous studies using Bacillus subtilis deletion mutants lacking SASP-alpha and -beta have shown that such mutations abolish the elevated resistance of spores to UV radiation. Analyses using circular dichroism and Fourier-transform infrared spectroscopy now demonstrate that binding alpha/beta-type SASPs to DNA in vitro causes a structural change in DNA, from the B to the A conformation. This may provide the basis whereby alpha/beta-type SASPs confer increased spore UV resistance in vivo--by changing spore DNA conformation, they alter DNA photochemistry such that UV irradiation produces spore photoproduct instead of the more lethal cyclobutane-type thymine dimers.

A-DNA accommodates adducts derived from diol epoxides of polycyclic aromatic hydrocarbons bound in a "side-stacking" mode.

The minor groove of undistorted A-DNA provides a good binding site for planar, hydrophobic moieties such as unmetabolized polycyclic aromatic hydrocarbons (PAHs), and the base pairs at the ends of short oligodeoxynucleotide helices. It also accommodates the chief adduct derived from the metabolically activated form of the carcinogen benzo[a]pyrene. B-DNA lacks such a site. Computerized models have been generated for the major (N2-guanine-linked) adducts formed at this site by both + and - enantiomers of anti-benzo[a]pyrene-7,8-dihydrodiol-9,10-epoxide (anti-BPDE) with poly(dG).poly(dC) in the A-DNA conformation. The BPDE adducts lie in the shallow, relatively hydrophobic minor groove of the A-DNA after empirical potential energy minimization using the program AMBER. We term this binding mode "side-stacking." The side-stacked + anti-BPDE may constitute the chief carcinogenic lesion derived from benzo[a]pyrene.