Internet-Connected Cortical Organoids for Project-Based Stem Cell and Neuroscience Education.

The introduction of Internet-connected technologies to the classroom has the potential to revolutionize STEM education by allowing students to perform experiments in complex models that are unattainable in traditional teaching laboratories. By connecting laboratory equipment to the cloud, we introduce students to experimentation in pluripotent stem cell (PSC)-derived cortical organoids in two different settings: using microscopy to monitor organoid growth in an introductory tissue culture course and using high-density (HD) multielectrode arrays (MEAs) to perform neuronal stimulation and recording in an advanced neuroscience mathematics course. We demonstrate that this approach develops interest in stem cell and neuroscience in the students of both courses. All together, we propose cloud technologies as an effective and scalable approach for complex project-based university training.

Internet-connected cortical organoids for project-based stem cell and neuroscience education.

The introduction of internet-connected technologies to the classroom has the potential to revolutionize STEM education by allowing students to perform experiments in complex models that are unattainable in traditional teaching laboratories. By connecting laboratory equipment to the cloud, we introduce students to experimentation in pluripotent stem cell-derived cortical organoids in two different settings: Using microscopy to monitor organoid growth in an introductory tissue culture course, and using high density multielectrode arrays to perform neuronal stimulation and recording in an advanced neuroscience mathematics course. We demonstrate that this approach develops interest in stem cell and neuroscience in the students of both courses. All together, we propose cloud technologies as an effective and scalable approach for complex project-based university training.

Cloud-controlled microscopy enables remote project-based biology education in underserved Latinx communities.

Project-based learning (PBL) has long been recognized as an effective way to teach complex biology concepts. However, not all institutions have the resources to facilitate effective project-based coursework for students. We have developed a framework for facilitating PBL using remote-controlled internet-connected microscopes. Through this approach, one lab facility can host an experiment for many students around the world simultaneously. Experiments on this platform can be run on long timescales and with materials that are typically unavailable to high school classrooms. This allows students to perform novel research projects rather than just repeating standard classroom experiments. To investigate the impact of this program, we designed and ran six user studies with students worldwide. All experiments were hosted in Santa Cruz and San Francisco, California, with observations and decisions made remotely by the students using their personal computers and cellphones. In surveys gathered after the experiments, students reported increased excitement for science and a greater desire to pursue a career in STEM. This framework represents a novel, scalable, and effective PBL approach that has the potential to democratize biology and STEM education around the world.

Picroscope: low-cost system for simultaneous longitudinal biological imaging.

Simultaneous longitudinal imaging across multiple conditions and replicates has been crucial for scientific studies aiming to understand biological processes and disease. Yet, imaging systems capable of accomplishing these tasks are economically unattainable for most academic and teaching laboratories around the world. Here, we propose the Picroscope, which is the first low-cost system for simultaneous longitudinal biological imaging made primarily using off-the-shelf and 3D-printed materials. The Picroscope is compatible with standard 24-well cell culture plates and captures 3D z-stack image data. The Picroscope can be controlled remotely, allowing for automatic imaging with minimal intervention from the investigator. Here, we use this system in a range of applications. We gathered longitudinal whole organism image data for frogs, zebrafish, and planaria worms. We also gathered image data inside an incubator to observe 2D monolayers and 3D mammalian tissue culture models. Using this tool, we can measure the behavior of entire organisms or individual cells over long-time periods.

Cell type-specific roles of PAR1 in Coxsackievirus B3 infection.

Protease-activated receptor 1 (PAR1) is widely expressed in humans and mice, and is activated by a variety of proteases, including thrombin. Recently, we showed that PAR1 contributes to the innate immune response to viral infection. Mice with a global deficiency of PAR1 expressed lower levels of CXCL10 and had increased Coxsackievirus B3 (CVB3)-induced myocarditis compared with control mice. In this study, we determined the effect of cell type-specific deletion of PAR1 in cardiac myocytes (CMs) and cardiac fibroblasts (CFs) on CVB3-induced myocarditis. Mice lacking PAR1 in either CMs or CFs exhibited increased CVB3 genomes, inflammatory infiltrates, macrophages and inflammatory mediators in the heart and increased CVB3-induced myocarditis compared with wild-type controls. Interestingly, PAR1 enhanced poly I:C induction of CXCL10 in rat CFs but not in rat neonatal CMs. Importantly, activation of PAR1 reduced CVB3 replication in murine embryonic fibroblasts and murine embryonic cardiac myocytes. In addition, we showed that PAR1 reduced autophagy in murine embryonic fibroblasts and rat H9c2 cells, which may explain how PAR1 reduces CVB3 replication. These data suggest that PAR1 on CFs protects against CVB3-induced myocarditis by enhancing the anti-viral response whereas PAR1 on both CMs and fibroblasts inhibits viral replication.

Cancer Therapy-Associated Thrombosis.

[Figure: see text].

Drug binding dynamics of the dimeric SARS-CoV-2 main protease, determined by molecular dynamics simulation.

We performed molecular dynamics simulation of the dimeric SARS-CoV-2 (severe acute respiratory syndrome corona virus 2) main protease (Mpro) to examine the binding dynamics of small molecular ligands. Seven HIV inhibitors, darunavir, indinavir, lopinavir, nelfinavir, ritonavir, saquinavir, and tipranavir, were used as the potential lead drugs to investigate access to the drug binding sites in Mpro. The frequently accessed sites on Mpro were classified based on contacts between the ligands and the protein, and the differences in site distributions of the encounter complex were observed among the ligands. All seven ligands showed binding to the active site at least twice in 28 simulations of 200 ns each. We further investigated the variations in the complex structure of the active site with the ligands, using microsecond order simulations. Results revealed a wide variation in the shapes of the binding sites and binding poses of the ligands. Additionally, the C-terminal region of the other chain often interacted with the ligands and the active site. Collectively, these findings indicate the importance of dynamic sampling of protein-ligand complexes and suggest the possibilities of further drug optimisations.

Quantification of Viral and Host Biomarkers in the Liver of Rhesus Macaques: A Longitudinal Study of Zaire Ebolavirus Strain Kikwit (EBOV/Kik).

Zaire ebolavirus (EBOV) causes Ebola virus disease (EVD), which carries a fatality rate between 25% and 90% in humans. Liver pathology is a hallmark of terminal EVD; however, little is known about temporal disease progression. We used multiplexed fluorescent immunohistochemistry and in situ hybridization in combination with whole slide imaging and image analysis (IA) to quantitatively characterize temporospatial signatures of viral and host factors as related to EBOV pathogenesis. Eighteen rhesus monkeys euthanized between 3 and 8 days post-infection, and 3 uninfected controls were enrolled in this study. Compared with semiquantitative histomorphologic ordinal scoring, quantitative IA detected subtle and progressive features of early and terminal EVD that was not feasible with routine approaches. Sinusoidal macrophages were the earliest cells to respond to infection, expressing proinflammatory cytokine interleukin 6 (IL6) mRNA, which was subsequently also observed in fibrovascular compartments. The mRNA of interferon-stimulated gene-15 (ISG-15), also known as ISG15 ubiquitin like modifier (ISG15), was observed early, with a progressive and ubiquitous hybridization signature involving mesenchymal and epithelial compartments. ISG-15 mRNA was prominent near infected cells, but not in infected cells, supporting the hypothesis that bystander cells produce a robust interferon gene response. This study contributes to our current understanding of early EVD progression and illustrates the value that digital pathology and quantitative IA serve in infectious disease research.

Circulating Extracellular Vesicle Tissue Factor Activity During Orthohantavirus Infection Is Associated With Intravascular Coagulation.

BACKGROUND: Puumala orthohantavirus (PUUV) causes hemorrhagic fever with renal syndrome (HFRS). Patients with HFRS have an activated coagulation system with increased risk of disseminated intravascular coagulation (DIC) and venous thromboembolism (VTE). The aim of the study was to determine whether circulating extracellular vesicle tissue factor (EVTF) activity levels associates with DIC and VTE (grouped as intravascular coagulation) in HFRS patients. METHODS: Longitudinal samples were collected from 88 HFRS patients. Patients were stratified into groups of those with intravascular coagulation (n = 27) and those who did not (n = 61). We measured levels of circulating EVTF activity, fibrinogen, activated partial prothrombin time, D-dimer, tissue plasminogen activator (tPA), plasminogen activator inhibitor 1 (PAI-1), and platelets. RESULTS: Plasma EVTF activity was transiently increased during HFRS. Levels of EVTF activity were significantly associated with plasma tPA and PAI-1, suggesting that endothelial cells could be a potential source. Patients with intravascular coagulation had significantly higher peak EVTF activity levels compared with those who did not, even after adjustment for sex and age. The peak EVTF activity value predicting intravascular coagulation was 0.51 ng/L with 63% sensitivity and 61% specificity with area under the curve = 0.63 (95% confidence interval, 0.51-0.76) and P = .046. CONCLUSIONS: Plasma EVTF activity during HFRS is associated with intravascular coagulation.

An average-case sublinear forward algorithm for the haploid Li and Stephens model.

BACKGROUND: Hidden Markov models of haplotype inheritance such as the Li and Stephens model allow for computationally tractable probability calculations using the forward algorithm as long as the representative reference panel used in the model is sufficiently small. Specifically, the monoploid Li and Stephens model and its variants are linear in reference panel size unless heuristic approximations are used. However, sequencing projects numbering in the thousands to hundreds of thousands of individuals are underway, and others numbering in the millions are anticipated. RESULTS: To make the forward algorithm for the haploid Li and Stephens model computationally tractable for these datasets, we have created a numerically exact version of the algorithm with observed average case sublinear runtime with respect to reference panel size k when tested against the 1000 Genomes dataset. CONCLUSIONS: We show a forward algorithm which avoids any tradeoff between runtime and model complexity. Our algorithm makes use of two general strategies which might be applicable to improving the time complexity of other future sequence analysis algorithms: sparse dynamic programming matrices and lazy evaluation.

Superbubbles, Ultrabubbles, and Cacti.

A superbubble is a type of directed acyclic subgraph with single distinct source and sink vertices. In genome assembly and genetics, the possible paths through a superbubble can be considered to represent the set of possible sequences at a location in a genome. Bidirected and biedged graphs are a generalization of digraphs that are increasingly being used to more fully represent genome assembly and variation problems. In this study, we define snarls and ultrabubbles, generalizations of superbubbles for bidirected and biedged graphs, and give an efficient algorithm for the detection of these more general structures. Key to this algorithm is the cactus graph, which, we show, encodes the nested decomposition of a graph into snarls and ultrabubbles within its structure. We propose and demonstrate empirically that this decomposition on bidirected and biedged graphs solves a fundamental problem by defining genetic sites for any collection of genomic variations, including complex structural variations, without need for any single reference genome coordinate system. Further, the nesting of the decomposition gives a natural way to describe and model variations contained within large variations, a case not currently dealt with by existing formats [e.g., variant cell format (VCF)].

Temperature-Sensitive Substrate and Product Binding Underlie Temperature-Compensated Phosphorylation in the Clock.

Temperature compensation is a striking feature of the circadian clock. Here we investigate biochemical mechanisms underlying temperature-compensated, CKIδ-dependent multi-site phosphorylation in mammals. We identify two mechanisms for temperature-insensitive phosphorylation at higher temperature: lower substrate affinity to CKIδ-ATP complex and higher product affinity to CKIδ-ADP complex. Inhibitor screening of ADP-dependent phosphatase activity of CKIδ identified aurintricarboxylic acid (ATA) as a temperature-sensitive kinase activator. Docking simulation of ATA and mutagenesis experiment revealed K224D/K224E mutations in CKIδ that impaired product binding and temperature-compensated primed phosphorylation. Importantly, K224D mutation shortens behavioral circadian rhythms and changes the temperature dependency of SCN's circadian period. Interestingly, temperature-compensated phosphorylation was evolutionary conserved in yeast. Molecular dynamics simulation and X-ray crystallography demonstrate that an evolutionally conserved CKI-specific domain around K224 can provide a structural basis for temperature-sensitive substrate and product binding. Surprisingly, this domain can confer temperature compensation on a temperature-sensitive TTBK1. These findings suggest the temperature-sensitive substrate- and product-binding mechanisms underlie temperature compensation.

Microvesicle-associated tissue factor procoagulant activity for the preoperative diagnosis of ovarian cancer.

BACKGROUND: Tissue factor (TF) is involved in tumor growth and metastasis and contributes to venous thromboembolism (VTE) in cancer, including gynecological malignancies. The diagnostic value of microvesicle-associated TF procoagulant activity (MV TF PCA) in women with suspected ovarian cancer, however, has not been studied. OBJECTIVE: To evaluate MV TF PCA as a diagnostic tool in women with an ovarian mass of unknown etiology and as a predictive biomarker for perioperative VTE. METHODS: Plasma MVs were isolated by high-speed centrifugation and analyzed for TF-specific PCA by single-stage clotting assay. In addition, plasma TF antigen and soluble P-selectin (sCD62P) were measured by ELISA. RESULTS: D-Dimer, MV TF PCA, and sCD62P, but not the tumor marker, CA-125, significantly differentiated patients with malignant (n=40) from those with benign tumors (n=15) and healthy controls (n=34). In cancer patients, only D-Dimer and CA-125 correlated with the FIGO stage. An abnormal D-dimer had the highest sensitivity for the diagnosis of cancer, while MV TF PCA above the ROC curve-derived cut-off value of 182U/mL had the highest specificity. By multivariate logistic regression analysis, addition of MV TF PCA conferred diagnostic benefit to the single variables, CA-125 (p=0.052) and D-dimer (p=0.019). Perioperative VTE occurred in 16% of cancer patients and was associated with an advanced FIGO stage, but not MV TF PCA. There was no difference in plasma TF antigen levels between study groups. CONCLUSIONS: MV TF PCA, but not plasma TF antigen, may provide valuable additional information for the diagnostic work-up of women with suspected ovarian cancer.

Perturbation analyses of intermolecular interactions.

Conformational fluctuations of a protein molecule are important to its function, and it is known that environmental molecules, such as water molecules, ions, and ligand molecules, significantly affect the function by changing the conformational fluctuations. However, it is difficult to systematically understand the role of environmental molecules because intermolecular interactions related to the conformational fluctuations are complicated. To identify important intermolecular interactions with regard to the conformational fluctuations, we develop herein (i) distance-independent and (ii) distance-dependent perturbation analyses of the intermolecular interactions. We show that these perturbation analyses can be realized by performing (i) a principal component analysis using conditional expectations of truncated and shifted intermolecular potential energy terms and (ii) a functional principal component analysis using products of intermolecular forces and conditional cumulative densities. We refer to these analyses as intermolecular perturbation analysis (IPA) and distance-dependent intermolecular perturbation analysis (DIPA), respectively. For comparison of the IPA and the DIPA, we apply them to the alanine dipeptide isomerization in explicit water. Although the first IPA principal components discriminate two states (the α state and PPII (polyproline II) + ß states) for larger cutoff length, the separation between the PPII state and the ß state is unclear in the second IPA principal components. On the other hand, in the large cutoff value, DIPA eigenvalues converge faster than that for IPA and the top two DIPA principal components clearly identify the three states. By using the DIPA biplot, the contributions of the dipeptide-water interactions to each state are analyzed systematically. Since the DIPA improves the state identification and the convergence rate with retaining distance information, we conclude that the DIPA is a more practical method compared with the IPA. To test the feasibility of the DIPA for larger molecules, we apply the DIPA to the ten-residue chignolin folding in explicit water. The top three principal components identify the four states (native state, two misfolded states, and unfolded state) and their corresponding eigenfunctions identify important chignolin-water interactions to each state. Thus, the DIPA provides the practical method to identify conformational states and their corresponding important intermolecular interactions with distance information.

Perturbational formulation of principal component analysis in molecular dynamics simulation.

Conformational fluctuations of a molecule are important to its function since such intrinsic fluctuations enable the molecule to respond to the external environmental perturbations. For extracting large conformational fluctuations, which predict the primary conformational change by the perturbation, principal component analysis (PCA) has been used in molecular dynamics simulations. However, several versions of PCA, such as Cartesian coordinate PCA and dihedral angle PCA (dPCA), are limited to use with molecules with a single dominant state or proteins where the dihedral angle represents an important internal coordinate. Other PCAs with general applicability, such as the PCA using pairwise atomic distances, do not represent the physical meaning clearly. Therefore, a formulation that provides general applicability and clearly represents the physical meaning is yet to be developed. For developing such a formulation, we consider the conformational distribution change by the perturbation with arbitrary linearly independent perturbation functions. Within the second order approximation of the Kullback-Leibler divergence by the perturbation, the PCA can be naturally interpreted as a method for (1) decomposing a given perturbation into perturbations that independently contribute to the conformational distribution change or (2) successively finding the perturbation that induces the largest conformational distribution change. In this perturbational formulation of PCA, (i) the eigenvalue measures the Kullback-Leibler divergence from the unperturbed to perturbed distributions, (ii) the eigenvector identifies the combination of the perturbation functions, and (iii) the principal component determines the probability change induced by the perturbation. Based on this formulation, we propose a PCA using potential energy terms, and we designate it as potential energy PCA (PEPCA). The PEPCA provides both general applicability and clear physical meaning. For demonstrating its power, we apply the PEPCA to an alanine dipeptide molecule in vacuum as a minimal model of a nonsingle dominant conformational biomolecule. The first and second principal components clearly characterize two stable states and the transition state between them. Positive and negative components with larger absolute values of the first and second eigenvectors identify the electrostatic interactions, which stabilize or destabilize each stable state and the transition state. Our result therefore indicates that PCA can be applied, by carefully selecting the perturbation functions, not only to identify the molecular conformational fluctuation but also to predict the conformational distribution change by the perturbation beyond the limitation of the previous methods.