ABSTRACT
Hemochorial placentas have evolved defense mechanisms to prevent the vertical transmission of viruses to the immunologically underdeveloped fetus. Unlike somatic cells that require pathogen-associated molecular patterns to stimulate interferon production, placental trophoblasts constitutively produce type III interferons (IFNL) through an unknown mechanism. We demonstrate that transcripts of short interspersed nuclear elements (SINEs) embedded in miRNA clusters within the placenta trigger a viral mimicry response that induces IFNL and confers antiviral protection. Alu SINEs within primate-specific chromosome 19 (C19MC) and B1 SINEs within rodent-specific microRNA cluster on chromosome 2 (C2MC) produce dsRNAs that activate RIG-I-like receptors (RLRs) and downstream IFNL production. Homozygous C2MC knockout mouse trophoblast stem (mTS) cells and placentas lose intrinsic IFN expression and antiviral protection, whereas B1 RNA overexpression restores C2MCΔ/Δ mTS cell viral resistance. Our results uncover a convergently evolved mechanism whereby SINE RNAs drive antiviral resistance in hemochorial placentas, placing SINEs as integral players in innate immunity.
Subject(s)
MicroRNAs , Animals , Mice , Female , Pregnancy , MicroRNAs/genetics , Placenta , Interferon Lambda , Antiviral Agents , Short Interspersed Nucleotide ElementsABSTRACT
Despite the central role of the placenta in supporting a pregnancy, relatively little is known about transcriptomic and immune-cell changes that occur across gestation. To generate a reference gene expression map of first (T1), second (T2) and third (T3) trimester human placenta, and assess differences in transcriptome in maternal versus fetal side tissues sections of full-term placenta, we performed RNA-Seq analysis on 64 biopsy samples from 18 placentas across all three gestations. We identified 1120 differentially expressed genes in placenta tissues from T1 and T3 samples using a generalized linear model within DESeq2. In total, 411 and 709 genes were positively associated with T1 and T3 placenta, respectively. Comparison of the top 200 differentially expressed genes in the T1 placenta with T3 showed that most of the top enriched biological processes were related to cell division and proliferation. T1 and T2 tissues shared expression of fibroblast-specific COL6A2, HGF, and SPP1 genes. In T3 samples, the expression of genes relating to vasculature development and regulation were highly enriched. Monocytes and NK cell population increased in T3 compared to T1 and T2, whereas Hofbauer cell proportion expanded significantly in T2 and then decreased in T3 samples. There were no significant gene expression differences in the maternal vs. fetal side in T3 placentas. Gene expression patterns shift temporally across trimesters but not spatially across the placenta, at least at the resolution of the biopsy samples. The genes and gene set we identified here represent a valuable resource for studying pathology in pregnancy-related disorders.
Subject(s)
Placenta , Transcriptome , Female , Humans , Placenta/metabolism , PregnancyABSTRACT
Sensitive and accurate RT-qPCR tests are the primary diagnostic tools to identify SARS-CoV-2-infected patients. While many SARS-CoV-2 RT-qPCR tests are available, there are significant differences in test sensitivity, workflow (e.g. hands-on-time), gene targets and other functionalities that users must consider. Several publicly available protocols shared by reference labs and public health authorities provide useful tools for SARS-CoV-2 diagnosis, but many have shortcomings related to sensitivity and laborious workflows. Here, we describe a series of SARS-CoV-2 RT-qPCR tests that are originally based on the protocol targeting regions of the RNA-dependent RNA polymerase (RdRp) and envelope (E) coding genes developed by the Charité Berlin. We redesigned the primers/probes, utilized locked nucleic acid nucleotides, incorporated dual probe technology and conducted extensive optimizations of reaction conditions to enhance the sensitivity and specificity of these tests. By incorporating an RNase P internal control and developing multiplexed assays for distinguishing SARS-CoV-2 and influenza A and B, we streamlined the workflow to provide quicker results and reduced consumable costs. Some of these tests use modified enzymes enabling the formulation of a room temperature-stable master mix and lyophilized positive control, thus increasing the functionality of the test and eliminating cold chain shipping and storage. Moreover, a rapid, RNA extraction-free version enables high sensitivity detection of SARS-CoV-2 in about an hour using minimally invasive, self-collected gargle samples. These RT-qPCR assays can easily be implemented in any diagnostic laboratory and can provide a powerful tool to detect SARS-CoV-2 and the most common seasonal influenzas during the vaccination phase of the pandemic.
Subject(s)
COVID-19 , Influenza, Human , COVID-19/diagnosis , COVID-19 Testing , Humans , Influenza, Human/diagnosis , Nucleotides , RNA, Viral/analysis , RNA, Viral/genetics , Real-Time Polymerase Chain Reaction/methods , SARS-CoV-2/genetics , Sensitivity and Specificity , TechnologyABSTRACT
BACKGROUND: Extracellular RNAs (exRNAs) in biofluids are amenable to quantitative analysis and proposed as noninvasive biomarkers for monitoring organ function. Cell-lineage-specific microRNAs (miRNAs) are present in plasma as soluble ribonucleoproteins or enclosed in exRNA carriers and transported through the vasculature. However, more extensive studies of healthy individuals are needed to gain insights into the variability of plasma miRNA abundance and composition. METHODS: The exRNA composition of platelet-depleted plasma collected twice from 236 healthy individuals was characterized by small RNA sequencing. Plasma of pregnant women featuring dramatically increased placental miRNAs and samples from subject P12 with noticeably increased epithelial- and neuroendocrine-origin miRNAs were included for comparison. The miRNA content of 10 000g and 100 000g pellet fractions of plasma generated by ultracentrifugation was also determined. Data analysis methods included Pearson correlation, differential gene expression, and unsupervised clustering. RESULTS: The abundance changes for more variable miRNAs in plasma of normal individuals correlated between coexpressed cell-lineage-specific miRNAs of the liver, neuroendocrine organs, epithelial cells, and muscle. ExRNA of pellet fractions contained <2% of total plasma miRNA with modest enrichment of lineage-specific and variable miRNAs compared to supernatant. The abundance fold changes of miRNAs observed in pregnancy and P12 compared to normal exceeded interquartile variability of healthy individuals. The neuroendocrine miRNA signature of P12 persisted for more than 4 years and was absent in other individuals. CONCLUSIONS: This study defines the framework and effect size for screening of extensive plasma collections for miRNA phenotypes and biomarker discovery.
Subject(s)
MicroRNAs , Sequence Analysis, RNA , Biomarkers , Female , Humans , MicroRNAs/blood , MicroRNAs/genetics , Phenotype , Placenta , Pregnancy , Pregnant Women , Sequence Analysis, RNA/methodsABSTRACT
AIMS: Allograft rejection following heart transplantation (HTx) is a serious complication even in the era of modern immunosuppressive regimens and causes up to a third of early deaths after HTx. Allograft rejection is mediated by a cascade of immune mechanisms leading to acute cellular rejection (ACR) and/or antibody-mediated rejection (AMR). The gold standard for monitoring allograft rejection is invasive endomyocardial biopsy that exposes patients to complications. Little is known about the potential of circulating miRNAs as biomarkers to detect cardiac allograft rejection. We here present a systematic analysis of circulating miRNAs as biomarkers and predictors for allograft rejection after HTx using next-generation small RNA sequencing. METHODS AND RESULTS: We used next-generation small RNA sequencing to investigate circulating miRNAs among HTx recipients (10 healthy controls, 10 heart failure patients, 13 ACR, and 10 AMR). MiRNA profiling was performed at different time points before, during, and after resolution of the rejection episode. We found three miRNAs with significantly increased serum levels in patients with biopsy-proven cardiac rejection when compared with patients without rejection: hsa-miR-139-5p, hsa-miR-151a-5p, and hsa-miR-186-5p. We identified miRNAs that may serve as potential predictors for the subsequent development of ACR: hsa-miR-29c-3p (ACR) and hsa-miR-486-5p (AMR). Overall, hsa-miR-486-5p was most strongly associated with acute rejection episodes. CONCLUSIONS: Monitoring cardiac allograft rejection using circulating miRNAs might represent an alternative strategy to invasive endomyocardial biopsy.
Subject(s)
Heart Transplantation , MicroRNAs , Allografts , Biomarkers , Graft Rejection/diagnosis , Humans , MicroRNAs/geneticsABSTRACT
We now know RNA can survive the harsh environment of biofluids when encapsulated in vesicles or by associating with lipoproteins or RNA binding proteins. These extracellular RNA (exRNA) play a role in intercellular signaling, serve as biomarkers of disease, and form the basis of new strategies for disease treatment. The Extracellular RNA Communication Consortium (ERCC) hosted a two-day online workshop (April 19-20, 2021) on the unique challenges of exRNA data analysis. The goal was to foster an open dialog about best practices and discuss open problems in the field, focusing initially on small exRNA sequencing data. Video recordings of workshop presentations and discussions are available (https://exRNA.org/exRNAdata2021-videos/). There were three target audiences: experimentalists who generate exRNA sequencing data, computational and data scientists who work with those groups to analyze their data, and experimental and data scientists new to the field. Here we summarize issues explored during the workshop, including progress on an effort to develop an exRNA data analysis challenge to engage the community in solving some of these open problems.
ABSTRACT
Neuroendocrine neoplasms (NENs) are clinically diverse and incompletely characterized cancers that are challenging to classify. MicroRNAs (miRNAs) are small regulatory RNAs that can be used to classify cancers. Recently, a morphology-based classification framework for evaluating NENs from different anatomical sites was proposed by experts, with the requirement of improved molecular data integration. Here, we compiled 378 miRNA expression profiles to examine NEN classification through comprehensive miRNA profiling and data mining. Following data preprocessing, our final study cohort included 221 NEN and 114 non-NEN samples, representing 15 NEN pathological types and 5 site-matched non-NEN control groups. Unsupervised hierarchical clustering of miRNA expression profiles clearly separated NENs from non-NENs. Comparative analyses showed that miR-375 and miR-7 expression is substantially higher in NEN cases than non-NEN controls. Correlation analyses showed that NENs from diverse anatomical sites have convergent miRNA expression programs, likely reflecting morphological and functional similarities. Using machine learning approaches, we identified 17 miRNAs to discriminate 15 NEN pathological types and subsequently constructed a multilayer classifier, correctly identifying 217 (98%) of 221 samples and overturning one histological diagnosis. Through our research, we have identified common and type-specific miRNA tissue markers and constructed an accurate miRNA-based classifier, advancing our understanding of NEN diversity.
ABSTRACT
Extracellular mRNAs (ex-mRNAs) potentially supersede extracellular miRNAs (ex-miRNAs) and other RNA classes as biomarkers. We performed conventional small-RNA-sequencing (sRNA-seq) and sRNA-seq with T4 polynucleotide kinase (PNK) end-treatment of total exRNA isolated from serum and platelet-poor EDTA, ACD, and heparin plasma to study the effect on ex-mRNA capture. Compared to conventional sRNA-seq PNK-treatment increased the detection of informative ex-mRNAs reads up to 50-fold. The exRNA pool was dominated by hematopoietic cells and platelets, with additional contribution from the liver. About 60% of the 15- to 42-nt reads originated from the coding sequences, in a pattern reminiscent of ribosome-profiling. Blood sample type had a considerable influence on the exRNA profile. On average approximately 350 to 1,100 distinct ex-mRNA transcripts were detected depending on plasma type. In serum, additional transcripts from neutrophils and hematopoietic cells increased this number to near 2,300. EDTA and ACD plasma showed a destabilizing effect on ex mRNA and non-coding RNA ribonucleoprotein complexes compared to other plasma types. In a proof-of-concept study, we investigated differences between the exRNA profiles of patients with acute coronary syndrome (ACS) and healthy controls. The improved tissue resolution of ex mRNAs after PNK-treatment enabled us to detect a neutrophil-signature in ACS that escaped detection by ex miRNA analysis.
Subject(s)
Acute Coronary Syndrome/genetics , Blood Cells/metabolism , Cell-Free Nucleic Acids/genetics , Liver/metabolism , Muscle, Skeletal/metabolism , RNA-Seq/methods , Acute Coronary Syndrome/blood , Adult , Aged , Biomarkers/blood , Case-Control Studies , Cell-Free Nucleic Acids/blood , Citric Acid , Edetic Acid , Erythrocytes/metabolism , Female , Glucose/analogs & derivatives , Heparin , High-Throughput Nucleotide Sequencing , Humans , Leukocytes, Mononuclear/metabolism , Male , Middle Aged , Monocytes/metabolism , Neutrophils/metabolism , Plasma , Polynucleotide 5'-Hydroxyl-Kinase , Proof of Concept Study , Sequence Analysis, RNA/methods , Serum , Specimen HandlingABSTRACT
The placenta and decidua interact dynamically to enable embryonic and fetal development. Here, we report single-cell RNA sequencing of 14,341 and 6754 cells from first-trimester human placental villous and decidual tissues, respectively. Bioinformatic analysis identified major cell types, many known and some subtypes previously unknown in placental villi and decidual context. Further detailed analysis revealed proliferating subpopulations, enrichment of cell type-specific transcription factors, and putative intercellular communication in the fetomaternal microenvironment. This study provides a blueprint to further the understanding of the roles of these cells in the placenta and decidua for maintenance of early gestation as well as pathogenesis in pregnancy-related disorders.
Subject(s)
Biomarkers/analysis , Chorionic Villi/metabolism , Decidua/metabolism , Placenta/metabolism , Pregnancy Trimester, First/genetics , Single-Cell Analysis/methods , Trophoblasts/metabolism , Decidua/cytology , Female , High-Throughput Nucleotide Sequencing , Humans , Placenta/cytology , Pregnancy , Trophoblasts/cytologyABSTRACT
Circulating extracellular RNAs (exRNAs) have the potential to serve as biomarkers for a wide range of medical conditions. However, limitations in existing exRNA isolation methods and a lack of knowledge on parameters affecting exRNA variability in human samples may hinder their successful discovery and clinical implementation. Using combinations of denaturants, reducing agents, proteolysis, and revised organic extraction, we developed an automated, high-throughput approach for recovery of exRNAs and exDNA from the same biofluid sample. We applied this method to characterize exRNAs from 312 plasma and serum samples collected from 13 healthy volunteers at 12 time points over a 2-month period. Small RNA cDNA library sequencing identified nearly twofold increased epithelial-, muscle-, and neuroendocrine-cell-specific miRNAs in females, while fasting and hormonal cycle showed little effect. External standardization helped to detect quantitative differences in erythrocyte and platelet-specific miRNA contributions and in miRNA concentrations between biofluids. It also helped to identify a study participant with a unique exRNA phenotype featuring a miRNA signature of up to 20-fold elevated endocrine-cell-specific miRNAs and twofold elevated total miRNA concentrations stable for over 1 year. Collectively, these results demonstrate an efficient and quantitative method to discern exRNA phenotypes and suggest that plasma and serum RNA profiles are stable over months and can be routinely monitored in long-term clinical studies.
Subject(s)
Cell-Free Nucleic Acids/blood , Adult , Biomarkers/blood , Cell-Free Nucleic Acids/genetics , Cell-Free Nucleic Acids/isolation & purification , Female , Gene Library , High-Throughput Nucleotide Sequencing , Humans , Male , MicroRNAs/blood , MicroRNAs/geneticsABSTRACT
OBJECTIVE: MicroRNA (miRNA) are short nucleotide strands with a regulatory function in the cell. Several miRNAs have been shown to be useful as biomarkers for different neoplasms. The aim of this project was to assess whether levels of miRNA in cell free urine could be used as a biomarker in transitional cell carcinoma (TCC). RESULTS: cDNA libraries were produced based on small RNAs in urine samples of fourteen TCC patients and twenty healthy volunteers. Resulting reads were deep sequenced on Illumina HiSeq sequencer with the intent of characterizing cell free urine miRNA profiles. A statistically significant difference was found for a single miRNA; miR-210 was > sixfold higher in the TCC group compared to the control group. Furthermore, we were able to produce a diagnostic score by summing of standardized levels of overexpressed miRNA. This score was considerably higher in TCC patients with a sensitivity of 0.93, specificity of 0.76 and negative predictive value > 0.97.
Subject(s)
Biomarkers, Tumor/urine , Carcinoma, Transitional Cell/diagnosis , MicroRNAs/urine , Urologic Neoplasms/diagnosis , Carcinoma, Transitional Cell/urine , Case-Control Studies , Cell-Free System , DNA, Complementary/genetics , High-Throughput Nucleotide Sequencing , Humans , Urologic Neoplasms/urineABSTRACT
BACKGROUND: Urine is a potential source of biomarkers for diseases of the kidneys and urinary tract. RNA, including microRNA, is present in the urine enclosed in detached cells or in extracellular vesicles (EVs) or bound and protected by extracellular proteins. Detection of cell- and disease-specific microRNA in urine may aid early diagnosis of organ-specific pathology. In this study, we applied barcoded deep sequencing to profile microRNAs in urine of healthy volunteers, and characterized the effects of sex, urine fraction (cells vs. EVs) and repeated voids by the same individuals. RESULTS: Compared to urine-cell-derived small RNA libraries, urine-EV-derived libraries were relatively enriched with miRNA, and accordingly had lesser content of other small RNA such as rRNA, tRNA and sn/snoRNA. Unsupervised clustering of specimens in relation to miRNA expression levels showed prominent bundling by specimen type (urine cells or EVs) and by sex, as well as a tendency of repeated (first and second void) samples to neighbor closely. Likewise, miRNA profile correlations between void repeats, as well as fraction counterparts (cells and EVs from the same specimen) were distinctly higher than correlations between miRNA profiles overall. Differential miRNA expression by sex was similar in cells and EVs. CONCLUSIONS: miRNA profiling of both urine EVs and sediment cells can convey biologically important differences between individuals. However, to be useful as urine biomarkers, careful consideration is needed for biofluid fractionation and sex-specific analysis, while the time of voiding appears to be less important.
Subject(s)
Biomarkers/urine , Cell-Derived Microparticles/genetics , Extracellular Vesicles/genetics , Gene Expression Profiling , High-Throughput Nucleotide Sequencing/methods , MicroRNAs/genetics , MicroRNAs/urine , Adult , Female , Healthy Volunteers , Humans , Male , Principal Component Analysis , Young AdultABSTRACT
Recent studies implicated the RNA-binding protein with multiple splicing (RBPMS) family of proteins in oocyte, retinal ganglion cell, heart, and gastrointestinal smooth muscle development. These RNA-binding proteins contain a single RNA recognition motif (RRM), and their targets and molecular function have not yet been identified. We defined transcriptome-wide RNA targets using photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) in HEK293 cells, revealing exonic mature and intronic pre-mRNA binding sites, in agreement with the nuclear and cytoplasmic localization of the proteins. Computational and biochemical approaches defined the RNA recognition element (RRE) as a tandem CAC trinucleotide motif separated by a variable spacer region. Similar to other mRNA-binding proteins, RBPMS family of proteins relocalized to cytoplasmic stress granules under oxidative stress conditions suggestive of a support function for mRNA localization in large and/or multinucleated cells where it is preferentially expressed.
Subject(s)
Gene Expression Profiling , Protein Interaction Domains and Motifs , RNA-Binding Proteins/chemistry , RNA-Binding Proteins/metabolism , RNA/metabolism , Adult , Amino Acid Sequence , Animals , Binding Sites , Cells, Cultured , Female , HEK293 Cells , Humans , Male , Mice , Molecular Sequence Data , Multigene Family , Protein Binding , Protein Interaction Domains and Motifs/genetics , RNA/genetics , Sequence Homology, Amino Acid , TranscriptomeABSTRACT
Human LIN28A and LIN28B are RNA-binding proteins (RBPs) conserved in animals with important roles during development and stem cell reprogramming. We used Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation (PAR-CLIP) in HEK293 cells and identified a largely overlapping set of â¼3000 mRNAs at â¼9500 sites located in the 3' UTR and CDS. In vitro and in vivo, LIN28 preferentially bound single-stranded RNA containing a uridine-rich element and one or more flanking guanosines and appeared to be able to disrupt base-pairing to access these elements when embedded in predicted secondary structure. In HEK293 cells, LIN28 protein binding mildly stabilized target mRNAs and increased protein abundance. The top targets were its own mRNAs and those of other RBPs and cell cycle regulators. Alteration of LIN28 protein levels also negatively regulated the abundance of some but not all let-7 miRNA family members, indicating sequence-specific binding of let-7 precursors to LIN28 proteins and competition with cytoplasmic miRNA biogenesis factors.
Subject(s)
MicroRNAs/metabolism , RNA, Messenger , RNA-Binding Proteins , 3' Untranslated Regions , Binding Sites , Cytoplasm/metabolism , Gene Expression Regulation , HEK293 Cells , Humans , Immunoprecipitation , Protein Binding , RNA, Messenger/genetics , RNA, Messenger/metabolism , RNA-Binding Proteins/genetics , RNA-Binding Proteins/metabolism , Stem Cells/metabolismABSTRACT
Bacterial cold shock proteins (CSPs) regulate the cellular response to temperature downshift. Their general principle of function involves RNA chaperoning and transcriptional antitermination. Here we present two crystal structures of cold shock protein B from Bacillus subtilis (Bs-CspB) in complex with either a hexanucleotide (5'-UUUUUU-3') or heptanucleotide (5'-GUCUUUA-3') single-stranded RNA (ssRNA). Hydrogen bonds and stacking interactions between RNA bases and aromatic sidechains characterize individual binding subsites. Additional binding subsites which are not occupied by the ligand in the crystal structure were revealed by NMR spectroscopy in solution on Bs-CspB·RNA complexes. Binding studies demonstrate that Bs-CspB associates with ssDNA as well as ssRNA with moderate sequence specificity. Varying affinities of oligonucleotides are reflected mainly in changes of the dissociation rates. The generally lower binding affinity of ssRNA compared to its ssDNA analog is attributed solely to the substitution of thymine by uracil bases in RNA.
Subject(s)
Bacterial Proteins/chemistry , Heat-Shock Proteins/chemistry , RNA/chemistry , Amino Acid Sequence , Crystallography, X-Ray , DNA, Single-Stranded/chemistry , Hydrogen Bonding , Kinetics , Ligands , Molecular Sequence Data , Oligonucleotides/chemistry , Protein Conformation , Solutions , Thermodynamics , Thymine/chemistryABSTRACT
The coxsackievirus-adenovirus receptor (CAR) is a member of the Ig superfamily strongly expressed in the developing nervous system. Our histological investigations during development reveal an initial uniform distribution of CAR on all neural cells with a concentration on membranes that face the margins of the nervous system (e.g., the basal laminae and the ventricular side). At more advanced stages, CAR becomes downregulated and restricted to specific regions including areas rich in axonal and dendritic surfaces. To study the function of CAR on neural cells, we used the fiber knob of the adenovirus, extracellular CAR domains, blocking antibodies to CAR, as well as CAR-deficient neural cells. Blocking antibodies were found to inhibit neurite extension in retina organ and retinal explant cultures, whereas the application of the recombinant fiber knob of the adenovirus subtype Ad2 or extracellular CAR domains promoted neurite extension and adhesion to extracellular matrices. We observed a promiscuous interaction of CAR with extracellular matrix glycoproteins, which was deduced from analytical ultracentrifugation experiments, affinity chromatography, and adhesion assays. The membrane proximal Ig domain of CAR, termed D2, was found to bind to a fibronectin fragment, including the heparin-binding domain 2, which promotes neurite extension of wild type, but not of CAR-deficient neural cells. In contrast to heterophilic interactions, homophilic association of CAR involves both Ig domains, as was revealed by ultracentrifugation, chemical cross-linking, and adhesion studies. The results of these functional and binding studies are correlated to a U-shaped homodimer of the complete extracellular domains of CAR detected by x-ray crystallography.
Subject(s)
Cell Differentiation/genetics , Central Nervous System/embryology , Central Nervous System/metabolism , Gene Expression Regulation, Developmental/genetics , Neurogenesis/genetics , Neurons/metabolism , Receptors, Virus/metabolism , Animals , Antibodies, Blocking/pharmacology , CHO Cells , Cell Adhesion/genetics , Cells, Cultured , Central Nervous System/cytology , Chick Embryo , Coxsackie and Adenovirus Receptor-Like Membrane Protein , Cricetinae , Cricetulus , Crystallography, X-Ray , Dimerization , Extracellular Matrix Proteins/metabolism , Fibronectins/metabolism , Humans , Mice , Mice, Knockout , NIH 3T3 Cells , Neurites/metabolism , Neurites/ultrastructure , Neurons/cytology , Organ Culture Techniques , Protein Structure, Tertiary/physiology , Receptors, Virus/chemistry , Receptors, Virus/genetics , Retina/cytology , Retina/embryology , Retina/metabolismABSTRACT
In previous work, a strongly stabilized variant of the beta1 domain of streptococcal protein G (Gbeta1) was obtained by an in vitro selection method. This variant, termed Gbeta1-M2, contains the four substitutions E15V, T16L, T18I, and N37L. Here we elucidated the molecular basis of the observed strong stabilizations. The contributions of these four residues were analyzed individually and in various combinations, additional selections with focused Gbeta1 gene libraries were performed, and the crystal structure of Gbeta1-M2 was determined. All single substitutions (E15V, T16L, T18I, and N37L) stabilize wild-type Gbeta1 by contributions of between 1.6 and 6.0 kJ mol(-1) (at 70 degrees C). Hydrophobic residues at positions 16 and 37 provide the major contribution to stabilization by enlarging the hydrophobic core of Gbeta1. They also increase the tendency to form dimers, as shown by dependence on the concentration of apparent molecular mass in analytical ultracentrifugation, by concentration-dependent stability, and by a strongly increased van't Hoff enthalpy of unfolding. The 0.88-A crystal structure of Gbeta1-M2 and NMR measurements in solution provide the explanation for the observed dimer formation. It involves a head-to-head arrangement of two Gbeta1-M2 molecules via six intermolecular hydrogen bonds between the two beta strands 2 and 2' and an adjacent self-complementary hydrophobic surface area, which is created by the T16L and N37L substitutions and a large 120 degrees rotation of the Tyr33 side chain. This removal of hydrophilic groups and the malleability of the created hydrophobic surface provide the basis for the dimer formation of stabilized Gbeta1 variants.
Subject(s)
Bacterial Proteins/chemistry , Protein Structure, Quaternary , Protein Subunits/chemistry , Amino Acid Substitution , Bacterial Proteins/genetics , Crystallography, X-Ray , Dimerization , Models, Molecular , Protein Denaturation , Protein Subunits/geneticsABSTRACT
Computational design and in vitro evolution are major strategies for stabilizing proteins. For the four critical positions 16, 18, 25, and 29 of the B domain of the streptococcal protein G (Gbeta1), they identified the same optimal residues at positions 16 and 25, but not at 18 and 29. Here we analyzed the energetic contributions of the residues from these two approaches by single and double mutant analyses and determined crystal structures for a variant from the calculation (I16/L18/E25/K29) and from the selection (I16/I18/E25/F29). The structural analysis explains the observed differences in stabilization. Residues 16, 18, and 29 line an invagination, which results from a packing defect between the helix and the beta-sheet of Gbeta1. In all stabilized variants, residues with larger side-chains occur at these positions and packing is improved. In the selected variant, packing is better optimized than in the computed variant. Such differences in side-chain packing strongly affect stability but are difficult to evaluate by computation.
Subject(s)
Antigens, Bacterial/chemistry , Bacterial Proteins/chemistry , Computational Biology , Evolution, Molecular , Antigens, Bacterial/genetics , Antigens, Bacterial/metabolism , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Computer Simulation , Crystallography, X-Ray , Models, Molecular , Protein Conformation , Protein Folding , Protein Structure, Secondary , ThermodynamicsABSTRACT
The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.
Subject(s)
Bacterial Proteins/chemistry , Bacterial Proteins/genetics , Cold Temperature , Genetic Variation , Protein Structure, Tertiary , Amino Acid Sequence , Bacterial Proteins/metabolism , Crystallography, X-Ray , Directed Molecular Evolution , Models, Molecular , Molecular Sequence Data , Mutation , Protein Denaturation , Sequence Alignment , Surface PropertiesABSTRACT
We have developed a novel DNA microarray-based approach for identification of the sequence-specificity of single-stranded nucleic-acid-binding proteins (SNABPs). For verification, we have shown that the major cold shock protein (CspB) from Bacillus subtilis binds with high affinity to pyrimidine-rich sequences, with a binding preference for the consensus sequence, 5'-GTCTTTG/T-3'. The sequence was modelled onto the known structure of CspB and a cytosine-binding pocket was identified, which explains the strong preference for a cytosine base at position 3. This microarray method offers a rapid high-throughput approach for determining the specificity and strength of ss DNA-protein interactions. Further screening of this newly emerging family of transcription factors will help provide an insight into their cellular function.