Structural basis of membrane skeleton organization in red blood cells.

The spectrin-based membrane skeleton is a ubiquitous membrane-associated two-dimensional cytoskeleton underneath the lipid membrane of metazoan cells. Mutations of skeleton proteins impair the mechanical strength and functions of the membrane, leading to several different types of human diseases. Here, we report the cryo-EM structures of the native spectrin-actin junctional complex (from porcine erythrocytes), which is a specialized short F-actin acting as the central organizational unit of the membrane skeleton. While an α-/ß-adducin hetero-tetramer binds to the barbed end of F-actin as a flexible cap, tropomodulin and SH3BGRL2 together create an absolute cap at the pointed end. The junctional complex is strengthened by ring-like structures of dematin in the middle actin layers and by patterned periodic interactions with tropomyosin over its entire length. This work serves as a structural framework for understanding the assembly and dynamics of membrane skeleton and offers insights into mechanisms of various ubiquitous F-actin-binding factors in other F-actin systems.

Structurally Resolved SARS-CoV-2 Antibody Shows High Efficacy in Severely Infected Hamsters and Provides a Potent Cocktail Pairing Strategy.

Understanding how potent neutralizing antibodies (NAbs) inhibit SARS-CoV-2 is critical for effective therapeutic development. We previously described BD-368-2, a SARS-CoV-2 NAb with high potency; however, its neutralization mechanism is largely unknown. Here, we report the 3.5-Å cryo-EM structure of BD-368-2/trimeric-spike complex, revealing that BD-368-2 fully blocks ACE2 recognition by occupying all three receptor-binding domains (RBDs) simultaneously, regardless of their "up" or "down" conformations. Also, BD-368-2 treats infected adult hamsters at low dosages and at various administering windows, in contrast to placebo hamsters that manifested severe interstitial pneumonia. Moreover, BD-368-2's epitope completely avoids the common binding site of VH3-53/VH3-66 recurrent NAbs, evidenced by tripartite co-crystal structures with RBDs. Pairing BD-368-2 with a potent recurrent NAb neutralizes SARS-CoV-2 pseudovirus at pM level and rescues mutation-induced neutralization escapes. Together, our results rationalized a new RBD epitope that leads to high neutralization potency and demonstrated BD-368-2's therapeutic potential in treating COVID-19.

Potent Neutralizing Antibodies against SARS-CoV-2 Identified by High-Throughput Single-Cell Sequencing of Convalescent Patients' B Cells.

The COVID-19 pandemic urgently needs therapeutic and prophylactic interventions. Here, we report the rapid identification of SARS-CoV-2-neutralizing antibodies by high-throughput single-cell RNA and VDJ sequencing of antigen-enriched B cells from 60 convalescent patients. From 8,558 antigen-binding IgG1+ clonotypes, 14 potent neutralizing antibodies were identified, with the most potent one, BD-368-2, exhibiting an IC50 of 1.2 and 15 ng/mL against pseudotyped and authentic SARS-CoV-2, respectively. BD-368-2 also displayed strong therapeutic and prophylactic efficacy in SARS-CoV-2-infected hACE2-transgenic mice. Additionally, the 3.8 Å cryo-EM structure of a neutralizing antibody in complex with the spike-ectodomain trimer revealed the antibody's epitope overlaps with the ACE2 binding site. Moreover, we demonstrated that SARS-CoV-2-neutralizing antibodies could be directly selected based on similarities of their predicted CDR3H structures to those of SARS-CoV-neutralizing antibodies. Altogether, we showed that human neutralizing antibodies could be efficiently discovered by high-throughput single B cell sequencing in response to pandemic infectious diseases.

Cryo-EM Structure and Assembly of an Extracellular Contractile Injection System.

Contractile injection systems (CISs) are cell-puncturing nanodevices that share ancestry with contractile tail bacteriophages. Photorhabdus virulence cassette (PVC) represents one group of extracellular CISs that are present in both bacteria and archaea. Here, we report the cryo-EM structure of an intact PVC from P. asymbiotica. This over 10-MDa device resembles a simplified T4 phage tail, containing a hexagonal baseplate complex with six fibers and a capped 117-nanometer sheath-tube trunk. One distinct feature of the PVC is the presence of three variants for both tube and sheath proteins, indicating a functional specialization of them during evolution. The terminal hexameric cap docks onto the topmost layer of the inner tube and locks the outer sheath in pre-contraction state with six stretching arms. Our results on the PVC provide a framework for understanding the general mechanism of widespread CISs and pave the way for using them as delivery tools in biological or therapeutic applications.

Structure of a Pancreatic ATP-Sensitive Potassium Channel.

ATP-sensitive potassium channels (KATP) couple intracellular ATP levels with membrane excitability. These channels play crucial roles in many essential physiological processes and have been implicated extensively in a spectrum of metabolic diseases and disorders. To gain insight into the mechanism of KATP, we elucidated the structure of a hetero-octameric pancreatic KATP channel in complex with a non-competitive inhibitor glibenclamide by single-particle cryoelectron microscopy to 5.6-Å resolution. The structure shows that four SUR1 regulatory subunits locate peripherally and dock onto the central Kir6.2 channel tetramer through the SUR1 TMD0-L0 fragment. Glibenclamide-bound SUR1 uses TMD0-L0 fragment to stabilize Kir6.2 channel in a closed conformation. In another structural population, a putative co-purified phosphatidylinositol 4,5-bisphosphate (PIP2) molecule uncouples Kir6.2 from glibenclamide-bound SUR1. These structural observations suggest a molecular mechanism for KATP regulation by anti-diabetic sulfonylurea drugs, intracellular adenosine nucleotide concentrations, and PIP2 lipid.

Structural insights into the assembly and mechanism of mpox virus DNA polymerase complex F8-A22-E4-H5.

The DNA replication of mpox virus is performed by the viral polymerase F8 and also requires other viral factors, including processivity factor A22, uracil DNA glycosylase E4, and phosphoprotein H5. However, the molecular roles of these viral factors remain unclear. Here, we characterize the structures of F8-A22-E4 and F8-A22-E4-H5 complexes in the presence of different primer-template DNA substrates. E4 is located upstream of F8 on the template single-stranded DNA (ssDNA) and is catalytically active, highlighting a functional coupling between DNA base-excision repair and DNA synthesis. Moreover, H5, in the form of tetramer, binds to the double-stranded DNA (dsDNA) region downstream of F8 in a similar position as PCNA (proliferating cell nuclear antigen) does in eukaryotic polymerase complexes. Omission of H5 or disruption of its DNA interaction showed a reduced synthesis of full-length DNA products. These structures provide snapshots for the working cycle of the polymerase and generate insights into the mechanisms of these essential factors in viral DNA replication.

Parental histone transfer caught at the replication fork.

In eukaryotes, DNA compacts into chromatin through nucleosomes1,2. Replication of the eukaryotic genome must be coupled to the transmission of the epigenome encoded in the chromatin3,4. Here we report cryo-electron microscopy structures of yeast (Saccharomyces cerevisiae) replisomes associated with the FACT (facilitates chromatin transactions) complex (comprising Spt16 and Pob3) and an evicted histone hexamer. In these structures, FACT is positioned at the front end of the replisome by engaging with the parental DNA duplex to capture the histones through the middle domain and the acidic carboxyl-terminal domain of Spt16. The H2A-H2B dimer chaperoned by the carboxyl-terminal domain of Spt16 is stably tethered to the H3-H4 tetramer, while the vacant H2A-H2B site is occupied by the histone-binding domain of Mcm2. The Mcm2 histone-binding domain wraps around the DNA-binding surface of one H3-H4 dimer and extends across the tetramerization interface of the H3-H4 tetramer to the binding site of Spt16 middle domain before becoming disordered. This arrangement leaves the remaining DNA-binding surface of the other H3-H4 dimer exposed to additional interactions for further processing. The Mcm2 histone-binding domain and its downstream linker region are nested on top of Tof1, relocating the parental histones to the replisome front for transfer to the newly synthesized lagging-strand DNA. Our findings offer crucial structural insights into the mechanism of replication-coupled histone recycling for maintaining epigenetic inheritance.

Nuclear export of pre-60S particles through the nuclear pore complex.

The nuclear pore complex (NPC) is the bidirectional gate that mediates the exchange of macromolecules or their assemblies between nucleus and cytoplasm1-3. The assembly intermediates of the ribosomal subunits, pre-60S and pre-40S particles, are among the largest cargoes of the NPC and the export of these gigantic ribonucleoproteins requires numerous export factors4,5. Here we report the cryo-electron microscopy structure of native pre-60S particles trapped in the channel of yeast NPCs. In addition to known assembly factors, multiple factors with export functions are also included in the structure. These factors in general bind to either the flexible regions or subunit interface of the pre-60S particle, and virtually form many anchor sites for NPC binding. Through interactions with phenylalanine-glycine (FG) repeats from various nucleoporins of NPC, these factors collectively facilitate the passage of the pre-60S particle through the central FG repeat network of the NPC. Moreover, in silico analysis of the axial and radial distribution of pre-60S particles within the NPC shows that a single NPC can take up to four pre-60S particles simultaneously, and pre-60S particles are enriched in the inner ring regions close to the wall of the NPC with the solvent-exposed surface facing the centre of the nuclear pore. Our data suggest a translocation model for the export of pre-60S particles through the NPC.

Structure, function and pharmacology of human itch receptor complexes.

In the clades of animals that diverged from the bony fish, a group of Mas-related G-protein-coupled receptors (MRGPRs) evolved that have an active role in itch and allergic signals1,2. As an MRGPR, MRGPRX2 is known to sense basic secretagogues (agents that promote secretion) and is involved in itch signals and eliciting pseudoallergic reactions3-6. MRGPRX2 has been targeted by drug development efforts to prevent the side effects induced by certain drugs or to treat allergic diseases. Here we report a set of cryo-electron microscopy structures of the MRGPRX2-Gi1 trimer in complex with polycationic compound 48/80 or with inflammatory peptides. The structures of the MRGPRX2-Gi1 complex exhibited shallow, solvent-exposed ligand-binding pockets. We identified key common structural features of MRGPRX2 and describe a consensus motif for peptidic allergens. Beneath the ligand-binding pocket, the unusual kink formation at transmembrane domain 6 (TM6) and the replacement of the general toggle switch from Trp6.48 to Gly6.48 (superscript annotations as per Ballesteros-Weinstein nomenclature) suggest a distinct activation process. We characterized the interfaces of MRGPRX2 and the Gi trimer, and mapped the residues associated with key single-nucleotide polymorphisms on both the ligand and G-protein interfaces of MRGPRX2. Collectively, our results provide a structural basis for the sensing of cationic allergens by MRGPRX2, potentially facilitating the rational design of therapies to prevent unwanted pseudoallergic reactions.

Structure of GPR101-Gs enables identification of ligands with rejuvenating potential.

GPR101 is an orphan G protein-coupled receptor actively participating in energy homeostasis. Here we report the cryo-electron microscopy structure of GPR101 constitutively coupled to Gs heterotrimer, which reveals unique features of GPR101, including the interaction of extracellular loop 2 within the 7TM bundle, a hydrophobic chain packing-mediated activation mechanism and the structural basis of disease-related mutants. Importantly, a side pocket is identified in GPR101 that facilitates in silico screening to identify four small-molecule agonists, including AA-14. The structure of AA-14-GPR101-Gs provides direct evidence of the AA-14 binding at the side pocket. Functionally, AA-14 partially restores the functions of GH/IGF-1 axis and exhibits several rejuvenating effects in wild-type mice, which are abrogated in Gpr101-deficient mice. In summary, we provide a structural basis for the constitutive activity of GPR101. The structure-facilitated identification of GPR101 agonists and functional analysis suggest that targeting this orphan receptor has rejuvenating potential.

Comprehensive structural characterization of the human AAA+ disaggregase CLPB in the apo- and substrate-bound states reveals a unique mode of action driven by oligomerization.

The human AAA+ ATPase CLPB (SKD3) is a protein disaggregase in the mitochondrial intermembrane space (IMS) and functions to promote the solubilization of various mitochondrial proteins. Loss-of-function CLPB mutations are associated with a few human diseases with neutropenia and neurological disorders. Unlike canonical AAA+ proteins, CLPB contains a unique ankyrin repeat domain (ANK) at its N-terminus. How CLPB functions as a disaggregase and the role of its ANK domain are currently unclear. Herein, we report a comprehensive structural characterization of human CLPB in both the apo- and substrate-bound states. CLPB assembles into homo-tetradecamers in apo-state and is remodeled into homo-dodecamers upon substrate binding. Conserved pore-loops (PLs) on the ATPase domains form a spiral staircase to grip and translocate the substrate in a step-size of 2 amino acid residues. The ANK domain is not only responsible for maintaining the higher-order assembly but also essential for the disaggregase activity. Interactome analysis suggests that the ANK domain may directly interact with a variety of mitochondrial substrates. These results reveal unique properties of CLPB as a general disaggregase in mitochondria and highlight its potential as a target for the treatment of various mitochondria-related diseases.

Structural basis of assembly of the human T cell receptor-CD3 complex.

The αß T cell receptor (TCR), in association with the CD3γε-CD3δε-CD3ζζ signalling hexamer, is the primary determinant of T cell development and activation, and of immune responses to foreign antigens. The mechanism of assembly of the TCR-CD3 complex remains unknown. Here we report a cryo-electron microscopy structure of human TCRαß in complex with the CD3 hexamer at 3.7 Å resolution. The structure contains the complete extracellular domains and all the transmembrane helices of TCR-CD3. The octameric TCR-CD3 complex is assembled with 1:1:1:1 stoichiometry of TCRαß:CD3γε:CD3δε:CD3ζζ. Assembly of the extracellular domains of TCR-CD3 is mediated by the constant domains and connecting peptides of TCRαß that pack against CD3γε-CD3δε, forming a trimer-like structure proximal to the plasma membrane. The transmembrane segment of the CD3 complex adopts a barrel-like structure formed by interaction of the two transmembrane helices of CD3ζζ with those of CD3γε and CD3δε. Insertion of the transmembrane helices of TCRαß into the barrel-like structure via both hydrophobic and ionic interactions results in transmembrane assembly of the TCR-CD3 complex. Together, our data reveal the structural basis for TCR-CD3 complex assembly, providing clues to TCR triggering and a foundation for rational design of immunotherapies that target the complex.

Unique Roles of the Non-identical MCM Subunits in DNA Replication Licensing.

A family of six homologous subunits, Mcm2, -3, -4, -5, -6, and -7, each with its own unique features, forms the catalytic core of the eukaryotic replicative helicase. The necessity of six similar but non-identical subunits has been a mystery since its initial discovery. Recent cryo-EM structures of the Mcm2-7 (MCM) double hexamer, its precursors, and the origin recognition complex (ORC)-Cdc6-Cdt1-Mcm2-7 (OCCM) intermediate showed that each of these subunits plays a distinct role in orchestrating the assembly of the pre-replication complex (pre-RC) by ORC-Cdc6 and Cdt1.

Engineering of the DNA replication and repair machinery to develop binary mutators for rapid genome evolution of Corynebacterium glutamicum.

Corynebacterium glutamicum is an important industrial workhorse for production of amino acids and chemicals. Although recently developed genome editing technologies have advanced the rational genetic engineering of C. glutamicum, continuous genome evolution based on genetic mutators is still unavailable. To address this issue, the DNA replication and repair machinery of C. glutamicum was targeted in this study. DnaQ, the homolog of ϵ subunit of DNA polymerase III responsible for proofreading in Escherichia coli, was proven irrelevant to DNA replication fidelity in C. glutamicum. However, the histidinol phosphatase (PHP) domain of DnaE1, the α subunit of DNA polymerase III, was characterized as the key proofreading element and certain variants with PHP mutations allowed elevated spontaneous mutagenesis. Repression of the NucS-mediated post-replicative mismatch repair pathway or overexpression of newly screened NucS variants also impaired the DNA replication fidelity. Simultaneous interference with the DNA replication and repair machinery generated a binary genetic mutator capable of increasing the mutation rate by up to 2352-fold. The mutators facilitated rapid evolutionary engineering of C. glutamicum to acquire stress tolerance and protein overproduction phenotypes. This study provides efficient tools for evolutionary engineering of C. glutamicum and could inspire the development of mutagenesis strategy for other microbial hosts.

Engineering and functional analysis of yeast with a monotypic 40S ribosome subunit.

Emerging evidence reveals that ribosomes are not monolithic but dynamic machines with heterogeneous protein compositions that can reshape ribosomal translational abilities and cellular adaptation to environmental changes. Duplications of ribosomal protein (RP) genes are ubiquitous among organisms and are believed to affect cell function through paralog-specific regulation (e.g., by generating heterogeneous ribosomes) and/or gene dose amplification. However, direct evaluations of their impacts on cell function remain elusive due to the highly heterogeneous cellular RP pool. Here, we engineered a yeast with homogeneous 40S RP paralog compositions, designated homo-40S, by deleting the entire set of alternative duplicated genes encoding yeast 40S RP paralogs. Homo-40S displayed mild growth defects along with high sensitivity to the translation inhibitor paromomycin and a significantly increased stop codon readthrough. Moreover, doubling of the remaining RP paralogous genes in homo-40S rescued these phenotypes markedly, although not fully, compared to the wild-type phenotype, indicating that the dose of 40S RP genes together with the heterogeneity of the contents was vital for maintaining normal translational functionalities and growth robustness. Additional experiments revealed that homo-40S improved paromomycin tolerance via acquisition of bypass mutations or evolved to be diploid to generate fast-growing derivatives, highlighting the mutational robustness of engineered yeast to accommodate environmental and genetic changes. In summary, our work demonstrated that duplicated RP paralogs impart robustness and phenotypic plasticity through both gene dose amplification and paralog-specific regulation, paving the way for the direct study of ribosome biology through monotypic ribosomes with a homogeneous composition of specific RP paralogs.

Machine learning versus multivariate logistic regression for predicting severe COVID-19 in hospitalized children with Omicron variant infection.

With the emergence of the Omicron variant, the number of pediatric Coronavirus Disease 2019 (COVID-19) cases requiring hospitalization and developing severe or critical illness has significantly increased. Machine learning and multivariate logistic regression analysis were used to predict risk factors and develop prognostic models for severe COVID-19 in hospitalized children with the Omicron variant in this study. Of the 544 hospitalized children including 243 and 301 in the mild and severe groups, respectively. Fever (92.3%) was the most common symptom, followed by cough (79.4%), convulsions (36.8%), and vomiting (23.2%). The multivariate logistic regression analysis showed that age (1-3 years old, odds ratio (OR): 3.193, 95% confidence interval (CI): 1.778-5.733], comorbidity (OR: 1.993, 95% CI:1.154-3.443), cough (OR: 0.409, 95% CI:0.236-0.709), and baseline neutrophil-to-lymphocyte ratio (OR: 1.108, 95% CI: 1.023-1.200), lactate dehydrogenase (OR: 1.993, 95% CI: 1.154-3.443), blood urea nitrogen (OR: 1.002, 95% CI: 1.000-1.003) and total bilirubin (OR: 1.178, 95% CI: 1.005-3.381) were independent risk factors for severe COVID-19. The area under the curve (AUC) of the prediction models constructed by multivariate logistic regression analysis and machine learning (RandomForest + TomekLinks) were 0.7770 and 0.8590, respectively. The top 10 most important variables of random forest variables were selected to build a prediction model, with an AUC of 0.8210. Compared with multivariate logistic regression, machine learning models could more accurately predict severe COVID-19 in children with Omicron variant infection.

Macrophage-derived PDGF-BB modulates glycolytic enzymes expression and pyroptosis in nucleus pulposus cells via PDGFR-ß/TXNIP pathway.

OBJECTIVE: Intervertebral Disc Degeneration (IVDD) is one of the leading causes of low back pain, significantly impacting both individuals and society. This study aimed to investigate the significance of macrophage infiltration and the role of macrophage-secreted platelet-derived growth factor-BB (PDGF-BB) in IVDD progression. METHODS: To confirm the protective function of macrophage-derived PDGF-BB on nucleus pulposus cells (NPCs), we employed Lysm-Cre transgenic mice to genetically ablate PDGF-B within the myeloid cells. Immunohistochemistry was utilized to detect the expression of glycolytic enzymes and pyroptosis-related proteins during the process of IVDD. Western blot, RT-PCR, ELISA and immunofluorescence were used to detect the protective effect of recombinant PDGF-BB on NPCs. RESULTS: Macrophage-derived PDGF-BB deficiency resulted in the loss of NPCs and the increased ossification of cartilage endplates during lumbar disc degeneration. Also, PDGF-BB deficiency triggered the inhibition of glycolytic enzymes' expression and the activation of pathways related to pyroptosis in the nucleus pulposus. Mechanistically, our results suggest that PDGF-BB predominantly conveys its protective influence on NPCs through the PDGF receptor- beta (PDGFR-ß)/ thioredoxin-interacting protein pathway. CONCLUSIONS: The absence of PDGF-BB originating from macrophages expedites the advancement of IVDD, whereas the application of PDGF-BB treatment holds the potential for retarding intervertebral disc degeneration in the human body.

Membranized Coacervate Microdroplets: from Versatile Protocell Models to Cytomimetic Materials.

Although complex coacervate microdroplets derived from associative phase separation of counter-charged electrolytes have emerged as a broad platform for the bottom-up construction of membraneless, molecularly crowded protocells, the absence of an enclosing membrane limits the construction of more sophisticated artificial cells and their use as functional cytomimetic materials. To address this problem, we and others have recently developed chemical-based strategies for the membranization of preformed coacervate microdroplets. In this Account, we review our recent work on diverse coacervate systems using a range of membrane building blocks and assembly processes. First, we briefly introduce the unusual nature of the coacervate/water interface, emphasizing the ultralow interfacial tension and broad interfacial width as physiochemical properties that require special attention in the judicious design of membranized coacervate microdroplets. Second, we classify membrane assembly into two different approaches: (i) interfacial self-assembly by using diverse surface-active building blocks such as molecular amphiphiles (fatty acids, phospholipids, block copolymers, protein-polymer conjugates) or nano- and microscale objects (liposomes, nanoparticle surfactants, cell fragments, living cells) with appropriate wettability; and (ii) coacervate droplet-to-vesicle reconfiguration by employing auxiliary surface reconstruction agents or triggering endogenous transitions (self-membranization) under nonstoichiometric (charge mismatched) conditions. We then discuss the key cytomimetic behaviors of membranized coacervate-based model protocells. Customizable permeability is achieved by synergistic effects operating between the molecularly crowded coacervate interior and surrounding membrane. In contrast, metabolic-like endogenous reactivity, diffusive chemical signaling, and collective chemical operations occur specifically in protocell networks comprising diverse populations of membranized coacervate microdroplets. In each case, these cytomimetic behaviors can give rise to functional microscale materials capable of promising cell-like applications. For example, immobilizing spatially segregated enzyme-loaded phospholipid-coated coacervate protocells in concentrically tubular hydrogels delivers prototissue-like bulk materials that generate nitric oxide in vitro, enabling platelet deactivation and inhibition of blood clot formation. Alternatively, therapeutic protocells with in vivo vasoactivity, high hemocompatibility, and increased blood circulation times are constructed by spontaneous assembly of hemoglobin-containing cell-membrane fragments on the surface of enzyme-loaded coacervate microdroplets. Higher-order properties such as artificial endocytosis are achieved by using nanoparticle-caged coacervate protocell hosts that selectively and actively capture guest nano- and microscale objects by responses to exogenous stimuli or via endogenous enzyme-mediated reactions. Finally, we discuss the current limitations in the design and programming of membranized coacervate microdroplets, which may help to guide future directions in this emerging research area. Taken together, we hope that this Account will inspire new advances in membranized coacervate microdroplets and promote their application in the development of integrated protocell models and functional cytomimetic materials.

Evaluation of the FPMC respiratory panel for detection of respiratory tract pathogens in nasopharyngeal swab and sputum specimens.

OBJECTIVES: The performance of the new Respiratory Pathogen panel (fluorescent probe melting curve, FPMC) for the qualitative detection of 12 organisms (chlamydia pneumoniae, mycoplasma pneumoniae, adenovirus, influenza A virus, influenza B virus, parainfluenza virus, rhinovirus, etc.) was assessed. METHODS: Prospectively collected nasopharyngeal swab (NPS) and sputum specimens (n = 635) were detected by using the FPMC panel, with the Sanger sequencing method as the comparative method. RESULTS: The overall percent concordance between the FPMC analysis method and the Sanger sequencing method was 100% and 99.66% for NPS and sputum specimens, respectively. The FPMC testified an overall positive percent concordance of 100% for both NPS and sputum specimens. The FPMC analysis method also testified an overall negative percent concordance of 100% and 99.38% for NPS and sputum specimens, respectively. CONCLUSIONS: The FPMC analysis method is a stable and accurate assay for rapid, comprehensive detecting for respiratory pathogens.

Predicting the Metal Mixture Toxicity with a Toxicokinetic-Toxicodynamic Model Considering the Time-Dependent Adverse Outcome Pathways.

Chemicals mainly exist in ecosystems as mixtures, and understanding and predicting their effects are major challenges in ecotoxicology. While the adverse outcome pathway (AOP) and toxicokinetic-toxicodynamic (TK-TD) models show promise as mechanistic approaches in chemical risk assessment, there is still a lack of methodology to incorporate the AOP into a TK-TD model. Here, we describe a novel approach that integrates the AOP and TK-TD models to predict mixture toxicity using metal mixtures (specifically Cd-Cu) as a case study. We preliminarily constructed an AOP of the metal mixture through temporal transcriptome analysis together with confirmatory bioassays. The AOP revealed that prolonged exposure time activated more key events and adverse outcomes, indicating different modes of action over time. We selected a potential key event as a proxy for damage and used it as a measurable parameter to replace the theoretical parameter (scaled damage) in the TK-TD model. This refined model, which connects molecular responses to organism outcomes, effectively predicts Cd-Cu mixture toxicity over time and can be extended to other metal mixtures and even multicomponent mixtures. Overall, our results contribute to a better understanding of metal mixture toxicity and provide insights for integrating the AOP and TK-TD models to improve risk assessment for chemical mixtures.