Subnanometer-resolution structures of the grass carp reovirus core and virion.

Grass carp reovirus (GCRV) is a member of the Aquareovirus genus of the family Reoviridae, a large family of double-stranded RNA (dsRNA) viruses infecting plants, insects, fishes and mammals. We report the first subnanometer-resolution three-dimensional structures of both GCRV core and virion by cryoelectron microscopy. These structures have allowed the delineation of interactions among the over 1000 molecules in this enormous macromolecular machine and a detailed comparison with other dsRNA viruses at the secondary-structure level. The GCRV core structure shows that the inner proteins have strong structural similarities with those of orthoreoviruses even at the level of secondary-structure elements, indicating that the structures involved in viral dsRNA interaction and transcription are highly conserved. In contrast, the level of similarity in structures decreases in the proteins situated in the outer layers of the virion. The proteins involved in host recognition and attachment exhibit the least similarities to other members of Reoviridae. Furthermore, in GCRV, the RNA-translocating turrets are in an open state and lack a counterpart for the sigma1 protein situated on top of the close turrets observed in mammalian orthoreovirus. Interestingly, the distribution and the organization of GCRV core proteins resemble those of the cytoplasmic polyhedrosis virus, a cypovirus and the structurally simplest member of the Reoviridae family. Our results suggest that GCRV occupies a unique structure niche between the simpler cypoviruses and the considerably more complex mammalian orthoreovirus, thus providing an important model for understanding the structural and functional conservation and diversity of this enormous family of dsRNA viruses.

Characterization of NARJ: A molecular chaperone involved in assembly of nitrate reductase

The major goal of this work was to define the role of accessory protein, NARJ, in assembly of nitrate reductase which is a membrane-bound multisubunit enzyme that can catalyze the reduction of nitrate to nitrite under anaerobic growth in E. coli. Nitrate reductase is encoded by the nar GHJI operon under control of the narG promoter. The purified nitrate reductase is composed of three subunits: $\alpha,\ \beta,$ and $\gamma.$ The NARJ protein which is encoded by the third gene (narJ) is not found to be associated with any of the purified preparations of the enzyme, but is required for active nitrate reductase. In this study the product of the narJ gene was identified. NARJ appeared to be produced at a reduced level, compared to the other proteins encoded by the nar operon. Since NARJ could not be overexpressed to a level for an efficient purification, NARJ was expressed and purified as a recombinant protein with polyhistidine tag. The recombinant protein NARJ-6His could functionally replace native NARJ. Purified NARJ-6His is a dimeric protein which contains no identifiable cofactors or unique secondary structure. NARJ was localized in the cytoplasm, and was not associated with nitrate reductase in the membrane. In vivo NARJ activated the $\alpha\beta$ complex and stabilized the $\alpha$ subunit against protease degradation. In the absence of the membrane-bound $\gamma$ subunit, NARJ formed an intermediate complex with $\alpha\beta$ in the cytosol. Based on these studies, NARJ fits the formal definition of a molecular chaperone. It appears to be required only for the biogenesis of nitrate reductase and, therefore, is defined as a private chaperone specifically involved in the assembly of nitrate reductase system. ^

Study of retinoic acid signaling in cancer cells: Cloning and characterization of retinoic acid responsive genes

Retinoids such as all-trans-retinoic acid (ATRA) are promising agents for cancer chemoprevention and therapy. ATRA can cause growth inhibition, induction of differentiation and apoptosis of a variety of cancer cells. These effects are thought to be mediated by nuclear retinoids receptors which are involved in ligand-dependent transcriptional activation of downstream target genes. Using differential display, we identified several retinoic acid responsive genes in the head and neck squamous carcinoma cells and lung cancer cells, including tissue type transglutaminase, cytochrome P450-related retinoic acid hydroxylase, and a novel gene, designated RAIG1. RAIG1 has two transcripts of 2.4 and 6.8 kbp, respectively, that are generated by alternative selection of polyadenylation sites. Both transcripts have the same open reading frame that encodes a protein comprised of 357 amino acid residues. The deduced RAIG1 protein sequence contains seven transmembrane domains, a signature structure of G protein-coupled receptors. RAIG1 mRNA is expressed at high level in fetal and adult lung tissues. Induction of RAIG1 expression by ATRA is rapid and dose-dependent. A fusion protein of RAIG1 and the green fluorescent protein was localized in the cell surface membrane and perinuclear vesicles in transiently transfected cells. The locus for RAIG1 gene was mapped to a region between D12S358 and D12S847 on chromosome 12p12.3-p13. Our study of the novel retinoic acid induced gene RAIG1 provide evidence for a possible interaction between retinoid and G protein signaling pathways.^ We further examined RAIG1 expression pattern in a panel of 84 cancer cell lines of different origin. The expression level varies greatly from very high to non-detectable. We selected a panel of different cancer cells to study the effects of retinoids and other differentiation agents. We observed: (1) In most cases, retinoids (including all-trans retinoic acid, 4HPR, CD437) could induce the expression of RAIG-1 in cells from cancers of the breast, colon, head and neck, lung, ovarian and prostate. (2) Compare to retinoids, butyrate is often a more potent inducer of RAIG-1 expression in many cancer cells. (3) Butyrate, Phenylacetate butyrate, (R)P-Butyrate and (S)P-Butyrate have different impact on RAIG1 expression which varies among different cell lines. Our results indicate that retinoids could restore RAIG1 expression that is down-regulated in many cancer cells.^ A mouse homologous gene, mRAIG1, was cloned by 5$\sp\prime$ RACE reaction. mRAIG1 cDNA has 2105 bp and shares 63% identity with RAIG1 cDNA. mRAIG1 encodes a polypeptide of 356 amino acid which is 76% identity with RAIG1 protein. mRAIG1 protein also has seven transmembrane domains which are structurally identical to those of RAIG1 protein. Only one 2.2 kbp mRAIG1 transcript could be detected. The mRAIG1 mRNA is also highly expressed in lung tissue. The expression of mRAIG1 gene could be induced by ATRA in several mouse embryonal carcinoma cells. The induction of mRAIG1 expression is associated with retinoic acid-induced neuroectoderm differentiation of P19 cells. Similarity in cDNA and protein sequence, secondary structure, tissue distribution and inducible expression by retinoic acid strongly suggest that the mouse gene is the homologue of the human RAIG1 gene. ^

VERIFICATION OF DNA PREDICTED PROTEIN SEQUENCES BY ENZYME HYDROLYSIS AND MASS SPECTROMETRIC ANALYSIS

The focus of this thesis lies in the development of a sensitive method for the analysis of protein primary structure which can be easily used to confirm the DNA sequence of a protein's gene and determine the modifications which are made after translation. This technique involves the use of dipeptidyl aminopeptidase (DAP) and dipeptidyl carboxypeptidase (DCP) to hydrolyze the protein and the mass spectrometric analysis of the dipeptide products.^ Dipeptidyl carboxypeptidase was purified from human lung tissue and characterized with respect to its proteolytic activity. The results showed that the enzyme has a relatively unrestricted specificity, making it useful for the analysis of the C-terminal of proteins. Most of the dipeptide products were identified using gas chromatography/mass spectrometry (GC/MS). In order to analyze the peptides not hydrolyzed by DCP and DAP, as well as the dipeptides not identified by GC/MS, a FAB ion source was installed on a quadrupole mass spectrometer and its performance evaluated with a variety of compounds.^ Using these techniques, the sequences of the N-terminal and C-terminal regions and seven fragments of bacteriophage P22 tail protein have been verified. All of the dipeptides identified in these analysis were in the same DNA reading frame, thus ruling out the possibility of a single base being inserted or deleted from the DNA sequence. The verification of small sequences throughout the protein sequence also indicates that no large portions of the protein have been removed after translation. ^

Purification and cloning of the novel phosphoprotein copine III and characterization of its associated kinase activity

The copines, named and first described by Creutz et al. (1998), comprise a two C2 domain-containing protein family that can aggregate phosphatidylserine membranes in a calcium-dependent manner. Although no enzymatic function has been attributed to copines, their carboxyl terminus shows homology to the A domain found in integrins that allows binding of magnesium ions. The secondary structure of A domains resembles a Rossmann fold, which can bind dinucleotides and is present in a number of intracellular enzymes. Due to a crossreacting activity of Mik b 1, an antibody to the IL-2R b chain, we were able to serendipitously clone human copine III (CIII). CIII is 65% identical to copine I (CI) and the 5 kb CIII transcript is expressed ubiquitously as determined by a multitissue Northern blot. A polyclonal antibody generated against the carboxyl terminus of CIII recognized CIII in immunoblots and immunoprecipitations. Phosphorylation of CIII was observed on serine and threonine residues, as determined by phosphoamino acid analysis. ^ Experiments were designed to determine whether or not any enzymatic activity, specifically kinase activity, was intrinsic to or associated with CIII. In vitro and in gel kinase assays were performed using transfected HA-tagged CI and CIII, immunoprecipitated endogenous CIII and purified endogenous CIII. The exogenous substrate MBP was phosphorylated in all in vitro kinase assays containing CIII protein purification and column chromatography expertise with me. ^

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Isolation of a myoglobin molten globule by selective cobalt(III)-induced unfolding

Reaction of the Schiff-base complex [Co(acetylacetonate-ethylenediimine)(NH3)2]+ with metmyoglobin at pH 6.5 yields a partially folded protein containing six Co(III) complexes. Although half of its α-helical secondary structure is retained, absorption and CD spectra indicate that the tertiary structure in both B-F and AGH domains is disrupted in the partially folded protein. In analogy to proton-induced unfolding, it is likely that the loss of tertiary structure is triggered by metal-ion binding to histidines. Cobalt(III)-induced unfolding of myoglobin is unique in its selectivity (other proteins are unaffected) and in allowing the isolation of the partially folded macromolecule (the protein does not refold or aggregate upon removal of free denaturant).

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.

Molecular and functional characterization of a novel low-affinity cation transporter (LCT1) in higher plants

The transport of cations across membranes in higher plants plays an essential role in many physiological processes including mineral nutrition, cell expansion, and the transduction of environmental signals. In higher plants the coordinated expression of transport mechanisms is essential for specialized cellular processes and for adaptation to variable environmental conditions. To understand the molecular basis of cation transport in plant roots, a Triticum aestivum cDNA library was used to complement a yeast mutant deficient in potassium (K+) uptake. Two genes were cloned that complemented the mutant: HKT1 and a novel cDNA described in this report encoding a cation transporter, LCT1 (low-affinity cation transporter). Analysis of the secondary structure of LCT1 suggests that the protein contains 8–10 transmembrane helices and a hydrophilic amino terminus containing sequences enriched in Pro, Ser, Thr, and Glu (PEST). The transporter activity was assayed using radioactive isotopes in yeast cells expressing the cDNA. LCT1 mediated low-affinity uptake of the cations Rb+ and Na+, and possibly allowed Ca2+ but not Zn2+ uptake. LCT1 is expressed in low abundance in wheat roots and leaves. The precise functional role of this cation transporter is not known, although the competitive inhibition of cation uptake by Ca2+ has parallels to whole plant and molecular studies that have shown the important role of Ca2+ in reducing Na+ uptake and ameliorating Na+ toxicity. The structure of this higher plant ion transport protein is unique and contains PEST sequences.

Antibody C219 recognizes an α-helical epitope on P-glycoprotein

The ABC transporter, P-glycoprotein, is an integral membrane protein that mediates the ATP-driven efflux of drugs from multidrug-resistant cancer and HIV-infected cells. Anti-P-glycoprotein antibody C219 binds to both of the ATP-binding regions of P-glycoprotein and has been shown to inhibit its ATPase activity and drug binding capacity. C219 has been widely used in a clinical setting as a tumor marker, but recent observations of cross-reactivity with other proteins, including the c-erbB2 protein in breast cancer cells, impose potential limitations in detecting P-glycoprotein. We have determined the crystal structure at a resolution of 2.4 Å of the variable fragment of C219 in complex with an epitope peptide derived from the nucleotide binding domain of P-glycoprotein. The 14-residue peptide adopts an amphipathic α-helical conformation, a secondary structure not previously observed in structures of antibody–peptide complexes. Together with available biochemical data, the crystal structure of the C219-peptide complex indicates the molecular basis of the cross-reactivity of C219 with non-multidrug resistance-associated proteins. Alignment of the C219 epitope with the recent crystal structure of the ATP-binding subunit of histidine permease suggests a structural basis for the inhibition of the ATP and drug binding capacity of P-glycoprotein by C219. The results provide a rationale for the development of C219 mutants with improved specificity and affinity that could be useful in antibody-based P-glycoprotein detection and therapy in multidrug resistant cancers.

Comparative analyses of the secondary structures of synthetic and intracellular yeast MFA2 mRNAs

The overall folded (global) structure of mRNA may be critical to translation and turnover control mechanisms, but it has received little experimental attention. Presented here is a comparative analysis of the basic features of the global secondary structure of a synthetic mRNA and the same intracellular eukaryotic mRNA by dimethyl sulfate (DMS) structure probing. Synthetic MFA2 mRNA of Saccharomyces cerevisiae first was examined by using both enzymes and chemical reagents to determine single-stranded and hybridized regions; RNAs with and without a poly(A) tail were compared. A folding pattern was obtained with the aid of the mfold program package that identified the model that best satisfied the probing data. A long-range structural interaction involving the 5′ and 3′ untranslated regions and causing a juxtaposition of the ends of the RNA, was examined further by a useful technique involving oligo(dT)-cellulose chromatography and antisense oligonucleotides. DMS chemical probing of A and C nucleotides of intracellular MFA2 mRNA was then done. The modification data support a very similar intracellular structure. When low reactivity of A and C residues is found in the synthetic RNA, ≈70% of the same sites are relatively more resistant to DMS modification in vivo. A slightly higher sensitivity to DMS is found in vivo for some of the A and C nucleotides predicted to be hybridized from the synthetic structural model. With this small mRNA, the translation process and mRNA-binding proteins do not block DMS modifications, and all A and C nucleotides are modified the same or more strongly than with the synthetic RNA.

Characterization of the KIF3C Neural Kinesin-like Motor from Mouse

Proteins of the kinesin superfamily define a class of microtubule-dependent motors that play crucial roles in cell division and intracellular transport. To study the molecular mechanism of axonal transport, a cDNA encoding a new kinesin-like protein called KIF3C was cloned from a mouse brain cDNA library. Sequence and secondary structure analysis revealed that KIF3C is a member of the KIF3 family. In contrast to KIF3A and KIF3B, Northern and Western analysis indicated that KIF3C expression is highly enriched in neural tissues such as brain, spinal cord, and retina. When anti-KIF3C antibodies were used to stain the cerebellum, the strongest signal came from the cell bodies and dendrites of Purkinje cells. In retina, anti-KIF3C mainly stains the ganglion cells. Immunolocalization showed that the KIF3C motor in spinal cord and sciatic nerve is mainly localized in cytoplasm. In spinal cord, the KIF3C staining was punctate; double labeling with anti-giantin and anti-KIF3C showed a clear concentration of the motor protein in the Golgi complex. Staining of ligated sciatic nerves demonstrated that the KIF3C motor accumulated at the proximal side of the ligated nerve, which suggests that KIF3C is an anterograde motor. Immunoprecipitation experiments revealed that KIF3C and KIF3A, but not KIF3B, were coprecipitated. These data, combined with previous data from other labs, indicate that KIF3C and KIF3B are “variable” subunits that associate with a common KIF3A subunit, but not with each other. Together these results suggest that KIF3 family members combinatorially associate to power anterograde axonal transport.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

Did homeodomain proteins duplicate before the origin of angiosperms, fungi, and metazoa?

Homeodomain proteins are transcription factors that play a critical role in early development in eukaryotes. These proteins previously have been classified into numerous subgroups whose phylogenetic relationships are unclear. Our phylogenetic analysis of representative eukaryotic sequences suggests that there are two major groups of homeodomain proteins, each containing sequences from angiosperms, metazoa, and fungi. This result, based on parsimony and neighbor-joining analyses of primary amino acid sequences, was supported by two additional features of the proteins. The two protein groups are distinguished by an insertion/deletion in the homeodomain, between helices I and II. In addition, an amphipathic alpha-helical secondary structure in the region N terminal of the homeodomain is shared by angiosperm and metazoan sequences in one group. These results support the hypothesis that there was at least one duplication of homeobox genes before the origin of angiosperms, fungi, and metazoa. This duplication, in turn, suggests that these proteins had diverse functions early in the evolution of eukaryotes. The shared secondary structure in angiosperm and metazoan sequences points to an ancient conserved functional domain.

Protein–RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.