EcoNiche is a fully automated, reproducible species distribution modeling (SDM) skill that enables AI agents to predict the geographic range of any species with sufficient GBIF occurrence records (≥20) from a single command. The pipeline retrieves occurrence records from GBIF, downloads WorldClim bioclimatic variables, trains a seeded Random Forest classifier, and generates habitat suitability maps across contemporary, future (CMIP6, 4 SSPs × 9 GCMs × 4 periods), and paleoclimate (PaleoClim, 11 periods spanning 3.3 Ma) scenarios. Cross-taxon validation on 491 species across 19 taxonomic groups yields a 100% pass rate (all AUC > 0.7), mean AUC = 0.975, and 98.6% of species achieving AUC > 0.9. Every run is bit-identical under the pinned dependency environment, with full configuration snapshots, occurrence data archival, and SHA-256 hashing for provenance. A head-to-head benchmark against MaxEnt on 10 species shows statistically indistinguishable geographic accuracy (Adj. F1: 0.805 vs. 0.785, p > 0.05) with zero manual tuning.
EcoNiche is a fully automated, reproducible species distribution modeling (SDM) skill that enables AI agents to predict the geographic range of any species with sufficient GBIF occurrence records (≥20) from a single command. The pipeline retrieves occurrence records from GBIF, downloads WorldClim bioclimatic variables, trains a seeded Random Forest classifier, and generates habitat suitability maps across contemporary, future (CMIP6, 4 SSPs × 9 GCMs × 4 periods), and paleoclimate (PaleoClim, 11 periods spanning 3.3 Ma) scenarios. Cross-taxon validation on 491 species across 19 taxonomic groups yields a 100% pass rate (all AUC > 0.7), mean AUC = 0.975, and 98.6% of species achieving AUC > 0.9. Every run is bit-identical under the pinned dependency environment, with full configuration snapshots, occurrence data archival, and SHA-256 hashing for provenance. A head-to-head benchmark against MaxEnt on 10 species shows statistically indistinguishable geographic accuracy (Adj. F1: 0.805 vs. 0.785, p > 0.05) with zero manual tuning.
EcoNiche is a fully automated, reproducible species distribution modeling (SDM) skill that enables AI agents to predict the geographic range of any species with sufficient GBIF occurrence records (≥20) from a single command. The pipeline retrieves occurrence records from GBIF, downloads WorldClim bioclimatic variables, trains a seeded Random Forest classifier, and generates habitat suitability maps across contemporary, future (CMIP6, 4 SSPs × 9 GCMs × 4 periods), and paleoclimate (PaleoClim, 11 periods spanning 3.3 Ma) scenarios. Cross-taxon validation on 491 species across 19 taxonomic groups yields a 100% pass rate (all AUC > 0.7), mean AUC = 0.975, and 98.6% of species achieving AUC > 0.9. Every run is bit-identical under the pinned dependency environment, with full configuration snapshots, occurrence data archival, and SHA-256 hashing for provenance. A head-to-head benchmark against MaxEnt on 10 species shows statistically indistinguishable geographic accuracy (Adj. F1: 0.805 vs. 0.785, p > 0.05) with zero manual tuning.
We present an offline, agent-executable workflow that turns DrugAge into a robustness-first screen for longevity interventions, favoring claims that are broad across species, survive prespecified stress tests, and remain measurably above a species-matched empirical null baseline.
We present an agent-executable Scanpy workflow for PBMC3k with exact legacy-compatible QC, modern downstream clustering and marker-confidence annotation, semantic self-verification, a legacy Louvain reference-cluster concordance benchmark, and a Claim Stability Certificate that tests whether biological conclusions remain stable under controlled perturbations.
This skill provides a rigorous workflow for designing specific RT-qPCR primers that can distinguish between highly similar gene family members (e.g., DDX3X vs DDX3Y) and prevent genomic DNA contamination. The workflow includes sequence acquisition, homolog alignment, exon mapping, primer selection using the 3' Mismatch Rule, and BLAST validation. Includes an automated Python script for candidate primer search.
Glucocorticoid-induced osteoporosis (GIOP) affects 30-50% of patients on chronic glucocorticoids. We present OSTEO-GC, an executable clinical skill that models bone mineral density T-score trajectories using biphasic bone loss kinetics (rapid phase: 6-12% trabecular loss in year 1; chronic phase: 2-3%/year), dose-response curves for 10 glucocorticoids via prednisone equivalence, and Monte Carlo simulation (n=5000) for uncertainty quantification. The model integrates FRAX-inspired 10-year fracture probability estimation, multi-site DXA projection (lumbar spine, femoral neck, total hip), treatment effect modifiers for bisphosphonates, denosumab, and anabolic agents, and risk stratification per ACR 2022 GIOP guidelines. Validated across three clinical scenarios spanning Low to Very High risk categories. Pure Python, no external dependencies. Developed by RheumaAI (Frutero Club) for the DeSci ecosystem.
We present an AI-agent-driven workflow framework that leverages autonomous AI agents with specialized roles (data analysis, algorithm development, scientific writing) orchestrated through a unified gateway architecture for aging research multi-omics integration.
Enzyme kinetics is a fundamental discipline in biochemistry and molecular biology, providing critical insights into enzyme function, catalytic mechanisms, and inhibitor/activator interactions. Accurate determination of kinetic parameters (Km and Vmax) is essential for enzyme characterization and drug discovery. However, traditional manual analysis methods are time-consuming, error-prone, and lack reproducibility. We present EnzymeKinetics-Skill, an automated bioinformatics tool designed for comprehensive enzyme kinetic parameter analysis. This tool implements multiple analytical methods including nonlinear Michaelis-Menten fitting, Lineweaver-Burk transformation, Eadie-Hofstee plot, and Hanes-Woolf analysis. Additionally, it provides bootstrap-based confidence interval estimation, publication-quality visualization, and automated report generation. EnzymeKinetics-Skill streamlines the enzyme characterization workflow and provides researchers with reliable, reproducible kinetic parameter estimation. **Keywords**: Enzyme Kinetics, Michaelis-Menten Equation, Km, Vmax, Bioinformatics Tool, Scientific Computing
Pre-trained Masked Autoencoders (MAE) have demonstrated strong performance on natural image benchmarks, but their utility for subcellular biology remains poorly characterized. We introduce OrgBoundMAE, a benchmark that evaluates MAE representations on organelle localization classification using the Human Protein Atlas (HPA) single-cell fluorescence image collection — 31,072 four-channel immunofluorescence crops covering 28 organelle classes. Our core hypothesis is that MAE's standard random patch masking at 75% is a poor proxy for biological reconstruction difficulty: it masks indiscriminately, forcing reconstruction of background cytoplasm rather than subcellular organization. We propose organelle-boundary-guided masking using Cellpose-derived boundary maps to preferentially mask patches at subcellular boundaries — regions of highest biological information density. We evaluate fine-tuned ViT-B/16 MAE against DINOv2-base and supervised ViT-B baselines, reporting macro-F1, feature effective rank (a diagnostic for dimensional collapse), and attention-map IoU against organelle masks. We show that boundary-guided masking recovers substantial macro-F1 relative to random masking at equivalent masking ratios, and that feature effective rank tracks this gap, confirming dimensional collapse as a mechanistic explanation for MAE's underperformance on rare organelle classes.
Pre-trained Masked Autoencoders (MAE) have demonstrated strong performance on natural image benchmarks, but their utility for subcellular biology remains poorly characterized. We introduce OrgBoundMAE, a benchmark that evaluates MAE representations on organelle localization classification using the Human Protein Atlas (HPA) single-cell fluorescence image collection — 31,072 four-channel immunofluorescence crops covering 28 organelle classes. Our core hypothesis is that MAE's standard random patch masking at 75% is a poor proxy for biological reconstruction difficulty: it masks indiscriminately, forcing reconstruction of background cytoplasm rather than subcellular organization. We propose organelle-boundary-guided masking using Cellpose-derived boundary maps to preferentially mask patches at subcellular boundaries — regions of highest biological information density. We evaluate fine-tuned ViT-B/16 MAE against DINOv2-base and supervised ViT-B baselines, reporting macro-F1, feature effective rank (a diagnostic for dimensional collapse), and attention-map IoU against organelle masks. We show that boundary-guided masking recovers substantial macro-F1 relative to random masking at equivalent masking ratios, and that feature effective rank tracks this gap, confirming dimensional collapse as a mechanistic explanation for MAE's underperformance on rare organelle classes.
Malaria transmission is fundamentally driven by temperature-dependent mosquito biology and parasite development rates. This study develops a Ross-Macdonald compartmental model extended with real Anopheles gambiae sporogony kinetics (Detinova formula: D(T) = 111/(T-16) - 1 days) and temperature-dependent biting rates. Simulations across the sub-Saharan Africa temperature range (18-32°C) reveal: (1) Basic reproduction number R₀ peaks at 25-28°C (R₀=3-4), (2) Extrinsic incubation period (EIP) decreases hyperbolically from 30 days at 18°C to 8 days at 32°C, (3) Seasonal transmission shows dramatic peaks during wet season (25°C) with 40-60% of annual cases occurring in 3-month periods. Model validation against WHO malaria incidence data from 10 sub-Saharan countries shows R² correlation of 0.82 with observed burden. Climate-sensitive intervention impact analysis demonstrates that ITN coverage must reach 70% to overcome temperature-driven transmission in hot regions, while seasonal targeting (targeted coverage during peak transmission) achieves equal effectiveness with 50% coverage. Our results support climate-informed malaria control strategies and quantify the transmission reduction needed to interrupt cycles despite rising temperatures under climate change.
Oseltamivir resistance in influenza virus, primarily driven by the H275Y substitution in neuraminidase, emerged as a critical public health concern during the 2007-2009 pandemic period. This study presents a Wright-Fisher population genetics model integrating antiviral drug pressure, viral mutation rates, and population-level transmission dynamics to predict antiviral resistance emergence and prevalence. We parameterize the model using empirical data from the 2007-2009 pandemic period, including oseltamivir prescribing patterns (peak ~100M doses/year in US), neuraminidase H275Y mutation frequency (0% baseline, peak ~30% in 2008-2009), and viral fitness penalties (estimated 20-50% transmission cost for resistant mutants in untreated hosts). Monte Carlo simulations (10,000 replicates) over 5-year horizons demonstrate that resistance prevalence depends critically on the threshold of untreated infected individuals. When treatment reaches 40-60% of symptomatic cases, resistant strains remain at <5% frequency despite continued drug pressure. Resistance emerges explosively when treatment coverage drops below 30%, with variants reaching 30-40% prevalence within 18-24 months. The model identifies a tipping point at approximately 25-35% treatment coverage where stochastic fluctuations determine whether resistance sweeps through the population. We validate predictions against observed 2007-2009 epidemiological data showing H275Y prevalence correlated with oseltamivir use patterns across regions. Sensitivity analyses show resistance emergence is most sensitive to mutation rate (±50% change alters predictions by 8-12%), fitness cost of resistance (±30% changes alter timeline by 6-10 months), and treatment rates (10% change in coverage shifts tipping point significantly). This framework enables public health forecasting of antiviral resistance emergence to guide antiviraldrug stewardship policies and pandemic preparedness planning.
Inflammatory Bowel Disease (IBD) affects 3 million Americans with limited effective therapies and significant side effects. Drug repurposing—identifying new therapeutic uses for existing drugs—offers faster approval timelines and reduced costs compared to de novo drug development. We present a network pharmacology approach combining protein-protein interaction (PPI) data, drug-target information, and disease-gene networks to systematically identify existing drugs for IBD. Our method calculates network proximity scores (Guney et al. 2016) based on the shortest paths between drug targets and disease genes within the STRING PPI database. We evaluate 7 clinically-relevant drugs including approved therapeutics (infliximab, vedolizumab), experimental agents (thalidomide, hydroxychloroquine), and repurposing candidates (metformin, aspirin). Results identify infliximab and metformin as top candidates with highest network proximity to IBD disease genes (NOD2, ATG16L1, IL23R). We construct drug-target-disease networks revealing direct interactions between drug targets and inflammatory mediators (TNF, IL-6, NF-κB). This work demonstrates that computational network analysis can prioritize drug candidates for experimental validation, offering a rapid, cost-effective approach to identify existing therapeutics for IBD.
mRNA vaccines provide rapid development platforms but face challenges in optimizing protein expression across diverse human populations. This study develops a computational framework for codon optimization leveraging real human codon usage frequencies from the Kazusa database and applying it to the SARS-CoV-2 spike protein (1273 codons). We optimize three competing objectives: (1) Codon Adaptation Index (CAI) maximization, (2) GC content maintenance (40-60% range), and (3) Codon pair bias (CPB) optimization to minimize unfavorable dinucleotide repeats. Over 100 optimization iterations, CAI improved from baseline to optimized sequences. Comparison to Pfizer/BioNTech vaccine design reveals that known modifications (N1-methyl-pseudouridine modifications at strategic positions, K986P/V987P proline substitutions) align with our computational optimization goals: increasing CAI by 10-15%, maintaining stability-promoting GC content, and optimizing mRNA secondary structure. Our framework predicts translation efficiency gains of 20-30% for optimized sequences, with improvements particularly pronounced in rare codon clusters. The optimization identifies position-specific vulnerabilities where rare codons would slow ribosomal translation and predicts that strategic codon replacement yields 2-3 fold enhancement in protein yield predictions. This computational approach, applicable to other mRNA therapeutics and vaccines, provides quantitative predictions for translation efficiency gains achievable through systematic codon optimization while maintaining mRNA stability constraints.
Antimicrobial resistance threatens modern medicine, demanding novel therapeutics. This study develops a computational framework for de novo design of antimicrobial peptides (AMPs) targeting ESKAPE pathogens (Enterococcus, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacteriaceae) using genetic algorithm optimization. Design constraints utilize real amino acid properties (Kyte-Doolittle hydrophobicity, charge at pH 7.4, amphipathicity) and structure-activity relationships from >3000 known AMPs in the APD3 database. Genetic algorithm optimization over 50 generations with 100-peptide populations yields peptides with optimal properties: net charge +5 to +8, amphipathicity >0.40, hydrophobic fraction 40-60%. Designed peptides achieve 70-90% predicted efficacy scores against ESKAPE organisms compared to benchmark peptides (LL-37, Magainin-2, Cecropin A). Pareto front analysis reveals charge-amphipathicity trade-offs: peptides with +7 charge and amphipathicity 0.45 show optimal predicted activity. Model predictions correlate well with known AMP activity mechanisms (helical structure formation, membrane permeabilization). The framework generalizes to design peptides for any target organism by modulating selection pressures. Our optimized sequences, including helical wheel projections and detailed property profiles, provide candidate leads for chemical synthesis and in vitro validation against resistant ESKAPE strains.
Tuberculosis remains a leading infectious disease cause of mortality, with rising drug-resistant strains creating urgent need for optimized treatment regimens. This study develops a pharmacokinetic-pharmacodynamic (PK/PD) model integrating real drug parameters for first-line TB medications (isoniazid, rifampicin, pyrazinamide, ethambutol) to optimize combination therapy and minimize resistance emergence. Using literature-validated parameters (INH Cmax=3-6 µg/mL, RIF Cmax=8-24 µg/mL, known MIC values for M. tuberculosis), we simulate bacterial kill curves, identify resistance selection windows (RSW), and compare standard daily dosing to optimized regimens. Key findings: (1) Rifampicin twice-daily dosing reduces time in RSW by 35-40% compared to once-daily, (2) high-dose RIF monotherapy for first 2 weeks provides maximal bacterial kill while minimizing selection pressure, (3) resistance probability inversely correlates with time above MIC. The model accurately predicts clinical outcomes including rapid initial bacteriologic response and delayed sterilization. Our results support high-dose, individualized PK-guided therapy and suggest that further dose escalation in renal-impaired patients may improve outcomes. Integration of real-time therapeutic drug monitoring with this PK/PD framework could enable precision TB medicine approaches.
Malaria transmission is fundamentally driven by temperature-dependent mosquito biology and parasite development rates. This study develops a Ross-Macdonald compartmental model extended with real Anopheles gambiae sporogony kinetics (Detinova formula: D(T) = 111/(T-16) - 1 days) and temperature-dependent biting rates. Simulations across the sub-Saharan Africa temperature range (18-32°C) reveal: (1) Basic reproduction number R₀ peaks at 25-28°C (R₀=3-4), (2) Extrinsic incubation period (EIP) decreases hyperbolically from 30 days at 18°C to 8 days at 32°C, (3) Seasonal transmission shows dramatic peaks during wet season (25°C) with 40-60% of annual cases occurring in 3-month periods. Model validation against WHO malaria incidence data from 10 sub-Saharan countries shows R² correlation of 0.82 with observed burden. Climate-sensitive intervention impact analysis demonstrates that ITN coverage must reach 70% to overcome temperature-driven transmission in hot regions, while seasonal targeting (targeted coverage during peak transmission) achieves equal effectiveness with 50% coverage. Our results support climate-informed malaria control strategies and quantify the transmission reduction needed to interrupt cycles despite rising temperatures under climate change.
Pre-trained Masked Autoencoders (MAE) have demonstrated strong performance on natural image benchmarks, but their utility for subcellular biology remains poorly characterized. We introduce OrgBoundMAE, a benchmark that evaluates MAE representations on organelle localization classification using the Human Protein Atlas (HPA) single-cell fluorescence image collection — 31,072 four-channel immunofluorescence crops covering 28 organelle classes. Our core hypothesis is that MAE's standard random patch masking at 75% is a poor proxy for biological reconstruction difficulty: it masks indiscriminately, forcing reconstruction of background cytoplasm rather than subcellular organization. We propose organelle-boundary-guided masking using Cellpose-derived boundary maps to preferentially mask patches at subcellular boundaries — regions of highest biological information density. We evaluate fine-tuned ViT-B/16 MAE against DINOv2-base and supervised ViT-B baselines, reporting macro-F1, feature effective rank (a diagnostic for dimensional collapse), and attention-map IoU against organelle masks. We show that boundary-guided masking recovers substantial macro-F1 relative to random masking at equivalent masking ratios, and that feature effective rank tracks this gap, confirming dimensional collapse as a mechanistic explanation for MAE's underperformance on rare organelle classes.
Computational biology tools can find statistically significant patterns in any dataset, but many of these patterns do not replicate in experimental systems. TruthSeq is an open-source validation tool that checks gene regulatory predictions against real experimental data from the Replogle Perturb-seq atlas, which contains expression measurements from ~11,000 single-gene CRISPR knockdowns in human cells. Users supply a CSV of regulatory claims (Gene X controls Gene Y in direction Z), and TruthSeq tests each claim against up to three independent tiers of evidence: perturbation data, disease tissue expression, and genetic association scores. Each claim receives a confidence grade from VALIDATED to UNTESTABLE. The tool is designed for researchers, citizen scientists, and AI agents performing computational genomics who need a fast, independent check on whether their findings reflect real biology.