Majority Of Substitutions Research Articles

Recurrent somatic mutations reveal new insights into consequences of mutagenic processes in cancer.

The sheer size of the human genome makes it improbable that identical somatic mutations at the exact same position are observed in multiple tumours solely by chance. The scarcity of cancer driver mutations also precludes positive selection as the sole explanation. Therefore, recurrent mutations may be highly informative of characteristics of mutational processes. To explore the potential, we use recurrence as a starting point to cluster >2,500 whole genomes of a pan-cancer cohort. We describe each genome with 13 recurrence-based and 29 general mutational features. Using principal component analysis we reduce the dimensionality and create independent features. We apply hierarchical clustering to the first 18 principal components followed by k-means clustering. We show that the resulting 16 clusters capture clinically relevant cancer phenotypes. High levels of recurrent substitutions separate the clusters that we link to UV-light exposure and deregulated activity of POLE from the one representing defective mismatch repair, which shows high levels of recurrent insertions/deletions. Recurrence of both mutation types characterizes cancer genomes with somatic hypermutation of immunoglobulin genes and the cluster of genomes exposed to gastric acid. Low levels of recurrence are observed for the cluster where tobacco-smoke exposure induces mutagenesis and the one linked to increased activity of cytidine deaminases. Notably, the majority of substitutions are recurrent in a single tumour type, while recurrent insertions/deletions point to shared processes between tumour types. Recurrence also reveals susceptible sequence motifs, including TT[C>A]TTT and AAC[T>G]T for the POLE and ‘gastric-acid exposure’ clusters, respectively. Moreover, we refine knowledge of mutagenesis, including increased C/G deletion levels in general for lung tumours and specifically in midsize homopolymer sequence contexts for microsatellite instable tumours. Our findings are an important step towards the development of a generic cancer diagnostic test for clinical practice based on whole-genome sequencing that could replace multiple diagnostics currently in use.

P081 Intron variation: The new allele pattern of 66 novel HLA class I alleles

Aim The aim of this study was to characterize the complete class I gene sequences of new class I alleles found in a population study of 1472 volunteers from the Argentina unrelated donor registry. Methods HLA-A, -B, -C alleles were amplified with primers annealing in the 5’ and 3’ UTRs and the sequences of exons and introns determined by next generation sequencing (NGS) with an Illumina Miseq using an in–house protocol and Conexio Assign software analysis. Initially class I and II amplicons were combined for sequencing; single locus sequencing of a second PCR amplicon was used to confirm variation. Allele sequences were submitted to the IPD-IMGT/HLA database. Results 1% of the 2944 alleles identified at each class I locus were novel. Three novel alleles were found multiple times. The majority (92%) differed by only a single nucleotide substitution from the known allele. Variants in noncoding regions predominated (79%) (T1). None of the intron variants appeared to impact RNA splicing motifs. The distribution of novel substitutions among introns is listed in T2. Unique substitutions, not shared by any other known allele at the locus, were found in 71% HLA-A, 58% HLA-B and 67% HLA-C novel alleles. Conclusions Sixty-six novel class I alleles were identified. The majority of substitutions were found in introns and appear to be random rather than resulting from gene conversion events. Download : Download high-res image (295KB) Download : Download full-size image Download : Download high-res image (73KB) Download : Download full-size image

Unraveling the structural basis of grazoprevir potency against clinically relevant substitutions in hepatitis C virus NS3/4A protease from genotype 1a

Grazoprevir is a potent pan-genotype and macrocyclic inhibitor of hepatitis C virus (HCV) NS3/4A protease and was developed for treating chronic HCV infection. In HCV genotype (GT) 1a, grazoprevir maintains potent activity against a majority of NS3 resistance-associated amino acid substitutions, including the highly prevalent and naturally occurring Q80K polymorphism that impacts simeprevir, another NS3/4A protease inhibitor. The basis for an unexpected difference in the clinical impact of some NS3 substitutions was investigated. Phenotypic analysis of resistance-associated substitutions identified in NS3 from GT1a-infected patients who failed therapy with grazoprevir (in combination with elbasvir, an inhibitor of HCV NS5A protein) showed that positions 56, 156, and 168 in NS3 were most impactful because they diminished protein-inhibitor interactions. Although an amino acid substitution from aspartic acid to alanine at position 168 (D168A) reduced the potency of grazoprevir, its combination with R155K unexpectedly nullified this effect. Molecular dynamics and free-energy surface studies indicated that Asp-168 is important in anchoring Arg-155 for ligand binding but is not critical for Lys-155 because of the inherent flexibility of its side chain. Moreover, modeling studies supported a strong direct cation-heterocycle interaction between the Lys-155 side chain of the double substitution, R155K/D168A, and the lone pair on the quinoxaline in grazoprevir. This unique interaction provides a structural basis for grazoprevir's higher potency than simeprevir, an inhibitor to which the double substitution confers a significant reduction in potency. Our findings are consistent with the detection of R155K/D168A in NS3 from virologic failures treated with simeprevir but not grazoprevir.

Molecular epidemiological and serological studies of bovine leukemia virus (BLV) infection in Thailand cattle

BLV is the etiological agent of enzootic bovine leucosis. BLV has negative effects on animal health and causes economic losses worldwide. However, epidemiological studies on BLV are relatively unknown in many parts of Asian countries. Thus, this study sought to explore BLV infections in cattle in Thailand to determine the extent of the geographic distribution of BLV and to measure its prevalence rates. For this study, 744 cattle from 11 farms in 9 provinces of Thailand were screened in 2013 and 2014 by ELISA and nested PCR. Of those cattle, 41 BLVs were genetically characterized using 188 BLV gp51 env gene sequences available in GenBank. The BLV prevalence in Thailand was high, ranging from 5.3% to 87.8%, as determined by PCR and 11.0% to 100% as determined by ELISA, according to geographical region. Phylogenetic analysis showed that Thailand BLVs belonged to genotypes 1 and 6 and a new genotype 10, which are sporadically observed across Thailand with a prevalence of 31.7%, 19.5%, and 48.8%, respectively. A significant number of amino acid substitutions were also found in the gp51 sequences, of which unique changes in genotype 10 have not been reported previously. Briefly, the majority of substitutions were confined to CD4+/CD8+ T-cell epitopes, neutralizing domains, and E-D-A epitopes. Those observations indicate that BLV infections in Thailand cattle are prevalent and that the geographic distribution of BLV is dynamic, with a high level of genetic diversity. This distribution implies a long-term BLV infection in cattle populations and the movement of infected cattle. In sum, this study suggests that intensive surveillance and effective prevention strategies are required to determine the prevalence of BLV in Thailand and control continuous infections with BLVs.

Genetic drift of human coronavirus OC43 spike gene during adaptive evolution.

Coronaviruses (CoVs) continuously threaten human health. However, to date, the evolutionary mechanisms that govern CoV strain persistence in human populations have not been fully understood. In this study, we characterized the evolution of the major antigen-spike (S) gene in the most prevalent human coronavirus (HCoV) OC43 using phylogenetic and phylodynamic analysis. Among the five known HCoV-OC43 genotypes (A to E), higher substitution rates and dN/dS values as well as more positive selection sites were detected in the S gene of genotype D, corresponding to the most dominant HCoV epidemic in recent years. Further analysis showed that the majority of substitutions were located in the S1 subunit. Among them, seven positive selection sites were chronologically traced in the temporal evolution routes of genotype D, and six were located around the critical sugar binding region in the N-terminal domain (NTD) of S protein, an important sugar binding domain of CoV. These findings suggest that the genetic drift of the S gene may play an important role in genotype persistence in human populations, providing insights into the mechanisms of HCoV-OC43 adaptive evolution.

Western equine encephalitis virus: evolutionary analysis of a declining alphavirus based on complete genome sequences.

Western equine encephalitis virus (WEEV) is an arbovirus from the genus Alphavirus, family Togaviridae, which circulates in North America between birds and mosquitoes, occasionally causing disease in humans and equids. In recent decades, human infection has decreased dramatically; the last documented human case in North America occurred in 1994, and the virus has not been detected in mosquito pools since 2008. Because limited information exists regarding the evolution of WEEV, we analyzed the genomic sequences of 33 low-passage-number strains with diverse geographic and temporal distributions and performed comprehensive phylogenetic analyses. Our results indicated that WEEV is a highly conserved alphavirus with only approximately 5% divergence in its most variable genes. We confirmed the presence of the previously determined group A and B lineages and further resolved group B into three sublineages. We also observed an increase in relative genetic diversity during the mid-20th century, which correlates with the emergence and cocirculation of several group B sublineages. The estimated WEEV population size dropped in the 1990s, with only the group B3 lineage being sampled in the past 20 years. Structural mapping showed that the majority of substitutions in the envelope glycoproteins occurred at the E2-E2 interface. We hypothesize that an event occurred in the mid-20th century that resulted in the increased genetic diversity of WEEV in North America, followed by genetic constriction due to either competitive displacement by the B3 sublineage or stochastic events resulting from a population decline. Western equine encephalitis virus (WEEV) has caused several epidemics that resulted in the deaths of thousands of humans and hundreds of thousands of equids during the past century. During recent decades, human infection decreased drastically and the virus has not been found in mosquito pools since 2008. Because limited information exists regarding the evolution of WEEV, we analyzed 33 complete genome sequences and conducted comprehensive phylogenetic analyses. We confirmed the presence of two major lineages, one of which diverged into three sublineages. Currently, only one of those sublineages is found circulating in nature. Understanding the evolution of WEEV over the past century provides a unique opportunity to observe an arbovirus that is in decline and to better understand what factors can cause said decline.

Screening ofRB1Alterations in Brazilian Patients With Retinoblastoma and Relatives With Retinoma: Phenotypic and Genotypic Associations

To identify constitutional alterations of the retinoblastoma 1 gene (RB1) in two cohorts of Brazilian patients with retinoblastoma and to analyze genotype-phenotype associations. Molecular screening was carried out by direct sequencing of the 27 RB1 exons and flanking regions in blood DNA of 71 patients with retinoblastoma and 4 relatives with retinoma, and with multiplex ligation-dependent probe amplification (MLPA) in 21 patients. The presumed impact of nucleotide substitutions on the structure of the retinoblastoma protein (pRB) was predicted by Polymorphism Phenotyping-2 (PolyPhen-2). Kaplan-Meier and log-rank test were used for estimating 60-month survival rates. One hundred two nucleotide substitutions were detected, 92 substitutions in 59 patients with retinoblastoma and 10 substitutions in 4 individuals with retinoma. Eight substitutions were novel. The majority of substitutions were intronic (86.2%). More than one substitution was present in 37.3% of patients. Twenty-one duplications and 11 deletions were found in 12 patients; some of which with both types of alterations. Duplications/deletions were found in four patients lacking constitutional alterations when analyzed by sequencing, and in eight patients carrying one or more polymorphic intronic substitutions. The global 60-month survival rate in patients was 91.8% (Confidence Interval95% = 85.0 - 99.1). Significant, lower survival rates were found in extraocular presentation (81.0%) versus intraocular tumors (P = 0.014), first enucleation after 1 month following diagnosis (80.9%) versus earlier first enucleation (P = 0.020), and relapse (100.0%) versus absence of relapse (P = 0.0005). Fifteen substitutions (4 intronic and 11 exonic) were identified as probably or likely pathogenic. Four of these 11 exonic substitutions were novel. Survival rates, however, were not affected by presence of these probably or likely pathogenic alterations, most of which not found in patients with retinoblastoma from other Latin American countries. These differences might be related to the different ethnic composition of the Latin American cohorts. Portuguese Abstract.

Human parainfluenza virus type 2 hemagglutinin-neuramindase gene: sequence and phylogenetic analysis of the Saudi strain Riyadh 105/2009

BackgroundAlthough human parainfluenza type 2 (HPIV-2) virus is an important respiratory pathogen, a little is known about strains circulating in Saudi Arabia.FindingsAmong 180 nasopharyngeal aspirates collected from suspected cases in Riyadh, only one sample (0.56%) was confirmed HPIV-2 positive by nested RT-PCR. The sample that was designated Riyadh 105/2009 was used for sequencing and phylogenetic analysis of the most variable virus gene; the haemagglutinin-neuramindase (HN). Comparison of HN gene of Riyadh 105/2009 strain and the relevant sequences available in GenBank revealed a strong relationship with Oklahoma-94-2009 strain. Phylogenetic analysis indicated four different clusters of HPIV-2 strains (G1-4). Twenty-three amino acid substitutions were recorded for Riyadh 105/2009, from which four are unique. The majority of substitutions (n=18) had changed their amino acids characteristics. By analyzing the effect of the recorded substitutions on the protein function using SIFT program, only two located at positions 360 and 571 were predicted to be deleterious.ConclusionsThe presented changes of Riyadh 105/2009 strain may possess potential effect on the protein structure and/or function level. This is the first report that describes partial characterization of Saudi HPIV-2 strain.

Mutagenesis of the West Nile virus NS2B cofactor domain reveals two regions essential for protease activity

The flavivirus NS2B/NS3 protease has received considerable attention as a target for the development of antiviral compounds. While substrate based inhibitors have been the primary focus to date, an approach focussing on NS2B cofactor displacement could prove to be an effective alternative. To understand better the role of the NS2B cofactor in protease activation, we conducted an alanine mutagenesis screen throughout the 42-residue central cofactor domain (NS2B(51-92)) of West Nile virus (WNV). Two sites critical for proteolytic activity were identified (NS2B(59-62) and NS2B(75-87)), where the majority of substitutions were found to significantly decrease proteolytic activity of a recombinant WNV NS2B/NS3 protease. These findings provide mechanistic insights into the structural and functional role that the cofactor may play in the substrate-bound and free protease complexes as well as providing novel sites for targeting new antiviral inhibitors.

Complete mitochondrial genome sequence of three Tetrahymena species reveals mutation hot spots and accelerated nonsynonymous substitutions in Ymf genes.

The ciliate Tetrahymena, a model organism, contains divergent mitochondrial (Mt) genome with unusual properties, where half of its 44 genes still remain without a definitive function. These genes could be categorized into two major groups of KPC (known protein coding) and Ymf (genes without an identified function). To gain insights into the mechanisms underlying gene divergence and molecular evolution of Tetrahymena (T.) Mt genomes, we sequenced three Mt genomes of T.paravorax, T.pigmentosa, and T.malaccensis. These genomes were aligned and the analyses were carried out using several programs that calculate distance, nucleotide substitution (dn/ds), and their rate ratios (ω) on individual codon sites and via a sliding window approach. Comparative genomic analysis indicated a conserved putative transcription control sequence, a GC box, in a region where presumably transcription and replication initiate. We also found distinct features in Mt genome of T.paravorax despite similar genome organization among these ∼47 kb long linear genomes. Another significant finding was the presence of at least one or more highly variable regions in Ymf genes where majority of substitutions were concentrated. These regions were mutation hotspots where elevated distances and the dn/ds ratios were primarily due to an increase in the number of nonsynonymous substitutions, suggesting relaxed selective constraint. However, in a few Ymf genes, accelerated rates of nonsynonymous substitutions may be due to positive selection. Similarly, on protein level the majority of amino acid replacements occurred in these regions. Ymf genes comprise half of the genes in Tetrahymena Mt genomes, so understanding why they have not been assigned definitive functions is an important aspect of molecular evolution. Importantly, nucleotide substitution types and rates suggest possible reasons for not being able to find homologues for Ymf genes. Additionally, comparative genomic analysis of complete Mt genomes is essential in identifying biologically significant motifs such as control regions.

Spatial Covariation of Mutation and Nonsynonymous Substitution Rates in Vertebrate Mitochondrial Genomes

Mitochondrial genomes encode fundamental subunits of the basic energy producing machinery of eukaryotic cells that are under strong functional constraint. Paradoxically, these genes evolve rapidly in general, and there is substantial variation in evolutionary rates among genes within genomes. In order to investigate spatial variation in selection intensity, we conducted tests of neutrality using ratios of synonymous to nonsynonymous substitutions (dN/dS = omega) on numerous protein gene segments from fishes and mammals. Values of omega were very low for nearly all genomic regions. However, values of both omega and dN varied in a clinal pattern with increasing distance from the light-strand origin of replication. Spatial heterogeneity of nonsynonymous substitution rates exhibits a significantly positive correlation with variation in mutation rates that are related to the mode of mitochondrial DNA replication. The finding that nonsynonymous substitution rates are proportional to mutation rates is expected if a majority of substitutions are selectively neutral or slightly deleterious. Spatial patterns of among-gene variation in nonsynonymous rates were highly similar between fishes and mammals, suggesting that forces governing mitochondrial gene evolution have remained relatively constant over 450 Myr of vertebrate evolution. Conservation of substitution patterns despite major shifts in thermal habit and metabolic demands among taxa implicates a conserved replication mechanism controlling relative mutation rates as a major determinant of mitochondrial protein evolution.

Human MutL homolog (MLH1) function in DNA mismatch repair: a prospective screen for missense mutations in the ATPase domain.

Germline mutations in the DNA mismatch repair (MMR) genes MSH2 and MLH1 are responsible for the majority of hereditary non-polyposis colorectal cancer (HNPCC), an autosomal-dominant early-onset cancer syndrome. Genetic testing of both MSH2 and MLH1 from individuals suspected of HNPCC has revealed a considerable number of missense codons, which are difficult to classify as either pathogenic mutations or silent polymorphisms. To identify novel MLH1 missense codons that impair MMR activity, a prospective genetic screen in the yeast Saccharomyces cerevisiae was developed. The screen utilized hybrid human-yeast MLH1 genes that encode proteins having regions of the yeast ATPase domain replaced by homologous regions from the human protein. These hybrid MLH1 proteins are functional in MMR in vivo in yeast. Mutagenized MLH1 fragments of the human coding region were synthesized by error-prone PCR and cloned directly in yeast by in vivo gap repair. The resulting yeast colonies, which constitute a library of hybrid MLH1 gene variants, were initially screened by semi-quantitative in vivo MMR assays. The hybrid MLH1 genes were recovered from yeast clones that exhibited a MMR defect and sequenced to identify alterations in the mutagenized region. This investigation identified 117 missense codons that conferred a 2-fold or greater decreased efficiency of MMR in subsequent quantitative MMR assays. Notably, 10 of the identified missense codons were equivalent to codon changes previously observed in the human population and implicated in HNPCC. To investigate the effect of all possible codon alterations at single residues, a comprehensive mutational analysis of human MLH1 codons 43 (lysine-43) and 44 (serine-44) was performed. Several amino acid replacements at each residue were silent, but the majority of substitutions at lysine-43 (14/19) and serine-44 (18/19) reduced the efficiency of MMR. The assembled data identifies amino acid substitutions that disrupt MLH1 structure and/or function, and should assist the interpretation of MLH1 genetic tests.

Patterns of Polymorphism in Wheat streak mosaic virus: Sequence Space Explored by a Clade of Closely Related Viral Genotypes Rivals That between the Most Divergent Strains

Nucleotide (nt) sequence polymorphism within a collection of Wheat streak mosaic virus (WSMV) isolates was examined. An ∼1267-nt region encompassing the coat protein (CP) cistron and flanking sequences was amplified by reverse transcription–polymerase chain reaction (RT-PCR) for each of 49 isolates not previously sequenced. Consensus sequences were compiled for each isolate based on sequences derived from three clones per RT-PCR product. Among 59 consensus sequences examined, only two were identical. Clades A–C contained divergent isolates from Mexico (Clade A); the Czech Republic, Hungary, and Russia (Clade B); and Iran (Clade C). Fifty-four closely related consensus sequences of isolates from the U.S. (51 sequences), Canada (1 sequence), and Turkey (2 sequences) comprised Clade D. Pair-wise nt divergence between two of the most distantly related sequences (Sidney 81 of Clade D and El Batán 3 of Clade A) was 20%, representing over half of the variable sites (34.1%) in the entire WSMV data set. Maximum pairwise nt divergence within Clade D was 3.6%, yet the proportion of all variable sites within Clade D (21.4%) was similar to that of the Sidney 81–El Batán 3 pair. Patterns of polymorphism within Clade D and the Sidney 81-El Batán 3 pair were remarkably similar with respect to synonymous, nonsynonymous, and noncoding substitutions, as were the proportions of substitutions as a function of nt position within codons. The majority of substitutions within Clade D were synonymous and randomly distributed throughout the coding region examined, whereas nonsynonymous substitutions exhibited a clumped distribution and mostly occurred within the 5′-proximal portion of the CP cistron. Because over half of the polymorphic sites within Clade D were of allele size class 1, the isolates appear to be evolving independently and in a nondeterministic manner, within the constraints of selection. These results indicate that Clade D has undergone substantial and, most likely, recent divergence with the majority of consensus sequence substitutions potentially neutral with respect to fitness. An estimate of evolution rate suggests that the present diversity within the U.S. population arose in about a century, a timeframe corresponding to the establishment of wheat monoculture in the Great Plains.

Mutation Spectra of Herpes Simplex Virus Type 1 Thymidine Kinase Mutants

To examine whether the exonuclease activity intrinsic to the polymerase (Pol) of herpes simplex virus type 1 can influence the mutational spectra, we applied the denaturing gradient gel electrophoresis (DGGE) system combined with sequencing to characterize thymidine kinase mutants isolated from both the wild-type virus and a mutant deficient in exonuclease activity, Y7. Wild-type viruses produced predominantly frameshift mutations (67%), whereas Y7 replicated a significantly lower proportion of frameshifts (21%; P < 0.005). Furthermore, the majority of substitutions were transitional changes in both groups, although they distributed differently. The implications of these findings are discussed.

Are RNA viruses adapting or merely changing?

. RNA viruses and retroviruses fix substitutions approximately 1 million-fold faster than their hosts. This diversification could represent an inevitable drift under purifying selection, the majority of substitutions being phenotypically neutral. The alternative is to suppose that most fixed mutations are beneficial to the virus, allowing it to keep ahead of the host and/or host population. Here, relative sequence diversification of different proteins encoded by viral genomes is found to be linear. The examples encompass a wide variety of retroviruses and RNA viruses. The smoothness of relative divergence spans quasispeciation following clonal infection, to variation among different isolates of the same virus, to viruses from different species or those associated with different diseases, indicating that the majority of fixed mutations likely reflects drift. This held for both mammalian and plant viruses, indicating that adaptive immunity doesn't necessarily shape the relative accumulation of amino acid substitutions. When compared to their hosts RNA viruses evolution appears conservative.

A polar, solvent-exposed residue can be essential for native protein structure

Background: A large energy gap between the native state and the non-native folded states is required for folding into a unique three-dimensional structure. The features that define this energy gap are not well understood, but can be addressed using de novo protein design. Previously, α 2D, a dimeric four-helix bundle, was designed and shown to adopt a native-like conformation. The high-resolution solution structure revealed that this protein adopted a bisecting U motif. Glu7, a solvent-exposed residue that adopts many conformations in solution, might be involved in defining the unique three-dimensional structure of α 2D. Results: A variety of hydrophobic and polar residues were substituted for Glu7 and the dynamic and thermodynamic properties of the resulting proteins were characterized by analytical ultracentrifugation, circular dichroism spectroscopy, and nuclear magnetic resonance spectroscopy. The majority of substitutions at this solvent-exposed position had little affect on the ability to fold into a dimeric four-helix bundle. The ability to adopt a unique conformation, however, was profoundly modulated by the residue at this position despite the similar free energies of folding of each variant. Conclusions: Although Glu7 is not involved directly in stabilizing the native state of α 2D, it is involved indirectly in specifying the observed fold by modulating the energy gap between the native state and the non-native folded states. These results provide experimental support for hypothetical models arising from lattice simulations of protein folding, and underscore the importance of polar interfacial residues in defining the native conformations of proteins.

Sequence replacements in the central β‐turn of plastocyanin

The role of beta-turns in dictating the structure of a beta-barrel protein is assessed by probing the tolerance of the central beta-turn of poplar plastocyanin to substitution by arbitrary sequences. Native plastocyanin binds copper and is colored bright blue. However, when the wild-type Pro47-Ser48-Gly49-Val50 turn sequence is replaced by arbitrary tetrapeptides, the vast majority (92/98 = 94%) of mutant proteins cannot fold into the native blue structure. Characterization of the colorless mutant proteins demonstrates that the majority of substitutions in this type II beta-turn disrupt the native structure severely. Gross structural changes are indicated by major differences in the CD spectra of the mutants relative to the wild-type protein, and by the much larger apparent size of mutant proteins in gel filtration experiments. These mutant proteins do not bind copper. Furthermore, Cys84 forms a disulfide bond readily in the colorless mutant proteins, indicating that it has moved away from the buried position it occupies in the native copper binding site and has become exposed. These results indicate that the central beta-turn in plastocyanin is not merely a default structure arising in response to the surrounding context; rather, sequence information in this turn plays an active role in dictating the location of a chain reversal in the beta-barrel structure. These findings are discussed in terms of their implications for the folding of natural proteins, as well as the design of de novo proteins.

Control of cationic amino acid transport and retroviral receptor functions in a membrane protein family.

A partial cDNA sequence indicated that the T lymphocyte early-activation gene (Tea) encodes a protein related to the dual-function ecotropic retrovirus receptor/cationic amino acid transporter (ecoR/CAT1), and RNA blots suggested highest Tea expression in T lymphocytes and liver (MacLeod, C.L., Finley, K., Kakuda, D. Kozad, C.A., and Wilkinson, M.F. (1990) Mol. Cell. Biol. 7, 3663-3674). The sequence of full-length Tea cDNA from liver (3683 bases) predicts a 657-amino-acid protein (CAT2 alpha) with 12-14 transmembrane domains. A long (515 base) region with six initiation codons and termination codons precedes the translation start codon. The liver Tea cDNA is identical to Tea cDNA from T lymphocytes (encoding CAT2 beta) with the exception of an apparent alternatively spliced sequence encoding a hydrophilic loop of 43 amino acids. The liver-specific sequence contains unique consensus sites for phosphorylation by cyclic AMP-dependent protein kinase and by protein kinase C. Injection of Xenopus oocytes with CAT2 alpha or CAT2 beta messenger RNA resulted in expression of Na(+)-independent cationic amino acid transport that was detected by current measurements under voltage-clamp. Although the amino acid sequences of the isoforms differ in only 21 of 43 residues with the majority of substitutions being conservative, the apparent affinity of CAT2 beta for arginine uptake was 70-fold higher than the CAT2 alpha isoform (Km 38 microM versus 2.7 mM). Neither isoform functioned as a receptor for ecotropic or amphotropic murine retroviruses. However, CAT1-CAT2 chimeric proteins that contain the first three putative extracellular loops of ecoR/CAT1 functioned as ecotropic receptors despite a diminished capacity to bind the viral envelope glycoprotein. The chimeric proteins also functioned as basic amino acid transporters with substrate affinities corresponding to the CAT2 isoform constituting the carboxyl-terminal portion. These results demonstrate that domains of these transporters can function in chimeric combinations to control viral receptor and transport functions.

Amino acid substitutions in pilin of Pseudomonas aeruginosa. Effect on leader peptide cleavage, amino-terminal methylation, and pilus assembly

A total of 37 separate mutants containing single and multiple amino acid substitutions in the leader and amino-terminal conserved region of the Type IV pilin from Pseudomonas aeruginosa were generated by oligonucleotide-directed mutagenesis. The effect of these substitutions on the secretion, processing, and assembly of the pilin monomers into mature pili was examined. The majority of substitutions in the highly conserved amino-terminal region of the pilin monomer had no effect on piliation. Likewise, substitution of several of the residues within the six amino acid leader sequence did not affect secretion and leader cleavage (processing), including replacement of one or both of the positively charged lysine residues with uncharged or negatively charged amino acids. One characteristic of the Type IV pili is the presence of an amino-terminal phenylalanine after leader peptide cleavage which is N-methylated prior to assembly of pilin monomers into pili. Substitution of the amino-terminal phenylalanine with a number of other amino acids, including polar, hydrophobic, and charged residues, did not affect proper leader cleavage and subsequent assembly into pili. Amino-terminal sequencing showed that the majority of substitute residues were also methylated. Substitution of the glycine residue at the -1 position to the cleavage site resulted in the inability to cleave the prepilin monomers and blocked the subsequent assembly of monomers into pili. These results indicate that despite the high degree of conservation in the amino-terminal sequences of the Type IV pili, N-methylphenylalanine at the +1 position relative to the leader peptide cleavage site is not strictly required for pilin assembly. N-Methylation of the amino acids substituted for phenylalanine was shown to have taken place in four of the five mutants tested, but it remains unclear as to whether pilin assembly is dependent on this modification. Recognition and proper cleavage of the prepilin by the leader peptidase appears to be dependent only on the glycine residue at the -1 position. Cell fractionation experiments demonstrated that pilin isolated from mutants deficient in prepilin processing and/or assembly was found in both inner and outer membrane fractions, indistinguishable from the results seen with the wild type.

Species specificity of human interleukin-3 demonstrated by cloning and expression of the homologous rhesus monkey (Macaca mulatta) gene

To enable preclinical studies on homologous interleukin-3 (IL-3) in primate species, we isolated the gene encoding Rhesus monkey IL-3 (RhIL- 3). The nucleotide sequence of the RhIL-3 gene displayed 92.9% homology with that of the human IL-3 (hIL-3) gene. The isolated RhIL-3 gene encodes a 143-amino acid (aa) precursor polypeptide, nine C-terminal residues shorter than the human protein. Protein homology was found to be 89.5% for the signal peptide (19 aa) and 80.5% for the mature protein (124 aa). Comparison of the human and RhIL-3 coding sequences showed that the majority of substitutions had occurred at amino acid replacement sites indicating a rapid evolution of the IL-3 protein. After expression of a genomic fragment in COS cells, RhIL-3 cDNA was constructed, which enabled large-scale production of the RhIL-3 polypeptide, RhIL-3 produced by Bacillus licheniformis and purified to homogeneity appeared to be approximately 100-fold more effective in stimulating Rhesus monkey hematopoietic progenitors than hIL-3, whereas RhIL-3 and hIL-3 showed comparable stimulatory activity on normal as well as malignant human hematopoietic cells. Thus, the rapid evolution of hIL-3 has resulted in a unidirectional species specificity, which most likely restricts the in vivo effects of hIL-3 in Macaca species.