Autonomous self-burying seed carriers for aerial seeding.

Aerial seeding can quickly cover large and physically inaccessible areas1 to improve soil quality and scavenge residual nitrogen in agriculture2, and for postfire reforestation3-5 and wildland restoration6,7. However, it suffers from low germination rates, due to the direct exposure of unburied seeds to harsh sunlight, wind and granivorous birds, as well as undesirable air humidity and temperature1,8,9. Here, inspired by Erodium seeds10-14, we design and fabricate self-drilling seed carriers, turning wood veneer into highly stiff (about 4.9 GPa when dry, and about 1.3 GPa when wet) and hygromorphic bending or coiling actuators with an extremely large bending curvature (1,854 m-1), 45 times larger than the values in the literature15-18. Our three-tailed carrier has an 80% drilling success rate on flat land after two triggering cycles, due to the beneficial resting angle (25°-30°) of its tail anchoring, whereas the natural Erodium seed's success rate is 0%. Our carriers can carry payloads of various sizes and contents including biofertilizers and plant seeds as large as those of whitebark pine, which are about 11 mm in length and about 72 mg. We compare data from experiments and numerical simulation to elucidate the curvature transformation and actuation mechanisms to guide the design and optimization of the seed carriers. Our system will improve the effectiveness of aerial seeding to relieve agricultural and environmental stresses, and has potential applications in energy harvesting, soft robotics and sustainable buildings.

Proteome-wide structural analysis quantifies structural conservation across distant species.

Traditional evolutionary biology research mainly relies on sequence information to infer evolutionary relationships between genes or proteins. In contrast, protein structural information has long been overlooked, although structures are more conserved and closely linked to the functions than the sequences. To address this gap, we conducted a proteome-wide structural analysis using experimental and computed protein structures for organisms from the three distinct domains, including Homo sapiens (eukarya), Escherichia coli (bacteria), and Methanocaldococcus jannaschii (archaea). We reveal the distribution of structural similarity and sequence identity at the genomic level and characterize the twilight zone, where signals obtained from sequence alignment are blurred and evolutionary relationships cannot be inferred unambiguously. We find that structurally similar homologous protein pairs in the twilight zone account for ∼0.004%-0.021% of all possible protein pair combinations, which translates to ∼8%-32% of the protein-coding genes, depending on the species under comparison. In addition, by comparing the structural homologs, we show that human proteins involved in the energy supply are more similar to their E. coli homologs, whereas proteins relating to the central dogma are more similar to their M. jannaschii homologs. We also identify a bacterial GPCR homolog in the E. coli proteome that displays distinctive domain architecture. Our results shed light on the characteristics of the twilight zone and the origin of different pathways from a protein structure perspective, highlighting an exciting new frontier in evolutionary biology.

Diversity of woodland strawberry inflorescences arises from heterochrony regulated by TERMINAL FLOWER 1 and FLOWERING LOCUS T.

A vast variety of inflorescence architectures have evolved in angiosperms. Here, we analyze the diversity and development of the woodland strawberry (Fragaria vesca) inflorescence. Contrary to historical classifications, we show that it is a closed thyrse: a compound inflorescence with determinate primary monopodial axis and lateral sympodial branches, thus combining features of racemes and cymes. We demonstrate that this architecture is generated by 2 types of inflorescence meristems differing in their geometry. We further show that woodland strawberry homologs of TERMINAL FLOWER 1 (FvTFL1) and FLOWERING LOCUS T (FvFT1) regulate the development of both the racemose and cymose components of the thyrse. Loss of functional FvTFL1 reduces the number of lateral branches of the main axis and iterations in the lateral branches but does not affect their cymose pattern. These changes can be enhanced or compensated by altering FvFT1 expression. We complement our experimental findings with a computational model that captures inflorescence development using a small set of rules. The model highlights the distinct regulation of the fate of the primary and higher-order meristems, and explains the phenotypic diversity among inflorescences in terms of heterochrony resulting from the opposite action of FvTFL1 and FvFT1 within the thyrse framework. Our results represent a detailed analysis of thyrse architecture development at the meristematic and molecular levels.

Glycolysis Induced by METTL14 Is Essential for Macrophage Phagocytosis and Phenotype in Cervical Cancer.

N 6-methyladenosine (m6A) is the most abundant mRNA modification in mammals and it plays a vital role in various biological processes. However, the roles of m6A on cervical cancer tumorigenesis, especially macrophages infiltrated in the tumor microenvironment of cervical cancer, are still unclear. We analyzed the abnormal m6A methylation in cervical cancer, using CaSki and THP-1 cell lines, that might influence macrophage polarization and/or function in the tumor microenvironment. In addition, C57BL/6J and BALB/c nude mice were used for validation in vivo. In this study, m6A methylated RNA immunoprecipitation sequencing analysis revealed the m6A profiles in cervical cancer. Then, we discovered that the high expression of METTL14 (methyltransferase 14, N6-adenosine-methyltransferase subunit) in cervical cancer tissues can promote the proportion of programmed cell death protein 1 (PD-1)-positive tumor-associated macrophages, which have an obstacle to devour tumor cells. Functionally, changes of METTL14 in cervical cancer inhibit the recognition and phagocytosis of macrophages to tumor cells. Mechanistically, the abnormality of METTL14 could target the glycolysis of tumors in vivo and vitro. Moreover, lactate acid produced by tumor glycolysis has an important role in the PD-1 expression of tumor-associated macrophages as a proinflammatory and immunosuppressive mediator. In this study, we revealed the effect of glycolysis regulated by METTL14 on the expression of PD-1 and phagocytosis of macrophages, which showed that METTL14 was a potential therapeutic target for treating advanced human cancers.

A Histopathologic Image Analysis for the Classification of Endocervical Adenocarcinoma Silva Patterns Depend on Weakly Supervised Deep Learning.

Twenty-five percent of cervical cancers are classified as endocervical adenocarcinomas (EACs), which comprise a highly heterogeneous group of tumors. A histopathologic risk stratification system known as the Silva pattern system was developed based on morphology. However, accurately classifying such patterns can be challenging. The study objective was to develop a deep learning pipeline (Silva3-AI) that automatically analyzes whole slide image-based histopathologic images and identifies Silva patterns with high accuracy. Initially, a total of 202 patients with EACs and histopathologic slides were obtained from Qilu Hospital of Shandong University for developing and internally testing the Silva3-AI model. Subsequently, an additional 161 patients and slides were collected from seven other medical centers for independent testing. The Silva3-AI model was developed using a vision transformer and recurrent neural network architecture, utilizing multi-magnification patches, and its performance was evaluated based on a class-specific area under the receiver-operating characteristic curve. Silva3-AI achieved a class-specific area under the receiver-operating characteristic curve of 0.947 for Silva A, 0.908 for Silva B, and 0.947 for Silva C on the independent test set. Notably, the performance of Silva3-AI was consistent with that of professional pathologists with 10 years' diagnostic experience. Furthermore, the visualization of prediction heatmaps facilitated the identification of tumor microenvironment heterogeneity, which is known to contribute to variations in Silva patterns.

m6Aexpress-enet: Predicting the regulatory expression m6A sites by an enet-regularization negative binomial regression model.

As the most abundant mRNA modification, m6A controls and influences many aspects of mRNA metabolism including the mRNA stability and degradation. However, the role of specific m6A sites in regulating gene expression still remains unclear. In additional, the multicollinearity problem caused by the correlation of methylation level of multiple m6A sites in each gene could influence the prediction performance. To address the above challenges, we propose an elastic-net regularized negative binomial regression model (called m6Aexpress-enet) to predict which m6A site could potentially regulate its gene expression. Comprehensive evaluations on simulated datasets demonstrate that m6Aexpress-enet could achieve the top prediction performance. Applying m6Aexpress-enet on real MeRIP-seq data from human lymphoblastoid cell lines, we have uncovered the complex regulatory pattern of predicted m6A sites and their unique enrichment pathway of the constructed co-methylation modules. m6Aexpress-enet proves itself as a powerful tool to enable biologists to discover the mechanism of m6A regulatory gene expression. Furthermore, the source code and the step-by-step implementation of m6Aexpress-enet is freely accessed at https://github.com/tengzhangs/m6Aexpress-enet.

Temperature-Driven Rotation Symmetry-Breaking States in an Atomic Kagome Metal KV3Sb5.

Kagome lattice AV3Sb5 has attracted tremendous interest because it hosts correlated and topological physics. However, an in-depth understanding of the temperature-driven electronic states in AV3Sb5 is elusive. Here we use scanning tunneling microscopy to directly capture the rotational symmetry-breaking effect in KV3Sb5. Through both topography and spectroscopic imaging of defect-free KV3Sb5, we observe a charge density wave (CDW) phase transition from an a0 × a0 atomic lattice to a robust 2a0 × 2a0 superlattice upon cooling the sample to 60 K. An individual Sb-atom vacancy in KV3Sb5 further gives rise to the local Friedel oscillation (FO), visible as periodic charge modulations in spectroscopic maps. The rotational symmetry of the FO tends to break at the temperature lower than 40 K. Moreover, the FO intensity shows an obvious competition against the intensity of the CDW. Our results reveal a tantalizing electronic nematicity in KV3Sb5, highlighting the multiorbital correlation in the kagome lattice framework.

Adenosine 2A receptor antagonists promote lymphocyte proliferation in sepsis by inhibiting Treg expression of PD-L1 in spleen.

The spleen is essential for lymphocyte proliferation, which is associated to sepsis prognosis. Adenosine 2A receptor (A2AR) blocking promotes lymphocyte proliferation in sepsis, however the mechanism is uncertain. Our sepsis cecum ligation perforation model showed that blocking A2AR increased survival and CD4+ cell numbers in a spleen-dependent mechanism. The sequencing of the transcriptome of the spleen indicated alterations in the expression of genes involved in the control of lymphocyte proliferation by inhibiting A2AR, including a reduction in the expression of PD-L1. Flow cytometry analysis of PD-L1 expression intensity in splenic cell subpopulations revealed that the Treg cell subpopulation was the strongest PD-L1-expressing cell population, and Treg PD-L1 expression decreased after blocking A2AR. In vitro activation of A2AR was able to upregulate PD-L1 expression of Treg and boost Treg capacity to limit lymphocyte proliferation, while blockage of PD-L1 partly reduced A2AR-activated Treg's ability to inhibit lymphocyte proliferation. In addition, blocking CREB phosphorylation significantly inhibited A2AR-induced PD-L1 expression. According to the findings of our research, inhibiting A2AR improves the prognosis of sepsis by lowering the level of PD-L1 expression by Treg in the spleen and reducing the inhibition of lymphocyte proliferation.

A Bioinspired Copper-Pair Catalyst in Metal-Organic Frameworks for Molecular Dioxygen Activation and Aerobic Oxidative C-N Coupling.

Open Cu sites were loaded to the UiO-67 metal-organic framework (MOF) skeleton by introduction of flexible Cu-binding pyridylmethylamine (pyma) side chains to the biphenyldicarboxylate linkers. Distance between Cu centers in the MOF pores was tuned by controlling the density of metal-binding side chains. "Interacted" Cu-pair or "isolated" monomeric Cu sites were achieved with high and low (pyma)Cu side chain loading, respectively. Spectroscopic and theoretical studies indicate that "interacted" Cu pairs can effectively bind and activate molecular dioxygen to form Cu2O2 clusters, which showed high catalytic activity for aerobic oxidative C-N coupling. On the contrary, MOF catalyst bearing isolated monomeric Cu sites only showed modest catalytic activity. Enhancement in catalytic performance for the Cu-pair catalyst is attributed to the remote synergistic effect of the paired Cu site, which binds molecular dioxygen and cleaves the O═O bond in a collaborative manner. This work demonstrates that noncovalently interacted metal-pair sites can effectively activate inert small molecules and promote heterogeneous catalytic processes.

Steering CO2 Electroreduction Selectivity U-Turn to Ethylene by Cu-Si Bonded Interface.

Copper (Cu), with the advantage of producing a deep reduction product, is a unique catalyst for the electrochemical reduction of CO2 (CO2RR). Designing a Cu-based catalyst to trigger CO2RR to a multicarbon product and understanding the accurate structure-activity relationship for elucidating reaction mechanisms still remain a challenge. Herein, we demonstrate a rational design of a core-shell structured silica-copper catalyst (p-Cu@m-SiO2) through Cu-Si direct bonding for efficient and selective CO2RR. The Cu-Si interface fulfills the inversion in CO2RR product selectivity. The product ratio of C2H4/CH4 changes from 0.6 to 14.4 after silica modification, and the current density reaches a high of up to 450 mA cm-2. The kinetic isotopic effect, in situ attenuated total reflection Fourier-transform infrared spectra, and density functional theory were applied to elucidate the reaction mechanism. The SiO2 shell stabilizes the *H intermediate by forming Si-O-H and inhibits the hydrogen evolution reaction effectively. Moreover, the direct-bonded Cu-Si interface makes bare Cu sites with larger charge density. Such bare Cu sites and Si-O-H sites stabilized the *CHO and activated the *CO, promoting the coupling of *CHO and *CO intermediates to form C2H4. This work provides a promising strategy for designing Cu-based catalysts with high C2H4 catalytic activity.

m 6 Aexpress-BHM: predicting m6A regulation of gene expression in multiple-groups context by a Bayesian hierarchical mixture model.

As the most abundant RNA modification, N6-methyladenosine (m6A) plays an important role in various RNA activities including gene expression and translation. With the rapid application of MeRIP-seq technology, samples of multiple groups, such as the involved multiple viral/ bacterial infection or distinct cell differentiation stages, are extracted from same experimental unit. However, our current knowledge about how the dynamic m6A regulating gene expression and the role in certain biological processes (e.g. immune response in this complex context) is largely elusive due to lack of effective tools. To address this issue, we proposed a Bayesian hierarchical mixture model (called m6Aexpress-BHM) to predict m6A regulation of gene expression (m6A-reg-exp) in multiple groups of MeRIP-seq experiment with limited samples. Comprehensive evaluations of m6Aexpress-BHM on the simulated data demonstrate its high predicting precision and robustness. Applying m6Aexpress-BHM on three real-world datasets (i.e. Flaviviridae infection, infected time-points of bacteria and differentiation stages of dendritic cells), we predicted more m6A-reg-exp genes with positive regulatory mode that significantly participate in innate immune or adaptive immune pathways, revealing the underlying mechanism of the regulatory function of m6A during immune response. In addition, we also found that m6A may influence the expression of PD-1/PD-L1 via regulating its interacted genes. These results demonstrate the power of m6Aexpress-BHM, helping us understand the m6A regulatory function in immune system.

Development and evolution of the Asteraceae capitulum.

Asteraceae represent one of the largest and most diverse families of plants. The evolutionary success of this family has largely been contributed to their unique inflorescences, capitula that mimic solitary flowers but are typically aggregates of multiple florets. Here, we summarize the recent molecular and genetic level studies that have promoted our understanding of the development and evolution of capitula. We focus on new results on patterning of the enlarged meristem resulting in the iconic phyllotactic arrangement of florets in Fibonacci numbers of spirals. We also summarize the current understanding of the genetic networks regulating the characteristic reproductive traits in the family such as floral dimorphism and differentiation of highly specialized floral organs. So far, developmental studies in Asteraceae are still limited to a very narrow selection of model species. Along with the recent advancements in genomics and phylogenomics, Asteraceae and its relatives provide an outstanding model clade for extended evo-devo studies to exploit the morphological diversity and the underlying molecular networks and to translate this knowledge to the breeding of the key crops in the family.

The hidden diversity of vascular patterns in flower heads.

Vascular systems are intimately related to the shape and spatial arrangement of the plant organs they support. We investigate the largely unexplored association between spiral phyllotaxis and the vascular system in Asteraceae flower heads. We imaged heads of eight species using synchrotron-based X-ray micro-computed tomography and applied original virtual reality and haptic software to explore head vasculature in three dimensions. We then constructed a computational model to infer a plausible patterning mechanism. The vascular system in the head of the model plant Gerbera hybrida is qualitatively different from those of Bellis perennis and Helianthus annuus, characterized previously. Cirsium vulgare, Craspedia globosa, Echinacea purpurea, Echinops bannaticus, and Tanacetum vulgare represent variants of the Bellis and Helianthus systems. In each species, the layout of the main strands is stereotypical, but details vary. The observed vascular patterns can be generated by a common computational model with different parameter values. In spite of the observed differences of vascular systems in heads, they may be produced by a conserved mechanism. The diversity and irregularities of vasculature stand in contrast with the relative uniformity and regularity of phyllotactic patterns, confirming that phyllotaxis in heads is not driven by the vasculature.

Terahertz metalens for generating multi-polarized focal points and images with uniform intensity distributions.

Metasurfaces have provided a flexible platform for designing ultracompact metalenses with unusual functionalities. However, traditional multi-foci metalenses are limited to generating circularly polarized (CP) or linearly polarized (LP) focal points, and the intensity distributions are always inhomogeneous/chaotical between the multiple focal points. Here, an inverse design approach is proposed to optimize the in-plane orientation of each meta-atom in a terahertz (THz) multi-foci metalens that can generate multi-polarized focal points with nearly uniform intensity distributions. As a proof-of-principle example, we numerically and experimentally demonstrate an inversely designed metalens for simultaneously generating multiple CP- and LP-based focal points with homogeneous intensity distributions, leading to a multi-polarized image (rather than the holography). Furthermore, the multi-channel and multi-polarized images consisting of multiple focal points with homogeneous intensity distributions are also numerically demonstrated. The unique approach for inversely designing multi-foci metalens that can generate multi-polarized focal points and images with uniform intensity distributions will enable potential applications in imaging and sensing.

Beam-shaped femtosecond laser printing of quasi-capsule-shaped holographic terahertz metasurfaces.

Terahertz (THz) metasurfaces have opened up a new avenue for the THz wavefront modulation. However, high-efficient and low-cost fabrication of THz metasurfaces remains a great challenge today. Here, quasi-capsule-shaped polarization-multiplexed holographic THz metasurfaces were printed by a beam-shaped femtosecond laser. The laser beam was spatially modulated by holograms of optimized cylindrical lens loaded on a spatial light modulator (SLM). The size of quasi-capsule apertures can be exquisitely and flexibly controlled by adjusting the focal length in holograms, pulse energy, and pulse number. Based on near-field diffraction and Burch encoding, an array of 100 × 100 basic unit apertures were initially designed, and a polarization-multiplexed THz metasurface was finally printed with a dimension of 8 mm × 8 mm. The function of polarization multiplexing was demonstrated, by which two kinds of images were reconstructed in response to X and Y-polarization THz waves, respectively. The present work highlights a great leap in fabrication method for THz metasurfaces and hopefully stimulates the development of miniaturized and integrated THz systems.

APLN promotes the proliferation, migration, and glycolysis of cervical cancer through the PI3K/AKT/mTOR pathway.

Apelin (APLN) is an endogenous ligand of the G protein-coupled receptor APJ (APLNR). APLN has been implicated in the development of multiple tumours. Herein, we determined the effect of APLN on the biological behaviour and underlying mechanisms of cervical cancer. The expression and survival curves of APLN were determined using Gene Expression Profiling Interactive Analysis. The cellular functions of APLN were detected using CCK-8, clone formation, EdU, Transwell assays, flow cytometry, and seahorse metabolic analysis. The underlying mechanisms were elucidated using gene set enrichment analysis and Western blotting. APLN was upregulated in the samples of patients with cervical cancer and is associated with poor prognosis. APLN knockdown decreased the proliferation, migration, and glycolysis of cervical cancer cells. The opposite results were observed when APLN was overexpressed. Mechanistically, we determined that APLN was critical for activating the PI3K/AKT/mTOR pathway via APLNR. APLN receptor inhibitor ML221 reversed the effect of APLN overexpression on cervical cancer cells. Treatment with LY294002, the PI3K inhibitor, drastically reversed the oncological behaviour of APLN-overexpressing C-33A cells. APLN promoted the proliferation, migration, and glycolysis of cervical cancer cells via the PI3K/AKT/mTOR pathway.

MicroRNA transcriptome analysis reveals the immune regulatory mechanism of Crassostrea hongkongesis against Vibrio harveyi infection.

MicroRNAs (miRNAs) are small non-coding RNA molecules that modulate target-genes expression and play crucial roles in post-transcriptional regulation and immune system regulation. The Hong Kong oyster (Crassostrea hongkongesis), as the main marine aquaculture shellfish in the South China Sea, not only has high economic and ecological value, but also is an ideal model for conducting research on pathogen host interaction. Vibrio harveyi, a Gram negative luminescent marine bacterium, is widely distributed in coastal water environments and can cause large-scale death of C. hongkongesis. However, little in formation is available on the immune regulatory mechanisms of C. hongkongesis infected with V. harveyi. Therefore, we performed microRNA transcriptome analysis for elucidating the immunoregulation mechanism of C. hongkongesis infected with V. harveyi. The results show that a total of 308468208 clean reads and 288371159 clean tags were obtained. 222 differentially expressed miRNAs were identified. A total of 388 target genes that were differentially expressed and negatively correlated with miRNA expression were predicted by 222 DEmiRs. GO enrichment analysis of 388 DETGs showed that they were mainly enriched in the immune-related term of membrane-bounded vesicle, endocytic vesicle lumen, antigen processing and presentation of exogenous peptide antigen via MHC class I, antigen processing and presentation of peptide antigen via MHC class I, and other immune-related term. KEGG enrichment analysis showed that DETGs were mainly enriched in the Complement and coagulation cascades, Herpes simplex virus 1 infection, Bacterial invasion of epithelial cells, Antigen processing and presentation and NOD-like receptor signaling pathway. The 16 key DEmiRs and their target genes form a regulatory network for seven immune-related pathways. These results suggest that V. harveyi infection induces a complex miRNA response with wide-ranging effects on immune gene expression in the C. hongkongesis. This study explored the immune response of C. hongkongesis to V. harveyi infection at the level of miRNAs, which provides new ideas for the healthy culture and selective breeding of C. hongkongesis.

Pt-S Bond-Enabled Temperature-Dependent Phosphorescence in S-Heteroaryl Tetradentate Pt(S

This study presents the rare examples of S-heteroaryl tetradentate Pt(S^C^N^O) luminescent complexes (PtSZ and PtSZtBu) containing a Pt-S bond. The presence of the Pt-S bond allows the novel Pt(S^C^N^O) complexes to exhibit temperature-dependent phosphorescent emission behavior. The PtSZtBu exhibits dual-emission phenomena and biexponential transient decay spectra above 250 K, indicating the presence of two minimal excited states in the potential energy surface (PES) of the T1 state. Through complementary experimental and computational studies, we have identified changes in orbital composition between Pt(dxy)-S(px) and Pt(dyz)-S(pz) in excited states with increasing temperature. This results in two energy minima, enabling the excited states to decay selectively and radiatively at different temperatures. Consequently, this leads to remarkable steady-state and transient emission spectra changes. Our work not only provides valuable insights for the development of novel Pt-S bond-based tetradentate Pt(II) complexes but also enhances our understanding of the distinctive properties governed by the Pt-S bond.

Inhibition of RXRA-mediated PLA2G2A improves delirium in COPD mice by regulating endoplasmic reticulum stress pathway and inhibiting cell apoptosis.

Delirium is a common psychiatric complication of chronic obstructive pulmonary disease (COPD). The relief of delirium is considered one of the beneficial ways to treat COPD. However, there are currently no specific drugs that alleviate delirium in COPD patients. Our research aimed to elucidate the specific mechanisms underlying delirium in COPD mice, while also seeking more effective therapeutic targets. In our study, bioinformatics analysis and qRT PCR were used to identify key factors in the development of delirium in COPD animal models. Open field and elevated plus maze tests were used to detect delirium in mice. Tunel staining and HE staining were used to analyze the apoptosis of mouse hippocampus cells. EdU and CCK-8 experiments were used to analyze PC-12 cells vitality and proliferation. JASPAR online database, dual luciferase reporting experiments, ChIP experiments, and IF staining were used to analyze the interaction between RXRA and PLA2G2A. RXRA is highly expressed in the brain tissue of COPD mice with delirium symptoms. The downregulation of RXRA inhibits the delirium state in COPD mice. This is mainly due to the reduction of endoplasmic reticulum stress and cell apoptosis by inhibiting the expression of RXRA. In addition, we also confirmed that RXRA is a transcription factor of PLA2G2A. RXRA has an inhibitory effect on the expression of PLA2G2A. In vitro experiments have confirmed that inhibition of the RXRA/PLA2G2A axis reduces cell apoptosis, thereby alleviating the occurrence and development of delirium in COPD mice. Inhibition of the RXRA/PLA2G2A axis reduces endoplasmic reticulum stress and cell apoptosis. This process alleviates the development of delirium in COPD mice.

Reaction kinetics of lithium-sulfur batteries with a polar Li-ion electrolyte: modeling of liquid phase and solid phase processes.

The present investigation fits the reaction kinetics of a lithium-sulfur (Li-S) battery with polar electrolyte employing a novel two-phase continuum multipore model. The continuum two-phase model considers processes in both the liquid electrolyte phase and the solid precipitates phase, where the diffusion coefficients of the Li+ ions in a solvent-softened solid state are determined from molecular dynamics simulations. Solubility experiments yield the saturation concentration of sulfur and lithium sulfides in the polar electrolyte employed in this study. The model describes the transport of dissolved molecular and ion species in pores of different size in solvated or desolvated form, depending on pore size. The Li-S reaction model in this study is validated for electrolyte 1 M LiPF6 in EC/DMC. It includes seven redox reactions and two cyclic non-electrochemical reactions in the cathode, and the lithium redox reaction at the anode. Electrochemical reactions are assumed to take place in the electrolyte solution or the solid state and cyclic reactions are assumed to take place in the liquid electrolyte phase only. The determination of the reaction kinetics parameters takes place via fitting the model predictions with experimental data of a cyclic voltammetry cycle with in operando UV-vis spectroscopy.