PyRETIS 3: Conquering rare and slow events without boundaries

We present and discuss the advancements made in PyRETIS 3, the third instalment of our Python library for an efficient and user-friendly rare event simulation, focused to execute molecular simulations with replica exchange transition interface sampling (RETIS) and its variations. Apart from a general rewiring of the internal code towards a more modular structure, several recently developed sampling strategies have been implemented. These include recently developed Monte Carlo moves to increase path decorrelation and convergence rate, and new ensemble definitions to handle the challenges of long-lived metastable states and transitions with unbounded reactant and product states. Additionally, the post-analysis software PyVisa is now embedded in the main code, allowing fast use of machine-learning algorithms for clustering and visualising collective variables in the simulation data.

CMRsim-A python package for cardiovascular MR simulations incorporating complex motion and flow

PURPOSE

To present an open-source MR simulation framework that facilitates the incorporation of complex motion and flow for studying cardiovascular MR (CMR) acquisition and reconstruction.

METHODS

CMRsim is a Python package that allows simulation of CMR images using dynamic digital phantoms with complex motion as input. Two simulation paradigms are available, namely, numerical and analytical solutions to the Bloch equations, using a common motion representation. Competitive simulation speeds are achieved using TensorFlow for GPU acceleration. To demonstrate the capability of the package, one introductory and two advanced CMR simulation experiments are presented. The latter showcase phase-contrast imaging of turbulent flow downstream of a stenotic section and cardiac diffusion tensor imaging on a contracting left ventricle. Additionally, extensive documentation and example resources are provided.

RESULTS

The Bloch simulation with turbulent flow using approximately 1.5 million particles and a sequence duration of 710 ms for each of the seven different velocity encodings took a total of 29 min on a NVIDIA Titan RTX GPU. The results show characteristic phase contrast and magnitude modulation present in real data. The analytical simulation of cardiac diffusion tensor imaging with bulk-motion phase sensitivity took approximately 10 s per diffusion-weighted image, including preparation and loading steps. The results exhibit the expected alteration of diffusion metrics due to strain.

CONCLUSION

CMRsim is the first simulation framework that allows one to feasibly incorporate complex motion, including turbulent flow, to systematically study advanced CMR acquisition and reconstruction approaches. The open-source package features modularity and transparency, facilitating maintainability and extensibility in support of reproducible research.

Paste, aggregate, or air? That is the question

The Ambassador Bridge between Detroit, Michigan, and Windsor, Ontario, has served for almost 100 years as North America's busiest international border crossing. But in 2025, the Ambassador will be replaced by the new Gordie Howe International Bridge. The Gordie Howe is a cable-stayed bridge, with two massive 220 m tall concrete piers on opposite banks of the St. Claire River, a single clear span of 853 m, and 42 m of clearance over this busy waterway. To ensure durability in this harsh freeze-thaw environment, air-entrained concrete is specified throughout. And, to ensure the quality of air entrainment, the ASTM C 457 Procedure C, Contrast Enhanced Method is employed. While a similar automated microscopic approach has been in use for well over a decade according to EN 480-11 Determination of air void characteristics in hardened concrete, this is the first large-scale application of automated air void assessment in North American infrastructure. According to the ASTM Procedure C, the air void characteristics are determined through digital image processing, while the paste content may be determined by either mix design parameters, manual point count, or 'other means'. Of these three options, point counting is used for Gordie Howe; but in parallel, during each point count, the digital image coordinates and phase identifications for each evaluated stop are recorded. This allows for training of a neural network, for automated determination of paste content, as demonstrated here.

Ontological approach for competency-based curriculum analysis

This article is dedicated to the development of a model for competencies within an educational program and its implementation through the use of semantic technologies. The model proposed by the authors is distinctive in that competencies are organized into a hierarchical data structure with arbitrary levels of nesting. Furthermore, the article presents an original solution for modelling the input requirements for studying a course, which is defined in the form of dependencies between the competencies generated by the course and the competencies of other courses. The outcome of this work is an ontological model of a competency-based curriculum, for which the authors have developed and implemented algorithms for data addition and retrieval, as well as for analyzing the consistency of the curriculum in terms of the input requirements for studying a discipline and the learning outcomes from previous periods. The findings presented in the article will prove to be valuable in the development of educational process management information systems and educational program constructors. They will also be instrumental in aligning diverse educational programs within the context of academic mobility.

Python workflow for the selection and identification of marker peptides-proof-of-principle study with heated milk

The analysis of almost holistic food profiles has developed considerably over the last years. This has also led to larger amounts of data and the ability to obtain more information about health-beneficial and adverse constituents in food than ever before. Especially in the field of proteomics, software is used for evaluation, and these do not provide specific approaches for unique monitoring questions. An additional and more comprehensive way of evaluation can be done with the programming language Python. It offers broad possibilities by a large ecosystem for mass spectrometric data analysis, but needs to be tailored for specific sets of features, the research questions behind. It also offers the applicability of various machine-learning approaches. The aim of the present study was to develop an algorithm for selecting and identifying potential marker peptides from mass spectrometric data. The workflow is divided into three steps: (I) feature engineering, (II) chemometric data analysis, and (III) feature identification. The first step is the transformation of the mass spectrometric data into a structure, which enables the application of existing data analysis packages in Python. The second step is the data analysis for selecting single features. These features are further processed in the third step, which is the feature identification. The data used exemplarily in this proof-of-principle approach was from a study on the influence of a heat treatment on the milk proteome/peptidome.

Unlocking the surface chemistry of ionic minerals: a high-throughput pipeline for modeling realistic interfaces

A systematic procedure is introduced for modeling charge-neutral non-polar surfaces of ionic minerals containing polyatomic anions. By integrating distance- and charge-based clustering to identify chemical species within the mineral bulk, our pipeline, PolyCleaver, renders a variety of theoretically viable surface terminations. As a demonstrative example, this approach was applied to forsterite (Mg2SiO4), unveiling a rich interface landscape based on interactions with formaldehyde, a relevant multifaceted molecule, and more particularly in prebiotic chemistry. This high-throughput method, going beyond techniques traditionally applied in the modeling of minerals, offers new insights into the potential catalytic properties of diverse surfaces, enabling a broader exploration of synthetic pathways in complex mineral systems.

py_ped_sim - A flexible forward genetic simulator for complex family pedigree analysis

Background

Large-scale family pedigrees are commonly used across medical, evolutionary, and forensic genetics. These pedigrees are tools for identifying genetic disorders, tracking evolutionary patterns, and establishing familial relationships via forensic genetic identification. However, there is a lack of software to accurately simulate different pedigree structures along with genomes corresponding to those individuals in a family pedigree. This limits simulation-based evaluations of methods that use pedigrees.

Results

We have developed a python command-line-based tool called py_ped_sim that facilitates the simulation of pedigree structures and the genomes of individuals in a pedigree. py_ped_sim represents pedigrees as directed acyclic graphs, enabling conversion between standard pedigree formats and integration with the forward population genetic simulator, SLiM. Notably, py_ped_sim allows the simulation of varying numbers of offspring for a set of parents, with the capacity to shift the distribution of sibship sizes over generations. We additionally add simulations for events of misattributed paternity, which offers a way to simulate half-sibling relationships. We validated the accuracy of our software by simulating genomes onto diverse family pedigree structures, showing that the estimated kinship coefficients closely approximated expected values.

Conclusions

py_ped_sim is a user-friendly and open-source solution for simulating pedigree structures and conducting pedigree genome simulations. It empowers medical, forensic, and evolutionary genetics researchers to gain deeper insights into the dynamics of genetic inheritance and relatedness within families.

SEAOP: a statistical ensemble approach for outlier detection in quantitative proteomics data

Quality control in quantitative proteomics is a persistent challenge, particularly in identifying and managing outliers. Unsupervised learning models, which rely on data structure rather than predefined labels, offer potential solutions. However, without clear labels, their effectiveness might be compromised. Single models are susceptible to the randomness of parameters and initialization, which can result in a high rate of false positives. Ensemble models, on the other hand, have shown capabilities in effectively mitigating the impacts of such randomness and assisting in accurately detecting true outliers. Therefore, we introduced SEAOP, a Python toolbox that utilizes an ensemble mechanism by integrating multi-round data management and a statistics-based decision pipeline with multiple models. Specifically, SEAOP uses multi-round resampling to create diverse sub-data spaces and employs outlier detection methods to identify candidate outliers in each space. Candidates are then aggregated as confirmed outliers via a chi-square test, adhering to a 95% confidence level, to ensure the precision of the unsupervised approaches. Additionally, SEAOP introduces a visualization strategy, specifically designed to intuitively and effectively display the distribution of both outlier and non-outlier samples. Optimal hyperparameter models of SEAOP for outlier detection were identified by using a gradient-simulated standard dataset and Mann-Kendall trend test. The performance of the SEAOP toolbox was evaluated using three experimental datasets, confirming its reliability and accuracy in handling quantitative proteomics.

suMRak: a multi-tool solution for preclinical brain MRI data analysis

Introduction

Magnetic resonance imaging (MRI) is invaluable for understanding brain disorders, but data complexity poses a challenge in experimental research. In this study, we introduce suMRak, a MATLAB application designed for efficient preclinical brain MRI analysis. SuMRak integrates brain segmentation, volumetry, image registration, and parameter map generation into a unified interface, thereby reducing the number of separate tools that researchers may require for straightforward data handling.

Methods and implementation

All functionalities of suMRak are implemented using the MATLAB App Designer and the MATLAB-integrated Python engine. A total of six helper applications were developed alongside the main suMRak interface to allow for a cohesive and streamlined workflow. The brain segmentation strategy was validated by comparing suMRak against manual segmentation and ITK-SNAP, a popular open-source application for biomedical image segmentation.

Results

When compared with the manual segmentation of coronal mouse brain slices, suMRak achieved a high Sørensen-Dice similarity coefficient (0.98 ± 0.01), approaching manual accuracy. Additionally, suMRak exhibited significant improvement (p = 0.03) when compared to ITK-SNAP, particularly for caudally located brain slices. Furthermore, suMRak was capable of effectively analyzing preclinical MRI data obtained in our own studies. Most notably, the results of brain perfusion map registration to T2-weighted images were shown, improving the topographic connection to anatomical areas and enabling further data analysis to better account for the inherent spatial distortions of echoplanar imaging.

Discussion

SuMRak offers efficient MRI data processing of preclinical brain images, enabling researchers' consistency and precision. Notably, the accelerated brain segmentation, achieved through K-means clustering and morphological operations, significantly reduces processing time and allows for easier handling of larger datasets.

Bibliometric insights into the evolution of uranium contamination reduction research topics: Focus on microbial reduction of uranium

Confronting the threat of environment uranium pollution, decades of research have yielded advanced and significant findings in uranium bioremediation, resulting in the accumulation of tremendous amount of high-quality literature. In this study, we analyzed over 10,000 uranium reduction-related papers published from 1990 to the present in the Web of Science based on bibliometrics, and revealed some critical information on knowledge structure, thematic evolution and additional attention. Methods including contribution comparison, co-occurrence and temporal evolution analysis are applied. The results of the distribution and impact analysis of authors, sources, and journals indicated that the United States is a leader in this field of research and China is on the rise. The top keywords remained stable, primarily focused on chemicals (uranium, iron, plutonium, nitrat, carbon), characters (divers, surfac, speciat), and microbiology (microbial commun, cytochrome, extracellular polymeric subst). Keywords related to new strains, reduction mechanisms and product characteristics demonstrated the strongest uptrend, while some keywords related to mechanism and performance were clearly emerging in the past 5 years. Furthermore, the evolution of the thematic progression can be categorized into three stages, commencing with the discovery of the enzymatic reduction of hexavalent uranium to tetravalent uranium, developing in the groundwater remediation process at uranium-contaminated sites, and delving into the research on microbial reduction mechanisms of uranium. For future research, enhancing the understanding of mechanisms, improving uranium removal performance, and exploring practical applications can be considered. This study provides unique insights into microbial uranium reduction research, providing valuable references for related studies in this field.

Machine learning regression algorithms to predict emissions from steam boilers

Currently, the modeling of complex chemical-physical processes is drastically influencing industrial development. Therefore, the analysis and study of the combustion process of the boilers using machine learning (ML) techniques are vital to increase the efficiency with which this equipment operates and reduce the pollution load they contribute to the environment. This work aims to predict the emissions of CO, CO2, NOx, and the temperature of the exhaust gases of industrial boilers from real data. Different ML algorithms for regression analysis are discussed. The following are input variables: ambient temperature, working pressure, steam production, and the type of fuel used in around 20 industrial boilers. Each boiler's emission data was collected using a TESTO 350 Combustion Gas Analyzer. The modeling, with a machine learning approach using the Gradient Boosting Regression algorithm, showed better performance in the predictions made on the test data, outperforming all other models studied. It was achieved with predicted values showing a mean absolute error of 0.51 and a coefficient of determination of 99.80%. Different regression models (DNN, MLR, RFR, GBR) were compared to select the most optimal. Compared to models based on Linear Regression, the DNN model has better prediction performance. The proposed model provides a new method to predict CO2, CO, NOx emissions, and exhaust gas outlet temperature.

Demystifying dimensionality reduction techniques in the 'omics' era: A practical approach for biological science students

Dimensionality reduction techniques are essential in analyzing large 'omics' datasets in biochemistry and molecular biology. Principal component analysis, t-distributed stochastic neighbor embedding, and uniform manifold approximation and projection are commonly used for data visualization. However, these methods can be challenging for students without a strong mathematical background. In this study, intuitive examples were created using COVID-19 data to help students understand the core concepts behind these techniques. In a 4-h practical session, we used these examples to demonstrate dimensionality reduction techniques to 15 postgraduate students from biomedical backgrounds. Using Python and Jupyter notebooks, our goal was to demystify these methods, typically treated as "black boxes", and empower students to generate and interpret their own results. To assess the impact of our approach, we conducted an anonymous survey. The majority of the students agreed that using computers enriched their learning experience (67%) and that Jupyter notebooks were a valuable part of the class (66%). Additionally, 60% of the students reported increased interest in Python, and 40% gained both interest and a better understanding of dimensionality reduction methods. Despite the short duration of the course, 40% of the students reported acquiring research skills necessary in the field. While further analysis of the learning impacts of this approach is needed, we believe that sharing the examples we generated can provide valuable resources for others to use in interactive teaching environments. These examples highlight advantages and limitations of the major dimensionality reduction methods used in modern bioinformatics analysis in an easy-to-understand way.

A spatio-temporal image analysis for growth of indeterminate pulmonary nodules detected by CT scan

The objective is to evaluate the performance of computational image classification for indeterminate pulmonary nodules (IPN) chronologically detected by CT scan. Total 483 patients with 670 abnormal pulmonary nodules, who were taken chest thin-section CT (TSCT) images at least twice and resected as suspicious nodules in our hospital, were enrolled in this study. Nodular regions from the initial and the latest TSCT images were cut manually for each case, and approached by Python development environment, using the open-source cv2 library, to measure the nodular change rate (NCR). These NCRs were statistically compared with clinico-pathological factors, and then, this discriminator was evaluated for clinical performance. NCR showed significant differences among the nodular consistencies. In terms of histological subtypes, NCR of invasive adenocarcinoma (ADC) were significantly distinguishable from other lesions, but not from minimally invasive ADC. Only for cancers, NCR was significantly associated with loco-regional invasivity, p53-immunoreactivity, and Ki67-immunoreactivity. Regarding Epidermal Growth Factor Receptor gene mutation of ADC-related nodules, NCR showed a significant negative correlation. On staging of lung cancer cases, NCR was significantly increased with progression from pTis-stage 0 up to pT1b-stage IA2. For clinical shared decision-making (SDM) whether urgent resection or watchful-waiting, receiver operating characteristic (ROC) analysis showed that area under the ROC curve was 0.686. For small-sized IPN detected by CT scan, this approach shows promise as a potential navigator to improve work-up for life-threatening cancer screening and assist SDM before surgery.

Finite element-based optimization procedure for an irregular domain with unstructured mesh

At present, structural optimization is a highly demanding area of research in engineering. Engineers aim to minimize material in a body while maintaining its usability and safety at the same time. Developing a user-friendly program to optimize a structure using the finite element method (FEM) is the goal of the current study. With the advent of additive manufacturing, the production of complex-shaped designs is showing promise. A detailed optimization algorithm based on solid isotropic material with penalization (SIMP) is presented in this paper. UnTop2D: An object-oriented Python program with a graphical user interface (GUI) has been developed, which can be applied to structures with both structured and unstructured meshes. The mesh is not required to be topologically ball and can be imported from professional meshing software. Any selected element can be frozen to prevent its removal during optimization, and wall elements can also be frozen for real-world scenarios. The optimized structure can be exported as an Abaqus input file for structural analysis and STL file for 3D printing. This paper presents several examples to demonstrate the effectiveness of the proposed procedure.

Scholar Metrics Scraper (SMS): automated retrieval of citation and author data

Academic departments, research clusters and evaluators analyze author and citation data to measure research impact and to support strategic planning. We created Scholar Metrics Scraper (SMS) to automate the retrieval of bibliometric data for a group of researchers. The project contains Jupyter notebooks that take a list of researchers as an input and exports a CSV file of citation metrics from Google Scholar (GS) to visualize the group's impact and collaboration. A series of graph outputs are also available. SMS is an open solution for automating the retrieval and visualization of citation data.

ReTimeML: a retention time predictor that supports the LC-MS/MS analysis of sphingolipids

The analysis of ceramide (Cer) and sphingomyelin (SM) lipid species using liquid chromatography-tandem mass spectrometry (LC-MS/MS) continues to present challenges as their precursor mass and fragmentation can correspond to multiple molecular arrangements. To address this constraint, we developed ReTimeML, a freeware that automates the expected retention times (RTs) for Cer and SM lipid profiles from complex chromatograms. ReTimeML works on the principle that LC-MS/MS experiments have pre-determined RTs from internal standards, calibrators or quality controls used throughout the analysis. Employed as reference RTs, ReTimeML subsequently extrapolates the RTs of unknowns using its machine-learned regression library of mass-to-charge (m/z) versus RT profiles, which does not require model retraining for adaptability on different LC-MS/MS pipelines. We validated ReTimeML RT estimations for various Cer and SM structures across different biologicals, tissues and LC-MS/MS setups, exhibiting a mean variance between 0.23 and 2.43% compared to user annotations. ReTimeML also aided the disambiguation of SM identities from isobar distributions in paired serum-cerebrospinal fluid from healthy volunteers, allowing us to identify a series of non-canonical SMs associated between the two biofluids comprised of a polyunsaturated structure that confers increased stability against catabolic clearance.

Digital measuring the ocular morphological parameters of guinea pig eye in vivo with Python

AIM

To quantitatively measure ocular morphological parameters of guinea pig with Python technology.

METHODS

Thirty-six eyeballs of eighteen 3-week-old guinea pigs were measured with keratometer and photographed to obtain the horizontal, coronal, and sagittal planes respectively. The corresponding photo pixels-actual length ratio was acquired by a proportional scale. The edge coordinates were identified artificially by ginput function. Circle and conic curve fitting were applied to fit the contour of the eyeball in the sagittal, coronal and horizontal view. The curvature, curvature radius, eccentricity, tilt angle, corneal diameter, and binocular separation angle were calculated according to the geometric principles. Next, the eyeballs were removed, canny edge detection was applied to identify the contour of eyeball in vitro. The results were compared between in vivo and in vitro.

RESULTS

Regarding the corneal curvature and curvature radius on the horizontal and sagittal planes, no significant differences were observed among results in vivo, in vitro, and the keratometer. The horizontal and vertical binocular separation angles were 130.6°±6.39° and 129.8°±9.58° respectively. For the corneal curvature radius and eccentricity in vivo, significant differences were observed between horizontal and vertical planes.

CONCLUSION

The Graphical interface window of Python makes up the deficiency of edge detection, which requires too much definition in Matlab. There are significant differences between guinea pig and human beings, such as exotropic eye position, oblique oval eyeball, and obvious discrepancy of binoculus. This study helps evaluate objectively the ocular morphological parameters of small experimental animals in emmetropization research.

An Ogden hyperelastic 3D micromechanical model to depict Poynting effect in brain white matter

Shear and torsional load on soft solids such as brain white matter purportedly exhibits the Poynting Effect. It is a typical nonlinear phenomenon associated with soft materials whereby they tend to elongate (positive Poynting effect) or contract (negative Poynting effect) in a direction perpendicular to the shearing or twisting plane. In this research, a novel 3D micromechanical Finite Element Model (FEM) has been formulated to describe the Poynting effect in bi-phasic modeled brain white matter (BWM) representative volume element (RVE) with axons tracts embedded in surrounding extracellular matrix (ECM) for simulating brain matter's response to pure and simple shear. In the presented BWM 3D FEM, nonlinear Ogden hyper-elastic material model is deployed to interpret axons and ECM material phases. The modeled bi-phasic RVEs have axons tied to the surrounding ECM. In this proof-of-concept (POC) FEM, three simple shear loading configurations and a pure shear case were analyzed. Root mean square deviation (RMSD) was calculated for stress and deformation response plots to understand the effect of axon-ECM orientations and loading conditions on the degree of Poynting behavior. Variations in normal stresses (S11, S22, or S33) perpendicular to the shear plane underscored the significance of axonal fiber-matrix interactions. From the simulated ensemble of cases, a transitional dominance trend was noticed, as simple sheared axons showed pronounced Poynting behavior, but shear deformation build-up in the purely sheared brain model exhibited the highest Poynting behavior at higher strain % limits. At lower strain limits, simple shear imparted across and perpendicular to axonal tract directions emerged as the dominant Poynting effect configurations. At high strains, the stress-strain% plots manifested mild strain stiffening effects and bending stresses in purely sheared axons, substantiated the strong non-linearity in brain tissues' response.

thebeat: A Python package for working with rhythms and other temporal sequences

thebeat is a Python package for working with temporal sequences and rhythms in the behavioral and cognitive sciences, as well as in bioacoustics. It provides functionality for creating experimental stimuli, and for visualizing and analyzing temporal data. Sequences, sounds, and experimental trials can be generated using single lines of code. thebeat contains functions for calculating common rhythmic measures, such as interval ratios, and for producing plots, such as circular histograms. thebeat saves researchers time when creating experiments, and provides the first steps in collecting widely accepted methods for use in timing research. thebeat is an open-source, on-going, and collaborative project, and can be extended for use in specialized subfields. thebeat integrates easily with the existing Python ecosystem, allowing one to combine our tested code with custom-made scripts. The package was specifically designed to be useful for both skilled and novice programmers. thebeat provides a foundation for working with temporal sequences onto which additional functionality can be built. This combination of specificity and plasticity should facilitate research in multiple research contexts and fields of study.

IFM calculator: An algorithm for interfragmentary motion calculation in finite element analysis

BACKGROUND

Interfragmentary motion (IFM) is a complex state that significantly impacts the healing process of fractures following implant placement. It is crucial to fully consider the IFM state after implantation in the design and biomechanical testing of implants. However, current finite element analysis software lacks direct tools for calculating IFM, and existing IFM tools do not offer a comprehensive solution in terms of accuracy, functionality, and visualization.

METHODS

In our study, we developed a Python-based algorithm for calculating IFM that addresses limitations. Our algorithm automatically calculated IFM distances, sliding distances, gaps, as well as the angles and rotation of the two fracture surfaces using all nodes on both sides of the fracture ends. Researchers could input data and selected desired parameters in the interface. The algorithm then performed the necessary calculations and presented the results in a clear and concise manner. The algorithm also provided comprehensive data export capabilities, allowing researchers to customize analyses based on specific needs.To provide a more intuitive demonstration of the calculation process and usage of IFM-Cal, we conducted simulations in Ansys using two rectangular blocks to compare the accuracy and function of three different methods (Point based method, contact tool and IFM-Cal).

RESULTS

The point-based method and the contact tool could not accurately calculate IFA, while IFM-Cal could provide a comprehensive evaluation of IFA. In simulation 1, the IFM distances calculated using the point sampling method, contact tool, and IFM-Cal were 2.00 mm, 3.15 mm, and 2.00 mm, respectively. In simulation 2, both the point sampling method and contact tool failed to calculate the interfragmentary angle (IFA), while the IFM-Cal algorithm estimated an angle of -7.87°, which had a small error compared to the ground-truth value of 7.9°.

CONCLUSION

We have developed an algorithm for computing IFM which can be utilized in finite element analysis and biomechanical experiments. By conducting comparative simulations with other existing algorithms, we have demonstrated the superior accuracy and expanded evaluation capabilities of our algorithm.

Tutorial: Analysis of central and peripheral motor unit properties from decomposed High-Density surface EMG signals with openhdemg

High-Density surface Electromyography (HD-sEMG) is the most established technique for the non-invasive analysis of single motor unit (MU) activity in humans. It provides the possibility to study the central properties (e.g., discharge rate) of large populations of MUs by analysis of their firing pattern. Additionally, by spike-triggered averaging, peripheral properties such as MUs conduction velocity can be estimated over adjacent regions of the muscles and single MUs can be tracked across different recording sessions. In this tutorial, we guide the reader through the investigation of MUs properties from decomposed HD-sEMG recordings by providing both the theoretical knowledge and practical tools necessary to perform the analyses. The practical application of this tutorial is based on openhdemg, a free and open-source community-based framework for the automated analysis of MUs properties built on Python 3 and composed of different modules for HD-sEMG data handling, visualisation, editing, and analysis. openhdemg is interfaceable with most of the available recording software, equipment or decomposition techniques, and all the built-in functions are easily adaptable to different experimental needs. The framework also includes a graphical user interface which enables users with limited coding skills to perform a robust and reliable analysis of MUs properties without coding.

Indonesian presidential election sentiment: Dataset of response public before 2024

Indonesia is one of the countries that is currently entering the political year for the election of President, Regional Heads, and Members of the Legislative in 2024. This has become a hot topic on social media, especially about the Presidential Election. Twitter is one of the platforms with the largest users in Indonesia. It is interesting to see the alignment of Twitter users towards presidential candidates who already have a carrying party, namely Ganjar Pranowo, Prabowo Subianto, and Anies Baswedan based on a sentiment analysis approach. User feedback data about Indonesian Presidential candidates are obtained from the Twitter platform using Twitter API with Python programming language. The data obtained was 30,000 data with each candidate as many as 10,000 data. Data is pulled in April 2023 with specific keywords. The time for data withdrawal is chosen based on the announcement of Presidential Candidates carried by political parties before the schedule for determining or campaigning for Presidential candidates. Current data can potentially be used again as a comparison of analysis of presidential candidates on campaign time spans and after campaigns or actual calculation results. The data that can be accessed is in CSV format and has gone through several stages such as labelling using Language experts, removing spam Tweets & empty cells and preprocessing.

Convolutional Neural Network (CNN) Model for the Classification of Varieties of Date Palm Fruits (Phoenix dactylifera L.)

The popularity and demand for high-quality date palm fruits (Phoenix dactylifera L.) have been growing, and their quality largely depends on the type of handling, storage, and processing methods. The current methods of geometric evaluation and classification of date palm fruits are characterised by high labour intensity and are usually performed mechanically, which may cause additional damage and reduce the quality and value of the product. Therefore, non-contact methods are being sought based on image analysis, with digital solutions controlling the evaluation and classification processes. The main objective of this paper is to develop an automatic classification model for varieties of date palm fruits using a convolutional neural network (CNN) based on two fundamental criteria, i.e., colour difference and evaluation of geometric parameters of dates. A CNN with a fixed architecture was built, marked as DateNET, consisting of a system of five alternating Conv2D, MaxPooling2D, and Dropout classes. The validation accuracy of the model presented in this study depended on the selection of classification criteria. It was 85.24% for fruit colour-based classification and 87.62% for the geometric parameters only; however, it increased considerably to 93.41% when both the colour and geometry of dates were considered.

PyRINEX: a new multi-purpose Python package for GNSS RINEX data

Since the first receiver independent exchange format (RINEX) version was released in 1989, it has gone through several versions, making the existing software, such as TEQC, incompatible with certain later versions. This study proposes a new Python package named PyRINEX, which is developed to batch process the most generally used versions of RINEX files, namely 2.0 and 3.0. The proposed package can be used to manage and edit numerous RINEX files as well as perform a data quality check function. PyRINEX can be easily imported into any Python IDE similar to any other open-source Python package, it also makes secondary development easy for users.

DARTS: an open-source Python pipeline for Ca2+ microdomain analysis in live cell imaging data

Ca2+ microdomains play a key role in intracellular signaling processes. For instance, they mediate the activation of T cells and, thus, the initial adaptive immune system. They are, however, also of utmost importance for activation of other cells, and a detailed understanding of the dynamics of these spatially localized Ca2+ signals is crucial for a better understanding of the underlying signaling processes. A typical approach to analyze Ca2+ microdomain dynamics is live cell fluorescence microscopy imaging. Experiments usually involve imaging a larger number of cells of different groups (for instance, wild type and knockout cells), followed by a time consuming image and data analysis. With DARTS, we present a modular Python pipeline for efficient Ca2+ microdomain analysis in live cell imaging data. DARTS (Deconvolution, Analysis, Registration, Tracking, and Shape normalization) provides state-of-the-art image postprocessing options like deep learning-based cell detection and tracking, spatio-temporal image deconvolution, and bleaching correction. An integrated automated Ca2+ microdomain detection offers direct access to global statistics like the number of microdomains for cell groups, corresponding signal intensity levels, and the temporal evolution of the measures. With a focus on bead stimulation experiments, DARTS provides a so-called dartboard projection analysis and visualization approach. A dartboard projection covers spatio-temporal normalization of the bead contact areas and cell shape normalization onto a circular template that enables aggregation of the spatiotemporal information of the microdomain detection results for the individual cells of the cell groups of interest. The dartboard visualization allows intuitive interpretation of the spatio-temporal microdomain dynamics at the group level. The application of DARTS is illustrated by three use cases in the context of the formation of initial Ca2+ microdomains after cell stimulation. DARTS is provided as an open-source solution and will be continuously extended upon the feedback of the community. Code available at: 10.5281/zenodo.10459243.

The probability of edge existence due to node degree: a baseline for network-based predictions

Important tasks in biomedical discovery such as predicting gene functions, gene-disease associations, and drug repurposing opportunities are often framed as network edge prediction. The number of edges connecting to a node, termed degree, can vary greatly across nodes in real biomedical networks, and the distribution of degrees varies between networks. If degree strongly influences edge prediction, then imbalance or bias in the distribution of degrees could lead to nonspecific or misleading predictions. We introduce a network permutation framework to quantify the effects of node degree on edge prediction. Our framework decomposes performance into the proportions attributable to degree and the network's specific connections using network permutation to generate features that depend only on degree. We discover that performance attributable to factors other than degree is often only a small portion of overall performance. Researchers seeking to predict new or missing edges in biological networks should use our permutation approach to obtain a baseline for performance that may be nonspecific because of degree. We released our methods as an open-source Python package (https://github.com/hetio/xswap/).

Computational reproducibility of Jupyter notebooks from biomedical publications

BACKGROUND

Jupyter notebooks facilitate the bundling of executable code with its documentation and output in one interactive environment, and they represent a popular mechanism to document and share computational workflows, including for research publications. The reproducibility of computational aspects of research is a key component of scientific reproducibility but has not yet been assessed at scale for Jupyter notebooks associated with biomedical publications.

APPROACH

We address computational reproducibility at 2 levels: (i) using fully automated workflows, we analyzed the computational reproducibility of Jupyter notebooks associated with publications indexed in the biomedical literature repository PubMed Central. We identified such notebooks by mining the article's full text, trying to locate them on GitHub, and attempting to rerun them in an environment as close to the original as possible. We documented reproduction success and exceptions and explored relationships between notebook reproducibility and variables related to the notebooks or publications. (ii) This study represents a reproducibility attempt in and of itself, using essentially the same methodology twice on PubMed Central over the course of 2 years, during which the corpus of Jupyter notebooks from articles indexed in PubMed Central has grown in a highly dynamic fashion.

RESULTS

Out of 27,271 Jupyter notebooks from 2,660 GitHub repositories associated with 3,467 publications, 22,578 notebooks were written in Python, including 15,817 that had their dependencies declared in standard requirement files and that we attempted to rerun automatically. For 10,388 of these, all declared dependencies could be installed successfully, and we reran them to assess reproducibility. Of these, 1,203 notebooks ran through without any errors, including 879 that produced results identical to those reported in the original notebook and 324 for which our results differed from the originally reported ones. Running the other notebooks resulted in exceptions.

CONCLUSIONS

We zoom in on common problems and practices, highlight trends, and discuss potential improvements to Jupyter-related workflows associated with biomedical publications.

MetaBCI: An open-source platform for brain-computer interfaces

BACKGROUND

Recently, brain-computer interfaces (BCIs) have attracted worldwide attention for their great potential in clinical and real-life applications. To implement a complete BCI system, one must set up several links to translate the brain intent into computer commands. However, there is not an open-source software platform that can cover all links of the BCI chain.

METHOD

This study developed a one-stop open-source BCI software, namely MetaBCI, to facilitate the construction of a BCI system. MetaBCI is written in Python, and has the functions of stimulus presentation (Brainstim), data loading and processing (Brainda), and online information flow (Brainflow). This paper introduces the detailed information of MetaBCI and presents four typical application cases.

RESULTS

The results showed that MetaBCI was an extensible and feature-rich software platform for BCI research and application, which could effectively encode, decode, and feedback brain activities.

CONCLUSIONS

MetaBCI can greatly lower the BCI's technical threshold for BCI beginners and can save time and cost to build up a practical BCI system. The source code is available at https://github.com/TBC-TJU/MetaBCI, expecting new contributions from the BCI community.

LAMPrimers iQ: New primer design software for loop-mediated isothermal amplification (LAMP)

Nucleic acids amplification is a widely used technique utilized for different manipulations with DNA and RNA. Although, polymerase chain reaction (PCR) remains the most popular amplification method, isothermal approaches are gained more attention last decades. Among these, loop-mediated isothermal amplification (LAMP) became an excellent alternative to PCR. LAMP requires an increased number of primers and, therefore, is considered a highly specific amplification reaction compared to PCR. LAMP primers design is still a non-trivial task, and all niceties should be taken into account during their selection. Here, we report on a new program called LAMPrimers iQ destined for high-quality LAMP primers design. LAMPrimers iQ is based on an original algorithm considering rigorous criteria for primers selection. Unlike alternative programs, LAMPrimers iQ can process long DNA or RNA sequences, and completely avoid primers that can form homo- and heterodimers. The quality of the primers designed was checked using SARS-CoV-2 coronavirus RNA as a model target. It was shown that primers selected with LAMPrimers iQ provide higher specificity and reliable detection of viral RNA compared to those obtained by alternative programs. The program is available at https://github.com/Restily/LAMPrimers-iQ.

Big Data Analysis in Computational Biology and Bioinformatics

Advancements in high-throughput technologies, genomics, transcriptomics, and metabolomics play an important role in obtaining biological information about living organisms. The field of computational biology and bioinformatics has experienced significant growth with the advent of high-throughput sequencing technologies and other high-throughput techniques. The resulting large amounts of data present both opportunities and challenges for data analysis. Big data analysis has become essential for extracting meaningful insights from the massive amount of data. In this chapter, we provide an overview of the current status of big data analysis in computational biology and bioinformatics. We discuss the various aspects of big data analysis, including data acquisition, storage, processing, and analysis. We also highlight some of the challenges and opportunities of big data analysis in this area of research. Despite the challenges, big data analysis presents significant opportunities like development of efficient and fast computing algorithms for advancing our understanding of biological processes, identifying novel biomarkers for breeding research and developments, predicting disease, and identifying potential drug targets for drug development programs.

An Approach to Track and Analyze the Trend of Antimicrobial Resistance Using Python: A Pilot Study for Anand, Gujarat, India-May 2022-August 2023

The present work deals with the analysis and monitoring of bacterial resistance in using Python for the state of Gujarat, India, where occurrences of drug-resistant bacteria are prevalent. This will provide an insight into the portfolio of drug-resistant bacteria reported, which can be used to track resistance behavior and to suggest a treatment regime for the particular bacteria. The present analysis has been done using Python on Jupyter Notebook as the integrated development environment and its data analysis libraries such as Pandas, Seaborn, and Matplotlib. The data have been loaded from excel file using Pandas and cleaned to transform features into required format. Seaborn and Matplotlib have been used to create data visualizations and represent the data inexplicable manner using graphs, plots, and tables. This program can be used to study disaster epidemiology, tracking, analyzing, and surveillance of antimicrobial resistance with a proper system integration approach.

Control Software Design for a Multisensing Multicellular Microscale Gas Chromatography System

Microscale gas chromatography (μGC) systems are miniaturized instruments that typically incorporate one or several microfabricated fluidic elements; such systems are generally well suited for the automated sampling and analysis of gas-phase chemicals. Advanced μGC systems may incorporate more than 15 elements and operate these elements in different coordinated sequences to execute complex operations. In particular, the control software must manage the sampling and analysis operations of the μGC system in a time-sensitive manner; while operating multiple control loops, it must also manage error conditions, data acquisition, and user interactions when necessary. To address these challenges, this work describes the investigation of multithreaded control software and its evaluation with a representative μGC system. The μGC system is based on a progressive cellular architecture that uses multiple μGC cells to efficiently broaden the range of chemical analytes, with each cell incorporating multiple detectors. Implemented in Python language version 3.7.3 and executed by an embedded single-board computer, the control software enables the concurrent control of heaters, pumps, and valves while also gathering data from thermistors, pressure sensors, capacitive detectors, and photoionization detectors. A graphical user interface (UI) that operates on a laptop provides visualization of control parameters in real time. In experimental evaluations, the control software provided successful operation and readout for all the components, including eight sets of thermistors and heaters that form temperature feedback loops, two sets of pressure sensors and tunable gas pumps that form pressure head feedback loops, six capacitive detectors, three photoionization detectors, six valves, and an additional fixed-flow gas pump. A typical run analyzing 18 chemicals is presented. Although the operating system does not guarantee real-time operation, the relative standard deviations of the control loop timings were <0.5%. The control software successfully supported >1000 μGC runs that analyzed various chemical mixtures.

Water Effective Diffusion Coefficient in Dairy Powder Calculated by Digital Image Processing and through Machine Learning Algorithms of CLSM Micrographs

Rehydration of dairy powders is a complex and essential process. A relatively new quantitative mechanism for monitoring powders' rehydration process uses the effective diffusion coefficient. This research focused on modifying a previously used labor-intensive method that will be able to automatically measure the real-time water diffusion coefficient in dairy powders based on confocal microscopy techniques. Furthermore, morphological characteristics and local hydration of individual particles were identified using an imaging analysis procedure written in Matlab©-R2023b and image analysis through machine learning algorithms written in Python™-3.11. The first model includes segmentation into binary images and labeling particles during water diffusion. The second model includes the expansion of data set selection, neural network training and particle markup. For both models, the effective diffusion follows Fick's second law for spherical geometry. The effective diffusion coefficient on each particle was computed from the dye intensity during the rehydration process. The results showed that effective diffusion coefficients for water increased linearly with increasing powder particle size and are in agreement with previously used methods. In summary, the models provide two independent machine measurements of effective diffusion coefficient based on the same set of micrographs and may be useful in a wide variety of high-protein powders.

Structural Outlier Detection and Zernike-Canterakis Moments for Molecular Surface Meshes-Fast Implementation in Python

Object retrieval systems measure the degree of similarity of the shape of 3D models. They search for the elements of the 3D model databases that resemble the query model. In structural bioinformatics, the query model is a protein tertiary/quaternary structure and the objective is to find similarly shaped molecules in the Protein Data Bank. With the ever-growing size of the PDB, a direct atomic coordinate comparison with all its members is impractical. To overcome this problem, the shape of the molecules can be encoded by fixed-length feature vectors. The distance of a protein to the entire PDB can be measured in this low-dimensional domain in linear time. The state-of-the-art approaches utilize Zernike-Canterakis moments for the shape encoding and supply the retrieval process with geometric data of the input structures. The BioZernike descriptors are a standard utility of the PDB since 2020. However, when trying to calculate the ZC moments locally, the issue of the deficiency of libraries readily available for use in custom programs (i.e., without relying on external binaries) is encountered, in particular programs written in Python. Here, a fast and well-documented Python implementation of the Pozo-Koehl algorithm is presented. In contrast to the more popular algorithm by Novotni and Klein, which is based on the voxelized volume, the PK algorithm produces ZC moments directly from the triangular surface meshes of 3D models. In particular, it can accept the molecular surfaces of proteins as its input. In the presented PK-Zernike library, owing to Numba's just-in-time compilation, a mesh with 50,000 facets is processed by a single thread in a second at the moment order 20. Since this is the first time the PK algorithm is used in structural bioinformatics, it is employed in a novel, simple, but efficient protein structure retrieval pipeline. The elimination of the outlying chain fragments via a fast PCA-based subroutine improves the discrimination ability, allowing for this pipeline to achieve an 0.961 area under the ROC curve in the BioZernike validation suite (0.997 for the assemblies). The correlation between the results of the proposed approach and of the 3D Surfer program attains values up to 0.99.

Mastering data visualization with Python: practical tips for researchers

Big data have revolutionized the way data are processed and used across all fields. In the past, research was primarily conducted with a focus on hypothesis confirmation using sample data. However, in the era of big data, this has shifted to gaining insights from the collected data. Visualizing vast amounts of data to derive insights is crucial. For instance, leveraging big data for visualization can help identify and predict characteristics and patterns related to various infectious diseases. When data are presented in a visual format, patterns within the data become clear, making it easier to comprehend and provide deeper insights. This study aimed to comprehensively discuss data visualization and the various techniques used in the process. It also sought to enable researchers to directly use Python programs for data visualization. By providing practical visualization exercises on GitHub, this study aimed to facilitate their application in research endeavors.

Mapping the Bibliometrics Landscape of AI in Medicine: Methodological Study

BACKGROUND

Artificial intelligence (AI), conceived in the 1950s, has permeated numerous industries, intensifying in tandem with advancements in computing power. Despite the widespread adoption of AI, its integration into medicine trails other sectors. However, medical AI research has experienced substantial growth, attracting considerable attention from researchers and practitioners.

OBJECTIVE

In the absence of an existing framework, this study aims to outline the current landscape of medical AI research and provide insights into its future developments by examining all AI-related studies within PubMed over the past 2 decades. We also propose potential data acquisition and analysis methods, developed using Python (version 3.11) and to be executed in Spyder IDE (version 5.4.3), for future analogous research.

METHODS

Our dual-pronged approach involved (1) retrieving publication metadata related to AI from PubMed (spanning 2000-2022) via Python, including titles, abstracts, authors, journals, country, and publishing years, followed by keyword frequency analysis and (2) classifying relevant topics using latent Dirichlet allocation, an unsupervised machine learning approach, and defining the research scope of AI in medicine. In the absence of a universal medical AI taxonomy, we used an AI dictionary based on the European Commission Joint Research Centre AI Watch report, which emphasizes 8 domains: reasoning, planning, learning, perception, communication, integration and interaction, service, and AI ethics and philosophy.

RESULTS

From 2000 to 2022, a comprehensive analysis of 307,701 AI-related publications from PubMed highlighted a 36-fold increase. The United States emerged as a clear frontrunner, producing 68,502 of these articles. Despite its substantial contribution in terms of volume, China lagged in terms of citation impact. Diving into specific AI domains, as the Joint Research Centre AI Watch report categorized, the learning domain emerged dominant. Our classification analysis meticulously traced the nuanced research trajectories across each domain, revealing the multifaceted and evolving nature of AI's application in the realm of medicine.

CONCLUSIONS

The research topics have evolved as the volume of AI studies increases annually. Machine learning remains central to medical AI research, with deep learning expected to maintain its fundamental role. Empowered by predictive algorithms, pattern recognition, and imaging analysis capabilities, the future of AI research in medicine is anticipated to concentrate on medical diagnosis, robotic intervention, and disease management. Our topic modeling outcomes provide a clear insight into the focus of AI research in medicine over the past decades and lay the groundwork for predicting future directions. The domains that have attracted considerable research attention, primarily the learning domain, will continue to shape the trajectory of AI in medicine. Given the observed growing interest, the domain of AI ethics and philosophy also stands out as a prospective area of increased focus.

PDBe CCDUtils: an RDKit-based toolkit for handling and analysing small molecules in the Protein Data Bank

While the Protein Data Bank (PDB) contains a wealth of structural information on ligands bound to macromolecules, their analysis can be challenging due to the large amount and diversity of data. Here, we present PDBe CCDUtils, a versatile toolkit for processing and analysing small molecules from the PDB in PDBx/mmCIF format. PDBe CCDUtils provides streamlined access to all the metadata for small molecules in the PDB and offers a set of convenient methods to compute various properties using RDKit, such as 2D depictions, 3D conformers, physicochemical properties, scaffolds, common fragments, and cross-references to small molecule databases using UniChem. The toolkit also provides methods for identifying all the covalently attached chemical components in a macromolecular structure and calculating similarity among small molecules. By providing a broad range of functionality, PDBe CCDUtils caters to the needs of researchers in cheminformatics, structural biology, bioinformatics and computational chemistry.

Development of an Arduino-based, open-control interface for hardware in the loop applications

This article presents a flexible control interface based on low-cost hardware solutions for electric drives which classically come either with a proprietary hardware solution or a high-cost interface solution. The interface presented can be used to connect a standard PC with an electric drive to enable testing simulation and control applications. The control interface is developed based on the open-source Python scripting language and Arduino's open-source and accessible hardware. The new interface communicates with the test stand through its I/O terminals via developed electronic amplifiers and creates a solid base for further development towards more extensive hardware in the loop simulations.

MASSA Algorithm: an automated rational sampling of training and test subsets for QSAR modeling

QSAR models capable of predicting biological, toxicity, and pharmacokinetic properties were widely used to search lead bioactive molecules in chemical databases. The dataset's preparation to build these models has a strong influence on the quality of the generated models, and sampling requires that the original dataset be divided into training (for model training) and test (for statistical evaluation) sets. This sampling can be done randomly or rationally, but the rational division is superior. In this paper, we present MASSA, a Python tool that can be used to automatically sample datasets by exploring the biological, physicochemical, and structural spaces of molecules using PCA, HCA, and K-modes. The proposed algorithm is very useful when the variables used for QSAR are not available or to construct multiple QSAR models with the same training and test sets, producing models with lower variability and better values for validation metrics. These results were obtained even when the descriptors used in the QSAR/QSPR were different from those used in the separation of training and test sets, indicating that this tool can be used to build models for more than one QSAR/QSPR technique. Finally, this tool also generates useful graphical representations that can provide insights into the data.

GUI development for SSDL to calibrate photon measuring equipment

Current legislation mandates the inspection and calibration of operational survey radiation monitoring instruments used in nuclear medicine, radiotherapy departments, and other fields utilizing ionizing radiation sources. To comply with national and international radiation protection standards, Morocco's National Secondary Standard Dosimetry Laboratory provides reliable calibration results with high accuracy and covers various measurement ranges using attenuators provided by the automated Gamma G10 irradiator or validated beam qualities produced by the X-ray irradiator type X80-320 kV. This study aims to develop a digital graphical user interface using Python programming language, designed for calibrating radiation protection measuring instruments . The interface is intended to facilitate all operations and calculations related to determining calibration factors and measurement uncertainties in accordance with the ISO 4037 standard, ensuring minimal processing time and minimizing potential error sources . The interface enables calculations to be recorded, as well as the establishment and electronic archiving of the calibration certificate and the report in PDF format using the Hypertext Preprocessor FPDF library (PHP FPDF). With the development of this interface, multiple instruments can be processed per day with high accuracy, streamlining the calibration process and improving efficiency.•The importance of compliance with international standards to ensure the quality and reliability of measurements in radiation protection was examined.•Description of X-ray and Gamma-ray irradiators designed for the calibration of radiation protection measuring instruments within the Secondary Dosimetry Calibration Laboratory (SSDL) which is a member of the WHO/IAEA network within the National Center for Radiation Protection of Morocco•Graphical User Interface using python for the calibration of photon measurement instruments for radiation protection purposes was developped.

Shennong: A Python toolbox for audio speech features extraction

We introduce Shennong, a Python toolbox and command-line utility for audio speech features extraction. It implements a wide range of well-established state-of-the-art algorithms: spectro-temporal filters such as Mel-Frequency Cepstral Filterbank or Predictive Linear Filters, pre-trained neural networks, pitch estimators, speaker normalization methods, and post-processing algorithms. Shennong is an open source, reliable and extensible framework built on top of the popular Kaldi speech processing library. The Python implementation makes it easy to use by non-technical users and integrates with third-party speech modeling and machine learning tools from the Python ecosystem. This paper describes the Shennong software architecture, its core components, and implemented algorithms. Then, three applications illustrate its use. We first present a benchmark of speech features extraction algorithms available in Shennong on a phone discrimination task. We then analyze the performances of a speaker normalization model as a function of the speech duration used for training. We finally compare pitch estimation algorithms on speech under various noise conditions.

kakapo: easy extraction and annotation of genes from raw RNA-seq reads

kakapo (kākāpō) is a Python-based pipeline that allows users to extract and assemble one or more specified genes or gene families. It flexibly uses original RNA-seq read or GenBank SRA accession inputs without performing global assembly of entire transcriptomes or metatranscriptomes. The pipeline identifies open reading frames in the assembled gene transcripts and annotates them. It optionally filters raw reads for ribosomal, plastid, and mitochondrial reads, or reads belonging to non-target organisms (e.g., viral, bacterial, human). kakapo can be employed for targeted assembly, to extract arbitrary loci, such as those commonly used for phylogenetic inference in systematics or candidate genes and gene families in phylogenomic and metagenomic studies. We provide example applications and discuss how its use can offset the declining value of GenBank's single-gene databases and help assemble datasets for a variety of phylogenetic analyses.

FLUTE: A Python GUI for interactive phasor analysis of FLIM data

Fluorescence lifetime imaging microscopy (FLIM) is a powerful technique used to probe the local environment of fluorophores. The fit-free phasor approach to FLIM data is increasingly being used due to its ease of interpretation. To date, no open-source graphical user interface (GUI) for phasor analysis of FLIM data is available in Python, thus limiting the widespread use of phasor analysis in biomedical research. Here, we present Fluorescence Lifetime Ultimate Explorer (FLUTE), a Python GUI that is designed to fill this gap. FLUTE simplifies and automates many aspects of the analysis of FLIM data acquired in the time domain, such as calibrating the FLIM data, performing interactive exploration of the phasor plot, displaying phasor plots and FLIM images with different lifetime contrasts simultaneously, and calculating the distance from known molecular species. After applying desired filters and thresholds, the final edited datasets can be exported for further user-specific analysis. FLUTE has been tested using several FLIM datasets including autofluorescence of zebrafish embryos and in vitro cells. In summary, our user-friendly GUI extends the advantages of phasor plotting by making the data visualization and analysis easy and interactive, allows for analysis of large FLIM datasets, and accelerates FLIM analysis for non-specialized labs.

A python system for regional landslide susceptibility assessment by integrating machine learning models and its application

Landslide susceptibility assessment is considered the first step in landslide risk assessment, but current studies mostly rely on GIS platforms or other software for data preprocessing. The modeling process is relatively complicated and multi-models cannot be integrated. With regard to this issue, this study develops a Python system for automatic assessment of regional landslide susceptibility. The Python system implements landslide susceptibility assessment through three modules: geographic data processing, machine learning modeling and result evaluation analysis. For geographic data processing, ten landslide influencing factors can be used to construct an evaluation factor dataset and reclassify the thematic maps based on the frequency ratio method. Four built-in machine learning models (logistic regression (LR), multi-layer perceptron (MLP), support vector machine (SVM) and extreme gradient boosting (XGBoost)) are integrated into the system to complete susceptibility modeling and calculation. Additionally, receiver operating characteristic (ROC) curves can be automatically generated to evaluate the accuracy. The system was then applied into Lantian County in Shaanxi Province as a demonstration example. The results show that the areas under the ROC curve (AUC) of the four models are 0.838 (LR)、0.882 (SVM)、0.809 (MLP) and 0.812 (XGBoost), respectively, indicating that the SVM model was the most suitable model for landslide susceptibility assessment in Lantian County in the Loess Plateau of China. The system has now been made open source on Github, which can effectively improve the efficiency of regional landslide susceptibility assessment, especially provide tools for data processing and modeling for non-professionals.

Using PharmaPy with Jupyter Notebook to teach digital design in pharmaceutical manufacturing

The use of digital tools in pharmaceutical manufacturing has gained traction over the past two decades. Whether supporting regulatory filings or attempting to modernize manufacturing processes to adopt new and quickly evolving Industry 4.0 standards, engineers entering the workforce must exhibit proficiency in modeling, simulation, optimization, data processing, and other digital analysis techniques. In this work, a course that addresses digital tools in pharmaceutical manufacturing for chemical engineers was adjusted to utilize a new tool, PharmaPy, instead of traditional chemical engineering simulation tools. Jupyter Notebook was utilized as an instructional and interactive environment to teach students to use PharmaPy, a new, open-source pharmaceutical manufacturing process simulator. Students were then surveyed to see if PharmaPy was able to meet the learning objectives of the course. During the semester, PharmaPy's model library was used to simulate both individual unit operations as well as multiunit pharmaceutical processes. Through the initial survey results, students indicated that: (i) through Jupyter Notebook, learning Python and PharmaPy was approachable from varied coding experience backgrounds and (ii) PharmaPy strengthened their understanding of pharmaceutical manufacturing through active pharmaceutical ingredient process design and development.

medigan: a Python library of pretrained generative models for medical image synthesis

Purpose

Deep learning has shown great promise as the backbone of clinical decision support systems. Synthetic data generated by generative models can enhance the performance and capabilities of data-hungry deep learning models. However, there is (1) limited availability of (synthetic) datasets and (2) generative models are complex to train, which hinders their adoption in research and clinical applications. To reduce this entry barrier, we explore generative model sharing to allow more researchers to access, generate, and benefit from synthetic data.

Approach

We propose medigan, a one-stop shop for pretrained generative models implemented as an open-source framework-agnostic Python library. After gathering end-user requirements, design decisions based on usability, technical feasibility, and scalability are formulated. Subsequently, we implement medigan based on modular components for generative model (i) execution, (ii) visualization, (iii) search & ranking, and (iv) contribution. We integrate pretrained models with applications across modalities such as mammography, endoscopy, x-ray, and MRI.

Results

The scalability and design of the library are demonstrated by its growing number of integrated and readily-usable pretrained generative models, which include 21 models utilizing nine different generative adversarial network architectures trained on 11 different datasets. We further analyze three medigan applications, which include (a) enabling community-wide sharing of restricted data, (b) investigating generative model evaluation metrics, and (c) improving clinical downstream tasks. In (b), we extract Fréchet inception distances (FID) demonstrating FID variability based on image normalization and radiology-specific feature extractors.

Conclusion

medigan allows researchers and developers to create, increase, and domain-adapt their training data in just a few lines of code. Capable of enriching and accelerating the development of clinical machine learning models, we show medigan's viability as platform for generative model sharing. Our multimodel synthetic data experiments uncover standards for assessing and reporting metrics, such as FID, in image synthesis studies.

Revolutionizing Diabetes Diagnosis: Machine Learning Techniques Unleashed

The intricate and multifaceted nature of diabetes disrupts the body's crucial glucose processing mechanism, which serves as a fundamental energy source for the cells. This research aims to predict the occurrence of diabetes in individuals by harnessing the power of machine learning algorithms, utilizing the PIMA diabetes dataset. The selected algorithms employed in this study encompass Decision Tree, K-Nearest Neighbor, Random Forest, Logistic Regression, and Support Vector Machine. To execute the experiments, two software tools, namely Waikato Environment for Knowledge Analysis (WEKA) version 3.8.1 and Python version 3.10, were utilized. To evaluate the performance of the algorithms, several metrics were employed, including true positive rate, false positive rate, precision, recall, F-measure, Matthew's correlation coefficient, receiver operating characteristic area, and precision-recall curves area. Furthermore, various errors such as Mean Absolute Error, Root Mean Squared Error, Relative Absolute Error, and Root Relative Squared Error were examined to assess the accuracy of the models. Upon conducting the experiments, it was observed that Logistic Regression outperformed the other techniques, exhibiting the highest precision of 81 percent using Python and 80.43 percent using WEKA. These findings shed light on the efficacy of machine learning in predicting diabetes and highlight the potential of Logistic Regression as a valuable tool in this domain.

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled60Co (Type M) and131I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method (expm). The differences generally observed between the implicit solvers andexpmare attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

Improving the usability of open health service delivery simulation models using Python and web apps

One aim of Open Science is to increase the accessibility of research. Within health services research that uses discrete-event simulation, Free and Open Source Software (FOSS), such as Python, offers a way for research teams to share their models with other researchers and NHS decision makers. Although the code for healthcare discrete-event simulation models can be shared alongside publications, it may require specialist skills to use and run. This is a disincentive to researchers adopting Free and Open Source Software and open science practices. Building on work from other health data science disciplines, we propose that web apps offer a user-friendly interface for healthcare models that increase the accessibility of research to the NHS, and researchers from other disciplines. We focus on models coded in Python deployed as streamlit web apps. To increase uptake of these methods, we provide an approach to structuring discrete-event simulation model code in Python so that models are web app ready. The method is general across discrete-event simulation Python packages, and we include code for both simpy and ciw implementations of a simple urgent care call centre model. We then provide a step-by-step tutorial for linking the model to a streamlit web app interface, to enable other health data science researchers to reproduce and implement our method.

A weather features dataset for prediction of short-term rainfall quantities in Uganda

Weather data is of great importance to the development of weather prediction models. However, the availability and quality of this data remains a significant challenge for most researchers around the world. In Uganda, obtaining observational weather data is very challenging due to the sparse distribution of weather stations and inconsistent data records. This has created critical gaps in data availability to run and develop efficient weather prediction models. To bridge this gap, we obtained country-specific time series hourly observational weather data. The data period is from 2020 to 2022 of 11 weather stations distributed in the four regions of Uganda. The data was accessed from the Ogimet data repository using the "climate" R-package. The automated procedures in the R-programming language environment allowed us to download user-defined data at a time resolution from an hourly to an annual basis. However, the raw data acquired cannot be used to learn rainfall patterns because it includes duplicates and non-uniform data. Therefore, this article presents a prepared and cleaned dataset that can be used for the prediction of short-term rainfall quantities in Uganda.