Expression of Atoh1, Gfi1, and Pou4f3 in the mature cochlea reprograms nonsensory cells into hair cells

Mechanosensory hair cells of the mature mammalian organ of Corti do not regenerate; consequently, loss of hair cells leads to permanent hearing loss. Although nonmammalian vertebrates can regenerate hair cells from neighboring supporting cells, many humans with severe hearing loss lack both hair cells and supporting cells, with the organ of Corti being replaced by a flat epithelium of nonsensory cells. To determine whether the mature cochlea can produce hair cells in vivo, we reprogrammed nonsensory cells adjacent to the organ of Corti with three hair cell transcription factors: Gfi1, Atoh1, and Pou4f3. We generated numerous hair cell-like cells in nonsensory regions of the cochlea and new hair cells continued to be added over a period of 9 wk. Significantly, cells adjacent to reprogrammed hair cells expressed markers of supporting cells, suggesting that transcription factor reprogramming of nonsensory cochlear cells in adult animals can generate mosaics of sensory cells like those seen in the organ of Corti. Generating such sensory mosaics by reprogramming may represent a potential strategy for hearing restoration in humans.

Characterization of Cancer Evolution Landscape Based on Accurate Detection of Somatic Mutations in Single Tumor Cells

Accurate detection of somatic mutations in single tumor cells is greatly desired as it allows us to quantify the single-cell mutation burden and construct the mutation-based phylogenetic tree. Here we developed scNanoSeq chemistry and profiled 842 single cells from 21 human breast cancer samples. The majority of the mutation-based phylogenetic trees comprise a characteristic stem evolution followed by the clonal sweep. We observed the subtype-dependent lengths in the stem evolution. To explain this phenomenon, we propose that the differences are related to different reprogramming required for different subtypes of breast cancer. Furthermore, we reason that the time that the tumor-initiating cell took to acquire the critical clonal-sweep-initiating mutation by random chance set the time limit for the reprogramming process. We refer to this model as a reprogramming and critical mutation co-timing (RCMC) subtype model. Next, in the sweeping clone, we observed that tumor cells undergo a branched evolution with rapidly decreasing selection. In the most recent clades, effectively neutral evolution has been reached, resulting in a substantially large number of mutational heterogeneities. Integrative analysis with 522-713X ultra-deep bulk whole genome sequencing (WGS) further validated this evolution mode. Mutation-based phylogenetic trees also allow us to identify the early branched cells in a few samples, whose phylogenetic trees support the gradual evolution of copy number variations (CNVs). Overall, the development of scNanoSeq allows us to unveil novel insights into breast cancer evolution.

Droplet-based transcriptome profiling of individual synapses

Synapses are crucial structures that mediate signal transmission between neurons in complex neural circuits and display considerable morphological and electrophysiological heterogeneity. So far we still lack a high-throughput method to profile the molecular heterogeneity among individual synapses. In the present study, we develop a droplet-based single-cell (sc) total-RNA-sequencing platform, called Multiple-Annealing-and-Tailing-based Quantitative scRNA-seq in Droplets, for transcriptome profiling of individual neurites, primarily composed of synaptosomes. In the synaptosome transcriptome, or 'synaptome', profiling of both mouse and human brain samples, we detect subclusters among synaptosomes that are associated with neuronal subtypes and characterize the landscape of transcript splicing that occurs within synapses. We extend synaptome profiling to synaptopathy in an Alzheimer's disease (AD) mouse model and discover AD-associated synaptic gene expression changes that cannot be detected by single-nucleus transcriptome profiling. Overall, our results show that this platform provides a high-throughput, single-synaptosome transcriptome profiling tool that will facilitate future discoveries in neuroscience.

[The correlation between plasma heat shock proteins 90α levels and white matter hyperintensity in patients with cerebral small vessel disease]

Objective: To investigate the relationship between plasma heat shock proteins 90α(Hsp90α) levels and the white matter hyperintensity(WMH) in patients with cerebral small vessel disease(SVD). Methods: Patients admitted to the Department of Neurology, the First Affiliated Hospital of Zhengzhou University from March to August 2021 and diagnosed with WMH by magnetic resonance examination (MRI) were selected as the case group, matched with physical examination patients who visited the Department of Medical Examination during the same period and showed no WMH on MRI and no history of neurological diseases as the control group, and the level of plasma Hsp90α was quantitatively detected by enzyme-linked immunosorbent assay. Mann-Whitney U test was used to compare whether there was a difference in plasma Hsp90α levels between the control group and the case group.Multivariate logistic regression analysis was used to explore the related factors of WMH in patients with SVD. Results: Of the 183 subjects, the control group (n=73) consisted of 28 males and 45 females, aged (54±10) years, while the case group (n=110) consisted of 71 males and 39 females, aged (64±10) years old. Plasma Hsp90α level was higher in the case group than that of the control group [53.33(35.33, 70.09) ng/ml vs 35.02(18.51, 54.95) ng/ml, P<0.001]. After adjusting for confounding factors by multivariate analysis, the results showed that plasma Hsp90α levels greater than 58.34 ng/ml was associated with WMH (P=0.002, OR=5.931, 95%CI:1.955-17.995). Conclusion: Higher level of plasma Hsp90α is associated with WMH in patients with SVD.

Single-cell Damagenome Profiling by Linear Copying and Splitting based Whole Genome Amplification (LCS-WGA)

Spontaneous DNA damage frequently occurs on the human genome, and it could alter gene expression by inducing mutagenesis or epigenetic changes. Therefore, it is highly desired to profile DNA damage distribution on the human genome and identify the genes that are prone to DNA damage. Here, we present a novel single-cell whole-genome amplification method which employs linear-copying followed by a split-amplification scheme, to efficiently remove amplification errors and achieve accurate detection of DNA damage in individual cells. In comparison to previous methods that measure DNA damage, our method uses a next-generation sequencing platform to detect misincorporated bases derived from spontaneous DNA damage with single-cell resolution.

Single-cell damagenome profiling unveils vulnerable genes and functional pathways in human genome toward DNA damage

We report a novel single-cell whole-genome amplification method (LCS-WGA) that can efficiently capture spontaneous DNA damage existing in single cells. We refer to these damage-associated single-nucleotide variants as "damSNVs," and the whole-genome distribution of damSNVs as the damagenome. We observed that in single human neurons, the damagenome distribution was significantly correlated with three-dimensional genome structures. This nonuniform distribution indicates different degrees of DNA damage effects on different genes. Next, we identified the functionals that were significantly enriched in the high-damage genes. Similar functionals were also enriched in the differentially expressed genes (DEGs) detected by single-cell transcriptome of both Alzheimer's disease (AD) and autism spectrum disorder (ASD). This result can be explained by the significant enrichment of high-damage genes in the DEGs of neurons for both AD and ASD. The discovery of high-damage genes sheds new lights on the important roles of DNA damage in human diseases and disorders.

Single-PanIN-seq unveils that ARID1A deficiency promotes pancreatic tumorigenesis by attenuating KRAS-induced senescence

ARID1A is one of the most frequently mutated epigenetic regulators in a wide spectrum of cancers. Recent studies have shown that ARID1A deficiency induces global changes in the epigenetic landscape of enhancers and promoters. These broad and complex effects make it challenging to identify the driving mechanisms of ARID1A deficiency in promoting cancer progression. Here, we identified the anti-senescence effect of Arid1a deficiency in the progression of pancreatic intraepithelial neoplasia (PanIN) by profiling the transcriptome of individual PanINs in a mouse model. In a human cell line model, we found that ARID1A deficiency upregulates the expression of aldehyde dehydrogenase 1 family member A1 (ALDH1A1), which plays an essential role in attenuating the senescence induced by oncogenic KRAS through scavenging reactive oxygen species. As a subunit of the SWI/SNF chromatin remodeling complex, our ATAC sequencing data showed that ARID1A deficiency increases the accessibility of the enhancer region of ALDH1A1. This study provides the first evidence that ARID1A deficiency promotes pancreatic tumorigenesis by attenuating KRAS-induced senescence through the upregulation of ALDH1A1 expression.

Contextual cues from cancer cells govern cancer-associated fibroblast heterogeneity

Cancer cells function as primary architects of the tumor microenvironment. However, the molecular features of cancer cells that govern stromal cell phenotypes remain unclear. Here, we show that cancer-associated fibroblast (CAF) heterogeneity is driven by lung adenocarcinoma (LUAD) cells at either end of the epithelial-to-mesenchymal transition (EMT) spectrum. LUAD cells that have high expression of the EMT-activating transcription factor ZEB1 reprogram CAFs through a ZEB1-dependent secretory program and direct CAFs to the tips of invasive projections through a ZEB1-driven CAF repulsion process. The EMT, in turn, sensitizes LUAD cells to pro-metastatic signals from CAFs. Thus, CAFs respond to contextual cues from LUAD cells to promote metastasis.

Bota-Rabassedas et al. show that EMT in lung adenocarcinoma cells activates a secretory process that governs CAF heterogeneity and, in turn, sensitizes lung adenocarcinoma cells to pro-metastatic signals from CAFs. Thus, EMT positions lung adenocarcinoma cells at the apex of a signaling hierarchy in the tumor microenvironment.

A study on the protective effect of molecular hydrogen on osteoradionecrosis of the jaw in rats

The aim of this study was to investigate the protective effect of hydrogen in a rat model of osteoradionecrosis of the jaw (ORNJ). The rats and bone marrow-derived mesenchymal stem cells (BMSCs) were pre-treated with hydrogen before receiving irradiation (7Gy per fraction, five fractions in total once a day for rats, 4Gy for BMSCs). Reactive oxygen species (ROS) and cell differentiation were measured in the BMSCs. Also, the radioprotective effect of hydrogen for ORNJ in Sprague-Dawley rats was examined by gross clinical manifestations, micro-computed tomography, and histology. Hydrogen significantly reduced the production of ROS in BMSCs after irradiation. The cell viability was significantly decreased after irradiation (P= 0.001), but pre-treatment with hydrogen before irradiation increased the cell viability (P= 0.025). Hydrogen considerably increased the cellular differentiation potential of the irradiated cells. Comparing with the rats underwent irradiaton only, those rats treated by hydrogen-rich saline significantly appeared improved occlusion, salivation, alopecia, oral ulcer, and less bone necrosis. Myofibroblasts accumulated overwhelmingly in the fibrosis medulla and around the sequestrum after irradiation, and this was decreased in the group pre-treated with hydrogen. Hydrogen may represent a strategy for the prevention and treatment of ORNJ. Its high efficacy and low toxicity suggest possible therapeutic application.

Fibroblast heterogeneity and its impact on extracellular matrix and immune landscape remodeling in cancer

Tumor progression is marked by dense collagenous matrix accumulations that dynamically reorganize to accommodate a growing and invasive tumor mass. Cancer-associated fibroblasts (CAFs) play an essential role in matrix remodeling and influence other processes in the tumor microenvironment, including angiogenesis, immunosuppression, and invasion. These findings have spawned efforts to elucidate CAF functionality at the single-cell level. Here, we will discuss how those efforts have impacted our understanding of the ways in which CAFs govern matrix remodeling and the influence of matrix remodeling on the development of an immunosuppressive tumor microenvironment.

Dissolvable Polyacrylamide Beads for High-Throughput Droplet DNA Barcoding

Droplet-based single cell sequencing technologies, such as inDrop, Drop-seq, and 10X Genomics, are catalyzing a revolution in the understanding of biology. Barcoding beads are key components for these technologies. What is limiting today are barcoding beads that are easy to fabricate, can efficiently deliver primers into drops, and thus achieve high detection efficiency. Here, this work reports an approach to fabricate dissolvable polyacrylamide beads, by crosslinking acrylamide with disulfide bridges that can be cleaved with dithiothreitol. The beads can be rapidly dissolved in drops and release DNA barcode primers. The dissolvable beads are easy to synthesize, and the primer cost for the beads is significantly lower than that for the previous barcoding beads. Furthermore, the dissolvable beads can be loaded into drops with >95% loading efficiency of a single bead per drop and the dissolution of beads does not influence reverse transcription or the polymerase chain reaction (PCR) in drops. Based on this approach, the dissolvable beads are used for single cell RNA and protein analysis.

Transcriptomic and epigenetic regulation of hair cell regeneration in the mouse utricle and its potentiation by Atoh1

The mammalian cochlea loses its ability to regenerate new hair cells prior to the onset of hearing. In contrast, the adult vestibular system can produce new hair cells in response to damage, or by reprogramming of supporting cells with the hair cell transcription factor Atoh1. We used RNA-seq and ATAC-seq to probe the transcriptional and epigenetic responses of utricle supporting cells to damage and Atoh1 transduction. We show that the regenerative response of the utricle correlates with a more accessible chromatin structure in utricle supporting cells compared to their cochlear counterparts. We also provide evidence that Atoh1 transduction of supporting cells is able to promote increased transcriptional accessibility of some hair cell genes. Our study offers a possible explanation for regenerative differences between sensory organs of the inner ear, but shows that additional factors to Atoh1 may be required for optimal reprogramming of hair cell fate.

Endocrine lineage biases arise in temporally distinct endocrine progenitors during pancreatic morphogenesis

Decoding the molecular composition of individual Ngn3 + endocrine progenitors (EPs) during pancreatic morphogenesis could provide insight into the mechanisms regulating hormonal cell fate. Here, we identify population markers and extensive cellular diversity including four EP subtypes reflecting EP maturation using high-resolution single-cell RNA-sequencing of the e14.5 and e16.5 mouse pancreas. While e14.5 and e16.5 EPs are constantly born and share select genes, these EPs are overall transcriptionally distinct concomitant with changes in the underlying epithelium. As a consequence, e16.5 EPs are not the same as e14.5 EPs: e16.5 EPs have a higher propensity to form beta cells. Analysis of e14.5 and e16.5 EP chromatin states reveals temporal shifts, with enrichment of beta cell motifs in accessible regions at later stages. Finally, we provide transcriptional maps outlining the route progenitors take as they make cell fate decisions, which can be applied to advance the in vitro generation of beta cells.

Endocrine progenitors form early in pancreatic development but the diversity of this cell population is unclear. Here, the authors use single cell RNA sequencing of the mouse pancreas at e14.5 and e16.5 to show that endocrine progenitors are temporally distinct and those formed later are more likely to become beta cells

Treatment of a Pediatric Case of Severe Hemorrhagic Cystitis: Case Report and Review of Literature

Hemorrhagic cystitis is one of the complications of allogeneic hematopoietic stem cell transplantation. Treatment of hemorrhagic cystitis is difficult, especially in pediatric patients. A pediatric case of severe hemorrhagic cystitis after hematopoietic stem cell transplantation was treated in our hospital with arterial embolization combined with corticosteroid therapy because the conventional therapy was invalid for him. After the treatment, hemorrhagic cystitis was cured. During follow-up, the patient was in stable condition, with normal urine, blood cells returned to normal, bone marrow was in complete remission state, and disease-free survival for more than 8 months. Selective bladder arterial embolism followed by corticosteroid therapy successfully treated the patient.

Mutual regulation of tumour vessel normalization and immunostimulatory reprogramming

Blockade of angiogenesis can retard tumour growth, but may also paradoxically increase metastasis^1,2. Vessel normalization (VN) may resolve this paradox³. VN involves increased pericyte coverage, improved tumour vessel perfusion, reduced vascular permeability, and consequently mitigated hypoxia³. While these processes alter tumour progression, their regulation is poorly understood. Here we show that Type 1 T helper (Th1) cells play a crucial role in VN. Bioinformatic analyses revealed that gene expression features related to VN correlate with immunostimulatory pathways, especially T lymphocyte (TL) infiltration/activities. To delineate the causal relationship, we employed various mouse models with VN or TL deficiencies. While VN disruption reduced TL infiltration as expected⁴, reciprocal depletion or inactivation of CD4⁺-TLs decreased VN, indicating a mutually-regulatory loop. Additionally, CD4⁺-TL activation by immune checkpoint blockade (ICB) increased VN. IFNγ⁺ Th1 cells are the major population associated with VN. Patient-derived xenograft (PDX) tumours growing in immunodeficient animal hosts exhibited enhanced hypoxia compared to the original tumours in immunocompetent human hosts, which was reduced by adoptive Th1 transfer. Our findings elucidate an unexpected role of Th1 in vasculature and immune reprogramming. Th1 cells may be a marker and a determinant of both ICB and anti-angiogenesis efficacies.

Protein Structure Prediction: The Next Generation

Over the last 10-15 years a general understanding of the chemical reaction of protein folding has emerged from statistical mechanics. The lessons learned from protein folding kinetics based on energy landscape ideas have benefited protein structure prediction, in particular the development of coarse grained models. We survey results from blind structure prediction. We explore how second generation prediction energy functions can be developed by introducing information from an ensemble of previously simulated structures. This procedure relies on the assumption of a funneled energy landscape keeping with the principle of minimal frustration. First generation simulated structures provide an improved input for associative memory energy functions in comparison to the experimental protein structures chosen on the basis of sequence alignment.

Probing meiotic recombination and aneuploidy of single sperm cells by whole-genome sequencing

Meiotic recombination creates genetic diversity and ensures segregation of homologous chromosomes. Previous population analyses yielded results averaged among individuals and affected by evolutionary pressures. We sequenced 99 sperm from an Asian male by using the newly developed amplification method-multiple annealing and looping-based amplification cycles-to phase the personal genome and map recombination events at high resolution, which are nonuniformly distributed across the genome in the absence of selection pressure. The paucity of recombination near transcription start sites observed in individual sperm indicates that such a phenomenon is intrinsic to the molecular mechanism of meiosis. Interestingly, a decreased crossover frequency combined with an increase of autosomal aneuploidy is observable on a global per-sperm basis.

Genome-wide detection of single-nucleotide and copy-number variations of a single human cell

Kindred cells can have different genomes because of dynamic changes in DNA. Single-cell sequencing is needed to characterize these genomic differences but has been hindered by whole-genome amplification bias, resulting in low genome coverage. Here, we report on a new amplification method-multiple annealing and looping-based amplification cycles (MALBAC)-that offers high uniformity across the genome. Sequencing MALBAC-amplified DNA achieves 93% genome coverage ≥1x for a single human cell at 25x mean sequencing depth. We detected digitized copy-number variations (CNVs) of a single cancer cell. By sequencing three kindred cells, we were able to identify individual single-nucleotide variations (SNVs), with no false positives detected. We directly measured the genome-wide mutation rate of a cancer cell line and found that purine-pyrimidine exchanges occurred unusually frequently among the newly acquired SNVs.

General properties of transcriptional time series in Escherichia coli

Gene activity is described by the time-series of discrete, stochastic mRNA production events. This transcriptional time-series exhibits intermittent, bursty behavior. One consequence of this temporal intricacy is that gene expression can be tuned by varying different features of the time-series. What schemes for varying the transcriptional time-series are observed in the cell? Are the observed properties of these time-series optimized for cellular function? To address these questions, we characterize mRNA copy-number statistics at single-molecule resolution from multiple Escherichia coli promoters. We find that the degree of burstiness depends only on the gene expression level, while being independent of the details of gene regulation. The observed behavior is explained by the underlying variation in the duration of bursting events. Using information theory, we find that the properties of the transcriptional time series allow the cell to efficiently map the extracellular concentration of inducer molecules to intracellular levels of mRNA and proteins.

Lysogen stability is determined by the frequency of activity bursts from the fate-determining gene

Bacterial lysogeny serves as a simple paradigm for cell differentiation.

We characterize the activity of the fate-determining genes, cI and cro, with single-molecule resolution.

Stability of the lysogenic state is found to depend in a simple manner on the frequency of activity bursts from cI.

The ability of living cells to maintain an inheritable memory of their gene-expression state is key to cellular differentiation. Bacterial lysogeny serves as a simple paradigm for long-term cellular memory. In this study, we address the following question: in the absence of external perturbation, how long will a cell stay in the lysogenic state before spontaneously switching away from that state? We show by direct measurement that lysogen stability exhibits a simple exponential dependence on the frequency of activity bursts from the fate-determining gene, cI. We quantify these gene-activity bursts using single-molecule-resolution mRNA measurements in individual cells, analyzed using a stochastic mathematical model of the gene-network kinetics. The quantitative relation between stability and gene activity is independent of the fine details of gene regulation, suggesting that a quantitative prediction of cell-state stability may also be possible in more complex systems.

Protein structure prediction: do hydrogen bonding and water-mediated interactions suffice?

The many-body physics of hydrogen bond formation in alpha-helices of globular proteins was investigated using a simple physics-based model. Specifically, a context-sensitive hydrogen bond potential, which depends on residue identity and degree of solvent exposure, was used in the framework of the Associated Memory Hamiltonian codes developed previously but without using local-sequence structure matches ("memories"). Molecular dynamics simulations employing the energy function using the context-sensitive hydrogen bond potential alone (the "amnesiac" model) were used to generate low energy structures for three alpha-helical test proteins. The resulting structures were compared to both the X-ray crystal structures of the test proteins and the results obtained using the full Associated Memory Hamiltonian previously used. Results show that the amnesiac Hamiltonian was able to generate structures with reasonably high structural similarity (Q approximately 0.4) to that of the native protein but only with the use of predicted secondary structure information encoding local steric signals. Low energy structures obtained using the amnesiac Hamiltonian without any a priori secondary structure information had considerably less similarity to the native protein structures (Q approximately 0.3). Both sets of results utilizing the amnesiac Hamiltonian are poorer than when local-sequence structure matches are used.

Variationally determined free energy profiles for structural models of proteins: characteristic temperatures for folding and trapping

Characterizing the phase diagram for proteins is important both for laboratory studies and for the development of structure prediction algorithms. Using a variational scheme, we calculated the generic features of the protein thermostability over a large range of temperatures for a set of more than 50 different proteins using a model based on native structure alone. Focusing on a specific system, protein G, we further examined, using a more realistic model that includes the nonnative interaction, the thermostability of both the native state and a collection of trap structures. By surveying the native structures for many proteins and by paying closer attention to the various trap structures of protein G, we obtained an overall understanding of the folding dynamics far from the conditions usually focused on; namely, those near the folding temperature alone. Two characteristic temperatures (shown to scale with folding temperature in general) signal drastic changes in the folding mechanism. The variational calculations suggest that most proteins would, indeed, fold in a barrierless manner below a critical temperature analogous to a spinodal in crystallization. For fixed interaction strengths, this temperature, however, seems to be generally very low, approximately 50% of the equilibrium folding temperature. Likewise, native proteins, in general, would unfold in a completely barrierless way at a temperature 25% above folding temperature according to these variational calculations. We also studied the distribution of free energy profiles for escape from a set of trap structures generated by simulations.

Establishing the entatic state in folding metallated Pseudomonas aeruginosa azurin

Understanding how the folding of proteins establishes their functional characteristics at the molecular level challenges both theorists and experimentalists. The simplest test beds for confronting this issue are provided by electron transfer proteins. The environment provided by the folded protein to the cofactor tunes the metal's electron transport capabilities as envisioned in the entatic hypothesis. To see how the entatic state is achieved one must study how the folding landscape affects and in turn is affected by the metal. Here, we develop a coarse-grained functional to explicitly model how the coordination of the metal (which results in a so-called entatic or rack-induced state) modifies the folding of the metallated Pseudomonas aeruginosa azurin. Our free-energy functional-based approach directly yields the proper nonlinear extra-thermodynamic free energy relationships for the kinetics of folding the wild type and several point-mutated variants of the metallated protein. The results agree quite well with corresponding laboratory experiments. Moreover, our modified free-energy functional provides a sufficient level of detail to explicitly model how the geometric entatic state of the metal modifies the dynamic folding nucleus of azurin.

Statistics of cellular signal transduction as a race to the nucleus by multiple random walkers in compartment/phosphorylation space

Cellular signal transduction often involves a reaction network of phosphorylation and transport events arranged with a ladder topology. If we keep track of the location of the phosphate groups describing an abstract state space, a simple model of signal transduction involving enzymes can be mapped on to a problem of how multiple biased random walkers compete to reach their target in the nucleus yielding a signal. Here, the first passage time probability and the survival probability for multiple walkers can be used to characterize the response of the network. The statistics of the first passage through the network has an asymmetric distribution with a long tail arising from the hierarchical structure of the network. This distribution implies a significant difference between the mean and the most probable signal transduction time. The response patterns for various external inputs generated by our model agree with recent experiments. In addition, the model predicts that there is an optimal phosphorylation enzyme concentration for rapid signal transduction.

Phi-value analysis of apo-azurin folding: comparison between experiment and theory

Pseudomonas aeruginosa azurin is a 128-residue beta-sandwich metalloprotein; in vitro kinetic experiments have shown that it folds in a two-state reaction. Here, we used a variational free energy functional to calculate the characteristics of the transition state ensemble (TSE) for folding of the apo-form of P. aeruginosa azurin and investigate how it responds to thermal and mutational changes. The variational method directly yields predicted chevron plots for wild-type and mutant apo-forms of azurin. In parallel, we performed in vitro kinetic-folding experiments on the same set of azurin variants using chemical perturbation. Like the wild-type protein, all apo-variants fold in apparent two-state reactions both in calculations and in stopped-flow mixing experiments. Comparisons of phi (phi) values determined from the experimental and theoretical chevron parameters reveal an excellent agreement for most positions, indicating a polarized, highly structured TSE for folding of P. aeruginosa apo-azurin. We also demonstrate that careful analysis of side-chain interactions is necessary for appropriate theoretical description of core mutants.

Nonequilibrium self-assembly of linear fibers: microscopic treatment of growth, decay, catastrophe and rescue

Many of the large structures of cells are constructed from fibers. These fibers self-assemble from individual proteins in a far-from-equilibrium fashion. Nonequilibrium self-assembly results in a highly dynamic process at the subcellular level that can be regulated and tuned to carry out many of the biological functions of the cell: growth, division and locomotion. We construct and analyze a nonequilibrium model of the dynamic end of a biological fiber that possesses site-resolved resolution. We solve for the steady states of this nonequilibrium system using a variational method. The results are compared to exact numerical solutions for systems with modest size. Using an effective reaction coordinate, we construct an effective potential from the steady-state distribution. The stochastic transitions of the system can be analyzed in this representation. We then apply this method to model microtubule systems. Predictions for macroscopic catastrophe, rescue and dynamic instability in the steady states are analyzed. We find that the length of the cap of the microtubule is small. The relations between the catastrophe/rescue rate and the growth rate are also discussed.

Role of Topology, Nonadditivity, and Water-Mediated Interactions in Predicting the Structures of α/β Proteins

The folding of alpha/beta proteins involves most of the commonly known structural and dynamic complexities of the protein energy landscapes. Thus, the interplay among different structural components, taking into account the cooperative interactions, is important in determining the success of protein structure prediction. In this work we present further developments of our knowledge-based force field for alpha/beta proteins, introducing more realistic modeling of many-body interactions governing the folding of beta-sheets. The model's innovations highlight both specific topological characteristics of secondary structures and the generic nonadditive interactions that are mediated by water. We also investigate how a coarse biasing of the protein morphology can be used to understand the role of heterogeneity in protein collapse. Analysis of the simulation results for three test alpha/beta proteins indicates that the addition of the topological and many-body ingredients to the model helps to greatly reduce the roughness in the energy landscape. Consequently, high quality candidate structures for alpha/beta proteins can be generated from simulated annealing runs, using very modest amounts of computer time.

The folding energy landscape and phosphorylation: modeling the conformational switch of the NFAT regulatory domain

An energy landscape approach predicts the conformational changes of the configurations of the regulatory domain of the protein nuclear factor of activated T cells (NFAT) caused by phosphorylation of specific multiple sites. Structurally local effects and secondary structural changes are modeled using all-atom Brownian dynamics to investigate the changes of the backbone torsional distributions upon phosphorylation. For tertiary and global changes, we employ a coarse-grained model to sample ensembles of conformations both with and without phosphorylation. At the secondary structure level, phosphorylation moderately increases the helical propensity and gives a more rigid local backbone conformation. The tertiary effects of phosphorylation caused by the extensive charge modification are more pronounced and collectively change the conformation of the regulatory domain of NFAT from a flexible globular ensemble to a rather rigid helical bundle, blocking access to the nuclear localization sequence. These studies give computational support to one scenario conjectured from experiments.

A funneled energy landscape for cytochrome c directly predicts the sequential folding route inferred from hydrogen exchange experiments

Proteins fold through a variety of mechanisms. For a given protein, folding routes largely depend on the protein's stability and its native-state geometry, because the landscape is funneled. These ideas are corroborated for cytochrome c by using a coarse-grained topology-based model with a perfect funnel landscape that includes explicit modeling of the heme. The results show the importance of the heme as a nucleation site and explain the observed hydrogen exchange patterns of cytochrome c within the context of energy landscape theory.

Biological effects of c-Mer receptor tyrosine kinase in hematopoietic cells depend on the Grb2 binding site in the receptor and activation of NF-kappaB

The c-Mer receptor tyrosine kinase (RTK) is most closely related to chicken c-Eyk and belongs to the Axl RTK subfamily. Although not detected in normal lymphocytes, c-Mer is expressed in B- and T-cell leukemia cell lines, suggesting an association with lymphoid malignancies. To gain an understanding of the role of this receptor in lymphoid cells, we expressed in murine interleukin-3 (IL-3)-dependent Ba/F3 pro-B-lymphocyte cells a constitutively active receptor, CDMer, formed from the CD8 extracellular domain and the c-Mer intracellular domain. Cells transfected with a plasmid encoding the CDMer receptor became IL-3 independent. When tyrosine (Y)-to-phenylalanine (F) mutations were introduced into c-Mer, only the Y867 change significantly reduced the IL-3-independent cell proliferation. The Y867 residue in the CDMer receptor mediated the binding of Grb2, which recruited the p85 phosphatidylinositol 3-kinase (PI 3-kinase). Despite the difference in promotion of proliferation, both the CDMer and mutant F867 receptors activated Erk in transfected cells. On the other hand, we found that both transcriptional activation of NF-kappaB and activation of PI 3-kinase were significantly suppressed with the F867 mutant receptor, suggesting that the activation of antiapoptotic pathways is the major mechanism for the observed phenotypic difference. Consistent with this notion, apoptosis induced by IL-3 withdrawal was strongly prevented by CDMer but not by the F867 mutant receptor.

Identification of the product of growth arrest-specific gene 6 as a common ligand for Axl, Sky, and Mer receptor tyrosine kinases

Axl, Sky, and Mer, members of an Axl/Sky receptor tyrosine kinase subfamily, are typified by the cell adhesion molecule-related extracellular domain. The product of growth arrest-specific gene 6 (Gas6), structurally homologous to the anticoagulant protein S, was recently identified as the ligand for Axl and Sky, but the ligand for Mer remained unknown. We have now obtained evidence that Gas6 can also function as a ligand for Mer. Co-precipitation analysis, using soluble receptors of Axl, Sky, and Mer (Axl-Fc, Sky-Fc, and Mer-Fc) composed of the extracellular domain of receptors fused to the Fc domain of immunoglobulin G1, clearly showed that Gas6, but not protein S, specifically bound to Axl-Fc, Sky-Fc, and Mer-Fc fusion proteins. Quantitative kinetic analyses using a BIAcore biosensor instrument revealed dissociation constants (Kd) of the binding of rat Gas6 to Axl-Fc, Sky-Fc, and Mer-Fc are 0.4, 2.7, and 29 nM, respectively. We also found that Gas6 stimulated tyrosine phosphorylation of Axl, Sky, and Mer receptors ectopically expressed in Chinese hamster ovary cells. Taken together, these findings suggest that Gas6 is a common ligand for Axl, Sky, and Mer, all known members of an Axl/Sky receptor subfamily.

Unique signal transduction of Eyk: constitutive stimulation of the JAK-STAT pathway by an oncogenic receptor-type tyrosine kinase

The proto-oncogene c-eyk, the cellular counterpart of a transforming oncogene, v-eyk, encodes a receptor protein tyrosine kinase with a distinctive extracellular region. We now demonstrate that c-Eyk can be constitutively activated through dimerization, and that the active Eyk displays a unique signaling pattern. When the kinase domain of c-Eyk was fused to the extracellular and transmembrane domains of CD8, the resulting chimera showed elevated kinase activity and caused cellular transformation. We found that the activated Eyk kinases, both v- and c-Eyk, constitutively stimulate the JAK-STAT pathway, while exerting little effect on other signaling routes such as the Ras-MAP kinase and the JNK pathways. The activated Eyk kinases specifically stimulate tyrosine phosphorylation of STAT1, STAT3 and JAK1. These downstream molecules also co-immunoprecipitate with the constitutively dimerized form of Eyk. The Eyk kinase activity is required for STAT1 stimulation. We found that the activation of STAT1 but not STAT3 correlates well with cellular transformation. In constitutively stimulating the JAK-STAT pathway, particularly STAT1, Eyk is unique in its downstream signaling and may be dependent on this pathway for cellular transformation.

Gene expression of human sperm component related to A4 amyloid precursor protein

Monoclonal antibodies were raised against a specific human sperm protein and designated as the YWK-II mAb. The partial cDNA encoding the protein was isolated from a rat testis lambda gt11 expression library and the amino acid sequence of the protein was deduced. The cytoplasmic-transmembrane domains of the deduced protein had high homology with the A4 amyloid precursor protein of Alzheimer's disease. To evaluate the stage of spermatogenesis when the gene was expressed, single-stranded 35S-labeled RNA probes were prepared from the cDNA. By an in situ hybridization technique the mRNA for the antigen was detected in germ cells at all stages of spermatogenesis. The finding that the gene is expressed in spermatogonia suggests possible involvement in the initiation of germ cell differentiation or in the detachment of spermatogonia from the basement membrane.

Expression of a sperm protein gene during spermatogenesis in mammalian testis: an in situ hybridization study

In previous work a specific membrane protein with an estimated Mr of 20.1 kDa was purified from rabbit sperm tails and designated as rSMP-B protein. Antibodies were raised against rSMP-B protein and used to isolate and identify the cDNA coding the rSMP-B protein from a rat testis lambda gt11 expression library. The nucleotide sequence of the cDNA was determined in a previous study. Single-stranded 35S-labeled RNA probes were prepared. With the techniques of in situ hybridization, rSMP-B mRNA was detected in spermatids of rat and rabbit testis. The present results support our previous observation that immunization of male rabbits with the rSMP-B protein results in the arrest of spermatogenesis at the spermatid stage. Overall, rSMP-B protein appears to be involved in spermiogenesis, and the synthesis of the mRNA encoding the protein occurs in germ cells during the postmeiotic haploid phase of spermatogenesis.