Mutational analysis and molecular modeling of the nonapeptide hormone binding domains of the [Arg8]vasotocin receptor.

To identify determinants that form nonapeptide hormone binding domains of the white sucker Catostomus commersoni [Arg8]vasotocin receptor, chimeric constructs encoding parts of the vasotocin receptor and parts of the isotocin receptor have been analyzed by [(3,5-3H)Tyr2, Arg8]vasotocin binding to membranes of human embryonic kidney cells previously transfected with the different cDNA constructs and by functional expression studies in Xenopus laevis oocytes injected with mutant cRNAs. The results indicate that the N terminus and a region spanning the second extracellular loop and its flanking transmembrane segments, which contains a number of amino acid residues that are conserved throughout the nonapeptide receptor family, contribute to the affinity of the receptor for its ligand. Nonapeptide selectivity, however, is mainly defined by transmembrane region VI and the third extracellular loop. These results are complemented by a molecular model of the vasotocin receptor obtained by aligning its sequence with those of other G-protein coupled receptors as well as that of bacteriorhodopsin. The model indicates that amino acid residues of transmembrane regions II-VII that are located close to the extracellular surface also contribute to the binding of vasotocin.

Molecular cloning and characterization of an invertebrate cellular retinoic acid binding protein

We have cloned a cDNA and gene from the tobacco hornworm, Manduca sexta, which is related to the vertebrate cellular retinoic acid binding proteins (CRABPs). CRABPs are members of the superfamily of lipid binding proteins (LBPs) and are thought to mediate the effects of retinoic acid (RA) on morphogenesis, differentiation, and homeostasis. This discovery of a Manduca sexta CRABP (msCRABP) demonstrates the presence of a CRABP in invertebrates. Compared with bovine/murine CRABP I, the deduced amino acid sequence of msCRABP is 71% homologous overall and 88% homologous for the ligand binding pocket. The genomic organization of msCRABP is conserved with other CRABP family members and the larger LBP superfamily. Importantly, the promoter region contains a motif that resembles an RA response element characteristic of the promoter region of most CRABPs analyzed. Three-dimensional molecular modeling based on postulated structural homology with bovine/murine CRABP I shows msCRABP has a ligand binding pocket that can accommodate RA. The existence of an invertebrate CRABP has significant evolutionary implications, suggesting CRABPs appeared during the evolution of the LBP superfamily well before vertebrate/invertebrate divergence, instead of much later in evolution in selected vertebrates.

Molecular basis for the recognition of two structurally different major histocompatibility complex/peptide complexes by a single T-cell receptor

2C is a typical alloreactive cytotoxic T lymphocyte clone that recognizes two different ligands. These ligands are adducts of the allo-major histocompatibility complex (MHC) molecule H-2Ld and an endogenous octapeptide, and of the self-MHC molecule H-2Kb and another peptide. MHC-binding and T-cell assays with synthetic peptides in combination with molecular modeling studies were employed to analyze the structural basis for this crossreactivity. The molecular surfaces of the two complexes differ greatly in densities and distributions of positive and negative charges. However, modifications of the peptides that increase similarity decrease the capacities of the resulting MHC peptide complexes to induce T-cell responses. Moreover, the roles of the peptides in ligand recognition are different for self- and allo-MHC-restricted T-cell responses. The self-MHC-restricted T-cell responses were finely tuned to recognition of the peptide. The allo-MHC-restricted responses, on the other hand, largely ignore modifications of the peptide. The results strongly suggest that adaptation of the T-cell receptor to the different ligand structures, rather than molecular mimicry by the ligands, is the basis for the crossreactivity of 2C. This conclusion has important implications for T-cell immunology and for the understanding of immunological disorders.

A combined kinetic and modeling study of the catalytic center subsites of human angiogenin.

Kinetic analysis and molecular modeling have been used to map the ribonucleolytic center of angiogenin (Ang). Pyrimidine nucleotides were found to interact very weakly with Ang, consistent with the inaccessible B1 pyrimidine binding site revealed by x-ray crystallography. Ang also lacks an effective phosphate binding site on the 5' side of B1. Although the B2 site that preferentially binds purines on the 3' side of B1 is also weak, its associated phosphate subsites make substantial contributions: both 3',5'-ADP and 5'-ADP have Ki values 6-fold lower than for 5'-AMP, and adding a 3'-phosphate to the substrate CpA increases Kcat/Km by 9-fold. Thus Ang has a functional P2 site on the 3' side of B2 and a site for a second phosphate on the 5' side of B2. Modeling of an Ang-d(ApTpApA) complex suggested that Arg-5 forms part of the P2 site and that a 2'-phosphate might bind more tightly than a 3'-phosphate. Both predictions were confirmed kinetically. The subsite map obtained by this combined approach indicated that 5'-diphosphoadenosine 2'-phosphate might be a more potent inhibitor than any of the nucleotides tested thus far. Indeed, its Ki value of 150 microM is 50-fold lower than that for the best nucleotide previously reported and 400-fold lower than the Km for the best dinucleotide substrate. This compound may serve as a suitable starting point for the eventual design of tight-binding inhibitors of Ang as antiangiogenic agents for human therapy.

Similar antigenic surfaces, rather than sequence homology, dictate T-cell epitope molecular mimicry.

Molecular mimicry, normally defined by the level of primary-sequence similarities between self and foreign antigens, has been considered a key element in the pathogenesis of autoimmunity. Here we describe an example of molecular mimicry between two overlapping peptides within a single self-antigen, both of which are recognized by the same human self-reactive T-cell clone. Two intervening peptides did not stimulate the T-cell clone, even though they share nine amino acids with the stimulatory peptides. Molecular modeling of major histocompatibility complex class II-peptide complexes suggests that both of the recognized peptides generate similar antigenic surfaces, although these are composed of different sets of amino acids. The molecular modeling of a peptide shifted one residue from the stimulatory peptide, which was recognized in the context of the same HLA molecule by another T-cell clone, generated a completely different antigenic surface. Functional studies using truncated peptides confirmed that the anchor residues of the two "mimicking" epitopes in the HLA groove differ. Our results show, for two natural epitopes, how molecular mimicry can occur and suggest that studies of potential antigenic surfaces, rather than sequence similarity, are necessary for analyzing suspected peptide mimicry.

A molecular anchor for stabilizing triple-helical DNA.

Molecular modeling has been used to predict that 2,6-disubstituted amidoanthraquinones, and not the 1,4 series, should preferentially interact with and stabilize triple-stranded DNA structures over duplex DNA. This is due to marked differences in the nature of chromophore-base stacking and groove accessibility for the two series. A DNA foot-printing method that monitors the extent of protection from DNase I cleavage on triplex formation has been used to examine the effects of a number of synthetic isomer compounds in the 1,4 and 2,6 series. The experimental results are in accord with the predicted behavior and confirm that the 1,4 series bind preferentially to double- rather than triple-stranded DNA, whereas the isomeric 2,6 derivatives markedly favor binding to triplex DNA.

Effect of mutating the regulatory phosphoserine and conserved threonine on the activity of the expressed catalytic domain of Acanthamoeba myosin I heavy chain kinase

Phosphorylation of Ser-627 is both necessary and sufficient for full activity of the expressed 35-kDa catalytic domain of myosin I heavy chain kinase (MIHCK). Ser-627 lies in the variable loop between highly conserved residues DFG and APE at a position at which a phosphorylated Ser/Thr also occurs in many other Ser/Thr protein kinases. The variable loop of MIHCK contains two other hydroxyamino acids: Thr-631, which is conserved in almost all Ser/Thr kinases, and Thr-632, which is not conserved. We determined the effects on the kinase activity of the expressed catalytic domain of mutating Ser-627, Thr-631, and Thr-632 individually to Ala, Asp, and Glu. The S627A mutant was substantially less active than wild type (wt), with a lower kcat and higher Km for both peptide substrate and ATP, but was more active than unphosphorylated wt. The S627D and S627E mutants were also less active than phosphorylated wt, i.e., acidic amino acids cannot substitute for phospho-Ser-627. The activity of the T631A mutant was as low as that of the S627A mutant, whereas the T632A mutant was as active as phosphorylated wt, indicating that highly conserved Thr-631, although not phosphorylated, is essential for catalytic activity. Asp and Glu substitutions for Thr-631 and Thr-632 were inhibitory to various degrees. Molecular modeling indicated that Thr-631 can hydrogen bond with conserved residue Asp-591 in the catalytic loop and that similar interactions are possible for other kinases whose activities also are regulated by phosphorylation in the variable loop. Thus, this conserved Thr residue may be essential for the activities of other Ser/Thr protein kinases as well as for the activity of MIHCK.

Activating mutations in the extracellular domain of the fibroblast growth factor receptor 2 function by disruption of the disulfide bond in the third immunoglobulin-like domain

Multiple human skeletal and craniosynostosis disorders, including Crouzon, Pfeiffer, Jackson–Weiss, and Apert syndromes, result from numerous point mutations in the extracellular region of fibroblast growth factor receptor 2 (FGFR2). Many of these mutations create a free cysteine residue that potentially leads to abnormal disulfide bond formation and receptor activation; however, for noncysteine mutations, the mechanism of receptor activation remains unclear. We examined the effect of two of these mutations, W290G and T341P, on receptor dimerization and activation. These mutations resulted in cellular transformation when expressed as FGFR2/Neu chimeric receptors. Additionally, in full-length FGFR2, the mutations induced receptor dimerization and elevated levels of tyrosine kinase activity. Interestingly, transformation by the chimeric receptors, dimerization, and enhanced kinase activity were all abolished if either the W290G or the T341P mutation was expressed in conjunction with mutations that eliminate the disulfide bond in the third immunoglobulin-like domain (Ig-3). These results demonstrate a requirement for the Ig-3 cysteine residues in the activation of FGFR2 by noncysteine mutations. Molecular modeling also reveals that noncysteine mutations may activate FGFR2 by altering the conformation of the Ig-3 domain near the disulfide bond, preventing the formation of an intramolecular bond. This allows the unbonded cysteine residues to participate in intermolecular disulfide bonding, resulting in constitutive activation of the receptor.

RNA-controlled polymorphism in the in vivo assembly of 180-subunit and 120-subunit virions from a single capsid protein

Repeated, specific interactions between capsid protein (CP) subunits direct virus capsid assembly and exemplify regulated protein–protein interactions. The results presented here reveal a striking in vivo switch in CP assembly. Using cryoelectron microscopy, three-dimensional image reconstruction, and molecular modeling, we show that brome mosaic virus (BMV) CP can assemble in vivo two remarkably distinct capsids that selectively package BMV-derived RNAs in the absence of BMV RNA replication: a 180-subunit capsid indistinguishable from virions produced in natural infections and a previously unobserved BMV capsid type with 120 subunits arranged as 60 CP dimers. Each such dimer contains two CPs in distinct, nonequivalent environments, in contrast to the quasi-equivalent CP environments throughout the 180-subunit capsid. This 120-subunit capsid utilizes most of the CP interactions of the 180-subunit capsid plus nonequivalent CP–CP interactions. Thus, the CP of BMV, and perhaps other viruses, can encode CP–CP interactions that are not apparent from mature virions and may function in assembly or disassembly. Shared structural features suggest that the 120- and 180-subunit capsids share assembly steps and that a common pentamer of CP dimers may be an important assembly intermediate. The ability of a single CP to switch between distinct capsids by means of alternate interactions also implies reduced evolutionary barriers between different capsid structures. The in vivo switch between alternate BMV capsids is controlled by the RNA packaged: a natural BMV genomic RNA was packaged in 180-subunit capsids, whereas an engineered mRNA containing only the BMV CP gene was packaged in 120-subunit capsids. RNA features can thus direct the assembly of a ribonucleoprotein complex between alternate structural pathways.

Stretched and overwound DNA forms a Pauling-like structure with exposed bases

We investigate structural transitions within a single stretched and supercoiled DNA molecule. With negative supercoiling, for a stretching force >0.3 pN, we observe the coexistence of B-DNA and denatured DNA from σ ≈ −0.015 down to σ = −1. Surprisingly, for positively supercoiled DNA (σ > +0.037) stretched by 3 pN, we observe a similar coexistence of B-DNA and a new, highly twisted structure. Experimental data and molecular modeling suggest that this structure has ≈2.62 bases per turn and an extension 75% larger than B-DNA. This structure has tightly interwound phosphate backbones and exposed bases in common with Pauling’s early DNA structure [Pauling, L. & Corey, R. B. (1953), Proc. Natl. Acad. Sci. USA 39, 84–97] and an unusual structure proposed for the Pf1 bacteriophage [Liu, D. J. & Day, L. A. (1994) Science 265, 671–674].

Nucleic acid duplexes incorporating a dissociable covalent base pair

We have used molecular modeling techniques to design a dissociable covalently bonded base pair that can replace a Watson-Crick base pair in a nucleic acid with minimal distortion of the structure of the double helix. We introduced this base pair into a potential precursor of a nucleic acid double helix by chemical synthesis and have demonstrated efficient nonenzymatic template-directed ligation of the free hydroxyl groups of the base pair with appropriate short oligonucleotides. The nonenzymatic ligation reactions, which are characteristic of base paired nucleic acid structures, are abolished when the covalent base pair is reduced and becomes noncoplanar. This suggests that the covalent base pair linking the two strands in the duplex is compatible with a minimally distorted nucleic acid double-helical structure.

A functional genetic screen identifies regions at the C-terminal tail and death-domain of death-associated protein kinase that are critical for its proapoptotic activity

Death-associated protein kinase (DAP-kinase) is a Ca+2/calmodulin-regulated serine/threonine kinase with a multidomain structure that participates in apoptosis induced by a variety of signals. To identify regions in this protein that are critical for its proapoptotic activity, we performed a genetic screen on the basis of functional selection of short DAP-kinase-derived fragments that could protect cells from apoptosis by acting in a dominant-negative manner. We expressed a library of randomly fragmented DAP-kinase cDNA in HeLa cells and treated these cells with IFN-γ to induce apoptosis. Functional cDNA fragments were recovered from cells that survived the selection, and those in the sense orientation were examined further in a secondary screen for their ability to protect cells from DAP-kinase-dependent tumor necrosis factor-α-induced apoptosis. We isolated four biologically active peptides that mapped to the ankyrin repeats, the “linker” region, the death domain, and the C-terminal tail of DAP-kinase. Molecular modeling of the complete death domain provided a structural basis for the function of the death-domain-derived fragment by suggesting that the protective fragment constitutes a distinct substructure. The last fragment, spanning the C-terminal serine-rich tail, defined a new regulatory region. Ectopic expression of the tail peptide (17 amino acids) inhibited the function of DAP-kinase, whereas removal of this region from the complete protein caused enhancement of the killing activity, indicating that the C-terminal tail normally plays a negative regulatory role. Altogether, this unbiased screen highlighted functionally important regions in the protein and revealed an additional level of regulation of DAP-kinase apoptotic function that does not affect the catalytic activity.

p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32–36°, whereas wt p53 bends it by 51–57°. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by ≈35° and ≈70°, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53–DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

Database resources of the National Center for Biotechnology Information

In addition to maintaining the GenBank® nucleic acid sequence database, the National Center for Biotechnology Information (NCBI) provides data analysis and retrieval resources that operate on the data in GenBank and a variety of other biological data made available through NCBI’s Web site. NCBI data retrieval resources include Entrez, PubMed, LocusLink and the Taxonomy Browser. Data analysis resources include BLAST, Electronic PCR, OrfFinder, RefSeq, UniGene, HomoloGene, Database of Single Nucleotide Polymorphisms (dbSNP), Human Genome Sequencing, Human MapViewer, GeneMap’99, Human–Mouse Homology Map, Cancer Chromosome Aberration Project (CCAP), Entrez Genomes, Clusters of Orthologous Groups (COGs) database, Retroviral Genotyping Tools, Cancer Genome Anatomy Project (CGAP), SAGEmap, Gene Expression Omnibus (GEO), Online Mendelian Inheri­tance in Man (OMIM), the Molecular Modeling Database (MMDB) and the Conserved Domain Database (CDD). Augmenting many of the Web applications are custom implementations of the BLAST program optimized to search specialized data sets. All of the resources can be accessed through the NCBI home page at: http://www.ncbi.nlm.nih.gov.

Selenoprotein oxidoreductase with specificity for thioredoxin and glutathione systems

Thioredoxin (Trx) and glutathione (GSH) systems are considered to be two major redox systems in animal cells. They are reduced by NADPH via Trx reductase (TR) or oxidized GSH (GSSG) reductase and further supply electrons for deoxyribonucleotide synthesis, antioxidant defense, and redox regulation of signal transduction, transcription, cell growth, and apoptosis. We cloned and characterized a pyridine nucleotide disulfide oxidoreductase, Trx and GSSG reductase (TGR), that exhibits specificity for both redox systems. This enzyme contains a selenocysteine residue encoded by the TGA codon. TGR can reduce Trx, GSSG, and a GSH-linked disulfide in in vitro assays. This unusual substrate specificity is achieved by an evolutionary conserved fusion of the TR and glutaredoxin domains. These observations, together with the biochemical probing and molecular modeling of the TGR structure, suggest a mechanism whereby the C-terminal selenotetrapeptide serves a role of a protein-linked GSSG and shuttles electrons from the disulfide center within the TR domain to either the glutaredoxin domain or Trx.