TsImpute: an accurate two-step imputation method for single-cell RNA-seq data.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) technology has enabled discovering gene expression patterns at single cell resolution. However, due to technical limitations, there are usually excessive zeros, called "dropouts," in scRNA-seq data, which may mislead the downstream analysis. Therefore, it is crucial to impute these dropouts to recover the biological information. RESULTS: We propose a two-step imputation method called tsImpute to impute scRNA-seq data. At the first step, tsImpute adopts zero-inflated negative binomial distribution to discriminate dropouts from true zeros and performs initial imputation by calculating the expected expression level. At the second step, it conducts clustering with this modified expression matrix, based on which the final distance weighted imputation is performed. Numerical results based on both simulated and real data show that tsImpute achieves favorable performance in terms of gene expression recovery, cell clustering, and differential expression analysis. AVAILABILITY AND IMPLEMENTATION: The R package of tsImpute is available at https://github.com/ZhengWeihuaYNU/tsImpute.

Are children with IgA nephropathy different from adult patients?

BACKGROUND: Previously, several studies have indicated that pediatric IgA nephropathy (IgAN) might be different from adult IgAN, and treatment strategies might be also different between pediatric IgAN and adult IgAN. METHODS: We analyzed two prospective cohorts established by pediatric and adult nephrologists, respectively. A comprehensive analysis was performed investigating the difference in clinical and pathological characteristics, treatment, and prognosis between children and adults with IgAN. RESULTS: A total of 1015 children and 1911 adults with IgAN were eligible for analysis. More frequent gross hematuria (88% vs. 20%, p < 0.0001) and higher proteinuria (1.8 vs. 1.3 g/d, p < 0.0001) were seen in children compared to adults. In comparison, the estimated glomerular filtration rate (eGFR) was lower in adults (80.4 vs. 163 ml/min/1.73 m2, p < 0.0001). Hypertension was more prevalent in adult patients. Pathologically, a higher proportion of M1 was revealed (62% vs. 39%, p < 0.0001) in children than in adults. S1 (62% vs. 28%, p < 0.0001) and T1-2 (34% vs. 8%, p < 0.0001) were more frequent in adults. Adjusted by proteinuria, eGFR, and hypertension, children were more likely to be treated with glucocorticoids than adults (87% vs. 45%, p < 0.0001). After propensity score matching, in IgAN with proteinuria > 1 g/d, children treated with steroids were 1.87 (95% CI 1.16-3.02, p = 0.01) times more likely to reach complete remission of proteinuria compared with adults treated with steroids. CONCLUSIONS: Children present significantly differently from adults with IgAN in clinical and pathological manifestations and disease progression. Steroid response might be better in children.

Working Condition Recognition Based on Transfer Learning and Attention Mechanism for a Rotary Kiln.

It is difficult to identify the working conditions of the rotary kilns due to the harsh environment in the kilns. The flame images of the firing zone in the kilns contain a lot of working condition information, but the flame image data sample size is too small to be used to fully extract the key features. In order to solve this problem, a method combining transfer learning and attention mechanism is proposed to extract key features of flame images, in which the deep residual network is used as the backbone network, the coordinate attention module is introduced to capture the position information and channel information on the branch of feature graphs, and the features of flame images obtained are further screened to improve the extraction ability. At the same time, migration learning is performed by the pre-trained ImageNet data set, and feature migration and parameter sharing are realized to cope with the training difficulty of a small sample data size. Moreover, an activation function Mish is introduced to reduce the loss of effective information. The experimental results show that, compared with traditional methods, the working condition recognition accuracy of rotary kilns is improved by about 5% with the proposed method.

LINC00346 promotes hepatocellular carcinoma progression via activating the JAK-STAT3 signaling pathway.

Hepatocellular carcinoma (HCC) remains the most common malignant tumor worldwide. Long noncoding RNAs can modulate various tumorigenic processes. In addition, growing evidence has indicated tha the Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathway is activated in multiple cancers, including HCC. Recently, it was found that LINC00346 can participate in several cancers. Nevertheless, the biological roles of LINC00346 in HCC have been barely investigated. In this study, the function of LINC00346 was specifically concentrated upon. We observed that LINC00346 was obviously elevated in HCC cells (Bel7404, Huh-6, HepG2, and QGY-7703 cells). Then, Bel7404 and HepG2 cells were overexpressed with LINC00346. Overexpression of LINC00346 repressed HCC cell survival and cell proliferation. In addition, apoptosis of Bel7404 and HepG2 cells was triggered by LINC00346 upregulation. Bel7404 and HepG2 cell cycle was arrested in the G1 phase by LINC00346. Meanwhile, we conducted wound-healing assay and Transwell invasion assays. As shown, we observed that the migratory and invasive capacities of Bel7404 and HepG2 cells were remarkably restrained by the increase of LINC00346. Moreover, we showed that LINC00346 overexpression activated the JAK-STAT3 pathway, which is involved in many cancers. Afterward, in vivo experiments were utilized and we proved that LINC00346 was able to induce HCC tumor growth via activating the JAK-STAT3 pathway. To conclude, we revealed the potential possibility of developing LINC00346 as an indicator for HCC.

Design of a free-form surface microlens array optical system with high efficiency and uniformity.

The microlens array has been widely applied in LED lighting source due to its special optical properties, but most of the research lacks the analysis and optimization of the complete mathematical models. Hence, the new design method of a free-form surface microlens array optical system is proposed in this paper. Based on the characteristics of TIR and the law of refraction, a complete mathematical model of the free-form microlens is established. By numerically solving a set of differential equations, the profile of the free-form surface microlens is obtained. Then we rotate the profile to get the free-form surface microlens. Finally, the proposed microlens array is simulated and analyzed in near-field and far-field situations, respectively. We also discuss the influence of microlens array characteristics on illumination performance. The result shows the uniformity and efficiency have been improved, both of which can reach more than 90%.

Multiple Binding Configurations of Fis Protein Pairs on DNA: Facilitated Dissociation versus Cooperative Dissociation.

As a master transcription regulator, the Fis protein influences over two hundred genes of E. coli. The Fis protein's nonspecific binding to DNA is widely acknowledged, and its kinetics of dissociation from DNA is strongly influenced by its surroundings: the dissociation rate increases as the concentration of the Fis protein in the solution phase increases. In this study, we use computational methods to explore the global binding energy landscape of the Fis1:Fis2:DNA ternary complex. The complex contains a binary-Fis molecular dyad whose formation relies on complex structural rearrangements. The simulations allow us to distinguish several different pathways for the dissociation of the protein from DNA with different functional outcomes and involving different protein stoichiometries: (1) simple exchange of proteins and (2) cooperative unbinding of two Fis proteins to yield bare DNA. In the case of exchange, the protein on the DNA is replaced by the solution-phase protein through competition for DNA binding sites. This process seen in fluorescence imaging experiments has been called facilitated dissociation. In the latter case of cooperative unbinding of pairs, two neighboring Fis proteins on DNA form a unique binary-Fis configuration via protein-protein interactions, which in turn leads to the codissociation of both molecules simultaneously, a process akin to the "molecular stripping" seen in the NFκB/IκB genetic broadcasting system. This simulation shows that the existence of multiple binding configurations of transcription factors can have a significant impact on the kinetics and outcome of transcription factor dissociation from DNA, with important implications for the systems biology of gene regulation by Fis.

Energy landscapes of a mechanical prion and their implications for the molecular mechanism of long-term memory.

Aplysia cytoplasmic polyadenylation element binding (CPEB) protein, a translational regulator that recruits mRNAs and facilitates translation, has been shown to be a key component in the formation of long-term memory. Experimental data show that CPEB exists in at least a low-molecular weight coiled-coil oligomeric form and an amyloid fiber form involving the Q-rich domain (CPEB-Q). Using a coarse-grained energy landscape model, we predict the structures of the low-molecular weight oligomeric form and the dynamics of their transitions to the ß-form. Up to the decamer, the oligomeric structures are predicted to be coiled coils. Free energy profiles confirm that the coiled coil is the most stable form for dimers and trimers. The structural transition from α to ß is shown to be concentration dependent, with the transition barrier decreasing with increased concentration. We observe that a mechanical pulling force can facilitate the α-helix to ß-sheet (α-to-ß) transition by lowering the free energy barrier between the two forms. Interactome analysis of the CPEB protein suggests that its interactions with the cytoskeleton could provide the necessary mechanical force. We propose that, by exerting mechanical forces on CPEB oligomers, an active cytoskeleton can facilitate fiber formation. This mechanical catalysis makes possible a positive feedback loop that would help localize the formation of CPEB fibers to active synapse areas and mark those synapses for forming a long-term memory after the prion form is established. The functional role of the CPEB helical oligomers in this mechanism carries with it implications for targeting such species in neurodegenerative diseases.

Molecular stripping in the NF-κB/IκB/DNA genetic regulatory network.

Genetic switches based on the [Formula: see text] system are master regulators of an array of cellular responses. Recent kinetic experiments have shown that [Formula: see text] can actively remove NF-κB bound to its genetic sites via a process called "molecular stripping." This allows the [Formula: see text] switch to function under kinetic control rather than the thermodynamic control contemplated in the traditional models of gene switches. Using molecular dynamics simulations of coarse-grained predictive energy landscape models for the constituent proteins by themselves and interacting with the DNA we explore the functional motions of the transcription factor [Formula: see text] and its various binary and ternary complexes with DNA and the inhibitor IκB. These studies show that the function of the [Formula: see text] genetic switch is realized via an allosteric mechanism. Molecular stripping occurs through the activation of a domain twist mode by the binding of [Formula: see text] that occurs through conformational selection. Free energy calculations for DNA binding show that the binding of [Formula: see text] not only results in a significant decrease of the affinity of the transcription factor for the DNA but also kinetically speeds DNA release. Projections of the free energy onto various reaction coordinates reveal the structural details of the stripping pathways.

Exploring the aggregation free energy landscape of the amyloid-ß protein (1-40).

A predictive coarse-grained protein force field [associative memory, water-mediated, structure, and energy model for molecular dynamics (AWSEM)-MD] is used to study the energy landscapes and relative stabilities of amyloid-ß protein (1-40) in the monomer and all of its oligomeric forms up to an octamer. We find that an isolated monomer is mainly disordered with a short α-helix formed at the central hydrophobic core region (L17-D23). A less stable hairpin structure, however, becomes increasingly more stable in oligomers, where hydrogen bonds can form between neighboring monomers. We explore the structure and stability of both prefibrillar oligomers that consist of mainly antiparallel ß-sheets and fibrillar oligomers with only parallel ß-sheets. Prefibrillar oligomers are polymorphic but typically take on a cylindrin-like shape composed of mostly antiparallel ß-strands. At the concentration of the simulation, the aggregation free energy landscape is nearly downhill. We use umbrella sampling along a structural progress coordinate for interconversion between prefibrillar and fibrillar forms to identify a conversion pathway between these forms. The fibrillar oligomer only becomes favored over its prefibrillar counterpart in the pentamer where an interconversion bottleneck appears. The structural characterization of the pathway along with statistical mechanical perturbation theory allow us to evaluate the effects of concentration on the free energy landscape of aggregation as well as the effects of the Dutch and Arctic mutations associated with early onset of Alzheimer's disease.

Comparing the Aggregation Free Energy Landscapes of Amyloid Beta(1-42) and Amyloid Beta(1-40).

Using a predictive coarse-grained protein force field, we compute and compare the free energy landscapes and relative stabilities of amyloid-ß protein (1-42) and amyloid-ß protein (1-40) in their monomeric and oligomeric forms up to the octamer. At the same concentration, the aggregation free energy profile of Aß42 is more downhill, with a computed solubility that is about 10 times smaller than that of Aß40. At a concentration of 40 µM, the clear free energy barrier between the pre-fibrillar tetramer form and the fibrillar pentamer in the Aß40 aggregation landscape disappears for Aß42, suggesting that the Aß42 tetramer has a more diverse structural range. To further compare the landscapes, we develop a cluster analysis based on the structural similarity between configurations and use it to construct an oligomerization map that captures the paths of easy interconversion between different but structurally similar states of oligomers for both species. A taxonomy of the oligomer species based on ß-sheet stacking topologies is proposed. The comparison of the two oligomerization maps highlights several key differences in the landscapes that can be attributed to the two additional C-terminal residues that Aß40 lacks. In general, the two terminal residues strongly stabilize the oligomeric structures for Aß42 relative to Aß40, and greatly facilitate the conversion from pre-fibrillar trimers to fibrillar tetramers.

Resolving the NFκB Heterodimer Binding Paradox: Strain and Frustration Guide the Binding of Dimeric Transcription Factors.

Many eukaryotic transcription factors function after forming oligomers. The choice of protein partners is a nonrandom event that has distinct functional consequences for gene regulation. In the present work we examine three dimers of transcription factors in the NFκB family: p50p50, p50p65, and p65p65. The NFκB dimers bind to a myriad of genomic sites and switch the targeted genes on or off with precision. The p65p50 heterodimer of NFκB is the strongest DNA binder, and its unbinding is controlled kinetically by molecular stripping from the DNA induced by IκB. In contrast, the homodimeric forms of NFκB, p50p50 and p65p65, bind DNA with significantly less affinity, which places the DNA residence of the homodimers under thermodynamic rather than kinetic control. It seems paradoxical that the heterodimer should bind more strongly than either of the symmetric homodimers since DNA is a nearly symmetric target. Using a variety of energy landscape analysis tools, here we uncover the features in the molecular architecture of NFκB dimers that are responsible for these drastically different binding free energies. We show that frustration in the heterodimer interface gives the heterodimer greater conformational plasticity, allowing the heterodimer to better accommodate the DNA. We also show how the elastic energy and mechanical strain in NFκB dimers can be found by extracting the principal components of the fluctuations in Cartesian coordinates as well as fluctuations in the space of physical contacts, which are sampled via simulations with a predictive energy landscape Hamiltonian. These energetic contributions determine the specific detailed mechanisms of binding and stripping for both homo- and heterodimers.

Binding of NFκB Appears to Twist the Ankyrin Repeat Domain of IκBα.

Total internal reflection fluorescence-based single-molecule Förster resonance energy transfer (FRET) measurements were previously carried out on the ankyrin repeat domain (ARD) of IκBα, the temporally regulated inhibitor of canonical NFκB signaling. Under native conditions, most of the IκBα molecules showed stable, high FRET signals consistent with distances between the fluorophores estimated from the crystal structures of the NFκB(RelA/p50)-IκBα complex. Similar high FRET efficiencies were found when the IκBα molecules were either free or in complex with NFκB(RelA/p50), and were interpreted as being consistent with the crystallographically observed ARD structure. An exception to this was observed when the donor and acceptor fluorophores were attached in AR3 (residue 166) and AR6 (residue 262). Surprisingly, the FRET efficiency was lower for the bound IκBα molecules (0.67) than for the free IκBα molecules (0.74), apparently indicating that binding of NFκB(RelA/p50) stretches the ARD of IκBα. Here, we conducted confocal-based single-molecule FRET studies to investigate this phenomenon in greater detail. The results not only recapitulated the apparent stretching of the ARD but also showed that the effect was more pronounced when the N-terminal domains (NTDs) of both RelA and p50 were present, even though the interface between NFκB(RelA/p50) and IκBα encompasses only the dimerization domains. We also performed mass spectrometry-detected amide hydrogen/deuterium exchange (HDXMS) experiments on IκBα as well as IκBα bound to dimerization-domain-only constructs or full-length NFκB(RelA/p50). Although we expected the stretched IκBα to have regions with increased exchange, instead the HDXMS experiments showed decreases in exchange in AR3 and AR6 that were more pronounced when the NFκB NTDs were present. Simulations of the interaction recapitulated the increased distance between residues 166 and 262, and also provide a plausible mechanism for a twisting of the IκBα ARD induced by interactions of the IκBα proline-glutamate-serine-threonine-rich sequence with positively charged residues in the RelA NTD.

Exploring the Free Energy Landscape of Nucleosomes.

The nucleosome is the fundamental unit for packaging the genome. A detailed molecular picture for its conformational dynamics is crucial for understanding transcription and gene regulation. We investigate the disassembly of single nucleosomes using a predictive coarse-grained protein DNA model with transferable force fields. This model quantitatively describes the thermodynamic stability of both the histone core complex and the nucleosome and predicts rates of transient nucleosome opening that match experimental measurements. Quantitative characterization of the free-energy landscapes reveals the mechanism of nucleosome unfolding in which DNA unwinding and histone protein disassembly are coupled. The interfaces between H2A-H2B dimers and the (H3-H4)2 tetramer are first lost when the nucleosome opens releasing a large fraction but not all of its bound DNA. For the short strands studied in single molecule experiments, the DNA unwinds asymmetrically from the histone proteins, with only one of its two ends preferentially exposed. The detailed molecular mechanism revealed in this work provides a structural basis for interpreting experimental studies of nucleosome unfolding.

Molecular Mechanism of Facilitated Dissociation of Fis Protein from DNA.

Fis protein is a nucleoid-associated protein that plays many roles in transcriptional regulation and DNA site-specific recombination. In contrast to the naïve expectation based on stoichiometry, recent single-molecule studies have shown that the dissociation of Fis protein from DNA is accelerated by increasing the concentration of the Fis protein. Because the detailed molecular mechanism of facilitated dissociation is still not clear, in this study, we employ computational methods to explore the binding landscapes of Fis:DNA complexes with various stoichiometries. When two Fis molecules are present, simulations uncover a ternary complex, where the originally bound Fis protein is partially dissociated from DNA. The simulations support a three-state sequential kinetic model (N ⇄ I → D) for facilitated dissociation, thus explaining the concentration-dependent dissociation.

The Aggregation Free Energy Landscapes of Polyglutamine Repeats.

Aggregates of proteins containing polyglutamine (polyQ) repeats are strongly associated with several neurodegenerative diseases. The length of the repeats correlates with the severity of the disease. Previous studies have shown that pure polyQ peptides aggregate by nucleated growth polymerization and that the size of the critical nucleus (n*) decreases from tetrameric to dimeric and monomeric as length increases from Q18 to Q26. Why the critical nucleus size changes with repeat-length has been unclear. Using the associative memory, water-mediated, structure and energy model, we construct the aggregation free energy landscapes for polyQ peptides of different repeat-lengths. These studies show that the monomer of the shorter repeat-length (Q20) prefers an extended conformation and that its aggregation indeed has a trimeric nucleus (n* ∼ 3), while a longer repeat-length monomer (Q30) prefers a ß-hairpin conformation which then aggregates in a downhill fashion at 0.1 mM. For an intermediate length peptide (Q26), there is an equal preference for hairpin and extended forms in the monomer which leads to a mixed inhomogeneous nucleation mechanism for fibrils. The predicted changes of monomeric structure and nucleation mechanism are confirmed by studying the aggregation free energy profile for a polyglutamine repeat with site-specific PG mutations that favor the hairpin form, giving results in harmony with experiments on this system.

Free energy landscapes for initiation and branching of protein aggregation.

Experiments on artificial multidomain protein constructs have probed the early stages of aggregation processes, but structural details of the species that initiate aggregation remain elusive. Using the associative-memory, water-mediated, structure and energy model known as AWSEM, a transferable coarse-grained protein model, we performed simulations of fused constructs composed of up to four copies of the Titin I27 domain or its mutant I27* (I59E). Free energy calculations enable us to quantify the conditions under which such multidomain constructs will spontaneously misfold. Consistent with experimental results, the dimer of I27 is found to be the smallest spontaneously misfolding construct. Our results show how structurally distinct misfolded states can be stabilized under different thermodynamic conditions, and this result provides a plausible link between the single-molecule misfolding experiments under native conditions and aggregation experiments under denaturing conditions. The conditions for spontaneous misfolding are determined by the interplay among temperature, effective local protein concentration, and the strength of the interdomain interactions. Above the folding temperature, fusing additional domains to the monomer destabilizes the native state, and the entropically stabilized amyloid-like state is favored. Because it is primarily energetically stabilized, the domain-swapped state is more likely to be important under native conditions. Both protofibril-like and branching structures are found in annealing simulations starting from extended structures, and these structures suggest a possible connection between the existence of multiple amyloidogenic segments in each domain and the formation of branched, amorphous aggregates as opposed to linear fibrillar structures.

Frustration in the energy landscapes of multidomain protein misfolding.

Frustration from strong interdomain interactions can make misfolding a more severe problem in multidomain proteins than in single-domain proteins. On the basis of bioinformatic surveys, it has been suggested that lowering the sequence identity between neighboring domains is one of nature's solutions to the multidomain misfolding problem. We investigate folding of multidomain proteins using the associative-memory, water-mediated, structure and energy model (AWSEM), a predictive coarse-grained protein force field. We find that reducing sequence identity not only decreases the formation of domain-swapped contacts but also decreases the formation of strong self-recognition contacts between ß-strands with high hydrophobic content. The ensembles of misfolded structures that result from forming these amyloid-like interactions are energetically disfavored compared with the native state, but entropically favored. Therefore, these ensembles are more stable than the native ensemble under denaturing conditions, such as high temperature. Domain-swapped contacts compete with self-recognition contacts in forming various trapped states, and point mutations can shift the balance between the two types of interaction. We predict that multidomain proteins that lack these specific strong interdomain interactions should fold reliably.

Predictive energy landscapes for protein-protein association.

We investigate protein-protein association using the associative-memory, water-mediated, structure, and energy model (AWSEM), a coarse-grained protein folding model that has been optimized using energy-landscape theory. The potential was originally parameterized by enforcing a funneled nature for a database of dimeric interfaces but was later further optimized to create funneled folding landscapes for individual monomeric proteins. The ability of the model to predict interfaces was not tested previously. The present results show that simulated annealing of the model indeed is able to predict successfully the native interfaces of eight homodimers and four heterodimers, thus amounting to a flexible docking algorithm. We go on to address the relative importance of monomer geometry, flexibility, and nonnative intermonomeric contacts in the association process for the homodimers. Monomer surface geometry is found to be important in determining the binding interface, but it is insufficient. Using a uniform binding potential rather than the water-mediated potential results in sampling of misbound structures that are geometrically preferred but are nonetheless energetically disfavored by AWSEM, as well as in nature. Depending on the stability of the unbound monomers, nonnative contacts play different roles in the association process. For unstable monomers, thermodynamic states stabilized by nonnative interactions correspond to productive, on-pathway intermediates and can, therefore, catalyze binding through a fly-casting mechanism. For stable monomers, in contrast, states stabilized by nonnative interactions generally correspond to traps that impede binding.

Simultaneous determination of neochlorogenic acid, chlorogenic acid, cryptochlorogenic acid and geniposide in rat plasma by UPLC-MS/MS and its application to a pharmacokinetic study after administration of Reduning injection.

A simple, specific and sensitive ultra-performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) method was established and validated for simultaneous determination of neochlorogenic acid, chlorogenic acid, cryptochlorogenic acid and geniposide in rat plasma using puerarin as an internal standard (IS). Plasma samples were pretreated by a one-step direct protein precipitation procedure with acetonitrile after acidified using as little as 50 µL plasma. Chromatographic separation was performed on an Acquity BEH C18 column (100 × 2.1 mm, 1.7 µm) at a flow rate of 0.2 mL/min by a gradient elution, using 0.2% acetic acid-methanol as mobile phase. The detection was performed on a triple quadrupole tandem mass spectrometer by multiple reaction monitoring via electrospray ionization source with negative ion mode. Calibration curves showed good linearity (r > 0.995) over wide concentration ranges. The intra- and inter-day precisions were <15%, and the accuracy was within ±8.0%. The validated method was successfully applied to a pharmacokinetic study of the four bioactive components in rats after intravenous administration of Reduning injection.