ProtoCode: Leveraging large language models (LLMs) for automated generation of machine-readable PCR protocols from scientific publications.

Protocol standardization and sharing are crucial for reproducibility in life sciences. In spite of numerous efforts for standardized protocol description, adherence to these standards in literature remains largely inconsistent. Curation of protocols are especially challenging due to the labor intensive process, requiring expert domain knowledge of each experimental procedure. Recent advancements in Large Language Models (LLMs) offer a promising solution to interpret and curate knowledge from complex scientific literature. In this work, we develop ProtoCode, a tool leveraging fine-tune LLMs to curate protocols into intermediate representation formats which can be interpretable by both human and machine interfaces. Our proof-of-concept, focused on polymerase chain reaction (PCR) protocols, retrieves information from PCR protocols at an accuracy ranging 69-100 % depending on the information content. In all tested protocols, we demonstrate that ProtoCode successfully converts literature-based protocols into correct operational files for multiple thermal cycler systems. In conclusion, ProtoCode can alleviate labor intensive curation and standardization of life science protocols to enhance research reproducibility by providing a reliable, automated means to process and standardize protocols. ProtoCode is freely available as a web server at https://curation.taxila.io/ProtoCode/.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event.

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Parallel nonfunctionalization of CK1δ/ε kinase ohnologs following a whole-genome duplication event.

Whole genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in S. cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post WGD. In S. cerevisiae, HRR25 encodes the casein kinase (CK) 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post and non-WGD species, and have non-conserved cellular localization, consistent with their ongoing loss of function. The analysis in N. castelli shows that the non-complementing ohnolog is expressed at a lower level and has become non-essential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

A universal sequencing read interpreter.

Massively parallel DNA sequencing has led to the rapid growth of highly multiplexed experiments in biology. These experiments produce unique sequencing results that require specific analysis pipelines to decode highly structured reads. However, no versatile framework that interprets sequencing reads to extract their encoded information for downstream biological analysis has been developed. Here, we report INTERSTELLAR (interpretation, scalable transformation, and emulation of large-scale sequencing reads) that decodes data values encoded in theoretically any type of sequencing read and translates them into sequencing reads of another structure of choice. We demonstrated that INTERSTELLAR successfully extracted information from a range of short- and long-read sequencing reads and translated those of single-cell (sc)RNA-seq, scATAC-seq, and spatial transcriptomics to be analyzed by different software tools that have been developed for conceptually the same types of experiments. INTERSTELLAR will greatly facilitate the development of sequencing-based experiments and sharing of data analysis pipelines.

Generation of functional hepatocytes by forward programming with nuclear receptors.

Production of large quantities of hepatocytes remains a major challenge for a number of clinical applications in the biomedical field. Directed differentiation of human pluripotent stem cells (hPSCs) into hepatocyte-like cells (HLCs) provides an advantageous solution and a number of protocols have been developed for this purpose. However, these methods usually follow different steps of liver development in vitro, which is time consuming and requires complex culture conditions. In addition, HLCs lack the full repertoire of functionalities characterising primary hepatocytes. Here, we explore the interest of forward programming to generate hepatocytes from hPSCs and to bypass these limitations. This approach relies on the overexpression of three hepatocyte nuclear factors (HNF1A, HNF6, and FOXA3) in combination with different nuclear receptors expressed in the adult liver using the OPTi-OX platform. Forward programming allows for the rapid production of hepatocytes (FoP-Heps) with functional characteristics using a simplified process. We also uncovered that the overexpression of nuclear receptors such as RORc can enhance specific functionalities of FoP-Heps thereby validating its role in lipid/glucose metabolism. Together, our results show that forward programming could offer a versatile alternative to direct differentiation for generating hepatocytes in vitro.

Barcode fusion genetics-protein-fragment complementation assay (BFG-PCA): tools and resources that expand the potential for binary protein interaction discovery.

Barcode fusion genetics (BFG) utilizes deep sequencing to improve the throughput of protein-protein interaction (PPI) screening in pools. BFG has been implemented in Yeast two-hybrid (Y2H) screens (BFG-Y2H). While Y2H requires test protein pairs to localize in the nucleus for reporter reconstruction, dihydrofolate reductase protein-fragment complementation assay (DHFR-PCA) allows proteins to localize in broader subcellular contexts and proves to be largely orthogonal to Y2H. Here, we implemented BFG to DHFR-PCA (BFG-PCA). This plasmid-based system can leverage ORF collections across model organisms to perform comparative analysis, unlike the original DHFR-PCA that requires yeast genomic integration. The scalability and quality of BFG-PCA were demonstrated by screening human and yeast interactions for >11 000 bait-prey pairs. BFG-PCA showed high-sensitivity and high-specificity for capturing known interactions for both species. BFG-Y2H and BFG-PCA capture distinct sets of PPIs, which can partially be explained based on the domain orientation of the reporter tags. BFG-PCA is a high-throughput protein interaction technology to interrogate binary PPIs that exploits clone collections from any species of interest, expanding the scope of PPI assays.

Myopalladin knockout mice develop cardiac dilation and show a maladaptive response to mechanical pressure overload.

Myopalladin (MYPN) is a striated muscle-specific immunoglobulin domain-containing protein located in the sarcomeric Z-line and I-band. MYPN gene mutations are causative for dilated (DCM), hypertrophic, and restrictive cardiomyopathy. In a yeast two-hybrid screening, MYPN was found to bind to titin in the Z-line, which was confirmed by microscale thermophoresis. Cardiac analyses of MYPN knockout (MKO) mice showed the development of mild cardiac dilation and systolic dysfunction, associated with decreased myofibrillar isometric tension generation and increased resting tension at longer sarcomere lengths. MKO mice exhibited a normal hypertrophic response to transaortic constriction (TAC), but rapidly developed severe cardiac dilation and systolic dysfunction, associated with fibrosis, increased fetal gene expression, higher intercalated disc fold amplitude, decreased calsequestrin-2 protein levels, and increased desmoplakin and SORBS2 protein levels. Cardiomyocyte analyses showed delayed Ca2+ release and reuptake in unstressed MKO mice as well as reduced Ca2+ spark amplitude post-TAC, suggesting that altered Ca2+ handling may contribute to the development of DCM in MKO mice.

Machine learning approach for discrimination of genotypes based on bright-field cellular images.

Morphological profiling is a combination of established optical microscopes and cutting-edge machine vision technologies, which stacks up successful applications in high-throughput phenotyping. One major question is how much information can be extracted from an image to identify genetic differences between cells. While fluorescent microscopy images of specific organelles have been broadly used for single-cell profiling, the potential ability of bright-field (BF) microscopy images of label-free cells remains to be tested. Here, we examine whether single-gene perturbation can be discriminated based on BF images of label-free cells using a machine learning approach. We acquired hundreds of BF images of single-gene mutant cells, quantified single-cell profiles consisting of texture features of cellular regions, and constructed a machine learning model to discriminate mutant cells from wild-type cells. Interestingly, the mutants were successfully discriminated from the wild type (area under the receiver operating characteristic curve = 0.773). The features that contributed to the discrimination were identified, and they included those related to the morphology of structures that appeared within cellular regions. Furthermore, functionally close gene pairs showed similar feature profiles of the mutant cells. Our study reveals that single-gene mutant cells can be discriminated from wild-type cells based on BF images, suggesting the potential as a useful tool for mutant cell profiling.

Complete Genome Sequence of Bacillus sp. Strain KH172YL63, Isolated from Deep-Sea Sediment.

Bacillus sp. strain KH172YL63 is a Gram-positive bacterium isolated from the deep-sea floor surface sediment at 3,308 m below sea level in the Nankai Trough in Japan. Here, we report the complete genome sequence of Bacillus sp. strain KH172YL63, which has a genome size of 4,251,700 bp and a G+C content of 44.8%.

Myopalladin promotes muscle growth through modulation of the serum response factor pathway.

BACKGROUND: Myopalladin (MYPN) is a striated muscle-specific, immunoglobulin-containing protein located in the Z-line and I-band of the sarcomere as well as the nucleus. Heterozygous MYPN gene mutations are associated with hypertrophic, dilated, and restrictive cardiomyopathy, and homozygous loss-of-function truncating mutations have recently been identified in patients with cap myopathy, nemaline myopathy, and congenital myopathy with hanging big toe. METHODS: Constitutive MYPN knockout (MKO) mice were generated, and the role of MYPN in skeletal muscle was studied through molecular, cellular, biochemical, structural, biomechanical, and physiological studies in vivo and in vitro. RESULTS: MKO mice were 13% smaller compared with wild-type controls and exhibited a 48% reduction in myofibre cross-sectional area (CSA) and significantly increased fibre number. Similarly, reduced myotube width was observed in MKO primary myoblast cultures. Biomechanical studies showed reduced isometric force and power output in MKO mice as a result of the reduced CSA, whereas the force developed by each myosin molecular motor was unaffected. While the performance by treadmill running was similar in MKO and wild-type mice, MKO mice showed progressively decreased exercise capability, Z-line damage, and signs of muscle regeneration following consecutive days of downhill running. Additionally, MKO muscle exhibited progressive Z-line widening starting from 8 months of age. RNA-sequencing analysis revealed down-regulation of serum response factor (SRF)-target genes in muscles from postnatal MKO mice, important for muscle growth and differentiation. The SRF pathway is regulated by actin dynamics as binding of globular actin to the SRF-cofactor myocardin-related transcription factor A (MRTF-A) prevents its translocation to the nucleus where it binds and activates SRF. MYPN was found to bind and bundle filamentous actin as well as interact with MRTF-A. In particular, while MYPN reduced actin polymerization, it strongly inhibited actin depolymerization and consequently increased MRTF-A-mediated activation of SRF signalling in myogenic cells. Reduced myotube width in MKO primary myoblast cultures was rescued by transduction with constitutive active SRF, demonstrating that MYPN promotes skeletal muscle growth through activation of the SRF pathway. CONCLUSIONS: Myopalladin plays a critical role in the control of skeletal muscle growth through its effect on actin dynamics and consequently the SRF pathway. In addition, MYPN is important for the maintenance of Z-line integrity during exercise and aging. These results suggest that muscle weakness in patients with biallelic MYPN mutations may be associated with reduced myofibre CSA and SRF signalling and that the disease phenotype may be aggravated by exercise.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs.

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

Complete Genome Sequence of Psychrobacter sp. Strain KH172YL61, Isolated from Deep-Sea Sediments in the Nankai Trough, Japan.

Psychrobacter sp. strain KH172YL61 is a Gram-negative bacterium isolated from deep-sea sediment in the Nankai Trough in Japan. Here, we report the complete genome sequence of this strain, which has a genome size of 3.19 Mb, with a G+C content of 44.0%.

Fast and global detection of periodic sequence repeats in large genomic resources.

Periodically repeating DNA and protein elements are involved in various important biological events including genomic evolution, gene regulation, protein complex formation, and immunity. Notably, the currently used genome editing tools such as ZFNs, TALENs, and CRISPRs are also all associated with periodically repeating biomolecules of natural organisms. Despite the biological importance of periodically repeating sequences and the expectation that new genome editing modules could be discovered from such periodical repeats, no software that globally detects such structured elements in large genomic resources in a high-throughput and unsupervised manner has been developed. We developed new software, SPADE (Search for Patterned DNA Elements), that exhaustively explores periodic DNA and protein repeats from large-scale genomic datasets based on k-mer periodicity evaluation. With a simple constraint, sequence periodicity, SPADE captured reported genome-editing-associated sequences and other protein families involving repeating domains such as tetratricopeptide, ankyrin and WD40 repeats with better performance than the other software designed for limited sets of repetitive biomolecular sequences, suggesting the high potential of this software to contribute to the discovery of new biological events and new genome editing modules.

The nebulin SH3 domain is dispensable for normal skeletal muscle structure but is required for effective active load bearing in mouse.

Nemaline myopathy (NM) is a congenital myopathy with an estimated incidence of 150,000 live births. It is caused by mutations in thin filament components, including nebulin, which accounts for about 50% of the cases. The identification of NM cases with nonsense mutations resulting in loss of the extreme C-terminal SH3 domain of nebulin suggests an important role of the nebulin SH3 domain, which is further supported by the recent demonstration of its role in IGF-1-induced sarcomeric actin filament formation through targeting of N-WASP to the Z-line. To provide further insights into the functional significance of the nebulin SH3 domain in the Z-disk and to understand the mechanisms by which truncations of nebulin lead to NM, we took two approaches: (1) an affinity-based proteomic screening to identify novel interaction partners of the nebulin SH3 domain; and (2) generation and characterization of a novel knockin mouse model with a premature stop codon in the nebulin gene, eliminating its C-terminal SH3 domain (NebΔSH3 mouse). Surprisingly, detailed analyses of NebΔSH3 mice revealed no structural or histological skeletal muscle abnormalities and no changes in gene expression or localization of interaction partners of the nebulin SH3 domain, including myopalladin, palladin, zyxin and N-WASP. Also, no significant effect on peak isometric stress production, passive tensile stress or Young's modulus was found. However, NebΔSH3 muscle displayed a slightly altered force-frequency relationship and was significantly more susceptible to eccentric contraction-induced injury, suggesting that the nebulin SH3 domain protects against eccentric contraction-induced injury and possibly plays a role in fine-tuning the excitation-contraction coupling mechanism.

ß1-Adrenergic receptors increase UCP1 in human MADS brown adipocytes and rescue cold-acclimated ß3-adrenergic receptor-knockout mice via nonshivering thermogenesis.

With the finding that brown adipose tissue is present and negatively correlated to obesity in adult man, finding the mechanism(s) of how to activate brown adipose tissue in humans could be important in combating obesity, type 2 diabetes, and their complications. In mice, the main regulator of nonshivering thermogenesis in brown adipose tissue is norepinephrine acting predominantly via ß(3)-adrenergic receptors. However, vast majorities of ß(3)-adrenergic agonists have so far not been able to stimulate human ß(3)-adrenergic receptors or brown adipose tissue activity, and it was postulated that human brown adipose tissue could be regulated instead by ß(1)-adrenergic receptors. Therefore, we have investigated the signaling pathways, specifically pathways to nonshivering thermogenesis, in mice lacking ß(3)-adrenergic receptors. Wild-type and ß(3)-knockout mice were either exposed to acute cold (up to 12 h) or acclimated for 7 wk to cold, and parameters related to metabolism and brown adipose tissue function were investigated. ß(3)-knockout mice were able to survive both acute and prolonged cold exposure due to activation of ß(1)-adrenergic receptors. Thus, in the absence of ß(3)-adrenergic receptors, ß(1)-adrenergic receptors are effectively able to signal via cAMP to elicit cAMP-mediated responses and to recruit and activate brown adipose tissue. In addition, we found that in human multipotent adipose-derived stem cells differentiated into functional brown adipocytes, activation of either ß(1)-adrenergic receptors or ß(3)-adrenergic receptors was able to increase UCP1 mRNA and protein levels. Thus, in humans, ß(1)-adrenergic receptors could play an important role in regulating nonshivering thermogenesis.

Mammary-derived growth inhibitor alters traffic of EGFR and induces a novel form of cetuximab resistance.

PURPOSE: Only few predictive factors for the clinical activity of anti-epidermal growth factor receptor (EGFR) therapy are available. Mammary-derived growth inhibitor (MDGI) is a small cytosolic protein suggested to play a role in the differentiation of epithelial cells. Here, we have investigated the effect of MDGI expression on the EGFR signaling and cetuximab responsiveness of cancer cells. EXPERIMENTAL DESIGN: MDGI mRNA expression was investigated in clinical breast and lung cancer samples and in nontransformed and malignant cell lines. The effect of ectopic expression of MDGI on EGFR, ErbB2, and integrin function and traffic was investigated in breast and lung cancer cell lines using multiple methods. The effect of anti-EGFR agents on these cells were tested by cell proliferation measurements and by assessing tumor growth of breast cancer cells in cetuximab treated and control athymic nude mice. RESULTS: Here, we show that although MDGI is absent in cultured cell lines because of epigenetic silencing, MDGI mRNA is expressed in 40% of clinical breast carcinomas and 85% of lung cancers. Ectopic expression of MDGI rendered breast and lung cancer cells resistant to the anti-EGFR antibody cetuximab in vitro and in an orthotopic breast cancer xenograft model in vivo. When expressed in cancer cells, MDGI induces intracellular accumulation of EGFR, but not ErbB2, and the internalized receptor is phosphorylated and not degraded. CONCLUSIONS: MDGI-driven inherent desensitization of cancer cells is a novel molecular mechanism for resistance to the anti-EGFR antibody therapy, and MDGI may be a biomarker for responsiveness to anti-EGFR antibody therapy.

Myotube formation on micro-patterned glass: intracellular organization and protein distribution in C2C12 skeletal muscle cells.

Proliferation and fusion of myoblasts are needed for the generation and repair of multinucleated skeletal muscle fibers in vivo. Studies of myocyte differentiation, cell fusion, and muscle repair are limited by an appropriate in vitro muscle cell culture system. We developed a novel cell culture technique [two-dimensional muscle syncytia (2DMS) technique] that results in formation of myotubes, organized in parallel much like the arrangement in muscle tissue. This technique is based on UV lithography-produced micro-patterned glass on which conventionally cultured C2C12 myoblasts proliferate, align, and fuse to neatly arranged contractile myotubes in parallel arrays. Combining this technique with fluorescent microscopy, we observed alignment of actin filament bundles and a perinuclear distribution of glucose transporter 4 after myotube formation. Newly formed myotubes contained adjacently located MyoD-positive and MyoD-negative nuclei, suggesting fusion of MyoD-positive and MyoD-negative cells. In comparison, the closely related myogenic factor Myf5 did not exhibit this pattern of distribution. Furthermore, cytoplasmic patches of MyoD colocalized with bundles of filamentous actin near myotube nuclei. At later stages of differentiation, all nuclei in the myotubes were MyoD negative. The 2DMS system is thus a useful tool for studies on muscle alignment, differentiation, fusion, and subcellular protein localization.

Diphenylene iodonium stimulates glucose uptake in skeletal muscle cells through mitochondrial complex I inhibition and activation of AMP-activated protein kinase.

NADPH oxidase inhibitors such as diphenylene iodonium (DPI) and apocynin lower whole body and blood glucose levels and improve diabetes when administered to rodents. Skeletal muscle has an important role in managing glucose homeostasis and we have used L6 cells, C(2)C(12) cells and primary muscle cells as model systems to investigate whether these drugs regulate glucose uptake in skeletal muscle cells. The data presented in this study show that apocynin does not affect glucose uptake in skeletal muscle cells in culture. Tat gp91ds, a chimeric peptide that inhibits NADPH oxidase activity, also failed to affect glucose uptake and we found no significant evidence of NADPH oxidase (subunits tested were Nox4, p22phox, gp91phox and p47phox mRNA) in skeletal muscle cells in culture. However, DPI increases basal and insulin-stimulated glucose uptake in L6 cells, C(2)C(12) cells and primary muscle cells. Detailed studies on L6 cells demonstrate that the increase of glucose uptake is via a mechanism independent of phosphoinositide-3 kinase (PI3K)/Akt but dependent on AMP-activated protein kinase (AMPK). We postulate that DPI through inhibition of mitochondrial complex 1 and decreases in oxygen consumption, leading to decreases of ATP and activation of AMPK, stimulates glucose uptake in skeletal muscle cells.